virtualMassForce::virtualMassForce ( const dictionary& multiphaseSurfaceForcesDict, const multiphase::transport& mtm, const PtrList<volScalarField>& alpha, const volScalarField& beta ) : multiphaseSurfaceForcesDict_(multiphaseSurfaceForcesDict), alpha_(alpha), beta_(beta), mtm_(mtm), virtualMassForce_ ( IOobject ( "virtualMassForce", beta.time().timeName(), beta.mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), beta.mesh(), dimensionedVector("zero", dimensionSet(0, 0, 0, 0, 0), vector(0.0, 0.0, 0.0)) ) {}
Foam::tmp<Foam::volScalarField> Foam::laminarFlameSpeedModels::Gulders::Su0pTphi ( const volScalarField& p, const volScalarField& Tu, scalar phi ) const { tmp<volScalarField> tSu0 ( new volScalarField ( IOobject ( "Su0", p.time().timeName(), p.db(), IOobject::NO_READ, IOobject::NO_WRITE ), p.mesh(), dimensionedScalar("Su0", dimVelocity, 0.0) ) ); volScalarField& Su0 = tSu0(); forAll(Su0, celli) { Su0[celli] = Su0pTphi(p[celli], Tu[celli], phi, 0.0); }
Foam::tmp<Foam::volScalarField> Foam::phaseEquationsOfState::adiabaticPerfectFluid::psi ( const volScalarField& p, const volScalarField& T ) const { return tmp<Foam::volScalarField> ( new volScalarField ( IOobject ( "psi", p.time().timeName(), p.mesh(), IOobject::NO_READ, IOobject::NO_WRITE, false ), (rho0_/(gamma_*(p0_ + B_))) *pow((p + B_)/(p0_ + B_), 1.0/gamma_ - 1.0) ) ); }
Foam::tmp<Foam::volScalarField> Foam::consumptionSpeed::omega0Sigma ( const volScalarField& sigma ) { tmp<volScalarField> tomega0 ( new volScalarField ( IOobject ( "omega0", sigma.time().timeName(), sigma.db(), IOobject::NO_READ, IOobject::NO_WRITE ), sigma.mesh(), dimensionedScalar ( "omega0", dimensionSet(1, -2, -1, 0, 0, 0, 0), 0 ) ) ); volScalarField& omega0 = tomega0(); volScalarField::InternalField& iomega0 = omega0.internalField(); forAll(iomega0, celli) { iomega0[celli] = omega0Sigma(sigma[celli], 1.0); }
void writeCellGraph ( const volScalarField& vsf, const word& graphFormat ) { fileName path(vsf.time().path()/"graphs"/vsf.time().timeName()); mkDir(path); graph ( vsf.name(), "x", vsf.name(), vsf.mesh().C().primitiveField().component(vector::X), vsf.primitiveField() ).write(path/vsf.name(), graphFormat); }
thermalModel::thermalModel(const volScalarField& T) : IOdictionary ( IOobject ( "thermalProperties", T.time().constant(), T.db(), IOobject::MUST_READ, IOobject::NO_WRITE ) ), T_(T), lawPtr_(thermalLaw::New("law", T_, subDict("thermal"))) {}
Foam::autoPtr<Foam::radiation::radiationModel> Foam::radiation::radiationModel::New ( const volScalarField& T ) { IOobject radIO ( "radiationProperties", T.time().constant(), T.mesh(), IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE, false ); word modelType("none"); if (radIO.typeHeaderOk<IOdictionary>(false)) { IOdictionary(radIO).lookup("radiationModel") >> modelType; }
Foam::tmp<Foam::volScalarField> Foam::phaseEquationsOfState::linear::rho ( const volScalarField& p, const volScalarField& T ) const { return tmp<Foam::volScalarField> ( new volScalarField ( IOobject ( "rho", p.time().timeName(), p.mesh(), IOobject::NO_READ, IOobject::NO_WRITE, false ), rho0_ + psi_*p ) ); }
Foam::XiModel::XiModel ( const dictionary& XiProperties, const hhuCombustionThermo& thermo, const compressible::RASModel& turbulence, const volScalarField& Su, const volScalarField& rho, const volScalarField& b, const surfaceScalarField& phi ) : XiModelCoeffs_ ( XiProperties.subDict ( word(XiProperties.lookup("XiModel")) + "Coeffs" ) ), thermo_(thermo), turbulence_(turbulence), Su_(Su), rho_(rho), b_(b), phi_(phi), Xi_ ( IOobject ( "Xi", b.time().timeName(), b.db(), IOobject::MUST_READ, IOobject::AUTO_WRITE ), b.mesh() ) {}