[newOFuser]

Hi

Could someone please translate the XiFoam code below into equation form? I want to check how they relate to the equations from Weller's paper:

Basically I want to see how Us(.)(grad (rho Xi)) in the transport equation for Xi equals
fvm::div(phi + phiXi, Xi, "div(phiXi,Xi)") - fvm::Sp(fvc::div(phiXi), Xi)


surfaceScalarField phiXi = phiSt –
fvc::interpolate(fvc::laplacian(turbulence->alphaEff(), b)/mgb)*nf
+ fvc::interpolate(rho)*fvc::interpolate(Su*(1.0/Xi - Xi))*nf;

fvScalarMatrix XiEqn
(
fvm::ddt(rho, Xi)
+ fvm::div(phi + phiXi, Xi, "div(phiXi,Xi)")
- fvm::Sp(fvc::div(phiXi), Xi)
==
rho*R
- fvm::Sp(rho*(R - G), Xi)
- fvm::Sp(rho*max(sigmat-sigmas,dimensionedScalar("0",sigmat.dimensions(), 0)),Xi)
);
XiEqn.relax();
XiEqn.solve();
--------------------
Thanks so much!
ak

[elvis]

Hi,

did you read http://www.tfd.chalmers.se/~hani/kur...sariReport.pdf

Tutorial XiFoam
Version: OpenFOAM-1.7.x
Author:

Ehsan Yasari