diff --git a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/turbulentTemperatureCoupledBaffle/turbulentTemperatureCoupledBaffleFvPatchScalarField.C b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/turbulentTemperatureCoupledBaffle/turbulentTemperatureCoupledBaffleFvPatchScalarField.C index 2d3981fdde590b28396716f44f9c5fdc195072d2..05a26301112735fc5ed26ef0ad0d53f3691e08c4 100644 --- a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/turbulentTemperatureCoupledBaffle/turbulentTemperatureCoupledBaffleFvPatchScalarField.C +++ b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/turbulentTemperatureCoupledBaffle/turbulentTemperatureCoupledBaffleFvPatchScalarField.C @@ -215,6 +215,8 @@ turbulentTemperatureCoupledBaffleFvPatchScalarField tmp<scalarField> turbulentTemperatureCoupledBaffleFvPatchScalarField::K() const { + const fvMesh& mesh = patch().boundaryMesh().mesh(); + if (KName_ == "none") { const compressible::RASModel& model = @@ -230,10 +232,32 @@ turbulentTemperatureCoupledBaffleFvPatchScalarField::K() const *thermo.rho()().boundaryField()[patch().index()] *thermo.Cp()().boundaryField()[patch().index()]; } - else + else if (mesh.objectRegistry::foundObject<volScalarField>(KName_)) { return patch().lookupPatchField<volScalarField, scalar>(KName_); } + else if (mesh.objectRegistry::foundObject<volSymmTensorField>(KName_)) + { + const symmTensorField& KWall = + patch().lookupPatchField<volSymmTensorField, scalar>(KName_); + + vectorField n = patch().nf(); + + return n & KWall & n; + } + else + { + FatalErrorIn + ( + "turbulentTemperatureCoupledBaffleFvPatchScalarField::K() const" + ) << "Did not find field " << KName_ + << " on mesh " << mesh.name() << " patch " << patch().name() + << endl + << "Please set 'K' to 'none', a valid volScalarField" + << " or a valid volSymmTensorField." << exit(FatalError); + + return scalarField(0); + } }