zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

updates

Browse Source
This commit is contained in:
andy
2009-02-20 19:12:27 +00:00
parent e2bad4ddac
commit 761ccd2bb3
2 changed files with 32 additions and 24 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -127,10 +127,10 @@ template<class CloudType>
Foam::scalar Foam::InjectionModel<CloudType>::setNumberOfParticles
(
const label parcels,
const scalar diameter,
const scalar volumeFraction,
const scalar rho,
const scalar volume
const scalar volumeFraction,
const scalar diameter,
const scalar rho,
)
{
scalar nP = 0.0;
@ -294,30 +294,38 @@ void Foam::InjectionModel<CloudType>::inject(TrackData& td)
vector pos = vector::zero;
setPositionAndCell(iParcel, timeInj, owner_.meshInfo(), pos, cellI);
// Diameter of parcels
scalar d = d0(iParcel, timeInj);
// Number of particles per parcel
scalar nP = setNumberOfParticles
(
newParcels,
d,
volFraction,
rho,
newVolume
);
// Velocity of parcels
vector U = velocity(iParcel, timeInj, owner_.meshInfo());
if (cellI >= 0)
{
scalar dt = time - timeInj;
td.cloud().addNewParcel(pos, cellI, d, U, nP, dt);
// Diameter of parcels
scalar d = d0(iParcel, timeInj);
// Number of particles per parcel
scalar nP = setNumberOfParticles
(
newParcels,
newVolume,
volFraction,
d,
rho
);
// Velocity of parcels
vector U = velocity(iParcel, timeInj, owner_.meshInfo());
// Lagrangian timestep
scalar dt = time - timeInj;
// Add the new parcel
td.cloud().addNewParcel(pos, cellI, d, U, nP, dt);
massInjected_ += nP*rho*mathematicalConstant::pi*pow3(d)/6.0;
parcelsAdded_++;
}
else
{
WarningIn("Foam::InjectionModel<CloudType>::inject(TrackData& td)")
<< "Failed to inject new parcel:" <<
<< " id = " << iParcel << ", position = " pos << nl << endl;
}
}
postInjectCheck();

6
src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.H
View File

@ -170,10 +170,10 @@ protected:
scalar setNumberOfParticles
(
const label parcels,
const scalar diameter,
const scalar volume,
const scalar volumeFraction,
const scalar rho,
const scalar volume
const scalar diameter,
const scalar rho
);
//- Post injection checks