Fact-checked by Grok 2 weeks ago

Unstructured grid

An unstructured grid, also known as an irregular or unstructured mesh, is a discrete representation of a computational domain in numerical methods such as computational fluid dynamics (CFD) and finite element analysis, composed of irregularly arranged cells—typically triangles in two dimensions or tetrahedra in three dimensions—with arbitrary connectivity between elements.^[1]^[2] This approach contrasts with structured grids, which use regular, ordered arrangements like Cartesian lattices, by requiring explicit storage of neighbor connectivity lists to navigate the mesh.^[2]^[3] Unstructured grids are particularly valued for their ability to conform to complex geometries, such as those in aerospace configurations, turbomachinery, or environmental modeling, where regular grids would require excessive refinement or distortion.^[1]^[3] Key generation techniques include advancing front methods, Delaunay triangulation, and adaptive refinement, which allow local adjustments in resolution to capture features like shocks, wakes, or boundaries with higher fidelity.^[1] Applications span hypersonic flows, rotor simulations, and multiphase fluid dynamics, enabling automated meshing that reduces preparation time from weeks to days for intricate domains.^[1]^[3] Despite their flexibility, unstructured grids present challenges, including higher memory demands due to connectivity data, increased computational complexity in solvers, and potential accuracy reductions from irregular element shapes, such as poor conditioning in tetrahedral cells lacking parallel faces.^[2] Recent advancements integrate machine learning, like graph neural networks and physics-informed neural networks, to handle the irregularity of these grids in modeling complex physics, improving efficiency for high-dimensional problems.^[3] Overall, unstructured grids remain essential for simulating real-world phenomena where geometric adaptability outweighs the added overhead.^[1]^[2]

Fundamentals

Definition

An unstructured grid is a tessellation of a domain in Euclidean space using simple polygonal or polyhedral cells, such as triangles in two dimensions or tetrahedra in three dimensions, arranged irregularly without a predefined topological ordering or uniform spacing.^[1] This irregularity enables the grid to adapt flexibly to boundaries and internal features of complex geometries, where nodes—representing vertices—are distributed non-uniformly, connected by edges to form faces, and ultimately bounded by cells or elements that fill the domain. Unlike structured grids with their predictable, layered connectivity, unstructured grids require explicit storage of these interconnections to navigate the mesh during computations.^[1] The primary purpose of an unstructured grid in numerical simulations is to discretize irregular domains into non-uniform control volumes, facilitating the approximation of partial differential equations (PDEs) that govern physical phenomena such as heat transfer, structural deformation, or fluid flow.^[4] By allowing variable cell sizes and shapes, these grids concentrate resolution in regions of high gradients or interest, such as near curved surfaces or flow discontinuities, while coarsening elsewhere to optimize computational efficiency without sacrificing accuracy in solving the underlying PDEs.^[1] Unstructured grids emerged in the 1970s and 1980s alongside advancements in finite element and finite volume methods, driven by the need to model real-world shapes like aircraft components or biological tissues that defy regular discretization.^[5] A seminal contribution was the 1970 paper by Frederick, Wong, and Edge, which introduced an early automatic triangulation algorithm for two-dimensional structural analysis using irregular triangular elements.^[6] For instance, in two dimensions, an unstructured triangular mesh can conform precisely to a curved boundary, such as an airfoil profile, providing better geometric fidelity than a uniform Cartesian grid that would require extensive refinement or stair-step approximations.^[1]

Comparison with structured grids

Structured grids are characterized by a regular, curvilinear array of cells organized in logical i-j-k indices, allowing implicit identification of neighboring cells through simple arithmetic operations on coordinates.^[7] This topology supports efficient tensor-product operations and predictable data access patterns, making them ideal for uniform domains but challenging for irregular boundaries.^[7] In contrast, unstructured grids feature arbitrary element connectivity without a fixed ordering, necessitating explicit storage of neighbor relationships, which enables adaptable meshing for diverse topologies.^[7] The primary differences arise in geometry handling and computational efficiency: structured grids provide uniformity that enhances cache performance and simplifies implementation for simple shapes, but they often require multi-block decompositions or coordinate transformations to fit complex geometries, potentially leading to distortions.^[8] Unstructured grids, commonly using triangles or tetrahedra, avoid such transformations by conforming directly to intricate surfaces, offering greater generality for problems like airfoils or multi-element wings.^[9] For instance, a structured O-grid around an airfoil uses concentric blocks for boundary layer resolution but may introduce singularities at junctions, whereas an unstructured triangular mesh can vary element density smoothly without blocking.^[7] Trade-offs highlight the balance between flexibility and overhead: unstructured grids simplify generation and adaptation for complex domains, reducing user effort in mesh creation, but they demand 2–3 times more memory per grid point for connectivity data and incur higher computational costs due to irregular access patterns.^[8] Structured grids excel in interpolation accuracy and solver efficiency for viscous flows but falter at boundaries, where grid quality degrades.^[7] Hybrid methods, such as Chimera or overset grids, mitigate these by overlaying structured blocks on unstructured backgrounds, combining the former's efficiency with the latter's adaptability for moving or multifaceted geometries.^[10]

