|
[Sponsors] |
November 8, 2009, 10:33 |
CFL number limit for InterFoam
|
#1 |
Senior Member
Suresh kumar Kannan
Join Date: Mar 2009
Location: Luxembourg, Luxembourg, Luxembourg
Posts: 129
Rep Power: 17 |
Hello everybody,
I have been using lesinterFoam for studying the breakup of liquid sheets. I have also been testing the solver for the CFL limit. Recently I read an article "Volume of fluid methods fro immiscible-fluid and free-surface flows" by Vinay R.Gopala, Berend G.M. van Wachem.Chemical Engineering Journal 141 (2008) 204-221 In this paper it has been specified a couple of times that for the inter gamma scheme that is implemented in interFoam it is required to have very small time steps, CFL<0.01 for keeping a sharp interFace. Based on tis comments I have done a comparative study of different CFL numbers for the same mesh and same boundary conditions. I studied CFL numbmers of 0.2,0.1,0.05,0.01 and 0.0075 for the liquid sheets. I found out that the cases are very much dependent on the CFL numbers, but the cases with CFL nummbers 0.01 and 0.0075 are closer to each other. Then I alos found in article by Henry Weller " Interface tracking capabilities of the Inter-Gamma Differencing scheme" that it ir required to use Co<1/2 for 2d cases and Co<1/3 for 3d cases. Then it is also specified in Gopalas paper that in order to use larger time steps subcycles are used to solve the VOF equation. For a given Courant number, the flow equations are solved initially and later the VOF equation is fractionally updated n times, where n represents the number of subcycles predefined in the code. Now my question is why has the CFL limit been specified differently in these two literatures, which one should I follow, should I follow the limit specified by Gopala in his paper to keep a sharp interface or the one specified by Henry in his article which looks quiet old to me. My second question what is the limit to which i can vary my CFL number if I can use more n gamma subcycles. Has anybody studied the accuracy of the solution by using higher CFL number basically higher time steps and increasing the ngamma subcycles. I mean if want to use a CFL number of 0.2 is it sufficent that I use 4 ngamma subcycles or shouldi use more. Is there any reference for this CFL limit and the ngamma subcycles. bye with regards K.Suresh kumar |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[mesh manipulation] Mesh Refinement | Luiz Eduardo Bittencourt Sampaio (Sampaio) | OpenFOAM Meshing & Mesh Conversion | 42 | January 8, 2017 13:55 |
DecomposePar unequal number of shared faces | maka | OpenFOAM Pre-Processing | 6 | August 12, 2010 10:01 |
[blockMesh] BlockMeshmergePatchPairs | hjasak | OpenFOAM Meshing & Mesh Conversion | 11 | August 15, 2008 08:36 |
Unaligned accesses on IA64 | andre | OpenFOAM | 5 | June 23, 2008 11:37 |
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues | michele | OpenFOAM Meshing & Mesh Conversion | 2 | July 15, 2005 05:15 |