# Presentation: 356559

## Geometry Modeling & Grid Generation

## The objectives of this discussion are to relate experiences and  
offer some practical advice with regard to: 
Modeling geometry for use with CFD analysis
Generating a grid for the CFD analysis

Geometry modeling and grid generation are often the most 
difficult and time-intensive aspects of a CFD analysis.

- offer some practical advice with regard to:
- Modeling geometry for use with CFD analysis
- Generating a grid for the CFD analysis
- Geometry modeling and grid generation are often the most
- difficult and time-intensive aspects of a CFD analysis.

- Objectives

- Focus will be on multi-zone, structured grids.

## The data describing the geometry of the body about which the 
flow is to be computed can be provided 
in a variety of means:
Coordinates of points
Blueprint mechanical drawings
CAD files
CFD analysis mainly needs the surfaces 
that “see” the flow.  
These surfaces need to be extracted from 
the geometry data in some manner.

- flow is to be computed can be provided
- in a variety of means:
- Coordinates of points
- Blueprint mechanical drawings
- CAD files
- CFD analysis mainly needs the surfaces
- that “see” the flow.
- These surfaces need to be extracted from
- the geometry data in some manner.

- Geometry Data

## Need to be reasonable with the level
of  geometric detail to be included.

Greater geometric detail leads to more 
grid points and longer computations. 

A sequential approach is best in which 
the detail of the geometry is increased 
as the CFD analyses proceed.

- of  geometric detail to be included.
- Greater geometric detail leads to more
- grid points and longer computations.
- A sequential approach is best in which
- the detail of the geometry is increased
- as the CFD analyses proceed.

- Level of Geometric Detail

- For the NASA VDC inlet, neglected VGs, struts, leaves on centerbody, and slot.

## Availability of CAD files for defining the geometry provides needed 
information in a ready form; however, issues are:
CAD files usually contain much more information than is needed by the CFD analysis (mechanical details).  Can take considerable  effort to extract out the curves and surfaces that are needed.
Many CAD packages are solid modelers.  They can output surfaces; however, the “quilt” of surfaces can have some strange trimmed surfaces that have gaps in their edges.
CAD files are usually generated by mechanical designers that are mindful of tolerances and clearances needed in the mechanical design; however, CFD desires no clearances or gaps.  CFD analyst may need to spend time to “close up” the surfaces.
CAD file format.   IGES file format is commonly used since it does well representing surfaces and curves.
Mechanical drawings contain all the fine features; however, CFD may not wish to capture such geometric resolution.  CFD analyst has to remove such features.

- information in a ready form; however, issues are:
- CAD files usually contain much more information than is needed by the CFD analysis (mechanical details).  Can take considerable  effort to extract out the curves and surfaces that are needed.
- Many CAD packages are solid modelers.  They can output surfaces; however, the “quilt” of surfaces can have some strange trimmed surfaces that have gaps in their edges.
- CAD files are usually generated by mechanical designers that are mindful of tolerances and clearances needed in the mechanical design; however, CFD desires no clearances or gaps.  CFD analyst may need to spend time to “close up” the surfaces.
- CAD file format.   IGES file format is commonly used since it does well representing surfaces and curves.
- Mechanical drawings contain all the fine features; however, CFD may not wish to capture such geometric resolution.  CFD analyst has to remove such features.

- CAD Geometry Data

## The geometry of the flow domain requires 
modeling much in the same manner as the 
body geometry.
Outer boundary of the flow domain needs to be out enough to avoid bad interactions between boundary conditions and the flow.
Constructing the flow domain usually requires sketching out the topology of the grid (especially for multi-zone, structured grids in which each zone must adhere to a structured topology).
May need to go to great lengths to fit a structured grid topology onto a geometry. ( “clamshells”, singular axis, nose cells) 

- modeling much in the same manner as the
- body geometry.
- Outer boundary of the flow domain needs to be out enough to avoid bad interactions between boundary conditions and the flow.
- Constructing the flow domain usually requires sketching out the topology of the grid (especially for multi-zone, structured grids in which each zone must adhere to a structured topology).
- May need to go to great lengths to fit a structured grid topology onto a geometry. ( “clamshells”, singular axis, nose cells)

- Flow Domain Geometry

## The geometry of the body and flow domain are then input to the grid 
generator for the generation of the grid. 
A surface “database” is used in the grid generator for projecting the surface grid so that the surface grid is true to the geometry.
Flow domain boundaries are used to construct zonal boundaries.
Most grid generators read IGES files.

- generator for the generation of the grid.
- A surface “database” is used in the grid generator for projecting the surface grid so that the surface grid is true to the geometry.
- Flow domain boundaries are used to construct zonal boundaries.
- Most grid generators read IGES files.

- Geometry Input to the Grid Generator

## Grid is “body-fitted” to follow the shape of the body (wing). 
Grid points (vertices) are arranged in an array structure with indices (i,j,k). 
Grids can end up with 10s or 100s of zones and millions of grid points.
Grid cell size varies to provide proper resolution of flow gradients.
Zones and cells are hexahedral.

- Grid points (vertices) are arranged in an array structure with indices (i,j,k).
- Grids can end up with 10s or 100s of zones and millions of grid points.
- Grid cell size varies to provide proper resolution of flow gradients.
- Zones and cells are hexahedral.

- Multi-Zone, Structured Grid Basics

- , i

- , j

- , k