Data structures

Node and element representation

In unstructured grids, nodes represent discrete points in the computational domain, typically stored as arrays of floating-point coordinates (x, y, z) in three dimensions or (x, y) in two dimensions, distinguishing between interior nodes within the domain and boundary nodes on the edges or surfaces.^[11] These coordinates enable the geometric positioning essential for defining the grid's spatial structure, with interior nodes often numbering in the millions for complex simulations to capture fine details.^[12] Elements, or cells, in unstructured grids are polygonal or polyhedral volumes that approximate the continuous domain, defined solely by indices referencing the connected nodes rather than explicit position data, allowing flexible topologies without regular patterning.^[11] In two dimensions, common elements include triangles formed by three node indices, while three-dimensional grids feature tetrahedra bounded by four nodes; more general convex elements such as quadrilaterals, prisms, or hexahedra may also be used, defined by their respective node counts.^[13] Simplicial elements like triangles and tetrahedra are preferred for their simplicity in generation and guaranteed convexity, though they often require more elements to achieve comparable approximation accuracy to convex polyhedra, which can provide better resolution with fewer cells in structured-like regions.^[13] Geometric properties of elements are computed from node coordinates to support numerical integration and analysis. The centroid of an element is the arithmetic mean of its vertex coordinates, serving as a reference point for interpolation.^[11] Element volumes are derived using determinant-based formulas; for a tetrahedron with vertices a, b, c, d, the volume V is given by

V = \frac{1}{6} \left| \det \begin{pmatrix} b_x - a_x & c_x - a_x & d_x - a_x \\ b_y - a_y & c_y - a_y & d_y - a_y \\ b_z - a_z & c_z - a_z & d_z - a_z \end{pmatrix} \right|,

or equivalently via the scalar triple product \frac{1}{6} |(b - a) \cdot ((c - a) \times (d - a))|, ensuring positive orientation for consistent integration.^[14] Jacobians, representing the local scaling of the mapping from reference to physical element space, are computed as the determinant of the derivative matrix of nodal shape functions, crucial for transforming integrals over reference elements to the actual geometry.^[14] Boundary representation in unstructured grids often employs dedicated structures to manage domain edges without altering interior connectivity. Boundary faces, shared by one interior element, define the surface, while ghost cells—fictitious elements adjacent to boundaries—extend the grid outward to enforce conditions like no-slip walls or inflows, populated via mirroring or extrapolation from real cells.^[15] This approach maintains the unstructured nature while facilitating accurate flux computations at interfaces.^[12]

Connectivity and storage

In unstructured grids, connectivity data primarily consists of element-to-node mappings, where each element records the indices of its constituent nodes to define its topology and geometry. For instance, a triangular element in a 2D mesh stores three node indices, while a tetrahedral element in 3D stores four, enabling the reconstruction of element shapes from global node coordinates. Complementing this, node-to-element adjacency lists are commonly employed to track all elements incident to a given node, supporting operations like gradient reconstruction or refinement.^[16]^[1] Storage formats for connectivity emphasize efficiency in memory and access speed, often leveraging sparse representations due to the irregular nature of unstructured grids. The compressed sparse row (CSR) format is widely used for storing face or element adjacencies, consisting of three arrays: one for non-zero values (e.g., neighbor indices), one for column indices, and one for row pointers that indicate the start of each row's data, thus compactly encoding the sparse graph of cell connections. Half-edge data structures offer a topology-focused alternative, representing the mesh through directed half-edges that link vertices, edges, and faces via references to twin, next, incident vertex, and incident face entities, facilitating rapid traversal without redundant storage. These formats prioritize contiguous memory layouts to minimize indirection overhead and improve cache locality during iterative solver operations.^[17]^[18] Boundary connectivity is managed through explicit pairings of boundary faces to their adjacent elements, typically including owner and neighbor identifiers to distinguish internal fluxes from boundary conditions, with boundary faces assigned a fictitious exterior neighbor. Orientation conventions, such as the right-hand rule—where curling the fingers of the right hand along the ordered node indices of a face points the thumb in the direction of the outward normal—ensure consistent normal vector directions for flux calculations across the domain.^[19]^[1] Access patterns in computational algorithms rely on these structures to efficiently iterate over local neighborhoods, such as traversing adjacency lists in CSR format to gather neighbor elements for finite difference stencils or flux assembly in solvers. Dual graphs, where cells correspond to nodes and shared faces to edges, enable Voronoi-like dual representations for cell-centered methods, allowing quick identification of adjacent cells without full geometric queries. In a finite volume method, for example, face areas and normals are stored as parallel vectors indexed to the connectivity arrays, permitting vectorized computation of interface fluxes between cells while respecting boundary pairings.^[17]^[1]

