solveFluid.H File Reference

Include dependency graph for solveFluid.H:

Go to the source code of this file.

Functions
fvOptions	correct (he)
thermo	correct ()
	if (Y.size())
volScalarField	Yt (IOobject("Yt", runTime.timeName(), mesh), mesh, dimensionedScalar("Yt", dimless, 0))
	forAll (Y, i)
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)
	if (coupled)
fvOptions	correct (he)
thermo	correct ()

Variables
Info<< "\Solving for fluid region "<< fluidRegions[i].name()<< endl;if(finalIter){ mesh.data().setFinalIteration(true);}if(frozenFlow){# 1 "/home/chef2/andy/OpenFOAM/release/v2606/OpenFOAM-v2606/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/EEqn.H" 1{ volScalarField &he=thermo.he();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);if(coupled) { fvMatrixAssemblyPtr->	addFvMatrix (EEqn)
	else
Info<< "Min/max T:"<< min(thermo.T()).value()<< ' '<< max(thermo.T()).value()<< endl;}}# 10 "/home/chef2/andy/OpenFOAM/release/v2606/OpenFOAM-v2606/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H" 2}else{ if(oCorr==0) { }# 1 "/home/chef2/andy/OpenFOAM/release/v2606/OpenFOAM-v2606/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/UEqn.H" 1 MRF.correctBoundaryVelocity(U);UEqn=(fvm::ddt(rho, U)+fvm::div(phi, U)+MRF.DDt(rho, U)+turbulence.divDevRhoReff(U)==fvOptions(rho, U));UEqn.relax();fvOptions.constrain(UEqn);if(momentumPredictor) { solve(UEqn==fvc::reconstruct((- ghf *fvc::snGrad(rho) - fvc::snGrad(p_rgh)) *mesh.magSf()), U.select(finalIter));fvOptions.correct(U);K=0.5 *magSqr(U);}# 19 "/home/chef2/andy/OpenFOAM/release/v2606/OpenFOAM-v2606/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H" 2# 1 "/home/chef2/andy/OpenFOAM/release/v2606/OpenFOAM-v2606/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/YEqn.H" 1tmp< fv::convectionScheme< scalar > >	mvConvection (nullptr)
	Qdot = reaction.Qdot()

Function Documentation

correct() [1/4]

fvOptions correct

(

he

)

correct() [2/4]

thermo correct

(

)

if() [1/2]

if

(

Y.

size()

)

Definition at line 106 of file solveFluid.H.

Yt()

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

forAll()

forAll	(	Y	,
		i
	)

Definition at line 130 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()

EEqn relax

(

)

constrain()

fvOptions constrain

(

EEqn

)

if() [2/2]

if

(

coupled

)

Definition at line 30 of file EEqn.H.

References EEqn(), and fvMatrixAssemblyPtr.

Here is the call graph for this function:

correct() [3/4]

fvOptions correct

(

he

)

correct() [4/4]

thermo correct

(

)

Variable Documentation

addFvMatrix

Info<< "\nSolving for fluid region " << fluidRegions[i].name() << endl;if (finalIter){ mesh.data().setFinalIteration(true);}if (frozenFlow){# 1 "/home/chef2/andy/OpenFOAM/release/v2606/OpenFOAM-v2606/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/EEqn.H" 1{ volScalarField& he = thermo.he(); 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); if (coupled) { fvMatrixAssemblyPtr-> addFvMatrix(EEqn)

Definition at line 42 of file solveFluid.H.

else

else

Initial value:

{
        EEqn.solve(he.select(finalIter))

Definition at line 45 of file solveFluid.H.

mvConvection

Info<< "Min/max T:" << min(thermo.T()).value() << ' ' << max(thermo.T()).value() << endl; }}# 10 "/home/chef2/andy/OpenFOAM/release/v2606/OpenFOAM-v2606/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H" 2}else{ if (oCorr == 0) { }# 1 "/home/chef2/andy/OpenFOAM/release/v2606/OpenFOAM-v2606/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/UEqn.H" 1 MRF.correctBoundaryVelocity(U); UEqn = ( fvm::ddt(rho, U) + fvm::div(phi, U) + MRF.DDt(rho, U) + turbulence.divDevRhoReff(U) == fvOptions(rho, U) ); UEqn.relax(); fvOptions.constrain(UEqn); if (momentumPredictor) { solve ( UEqn == fvc::reconstruct ( ( - ghf*fvc::snGrad(rho) - fvc::snGrad(p_rgh) )*mesh.magSf() ), U.select(finalIter) ); fvOptions.correct(U); K = 0.5*magSqr(U); }# 19 "/home/chef2/andy/OpenFOAM/release/v2606/OpenFOAM-v2606/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H" 2# 1 "/home/chef2/andy/OpenFOAM/release/v2606/OpenFOAM-v2606/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/YEqn.H" 1tmp<fv::convectionScheme<scalar> > mvConvection(nullptr)

Definition at line 104 of file solveFluid.H.

Referenced by if().

Qdot

Qdot = reaction.Qdot()

Definition at line 122 of file solveFluid.H.