Anisotropic mesh adaptation: towards user‐independent,mesh‐independent and solver‐independent CFD. Part I: general principles

The present paper is the lead article in a three‐part series on anisotropic mesh adaptation and its applications to structured and unstructured meshes. A flexible approach is proposed and tested on two‐dimensional, inviscid and viscous, finite volume and finite element flow solvers, over a wide range of speeds. The directional properties of an interpolation‐based error estimate, extracted from the Hessian of the solution, are used to control the size and orientation of mesh edges. The approach is encapsulated into an edge‐based anisotropic mesh optimization methodology (MOM), which uses a judicious sequence of four local operations: refinement, coarsening, edge swapping and point movement, to equi‐distribute the error estimate along all edges, without any recourse to remeshing. The mesh adaptation convergence of the MOM loop is carefully studied for a wide variety of test cases. The mesh optimization generic coupling of MOM with finite volume and finite element flow solvers is shown to yield the same final mesh no matter what the starting point is. It is also shown that on such optimized meshes, the need for computational fluid dynamics (CFD) stabilization artifices, such as upwinding or artificial viscosity, are drastically reduced, if not altogether eliminated, in most well‐posed formulations. These two conclusions can be considered significant steps towards mesh‐independent and solver‐independent CFD. The structure of the three‐part series is thus, 1, general principles; 2, methodology and applications to structured and unstructured grids; 3, applications to three‐dimensional flows. Copyright © 2000 John Wiley & Sons, Ltd.     

Anisotropic mesh adaptation: towards user‐independent,mesh‐independent and solver‐independent CFD. Part II. Structured grids

The present paper is the second article in a three‐part series on anisotropic mesh adaptation and its application to (2‐D) structured and unstructured meshes. In the first article, the theory was presented, the methodology detailed and brief examples given of the application of the method to both types of grids. The second part details the application of the mesh adaptation method to structured grids. The adaptation operations are restricted to mesh movement in order to avoid the creation of hanging nodes. Being based on a spring analogy with no restrictive orthogonality constraint, a wide grid motion is allowed. The adaptation process is first validated on analytical test cases and its high efficiency is shown on relevant transonic and supersonic benchmarks. These latter test cases are also solved on adapted unstructured grids to provide a reference for comparison studies. The third part of the series will demonstrate the capability of the methodology on 2‐D unstructured test cases. Copyright © 2002 John Wiley & Sons, Ltd.     

Anisotropic mesh adaption for time‐dependent problems

We propose a space–time adaptive procedure for a model parabolic problem based on a theoretically sound anisotropic a posteriori error analysis. A space–time finite element scheme (continuous in space but discontinuous in time) is employed to discretize this problem, thus allowing for non‐matching meshes at different time levels. Copyright © 2008 John Wiley & Sons, Ltd.     

Adjoint‐based error estimation and mesh adaptation for hybridized discontinuous Galerkin methods

We present a robust and efficient target‐based mesh adaptation methodology, building on hybridized discontinuous Galerkin schemes for (nonlinear) convection–diffusion problems, including the compressible Euler and Navier–Stokes equations. The hybridization of finite element discretizations has the main advantage that the resulting set of algebraic equations has globally coupled degrees of freedom (DOFs) only on the skeleton of the computational mesh. Consequently, solving for these DOFs involves the solution of a potentially much smaller system. This not only reduces storage requirements but also allows for a faster solution with iterative solvers. The mesh adaptation is driven by an error estimate obtained via a discrete adjoint approach. Furthermore, the computed target functional can be corrected with this error estimate to obtain an even more accurate value. The aim of this paper is twofold: Firstly, to show the superiority of adjoint‐based mesh adaptation over uniform and residual‐based mesh refinement and secondly, to investigate the efficiency of the global error estimate. Copyright © 2014 John Wiley & Sons, Ltd.     

A parallel cell‐based DSMC method on unstructured adaptive meshes

A parallel DSMC method based on a cell‐based data structure is developed for the efficient simulation of rarefied gas flows on PC‐clusters. Parallel computation is made by decomposing the computational domain into several subdomains. Dynamic load balancing between processors is achieved based on the number of simulation particles and the number of cells allocated in each subdomain. Adjustment of cell size is also made through mesh adaptation for the improvement of solution accuracy and the efficient usage of meshes. Applications were made for a two‐dimensional supersonic leading‐edge flow, the axi‐symmetric Rothe's nozzle, and the open hollow cylinder flare flow for validation. It was found that the present method is an efficient tool for the simulation of rarefied gas flows on PC‐based parallel machines. Copyright © 2004 John Wiley & Sons, Ltd.     