Generation techniques

Triangulation-based methods

Triangulation-based methods for generating unstructured grids primarily rely on Delaunay triangulation, which produces high-quality meshes by optimizing geometric properties of the resulting simplices. In two dimensions, a Delaunay triangulation of a set of points is a triangulation where the circumcircle of any triangle contains no other points from the set in its interior, a property that maximizes the minimum angle among all possible triangulations and thereby improves element quality for numerical simulations.^[13] This empty-circle criterion ensures well-shaped triangles, reducing aspect ratios and enhancing interpolation accuracy in applications like finite element methods.^[20] The incremental insertion algorithm, notably the Bowyer-Watson method, constructs the Delaunay triangulation by successively adding points to an initial triangulation and updating the mesh through cavity removal and retriangulation. In this approach, inserting a new point identifies the cavity as the union of triangles whose circumcircles contain the point, deletes those triangles, and reconnects the cavity boundary to the new point while preserving the Delaunay property. To extend this to higher dimensions or ensure conformity, Lawson flips are applied post-insertion: these edge-swapping operations locally optimize the triangulation by replacing diagonals that violate the empty-circle property until global Delaunay satisfaction is achieved. For robustness against floating-point errors in geometric tests like orientation and incircle predicates, implementations incorporate exact arithmetic via adaptive precision, as developed by Shewchuk, enabling reliable computation even for degenerate configurations.^[21] In three dimensions, Delaunay tetrahedralization faces additional challenges like sliver elements—poorly shaped tetrahedra with small dihedral angles. Refinement algorithms, such as extensions of Ruppert's method, address this by inserting Steiner points (e.g., at circumcenters or midpoints) to improve quality by bounding small angles and avoiding degenerate elements.^[22] Divide-and-conquer algorithms, which recursively partition point sets and merge triangulations, offer efficient alternatives to incremental methods, achieving optimal O(n log n) time complexity; Shewchuk's implementations exemplify this with robust predicates to handle numerical precision.^[23] To incorporate boundaries and constraints, constrained Delaunay triangulation enforces prescribed edges, such as those defining a domain outline, by inserting them without fully violating the empty-circle property, often through iterative flips or point perturbations that maintain proximity to the ideal Delaunay while ensuring conformity. This is crucial for meshing complex geometries, where boundary segments must be preserved exactly. Quality improvements in these methods often involve refinement techniques like longest-edge bisection, which iteratively bisects the longest edge of elements to propagate refinement and avoid degenerate slivers in 3D; in tetrahedral meshes, this can be combined with 8-tetrahedra partitions to ensure non-degeneracy and uniform angle improvement across levels. A practical example is generating a triangular mesh for an airfoil profile: boundary points along the NACA 0012 contour are first connected via constrained Delaunay triangulation to form the outline, followed by incremental insertion of interior points distributed based on a sizing function, yielding a graded mesh with refined elements near the leading edge for capturing boundary layer effects.^[24]

Front-based and other algorithms

The advancing front method is a boundary-driven technique for generating unstructured grids, particularly effective for creating layered meshes near complex surfaces. It begins with the discretization of the domain boundary into an initial "front" consisting of edges (in 2D) or faces (in 3D), from which new nodes and elements are progressively created and propagated inward until the domain is filled.^[25] As the front advances, element creation prioritizes ideal shapes (e.g., equilateral triangles or regular tetrahedra) by computing new node positions based on local front geometry, spacing criteria, and smoothness controls to minimize distortion and ensure boundary conformity.^[26] This method, first formulated for unstructured triangular mesh generation in 1971, excels in preserving boundary integrity and generating high-aspect-ratio cells suitable for viscous flows.^[25] Variants of the advancing front method incorporate local remeshing to maintain front quality during propagation, where intersecting or degenerate elements are detected and repaired through node insertion or edge swapping. Collision detection algorithms prevent overlaps by monitoring front advancement and adjusting node placement, often using geometric predicates to resolve intersections in real-time. These enhancements address early limitations in handling concave boundaries and ensure robustness for 3D tetrahedral grids.^[27] Other algorithms complement the advancing front approach, including octree-based subdivision, which starts with a cubic domain decomposition into an octree hierarchy and refines selectively to produce tetrahedral elements by splitting cubes into primitives while adapting to geometry and resolution requirements. Frontal-Delaunay hybrids combine the boundary propagation of advancing fronts with global Delaunay triangulation for interior filling, where the front provides initial points and the Delaunay step optimizes connectivity to avoid poor elements.^[28] These methods offer alternatives for hierarchical or hybrid meshing in complex geometries.^[29] Recent advancements include machine learning-based techniques for unstructured mesh generation, such as implicit geometry neural networks that predict mesh size functions and automate node placement for complex domains, improving efficiency and adaptability as of 2024.^[30] These approaches leverage neural networks to handle irregularity and reduce manual intervention in meshing workflows. For large-scale applications, parallel generation employs domain decomposition to distribute the front across processors, with space-filling curves (e.g., Hilbert curves) used for load balancing by recursively partitioning the domain into contiguous subregions that minimize inter-processor communication.^[31] This approach enables efficient scaling on distributed-memory systems, as demonstrated in consistent parallel advancing front techniques that synchronize fronts at subdomain interfaces.^[32] An illustrative example is the meshing of a turbine blade, where the advancing front starts from the blade's surface contours to extrude prismatic layers inward, capturing boundary layer resolution with high-aspect-ratio tetrahedra near the walls while filling the interior with isotropic elements.^[26] This produces grids suitable for computational fluid dynamics simulations of rotating machinery.^[33]

