
[Sponsors] 
October 18, 2001, 08:30 
why compact scheme?

#1 
Guest
Posts: n/a

A question: Is there any disadvantages with compact difference scheme? If the answer is `no', why is compact scheme not so popular?


October 19, 2001, 13:30 
Re: why compact scheme?

#2 
Guest
Posts: n/a

Compact schemes are known to provide spectrallike spatial accuracies. Spectral schemes (FFT mostly) are often used for uniform grids. For nonuniform grids, you can either use Chebyshevspectral or compact schemes.
Much like spectral methods, compact schemes are highly nonlocal. If there are narrow regions with strong gradients in the flow (like thin exothermic flames, multifluid interfaces) it is not a good idea to use these methods. They tend to create dispersive behavior. 

October 22, 2001, 08:46 
Re: why compact scheme?

#3 
Guest
Posts: n/a

However, most compact schemes involve few nodes thus are high local, that is indiated by it's name. In fact, what I want to know it that since compact schemes are high accuracy in nature, why there are so few applications of them? Many guyes in this forum said `i use upwind scheme for advective terms and secondorder center difference scheme for viscous terms ...'


October 22, 2001, 12:27 
Re: why compact scheme?

#4 
Guest
Posts: n/a

Compact schemes are very expensive as the derivatives are also taken as unknowns. Eventhough, there has been a publication of using this derivative for the viscous terms. It necessarily involves inverting matrices to calculate the fluxes. I would also extrapolate that it would become harder with arbitrary grids (as is true with other schemes as well).
I agree with Paul that compact doesn't necessarily mean central or upwind. A compact scheme can be either! chidu... 

October 22, 2001, 13:08 
Re: why compact scheme?

#5 
Guest
Posts: n/a

For a given field F, you compute the derivatives at all points along a line (in a vector F_x) using a linear system A F_x = b. Both the matrix elements in "A" and vector elements in "b" are computed using compact stencils but that does not mean the derivative computation has a local character.
F_x = inv{A}.b inv{A} most certainly does not have a compact support. As to why compact schemes are not used often, one likely reason may be that they are not as accurate on curvilinear grids as they are on cartesian meshes. Also, building in the BCs into the linear system (when using a implicit solver) is not too convenient although not impossible. They also have a centrallike nature which leads to dispersive errors whereas the dissipative oddorder upwind/upwindbiased are often considered more stable for advection. 

October 23, 2001, 02:22 
Re: why compact scheme?

#6 
Guest
Posts: n/a

Kalyan's explanation is very illustrative. In that sense, many implicit compact schemes may not `compact'. However, these do exist some explicit compact schemes which don't necessary lead to a linear system. This maybe another question.
Thank you, kalyan and chidu. 

October 23, 2001, 11:55 
Re: why compact scheme?

#7 
Guest
Posts: n/a

Can you point me to references for 'explicit' compact schemes? My readings include papers by Lele and Krishnan. I was under the impression that the compactness of the stencil for a given truncation error is acheived through the implicitness of the derivatives!!
chidu... 

October 23, 2001, 11:58 
Re: why compact scheme?

#8 
Guest
Posts: n/a

Hi Kalyan,
I agree with you about the nonlocality implicit in the inversion process. Also, like I said it would be cumbersome on arbitrary grids. But I think, although most of the use of compact schemes have tended to be centrallike, it is possible to derive compact schemes that upwindbased, no? chidu... 

October 23, 2001, 14:23 
Re: why compact scheme?

#9 
Guest
Posts: n/a

I am sure you can build in some upwind type nature into the compact schemes.
I am no expert in schemes let alone the compact schemes but when I do LES I do have to care a great deal about the numerical errors. From what I see, compact schemes have been very attractive for temporal simulations but for spatial problems, the outflow BC is hard to handle. This is yet another problem with these schemes. As Chidu suggested, if some sort of upwinding were to be builtin, then perhaps this problem can be handled more reasonably. 

October 25, 2001, 09:11 
Noncentered compact schemes

#10 
Guest
Posts: n/a

