Foam::scalar Foam::solidRegionDiffNo ( const fvMesh& mesh, const Time& runTime, const volScalarField& Cprho, const volScalarField& kappa ) { scalar DiNum = 0.0; scalar meanDiNum = 0.0; //- Take care: can have fluid domains with 0 cells so do not test for // zero internal faces. surfaceScalarField kapparhoCpbyDelta ( mesh.surfaceInterpolation::deltaCoeffs() * fvc::interpolate(kappa) / fvc::interpolate(Cprho) ); DiNum = gMax(kapparhoCpbyDelta.internalField())*runTime.deltaT().value(); meanDiNum = (average(kapparhoCpbyDelta)).value()*runTime.deltaT().value(); Info<< "Region: " << mesh.name() << " Diffusion Number mean: " << meanDiNum << " max: " << DiNum << endl; return DiNum; }
Foam::scalar Foam::compressibleCourantNo( const fvMesh& mesh, const Time& runTime, const volScalarField& rho, const surfaceScalarField& phi ) { scalar CoNum = 0.0; scalar meanCoNum = 0.0; //- Can have fluid domains with 0 cells so do not test. //if (mesh.nInternalFaces()) { surfaceScalarField SfUfbyDelta = mesh.surfaceInterpolation::deltaCoeffs() * mag(phi) / fvc::interpolate(rho); CoNum = max(SfUfbyDelta/mesh.magSf()) .value()*runTime.deltaT().value(); meanCoNum = (sum(SfUfbyDelta)/sum(mesh.magSf())) .value()*runTime.deltaT().value(); } Info<< "Region: " << mesh.name() << " Courant Number mean: " << meanCoNum << " max: " << CoNum << endl; return CoNum; }
Foam::scalar Foam::compressibleCourantNo ( const fvMesh& mesh, const Time& runTime, const volScalarField& rho, const surfaceScalarField& phi ) { scalarField sumPhi ( fvc::surfaceSum(mag(phi))().primitiveField() / rho.primitiveField() ); scalar CoNum = 0.5*gMax(sumPhi/mesh.V().field())*runTime.deltaTValue(); scalar meanCoNum = 0.5*(gSum(sumPhi)/gSum(mesh.V().field()))*runTime.deltaTValue(); Info<< "Region: " << mesh.name() << " Courant Number mean: " << meanCoNum << " max: " << CoNum << endl; return CoNum; }