Go to the documentation of this file.
43 const scalar trackTime
53 Pout<<
"Time = " <<
mesh().time().timeName()
54 <<
" trackTime = " << trackTime
55 <<
" steptFraction() = " << stepFraction() <<
endl;
59 const scalar sfrac = stepFraction();
61 const scalar
f = 1 - stepFraction();
62 trackToAndHitFace(
f*trackTime*U_,
f,
cloud, td);
64 const scalar dt = (stepFraction() - sfrac)*trackTime;
66 const tetIndices tetIs = this->currentTetIndices();
71 scalar rhop =
cloud.rhop();
72 scalar magUr =
mag(Uc - U_);
75 scalar
Re = magUr*d_/nuc;
79 ReFunc += 0.15*
pow(
Re, 0.687);
82 scalar Dc = (24.0*nuc/d_)*ReFunc*(3.0/4.0)*(rhoc/(d_*rhop));
84 U_ = (U_ + dt*(Dc*Uc + (1.0 - rhoc/rhop)*td.
g()))/(1.0 + dt*Dc);
109 const vector nw = normal();
116 U_ -= (1.0 +
cloud.e())*Un*nw;
int debug
Static debugging option.
A Cloud of solid particles.
void hitWallPatch(solidParticleCloud &cloud, trackingData &td)
Overridable function to handle the particle hitting a wallPatch.
defineTemplateTypeNameAndDebug(faScalarMatrix, 0)
Ostream & endl(Ostream &os)
Add newline and flush stream.
prefixOSstream Pout
An Ostream wrapper for parallel output to std::cout.
dimensionSet transform(const dimensionSet &ds)
Return the argument; transformations do not change the dimensions.
virtual void transformProperties(const tensor &T)
Transform the physical properties of the particle.
const interpolationCellPoint< vector > & UInterp() const
const interpolationCellPoint< scalar > & nuInterp() const
PtrList< coordinateSystem > coordinates(solidRegions.size())
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
bool keepParticle
Flag to indicate whether to keep particle (false = delete)
bool switchProcessor
Flag to switch processor.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Type interpolate(const cellPointWeight &cpw) const
Interpolate field for the given cellPointWeight.
Class used to pass tracking data to the trackToFace function.
Storage and named access for the indices of a tet which is part of the decomposition of a cell.
A cloud is a registry collection of lagrangian particles.
virtual void transformProperties(const tensor &T)
Transform the physical properties of the particle.
bool move(solidParticleCloud &, trackingData &, const scalar)
Move.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
scalarField Re(const UList< complex > &cf)
Extract real component.
const interpolationCellPoint< scalar > & rhoInterp() const
bool hitPatch(solidParticleCloud &cloud, trackingData &td)
Overridable function to handle the particle hitting a patch.
void hitProcessorPatch(solidParticleCloud &cloud, trackingData &td)
Overridable function to handle the particle hitting a.