Comparison of adjoint‐based and feature‐based grid adaptation for functional outputs

Adjoint‐based and feature‐based grid adaptive strategies are compared for their robustness and effectiveness in improving the accuracy of functional outputs such as lift and drag coefficients. The output‐based adjoint approach strives to improve the adjoint error estimates that relate the local residual errors to the global error in an output function via adjoint variables as weight functions. A conservative adaptive indicator that takes into account the residual errors in both the primal (flow) and dual (adjoint) solutions is implemented for the adjoint approach. The physics‐based feature approach strives to identify and resolve significant features of the flow to improve functional accuracy. Adaptive indicators that represent expansions and compressions in the flow direction and gradients normal to the flow direction are implemented for the feature approach. The adaptive approaches are compared for functional outputs of three‐dimensional arbitrary Mach number flow simulations on mixed‐element unstructured meshes. Grid adaptation is performed with h‐refinement and results are presented for inviscid, laminar and turbulent flows. Copyright © 2006 John Wiley & Sons, Ltd.     

Sequential feature‐based mesh movement and adjoint error‐based mesh refinement

Nowadays, aerodynamic computational modeling is carried out on a daily basis in an industrial setting. This is done with the aim of predicting the performance and flow characteristics of new components. However, limited resources in terms of time and hardware force the engineer to employ relatively coarse computational grids, thus achieving results with variable degree of inaccuracy. In this article, a novel combination of feature and adjoint‐based mesh adaptation methods is investigated and applied to typical three‐dimensional turbomachinery cases, such as compressor and fan blades. The proposed process starts by employing feature‐based mesh movement to improve the global flow solution and then adjoint refinement to tune the mesh for each quantity of interest. Comparison of this process with one utilizing only the adjoint refinement procedure shows significant benefits in terms of accuracy of the performance quantity.     

Drag prediction,decomposition and visualization in unstructured mesh CFD solver of TAS‐code

The accuracy of drag prediction in unstructured mesh CFD solver of TAS (Tohoku University Aerodynamic Simulation) code is discussed using a drag decomposition method. The drag decomposition method decomposes total drag into wave, profile, induced and spurious drag components, the latter resulting from numerical diffusion and errors. The mesh resolution analysis is conducted by the drag decomposition method. The effect of an advanced unstructured mesh scheme of U‐MUSCL reconstruction is also investigated by the drag decomposition method. The computational results show that the drag decomposition method reliably predicts drag and is capable of meaningful drag decomposition. The accuracy of drag prediction is increased by eliminating the spurious drag component from the total drag. It is also confirmed that the physical drag components are almost independent of the mesh resolution and scheme modification. Copyright © 2007 John Wiley & Sons, Ltd.     

A moving mesh algorithm for aero‐thermo‐mechanical modelling in turbomachinery

This paper describes the development of a mesh deformation method used for aero‐thermo‐mechanical coupling of turbo‐engine components. The method is based on the nonlinear solution of an elastic medium analogy, solved using finite element discretisation and modified to let the boundary nodes be free to slide over the deflected surfaces. This sliding technique relies on a B‐spline reconstruction of the moving boundary and increases the robustness of the method in situations where the boundary deflection field presents significant gradients or large relative motion between two distinct boundaries. The performance of the method is illustrated with the application to an interstage cavity of a turbine assembly, subjected to the deformations computed by a coupled thermo‐mechanical analysis of the engine component. Copyright © 2011 John Wiley & Sons, Ltd.     

A gas‐kinetic scheme for the simulation of turbulent flows on unstructured meshes

This paper presents a numerical method for simulating turbulent flows via coupling the Boltzmann BGK equation with Spalart–Allmaras one equation turbulence model. Both the Boltzmann BGK equation and the turbulence model equation are carried out using the finite volume method on unstructured meshes, which is different from previous works on structured grid. The application of the gas‐kinetic scheme is extended to the simulation of turbulent flows with arbitrary geometries. The adaptive mesh refinement technique is also adopted to reduce the computational cost and improve the efficiency of meshes. To organize the unstructured mesh data structure efficiently, a non‐manifold hybrid mesh data structure is extended for polygonal cells. Numerical experiments are performed on incompressible flow over a smooth flat plate and compressible turbulent flows around a NACA 0012 airfoil using unstructured hybrid meshes. These numerical results are found to be in good agreement with experimental data and/or other numerical solutions, demonstrating the applicability of the proposed method to simulate both subsonic and transonic turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.     

Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

