Grid classification

Last updated June 19, 2024

In applied mathematics, a grid or mesh is defined as the set of smaller shapes formed after discretisation of a geometric domain. Meshing has applications in the fields of geography, designing, computational fluid dynamics,^[1] and more generally in partial differential equations numerical solving. The geometric domain can be in any dimension. The two-dimensional meshing includes simple polygon, polygon with holes, multiple domain and curved domain. In three dimensions there are three types of inputs. They are simple polyhedron, geometrical polyhedron and multiple polyhedrons. Before defining the mesh type it is necessary to understand elements (their shape and size).

Elements

An element of a mesh is one of the smaller shapes evoked in the introduction. The shape of the elements is of great importance in solving problems in computational fluid dynamics. They are typically based on aspect ratio i.e. the aspect ratio of element decide whether a particular element would be good to use or we should go for another element with different aspect ratio. For example, if the aspect ratio is large the speed of solver reduces while if this ratio is small the solver speed increases. Large aspect ratio has another limitation of leading to interpolation errors. But if the results vary with direction then we use large aspect ratio.

Fluid flow equation and coordinate system

Most of the fluid flow equations are easily solved by discretizing procedures using the Cartesian coordinate system.^[2] In this system the implementation of finite volume method is simpler and easier to understand. But most of the engineering problems deal with complex geometries that don’t work well in the Cartesian coordinate system. When the boundary region of the flow does not coincide with the coordinate lines of the structured grid then we can solve the problem by geometry approximation. Figures 1a. and 1b. shows how a cylinder can be approximated with the Cartesian coordinate system.

The curve geometry of cylinder in Cartesian coordinate system is approximated by using stepwise approximation. But this method requires large time and is very tedious to work with. Other than this problem there is one more problem which is the cells inside the solid part of the cylinder, which are called dead cells, are not involved in the calculations so they should be removed, otherwise they would consume extra space in computer or other resources. Stepwise approximation is not smooth and thus leads to significant error, though the grid can be refined by using a fine mesh to cover the wall region but this leads to waste of computer memory resources.

Therefore, there are limitation in using methods in computational fluid dynamics based on simple coordinate system (Cartesian or cylindrical) as these systems fails while modeling of complex geometries like that of an aerofoil, furnaces, gas turbine combustors, IC-engine etc.

Classification of grids in computational fluid dynamics

a) Structured curvilinear grid arrangements (vertices having similar neighborhood).

b) Unstructured grid arrangements (vertices having variation in neighborhood).

Structured curvilinear grids

1) Grid points are identified at the intersection of coordinate line.

2) There are fixed number of neighboring grids for interior grid.

3) They can be arranged into an array and can be named by indices I, J, K f (In three dimensions).

These are also known as body-fitted grids and works on the principle of mapping the flow domain onto computational domain with simple shape. The mapping is quite tedious if it involves Complex geometry. In order to model this type of geometry we divide the flow region into various smaller sub domains. All these regions are meshed separately and joined up correctly with the neighbors. This type of arrangement is known as block-structured grid. This type of system is more flexible than the previous one. Two-dimensional structured mesh use quadrilaterals elements, while three-dimension meshes use hexahedra. There are two types of body-fitted coordinate grids:

a) Orthogonal curvilinear coordinate.

In orthogonal mesh the grid lines are perpendicular to intersection. This is shown in Figure 2.

b) Non–orthogonal coordinate.

Figure 3 shows non-orthogonal grids. The figure shows the grid lines do not intersect at 90-degree angle. In both these cases the domain boundaries coincide with the coordinate lines; therefore all the geometrical details can be incorporated. Grids can be refined easily to capture important flow features.

Comparison between Cartesian and curvilinear grids

Comparison between Cartesian and curvilinear grids shows that in Cartesian grid cells are wasted in dealing with objects. The distribution of function is very fine in curvilinear grid. The resources required in curvilinear grids are less as compared to Cartesian grids thus saving much memory. Therefore, we can say that coarse grids are able to capture flow details efficiently.

Disadvantages of curvilinear grids

Difficulties associated with the curvilinear grids are related to equations.^[3]

While in Cartesian system the equation can be solved easily with less difficulty, but in curvilinear coordinate system it is difficult to solve the complex equations. Difference between various techniques lies in the fact that what type of grid arrangement is required and the dependent variable that is required in momentum equation. To generate meshes so that it includes all the geometrical features mapping is very important. In mapping physical geometry is mapped with computational geometry.

There are difficulties which we face in generating the body-fitted grids in geometries like IC engine combustion chamber. For example, the valve mapping in internal combustion engine is done very carefully so that the region of one type is mapped carefully with another type of region. There are regions where dense mesh is done deliberately to accommodate complex features. But this results in unnecessary grid resolution which leads to local variation of solution domain.

Block-structured grid

In this type of grid the domain is divided into different regions. Each region has different type of mesh structure. It is also possible that different coordinate system can be used for different regions. This makes the grids far more flexible. This also makes the refinement in the region where the geometry is to be captured more precise. Figure 4 shows the use of block grid technique. The advantages of this technique are that its generation is easy, the equations that we work with are easily discretized and curved boundaries are easily accommodated. Different blocks can be handled with required degree of fineness in mesh. GridPro can create beautiful multiblock structured meshes.

Unstructured grids

In far more complex geometries it is logical to use large number of blocks and therefore it leads to unstructured grids.^[4] These are widely accepted in computational fluid dynamics because this gives us far more flexibility and computer resources are also efficiently utilized. In this case two-dimensional unstructured mesh uses triangle elements while three-dimensional uses tetrahedral elements. These are combination of small structured mesh arranged in unstructured pattern. In this type of grid each single cell is treated as a block. There is no structure of coordinate lines that is given by the grid. The advantage of this type of grids is that mesh can be refined wherever needed. This is based on the fact since the control volume can be of any shape therefore restriction on number of adjacent cell is lifted. This uses various combinations of cell shapes. Hybrid grid can also be used. Hybrid grids are those in which mixture of triangular and quadrilateral elements are used to construct grids. In three dimension combination of tetrahedral and hexahedral elements results in hybrid grid. Example of hybrid grid is shown in Figure 5. Unstructured grids reduces the time which is required for meshing and mapping. Therefore, grid generation is faster and easy. Various automatic techniques especially those associated with Finite Element Method also utilize unstructured grids. Mesh refinement and adaption is easy in unstructured grids.

Mesh improvement

This is done by two techniques namely smoothing and flipping. In mesh smoothing the location of mesh vertices are adjusted. While in Flipping the diagonals of triangulated quadrilateral are exchanged. Flipping improves the quality measure of triangles.

Related Research Articles

Linear algebra is the branch of mathematics concerning linear equations such as:

In geometry, a coordinate system is a system that uses one or more numbers, or coordinates, to uniquely determine the position of the points or other geometric elements on a manifold such as Euclidean space. The order of the coordinates is significant, and they are sometimes identified by their position in an ordered tuple and sometimes by a letter, as in "the x-coordinate". The coordinates are taken to be real numbers in elementary mathematics, but may be complex numbers or elements of a more abstract system such as a commutative ring. The use of a coordinate system allows problems in geometry to be translated into problems about numbers and vice versa; this is the basis of analytic geometry.

Computational fluid dynamics (CFD) is a branch of fluid mechanics that uses numerical analysis and data structures to analyze and solve problems that involve fluid flows. Computers are used to perform the calculations required to simulate the free-stream flow of the fluid, and the interaction of the fluid with surfaces defined by boundary conditions. With high-speed supercomputers, better solutions can be achieved, and are often required to solve the largest and most complex problems. Ongoing research yields software that improves the accuracy and speed of complex simulation scenarios such as transonic or turbulent flows. Initial validation of such software is typically performed using experimental apparatus such as wind tunnels. In addition, previously performed analytical or empirical analysis of a particular problem can be used for comparison. A final validation is often performed using full-scale testing, such as flight tests.

Numerical methods for partial differential equations is the branch of numerical analysis that studies the numerical solution of partial differential equations (PDEs).

In numerical analysis, adaptive mesh refinement (AMR) is a method of adapting the accuracy of a solution within certain sensitive or turbulent regions of simulation, dynamically and during the time the solution is being calculated. When solutions are calculated numerically, they are often limited to predetermined quantified grids as in the Cartesian plane which constitute the computational grid, or 'mesh'. Many problems in numerical analysis, however, do not require a uniform precision in the numerical grids used for graph plotting or computational simulation, and would be better suited if specific areas of graphs which needed precision could be refined in quantification only in the regions requiring the added precision. Adaptive mesh refinement provides such a dynamic programming environment for adapting the precision of the numerical computation based on the requirements of a computation problem in specific areas of multi-dimensional graphs which need precision while leaving the other regions of the multi-dimensional graphs at lower levels of precision and resolution.

Mesh generation is the practice of creating a mesh, a subdivision of a continuous geometric space into discrete geometric and topological cells. Often these cells form a simplicial complex. Usually the cells partition the geometric input domain. Mesh cells are used as discrete local approximations of the larger domain. Meshes are created by computer algorithms, often with human guidance through a GUI, depending on the complexity of the domain and the type of mesh desired. A typical goal is to create a mesh that accurately captures the input domain geometry, with high-quality (well-shaped) cells, and without so many cells as to make subsequent calculations intractable. The mesh should also be fine in areas that are important for the subsequent calculations.

A regular grid is a tessellation of n-dimensional Euclidean space by congruent parallelotopes. Its opposite is irregular grid.

<span class="mw-page-title-main">Computational electromagnetics</span> Branch of physics

Computational electromagnetics (CEM), computational electrodynamics or electromagnetic modeling is the process of modeling the interaction of electromagnetic fields with physical objects and the environment using computers.

In finite element analysis, the spatial twist continuum (STC) is a dual representation of a hexahedral mesh that defines the global connectivity constraint. Generation of an STC can simplify the automated generation of a mesh. The method was published in 1993 by a group led by Peter Murdoch.

Parallel mesh generation in numerical analysis is a new research area between the boundaries of two scientific computing disciplines: computational geometry and parallel computing. Parallel mesh generation methods decompose the original mesh generation problem into smaller subproblems which are solved (meshed) in parallel using multiple processors or threads. The existing parallel mesh generation methods can be classified in terms of two basic attributes:

the sequential technique used for meshing the individual subproblems and
the degree of coupling between the subproblems.

In computational fluid dynamics, the volume of fluid (VOF) method is a free-surface modelling technique, i.e. a numerical technique for tracking and locating the free surface. It belongs to the class of Eulerian methods which are characterized by a mesh that is either stationary or is moving in a certain prescribed manner to accommodate the evolving shape of the interface. As such, VOF is an advection scheme—a numerical recipe that allows the programmer to track the shape and position of the interface, but it is not a standalone flow solving algorithm. The Navier–Stokes equations describing the motion of the flow have to be solved separately. The same applies for all other advection algorithms.

In computational fluid dynamics, the immersed boundary method originally referred to an approach developed by Charles Peskin in 1972 to simulate fluid-structure (fiber) interactions. Treating the coupling of the structure deformations and the fluid flow poses a number of challenging problems for numerical simulations. In the immersed boundary method the fluid is represented in an Eulerian coordinate system and the structure is represented in Lagrangian coordinates. For Newtonian fluids governed by the Navier–Stokes equations, the fluid equations are

The finite element method (FEM) is a popular method for numerically solving differential equations arising in engineering and mathematical modeling. Typical problem areas of interest include the traditional fields of structural analysis, heat transfer, fluid flow, mass transport, and electromagnetic potential.

<span class="mw-page-title-main">TELEMAC</span> Open source software for computational fluid dynamics

In computational fluid dynamics, TELEMAC is short for the open TELEMAC-MASCARET system, or a suite of finite element computer program owned by the Laboratoire National d'Hydraulique et Environnement (LNHE), part of the R&D group of Électricité de France. After many years of commercial distribution, a Consortium was officially created in January 2010 to organize the open source distribution of the open TELEMAC-MASCARET system now available under GPLv3.

A mesh is a representation of a larger geometric domain by smaller discrete cells. Meshes are commonly used to compute solutions of partial differential equations and render computer graphics, and to analyze geographical and cartographic data. A mesh partitions space into elements over which the equations can be solved, which then approximates the solution over the larger domain. Element boundaries may be constrained to lie on internal or external boundaries within a model. Higher-quality (better-shaped) elements have better numerical properties, where what constitutes a "better" element depends on the general governing equations and the particular solution to the model instance.

Gerris is computer software in the field of computational fluid dynamics (CFD). Gerris was released as free and open-source software, subject to the requirements of the GNU General Public License (GPL), version 2 or any later.

Fluid motion is governed by the Navier–Stokes equations, a set of coupled and nonlinear partial differential equations derived from the basic laws of conservation of mass, momentum and energy. The unknowns are usually the flow velocity, the pressure and density and temperature. The analytical solution of this equation is impossible hence scientists resort to laboratory experiments in such situations. The answers delivered are, however, usually qualitatively different since dynamical and geometric similitude are difficult to enforce simultaneously between the lab experiment and the prototype. Furthermore, the design and construction of these experiments can be difficult, particularly for stratified rotating flows. Computational fluid dynamics (CFD) is an additional tool in the arsenal of scientists. In its early days CFD was often controversial, as it involved additional approximation to the governing equations and raised additional (legitimate) issues. Nowadays CFD is an established discipline alongside theoretical and experimental methods. This position is in large part due to the exponential growth of computer power which has allowed us to tackle ever larger and more complex problems.

In physics, a free surface flow is the surface of a fluid flowing that is subjected to both zero perpendicular normal stress and parallel shear stress. This can be the boundary between two homogeneous fluids, like water in an open container and the air in the Earth's atmosphere that form a boundary at the open face of the container.

<span class="mw-page-title-main">FEATool Multiphysics</span>

FEATool Multiphysics is a physics, finite element analysis (FEA), and partial differential equation (PDE) simulation toolbox. FEATool Multiphysics features the ability to model fully coupled heat transfer, fluid dynamics, chemical engineering, structural mechanics, fluid-structure interaction (FSI), electromagnetics, as well as user-defined and custom PDE problems in 1D, 2D (axisymmetry), or 3D, all within a graphical user interface (GUI) or optionally as script files. FEATool has been employed and used in academic research, teaching, and industrial engineering simulation contexts.

References

↑ Patankar, Suhas V. (1980). Numerical Heat Transfer and Fluid FLow. Hemisphere Publishing Corporation. ISBN 0891165223.
↑ Marshall, D., and Ruffin, S.M., " An Embedded Boundary Cartesian Grid Scheme for Viscous Flows using a New Viscous Wall Boundary Condition Treatment,” AIAA Paper 2004-0581, presented at the AIAA 42nd Aerospace Sciences Meeting, January 2004.
↑ Versteeg, H.K.; Malalasekera, W. (2007). An introduction to computational fluid dynamics : the finite volume method (2nd ed.). Harlow: Prentice Hall. ISBN 9780131274983.
↑ Karmna, Steve L. Jr., "Splitfor: A 3D Unstructured Cartesian Prismatic Grid CFD Code for Complex Geometries," AIAA paper 95-0343, presented at the 33rd Aerospace Sciences Meeting and Exhibit, Reno Nevada, January 1995.