Calculation of the CFL-number

Felix_Sp · December 16, 2021, 05:06

Hello everyone,

I was just wondering how the calculation of the CFL-number works in REEF3D.
It is stated in https://doi.org/10.1016/j.compfluid.2016.09.01 that the CFL-number is calculated with the contribution of advection, diffusion and the maximum of the source terms:

$\Delta t \leq 2\left(\left(\frac{\left | u \right |_{max}}{dx}+D\right)+\sqrt{\left(\frac{\left | u \right |_{max}}{dx}+D\right)^2+\frac{\left | S \right |_{max}}{dx}}\right)^{-1}$

As far as I understand is that if I now apply D 20 2 and T 11 1 (and do not have sediment or heat transport) the equation simplifies to:

$\Delta t \leq \frac{dx}{\left | u \right |_{max}}$

I then use N 47 as a scaling factor, let's say $\alpha$ , for the calculated time step, e.g., 0.3, which leads ultimately to:

$\Delta t_{applied} = \Delta t * \alpha = \Delta t * 0.3$

Is this the way the actual time step is calculated? Or did I misunderstand something? Thanks for your help in advance!

kamath · December 16, 2021, 06:35

That is correct

Felix_Sp · December 16, 2021, 08:57

Thanks for the quick reply! You guys are always really fast!

kamath · December 17, 2021, 04:07

Thank you Felix

Hope you continue enjoying working with REEF3D

December 16, 2021, 05:06	Calculation of the CFL-number	#1
Felix_Sp Member Felix S. Join Date: Feb 2021 Location: Germany, Braunschweig Posts: 86 Rep Power: 6	Hello everyone, I was just wondering how the calculation of the CFL-number works in REEF3D. It is stated in https://doi.org/10.1016/j.compfluid.2016.09.01 that the CFL-number is calculated with the contribution of advection, diffusion and the maximum of the source terms: $\Delta t \leq 2\left(\left(\frac{\left \| u \right \|_{max}}{dx}+D\right)+\sqrt{\left(\frac{\left \| u \right \|_{max}}{dx}+D\right)^2+\frac{\left \| S \right \|_{max}}{dx}}\right)^{-1}$ As far as I understand is that if I now apply D 20 2 and T 11 1 (and do not have sediment or heat transport) the equation simplifies to: $\Delta t \leq \frac{dx}{\left \| u \right \|_{max}}$ I then use N 47 as a scaling factor, let's say $\alpha$ , for the calculated time step, e.g., 0.3, which leads ultimately to: $\Delta t_{applied} = \Delta t * \alpha = \Delta t * 0.3$ Is this the way the actual time step is calculated? Or did I misunderstand something? Thanks for your help in advance!

December 16, 2021, 06:35		#2
kamath Senior Member Arun Kamath Join Date: Nov 2014 Location: Trondheim, Norway Posts: 265 Rep Power: 13	That is correct __________________ Arun X years with REEF3D

December 17, 2021, 04:07		#4
kamath Senior Member Arun Kamath Join Date: Nov 2014 Location: Trondheim, Norway Posts: 265 Rep Power: 13	Thank you Felix Hope you continue enjoying working with REEF3D __________________ Arun X years with REEF3D

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Foam::error::PrintStack	almir	OpenFOAM Running, Solving & CFD	92	May 21, 2024 07:56
AMI speed performance	danny123	OpenFOAM	21	October 24, 2020 04:13
simpleFoam parallel	AndrewMortimer	OpenFOAM Running, Solving & CFD	12	August 7, 2015 18:45
Compressor Simulation using rhoPimpleDyMFoam	Jetfire	OpenFOAM Running, Solving & CFD	107	December 9, 2014 13:38
[blockMesh] BlockMeshmergePatchPairs	hjasak	OpenFOAM Meshing & Mesh Conversion	11	August 15, 2008 07:36

December 16, 2021, 08:57		#3
Felix_Sp Member Felix S. Join Date: Feb 2021 Location: Germany, Braunschweig Posts: 86 Rep Power: 6	Thanks for the quick reply! You guys are always really fast!