Commit 6247f6a5 authored by Henry's avatar Henry
Browse files

Correct child parcel properties

Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1287
parent 38b60d68
...@@ -253,7 +253,7 @@ void Foam::SprayParcel<ParcelType>::calcBreakup ...@@ -253,7 +253,7 @@ void Foam::SprayParcel<ParcelType>::calcBreakup
const vector g = td.cloud().g().value(); const vector g = td.cloud().g().value();
scalar massChild = 0.0; scalar parcelMassChild = 0.0;
scalar dChild = 0.0; scalar dChild = 0.0;
if if
( (
...@@ -279,24 +279,26 @@ void Foam::SprayParcel<ParcelType>::calcBreakup ...@@ -279,24 +279,26 @@ void Foam::SprayParcel<ParcelType>::calcBreakup
Urmag, Urmag,
this->tMom(), this->tMom(),
dChild, dChild,
massChild parcelMassChild
) )
) )
{ {
scalar Re = rhoAv*Urmag*dChild/muAv; scalar Re = rhoAv*Urmag*dChild/muAv;
this->mass0() -= massChild;
// Add child parcel as copy of parent // Add child parcel as copy of parent
SprayParcel<ParcelType>* child = new SprayParcel<ParcelType>(*this); SprayParcel<ParcelType>* child = new SprayParcel<ParcelType>(*this);
child->mass0() = massChild;
child->d() = dChild; 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 = const forceSuSp Fcp =
forces.calcCoupled(*child, dt, massChild, Re, muAv); forces.calcCoupled(*child, dt, massChild, Re, muAv);
const forceSuSp Fncp = const forceSuSp Fncp =
forces.calcNonCoupled(*child, dt, massChild, Re, muAv); forces.calcNonCoupled(*child, dt, massChild, Re, muAv);
child->age() = 0.0;
child->liquidCore() = 0.0; child->liquidCore() = 0.0;
child->KHindex() = 1.0; child->KHindex() = 1.0;
child->y() = td.cloud().breakup().y0(); child->y() = td.cloud().breakup().y0();
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment