Enhancing control of grid distribution in algebraic grid generation

Three techniques are presented to enhance the control of grid-point distribution for a class of algebraic grid generation methods known as the two-, four- and six-boundary methods. First, multidimensional stretching functions are presented, and a technique is devised to construct them based on the desired distribution of grid points along certain boundaries. Second, a normalization procedure is proposed which allows more effective control over orthogonality of grid lines at boundaries and curvature of grid lines near boundaries. And third, interpolating functions based on tension splines are introduced to control curvature of grid lines in the interior of the spatial domain. In addition to these three techniques, consistency conditions are derived which must be satisfied by all user-specified data employed in the grid generation process to control grid-point distribution. The usefulness of the techniques developed in this study was demonstrated by using them in conjunction with the two- and four-boundary methods to generate several grid systems, including a three-dimensional grid system in the coolant passage of a radial turbine blade with serpentine channels and pin fins.     

Boundary conformed co-ordinate systems for selected two-dimensional fluid flow problems. Part I: Generation of BFGs

In this paper the generation of general curvilinear co-ordinate systems for use in selected two-dimensional fluid flow problems is presented. The curvilinear co-ordinate systems are obtained from the numerical solution of a system of Poisson equations. The computational grids obtained by this technique allow for curved grid lines such that the boundary of the solution domain coincides with a grid line. Hence, these meshes are called boundary fitted grids (BFG). The physical solution area is mapped onto a set of connected rectangles in the transformed (computational) plane which form a composite mesh. All numerical calculations are performed in the transformed plane. Since the computational domain is a rectangle and a uniform grid with mesh spacings Δξ = Δη = 1 (in two-dimensions) is used, the computer programming is substantially facilitated. By means of control functions, which form the r.h.s. of the Poisson equations, the clustering of grid lines or grid points is governed. This allows a very fine resolution at certain specified locations and includes adaptive grid generation. The first two sections outline the general features of BFGs, and in section 3 the general transformation rules along with the necessary concepts of differential geometry are given. In section 4 the transformed grid generation equations are derived and control functions are specified. Expressions for grid adaptation arc also presented. Section 5 briefly discusses the numerical solution of the transformed grid generation equations using sucessive overrelaxation and shows a sample calculation where the FAS (full approximation scheme) multigrid technique was employed. In the companion paper (Part II), the application of the BFG method to selected fluid flow problems is addressed.     

Surface grid generation with a linkage to geometric generation

This paper describes a new method to generate surface grids over complex configurations defined by a geometric generation system. The scheme is designed for direct utilization of the surface definition provided by a geometric modeller based on a boundary representation (the so-called B-rep modeller). Thus, the conversion of the geometric representation for the surface grid generator is not required. Consequently, this technique eliminates not only laborious tedium in the conversion of data, but also errors in the representation of the surface induced in the process of the conversion. The proposed method is accomplished over several stages. First, the triangulation is performed on the surface of the geometry, on which the area to be grided is laid. Then linear partial differential equations are mapped and solved on these triangular elements. Finally, the surface grid is constructed by searching for the contours inside the solution domain. After the co-ordinate values of the grid points are obtained by a linear interpolation within each triangular element, these values are mapped onto the surface of the geometry through surface parametric functions provided by the B-rep modeller. An example of generating surface grid over a car configuration is given to illustrate the capability of the method.     

Numerical solution of solidification in a square prism using an algebraic grid generation technique

The solidification of an infinitely long square prism was analyzed numerically. A front fixing technique along with an algebraic grid generation scheme was used, where the finite difference form of the energy equation is solved for the temperature distribution in the solid phase and the solid–liquid interface energy balance is integrated for the new position of the moving solidification front. Results are given for the moving solidification boundary with a circular phase change interface. An algebraic grid generation scheme was developed for two-dimensional domains, which generates grid points separated by equal distances in the physical domain. The current scheme also allows the implementation of a finer grid structure at desired locations in the domain. The method is based on fitting a constant arc length mesh in the two computational directions in the physical domain. The resulting simultaneous, nonlinear algebraic equations for the grid locations are solved using the Newton-Raphson method for a system of equations. The approach is used in a two-dimensional solidification problem, in which the liquid phase is initially at the melting temperature, solved by using a front-fixing approach. The difference of the current study lies in the fact that front fixing is applied to problems, where the solid–liquid interface is curved such that the position of the interface, when expressed in terms of one of the coordinates is a double valued function. This requires a coordinate transformation in both coordinate directions to transform the complex physical solidification domain to a Cartesian, square computational domain. Due to the motion of the solid–liquid interface in time, the computational grid structure is regenerated at every time step.     

