zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: streamLineParticle : was using GREAT for deltaT which triggers lots of problems.

Browse Source 
Now using sqrt(GREAT).
This commit is contained in:
mattijs
2010-09-29 13:05:38 +01:00
parent fdb1c0be4b
commit 5c0c1f63e0
1 changed files with 6 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
src/postProcessing/functionObjects/field/streamLine/streamLineParticle.C
View File

@ -157,7 +157,9 @@ bool Foam::streamLineParticle::move(streamLineParticle::trackData& td)
td.switchProcessor = false; td.switchProcessor = false;
td.keepParticle = true; td.keepParticle = true;
scalar deltaT = GREAT; //cloud().pMesh().time().deltaTValue(); // Set very large dt. Note: cannot use GREAT since 1/GREAT is SMALL
// which is a trigger value for the tracking...
scalar deltaT = Foam::sqrt(GREAT); //cloud().pMesh().time().deltaTValue();
scalar tEnd = (1.0 - stepFraction())*deltaT; scalar tEnd = (1.0 - stepFraction())*deltaT;
scalar dtMax = tEnd; scalar dtMax = tEnd;
@ -169,6 +171,9 @@ bool Foam::streamLineParticle::move(streamLineParticle::trackData& td)
&& tEnd > ROOTVSMALL && tEnd > ROOTVSMALL
) )
{ {
// TBD: implement subcycling so step through cells in more than
// one step.
// set the lagrangian time-step // set the lagrangian time-step
scalar dt = min(dtMax, tEnd); scalar dt = min(dtMax, tEnd);