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); }
forAll(Su0Bf[patchi], facei) { Su0Bf[patchi][facei] = Su0pTphi ( p.boundaryField()[patchi][facei], Tu.boundaryField()[patchi][facei], phi, 0.0 ); }
volScalarField& Su0 = tSu0(); forAll(Su0, celli) { Su0[celli] = Su0pTphi(p[celli], Tu[celli], phi, 0.0); } forAll(Su0.boundaryField(), patchi) { forAll(Su0.boundaryField()[patchi], facei) { Su0.boundaryField()[patchi][facei] = Su0pTphi ( p.boundaryField()[patchi][facei], Tu.boundaryField()[patchi][facei], phi, 0.0 ); } } return tSu0; } Foam::tmp<Foam::volScalarField> Foam::laminarFlameSpeedModels::Gulders::Su0pTphi ( const volScalarField& p, const volScalarField& Tu, const volScalarField& phi
forAll(Su0p, facei) { Su0p[facei] = Su0pTphi(pp[facei], Tup[facei], phi); }