# Discretization method for DNS

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 January 23, 2009, 09:15 Discretization method for DNS #1 chandra Guest   Posts: n/a Hi, I would like have some suggestion from those who has performed DNS for turbulent flows. (1). Which spatial discretization method is suitable and of what order? I looked into the 6th order compact schemes but it seems quite expensive to obtain all the derivatives after solving huge matrices. So, I'm making my mind to use a 6-th order Central-Difference scheme on non-uniform grids. However, I am not sure how should I derive all those 6-th order expressions for various derivatives. Any efficient idea? Also, if I use 6th order central-difference, I need to extrapolate the variables for three imaginary cells on the boundary. Would it bring any significant inaccuracy? (2). Same is true with the time-marching scheme. I'm currently using a special 3rd order Runge-Kutta scheme for Differential Algebraic Equations (DAE) proposed by M.Arnold: "Half-Explicit Runge-Kutta Methods with Explicit Stages for Differential-Algebraic Systems of Index 2". This is the link to download it: http://sim.mathematik.uni-halle.de/~...s/1998/A98b.ps I'm planning to change the time-marching scheme also, either with the general Runge-Kutta schemes for ODEs or if there is any other better and efficient scheme. Please suggest me with your opinions. Thank you very much in anticipation!

 January 24, 2009, 03:48 Re: Discretization method for DNS #2 Harish Guest   Posts: n/a The best way would be to use explicit central difference scheme of 6th order/ Spectral methods along the periodic direction and non uniform 6th order compact scheme in the wall normal direction. Do not use coordinate transformation in the wall normal direction as it would introduce numerical errors which might introduce numerical errorsn. Refer to the non uniform compact scheme paper by Gamet et.al in journal of computational physics. They performed DNS of compressible channel flow using non uniform compact scheme. For time marching use Runga-Kutta method implemented in low-storage form. This would help in the reduction of memory requirements by a factor of 2 or 3. Hope this helps. Ravindra Shende likes this.

 January 24, 2009, 09:15 Re: Discretization method for DNS #3 chandra Guest   Posts: n/a Thanks for your reply. I looked at Gamet's paper in "International Journal for Numerical Methods in Fluids" with title "Compact Finite Difference Schemes in Non-Uniform Meshes. Application to Direct Numerical Simulations of Compressible Flows". However, I couldn't find his paper in JCP. Could you please be more specific about it. In the Gamet's paper I mentioned, he has described the compact scheme of 4th-order on non-uniform mesh using 5-point stencils resulting in a tri-diagonal system. That looks attractive. However, if I want to use 6th-order non-uniform mesh using 3-point stencils, "Chu & Fan" have proposed some schemes on non-uniform grids. But, the problem with using that is the resulting system is not a band-diagonal system. Instead, it seems to be a general sparse matrix which is computationally very expensive. If I carry on Gamet's idea, I think, I need to use a 7-point stencils for 6th-order scheme. It is same as a 6th order CD-scheme. So, I still have to figure out which one to use Please suggest your opinions! It will be very helpful.

 January 24, 2009, 09:37 Re: Discretization method for DNS #4 chandra Guest   Posts: n/a I think, I don't need to use 7-point stencil in order to get 6th-order compact scheme from Gamet's idea. 5-point stencil is sufficient. Also, I read that Compact Schemes are much more accurate that ordinary CD-schemes. Perhaps, that's the reason why people prefer compact schemes over CD-schemes in DNS.

 January 24, 2009, 21:50 Re: Discretization method for DNS #5 Harish Guest   Posts: n/a I think that i gave the wrong journal. It was published in IJNMF. Compact schemes resolve higher range of numerical wavenumbers compared to ordinary CD schemes and hence are more accurate. The problems with using implicit scheme if you have periodic directions is that the tridiagonal solver becomes much more expensive than a non periodic tridiagonal solver. It would be best to use Upwind compact scheme of Zhong or DRP scheme of Tam and Web in the periodic direction to make the computations faster. Also explicit schemes in certain directions will make it easier to parallelize the code in slices. Ravindra Shende likes this.

 January 25, 2009, 09:44 Re: Discretization method for DNS #6 chandra Guest   Posts: n/a Thank you very much for the valuable information. Now, I can understand why not to use the compact scheme in periodic direction. I have one final question though: I am using the classic MAC method by Harlow and Welch for my computation. The grid is staggered i.e. velocity vectors are defined on the face centers and the scalars are on the cell-centers. The compact schemes what I have come across until now or what you have suggested, all are for collocated grids, not for the staggered grids. If u have any suggestion in this regard or any reference, could u please let me know? It will indeed be very helpful to me. Thank you very much once again.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Andyjoe Open Source Meshers: Gmsh, Netgen, CGNS, ... 13 March 14, 2012 05:35 D.S.Nasan Main CFD Forum 6 September 4, 2008 02:08 student Main CFD Forum 2 June 14, 2008 16:43 mneo FLUENT 1 June 7, 2005 10:26 CFDtoy Main CFD Forum 0 February 15, 2005 13:51

All times are GMT -4. The time now is 13:42.

 Contact Us - CFD Online - Privacy Statement - Top