|
[Sponsors] |
Implementation of advection in alphaEqn.H in settlingFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 17, 2015, 00:34 |
Implementation of advection in alphaEqn.H in settlingFoam
|
#1 |
New Member
Xavier Pivan
Join Date: May 2015
Posts: 10
Rep Power: 10 |
Dear Foamers,
I wish to add an advection term (div(U*Alpha)) in the settlingFoam solver, alphaEqn.H. I try different option but no one works. Does anyone of you have an idea how to do that ? The problem lies on the fact that U and Alpha are matrix I think. volScalarField divU(fvc::div(phi/fvc::interpolate(rho))); fvScalarMatrix AlphaEqn ( fvm::ddt(rho, Alpha) + fvm::div(phiAlpha, Alpha) + fvm::div(divU,Alpha) //I tried this as divU is showned just up - fvm::laplacian(mut, Alpha) ); Help would be appreciate, Xavier |
|
July 17, 2015, 02:09 |
|
#2 |
Senior Member
Hesam
Join Date: Feb 2015
Posts: 139
Rep Power: 11 |
Hi Xavier,
I have the same problem with a small difference.when I use phi/fvc::interpolate(rho), result of alpha become with high gradients in some locations and I don't know the reason. But in your case i think if your solver even make succesfully it has dimentional problem in divU. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
alphaEqn.H in twoPhaseEulerFoam | cheng1988sjtu | OpenFOAM Bugs | 15 | May 1, 2016 16:12 |
ENO advection implementation | mpash | Main CFD Forum | 0 | January 23, 2014 23:20 |
Use advection step and Dispersion step to solve advection diffusion equation | cooljd | Main CFD Forum | 0 | October 15, 2012 15:44 |
Total variation diminishing scheme for advection or advection dispersion equation | cooljd | Main CFD Forum | 7 | July 27, 2012 08:32 |
Implementation issues of fvSchemes / laplacianScheme, in particular gaussLaplacianSch | thomek | OpenFOAM Programming & Development | 0 | October 18, 2010 05:10 |