5.1 New software

6.1 Mesh glossary

Participants: Frédéric Alauzet [correspondant], Paul Louis George [correspondant], Adrien Loseille.

After the publication of the third volume on Meshing, Geometric Modeling and Numerical Simulation 3 was published in 2020 14, 31, 29, books also available in French 13, 30, 28, a glossary in French on mesh has been written 6. From A to Z, this glossary provides definitions and gives a number of comments. Moreover, subtilities are given probably not so easy to understand even by people familiar with the french dialect.

6.2 Numerical simulations on GPU with the GMlib v3.0 library

Participants: Loïc Maréchal [correspondant], Julien Vanharen.

The whole library was completely rewritten to implement an automatic finite-element shader generation that converts a simple user source code into an OpenCL source that is in compiled on the GPU at run time. The library handles all meshing data structures, from file reading, renumbering and vectorizing for efficient access on the GPU, and transfer to the graphic card, all automatically and transparently. With this framework, the user can focus on the calculation part of the code, known as kernel, as all the rest is taken care of by the library. The OpenCL language was chosen as it is hardware agnostic and runs on any GPU (Intel, Nvidia and AMD) and can also use the multicore and vector capacities of modern CPUs.

Julien Vanharen developed a basic heat solver using the v3.0 as a test case so we could validate the software with various boundary conditions, calculation scheme, unstructured meshes and different memory access patterns with success. Even with basic calculation which does not stress the full GPU's power, we achieved two orders of magnitude greater speed against a single CPU core and one order of magnitude compared to a multithreaded implementation. As Julien moved to ONERA, we plan on setting up a collaboration between the two teams to implement more complex HPC codes.

6.3 High Order Meshing: from a straight mesh to a curved one

Participants: Loïc Maréchal [correspondant].

Works continued on P1toPk, a software that transform any first order hybrid mesh (triangles, quads, tets, pyramids, prisms and hexes) into a second order one while respecting a prescribe surface curvature. Efforts were made on boundary layers curving, which was challenging because jacobian validity is harder to guarantee as the elements get highly stretched, and a lot effort were also made to speed up the code by optimizing mathematical operations and parallelizing them. The code is now mature enough to be sent to industrial users for real life usage and we are waiting for valuable feedback in the present year.

6.4 Pixel-exact rendering for high-order meshes and solutions

Participants: Adrien Loseille [correspondant], Rémi Feuillet, Matthieu Maunoury.

Classic visualization software like ParaView 40, TecPlot 41, FieldView 46, Ensight 39, Medit 26, Vizir (OpenGL legacy based version) 52, Gmsh 36, ...  historically rely on the display of linear triangles with linear solutions on it. More precisely, each element of the mesh is divided into a set of elementary triangles. At each vertex of the elementary triangle is attached a value and an associated color. The value and the color inside the triangle is then deduced by a linear interpolation inside the triangle. With the increase of high-order methods and high-order meshes, these softwares adapted their technology by using subdivision methods. If a mesh has high-order elements, these elements are subdivided into a set of linear triangles in order to approximate the shape of the high-order element 82. Likewise, if a mesh has a high-order solution on it, each element is subdivided into smaller linear triangles in order to approximate the rendering of the high-order solution on it. The subdivision process can be really expensive if it is done in a naive way. For this reason, mesh adaptation procedures 69, 57, 58 are used to efficiently render high-order solutions and high-order elements using the standard linear rendering approaches. Even when optimized these approaches do have a huge RAM memory footprint as the subdivision is done on CPU in a preprocessing step. Also the adaptive subdivision process can be dependent on the palette (i.e. the range of values where the solution is studied) as the color only vary when the associated value is in this range. In this case, a change of palette inevitably imposes a new adaptation process. Finally, the use of a non conforming mesh adaptation can lead to a discontinuous rendering for a continuous solution.

Other approaches are specifically devoted to high-order solutions and are based on ray casting 63, 64, 66. The idea is for a given pixel, to find exactly its color. To do so, for each pixel, rays are cast from the position of the screen in the physical space and their intersection with the scene determines the color for the pixel. If high-order features are taken into account, it determines the color exactly for this pixel. However, this method is based on two non-linear problems: the root-finding problem and the inversion of the geometrical mapping. These problems are really costly and do not compete with the interactivity of the standard linear rendering methods even when these are called with a subdivision process unless they are done conjointly on the GPU. However, synchronization between GPU and OpenGL buffer are non-trivial combination.

The proposed method intends to be a good compromise between both methods. It does guarantee pixel-exact rendering on linear elements without extra subdivision or ray casting and it keeps the interactivity of a classical method. Moreover, the subdivision of the curved entities is done on the fly on GPU which leaves the RAM memory footprint at the size of the loaded mesh.

