CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

query regarding deltaT in multiphaseEulerFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 21, 2013, 08:15
Question query regarding deltaT in multiphaseEulerFoam
  #1
New Member
 
Join Date: Mar 2011
Posts: 10
Rep Power: 6
bigboss.abhishek is on a distinguished road
I am trying to simulate a 3 phase flow by using multiphaseEulerFoam solver. I start with deltaT = 0.001. But during long runs its value decreases. And it is never constant. My second query is that whenever i try to restart the run from t=0 , Openfoam never takes the original deltaT = 0.001 value. I am posting the snippets of my terminal window and control dict file.


Courant Number mean: 0.0570183 max: 0.49767
deltaT = 0.000768923
Time = 0.475434

MULES: Solving for alphawater
water volume fraction, min, max = 0.237401 -1.39992e-09 1
MULES: Solving for alphaoil
oil volume fraction, min, max = 0.331176 -1.28595e-11 1
MULES: Solving for alphaair
air volume fraction, min, max = 0.431422 -2.89095e-11 1
Phase-sum volume fraction, min, max = 1 0.999999 1
MULES: Solving for alphawater
water volume fraction, min, max = 0.237401 -3.13086e-09 1
MULES: Solving for alphaoil
oil volume fraction, min, max = 0.331176 -1.26736e-11 1
MULES: Solving for alphaair
air volume fraction, min, max = 0.431422 -2.9056e-11 1
Phase-sum volume fraction, min, max = 1 0.999999 1
MULES: Solving for alphawater
water volume fraction, min, max = 0.237401 -3.0326e-09 1
MULES: Solving for alphaoil
oil volume fraction, min, max = 0.331176 -1.24904e-11 1
MULES: Solving for alphaair
air volume fraction, min, max = 0.431422 -1.09878e-10 1
Phase-sum volume fraction, min, max = 1 0.999999 1
GAMG: Solving for p, Initial residual = 0.00186816, Final residual = 3.4607e-05, No Iterations 2
time step continuity errors : sum local = 2.95067e-07, global = 1.07225e-10, cumulative = -4.99832e-07
GAMG: Solving for p, Initial residual = 0.00130052, Final residual = 2.03657e-05, No Iterations 2
time step continuity errors : sum local = 1.73768e-07, global = 7.56413e-11, cumulative = -4.99757e-07
GAMGPCG: Solving for p, Initial residual = 0.000449752, Final residual = 3.46072e-08, No Iterations 5
time step continuity errors : sum local = 2.95375e-10, global = 2.05615e-11, cumulative = -4.99736e-07
ExecutionTime = 8560.98 s

i started from deltaT=0.001 but it decreased after some iterations


below is my controldict file


/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application multiphaseEulerFoam;


startTime 0;

startFrom latestTime;

stopAt endTime;

endTime 3;

deltaT 0.001;

writeControl adjustableRunTime;

writeInterval 0.2;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

runTimeModifiable yes;

adjustTimeStep yes;

maxCo 0.5;
maxAlphaCo 0.5;

maxDeltaT 1;


// ************************************************** *********************** //
bigboss.abhishek is offline   Reply With Quote

Old   August 21, 2013, 10:37
Default
  #2
Senior Member
 
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 195
Rep Power: 10
kwardle is on a distinguished road
In this solver and any others of the interFoam family with interfaceSharpening, the time step is by default adapted to maintain the stability criterion of the Co number. In your controlDict you have maxCo = 0.5. and adjustTimeStep = yes. While it is always a good idea to use Co-controlled time step with this type of solver (unless you know exactly what the flow is doing at all times), if you want to keep the time step from changing from your deltaT setting in controlDict, set "adjustTimeStep no;"

Alternatively, depending on how many nAlphaSubCycles you are using (fvSolutions) you could increase the value of maxCo instead to optimize the deltaT and still keep things stable. You want to keep maxCo/nAlphaSubCycles > 0.3.
Hope this helps.
-Kent
kwardle is offline   Reply With Quote

Old   August 22, 2013, 00:08
Default query regarding deltaT in multiphaseEuler
  #3
New Member
 
Join Date: Mar 2011
Posts: 10
Rep Power: 6
bigboss.abhishek is on a distinguished road
i changed the controldict as follows

application multiphaseEulerFoam;


startTime 0;

startFrom latestTime;

stopAt endTime;

endTime 3;

deltaT 0.001;

writeControl runTime;

writeInterval 0.2;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

runTimeModifiable yes;

adjustTimeStep no;

maxCo 0.5;
maxAlphaCo 0.5;

maxDeltaT 1;


// ************************************************** *********************** //

Its working fine..... ...thanku
bigboss.abhishek is offline   Reply With Quote

Old   August 22, 2013, 00:12
Lightbulb query regarding deltaT in multiphaseEuler
  #4
New Member
 
Join Date: Mar 2011
Posts: 10
Rep Power: 6
bigboss.abhishek is on a distinguished road
whenever i try to restart the run from t=0 , Openfoam never takes the original deltaT = 0.001 value ??? i am not able to find out why it happened. Although what i do now is that i have made a separate copy of my 0 time directory. whenever i need to take a fresh run i use it and do setFields. I am still not able to figure out the problem. Few people say its a bug.
bigboss.abhishek is offline   Reply With Quote

Old   November 25, 2014, 05:43
Default
  #5
New Member
 
herman eselte
Join Date: Nov 2014
Posts: 3
Rep Power: 2
hereselte is on a distinguished road
Hello,

just remove the directory named UNIFORM from that point you will continue.

H
hereselte is offline   Reply With Quote

Reply

Tags
controldict, deltat small, openfoam 2.1.1

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
interDyMFoam deltaT drops non-stop Virtual-iCFD OpenFOAM Running, Solving & CFD 5 June 12, 2015 04:30
how to query node coordinates? phsieh2005 OpenFOAM Post-Processing 2 August 16, 2013 01:48
multiphaseEulerFoam kwardle OpenFOAM 3 March 25, 2013 11:31
same geometry,structured and unstructured mesh,different behaviour. sharonyue OpenFOAM Running, Solving & CFD 13 January 2, 2013 23:40
deltaT in calculation of time step continuity errors misakagan OpenFOAM Running, Solving & CFD 0 August 27, 2010 10:13


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