|
[Sponsors] |
December 15, 2011, 23:55 |
Lagrangian Particle Tracking Initial step
|
#1 |
New Member
Bruce Hartley
Join Date: Jun 2011
Posts: 9
Rep Power: 15 |
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?
|
|
December 16, 2011, 00:50 |
|
#2 |
Senior Member
Rasmus Hemph
Join Date: Mar 2009
Location: Sweden
Posts: 108
Rep Power: 17 |
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 |
|
December 16, 2011, 01:40 |
Lagrangian particle Tracking getting positions and other information
|
#3 |
New Member
Bruce Hartley
Join Date: Jun 2011
Posts: 9
Rep Power: 15 |
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 |
|
December 16, 2011, 03:06 |
|
#4 |
Senior Member
Rasmus Hemph
Join Date: Mar 2009
Location: Sweden
Posts: 108
Rep Power: 17 |
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 |
|
September 2, 2015, 05:32 |
|
#5 | |
New Member
Join Date: Apr 2015
Posts: 9
Rep Power: 11 |
Sorry for reviving such an old thread, but I have probably the same problem like described in this thread:
Lagragian particles get the right force but wrong velocity It is the same factor of ~100 for me, too. it seems like for some reason the particle movement gets called several times under certain conditions which I yet don't understand. I traced my function calls, but there seems to be no difference between a timestep that calls the movement of the particle once and one that does it several times. EDIT: In case it's of interest here is the stack of function calls: Quote:
Ok, the problem was just that the given density and particle diameter didn't fit to the chosen mass. OpenFoam does everything correctly. It checks the particle's diameter and with the given density it calculates its mass and uses that and not what is given in the ParticleProperties under massTotal which is not the particles mass but the mass of all particles in the injected parcel. Last edited by Kabra; September 2, 2015 at 11:53. |
||
Tags |
lagrangian particles, starting conditions, wall function |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
alphaEqn.H in twoPhaseEulerFoam | cheng1988sjtu | OpenFOAM Bugs | 15 | May 1, 2016 16:12 |
ForcesCoeffs | ronaldo | OpenFOAM | 4 | September 14, 2009 07:11 |
Error log | vw.cfd | OpenFOAM | 6 | August 7, 2009 05:44 |
Computational time | sunnysun | OpenFOAM Running, Solving & CFD | 5 | March 16, 2009 03:32 |
Modeling in micron scale using icoFoam | m9819348 | OpenFOAM Running, Solving & CFD | 7 | October 27, 2007 00:36 |