We are developing a software, ViZiR 4, with exact non linear solution rendering to address the high-order visualization challenge 1. ViZiR 4 is bundled as a light, simple and interactive high-order meshes and solutions visualization software. It is based on OpenGL 4 core graphic pipeline. The use of OpenGL Shading Language (GLSL) allows to perform pixel exact rendering of high order solutions on straight elements and almost pixel exact rendering on curved elements (high-order meshes). ViZiR 4 enables the representation of high order meshes (up to degree 4) and high order solutions (up to degree 10) with pixel exact rendering. Furthermore, in comparison with standard rendering techniques based on legacy OpenGL, the use of OpenGL 4 core version improves the speed of rendering, reduces the memory footprint and increases the flexibility. Many post-processing tools, such as picking, hidding surfaces, isolines, clipping, capping, are integrated to enable on the fly the analysis of the numerical results.

We added in ViZiR 4 the visualization of polygonal and polyhedral meshes. Such meshes offer flexibility as the number of vertices and faces are arbitrary. Dual meshes are examples of such meshes. Many visualization softwares do not handle polygons and when it is the case, the interactivity is often limited. We showed that our graphic pipeline can be used to render polygons. New keywords have been introduced in the libMeshb library to define these polygons and polyhedra. Furthermore, some functions have been introduced to ease the access of these data. New tessellation algorithms have been developed and do not add extra vertices in the tessellation as the aim is to minimize the number of triangles to maximize the rendering performances.

6.5 High-order mesh generation

Participants: Frédéric Alauzet [correspondant], Adrien Loseille, Rémi Feuillet, Dave Marcum, Lucien Rochery.

For years, the resolution of numerical methods has consisted in solving Partial Derivative Equations by means of a piecewise linear representation of the physical phenomenon on linear meshes. This choice was merely driven by computational limitations. With the increase of the computational capabilities, it became possible to increase the polynomial order of the solution while keeping the mesh linear. This was motivated by the fact that even if the increase of the polynomial order requires more computational resources per iteration of the solver, it yields a faster convergence of the approximation error 3  81 and it enables to keep track of unsteady features for a longer time and with a coarser mesh than with a linear approximation of the solution. However, in  17, 45, it was theoretically shown that for elliptic problems the optimal convergence rate for a high-order method was obtained with a curved boundary of the same order and in  12, evidence was given that without a high-order representation of the boundary the studied physical phenomenon was not exactly solved using a high-order method. In  85, it was even highlighted that, in some cases, the order of the mesh should be of a higher degree than the one of the solver. In other words, if the used mesh is not a high-order mesh, then the obtained high-order solution will never reliably represent the physical phenomenon.

Based on these issues, the development of high-order mesh generation procedures appears mandatory 22. To generate high-order meshes, several approaches exist. The first approach was tackled twenty years ago  20 for both surface and volume meshing. At this moment the idea was to use all the meshing tools to get a valid high-order mesh. The same problem was revisited a few years later in  75 for bio-medical applications. In these first approaches and in all the following, the underlying idea is to use a linear mesh and elevate it to the desired order. Some make use of a PDE or variational approach to do so  7, 67, 23, 60, 80, 83, 37, others are based on optimization and smoothing operations and start from a linear mesh with a constrained high-order curved boundary in order to generate a suitable high-order mesh 44, 27, 78. Also, when dealing with Navier-Stokes equations, the question of generating curved boundary layer meshes (also called viscous meshes) appears. Most of the time, dedicated approaches are set-up to deal with this problem  61, 43. In all these techniques, the key feature is to find the best deformation to be applied to the linear mesh and to optimize it. The prerequisite of these methods is that the initial boundary is curved and will be used as an input data. A natural question is consequently to study an optimal position of the high-order nodes on the curved boundary starting from an initial linear or high-order boundary mesh. This can be done in a coupled way with the volume  70, 79 or in a preprocessing phase  71, 72. In this process, the position of the nodes is set by projection onto the CAD geometry or by minimization of an error between the surface mesh and the CAD surface. Note that the vertices of the boundary mesh can move as well during the process. In the case of an initial linear boundary mesh with absence of a CAD geometry, some approaches based on normal reconstructions can be used to create a surrogate for the CAD model  82, 38. Finally, a last question remains when dealing with such high-order meshes: Given a set of degrees of freedom, is the definition of these objects always valid ? Until the work presented in  33, 42, 32, no real approach was proposed to deal in a robust way with the validity of high-order elements. The novelty of these approaches was to see the geometrical elements and their Jacobian as Bézier entities. Based on the properties of the Bézier representation, the validity of the element is concluded in a robust sense, while the other methods were only using a sampling of the Jacobian to conclude about its sign without any warranty on the whole validity of the elements.

