diff --git a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.C index 9f3394c28d..3773321506 100644 --- a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.C +++ b/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.C @@ -241,6 +241,8 @@ bool Foam::InteractingKinematicParcel::move(TrackData& td) p.U() += 0.5*deltaT*p.f()/p.mass(); + pi_ += 0.5*deltaT*tau_; + break; } diff --git a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.H b/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.H index 475c2db163..48ae97ed79 100644 --- a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.H +++ b/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.H @@ -34,6 +34,9 @@ Description - turbulent dispersion - wall interactions + Includes many-body collisions and rotational motion. Only + spherical particles currently supported. + SourceFiles InteractingKinematicParcelI.H InteractingKinematicParcel.C @@ -432,6 +435,9 @@ public: //- Particle moment of inertia around diameter axis inline scalar momentOfInertia() const; + //- Particle angular velocity + inline vector omega() const; + //- Particle volume inline scalar volume() const; diff --git a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcelI.H b/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcelI.H index 0f95bceb25..ee2df1df80 100644 --- a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcelI.H +++ b/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcelI.H @@ -437,6 +437,14 @@ Foam::InteractingKinematicParcel::momentOfInertia() const } +template +inline Foam::vector +Foam::InteractingKinematicParcel::omega() const +{ + return pi_/momentOfInertia(); +} + + template inline Foam::scalar Foam::InteractingKinematicParcel::volume() const {