CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   Lagrangian Particle Tracking Initial step (http://www.cfd-online.com/Forums/openfoam-programming-development/95405-lagrangian-particle-tracking-initial-step.html)

Bruce Hartley December 16, 2011 00:55

Lagrangian Particle Tracking Initial step
 
I have modified some particle tracking with OpenFoam to track particles in electromagnetic fields using the standard lagrangian solvers. I find that the first step from T=0 to the next time step is not consistent with the following steps. The particle's move is relatively large for the first time step by a factor of about 100 larger than the following ones and in a direction which is not compatible with the field. It seems to be dependent on the starting positions of the particles so may depend on the wall boundary conditions. Can anyone suggest what might be affecting the first step in the particle tracking and indicate how it might be overcome. Has anyone else noticed such a problem?

hemph December 16, 2011 01:50

Hi Bruce,
If you have done any restarts before, the first time step is governed by the deltaT that is in the uniform/time dictionary. This has tricked by this many times..

//Rasmus

Bruce Hartley December 16, 2011 02:40

Lagrangian particle Tracking getting positions and other information
 
Hi, Rasmus,

Thank you for the reply but I'm afraid that I don't exactly understand what is required. Do I need to do something to the Control Dictionary or do I need to add another dictionary which specifies other time intervals or what? I have looked at wall conditions but can't seem to affect anything through them.

Can you give me an example of what I may need to do?

Thanks

Bruce Hartley

hemph December 16, 2011 04:06

Sorry if I was not clear enough. If you are in your case-directory and change dir to a timestep other than 0, say 1.0, you will find a directory called uniform. Inside this directory you will find the time-dictionary, where deltaT is specified for this time step.

If you restart from 1.0 instead of 0, the first time step is read from this deltaT instead of the deltaT in system/controlDict.

//Rasmus


All times are GMT -4. The time now is 13:29.