In this context, several issues have been addressed : the analogy between high-order and Bézier elements, the development of high-order error estimates suitable for parametric high-order surface mesh generation and the generalization of mesh optimization operators and their applications to curved mesh generation, moving-mesh methods, boundary layer mesh generation and mesh adaptation.

Metric fields are the link between particular error estimates - be they for low-order 47, 48 or high-order methods 19, for the solution of a PDE 50 or a quantity of interest derived from it such as drag or lift 51 - and automatic mesh adaptation. In the case of linear meshes, a metric field locally distorts the measure or distance such that, when the mesh adaptation algorithm has constructed an uniform mesh in the induced Riemannian space, it is strongly anisotropic in the usual Euclidean (physicial) space. As such, anisotropy arises naturally, without it ever being explicitely sought by the (re)meshing algorithm.

We seek to extend these principles of metric-based $P^1$ adaptation to high-order meshes. In particular, we expect the meshing process to naturally recover curvature from the variations of the metric field, very much like $P^1$ remeshing recovers anisotropy from local values of the metric field. As such, curvature must be the consequence of a simple geometric property computed in the Riemannian space, like anisotropy is the consequence of unitness in a space where distances are distorted. Therefore, we propose Riemannian edge length minimization (or geodesic seeking as in 84) as the driver for metric field curvature recovery.

The metric field's own intrinsic curvature may derive from any error estimate, be it boundary approximation error 24, 21 or an interpolation error estimate. So far, interpolation error estimates on high-order elements are limited to isotropy (16 in $L^2$ and 62 in $L^1$ norms) or require that the curvature of the element be bounded, essentially establishing a range where it may be considered linear 10.

Robustness and modularity of the general remeshing algorithm may be derived from the use of a single topological operator such as the cavity operator 49, 53, 54. This operator remeshes element subsets (so-called cavities) by reconnecting cavity boundary nodes to a given point in space, already belonging to the mesh or otherwise. This very elementary operation can handle the most common topological operations: insertions, collapses, edge or face swaps. Therefore, it is central both to mesh adaptation (node creation, deletion) and to mesh optimization (mainly swaps).

A new $P^2$ cavity operator was devised and implemented in $3D$. Volume curvature is driven, first, by a purely metric-based curving procedure using a Riemannian edge length minimization algorithm. The considered cost-function is strongly non-linear, as it depends on the geometry of the edge as well as on the metric field along the edge which, itself, depends on the geometry of the edge. Analytical derivatives of this quantity with regards to edge shape were computed nonetheless, which was necessary to implement optimization in $3D$ with reasonable CPU costs. This first required the extension of the log-euclidean metric interpolation scheme to high-order elements. A version for simplices of all dimensions and degrees was proposed based on Lagrange interpolation of the log-metrics. This means that the metric field along the $P^2$ edge can now be described in terms of the metrics at the two extremities and at the Lagrange control point. The dependency of the metric at the Lagrange control point must in turn be identified. It is based on localization of the point in a so-called back-mesh and, again, log-euclidean metric interpolation in the host back element. At last, derivatives of Riemannian edge length could be computed, and were used to implement a Newton minimization algorithm that is capable of finding local minima of edge length at a rate of 4e4 edges per second while remaining entirely consistent with the log-euclidean framework used in our anisotropic remeshing algorithm. Secondly, final cavities are made valid or improved using a new Jacobian correction algorithm that can very quickly optimize entire shells of edges using a recently identified property of Jacobian determinant control coefficients, namely that they depend linearly on a given control point. Using the very powerful simplex algorithm for linear programs, we can find a global solution to maximizing the minimum Jacobian of all elements sharing an edge in a single pass (with average runtime of 5e-5 seconds on a laptop). Three communications have been made on the subject in the last year, with peer-reviewed proceedings 3, 4, 2. Applications on large meshes deriving from industrial cases were carried out, illustrating the ability of these optimization procedures to curved meshes of 20M elements in under 10min.

6.6 Unstructured anisotropic mesh adaptation for 3D RANS turbomachinery applications

Participants: Frédéric Alauzet [correspondant], Adrien Loseille [correspondant], Julien Vanharen.

