a884a44e47
Patch contributed by Tobias Holzmann Resolves patch request https://bugs.openfoam.org/view.php?id=2524
35 lines
901 B
C
35 lines
901 B
C
fvMesh& mesh = solidRegions[i];
|
|
solidThermo& thermo = thermos[i];
|
|
|
|
tmp<volScalarField> trho = thermo.rho();
|
|
const volScalarField& rho = trho();
|
|
|
|
tmp<volScalarField> tcp = thermo.Cp();
|
|
const volScalarField& cp = tcp();
|
|
|
|
tmp<volSymmTensorField> taniAlpha;
|
|
if (!thermo.isotropic())
|
|
{
|
|
volSymmTensorField& aniAlpha = aniAlphas[i];
|
|
tmp<volVectorField> tkappaByCp = thermo.Kappa()/cp;
|
|
const coordinateSystem& coodSys = coordinates[i];
|
|
|
|
aniAlpha.primitiveFieldRef() =
|
|
coodSys.R().transformVector(tkappaByCp());
|
|
aniAlpha.correctBoundaryConditions();
|
|
|
|
taniAlpha = tmp<volSymmTensorField>
|
|
(
|
|
new volSymmTensorField(aniAlpha)
|
|
);
|
|
}
|
|
|
|
|
|
volScalarField& h = thermo.he();
|
|
|
|
const volScalarField& betav = betavSolid[i];
|
|
|
|
fv::options& fvOptions = solidHeatSources[i];
|
|
|
|
bool& resReachedSolid = residualReachedSolid[i];
|
|
bool& residualControlUsed = residualControlUsedSolid[i];
|