
[Sponsors] 
November 12, 2005, 08:50 
Length scale for CFL in unstructured grids

#1 
Guest
Posts: n/a

Hi.
How do you calculate the length scale for CFL in unstructured grids? Since the grid is not always aligned with the flow, I'm using the square root of the control volume area  in 2D. Is that correct? TIA 

November 12, 2005, 13:00 
Re: Length scale for CFL in unstructured grids

#2 
Guest
Posts: n/a

Hi Marcio,
I'm not sure for CVFEM formulation but I think there isn't a more or less correct way for computing it. I've already implemented a 2D FEM incompressible flow program using: he=sqrt(2*Ae), where he is the element "length" and Ae is the area of the element. This element size is based on the side of the square with 2 times the area of the triangle. but I´ve also tested other approaches, for example: he=sqrt(4*Ae/pi), element size based on the diameter of the circle with the same area. and he = max(e_1, e_2, e_3), where e_i are the length of the edges of the element (3 for triangles, of course) Besides the CFL number, I use the element size to compute the stabilization terms of the formulation (SUPG and PSPG) and all these approaches produced similar results. Nowadays, I´ve been working in a 3D solver based on the same formulation (stabilized FEM) and using the element length based on the diameter of the sphere with the same volume of the element but, as I've cited before for the 2D case, there are also other ways of computing it (cubic square of 6 times the volume for tetrahedra for example). I know that I haven't helped you but it is always a good subject to discuss because I think that the research community has not found the best way of computing it yet. Thus, for the time step, sometimes I do a roughly approximation taking the minor element size from the mesh and computing the time step for the desired CFL. Regards Renato. 

November 13, 2005, 05:16 
Re: Length scale for CFL in unstructured grids

#3 
Guest
Posts: n/a

Dear Marcio,
There is no unique way of computing a length scale in unstructured meshes, and sqrt(A) is one way of doing it. Although, other alternatives are possible, the present one which you are using does give a reasonable measure, except in certain cases of heavy skewness etc..., so I presume that you can proceed with the same . Regards, Ganesh 

November 13, 2005, 23:03 
Re: Length scale for CFL in unstructured grids

#4 
Guest
Posts: n/a

Apply the finite volume scheme to a linear convection/diffusion equation and do a stability analysis to get the CFL condition. Then generalize this to your case by replacing the wavespeed of the linear equation with the eigenvalues of the system case. The resulting formula will turn out to be costly and people usually make approximations. Any representative measure of the element size like side of smallest length, smallest altitude, etc. will work but may not give the least restrictive condition. See Blazek for some discussions on this topic.


November 14, 2005, 07:26 
Re: Length scale for CFL in unstructured grids

#5 
Guest
Posts: n/a

Thank you all for your answers. I've tested this scheme with CVFEM this weekend and it worked. The difference from the usual finite element scheme is that here I use the control volume area instead of the element area  so the length scale is roughly 1/3 of the traditional FEM for a 9node finite element. Also, for skewed elements I'm going to test the smallest edge as the length scale, but I've had no problem with it when solving a laminar channel (Re = 500) flow with streamwiseskewed elements with a aspect ratio = 2 and 3.
Regards Márcio Ricardo 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Turbulence length scale and integral length scale  rizhang  CFX  2  April 22, 2016 07:22 
Integral Length Scale vs Turbulent Length Scale  Scott Nordsen  FLUENT  2  January 20, 2014 02:17 
LES inlet length scale  f.neilson  Main CFD Forum  0  January 31, 2005 14:26 
Characteristic Length Scale  Thomas P. Abraham  Main CFD Forum  2  April 20, 1999 17:16 