From 0ca9ddbd49c21f28ce3ccf79aaddc0fd7c7f1e3f Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Thu, 28 Jun 2012 16:24:22 +0100
Subject: [PATCH] BUG: spray parcel - corrected calculation of tMom - mantis
#565
---
.../spray/parcels/Templates/SprayParcel/SprayParcel.C | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
index 8e291b75b71..3381bbc087f 100644
--- a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
+++ b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -246,7 +246,7 @@ void Foam::SprayParcel<ParcelType>::calcBreakup
const scalar mass = p.mass();
const forceSuSp Fcp = forces.calcCoupled(p, dt, mass, Re, muAv);
const forceSuSp Fncp = forces.calcNonCoupled(p, dt, mass, Re, muAv);
- scalar tMom = 1.0/(Fcp.Sp() + Fncp.Sp());
+ scalar tMom = mass/(Fcp.Sp() + Fncp.Sp());
const vector g = td.cloud().g().value();
@@ -301,7 +301,7 @@ void Foam::SprayParcel<ParcelType>::calcBreakup
child->tc() = -GREAT;
child->ms() = 0.0;
child->injector() = this->injector();
- child->tMom() = 1.0/(Fcp.Sp() + Fncp.Sp());
+ child->tMom() = massChild/(Fcp.Sp() + Fncp.Sp());
child->user() = 0.0;
child->setCellValues(td, dt, cellI);
--
GitLab