From fffe67d45f03c1d15e8b2f243d35a4064774455f Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Wed, 30 Mar 2011 16:24:35 +0100
Subject: [PATCH] Porosity: Temporarily disable pZones.ddt(rho, U) while the
term is being investigated
---
.../compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/UEqn.H | 3 ++-
applications/solvers/lagrangian/LTSReactingParcelFoam/UEqn.H | 4 ++--
.../lagrangian/porousExplicitSourceReactingParcelFoam/UEqn.H | 4 ++--
.../solvers/multiphase/interFoam/porousInterFoam/UEqn.H | 3 ++-
4 files changed, 8 insertions(+), 6 deletions(-)
diff --git a/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/UEqn.H b/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/UEqn.H
index faae8152ca6..2a9b92d9f0c 100644
--- a/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/UEqn.H
+++ b/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/UEqn.H
@@ -2,7 +2,8 @@
tmp<fvVectorMatrix> UEqn
(
- pZones.ddt(rho, U)
+ //pZones.ddt(rho, U)
+ fvc::ddt(rho, U)
+ fvm::div(phi, U)
+ turbulence->divDevRhoReff(U)
);
diff --git a/applications/solvers/lagrangian/LTSReactingParcelFoam/UEqn.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/UEqn.H
index fee0fe1a689..35394576f88 100644
--- a/applications/solvers/lagrangian/LTSReactingParcelFoam/UEqn.H
+++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/UEqn.H
@@ -1,7 +1,7 @@
fvVectorMatrix UEqn
(
-// fvm::ddt(rho, U)
- pZones.ddt(rho, U)
+ //pZones.ddt(rho, U)
+ fvm::ddt(rho, U)
+ fvm::div(phi, U)
+ turbulence->divDevRhoReff(U)
==
diff --git a/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/UEqn.H b/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/UEqn.H
index fee0fe1a689..35394576f88 100644
--- a/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/UEqn.H
+++ b/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/UEqn.H
@@ -1,7 +1,7 @@
fvVectorMatrix UEqn
(
-// fvm::ddt(rho, U)
- pZones.ddt(rho, U)
+ //pZones.ddt(rho, U)
+ fvm::ddt(rho, U)
+ fvm::div(phi, U)
+ turbulence->divDevRhoReff(U)
==
diff --git a/applications/solvers/multiphase/interFoam/porousInterFoam/UEqn.H b/applications/solvers/multiphase/interFoam/porousInterFoam/UEqn.H
index 9d84a865f52..246a79d3c9d 100644
--- a/applications/solvers/multiphase/interFoam/porousInterFoam/UEqn.H
+++ b/applications/solvers/multiphase/interFoam/porousInterFoam/UEqn.H
@@ -10,7 +10,8 @@
fvVectorMatrix UEqn
(
- pZones.ddt(rho, U)
+ //pZones.ddt(rho, U)
+ fvc::ddt(rho, U)
+ fvm::div(rhoPhi, U)
- fvm::laplacian(muEff, U)
- (fvc::grad(U) & fvc::grad(muEff))
--
GitLab