From 6247f6a5d46c2d27dcaadb48bbb905e6c236a7b5 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Thu, 18 Dec 2014 13:07:17 +0000
Subject: [PATCH] Correct child parcel properties Resolves bug-report
http://www.openfoam.org/mantisbt/view.php?id=1287
---
.../parcels/Templates/SprayParcel/SprayParcel.C | 12 +++++++-----
1 file changed, 7 insertions(+), 5 deletions(-)
diff --git a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
index 062c4e45a07..991d25ae1ca 100644
--- a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
+++ b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
@@ -253,7 +253,7 @@ void Foam::SprayParcel<ParcelType>::calcBreakup
const vector g = td.cloud().g().value();
- scalar massChild = 0.0;
+ scalar parcelMassChild = 0.0;
scalar dChild = 0.0;
if
(
@@ -279,24 +279,26 @@ void Foam::SprayParcel<ParcelType>::calcBreakup
Urmag,
this->tMom(),
dChild,
- massChild
+ parcelMassChild
)
)
{
scalar Re = rhoAv*Urmag*dChild/muAv;
- this->mass0() -= massChild;
// Add child parcel as copy of parent
SprayParcel<ParcelType>* child = new SprayParcel<ParcelType>(*this);
- child->mass0() = massChild;
child->d() = dChild;
- child->nParticle() = massChild/(this->rho()*this->volume(dChild));
+ child->d0() = dChild;
+ const scalar massChild = child->mass();
+ child->mass0() = massChild;
+ child->nParticle() = parcelMassChild/massChild;
const forceSuSp Fcp =
forces.calcCoupled(*child, dt, massChild, Re, muAv);
const forceSuSp Fncp =
forces.calcNonCoupled(*child, dt, massChild, Re, muAv);
+ child->age() = 0.0;
child->liquidCore() = 0.0;
child->KHindex() = 1.0;
child->y() = td.cloud().breakup().y0();
--
GitLab