diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffNo.C b/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffNo.C index a93f19e199..bb395003a3 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffNo.C +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffNo.C @@ -68,14 +68,16 @@ Foam::scalar Foam::solidRegionDiffNo scalar DiNum = 0.0; scalar meanDiNum = 0.0; - volScalarField K = mag(Kdirectional); + volScalarField K(mag(Kdirectional)); //- Take care: can have fluid domains with 0 cells so do not test for // zero internal faces. - surfaceScalarField KrhoCpbyDelta = + surfaceScalarField KrhoCpbyDelta + ( mesh.surfaceInterpolation::deltaCoeffs() * fvc::interpolate(K) - / fvc::interpolate(Cprho); + / fvc::interpolate(Cprho) + ); DiNum = gMax(KrhoCpbyDelta.internalField())*runTime.deltaT().value(); diff --git a/applications/solvers/stressAnalysis/solidEquilibriumDisplacementFoam/kineticEnergyLimiter.H b/applications/solvers/stressAnalysis/solidEquilibriumDisplacementFoam/kineticEnergyLimiter.H index eae3c3d327..169c7f229f 100644 --- a/applications/solvers/stressAnalysis/solidEquilibriumDisplacementFoam/kineticEnergyLimiter.H +++ b/applications/solvers/stressAnalysis/solidEquilibriumDisplacementFoam/kineticEnergyLimiter.H @@ -1,6 +1,6 @@ if (!(runTime.timeIndex() % 5)) { - volScalarField kineticEnergy(magSqr(Dcorr); + volScalarField kineticEnergy(magSqr(Dcorr)); dimensionedScalar intKineticEnergy = fvc::domainIntegrate(kineticEnergy); Info<< "kineticEnergy = " << intKineticEnergy.value() << endl;