|February 23, 2015, 01:56||
continuity source term in momentum Eq
Join Date: Sep 2009
Location: Tehran, Iran
Blog Entries: 1Rep Power: 17
would you please help me to figure out following source codes:
in interPhaseChangeFoam, UEqn.H:
fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U) + turbulence->divDevRhoReff(rho, U) );
fvVectorMatrix UcEqn ( fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc) - fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc) + continuousPhaseTurbulence->divDevRhoReff(Uc) == (1.0/rhoc)*cloudSU );
as much as i understand both "Sp" lines should be mathematically zero , they are continuity terms which multiply with U, my question is why they add to momentum equation and why it considered as negative source terms?
|February 24, 2015, 02:09||
Join Date: Jul 2009
Blog Entries: 1Rep Power: 10
I guess the additional term is due to variable density (rho). Had it been incompressible fluid it would have been zero. Since, for compressible flows these variables cannot be ignored. Please rectify me if I miss anything as I am not well versed with compressible flows.
|March 5, 2015, 18:22||
Join Date: Jun 2012
Location: Torino, Italy
Posts: 743Rep Power: 9
Pls see the image:
The reason maybe in this way, its conservative.
This is like twoPhaseEulerFoam's code. In openfoam 22x, this kind of equation is implemented in non-conservative way, but in OpenFOAM 23x, its implemented in conservative way. U can compare the code to check it out.
|March 13, 2015, 05:39||
Join Date: Feb 2014
Posts: 25Rep Power: 4
I am not sure but perhaps this has something to do with the diagonal dominance of the system? Although this would only make sense if the Sp term would be negative before converging to zero... if it would be slightly positive before converging to zero it would negatively affect the diagonal dominance...
Again, I am not sure about this, one of the remaining mysteries I am trying to figure out.
|continuity source, momentum|
|Thread||Thread Starter||Forum||Replies||Last Post|
|implicit - scalar product source term in momentum equation||vinch||OpenFOAM Running, Solving & CFD||0||October 28, 2014 15:57|
|Trouble compiling utilities using source-built OpenFOAM||Artur||OpenFOAM Programming & Development||14||October 29, 2013 11:59|
|dynamic Mesh is faster than MRF????||sharonyue||OpenFOAM Running, Solving & CFD||14||August 26, 2013 07:47|
|"parabolicVelocity" in OpenFoam 2.1.0 ?||sawyer86||OpenFOAM Running, Solving & CFD||21||February 7, 2012 12:44|
|DecomposePar links against liblamso0 with OpenMPI||jens_klostermann||OpenFOAM Bugs||11||June 28, 2007 17:51|