solveFluid.H File Reference

Include dependency graph for solveFluid.H:

Go to the source code of this file.

Functions
	if (pimple.frozenFlow())
tmp< fvVectorMatrix >	tUEqn (fvm::ddt(rho, U)+fvm::div(phi, U)+MRF.DDt(rho, U)+turbulence.divDevRhoReff(U)==fvOptions(rho, U))
UEqn	relax ()
fvOptions	constrain (UEqn)
fvOptions	correct (U)
tmp< fv::convectionScheme< scalar > >	mvConvection (nullptr)
	if (Y.size())
tmp< volScalarField >	Yt (nullptr)
	if (Y.size())
	forAll (Y, i)
	if (Y.size())
fvScalarMatrix	EEqn (fvm::ddt(rho, he)+fvm::div(phi, he)+fvc::ddt(rho, K)+fvc::div(phi, K)+(he.name()=="e" ? fvc::div(fvc::absolute(phi/fvc::interpolate(rho), U), p, "div(phiv,p)") :-dpdt) - fvm::laplacian(turbulence.alphaEff(), he)==rho *(U &g)+rad.Sh(thermo, he)+reaction.Qdot()+fvOptions(rho, he))
EEqn	relax ()
fvOptions	constrain (EEqn)
EEqn	solve ()
fvOptions	correct (he)
thermo	correct ()
	while (pimple.correctNonOrthogonal())
pressureControl	limit (p)
	if (!mesh.steady())
	if (thermo.dpdt())
	if (!mesh.steady())
	if (pimples.pimpleTurbCorr(i))
	if (!mesh.steady() &&pimples.finalPimpleIter())

Variables
	else
fvVectorMatrix &	UEqn = tUEqn.ref()
Info<< "Min/max T:"<< min(thermo.T()).value()<< ' '<< max(thermo.T()).value()<< endl;}# 14 "/home/ubuntu/OpenFOAM-7/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H" 2 while(pimple.correct()) {# 1 "/home/ubuntu/OpenFOAM-7/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H" 1if(!mesh.steady() &&!pimple.simpleRho()){ rho=thermo.rho();}volScalarField rAU("rAU", 1.0/UEqn.A());surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho *rAU));volVectorField HbyA(constrainHbyA(rAU *UEqn.H(), U, p_rgh));if(pimple.nCorrPiso()<=1){ tUEqn.clear();}surfaceScalarField phig(-rhorAUf *ghf *fvc::snGrad(rho) *mesh.magSf());surfaceScalarField phiHbyA("phiHbyA", fvc::flux(rho *HbyA)+MRF.zeroFilter(rhorAUf *fvc::ddtCorr(rho, U, phi)));MRF.makeRelative(fvc::interpolate(rho), phiHbyA);const bool closedVolume=mesh.steady() &&adjustPhi(phiHbyA, U, p_rgh);const bool adjustMass=closedVolume &&!thermo.incompressible();phiHbyA+=phig;constrainPressure(p_rgh, rho, U, phiHbyA, rhorAUf, MRF);{ fvScalarMatrix p_rghEqnComp(fvc::ddt(rho)+psi *correction(fvm::ddt(p_rgh))==fvOptions(psi, p_rgh, rho.name()));if(pimple.transonic()) { surfaceScalarField phid("phid",(fvc::interpolate(psi)/fvc::interpolate(rho)) *phiHbyA);phiHbyA -=fvc::interpolate(psi *p_rgh) *phiHbyA/fvc::interpolate(rho);p_rghEqnComp+=fvm::div(phid, p_rgh);} tmp< volScalarField >	psip0 (mesh.steady() ? tmp< volScalarField >() :psi *p)
	p = p_rgh + rho*gh
	else

Function Documentation

if() [1/9]

if

(

pimple.

frozenFlow()

)

Definition at line 1 of file solveFluid.H.

tUEqn()

tmp<fvVectorMatrix> tUEqn

(

fvm::ddt(rho, U)+fvm::div(phi, U)+MRF.DDt(rho, U)+turbulence.divDevRhoReff(U)

= =fvOptions(rho, U)

)

relax() [1/2]

UEqn relax

(

)

constrain() [1/2]

fvOptions constrain

(

UEqn

)

correct() [1/3]

fvOptions correct

(

U

)

mvConvection()

tmp<fv::convectionScheme<scalar> > mvConvection

(

nullptr

)

if() [2/9]

if

(

Y.

size()

)

Definition at line 4 of file solveFluid.H.

Yt()

tmp<volScalarField> Yt

(

nullptr

)

if() [3/9]

if

(

Y.

size()

)

Definition at line 23 of file solveFluid.H.

References turbulence.

forAll()

forAll	(	Y	,
		i
	)

Definition at line 36 of file solveFluid.H.

if() [4/9]

if

(

Y.

size()

)