A stabilized Nitsche‐type extended embedding mesh approach for 3D low‐ and high‐Reynolds‐number flows

This paper presents a stabilized extended finite element method (XFEM) based fluid formulation to embed arbitrary fluid patches into a fixed background fluid mesh. The new approach is highly beneficial when it comes to computational grid generation for complex domains, as it allows locally increased resolutions independent from size and structure of the background mesh. Motivating applications for such a domain decomposition technique are complex fluid‐structure interaction problems, where an additional boundary layer mesh is used to accurately capture the flow around the structure. The objective of this work is to provide an accurate and robust XFEM‐based coupling for low‐ as well as high‐Reynolds‐number flows. Our formulation is built from the following essential ingredients: Coupling conditions on the embedded interface are imposed weakly using Nitsche's method supported by extra terms to guarantee mass conservation and to control the convective mass transport across the interface for transient viscous‐dominated and convection‐dominated flows. Residual‐based fluid stabilizations in the interior of the fluid subdomains and accompanying face‐oriented fluid and ghost‐penalty stabilizations in the interface zone stabilize the formulation in the entire fluid domain. A detailed numerical study of our stabilized embedded fluid formulation, including an investigation of variants of Nitsche's method for viscous flows, shows optimal error convergence for viscous‐dominated and convection‐dominated flow problems independent of the interface position. Challenging two‐dimensional and three‐dimensional numerical examples highlight the robustness of our approach in all flow regimes: benchmark computations for laminar flow around a cylinder, a turbulent driven cavity flow at Re = 10000 and the flow interacting with a three‐dimensional flexible wall. Copyright © 2016 John Wiley & Sons, Ltd.     

Internal wave reflections and transmissions arising from a non-uniform mesh. Part I: An analysis for the Crank–Nicolson linear finite element scheme

This is the first of a series of three related papers dealing with some of the consequences of non-uniform meshes in a numerical model. In this paper the accuracy of the Crank–Nicolson linear finite element scheme, which is applied to the linear shallow water equations, is examined in the context of a single abrupt change in nodal spacing. The (in)accuracy is quantified in terms of reflection and transmission coefficients. An incident wave impinging on the interface between two regions with different nodal spacings is shown to give rise to no reflected waves and two transmitted waves. The analysis is verified using three different wavelengths (2Δx, 4Δx 8Δx) in three ‘hot-start’ numerical experiments with a mesh expansion factor of 2 and three experiments with a mesh contraction factor of 1/2. An energy flux analysis based on the concept of group velocity shows that energy is conserved across the interface.     

Residual‐based turbulence models for moving boundary flows: hierarchical application of variational multiscale method and three‐level scale separation

This paper presents residual‐based turbulence models for problems with moving boundaries and interfaces. The method is developed via a hierarchical application of variational multiscale ideas and the models are cast in an arbitrary Lagrangian–Eulerian (ALE) frame to accommodate the deformation of domain boundaries. An overlapping additive decomposition of velocity and pressure fields into coarse and fine scale components leads to coarse and fine scale mixed‐field problems. The problem governing fine scales is subjected to a further decomposition of the fine scale velocity into overlapping components termed as fine scales level I and level II. In turn, in the bottom‐up integration of scales, the model for level II fine scales serves to stabilize the problem governing level I fine scales, and model for level I fields yields the turbulence models. From the computational perspective, the coarse scales are represented in terms of the standard Lagrange shape functions, whereas level I and level II scales are represented via quadratic and fourth order polynomial bubbles, respectively. Because of the bubble functions approach employed in the consistently derived fine scale models, the resulting method is free of any embedded or tunable parameters. The proposed turbulence models share a common feature with the LES models in that the largest scales in the flow are numerically resolved, whereas the subgrid scales are modeled. The method is applied to flow around a plunging airfoil at Re = 40,000, and results are compared with experimental and numerical data published in the literature. Also presented are the results for the plunging airfoil at Re = 60,000 to show the robustness and range of applicability of the method. Copyright © 2013 John Wiley & Sons, Ltd.