The scope of this paper is to demonstrate the viability and efficiency of unstructured anisotropic mesh adaptation techniques to turbomachinery applications. The main difficulty in turbomachinery is the periodicity of the domain that must be taken into account inthe solution mesh-adaptive process. The periodicity is strongly enforced in the flow solver using ghost cells to minimize the impact on the source code. For the mesh adaptation, the local remeshing is done in two steps. First, the inner domain is remeshed with frozen periodic frontiers, and, second, the periodic surfaces are remeshed after moving geometrice ntities from one side of the domain to the other. One of the main goal of this work is to demonstrate how mesh adaptation, thanks to its automation, is able to generate meshes that are extremely difficult to envision and almost impossible to generate manually. This study only considers feature-based error estimate based on the standard multi-scale Lpinterpolation error estimate. We presents all the specific modifications that have been introduced in the adaptive process to deal with periodic simulations used for turbomachinery applications. The periodic mesh adaptation strategy is then tested and validated on the LS89 high pressure axial turbine vane and the NASA Rotor 37 test cases.

6.7 Mixed-element mesh adaptation for CFD simulations

Participants: Frédéric Alauzet [correspondant], Lucille Marie Tenkès, Julien Vanharen, Adrien Loseille, Cosimo Tarsia Morisco.

Due to their various nature, physical phenomena that we seek to capture in CFD simulations may have specific specific mesh requirements. For example, to solve the boundary layer, some numerical schemes favor structured meshes respecting alignment with the boundary of the domain, while these constraints are not necessary elsewhere. Our approach is to use the techniques of metric-based mesh adaptation to generate a mixed-element mesh that can fulfill these different mesh requirements. This approach is based on the metric-orthogonal point-placement, creating some structured parts from the intrinsic directional information bore by the metric-field. Some unstructured areas may remain where structure is not needed. The main goals of this work are to improve the orthogonality of the output mesh and its alignment with the metric field. This work has three main axes. First, we have improved the preprocessing gradation step to smooth the metric field and improve the orthogonality of the final mesh. Then, we have studied two methods to obtain quadrilaterals: one using an a priori quadrilaterals recombination, the other detecting straightforwardly the orthogonal patterns during the re-meshing step. Finally, the work on the solver Wolf has been carried on and corrected to perform robust and accurate simulations on mixed-element meshes. These three developments were embodied in a mixed-element adaptation loop. The first two topics are detailed in what follows.

6.8 Coupled flow and adjoint solver

Participants: Francesco Clerici [correspondant], Frédéric Alauzet.

When solving the RANS equations, usually one decouples the equations relative to the mean-flow and the equations relative to turbulence. This division provides two separated systems to be solved at each time step, one relative to the mean-flow and the other relative to the turbulence. This presents two main drawbacks: in the flow solver, the Jacobian of the system lacks of the terms bounding the mean-flow and the turbulence, and this can slow down the residual convergence. The second drawback regards the adjoint problem, which consists into a linear system assembled with the transpose of the Jacobian matrix of the residuals and, on the right-hand side, the derivative of an aeronautical coefficient with respect to the flow variables. A Jacobian missing the coupling terms between the mean-flow and the turbulence provides a null adjoint turbulent viscosity, and this is a limitation in the development of more complex discretization error estimates. We have therefore developed a 2D version of the coupled flow and adjoint solver which includes in the Jacobian the coupling terms between the mean-flow and the turbulence. When have tested this method on the 2D geometry provided for the 4th CFD AIAA High Lift Prediction Workshop, and the result provided several features to emerge, such as a high mesh refinement inside the boundary layer of the leading edges, and inside the regions of high turbulence destruction. The work is pursued in collaboration with Philippe Spalart (Boeing), and is presented in 18.

6.9 Turbulent error estimate

Participants: Francesco Clerici [correspondant], Frédéric Alauzet.

Goal-oriented mesh adaptation is a methodology used to adapt the mesh in order to minimize the discretization error commited on a functional depending on the solution. As an intermediate step, one finds an upper bound to such discretization error taking the form of a weighted sum of the interpolation errors of the solution, and such upper bound is called error estimate. Regarding Wolf and the RANS equations, up to now we have focused only on the mean-flow part of such an error estimate, meaning that the terms coming from turbulence have been neglected. The scope of this work is to enrich the error estimate with the information coming from the turbulence. In particular, the methodology has been tested on the 2D geometry provided for the 4th CFD AIAA High Lift Prediction Workshop, providing high mesh refinements on the boundaries of the turbulent regions. Also this work is pursued in collaboration with Philippe Spalart (Boeing), and is presented in 18.

7.1 Bilateral contracts with industry

Participants: Frédéric Alauzet [correspondant], Adrien Loseille.

• Boeing
• Safran Tech

7.2 Bilateral grants with industry

Participants: Frédéric Alauzet [correspondant], Adrien Loseille.

