CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   MULESTemplates.C for compressible two phase flow (http://www.cfd-online.com/Forums/openfoam/92334-mulestemplates-c-compressible-two-phase-flow.html)

nimasam September 10, 2011 15:23

MULESTemplates.C for compressible two phase flow
 
i have a two phase flow with variable density which change with pressure, temperature and chemical conversion!
look following step to derive a formulation for alphal1

1)ddt(rho) + div.(rho,U) = 0

and rho = rho1*alpha1+rho2*(1-alpha1)=(rho1-rho2)*alpha1+rho2

so:

2)ddt(rho_*alpha1) + div.(rho_,U) + ddt(rho2) + div.(rho2,U) = 0

and rho_ = rho1 - rho2

from MULESTemplates.C :

fvScalarMatrix psiConvectionDiffusion
(
fvm::ddt(rho, psi)
+ fv::gaussConvectionScheme<scalar>(mesh, phi, UDs).fvmDiv(phi, psi)
//- fv::gaussLaplacianScheme<scalar, scalar>(mesh, CDs, snGrads)
//.fvmLaplacian(Dpsif, psi)
- fvm::Sp(Sp, psi)
- Su
);

and from incompressibleInterFoam or interPhaseChangeFoam :
MULES::implicitSolve(oneField(), alphalpha1, phi, phiAlphalpha1, Sp, Su, 1, 0);

can i use rho_ instead of oneField() in above structure?


All times are GMT -4. The time now is 01:27.