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

VoF-Lagrangian Particle Tracking interface issue

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

LinkBack Thread Tools Search this Thread Display Modes
Old   June 11, 2021, 11:20
Default VoF-Lagrangian Particle Tracking interface issue
New Member
Join Date: Jun 2021
Posts: 2
Rep Power: 0
Danny_ is on a distinguished road
I'm trying to simulate a multiphase system represented by a hexaedral vessel partially filled with water which is later to be charged by solid particles injected from the upper part. The remaining part of the vessel is obviously filled by air at t=0.
To do that, I coupled the interFoam native solver with the src/lagrangian/intermediate directory, as described in "" and managed to correct all the errors found.
The simulation runs fine; however, the particles injected once fallen below because of gravity, stop at the water-air initial horizontal interface and never cross it. Therefore, I am not able to simulate the mixing process I am interested at for this system.
So far I changed repeatedly different parameters I supposed relevant (U_O of particles injection, rho.air, rho.water, nu.air, nu.water, parcelsPerSecond) but nothing works.
I suppose that the problem is that a drag model also for the lower phase which is not interested by the initial injection event cannot be set in the kinematicCloudProperties dictionary (otherwise the simulation immediately crashes)
Any suggestion to solve this issue would be highly appreciated.
Kind Regards,
Attached Files
File Type: txt alpha.water.txt (1.0 KB, 4 views)
File Type: txt blockMeshDict.txt (1.4 KB, 1 views)
File Type: txt controlDict.txt (1.1 KB, 3 views)
File Type: txt fvSchemes.txt (1.2 KB, 3 views)
File Type: txt g.txt (768 Bytes, 2 views)
Danny_ is offline   Reply With Quote

Old   June 11, 2021, 11:21
New Member
Join Date: Jun 2021
Posts: 2
Rep Power: 0
Danny_ is on a distinguished road
Remaining files
Attached Files
File Type: txt kinematicCloudProperties.txt (4.7 KB, 9 views)
File Type: txt p_rgh.txt (1.3 KB, 2 views)
File Type: txt setFieldsDict.txt (947 Bytes, 2 views)
File Type: txt transportProperties.txt (973 Bytes, 5 views)
File Type: txt U.txt (1.1 KB, 2 views)
Danny_ is offline   Reply With Quote

Old   August 26, 2021, 17:49
Julio Pieri
Join Date: Sep 2017
Posts: 75
Rep Power: 7
JulioPieri is on a distinguished road
Dear Danny,

I haven't worked with interFoam coupled with LPT, however I'm using swak4Foam together with the twoPhaseEulerFoam - which is also multiphase. Since I don't have the correct solver, I'm unable to reproduce your case here. Do you have any pictures of the particles "sticking" to the interface?

Things I'd try if this were my case:
- change drag model to sphereDrag, or remove it completely
- change base pressure of your case to 1e5
- remove pairCollision model to simplify the solution and speed up the debugging.
- include particleTrap function
- introduce a disturbance on the alpha.water field and see if the particles at least oscillate together with the interface (suggesting a "crossing" difficulty) or if they are simply maintaining position (thus unrelated with alpha).
- add an interface force to pull particles into the water - this might be unrealistic but may be a way to bypass this issue.

I also wouldn't discart implementation issues, since this is a custom solver. Try running swak4Foam for debugging as well.

Lastly, I recently had a problem where my particles wouldn't move at all with the velocity field. After long hours of trying, I decided to simple reset the case and start from another one and that solved the issue. Sometimes, and I don't know why, this simply works.

Edit: also, double check you mu field in the kinematicCloud file. Since the drag is calculated from that, make sure you're getting the correct value. I faced and solved a similar problem yesterday, right after replying.

Best regards

Last edited by JulioPieri; August 27, 2021 at 11:56.
JulioPieri is offline   Reply With Quote

Old   August 30, 2021, 19:38
Francisco T
Join Date: Nov 2011
Location: Melbourne, Australia
Posts: 64
Blog Entries: 1
Rep Power: 13
frantov is on a distinguished road
are you aware of the solver mppicinterfoam ?

it is available on the version
frantov 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
Eulerian Multiphase Model vs Lagrangian Particle Tracking ajjadhav CFX 14 December 7, 2020 17:22
Wind turbine simulation Saturn CFX 58 July 3, 2020 02:13
Question about adaptive timestepping Guille1811 CFX 25 November 12, 2017 18:38
Particle Reynolds number calculation in Lagrangian tracking? jiejie OpenFOAM Running, Solving & CFD 5 July 6, 2012 05:47
LES + VOF + Particle tracking mittal OpenFOAM 0 June 29, 2010 07:41

All times are GMT -4. The time now is 04:01.