• Projet RAPID DGA

8.1 Publications of the year

8.2 Cited publications

• 7 articleR.R. Abgrall, C.C. Dobrzynski and A.A. Froehly. A method for computing curved meshes via the linear elasticity analogy, application to fluid dynamics problems.International Journal for Numerical Methods in Fluids7642014, 246--266
  back to text
• 8 articleF.F. Alauzet and A.A. Loseille. A decade of progress on anisotropic mesh adaptation for Computational Fluid Dynamics.722016, 13-39
  back to text
• 9 articleF.F. Alauzet and A.A. Loseille. High Order Sonic Boom Modeling by Adaptive Methods.2292010, 561-593
  back to text
• 10 bookT.Thomas Apel. Anisotropic finite elements: local estimates and applications.3Teubner Stuttgart1999
  back to text
• 11 articleN.N. Barral, G.G. Olivier and F.F. Alauzet. Metric-based anisotropic mesh adaptation for three-dimensional time-dependent problems involving moving geometries.3312017, 157-187
  back to text
• 12 articleF.F. Bassi and S.S. Rebay. High-order accurate discontinuous finite element solution of the 2D Euler equations.Journal of computational physics13821997, 251--285
  back to text
• 13 bookH.Houman Borouchaki and P. L.Paul Louis George. Maillage, modélisation géométrique et simulation numérique 1: Fonctions de forme, triangulations et modélisation géométrique.1ISTE Group2017
  back to text
• 14 bookH.Houman Borouchaki and P. L.Paul Louis George. Meshing, Geometric Modeling and Numerical Simulation 1: Form Functions, Triangulations and Geometric Modeling.John Wiley & Sons2017
  back to text
• 15 articleH.H. Borouchaki, P.P. Laug and P.P.L. George. Parametric surface meshing using a combined advancing-front -- generalized-Delaunay approach.International Journal for Numerical Methods in Engineering491-22000, 233-259
  back to text
• 16 articleL.Lorenzo Botti. Influence of Reference-to-Physical Frame Mappings on Approximation Properties of Discontinuous Piecewise Polynomial Spaces.Journal of Scientific Computing5209 2012
  DOIback to text
• 17 incollectionP. G.P. G. Ciarlet and P.-A.P.-A. Raviart. The combined effect of curved boundaries and numerical integration in isoparametric finite element methods.The mathematical foundations of the finite element method with applications to partial differential equationsElsevier1972, 409--474
  back to text
• 18 articleF.F. Clerici, P.P. Spalart and F.F. Alauzet. Coupled adjoint solver and turbulent error estimate for anisotropic mesh adaptation in high-fidelity RANS simulations..AIAA J.2022
  back to textback to text
• 19 articleO.Olivier Coulaud and A.Adrien Loseille. Very High Order Anisotropic Metric-Based Mesh Adaptation in 3D.Procedia Engineering16325th International Meshing Roundtable2016, 353 - 365URL: http://www.sciencedirect.com/science/article/pii/S187770581633380X
  DOIback to text
• 20 inproceedingsS.S. Dey, R. M.R. M. O'bara and M. S.M. S. Shephard. Curvilinear Mesh Generation in 3D.Proceedings of the 7th International Meshing Roundtable1999, 407--417
  back to text
• 21 inproceedingsR.Rémi Feuillet, O.Olivier Coulaud and A.Adrien Loseille. Anisotropic Error Estimate for High-order Parametric Surface Mesh Generation.28th International Meshing RoundtableBuffalo, NY, United StatesOctober 2019
  HALDOIback to text
• 22 articleR.Rémi Feuillet, A.Adrien Loseille and F.Frédéric Alauzet. Optimization of P2 meshes and applications.Computer-Aided Design124April 2020, 102846
  HALDOIback to text
• 23 articleM.M. Fortunato and P.-O.P.-O. Persson. High-order Unstructured Curved Mesh Generation Using the Winslow Equations.J. Comput. Phys.307February 2016, 1--14
  back to text
• 24 miscP. J.Pascal J. Frey. About Surface Remeshing.2000
  back to text
• 25 inproceedingsP.P.J. Frey. About surface remeshing.Proceedings of the 9th international meshing roundtableNew Orleans, LO, USA2000, 123-136
  back to text
• 26 miscP. J.P. J. Frey. Medit: An interactive mesh visualization software, INRIA Technical Report RT0253.2001
  back to text
• 27 incollectionA.A. Gargallo-Peiró, X.X. Roca, J.J. Peraire and J.J. Sarrate. Defining quality measures for mesh optimization on parameterized CAD surfaces.Proceedings of the 21st International Meshing RoundtableSpringer2013, 85--102
  back to text
