Commit a8d206f2 authored by sergio's avatar sergio

BUG: Correct div(u*p) for Teqn in compressibleInterFoam .Gitlab 1400.

parent 506990e0
......@@ -4,7 +4,7 @@
fvm::ddt(rho, T) + fvm::div(rhoPhi, T) - fvm::Sp(contErr, T)
- fvm::laplacian(turbulence.alphaEff(), T)
+ (
(divU*p)() - contErr/rho*p
divUp - contErr/rho*p
+ (fvc::ddt(rho, K) + fvc::div(rhoPhi, K))() - contErr*K
)
*(
......
......@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2004-2011, 2017 OpenCFD Ltd.
\\ / A nd | Copyright (C) 2004-2011, 2019 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
| Copyright (C) 2011-2017 OpenFOAM Foundation
......@@ -88,10 +88,11 @@ int main(int argc, char *argv[])
{
#include "readDyMControls.H"
// Store divU from the previous mesh so that it can be mapped
// Store divU and divUp from the previous mesh so that it can be mapped
// and used in correctPhi to ensure the corrected phi has the
// same divergence
volScalarField divU("divU0", fvc::div(fvc::absolute(phi, U)));
volScalarField divUp("divUp", fvc::div(fvc::absolute(phi, U), p));
if (LTS)
{
......
......@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2004-2011, OpenCFD Ltd.
\\ / A nd | Copyright (C) 2004-2019, OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
| Copyright (C) 2011-2017 OpenFOAM Foundation
......@@ -118,7 +118,7 @@ int main(int argc, char *argv[])
turbulence.correctPhasePhi();
#include "UEqn.H"
volScalarField divU(fvc::div(fvc::absolute(phi, U)));
volScalarField divUp("divUp", fvc::div(fvc::absolute(phi, U), p));
#include "TEqn.H"
// --- Pressure corrector loop
......
......@@ -34,6 +34,7 @@ divSchemes
div(rhoPhi,T) Gauss linear;
div(rhoPhi,K) Gauss linear;
div((phi+meshPhi),p) Gauss linear;
div(U,p) Gauss linear;
Please register or sign in to reply
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
......
......@@ -39,6 +39,7 @@ divSchemes
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
div((phi+meshPhi),p) Gauss linear;
div(phi,p) Gauss linear;
div((muEff*dev2(T(grad(U))))) Gauss linear;
}
......
Markdown is supported
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