zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Fixes #416

Browse Source
This commit is contained in:
Prashant
2017-03-02 10:22:24 +05:30
parent 05de53d086
commit 6253452dcf
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
src/functionObjects/solvers/scalarTransport/scalarTransport.C
View File

@ -277,13 +277,13 @@ bool Foam::functionObjects::scalarTransport::execute()
const volScalarField& alpha = const volScalarField& alpha =
mesh_.lookupObject<volScalarField>(phaseName_); mesh_.lookupObject<volScalarField>(phaseName_);
const surfaceScalarField& limitedPhiAlpa = const surfaceScalarField& limitedPhiAlpha =
mesh_.lookupObject<surfaceScalarField>(phasePhiCompressedName_); mesh_.lookupObject<surfaceScalarField>(phasePhiCompressedName_);
D *= pos(alpha - 0.99); D *= pos(alpha - 0.99);
// Reset D dimensions consistent with limitedPhiAlpa // Reset D dimensions consistent with limitedPhiAlpha
D.dimensions().reset(limitedPhiAlpa.dimensions()/dimLength); D.dimensions().reset(limitedPhiAlpha.dimensions()/dimLength);
// Solve // Solve
tmp<surfaceScalarField> tTPhiUD; tmp<surfaceScalarField> tTPhiUD;
@ -292,7 +292,7 @@ bool Foam::functionObjects::scalarTransport::execute()
fvScalarMatrix sEqn fvScalarMatrix sEqn
( (
fvm::ddt(s) fvm::ddt(s)
+ fvm::div(limitedPhiAlpa, s, divScheme) + fvm::div(limitedPhiAlpha, s, divScheme)
- fvm::laplacian(D, s, laplacianScheme) - fvm::laplacian(D, s, laplacianScheme)
== ==
alpha*fvOptions_(s) alpha*fvOptions_(s)