CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Wiki > Overset grids

Overset grids

From CFD-Wiki

Revision as of 01:11, 29 April 2006 by Jasond (Talk | contribs)
Jump to: navigation, search

Contents

Introduction

A common difficulty in simulating complex fluid flow problems is that not every geometry can be well represented using a single, contiguous (structured or unstructured) grid. In many cases, different geometrical features are best represented by different grid types. One approach these difficulties is the construction of a grid system made up of blocks of overlapping structured grids. This technique is referred to as the Chimera or overset grid approach. In a full Chimera grid system, a complex geometry is decomposed into a system of geometrically simple overlapping grids. Boundary information is exchanged between these grids via interpolation of the flow variables, and many gridpoints may not be used in the solution (these points are sometimes called hole points). Each block has boundary or fringe points, which lie in the interior of a neighboring block (or blocks) and will require information from that containing block. The data that must be generated to successfully complete a Chimera-type calculation is not insignificant, and thus has been automated to a high degree. In very general terms, there are three steps to setting up an overset simulation:

  1. Grid generation
  2. Hole cutting
  3. Determination of interpolation weights

Note that in some systems, one or more (or all) of these steps may be combined. Finally, while the Chimera technique is most often associated with traditional finite volume/difference CFD codes, it can in principle be applied with other discretization schemes.

Grid generation

Grids for an overset simulation are generally simple and structured, and are often generated hyperbolic or marching techniques. (Add more here)

Hole Cutting

Add example here

Chimera interpolation

Many of the tools required to compute the data required to do an overset simulation can be built upon the trilinear interpolation formula. The discussion here assumes hexahedral cells with nodal representation of flow variables available. Given a hexahedral cell formed by eight points \vec{x}_1,\vec{x}_2,\ldots\vec{x}_8 (see figure, which needs to be redrawn), the coordinates of any point in the interior of this cell may be written as


\begin{matrix}
\vec{x}=\vec{x}_1 & + &\left(-\vec{x}_1+\vec{x}_2\right)dj\\
                  & + &\left(-\vec{x}_1+\vec{x}_4\right)dk\\
                  & + &\left(-\vec{x}_1+\vec{x}_5\right)dl\\
                  & + &\left(\vec{x}_1-\vec{x}_2-\vec{x}_4+\vec{x}_3\right)dj\ dk\\
                  & + &\left(\vec{x}_1-\vec{x}_2-\vec{x}_4+\vec{x}_3\right)dj\ dk\\
                  & + &\left(\vec{x}_1-\vec{x}_2-\vec{x}_5+\vec{x}_6\right)dj\ dl\\
                  & + &\left(\vec{x}_1-\vec{x}_4-\vec{x}_5+\vec{x}_8\right)dk\ dl\\
                  & + &\left(-\vec{x}_1+\vec{x}_2-\vec{x}_3+\vec{x}_4
                    +\vec{x}_5-\vec{x}_6+\vec{x}_7-\vec{x}_8\right)dj\ dk\ dl,\\

\end{matrix}

(This equation is really ugly)

where dj, dk, and dl are the interpolation weights and are in the interval [0,1]. For any boundary point \vec{x}_b and any cell base point \vec{x}_1 we can solve the interpolation for for the weights. If the cell does not contain \vec{x}_b, then one or more of the weights will not be in the proper range. This can be used to construct a Newton-like gradient search technique. For more details, the reader is referred to the chapter by Meakin (1999) and the references contained therein. In practice, the initial condition for this solution procedure is quite important. It is usually best to perform some spatial partitioning to narrow the search range.

Software

To use the overset approach, one will generally need two software components: a preprocessing program to generate the overset data (cut holes, interpolation weights, etc.) and a compatible flow solver.

Preprocessors

A number of preprocessing packages have been developed, with varying levels of complexity. NASA's PEGASUS [Suhs et al 2002], which is primarily a grid joining code, takes an existing grids and prepares it for use in an overset simulation. The grid is generated separately. Another program is Chalmesh [Petersson (1999)], which generates grids and the interpolation data simultaneously. A comprehensive package, called Chimera Grid Tools (CGT), has been developed (primarily by NASA employees) that is intended cover all aspects of the preprocessing , including grid generation. The CGT package includes many utilities, including PEGASUS [Suhs et al (2002)], that automate most of the process of going from CAD model to simulation. Of these three packages, only Chalmesh is freely available. Unfortuately, the other two are subject to U.S. goverment export controls, and are thus not generally available (see here for more on this). However, since most of the algorithms and techniques used are described in the literature, it is possible to do this sort of simulation without access to these packages. (Commercial Packages?)

Solvers

  • INS2/3D
  • OVERFLOW
  • OVERTURE/OVERBLOWN
  • CFDSHIP-IOWA

References

Meakin, Robert L. (1999), "Composite Overset Structured Grids", Chapter 11, Handbook of Grid Generation, CRC Press.

Petersson, N. Anders (1999), "Hole-Cutting for Three-Dimensional Overlapping Grids", SIAM Journal on Scientific Computing, Vol. 21, No. 2, pp 646-665.

Suhs, Norman E. and Rogers, Stuart E. and Dietz, W. E. (2002), "PEGASUS 5: An Automatic Pre-Processor for Overset-Grid CFD", AIAA Paper 2002-3186, AIAA 32nd Fluid Dynamics Conference, St. Louis.

External links

My wiki