A numerical procedure for the generation of orthogonal body-fitted co-ordinate systems with direct determination of grid points on the boundary

The procedure proposed is based on the solution by finite difference means of a set of Laplace's equations, by the application of a relaxation method. The curvilinear orthogonal grid so generated is fitted to a 2-D physical domain with closed boundary and the contribution of the present work consists in the arbitrary choice of grid points on two adjacent boundaries, in order to achieve the desired density of grid points where the geometry of the boundaries varies rapidly. The method proposed is rapid and stable. Some characteristic examples are finally presented.     

A large-scale parallel hybrid grid generation technique for realistic complex geometry

High-Performance Computing (HPC) systems and Computational Fluid Dynamics (CFD) have made significant progress in recent years; however, as the basis of the large-scale parallel computing, the massive grid generation of billions of cells has become a bottleneck problem. In this study, a parallel grid generation technique is proposed to generate large-scale mixed grids with arbitrary cell types and scales. The basic idea of our method is analogous to the global mesh refinement technique. An initial coarse grid with arbitrary cell types is regarded as a background mesh which is partitioned into subzones, and subzones are assigned onto different CPU cores. After the cells and faces in each subzone are split, the inserted new points of the solid wall are projected onto the original CAD entities to preserve the geometry accurately. Finally, the tangled cells caused by the projection in the boundary layer are untangled by a local Radial Basis Function mesh deformation technique. Furthermore, a parallel partition approach and an efficient wall distance computing technique for massive grids are developed also to shorten the preprocessing time. The tests show that the preprocessing efficiency has been increased by two or three orders compared with traditional methods. Billions of grids are generated for the AIAA JSM high-lift model and the Chinese CHN-T1 transport model to test the ability of the parallel grid generation technique. The maximum scale up to 19 billion mixed elements is generated using 16 384 CPU cores in parallel, and the mesh quality is acceptable for CFD simulations.     

An algorithm of finite element grid generation for three-dimensional complex connected domain

An automated quasi three-dimensional finite element grid generation method is presented for particular three-dimensional complex connected domain, across which some are simply connected two-dimensional.regions and some are multiply connected two-dimensional regions. A subdivision algorithm based on the variational principle has been developed to ascertain the smoothness and orthogonality of the generated grid in any cross sections. Smooth transition between the simply and multiply connected regions is maintained. For illustration, the method is applied to generate a finite element three-dimensional grid for human above knee stump.     

Solving Multi-dimensional Evolution Problems with Localized Structures using Second Generation Wavelets

A dynamically adaptive numerical method for solving multi-dimensional evolution problems with localized structures is developed. The method is based on the general class of multi-dimensional second-generation wavelets and is an extension of the second-generation wavelet collocation method of Vasilyev and Bowman to two and higher dimensions and irregular sampling intervals. Wavelet decomposition is used for grid adaptation and interpolation, while O ( N ) hierarchical finite difference scheme, which takes advantage of wavelet multilevel decomposition, is used for derivative calculations. The prowess and computational efficiency of the method are demonstrated for the solution of a number of two-dimensional test problems.     

A method for generating irregular computational grids in multiply connected planar domains

A method for generating irregular triangular computational grids in two-dimensional multiply connected domains is described. A set of points around each body is defined using a simple grid generation technique appropriate to the geometry of each body. The Voronoi regions associated with the resulting global point distribution are constructed from which the Delaunay triangulation of the set of points is thus obtained. The definition of Voronoi regions ensures that the triangulation produces triangles of reasonable aspect ratios given a grid point distribution. The approach readily accommodates local clustering of grid points to facilitate variable resolution of the domain. The technique is generally applicable and has been used with success in computing triangular grids in multiply connected planar domains. The suitability of such grids for flow calculations is demonstrated using a finite element method for solution of the inviscid transonic flow over two- dimensional high-lift aerofoil configurations.     

Hybrid grid generation for viscous flow analysis

Cartesian grid with cut‐cell method has drawn attention of CFD researchers owing to its simplicity. However, it suffers from the accuracy near the boundary of objects especially when applied to viscous flow analysis. Hybrid grid consisting of Cartesian grid in the background, body‐fitted layer near the object, and transition layer connecting the two is an interesting alternative. In this paper, we propose a robust method to generate hybrid grid in two‐dimensional (2D) and three‐dimensional (3D) space for viscous flow analysis. In the first step, body‐fitted layer made of quadrangles (in 2D) or prisms (in 3D) is created near the object's boundary by extruding front nodes with a speed function depending on the minimum normal curvature obtained by quadric surface fitting. To solve global interferences effectively, a level set method is used to find candidates of colliding cells. Then, axis‐aligned Cartesian grid (quadtree in 2D or octree in 3D) is filled in the rest of the domain. Finally, the gap between body‐fitted layer and Cartesian grid is connected by transition layer composed of triangles (in 2D) or tetrahedrons (in 3D). Mesh in transition layer is initially generated by constrained Delaunay triangulation from sampled points based on size function and is further optimized to provide smooth connection. Our approach to automatic hybrid grid generation has been tested with many models including complex geometry and multi‐body cases, showing robust results in reasonable time. Copyright © 2012 John Wiley & Sons, Ltd.     

A two-step godunov-type scheme for the euler equations

A second-order Godunov-type scheme for the Euler equations in conservation form is derived. The method is based on the ENO formulation proposed by Harten et al. The fundamental difference lies in the use of a two-step scheme to compute the time evolution. The scheme is TVD in the linear scalar case, and gives oscillation-free solutions when dealing with nonlinear hyperbolic systems. The admissible time step is twice that of classical Godunovtype schemes. This feature makes it computationally cheaper than one-step schemes, while requiring the same computer storage.

---

Sommario Viene data una nuova estensione al secondo ordine del metodo di Godunov per la soluzione delle equazioni di Eulero in forma conservativa. Il metodo é basato sulla formulazione ENO proposta da Harten et al. La differenza fondamentale consiste nel calcolo dell'evoluzione temporale, ottenuta mediante uno schema a due passi. Questo consente l'uso di un passo di integrazione nel tempo doppio rispetto agli altri schemi alla Godunov ad un solo passo. Il metodo proposto risulta quindi piú efficiente e puó inoltre essere implementato senza alcun aumento dell'occupazione di memoria. Viene dimostrato che lo schema é TVD nel caso lineare, e che fornisce soluzioni prive di oscillazioni spurie nel caso di sistemi non-lineari.

     

An adaptive Cartesian grid generation method for ‘Dirty’ geometry

Traditional structured and unstructured grid generation methods need a ‘water‐tight’ boundary surface grid to start. Therefore, these methods are named boundary to interior (B2I) approaches. Although these methods have achieved great success in fluid flow simulations, the grid generation process can still be very time consuming if ‘non‐water‐tight’ geometries are given. Significant user time can be taken to repair or clean a ‘dirty’ geometry with cracks, overlaps or invalid manifolds before grid generation can take place. In this paper, we advocate a different approach in grid generation, namely the interior to boundary (I2B) approach. With an I2B approach, the computational grid is first generated inside the computational domain. Then this grid is intelligently ‘connected’ to the boundary, and the boundary grid is a result of this ‘connection’. A significant advantage of the I2B approach is that ‘dirty’ geometries can be handled without cleaning or repairing, dramatically reducing grid generation time. An I2B adaptive Cartesian grid generation method is developed in this paper to handle ‘dirty’ geometries without geometry repair. Comparing with a B2I approach, the grid generation time with the I2B approach for a complex automotive engine can be reduced by three orders of magnitude. Copyright © 2002 John Wiley & Sons, Ltd.     

A finite volume method for numerical grid generation

A novel method to generate body‐fitted grids based on the direct solution for three scalar functions is derived. The solution for scalar variables ξ, η and ζ is obtained with a conventional finite volume method based on a physical space formulation. The grid is adapted or re‐zoned to eliminate the residual error between the current solution and the desired solution, by means of an implicit grid‐correction procedure. The scalar variables are re‐mapped and the process is reiterated until convergence is obtained. Calculations are performed for a variety of problems by assuming combined Dirichlet–Neumann and pure Dirichlet boundary conditions involving the use of transcendental control functions, as well as functions designed to effect grid control automatically on the basis of boundary values. The use of dimensional analysis to build stable exponential functions and other control functions is demonstrated. Automatic procedures are implemented: one based on a finite difference approximation to the Cristoffel terms assuming local‐boundary orthogonality, and another designed to procure boundary orthogonality. The performance of the new scheme is shown to be comparable with that of conventional inverse methods when calculations are performed on benchmark problems through the application of point‐by‐point and whole‐field solution schemes. Advantages and disadvantages of the present method are critically appraised. Copyright © 1999 National Research Council of Canada.     

Numerical predictions for unsteady viscous flow past an array of cylinders

The unsteady two-dimensional flow around an array of circular cylinders submerged in a uniform onset flow is analysed. The fluid is taken to be viscous and incompressible. The array of cylinders consists of two horizontal rows extending to infinity in the upstream and downstream directions. The centre-to-centre distance between adjacent cylinders is fixed at three diameters, and the rows are staggered. Advantage is taken of spatially periodic boundary conditions in the flow direction. This reduces the computational domain to a rectangular region surrounding a single circular cylinder. Two cases, for Reynolds numbers of 1000 and 10,000, are presented.     

Orthogonal grid generation for Navier–Stokes computations

Body conforming orthogonal grids were generated using a fast hyperbolic method for aerofoils, and were used to solve the Navier–Stokes equation in the generalized orthogonal system for the first time for time accurate simulation of incompressible flow. For grid generation, the Beltrami equation and the definition equation for the orthogonality are solved using a finite difference method. The grids generated around aerofoils by this method have better orthogonality than the results published by earlier investigators. The Navier–Stokes equation at Reynolds numbers of 3000 and 35 000 for NACA 0012 and NACA 0015 respectively, have been solved as an application. The obtained results match quite well with the corresponding experimental results. © 1998 John Wiley & Sons, Ltd.     

An extended algebraic reconstruction technique (ART) for density-gradient projections: laser speckle photographic tomography

An extended algebraic reconstruction technique (ART) is presented for tomographic image reconstruction from the density-gradient projections, such as laser speckle photography. The essence of the extended ART is that the density-gradient projection data of speckle photography (Eq. (1)) are first numerically integrated to the algebraic representation of interferometric fringe number data (Eq. (12)), which ART can readily reconstruct into the cross-sectional field. The extended ART is numerically examined by using two computer synthesized phantom fields, and experimentally by using asymmetric single and double helium jets in air. The experimentally reconstructed images were also compared with the direct measurements of helium concentration using an oxygen analyzing probe. The extended ART method shows an improved accuracy and is proposed to use to tomographically reconstruct the density-gradient projections over the previous Fourier convolution (FC) method (Liu et al. 1989). Received: 26 June 1998/Accepted: 18 March 1999     

Constitutive model for third harmonic generation in elastic solids

In this article, we present a new constitutive model for studying ultrasonic third harmonic generation in elastic solids. The model is hyperelastic in nature with two parameters characterizing the linear elastic material response and two other parameters characterizing the nonlinear response. The limiting response of the model as the nonlinearity parameters tend to zero is shown to be the well-known St Venant–Kirchhoff model. Also, the symmetric response of the model in tension and compression and its role in third harmonic generation is shown. Numerical simulations are carried out to study third harmonic generation in materials characterized by the proposed constitutive model. Predicted third harmonic guided wave generation reveals an increasing third harmonic content with increasing nonlinearity. On the other hand, the second harmonics are independent of the nonlinearity parameters and are generated due to the geometric nonlinearity. The feasibility of determining the nonlinearity parameters from third harmonic measurements is qualitatively discussed.     

Solution of shallow water equations for complex flow domains via boundary-fitted co-ordinates

Many problems of applied oceanography and environmental science demand the solution of the momentum, mass and energy equations on physical domains having curving coastlines. Finite-difference calculations representing the boundary as a step function may give inaccurate results near the coastline where simulation results are of greatest interest for numerous applications. This suggests the use of methods which are capable of handling the problem of boundary curvature. This paper presents computational results for the shallow water equations on a circular ring of constant depth, employing the concept of boundary fitted grids (BFG) for an accurate representation of the boundary. All calculations are performed on a rectangle in the transformed plane using a mesh with square grid spacing. Comparisons of the simulations of transient normal mode oscillations and analytic solutions are shown, demonstrating that this technique yields accurate results in situations (provided that there is a reasonable choice of grid) involving a curved boundary. The software developed allows application to any two-dimensional area, regardless of the complexity of the geometry. Simulation runs were made with two co-ordinate systems. For the first system, the grid point distribution was obtained from polar co-ordinates. For the second one, grid point positions were calculated numerically, solving Poisson's equation. It was found that small variations in the metric coefficients do not deteriorate the accuracy of the simulation results. Moreover, comparisons of surface elevation and velocity components at grid points near the inner and outer radii obtained from an x?y Cartesian grid model with the BFG simulation were made. The former model produced inacccuracies at grid points near boundaries, and, owing to the large number of mesh points used to yield the necessary fine resolution, the computation time was found to be a factor of three higher.