diff --git a/src/functionObjects/solvers/scalarTransport/scalarTransport.C b/src/functionObjects/solvers/scalarTransport/scalarTransport.C index 641241b91e..ebee2bc436 100644 --- a/src/functionObjects/solvers/scalarTransport/scalarTransport.C +++ b/src/functionObjects/solvers/scalarTransport/scalarTransport.C @@ -277,13 +277,13 @@ bool Foam::functionObjects::scalarTransport::execute() const volScalarField& alpha = mesh_.lookupObject(phaseName_); - const surfaceScalarField& limitedPhiAlpa = + const surfaceScalarField& limitedPhiAlpha = mesh_.lookupObject(phasePhiCompressedName_); D *= pos(alpha - 0.99); - // Reset D dimensions consistent with limitedPhiAlpa - D.dimensions().reset(limitedPhiAlpa.dimensions()/dimLength); + // Reset D dimensions consistent with limitedPhiAlpha + D.dimensions().reset(limitedPhiAlpha.dimensions()/dimLength); // Solve tmp tTPhiUD; @@ -292,7 +292,7 @@ bool Foam::functionObjects::scalarTransport::execute() fvScalarMatrix sEqn ( fvm::ddt(s) - + fvm::div(limitedPhiAlpa, s, divScheme) + + fvm::div(limitedPhiAlpha, s, divScheme) - fvm::laplacian(D, s, laplacianScheme) == alpha*fvOptions_(s)