Skip to content
Snippets Groups Projects
Commit b3bfdd0b authored by Andrew Heather's avatar Andrew Heather
Browse files

Merge branch 'Fix_GL_416' into 'master'

Fixes #416

Closes #416

See merge request !96
parents 05de53d0 6253452d
Branches
Tags
No related merge requests found
......@@ -277,13 +277,13 @@ bool Foam::functionObjects::scalarTransport::execute()
const volScalarField& alpha =
mesh_.lookupObject<volScalarField>(phaseName_);
const surfaceScalarField& limitedPhiAlpa =
const surfaceScalarField& limitedPhiAlpha =
mesh_.lookupObject<surfaceScalarField>(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<surfaceScalarField> 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)
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment