Journal
COMPUTER AIDED GEOMETRIC DESIGN
Volume 65, Issue -, Pages 48-75Publisher
ELSEVIER SCIENCE BV
DOI: 10.1016/j.cagd.2018.03.023
Keywords
-
Funding
- European Union's Horizon 2020 research and innovation program [678727]
- H2020 Societal Challenges Programme [678727] Funding Source: H2020 Societal Challenges Programme
Ask authors/readers for more resources
The generation of an analysis-suitable computational grid from a description of no more than its boundaries is a common problem in numerical analysis. Most classical meshing techniques for finite-volume, finite-difference or finite-element applications such as the Advancing Front Method (Schoberl, 1997), Delaunay Triangulation (Triangle, 1996) and elliptic or hyperbolic meshing schemes (Thompson et al., 1998) operate with linear or multi-linear but straight-sided elements for the generation of structured and unstructured meshes, respectively, whereas the generation of high-quality curved meshes is still considered a major challenge. A recent development is the introduction of Isogeometric Analysis (IgA) (Hughes et al., 2005), which can be considered as a natural high-order generalisation of the finite-element method. A description of the geometry (Omega) over cap is accomplished via a mapping operator x : (Omega) over cap -> Omega that maps the unit hypercube in R-n onto an approximation Omega of (Omega) over cap utilizing a linear combination of higher-order spline functions. The numerical simulation is then carried out in the computational domain (Omega) over cap via a 'pull back' using the mapping operator x. The advantage is that the flexibility of higher-order spline-functions usually allows for an accurate description of (Omega) over cap with much fewer elements which can significantly reduce the computational effort required for this step compared to traditional low-order methods. Furthermore, an analytical description of the geometry can be turned back into a traditional (structured or unstructured) grid by performing a large number of function evaluations in x. This can, for instance, be utilized for local refinement without the need for remeshing. A potential drawback of curved instead of linear elements is that the meshing techniques required for the creation of folding-free mappings tend to be more sophisticated and that it is a less trivial task to verify that the resulting mapping is indeed bijective. For the purpose of creating folding-free mappings utilizing spline-functions, we will present an algorithm adopting the principles of Elliptic Grid Generation (EGG) whose basic principles have been adapted to meet the needs of IgA. In R-2, EGG has particularly appealing properties since bijectivity of the resulting mapping is guaranteed as long as the numerical accuracy to the computational approach of the mapping is sufficient. We will present an algorithm that is capable of generating folding-free mappings from a large number of geometry contours, including complicated geometries from industrial applications with extreme aspect ratios. This is accomplished by combining EGG with automatized reparameterization techniques and a sophisticated numerical approach for solving resulting governing (nonlinear) equations. The algorithm is equipped with the means to verify the bijectivity of the resulting mapping and with automatized defect-correction methods in case a violation of bijectivity is detected. Furthermore, we will present possible strategies for the generation of folding-free mappings for certain types of volumetric geometries by combining EGG with transfinite-interpolation as well as a number of other applications such as time-dependent settings. All applications are provided with example geometries. (C) 2018 Elsevier B.V. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available