• 28 bookP. L.Paul Louis George, F.Frédéric Alauzet, A.Adrien Loseille and L.Loïc Maréchal. Maillage, modélisation géométrique et simulation numérique 3: Stockage, transformation, utilisation et visualisation de maillage.4ISTE Group2020
  back to text
• 29 bookP. L.Paul Louis George, F.Frédéric Alauzet, A.Adrien Loseille and L.Loïc Maréchal. Meshing, Geometric Modeling and Numerical Simulation 3: Storage, Visualization and In Memory Strategies.John Wiley & Sons2020
  back to text
• 30 bookP. L.Paul Louis George, H.Houman Borouchaki, F.Frédéric Alauzet, P.Patrick Laug, A.Adrien Loseille and L.Loïc Maréchal. Maillage, modélisation géométrique et simulation numérique 2: Métriques, maillages et adaptation de maillages.2ISTE Group2018
  back to text
• 31 bookP. L.Paul Louis George, H.Houman Borouchaki, F.Frederic Alauzet, P.Patrick Laug, A.Adrien Loseille and L.Loïc Maréchal. Meshing, Geometric Modeling and Numerical Simulation, Volume 2: Metrics, Meshes and Mesh Adaptation.John Wiley & Sons2019
  back to text
• 32 articleP.P.~L. George, H.H. Borouchaki and N.N. Barral. Geometric validity (positive jacobian) of high-order Lagrange finite elements, theory and practical guidance.Engineering with computers3232016, 405--424
  back to text
• 33 articleP.P.~L. George and H.H. Borouchaki. Construction of tetrahedral meshes of degree two.International Journal for Numerical Methods in Engineering9092012, 1156,1182
  back to textback to text
• 34 articleP.P.L. George and H.H. Borouchaki. ``Ultimate'' robustness in meshing an arbitrary polyhedron.5872003, 1061-1089
  back to text
• 35 articleP.P.L. George, F.F. Hecht and E.E. Saltel. Automatic mesh generator with specified boundary.921991, 269-288
  back to text
• 36 articleC.C. Geuzaine and J.-F.J.-F. Remacle. Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities.International Journal for Numerical Methods in Engineering79112009, 1309-1331
  back to text
• 37 articleR.R. Hartmann and T.T. Leicht. Generation of unstructured curvilinear grids and high-order discontinuous Galerkin discretization applied to a 3D high-lift configuration.International Journal for Numerical Methods in Fluids8262016, 316-333
  back to text
• 38 articleJ.J. Ims and Z. J.Z. J. Wang. Automated low-order to high-order mesh conversion.Engineering with Computers351Jan 2019, 323--335
  back to text
• 39 miscA.Ansys Inc.. Ensight.
  back to text
• 40 miscK.KitWare Inc.. ParaView.
  back to text
• 41 miscT.TecPlot Inc.. TecPlot.
  back to text
• 42 articleA.A. Johnen, J.-F.J.-F. Remacle and C.C. Geuzaine. Geometrical validity of curvilinear finite elements.Journal of Computational Physics233152013, 359-372
  back to text
• 43 inproceedingsS. L.S. L. Karman. Curving for Viscous Meshes.27th International Meshing RoundtableChamSpringer International Publishing2019, 303--325
  back to text
• 44 incollectionS. L.S. L. Karman, J. T.J T. Erwin, R. S.R. S. Glasby and D.D. Stefanski. High-Order Mesh Curving Using WCN Mesh Optimization.46th AIAA Fluid Dynamics ConferenceAIAA AVIATION ForumAmerican Institute of Aeronautics and Astronautics2016
  back to text
• 45 articleM.M. Lenoir. Optimal isoparametric finite elements and error estimates for domains involving curved boundaries.SIAM journal on numerical analysis2331986, 562--580
  back to text
• 46 miscI.Intelligent Light. FieldView.
  back to text
• 47 articleA.Adrien Loseille and F.Frédéric Alauzet. Continuous mesh framework part I: well-posed continuous interpolation error.SIAM Journal on Numerical Analysis4912011, 38--60
  back to text
• 48 articleA.Adrien Loseille and F.Frédéric Alauzet. Continuous mesh framework part II: validations and applications.SIAM Journal on Numerical Analysis4912011, 61--86
  back to text
• 49 articleA.Adrien Loseille, F.Frédéric Alauzet and V.Victorien Menier. Unique cavity-based operator and hierarchical domain partitioning for fast parallel generation of anisotropic meshes.Computer-Aided Design852017, 53-67
  back to textback to text
