Adapt turbulent BCs to postFixedObjReg

This commit is contained in:
Henrik Rusche 2013-10-10 00:32:41 +02:00
parent ffe78c66d8
commit 9f8cf26779
6 changed files with 12 additions and 12 deletions

View file

@ -140,10 +140,10 @@ updateCoeffs()
}
const fvPatchVectorField& Up =
patch().lookupPatchField<volVectorField, vector>(UName_);
lookupPatchField<volVectorField, vector>(UName_);
const fvsPatchScalarField& phip =
patch().lookupPatchField<surfaceScalarField, scalar>(phiName_);
lookupPatchField<surfaceScalarField, scalar>(phiName_);
this->refValue() = SMALL + 1.5*sqr(intensity_)*magSqr(Up);
this->valueFraction() = neg(phip);

View file

@ -136,10 +136,10 @@ void turbulentMixingLengthDissipationRateInletFvPatchScalarField::updateCoeffs()
const scalar Cmu75 = pow(Cmu, 0.75);
const fvPatchScalarField& kp =
patch().lookupPatchField<volScalarField, scalar>(kName_);
lookupPatchField<volScalarField, scalar>(kName_);
const fvsPatchScalarField& phip =
patch().lookupPatchField<surfaceScalarField, scalar>(phiName_);
lookupPatchField<surfaceScalarField, scalar>(phiName_);
this->refValue() = Cmu75*kp*sqrt(kp)/mixingLength_;
this->valueFraction() = neg(phip);

View file

@ -136,10 +136,10 @@ void turbulentMixingLengthFrequencyInletFvPatchScalarField::updateCoeffs()
const scalar Cmu25 = pow(Cmu, 0.25);
const fvPatchScalarField& kp =
patch().lookupPatchField<volScalarField, scalar>(kName_);
lookupPatchField<volScalarField, scalar>(kName_);
const fvsPatchScalarField& phip =
patch().lookupPatchField<surfaceScalarField, scalar>(phiName_);
lookupPatchField<surfaceScalarField, scalar>(phiName_);
this->refValue() = sqrt(kp)/(Cmu25*mixingLength_);
this->valueFraction() = neg(phip);

View file

@ -136,10 +136,10 @@ void turbulentMixingLengthDissipationRateInletFvPatchScalarField::updateCoeffs()
const scalar Cmu75 = pow(Cmu, 0.75);
const fvPatchScalarField& kp =
patch().lookupPatchField<volScalarField, scalar>(kName_);
lookupPatchField<volScalarField, scalar>(kName_);
const fvsPatchScalarField& phip =
patch().lookupPatchField<surfaceScalarField, scalar>(phiName_);
lookupPatchField<surfaceScalarField, scalar>(phiName_);
this->refValue() = Cmu75*kp*sqrt(kp)/mixingLength_;
this->valueFraction() = neg(phip);

View file

@ -136,10 +136,10 @@ void turbulentMixingLengthFrequencyInletFvPatchScalarField::updateCoeffs()
const scalar Cmu25 = pow(Cmu, 0.25);
const fvPatchScalarField& kp =
patch().lookupPatchField<volScalarField, scalar>(kName_);
lookupPatchField<volScalarField, scalar>(kName_);
const fvsPatchScalarField& phip =
patch().lookupPatchField<surfaceScalarField, scalar>(phiName_);
lookupPatchField<surfaceScalarField, scalar>(phiName_);
this->refValue() = sqrt(kp)/(Cmu25*mixingLength_);
this->valueFraction() = neg(phip);

View file

@ -158,10 +158,10 @@ void RWallFunctionFvPatchSymmTensorField::evaluate
)
{
const scalarField& nutw =
patch().lookupPatchField<volScalarField, scalar>(nutName_);
lookupPatchField<volScalarField, scalar>(nutName_);
const fvPatchVectorField& Uw =
patch().lookupPatchField<volVectorField, vector>(UName_);
lookupPatchField<volVectorField, vector>(UName_);
// Old treatment: anisotropy not enforced
vectorField snGradU = Uw.snGrad();