43 const scalar trackTime
54 <<
" trackTime = " << trackTime
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;
Base cloud calls templated on particle type.
static word timeName(const scalar t, const int precision=precision_)
A cloud is a registry collection of lagrangian particles.
Type interpolate(const cellPointWeight &cpw) const
Interpolate field for the given cellPointWeight.
const Time & time() const noexcept
Return time registry.
bool switchProcessor
Flag to switch processor.
bool keepParticle
Flag to indicate whether to keep particle (false = delete)
virtual void transformProperties(const tensor &T)
Transform the physical properties of the particle.
tetIndices currentTetIndices() const noexcept
Return indices of the current tet that the particle occupies.
const polyMesh & mesh() const noexcept
Return the mesh database.
const barycentric & coordinates() const noexcept
Return current particle coordinates.
scalar stepFraction() const noexcept
Return the fraction of time-step completed.
void trackToAndHitFace(const vector &direction, const scalar fraction, TrackCloudType &cloud, trackingData &td)
Convenience function. Combines trackToFace and hitFace.
A Cloud of solid particles.
Class used to pass tracking data to the trackToFace function.
const interpolationCellPoint< scalar > & rhoInterp() const
const interpolationCellPoint< vector > & UInterp() const
const interpolationCellPoint< scalar > & nuInterp() const
virtual void transformProperties(const tensor &T)
Transform the physical properties of the particle.
void hitWallPatch(solidParticleCloud &cloud, trackingData &td)
Overridable function to handle the particle hitting a wallPatch.
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.
Storage and named access for the indices of a tet which is part of the decomposition of a cell.
#define defineTemplateTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information for templates, useful.
scalarField Re(const UList< complex > &cf)
Extract real component.
dimensionSet transform(const dimensionSet &ds)
Return the argument; transformations do not change the dimensions.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Ostream & endl(Ostream &os)
Add newline and flush stream.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
prefixOSstream Pout
OSstream wrapped stdout (std::cout) with parallel prefix.