- The grid is often provides the only geometric information.

## Single-Zone Topology 

- A zonal boundary of a zone may be topologically connected
- to itself or other boundaries of the same zone:

- Nose zone

- Zone with singular axis

- Self-closing zone

## Flow information is exchanged across zonal boundaries.
Connectivity defines how a zonal boundary is connected to other zonal boundaries.

- Flow information is exchanged across zonal boundaries.
- Connectivity defines how a zonal boundary is connected to other zonal boundaries.

- Multi-Zone Connectivity

- 1 & 3

- 2

- 4

- GMAN can automatically compute zonal connectivity for abutting (1 & 2).
- Overlapping (overset) grids use tri-linear interpolation to exchange information.  This works best when cells in overset region are approximately the same size.

- Types of Zonal Connectivity
-    1)  Abutting, Point-to-point match
-    2)  Abutting, Non point-to-point match
-    3)  Overlapping, Point-to-point match
-    4)  Overlapping, Non point-to-point match

- Outer boundary

- Interior boundary

## Overlapped Grids   

- Example of NLR Airfoil with flap
- Flap grid overlaps airfoil grid.
- Use flap grid to cut a hole in airfoil grid.
- Indicate edge points of cut are fringe points to interpolate boundary information from flap grid solution.
- Outer edge of flap grid is fringe boundary of flap grid and receives interpolated data from the airfoil grid.

- Fringe points

- Hole

- Fringe points

## Remarks on Zonal Connectivity

- Avoid placing zonal boundaries in regions with large flow gradients (shocks, boundary layers, shock / boundary layer interaction).  *This is often unavoidable.*
- Best if zonal boundaries are place *normal* to the flow direction.
- Placing zonal boundaries parallel to the flow in regions of high gradients is generally not good, but often unavoidable.
- Point-to-point matching across zones is best since inviscid flux information can be directly transferred and errors are less in interpolation of viscous and turbulent information.
- Best if zones overlap by at least two grid points, especially in turbulent flows.  This improves the interpolation.  However, this may be difficult for grids with lots of zones and complicated geometry.
- Tools are available for combining, splitting, and overlapping zones in a grid

## Grid Clustering at Walls

- Walls with boundary layers, require the grid to be
- “clustered” normal to the wall to resolve the large
- gradients through the boundary layer.
- Non-dimensional distance from a wall in a
- turbulent boundary layer is given by y+
- General rule is that first grid point off the wall
- should have y+ < 1 to accurately define a turbulent
- velocity profile.  Can go up to y+=5 for less
- accuracy.  Resolution of heat fluxes requires y+
- about 0.1.
- With a smooth grid about 5-10 grid points
- will be in the viscous sublayer (y+ < 30).
- Even inviscid wall needs a little clustering.

- y1 = 1.0E-05 ft
- for  c = 1.0 ft
- 6 points in sublayer
- 6 points in log layer

## Grids for Wall Functions

- A wall function computes properties in the
- log layer (15< y+ < 100) of an attached
- turbulent boundary layer at a wall.
- Can remove the grid points in the log layer.
- This reduces the overall number of grid points.
- Also removes the smallest grid cells that inhibit iterative convergence.
- The full grid is generated and then the log
- region points are removed (CFSUBSET
- utility does this).

- Full grid.

- Grid for use with a wall function.

## Grid Sequencing

- Grid sequencing in a CFD solver is when the solution is performed on only
- certain grid points, such as only every other grid point along each grid line.
- This essentially coarsens the grid.
- Fewer grid points reduces the number of computations, and so, CPU time.
- Larger grid spacings help damp higher frequencies that limit iterative convergence.
- This is a good approach for starting the solution, damp out initial transients, and
- get a coarse-grid solution for a grid convergence study.
- The requirement for grid generation is that the number of grid points along a grid
- line *N* be constrained by
- where *n* is the number of levels for sequencing and *m* is some integer.  One level
- of sequencing would mean every other grid point is solved, and so, require *n* = 1.

## Grid Stretching and Smoothness

- Grid stretching is the ratio of distances between adjacent grid points along a
- grid line.
- Stretching is defined so that *S*  1.
- Can express as a percentage as      *S% = (S-1)*100%*
- General rule is to keep stretching below 15-20% to keep numerical errors
- within bounds.  This is especially important in regions with strong flow
- gradients.  Stretching upwards of 30% to 35% could be used for simulations
- where accuracy demands are less.
- A smooth grid is one in which the stretching varies smoothly over the grid.

## Grid Orthogonality

- Grid orthogonality is the angle that a grid line makes with the other grid lines
- intersecting at a grid point.
- Orthogonality is defined so that **  90o.
- General rule is to keep orthogonality as close to 90o as possible, especially normal
- to solid surfaces to provide accurate estimates of the normal to the walls for wall
- boundary conditions.  A minimum angle is probably 45o.

## CAD and Grid Generation Tools

- CAD Tools available here at Glenn:
- AutoCad
- Pro/Engineer
- Unigraphics
- ICEM DDN
- Grid Generation tools available at Glenn
- GRIDGEN (available via setup on Unix system)
- ICEM CFD (licenses held by some people)
- SmaggIce (Iced airfoils)
- NPARC System provides tools for scaling, rotating, manipulating, translating grids and interpolating solutions onto new grids.
- NASA Langley has number of grid tools (VGRID, GeoLab)
- Variety of smaller, specialized codes

- Next 