songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

BUG: lagrangian injection model - add to addedMass if particle is immediately absorbed - mantis #1036

Browse Source
This commit is contained in:
andy 2013-10-10 11:02:35 +01:00
parent a787c35c48
commit f285247d1a
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C
View File

@ -550,7 +550,7 @@ void Foam::InjectionModel<CloudType>::inject(TrackData& td)
if (cellI > -1) if (cellI > -1)
{ {
// Lagrangian timestep // Lagrangian timestep
scalar dt = time - timeInj; const scalar dt = time - timeInj;
// Apply corrections to position for 2-D cases // Apply corrections to position for 2-D cases
meshTools::constrainToMeshCentre(mesh, pos); meshTools::constrainToMeshCentre(mesh, pos);
@ -586,8 +586,11 @@ void Foam::InjectionModel<CloudType>::inject(TrackData& td)
pPtr->rho() pPtr->rho()
); );
const scalar mParcel0 = pPtr->nParticle()*pPtr->mass();
if (!pPtr->move(td, dt)) if (!pPtr->move(td, dt))
{ {
massAdded += mParcel0;
delete pPtr; delete pPtr;
} }
else else
@ -595,7 +598,7 @@ void Foam::InjectionModel<CloudType>::inject(TrackData& td)
if (pPtr->nParticle() >= 1.0) if (pPtr->nParticle() >= 1.0)
{ {
td.cloud().addParticle(pPtr); td.cloud().addParticle(pPtr);
massAdded += pPtr->nParticle()*pPtr->mass(); massAdded += mParcel0;
parcelsAdded++; parcelsAdded++;
} }
else else