From fd350da163c953b20308925bf5ac2e18bb9f5a30 Mon Sep 17 00:00:00 2001 From: Hrvoje Jasak Date: Wed, 29 Sep 2010 19:58:46 +0100 Subject: [PATCH] Fixed mesh motion error for motion flux on a moving mesh --- src/finiteVolume/finiteVolume/fvc/fvcMeshPhi.C | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/finiteVolume/finiteVolume/fvc/fvcMeshPhi.C b/src/finiteVolume/finiteVolume/fvc/fvcMeshPhi.C index b056fd7bf..65327249e 100644 --- a/src/finiteVolume/finiteVolume/fvc/fvcMeshPhi.C +++ b/src/finiteVolume/finiteVolume/fvc/fvcMeshPhi.C @@ -92,7 +92,7 @@ void Foam::fvc::makeRelative { if (phi.mesh().moving()) { - phi -= fvc::meshPhi(rho, U); + phi -= rho*fvc::meshPhi(rho, U); } } @@ -105,7 +105,7 @@ void Foam::fvc::makeRelative { if (phi.mesh().moving()) { - phi -= fvc::meshPhi(rho, U); + phi -= fvc::interpolate(rho)*fvc::meshPhi(rho, U); } } @@ -131,7 +131,7 @@ void Foam::fvc::makeAbsolute { if (phi.mesh().moving()) { - phi += fvc::meshPhi(rho, U); + phi += rho*fvc::meshPhi(rho, U); } } @@ -144,7 +144,7 @@ void Foam::fvc::makeAbsolute { if (phi.mesh().moving()) { - phi += fvc::meshPhi(rho, U); + phi += fvc::interpolate(rho)*fvc::meshPhi(rho, U); } }