inclusion of point mass sources

2025-11-28 03:28:01 +00:00 · 2009-06-03 19:40:33 +01:00
parent 2bb2266872
commit 33e452f4d3
2 changed files with 3 additions and 0 deletions
--- 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();
--- 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()
    );
 }