Applications

Computational fluid dynamics

Unstructured grids play a pivotal role in computational fluid dynamics (CFD) by enabling the discretization of the Navier-Stokes equations using finite volume methods on complex geometries. These grids facilitate the integration of hybrid meshes, which combine prismatic elements near solid walls to resolve boundary layers with tetrahedral elements in the bulk flow region, thereby improving accuracy in capturing viscous effects without excessive computational cost.^[34] This approach is particularly effective for aerospace applications involving irregular shapes, such as aircraft fuselages or turbine blades, where structured grids would require excessive manual intervention.^[12] Solver adaptations for unstructured grids often employ edge-based schemes to compute fluxes efficiently across irregular stencils. In these schemes, fluxes are evaluated along grid edges rather than cell faces, reducing memory requirements and enabling scalable implementations for large-scale simulations of the Euler or Navier-Stokes equations.^[35] For convergence acceleration, agglomeration multigrid methods coarsen the grid by merging elements into larger control volumes, accommodating the irregular connectivity of unstructured meshes while preserving solution smoothness.^[36] These techniques have demonstrated robust performance in steady-state viscous flow predictions.^[36] Specific techniques leveraging unstructured grids include unsteady simulations via the arbitrary Lagrangian-Eulerian (ALE) formulation, which allows mesh motion to track deforming boundaries while maintaining grid quality. ALE methods on unstructured grids have been applied to simulate fluid-structure interactions, such as flapping wings or deploying control surfaces, by combining Lagrangian motion near interfaces with Eulerian fixed grids in the far field.^[37] Additionally, overset unstructured grids enable the modeling of rotor-stator interactions in turbomachinery by overlapping independent meshes around rotating and stationary components, facilitating relative motion without grid regeneration.^[38] This approach captures unsteady wake interactions with high fidelity, as validated in compressor stage simulations.^[39] A notable case study arises from the AIAA Drag Prediction Workshops, where unstructured grids with anisotropic refinement have improved drag predictions for transonic airfoils like the NASA Common Research Model. Anisotropic adaptation, focusing refinement along flow gradients, reduced discretization errors in near-wall regions, yielding drag coefficients consistent with experimental data across multiple workshop cycles.^[40] These results highlight the efficacy of unstructured methods in high-fidelity aerodynamic design.^[41] Unstructured grids address key challenges in aerospace CFD, such as capturing shocks and vortices in hypersonic flows. Shock-capturing techniques on tetrahedral meshes resolve bow shocks around reentry vehicles with minimal dissipation, while shock-fitting methods align the grid to the shock surface for sharper resolution.^[42] In vortex-dominated flows, like those in high-speed inlets, hybrid unstructured grids enhance vorticity preservation, enabling accurate prediction of separation and mixing phenomena.^[43] These capabilities have been demonstrated in simulations of Mach 5+ flows, where unstructured approaches outperform structured grids in handling multi-scale features.^[44]

Finite element analysis

Unstructured grids play a crucial role in finite element analysis (FEA) for solving partial differential equations in structural and multiphysics problems, particularly those involving complex geometries like elasticity and heat transfer. The Galerkin method is applied over unstructured elements, such as tetrahedra or prisms, to approximate solutions by projecting the weak form of the governing equations onto a finite-dimensional subspace spanned by basis functions associated with mesh nodes.^[45] This approach enables flexible discretization of irregular domains without the rigidity of structured grids, making it suitable for problems in solid mechanics where boundary conformity is essential. For curved elements, isoparametric mappings transform the reference element coordinates to the physical domain, ensuring accurate representation of geometry while maintaining computational efficiency in integration.^[46] The assembly process in FEA with unstructured grids begins with local computation of element stiffness matrices for each element in the mesh. For linear elasticity, the element stiffness matrix \mathbf{K}^e is obtained through numerical quadrature of the integral

\mathbf{K}^e = \int_{V^e} \mathbf{B}^T \mathbf{D} \mathbf{B} \, dV,

where \mathbf{B} is the strain-displacement matrix, \mathbf{D} is the material constitutive matrix, and V^e is the element volume; this integral is evaluated using the Jacobian from isoparametric mappings. These local matrices are then scattered into the global stiffness matrix \mathbf{K} using the element-node connectivity data, which maps local degrees of freedom to global ones, ensuring efficient sparse storage and solution of the resulting system \mathbf{K} \mathbf{u} = \mathbf{f}. This connectivity-driven assembly is particularly advantageous for unstructured grids, as it handles arbitrary topologies without predefined ordering.^[47] Higher-order elements enhance accuracy in FEA applications requiring precise stress fields, such as fracture or contact problems. Quadratic tetrahedral elements (TET10), employing Lagrange polynomials as shape functions, capture quadratic variations in displacements and thus provide superior resolution of stress gradients compared to linear elements, reducing locking effects in nearly incompressible materials.^[48] These elements are integrated via Gauss quadrature on the reference tetrahedron, with shape functions defined hierarchically to maintain C^0 continuity across element boundaries. In multiphysics contexts, unstructured grids enable seamless coupling of FEA with computational fluid dynamics (CFD) on shared meshes for fluid-structure interaction (FSI) simulations, where structural deformations influence fluid loads and vice versa. This shared discretization avoids interface projection errors, allowing monolithic or partitioned solvers to enforce continuity at fluid-solid interfaces.^[49] For example, in vehicle crash simulations, adaptive tetrahedral meshes refine resolution at impact zones to accurately model large deformations and stress concentrations in components like frames, improving predictions of energy absorption and failure modes.^[50]

Advantages and challenges

Benefits

Unstructured grids offer substantial geometric flexibility, allowing them to conform closely to irregular boundaries and complex geometries without requiring multi-block decomposition, which streamlines the preprocessing phase and reduces time for CAD model imports.^[51] This adaptability is particularly valuable in industrial applications where geometries are often non-uniform, enabling efficient meshing of real-world domains like aircraft fuselages or turbine blades.^[52] A key advantage lies in their support for adaptive refinement, including local h-adaptation (mesh size adjustment) or p-adaptation (polynomial order variation), which targets high-gradient regions such as boundary layers to enhance solution accuracy while minimizing computational overhead from unnecessary global refinements.^[52] Such techniques control discretization errors effectively in computational fluid dynamics (CFD) and finite element analysis (FEA), leading to higher-fidelity simulations on irregular domains.^[53] Unstructured grids also excel in parallel scalability, as their irregular connectivity lends itself to graph partitioning algorithms like those in the METIS library, which enable balanced load distribution across processors for efficient computation on large-scale irregular domains.^[54] This facilitates better performance in distributed computing environments compared to structured grids on non-uniform geometries. Hybrid meshing strategies integrate unstructured grids with structured boundary layers, providing high-resolution near viscous walls while maintaining flexibility elsewhere, thus optimizing accuracy and efficiency in flows with strong gradients.^[55] In CFD benchmarks on complex shapes, unstructured grids have shown quantitative benefits, including overall time savings through reduced meshing efforts compared to fully structured approaches in practical workflows.^[56]

Limitations and quality metrics

Unstructured grids, while flexible for complex geometries, impose significant computational limitations compared to structured grids. The primary drawback stems from indirect addressing in data structures, which requires additional integer operations to access neighboring elements, leading to higher CPU overhead during solver iterations. This indirect indexing can result in substantially larger memory requirements, as connectivity information must be explicitly stored for each element, increasing storage needs by factors that depend on mesh complexity but are generally higher than the compact array-based representations of structured grids. Furthermore, the irregular connectivity complicates debugging of solver instabilities, as tracing error propagation across non-uniform elements is more challenging than in predictable structured layouts. To assess and ensure the reliability of unstructured grids, several quality metrics are employed to evaluate element shape, size, and orientation. The aspect ratio measures the elongation of an element by comparing the maximum to minimum edge lengths, with values ideally close to 1; ratios exceeding 10:1 often lead to numerical instability in solvers due to ill-conditioned matrices. Skewness quantifies deviation from an equilateral or equiangular ideal shape, typically defined via the maximum cosine of angles between edges or faces, where high skewness indicates distortion that can amplify discretization errors. The Jacobian determinant, derived from the transformation matrix mapping reference to physical elements, ensures invertibility and positive volume; a determinant near zero signals inverted or degenerate elements, rendering the mesh invalid for analysis. Strategies to mitigate poor quality include post-generation improvement techniques such as smoothing and untangling. Laplacian smoothing repositions interior nodes to the average position of their neighbors, iteratively improving element shapes while preserving boundary conditions; adaptive variants apply changes only if quality metrics like aspect ratio enhance, avoiding inversion in concave regions and achieving up to 23-fold speedup on parallel hardware.^[57] For inverted elements, optimization-based untangling maximizes the minimum signed area or volume of local simplices through linear programming, ensuring convex level sets for convergence and effectively resolving tangles without connectivity changes. These methods can be combined, with untangling preceding smoothing to first guarantee validity before optimizing quality. Error estimation in unstructured grids relies on a posteriori metrics to guide adaptive refinement, focusing on solution residuals or inter-element gradient discontinuities. Residual-based estimators compute local imbalances in the governing equations, while gradient recovery techniques average discontinuous gradients across element interfaces to detect jumps, providing both upper and lower bounds on the true error with effectivity indices that approach unity for well-adapted meshes. These indicators enable targeted adaptation, concentrating resolution where errors are high without global remeshing. A notable challenge in three-dimensional tetrahedral meshes is the persistence of slivers—thin, needle-like elements with small dihedral angles that degrade conditioning. Sliver removal often employs variational geometry optimization, such as gradient-based shape matching energy that penalizes low-height tetrahedra by minimizing deviations from a regular template simplex, using Newton's method for vertex updates and edge flips for connectivity. This approach suppresses slivers by enforcing uniform face areas relative to volume, yielding robust meshes suitable for high-fidelity simulations.

References

[1]
[PDF] Unstructured Grid Generation
Unstructured grid generation in Computational. Fluid Dynamics. (CFD) is the discrete representation of flow domains by irregular arrangements of points.
[2]
Grid Systems | What Are the Advantages and Disadvantages?
In unstructured grids, volume elements can be joined in any manner, and special lists must be kept to identify neighboring elements. Structured Meshing of ...
[3]
Machine learning for modelling unstructured grid data in ...
Unstructured grid data are essential for modelling complex geometries and dynamics in computational physics. Yet, their inherent irregularity presents ...
[4]
Numerical Grids - Computational Fluid Dynamics Literature - CCC
Numerical grid generation plays a critical role in any scientific computing problem when the geometry of the underlying region is complex or when the solution ...<|separator|>
[5]
[PDF] Introduction - CS@Purdue
Automatic unstructured mesh generation for finite element methods began in 1970 with an article by C. O. Frederick, Y. C. Wong, and F. W. Edge entitled “Two ...Missing: origin | Show results with:origin
[6]
Two‐dimensional automatic mesh generation for structural analysis
Aug 6, 2025 · 1 The gift wrapping algorithm was proposed in 1970 to compute 2D Delaunay triangulation by Frederick, Wong, and Edge [31] and by Chand and Kapur ...
[7]
[PDF] Geometry Modeling, Flow Domain Modeling, and Grid Generation
Unstructured grids allow greater flexibility in generating and adapting grids at the expense of greater storage of cell information. A Hybrid Grid contains both ...Missing: differences | Show results with:differences
[8]
[PDF] 1 A MULTIBLOCK APPROACH FOR CALCULATING ... - NASA
Although the unstructured grid approach enjoys these advantages over structured grids, flow solvers that utilize it suffer several disadvantages. These ...
[9]
An ANN-based advancing double-front method for automatic ... - NIH
After several decades of effort by CFD researchers, unstructured grid generation techniques have made great progress. However, as we know, all the ...
[10]
[PDF] C_..-3
An overset grid approach is really an unstructured collection of over- lapping structured grid components. As such, the approach should enjoy most of the grid ...
[11]
[PDF] Meshing for the Finite Element Method - Florida State University
The finite element method begins by looking at a complicated region, and thinking of it as a mesh of smaller, simpler subregions. The subregions are simple, ...
[12]
[PDF] UNSTRUCTURED GRID METHODS FOR THE SIMULATION OF 3D ...
... flow algorithms, the adoption of the edge based representation for the unstructured grid has important, computational implications. This is illustrated in ...
[13]
[PDF] Chapter 10 - Unstructured Mesh Generation - People @EECS
One of the central tools of scientific computing is the fifty-year old finite element method—a numerical method for approximating solutions to partial.
[14]
[PDF] FEM BASIS FUNCTIONS FOR A TETRAHEDRON
Dec 23, 2010 · Exercise: Show that if the tetrahdron T=1a,b,c. d1 is actually the reference tetrahedron, then the determinant formulas for the basis functions ...
[15]
(PDF) Parallel Multi-Level Ghost Cell Generation for Distributed ...
Oct 28, 2017 · Ghost cells are important for distributed memory parallel operations that require neighborhood information, and are required for correctness ...
[16]
vtkUnstructuredGrid Class Reference - VTK
vtkUnstructuredGrid Class Reference dataset represents arbitrary combinations of all possible cell types.
[17]
Data Structure for Unstructured Meshes with GPU Support
A general multi-purpose data structure for an efficient representation of conforming unstructured homogeneous meshes for scientific computations on CPU and GPU ...Missing: half- | Show results with:half-
[18]
https://dl.acm.org/doi/pdf/10.1145/2614348.2614376
[19]
OpenFOAM v7 User Guide - 5.1 Mesh description - CFD Direct
Jul 10, 2019 · The direction of the face normal vector is defined by the right-hand rule, i.e. looking towards a face, if the numbering of the points follows ...Missing: unstructured | Show results with:unstructured<|control11|><|separator|>
[20]
[PDF] Delaunay Refinement Mesh Generation
May 18, 1997 · Delaunay refinement generates unstructured meshes by maintaining a Delaunay triangulation or tetrahedralization, refining it by inserting ...
[21]
[PDF] Robust Adaptive Floating-Point Geometric Predicates
The incircle test determines whether a point lies inside, outside, or on a circle or sphere, and is used for Delaunay triangulation [8]. Inexact versions of ...
[22]
[PDF] A Delaunay Re nement Algorithm for Quality 2-Dimensional Mesh ...
The algorithm triangulates polygons, ensuring bounded aspect ratio triangles and a number of triangles within a constant factor of optimal, using successive ...
[23]
Triangle: A Two-Dimensional Quality Mesh Generator and Delaunay ...
Triangle generates exact Delaunay triangulations, constrained Delaunay triangulations, conforming Delaunay triangulations, Voronoi diagrams, and high-quality ...Delaunay triangulations · Demonstration · Command line switches
[24]
TRIANGLE - Mesh generation and Delaunay triangulation
Oct 9, 2014 · TRIANGLE is a C program which generates meshes, Delaunay triangulations and Voronoi diagrams for 2D pointsets, by Jonathan Shewchuk.
[25]
Advancing Front Technique - an overview | ScienceDirect Topics
The advancing front technique, first formulated by George (1971), is an unstructured grid generation method which preserves boundary integrity.Missing: seminal | Show results with:seminal
[26]
[PDF] Unstructured Viscous Grid Generation by Advancing-Front Method
A new method of generating unstructured triangular/tetrahedral grids with high-aspect-ratio. ceUs is proposed. The method is based on a new grid-marching.Missing: seminal | Show results with:seminal
[27]
3D Delaunay mesh generation coupled with an advancing-front ...
This paper describes a fully automatic mesh generation method suitable for domains of any shape in R 3. Initially, a so-called boundary mesh is generated ...Missing: hybrid | Show results with:hybrid
[28]
[PDF] An Algorithm of Combining of Advancing Front and Delaunay ...
In this paper an algorithm of unstructured grid generation is presented. It is based on coupling of the Delaunay triangulation and Advancing Front Tech- nique.
[29]
[PDF] An Advancing Front Delaunay Triangulation Algorithm ... - DTIC
'Ifv. A new algorithm is described for generating an unstructured mesh about ar -tT'A twU- dimensional configuration. Mesh points are generated ...<|separator|>
[30]
Parallelizing an Unstructured Grid Generator with a Space-Filling ...
Aug 18, 2000 · This new approach is based on a space- filling curve. The space-filling curve's indices are calculated recursively and in parallel, thus leading ...
[31]
CPAFT: A Consistent Parallel Advancing Front Technique for ... - arXiv
May 31, 2024 · In this paper, we propose a novel, consistent parallel advancing front technique (CPAFT) by combining the advancing front technique, the domain decomposition ...<|separator|>
[32]
[PDF] Development of an Automatic TET10 Meshing program for rotating ...
Mar 5, 2020 · This thesis addresses mesh generation with the focus on automation. It targets the meshing of a complex turbine blade and disk assembly.
[33]
Unstructured Grid Approaches for Accurate Aeroheating Simulations
Jun 28, 2007 · The use of tetrahedral, prismatic, and hybrid hexahedral-prismatic-tetrahedral grids for the accurate prediction of aerodynamic heating at ...
[34]
An efficient implementation of edge-based discretization without ...
Aug 6, 2025 · The node-centered edge-based discretization method is widely used today in practical unstructured-grid computational fluid dynamics (CFD) ...
[35]
[PDF] On Convergence Acceleration Techniques for Unstructured Meshes
Abstract. A discussion of convergence acceleration techniques as they relate to computational fluid dynamics problems on unstructured meshes is given.
[36]
Multigrid solution of the two-dimensional Euler equations on ...
An aspect ratio agglomeration multigrid for unstructured grids. 16 January 2013 | International Journal for Numerical Methods in Fluids, Vol. 72, No. 10.
[37]
An unstructured mesh arbitrary Lagrangian-Eulerian unsteady ...
Jun 1, 2016 · In this paper, an unstructured mesh Arbitrary Lagrangian-Eulerian (ALE) incompressible flow solver is developed to investigate the ...
[38]
Numerical analysis of a high-order unstructured overset grid method ...
Jun 15, 2018 · Evidence is presented which shows that significant influence of the rotor-stator interaction is captured, and the pressure fluctuation ...
[39]
[PDF] Numerical analysis of a high-order unstructured overset grid ... - HAL
May 20, 2020 · overset grid method is intended to be used at each rotor-stator interface in a multi-row turbomachine. Its range of applications is ...
[40]
[PDF] Grid Quality and Resolution Issues from the Drag Prediction ...
Dec 5, 2006 · The drag prediction workshop series (DPW), held over the last six years, and sponsored by the AIAA Applied Aerodynamics Committee, ...
[41]
Summary Data from the Sixth AIAA CFD Drag Prediction Workshop
The focus of the workshop was to compare absolute drag predictions, including the variation due to grid type and turbulence model type.
[42]
Hypersonic Flow Computations on Unstructured Grids
Nov 13, 2012 · Hypersonic Flow Computations on Unstructured Grids: Shock-Capturing vs. Shock-Fitting Approach. Aldo Bonfiglioli and; Renato Paciorri.
[43]
A shock-fitting technique for 2D unstructured grids - ScienceDirect.com
This paper illustrates the algorithmic features of this original technique and the results obtained in the computation of the hypersonic flow past a circular ...
[44]
An Unstructured, Three-Dimensional, Shock-Fitting Solver for ...
Aug 10, 2025 · A novel unstructured shock-fitting algorithm for three-dimensional flows is presented in this paper. The fitted shock front is described ...
[45]
Elastic full waveform inversion on unstructured meshes by the finite ...
In this paper, we develop a new numerical method for elastic full waveform inversion (FWI) on unstructured triangular meshes.
[46]
[PDF] An Isoparametric Tangled Finite Element Method for Handling ...
This paper presents an isoparametric tangled finite element method (i-TFEM) method for handling tangled high order/curvilinear meshes.
[47]
CMES | Matrix-Free Higher-Order Finite Element Method for Parallel ...
Matrix-Free Higher-Order Finite Element Method for Parallel Simulation of Compressible and Nearly-Incompressible Linear Elasticity on Unstructured Meshes.
[48]
[PDF] Finite Element Algorithms and Data Structures on Graphical ...
The layout of stiffness matrix data is one of the most important factors that determine the speed of the finite element method on the GPU: in the assembly phase ...
[49]
Finite Element Simulation of Articular Contact Mechanics with ...
In this study we implemented both ten-node (TET10) and fifteen-node (TET15) quadratic tetrahedral elements in FEBio (www.febio.org) and compared their accuracy, ...
[50]
Advances on finite element discretization of fluid–structure ...
Sep 30, 2025 · We review the main features of an unfitted finite element method for interface and fluid–structure interaction problems based on a ...<|separator|>
[51]
(PDF) Finite Element Analysis of Car Frame Frontal Crash using ...
Aug 8, 2025 · A mesh convergence study was carried out using the Adaptive mesh refinement technique. This technique makes refinement in the critical area of ...
[52]
Unstructured Grid - an overview | ScienceDirect Topics
Unstructured grids refer to a type of grid generation in which the pattern of connections between grid points can vary from point to point, allowing for ...
[53]
[PDF] Unstructured Grid Adaptation: Status, Potential Impacts, and ...
Unstructured grid adaptation is a powerful tool to control discretization error for Com- putational Fluid Dynamics (CFD). It has enabled key increases in ...
[54]
First benchmark of the Unstructured Grid Adaptation Working Group
Unstructured grid adaptation is a technology that holds the potential to improve the automation and accuracy of computational fluid dynamics and other ...
[55]
METIS—A Software Package for Partitioning Unstructured Graphs ...
A software package for partitioning unstructured graphs, partitioning meshes, and computing fill-reducing orderings of sparse matrices.
[56]
[PDF] A hybrid structured / unstructured grid strategy for the CFD ... - HAL
Dec 20, 2019 · Structured zones may be kept for the sake of efficiency and of accuracy in viscous layers, whereas unstructured zones may enable an easier mesh ...
[57]
[PDF] development of fast unstructured cfd code “fastar”
shows four times faster convergence with global ... [7] Nakahashi, K., et al., “Some challenges of realistic flow simulations by unstructured grid CFD”,.