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

Lagrangian particle tracking with breakup/coalescence

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

Like Tree2Likes
  • 1 Post By A_Pete
  • 1 Post By

LinkBack Thread Tools Search this Thread Display Modes
Old   May 30, 2016, 05:16
Default Lagrangian particle tracking with breakup/coalescence
Andreas Weber
Join Date: Jun 2014
Posts: 37
Rep Power: 11 is on a distinguished road
Hi everyone,

lagrangian particles can be tracked using "particleTracks" in postprocessing, so far so good. But when taking a closer look at the code, i found that the tracks are produced with the "origId" and "origProc" of each parcel. This works allright for simple particle simulation where parcels cannot breakup or coalesce. The "origId" is assigned to the parcel on injection and will not change anymore, so it can be easily tracked.
But when you have breakup in you simulation, it gets messed up because the daughter parcels will be created as a copy of the motherparcel. The daughterparcel will therefore get the same "origId", so it will be possible to have different parcels with the same origId in your domain. Particle tracking will fail. (It still works, but the tracks are wrong).

All in all... parcels do not have a uniqe ID that could be used for tracking in such a complicated simulation. I was wondering how that could be achieved, but since you have to count all injected parcels + added daughter parcels from breakup it could get all messed up.

So 1st question: How can i create a counter, that is accessible from injectionModel and breakupModel?
2nd question: would that also work with parallel runs?

Best regards
Andreas is offline   Reply With Quote

Old   May 30, 2016, 10:59
Join Date: Jul 2011
Posts: 54
Rep Power: 14
A_Pete is on a distinguished road
Hi Andreas,

have you had a look at the "particleCount_" variable as well as the "getNewParticleID()" function in particle.H? It says particleCount_ may be used to provide a unique ID. And getNewParticleID() should generate a new unique particle creation ID.

How does the break-up process work? Can't you manipulate the code in a way that you copy the child particle from the mother object but giving it a new ID?

Best regards,
granzer likes this.
A_Pete is offline   Reply With Quote

Old   May 31, 2016, 07:22
Andreas Weber
Join Date: Jun 2014
Posts: 37
Rep Power: 11 is on a distinguished road

well, wow, that was surprisingly easy ;-)
Have not seen that getNewParticleID() before.

I have simply added this into the breakup part in SprayParcel.C:
child->origId() = this->getNewParticleID();
which gives every new daughterparcel a nice new origId()...
Evey parcel will have a uniqe ID now!

Great thanks to you!

BTW: This will also enable to change coalescence (stochastic collision) models to work more efficient. Working on a code, which will check particle collision with respect to neighbouring cells as well.
granzer likes this. is offline   Reply With Quote

Old   October 3, 2023, 08:55
New Member
Join Date: Oct 2022
Posts: 7
Rep Power: 3
erice is on a distinguished road
Hi Andreas,

How to use the coalesce model for lagrangian particlle.
I do not find it in submodels for lagrangian.
erice is offline   Reply With Quote


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

Similar Threads
Thread Thread Starter Forum Replies Last Post
[swak4Foam] Validation of swak4Foam for lagrangian particle tracking nero235 OpenFOAM Community Contributions 13 March 28, 2017 05:45
Request for Lagrangian Particle Tracking Validation or Verification Paper Mojtaba.a OpenFOAM Verification & Validation 6 May 23, 2016 01:47
Lagrangian Particle Tracking in Eulerian-Eulerian Multiphase Flow DarrenC CFX 5 April 7, 2016 14:50
Ubuntu 12.10 + openfoam2.2.0 ==> paraview error message peteryuan OpenFOAM Installation 6 August 18, 2013 18:00
injection problem Mark New FLUENT 0 August 4, 2013 01:30

All times are GMT -4. The time now is 02:42.