The man you want to investigate is Andrei Tolstykh. He is a Russian academist who has investigated noncentered compact schemes for fluid dynamins for much of his life. He has a book, translated to english, entitled "High accuracy noncentered compact difference schemes for fluid dynamics applications" which describe how to get 3rd and 5th orders of accuracy from 3 point noncentered stencils. Hope this helps.


October 25, 2001, 13:40 
Re: Noncentered compact schemes

#11 
Guest
Posts: n/a

One can't be expected to buy a book just for noncentered compact schemes (especially if he/she does not work full time with scheme development). Can you provide us with us references in the form of papers.
Thank You. 

October 25, 2001, 22:10 
Re: Noncentered compact schemes

#12 
Guest
Posts: n/a

Most of his papers are in Russian. These papers can be found in Russian journals such as Zh. Vychisl. Mat. Mat. Fiz., or other equally clear titles. However, I would recommend searching author databases for "A. I. Tolstykh". In scanning the index of the book, it appears that the english translations of some of his papers can be found in "USSR Comput. Math. Math. Phys." Try looking in the following Vol, year, page listings:
18, 1978, 134147 16, 1976, 92100 21, 1981, 339354 25, 1985, 127134 29, 1989, 183194 amongst others. Once again, the easiest and most cohesive method is to check out his book as I listed earlier. The book seems straight forward to this selftaught CFD student, and will provide the most useful instruction, as opposed to hunting through journals. 

October 31, 2001, 23:31 
Re: why compact scheme?

#13 
Guest
Posts: n/a

I think compact scheme is very accurate for hyperbolic equation (see Ghrist's dissertation, The University of Colorado). Obtaining the derivatives is simple, you just solve tridiagonal matrices by using Thomas algorithm (for 4th & 6th order).It is difficult to handle the BC for the derivatives.


November 2, 2001, 11:00 
Re: why compact scheme?

#14 
Guest
Posts: n/a

It's true that compact scheme can be central or upwind. (the former is popularised by Lele in a 1992 JCP paper, the later is dicussed briefly by Hirsch in "Numerical Computation of Internal and External Flows vol1, which should be found easily in univeristy libraries).
Looking at Lele description, one can see quite clearly that the uniformity (as well as regularity) of Cartesian grid plays important role in supressing/removing higher order error.Essentially in compact scheme the multiplier of neighbouring derivative is set in a way to minimise error of the original/base explicit scheme. Thus one get higher order of accuarcy for a given compact stencil. It means that its application for curvilinear or arbitrary structured grid is not easily done. (see Smirnov, Lacor and Baelmans in ECCOMAS 2000 for curvilinear case). I did some study in solving linear problem of convectiong a double sine wave profile (2D case) in unstructured grid with the following conclusions: 1. in term of cost the scheme is as expensive as implicit solver. (generally less than 10 SGS iteration is enough to obatin the derivative in one time level) 2. accurate result can be obtained if the grid is formed by triangulating uniform structured grid as the dispersive error is minimised. 3. the error revert back to original/base explicit scheme if the grid becomes completely irregular such as the one produced by Delaunay triangulation since it is difficult to chose multiplier that minimise error in all direction (notably the cross derivative term) In the light of the result it is possible to see why that compact scheme may not be so popular. The coding complexity is not justified for arbitrary grid. Having said that, it is not a final word by any means. Further research on this topic is still being carried out Cheers 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
2nd order upwind scheme (Fluent and CFX)  Far  FLUENT  0  May 22, 2011 02:50 
how to understand high resolution scheme and high order scheme  iilw1314  Main CFD Forum  3  September 30, 2010 11:43 
Definition of limiter function for central dirrerencing scheme  sebastian_vogl  OpenFOAM Running, Solving & CFD  0  January 5, 2009 12:08 
Compact scheme for nonuniform grid?  Dong  Main CFD Forum  3  July 2, 2008 00:44 
extrapolation in MUSCL scheme  Chandra  Main CFD Forum  6  February 14, 2007 12:21 