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())
volScalarField	Yt (IOobject("Yt", runTime.timeName(), mesh), mesh, dimensionedScalar("Yt", dimless, 0))
	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)+Qdot+fvOptions(rho, he))
EEqn	relax ()
fvOptions	constrain (EEqn)
EEqn	solve ()
fvOptions	correct (he)
thermo	correct ()
	while (pimple.correctNonOrthogonal())
	if (!mesh.steady())
	if (thermo.dpdt())
	if (!mesh.steady())
	if (adjustMass)
	if (adjustMass||pLimited)
	if (adjustMass||pLimited||mesh.steady()||pimple.simpleRho())
	if (mesh.steady() &&!pimple.transonic())

Variables
	else
fvVectorMatrix &	UEqn = tUEqn.ref()
	Qdot = reaction.Qdot()
Info<< "Min/max T:"<< min(thermo.T()).value()<< ' '<< max(thermo.T()).value()<< endl;}# 14 "/home/ubuntu/OpenFOAM-6/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H" 2 while(pimple.correct()) {# 1 "/home/ubuntu/OpenFOAM-6/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=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)));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/10]

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/10]

if

(

Y.

size()

)

Definition at line 4 of file solveFluid.H.

Yt()

volScalarField Yt	(	IOobject("Yt", runTime.timeName(), mesh)	,
		mesh	,
		dimensionedScalar("Yt", dimless, 0)
	)

forAll()

forAll	(	Y	,
		i
	)

Definition at line 28 of file solveFluid.H.

References rho, and thermo.

if() [3/10]

if

(

Y.

size()

)

Definition at line 57 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)+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 102 of file solveFluid.H.

if() [4/10]

if

(

!mesh.

steady()

)

Definition at line 141 of file solveFluid.H.

if() [5/10]

if

(

thermo.

dpdt()

)

Definition at line 148 of file solveFluid.H.

if() [6/10]

if

(

!mesh.

steady()

)

Definition at line 154 of file solveFluid.H.

if() [7/10]

if

(

adjustMass

)

Definition at line 169 of file solveFluid.H.

if() [8/10]

if

(

adjustMass||

pLimited

)

Definition at line 176 of file solveFluid.H.

if() [9/10]

if

(

adjustMass||pLimited||mesh.

steady)||pimple.simpleRho(

)

Definition at line 182 of file solveFluid.H.

if() [10/10]

if

(

mesh.

steady) &&!pimple.transonic(

)

Definition at line 187 of file solveFluid.H.

Variable Documentation

else [1/2]

else

Initial value:

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

    }

# 1 "/home/ubuntu/OpenFOAM-6/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.

Qdot

Qdot = reaction.Qdot()

Definition at line 20 of file solveFluid.H.

psip0

Info<< "Min/max T:" << min(thermo.T()).value() << ' ' << max(thermo.T()).value() << endl;}# 14 "/home/ubuntu/OpenFOAM-6/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H" 2 while (pimple.correct()) {# 1 "/home/ubuntu/OpenFOAM-6/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 = 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)) ); 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 100 of file solveFluid.H.

p

p = p_rgh + rho*gh

Definition at line 138 of file solveFluid.H.

else [2/2]

else

Initial value:

{

}


const bool pLimited = pressureControl.limit(p)

Definition at line 160 of file solveFluid.H.