BUG: kinematicParcelFoam: correct source-term value (fixes #3446)

- likely results in output changes for denser carrier fluids Co-authored-by: Kutalmis Bercin <kutalmis.bercin@esi-group.com>
2025-11-28 03:28:01 +00:00 · 2025-10-03 15:29:52 +01:00
parent 7a266b566f
commit 9c4d0cdeef
2 changed files with 6 additions and 1 deletions
--- a/applications/solvers/lagrangian/kinematicParcelFoam/UEqn.H
+++ b/applications/solvers/lagrangian/kinematicParcelFoam/UEqn.H
@ -6,7 +6,7 @@
      + MRF.DDt(U)
      + turbulence->divDevReff(U)
     ==
-        parcels.SU(U, true)
+        invRhoInf*parcels.SU(U)
      + fvOptions(U)
    );

--- a/applications/solvers/lagrangian/kinematicParcelFoam/createFields.H
+++ b/applications/solvers/lagrangian/kinematicParcelFoam/createFields.H
@ -39,6 +39,11 @@ dimensionedScalar rhoInfValue
    laminarTransport
 );

+const dimensionedScalar invRhoInf
+(
+    dimless/dimDensity, scalar(1)/rhoInfValue.value()
+);
+
 volScalarField rhoInf
 (
    IOobject