From e53d830742203c3201109d0349925294df91393c Mon Sep 17 00:00:00 2001 From: Henry Date: Tue, 1 Oct 2013 10:37:16 +0100 Subject: [PATCH 1/2] rhoPimpleDyMFoam: correct dpdt for mesh-motion --- .../solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H index 14ba97d1c2..fd892bb932 100644 --- a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H +++ b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H @@ -115,5 +115,5 @@ K = 0.5*magSqr(U); if (thermo.dpdt()) { - dpdt = fvc::ddt(p); + dpdt = fvc::ddt(p) - fvc::div(fvc::meshPhi(rho, U), p); } From e42af28f7dbca878a6edb6b2ff9bdf7ab6ff1937 Mon Sep 17 00:00:00 2001 From: Henry Date: Tue, 1 Oct 2013 11:29:32 +0100 Subject: [PATCH 2/2] Updated tutorial --- .../rhoPimpleDyMFoam/annularThermalMixer/system/fvSchemes | 1 + 1 file changed, 1 insertion(+) diff --git a/tutorials/compressible/rhoPimpleDyMFoam/annularThermalMixer/system/fvSchemes b/tutorials/compressible/rhoPimpleDyMFoam/annularThermalMixer/system/fvSchemes index ab1cc7da4a..4256201f44 100644 --- a/tutorials/compressible/rhoPimpleDyMFoam/annularThermalMixer/system/fvSchemes +++ b/tutorials/compressible/rhoPimpleDyMFoam/annularThermalMixer/system/fvSchemes @@ -32,6 +32,7 @@ divSchemes div(phi,U) Gauss linearUpwind grad(U); div(phi,h) Gauss linearUpwind grad(h); div(phi,K) Gauss linear; + div(meshPhi,p) Gauss linear; div(phi,k) Gauss upwind; div(phi,epsilon) Gauss upwind; div((muEff*dev2(T(grad(U))))) Gauss linear;