From 33e452f4d3a56235df0c7c111a27176a2ae706c8 Mon Sep 17 00:00:00 2001 From: andy Date: Wed, 3 Jun 2009 19:40:33 +0100 Subject: [PATCH] inclusion of point mass sources --- .../solvers/lagrangian/trackedReactingParcelFoam/pEqn.H | 2 ++ .../solvers/lagrangian/trackedReactingParcelFoam/rhoEqn.H | 1 + 2 files changed, 3 insertions(+) diff --git a/applications/solvers/lagrangian/trackedReactingParcelFoam/pEqn.H b/applications/solvers/lagrangian/trackedReactingParcelFoam/pEqn.H index b26a6594aa..b6207ce0c4 100644 --- a/applications/solvers/lagrangian/trackedReactingParcelFoam/pEqn.H +++ b/applications/solvers/lagrangian/trackedReactingParcelFoam/pEqn.H @@ -41,6 +41,7 @@ if (transonic) - lapTerm() == reactingParcels.Srho() + + pointMassSources.Su() ); pEqn.solve(); @@ -80,6 +81,7 @@ else - lapTerm() == reactingParcels.Srho() + + pointMassSources.Su() ); pEqn.solve(); diff --git a/applications/solvers/lagrangian/trackedReactingParcelFoam/rhoEqn.H b/applications/solvers/lagrangian/trackedReactingParcelFoam/rhoEqn.H index 03eed2dc82..9517d2d912 100644 --- a/applications/solvers/lagrangian/trackedReactingParcelFoam/rhoEqn.H +++ b/applications/solvers/lagrangian/trackedReactingParcelFoam/rhoEqn.H @@ -37,6 +37,7 @@ Description + fvc::div(phi) == reactingParcels.Srho() + + pointMassSources.Su() ); }