• 50 phdthesisA.Adrien Loseille. Anisotropic 3D hessian-based multi-scale and adjoint-based mesh adaptation for Computational fluid dynamics: Application to high fidelity sonic boom prediction..Université Pierre et Marie Curie - Paris VIDecember 2008
  HALback to text
• 51 articleA.A. Loseille, A.A. Dervieux and F.F. Alauzet. Fully anisotropic goal-oriented mesh adaptation for 3D steady Euler equations.Journal of Computational Physics22982010, 2866 - 2897URL: http://www.sciencedirect.com/science/article/pii/S0021999109007001
  DOIback to text
• 52 miscA.A. Loseille, H.H. Guillard and A.A. Loyer. An introduction to Vizir: an interactive mesh visualization and modification software.2016
  back to text
• 53 inbookA.Adrien Loseille and R.Rainald Lohner. Cavity-Based Operators for Mesh Adaptation.51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace ExpositionURL: https://arc.aiaa.org/doi/abs/10.2514/6.2013-152
  DOIback to text
• 54 inbookA.Adrien Loseille. Recent Improvements on Cavity-Based Operators for RANS Mesh Adaptation.2018 AIAA Aerospace Sciences MeetingURL: https://arc.aiaa.org/doi/abs/10.2514/6.2018-0922
  DOIback to text
• 55 inproceedingsL.L. Maréchal. A new approach to octree-based hexahedral meshing.2001, 209--221
  back to text
• 56 inproceedingsL.L. Maréchal. Advances in Octree-Based All-Hexahedral Mesh Generation: Handling Sharp Features.18Salt Lake City, UT, USA2009, 65-84
  back to text
• 57 articleM.M. Maunoury, C.C. Besse, V.V. Mouysset, S.S. Pernet and P.-A.P.-A. Haas. Well-suited and adaptive post-processing for the visualization of hp simulation results.Journal of Computational Physics3752018, 1179 - 1204
  back to text
• 58 phdthesisM.M. Maunoury. Méthode de visualisation adaptée aux simulations d'ordre élevé. Application à la compression-reconstruction de champs rayonnés pour des ondes harmoniques..Université de ToulouseFebruary 2019
  back to text
• 59 inproceedingsD.D.J. Mavriplis. Results from the 3rd Drag Prediction Workshop using NSU3D unstructured mesh solver.45AIAA-2007-0256, Reno, NV, USAJan 2007
  back to text
• 60 articleD.D. Moxey, D.D. Ekelschot, Ü.Ü. Keskin, S.S.J. Sherwin and J.J. Peirò. High-order curvilinear meshing using a thermo-elastic analogy.Computer-Aided Design722016, 130 - 139
  back to text
• 61 articleD.D. Moxey, M.M.D. Green, S.S.J. Sherwin and J.J. Peiró. An isoparametric approach to high-order curvilinear boundary-layer meshing.Computer Methods in Applied Mechanics and Engineering2832015, 636 - 650
  back to text
• 62 articleD.David Moxey, S. P.Shankar P Sastry and R. M.Robert M Kirby. Interpolation error bounds for curvilinear finite elements and their implications on adaptive mesh refinement.Journal of Scientific Computing7822019, 1045--1062
  back to text
• 63 articleB.B. Nelson, R.R. Haimes and R.R.~M. Kirby. GPU-Based Interactive Cut-Surface Extraction From High-Order Finite Element Fields.IEEE Transactions on Visualization and Computer Graphics17122011, 1803--11
  back to text
• 64 articleB.B. Nelson, E.E. Liu, R. M.R. M. Kirby and R.R. Haimes. ElVis: A System for the Accurate and Interactive Visualization of High-Order Finite Element Solutions.IEEE Transactions on Visualization and Computer Graphics18122012, 2325-2334
  back to text
• 65 inproceedingsM.M. Park, J.J. Krakos, T.T.R. Michal, A.A. Loseille and J.J. Alonso. Unstructured Grid Adaptation: Status, Potential Impacts, and Recommended Investments Toward CFD Vision 2030.462016-3323, Washington, D.C., USA2016
  back to text
• 66 incollectionJ.J. Peiró, D.D. Moxey, B.B. Jordi, S. J.S. J. Sherwin, B.B.W. Nelson, R. M.R. M. Kirby and R.R. Haimes. High-Order Visualization with ElVis.Notes on Numerical Fluid Mechanics and Multidisciplinary DesignSpringer International Publishing2015, 521--534
  back to text
• 67 inproceedingsP.-O.P.-O. Persson and J.J. Peraire. Curved mesh generation and mesh refinement using Lagrangian solid mechanics.47th AIAA Aerospace Sciences Meeting including The New Horizons Forum and Aerospace Exposition2009, 949
  back to text
