CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Lagrangian Particle, first particle integration step wrong dt and age?

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   February 13, 2020, 06:39
Default Lagrangian Particle, first particle integration step wrong dt and age?
New Member
Join Date: Feb 2020
Posts: 1
Rep Power: 0
Sverre Magnus is on a distinguished road

I am using reactingParcelFoam as LTS (local -time-stepping) for a steady state analysis. Using a precalculated flow result (reactingCloud1Properties -> active no, I do one steady particle injection (reactingCloud1Properties -> active yes; transient no. I use v1812 at the moment.

1. First thing I noticed is that deltaTMax from reactingCloud1Propetries is not used in algorithm.

2. In the particle integration loop in "Foam::KinematicParcel<ParcelType>::move" in KinematicParcel.C in the first step the particle velocity is not yet calculated so U_ is (0 0 0), calcVelocity was not yet called. The displacement is reported as:
Tracking from (0.02 0.005 0.003) along (0 0 0) to (0.02 0.005 0.003) which makes sense when particle velocity is zero. But what is somehow strange is that the stepfraction is increased as a result of factor f has size of maxCo. Threafter the time-step dt is calculated to a non-zero value due to non-zero stepfraction. Time-step is added to age, so the particle has a starting age which is proportional to maxCo and trackTime (trackTime is the maxTrackTime from reactingCloud1Properties).
First step of while loop:
f = 1- p.stepFraction (latter is zero at start) so f = 1
then f = min(f,maxCo), with maxCo of 0.01 f is 0.01.
The trackToFace (actually in particle::trackToStationaryTri) algorithm sets stepFraction to 0.01. So with this stepfraction, sfrac (is stepFraction at start=0) and trackTime of 10s dt is 0.1s and age increases to 0.1s without any movement. This is also visible in the age file and one can visualize that in paraview. The particle velocity is calculated afterwards with (I guess) the new particle position and flow conditions and is used in the next tracking step, which seams to be ok.

That would mean in my opinion that the particle age is incorrect. Did anybody observe the same behaviour?
Sverre Magnus is offline   Reply With Quote


lagrangian method, particle cloud

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

All times are GMT -4. The time now is 08:32.