
[Sponsors] 
November 15, 2005, 14:59 
Mesh Size

#1 
Guest
Posts: n/a

Hello Friends !!
How i calculate Local mesh size for 2D and 3D?? My project is Compressible flow analysis on unstructured meshes. Now i am working in Quadratic reconstruction. At present i am developing code for quadratic reconstruction. for that i have to calculate the gradients with 2nd order accurate. It accomplished by both GreenGauss and Leastsquare...I am confused in Constant coefficient matrix for Least square method. so plz give the clear view of quadratic reconstruction. One more help... I need this below mentioned paper for my further progress.. Baldwin and Lomax "Thin layer approximation and algebraic model for separated turbulent flows", AIAA Paper, No. 78257, 1978. 

November 17, 2005, 02:19 
Re: Mesh Size

#2 
Guest
Posts: n/a

Dear Mohamed Yusuf,
Both GG and QR can be used for obtaining the gradients. However, not that GG by itself cannot give secondorder accurate gradients on arbitrary meshes. As far as your question on QR is concerened, this is the concept. 1. Consider any cell 'i' surrounded by support cells , 'j', j running from 1 to N. Thus cell 'i' has N supports. 2. We need to find the gradients at 'i'. 3. For every cell 'i' and support 'j', the cellcentroidal values \phi_i and \phi_j are known. Morevoer, \phi_j can be expanded in Taylors series about 'i' and the gradients and hessians at 'i' 4. To get the gradients and hessians we need to minimise the error between \phi_j at cell centroid of cell 'j' and \phi_j obtained through TS expansion. We minimise the sum of square of the error w.r.t to gradients and hessians,ie d(\sum e^2)/d(\phi_x) = 0 etc .. where denominator represents the gradient in x for the qty. phi and the sum is over all the supports. 5. This would give us a system of equations Ax = B, where x would have the gradients and hessians ( hence a coloumn vector of 5 in 2d and 9 in 3d). A is the geometric matrix which would only be a fucntion of the distances between cell centroids, boyh in x and y coordinates, mmed over N supports in suitable fashion ( 5*5 and 9*9 in 2d/3d) and b is the RHS vector, depending on the difference in solution as well as the distance b/w cellcentroids. 6. To get x , the matrix A needs to be inverted. For QR, this could be a cumbersome operation in 3D. The inversion is achieved using Singular value decomposition or Householder's algorithm ( See Numerical Recipes in Fortran/C for details) Note that the number of supports in 2D for QR needs to be atleast 5 and 610 is a widely accepted number. Howvere, in 3d this number could be excessively large like 40, and then computational efforts could pick up. You can also look into papers regarding Defect Correction in this regard, for 3d applications. As far as the B/L turbulence model is concerned if you could kindly send me your mailid, I could send you the paper in .pdf format Hope this helps Regards, Ganesh 

November 17, 2005, 08:41 
Re: Mesh Size

#3 
Guest
Posts: n/a

Thanks for ur reply...
But in Delaneye and Essers 1997 AIAA paper.. they considered 2nd order accurate gradient.. Then one more help.. In viscous computation, I am facing problem with calculation of gradient at faces(edges). I used GG theorem.. but i am not getting exact zero velocity at Noslip wall. Somewhere i am making mistakes but i cant realize.. Any other approach to calculate gradients for viscous? Advance thanks... 

November 17, 2005, 11:43 
Re: Mesh Size

#4 
Guest
Posts: n/a

Dear Yousuf,
In Delanye and Essers 1997 AIAA paper, they have first used GG to obtain gradients to first order and then corrected it to second order. This is possible option you could try. The fact that GG cannot by itself give gradients to second order accuracy on arbitrary meshes is still valid. In calculating the gradients at the boundary, use the boundary condition (say noslip)and consider a covolume made up by the nodes of the face and the cell centroid of the cell containing the boundary face. The viscous gradient computations are calculated using GG/QR. The AIAA 97 paper you referred to is a very good reference. Hope this helps Regards, Ganesh 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
SnappyHexMesh for internal Flow  vishwa  OpenFOAM Native Meshers: snappyHexMesh and Others  24  June 27, 2016 08:54 
Mesh motion with Translation & Rotation  Doginal  CFX  2  January 12, 2014 07:21 
Gambit Unstructed (Tet) Mesh  size control  newbie  FLUENT  3  August 26, 2008 11:38 
Icemcfd 11: Loss of mesh from surface mesh option?  Joe  CFX  2  March 26, 2007 18:10 
How to decide the mesh size?  echo  CDadapco  1  May 23, 2006 05:34 