Definition at line 65 of file solveFluid.H.

EEqn()

fvScalarMatrix EEqn

(

fvm::ddt(rho, he)+fvm::div(phi, he)+fvc::ddt(rho, K)+fvc::div(phi, K)+(he.name()=="e" ? fvc::div(fvc::absolute(phi/fvc::interpolate(rho), U), p, "div(phiv,p)") :-dpdt) - fvm::laplacian(turbulence.alphaEff(), he)

= =rho *(U &g)+rad.Sh(thermo, he)+reaction.Qdot()+fvOptions(rho, he)

)

Initial value:

{
    volScalarField& he = thermo.he()

relax() [2/2]

EEqn relax

(

)

constrain() [2/2]

fvOptions constrain

(

EEqn

)

solve()

EEqn solve

(

)

correct() [2/3]

fvOptions correct

(

he

)

correct() [3/3]

thermo correct

(

)

while()

while

(

pimple.

correctNonOrthogonal()

)

Definition at line 104 of file solveFluid.H.

limit()

pressureControl limit

(

p

)

if() [5/9]

if

(

!mesh.

steady()

)

Definition at line 145 of file solveFluid.H.

if() [6/9]

if

(

thermo.

dpdt()

)

Definition at line 152 of file solveFluid.H.

if() [7/9]

if

(

!mesh.

steady()

)

Definition at line 158 of file solveFluid.H.

if() [8/9]

if

(

pimples.

pimpleTurbCorri

)

Definition at line 22 of file solveFluid.H.

if() [9/9]

if

(

!mesh.

steady) &&pimples.finalPimpleIter(

)

Definition at line 27 of file solveFluid.H.

References rho, and thermo.

Variable Documentation

else [1/2]

else

Initial value:

{
    if (!mesh.steady() && pimples.nCorrPimple() <= 1)
    {

    }

# 1 "/home/ubuntu/OpenFOAM-7/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/UEqn.H" 1
    

    MRF.correctBoundaryVelocity(U)

Definition at line 6 of file solveFluid.H.

UEqn

fvVectorMatrix& UEqn = tUEqn.ref()

Definition at line 25 of file solveFluid.H.

psip0

Info<< "Min/max T:" << min(thermo.T()).value() << ' ' << max(thermo.T()).value() << endl;}# 14 "/home/ubuntu/OpenFOAM-7/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H" 2 while (pimple.correct()) {# 1 "/home/ubuntu/OpenFOAM-7/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H" 1if (!mesh.steady() && !pimple.simpleRho()){ rho = thermo.rho();}volScalarField rAU("rAU", 1.0/UEqn.A());surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU));volVectorField HbyA(constrainHbyA(rAU*UEqn.H(), U, p_rgh));if (pimple.nCorrPiso() <= 1){ tUEqn.clear();}surfaceScalarField phig(-rhorAUf*ghf*fvc::snGrad(rho)*mesh.magSf());surfaceScalarField phiHbyA( "phiHbyA", fvc::flux(rho*HbyA) + MRF.zeroFilter(rhorAUf*fvc::ddtCorr(rho, U, phi)));MRF.makeRelative(fvc::interpolate(rho), phiHbyA);const bool closedVolume = mesh.steady() && adjustPhi(phiHbyA, U, p_rgh);const bool adjustMass = closedVolume && !thermo.incompressible();phiHbyA += phig;constrainPressure(p_rgh, rho, U, phiHbyA, rhorAUf, MRF);{ fvScalarMatrix p_rghEqnComp ( fvc::ddt(rho) + psi*correction(fvm::ddt(p_rgh)) == fvOptions(psi, p_rgh, rho.name()) ); if (pimple.transonic()) { surfaceScalarField phid ( "phid", (fvc::interpolate(psi)/fvc::interpolate(rho))*phiHbyA ); phiHbyA -= fvc::interpolate(psi*p_rgh)*phiHbyA/fvc::interpolate(rho); p_rghEqnComp += fvm::div(phid, p_rgh); } tmp<volScalarField> psip0(mesh.steady() ? tmp< volScalarField >() :psi *p)

Definition at line 102 of file solveFluid.H.

p

p = p_rgh + rho*gh

Definition at line 140 of file solveFluid.H.

else [2/2]

else

Initial value:

{

}




if (adjustMass)
{
    p += (initialMass - fvc::domainIntegrate(thermo.rho()))
        /fvc::domainIntegrate(psi);
    p_rgh = p - rho*gh;
    p.correctBoundaryConditions();
}


if (adjustMass || mesh.steady() || pimple.simpleRho())
{
    rho = thermo.rho();
}

if (mesh.steady() && !pimple.transonic())
{
    rho.relax();
}

Info<< "Min/max rho:" << min(rho).value() << ' '
    << max(rho).value() << endl

Definition at line 164 of file solveFluid.H.