CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Is it possible to add massless particles to an incompressible flow at solver level?

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

LinkBack Thread Tools Search this Thread Display Modes
Old   April 11, 2017, 11:01
Default Is it possible to add massless particles to an incompressible flow at solver level?
New Member
Laura Sumner
Join Date: Mar 2017
Posts: 13
Rep Power: 9
LauraSumner is on a distinguished road

I have created a solver based on pimpleDyMFoam which I can post process using paraview to find the pathlines of particles injected into the flow using particle Tracer etc.

However, this is incredibly slow and time consuming, and the main issue is the paraView just crashes.

I need to be able to find the paths of particles in my flow for long times and therefore large datasets. My paraView simply cannot handle this. Plus this is a pretty ugly way of finding such information.

My initial thought was to research how to post process this kind of thing using the command line instead of running paraView. However the most elegant solution would be some combination of solvers which have massless particles added to the flow, so that I can just output their coordinates at each timestep instead of post processing at all. The solution must run in parallel as well!

My questions are:
1) Is this possible?
2) If so, is my method below anywhere near the right kind of thing?
3) If it is, what am i doing wrong?!


I am trying to merge the icoUncoupledKinematicParcelFoam somehow with my solver, because I thought that because of this line:


it would read in the newly calculated velocity field (U) and evolve the cloud accordingly. I have looked at other people's solvers that I have found from researching this problem and it seems they do similar things, however the cases I have found have been for old versions of OpenFOAM (I use 3.0.1) so it's hard to compare/I cant compile their code.

So far I have added all the files I think I need to the options file of my new solver, and added these lines of code within my time loop:

Info<< "Evolving " << << endl;



mu =*rhoInfValue;


(i also added this at the beginning of my main:
"specify alternative cloud name. default is 'kinematicCloud'"


when i compile i get errors saying:

error: expected primary-expression before ‘.’ token
Info<< "Evolving " << << endl;

and that

error: ‘mu’ was not declared in this scope
mu =*rhoInfValue;

which is strange because I am sure all the header files called by the icoUncoupledKinematicParcelFoam are included and that solver returns no such errors.

Perhaps this is also a foolish idea to begin with, any insight would be appreciated!!


LauraSumner is offline   Reply With Quote

Old   September 30, 2019, 17:13
New Member
Join Date: Mar 2015
Posts: 24
Rep Power: 11
rasool_soofi is on a distinguished road

did you ever find out the solution? Having the same issue myself!
rasool_soofi is offline   Reply With Quote


dynamic mesh, incompressible flow, kinematiccloud, particle tracking

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
How is the accuracy using VOF (pressure based solver) on supersonic flow air_fun Fluent Multiphase 1 August 25, 2021 00:59
[ANSYS Meshing] Help with element size sandri_92 ANSYS Meshing & Geometry 14 November 14, 2018 07:54
dispersion model with lagragian particle tracking model for incompressible flows eelcovv OpenFOAM Running, Solving & CFD 54 April 10, 2018 09:36
How to simulate dilute solid particles in gas flow? chpjz0391 OpenFOAM Running, Solving & CFD 4 March 22, 2016 19:32
Incompressible flow solver (staggered grid) J. Ehrhard Main CFD Forum 1 October 8, 1998 19:47

All times are GMT -4. The time now is 10:06.