• 68 articleG.Guillaume Puigt, V.Valerie Auffray and J.-D.J-D Müller. Discretisation of diffusive fluxes on hybrid grids.Journal of Computational Physics22952010, 1425--1447
  back to text
• 69 articleJ.-F.J.-F. Remacle, N.N. Chevaugeon, E.E. Marchandise and C.C. Geuzaine. Efficient visualization of high-order finite elements.International Journal for Numerical Methods in Engineering6952007, 750-771
  back to text
• 70 articleE.E. Ruiz-Gironès, X.X. Roca and J.J. Sarrate. High-order mesh curving by distortion minimization with boundary nodes free to slide on a 3D CAD representation.Computer-Aided Design7223rd International Meshing Roundtable Special Issue: Advances in Mesh Generation2016, 52 - 64
  back to text
• 71 articleE.E. Ruiz-Gironès, J.J. Sarrate and X.X. Roca. Defining an $L^2$-disparity Measure to Check and Improve the Geometric Accuracy of Non-interpolating Curved High-order Meshes.Procedia Engineering1242015, 122--134
  back to text
• 72 articleE.E. Ruiz-Gironès, J.J. Sarrate and X.X. Roca. Generation of curved high-order meshes with optimal quality and geometric accuracy.Procedia engineering1632016, 315--327
  back to text
• 73 articleC.C.L. Rumsey, J.J.P. Slotnick, M.M. Long, R.R.A. Stuever and T.T.R. Wayman. Summary of the first AIAA CFD High-Lift Prediction Workshop.Journal of Aircraft4862011, 2068-2079
  back to text
• 74 articleC.C.L. Rumsey and J.J.P. Slotnick. Overview and Summary of the Second AIAA High Lift Prediction Workshop.Journal of Aircraft5242015, 1006-1025
  back to text
• 75 articleS. J.S. J. Sherwin and J.J. Peiró. Mesh generation in curvilinear domains using high-order elements.International Journal for Numerical Methods in Engineering5312002, 207-223
  back to text
• 76 techreportJ.J. Slotnick, A.A. Khodadoust, J.J. Alonso, D.D. Darmofal, W.W. Gropp, E.E. Lurie and D.D. Mavriplis. CFD Vision 2030 Study: A path to revolutionary computational aerosciences.NASAMarch 2014
  back to text
• 77 incollectionL.-M.Lucille-Marie Tenkes and F.Frédéric Alauzet. Size Gradation Control for Anisotropic Hybrid Meshes.Numerical Geometry, Grid Generation and Scientific Computing143Lecture Notes in Computational Science and EngineeringSpringer International PublishingMay 2021, 127-139
  HALDOIback to text
• 78 articleT.T. Toulorge, C.C. Geuzaine, J.-F.J.-F. Remacle and J.J. Lambrechts. Robust untangling of curvilinear meshes.Journal of Computational Physics2542013, 8 - 26
  back to text
• 79 articleT.T. Toulorge, J.J. Lambrechts and J.-F.J.-F. Remacle. Optimizing the geometrical accuracy of curvilinear meshes.Journal of Computational Physics3102016, 361 - 380
  back to text
• 80 articleM.M. Turner, J.J. Peirò and D.D. Moxey. A Variational Framework for High-order Mesh Generation.Procedia Engineering163Supplement C25th International Meshing Roundtable2016, 340 - 352
  back to text
• 81 phdthesisJ.J. Vanharen. High-order numerical methods for unsteady flows around complex geometries.Université de Toulouse2017
  back to text
• 82 inproceedingsA.A. Vlachos, P.P. Jörg, C.C. Boyd and J.J.~L. Mitchell. Curved PN Triangles.Proceedings of the 2001 Symposium on Interactive 3D Graphics2001, 159-166
  back to textback to text
• 83 articleZ. Q.Z. Q. Xie, R.R. Sevilla, O.O. Hassan and K.K. Morgan. The generation of arbitrary order curved meshes for 3D finite element analysis.Computational Mechanics513Mar 2013, 361--374
  back to text
• 84 inproceedingsR.Ruili Zhang, A.Amaury Johnen and J.-F.Jean-François Remacle. Curvilinear mesh adaptation.International Meshing RoundtableSpringer2018, 57--69
  back to text
• 85 inproceedingsP.P. Zwanenburg and S.S. Nadarajah. On the Necessity of Superparametric Geometry Representation for Discontinuous Galerkin Methods on Domains with Curved Boundaries.23rd AIAA Computational Fluid Dynamics Conference2017
  back to text