Feature‐based adaptive mesh refinement for wingtip vortices

Feature‐based solution‐adaptive mesh refinement is an attractive strategy when it is known a priori that the resolution of certain key features is critical to achieving the objectives of a simulation. In this paper, we apply vortex characterization techniques, which are typically employed to visualize vortices, to identify regions of the computational domain for mesh refinement. We investigate different refinement strategies that are facilitated by these vortex characterization techniques to simulate the flow past a wing in a wind tunnel. Our results, which we compare with experimental data, indicate that it is necessary to refine the region within and near the vortex extent surface to obtain an accurate prediction. Application of the identified mesh refinement strategy also produced observed improvement in the results predicted for a spinning missile with deflected canards. Copyright © 2010 John Wiley & Sons, Ltd.     

On the use of anisotropic a posteriori error estimators for the adaptative solution of 3D inviscid compressible flows

This paper describes the use of an a posteriori error estimator to control anisotropic mesh adaptation for computing inviscid compressible flows. The a posteriori error estimator and the coupling strategy with an anisotropic remesher are first introduced. The mesh adaptation is controlled by a single‐parameter tolerance (TOL) in regions where the solution is regular, whereas a condition on the minimal element size h_min is enforced across solution discontinuities. This h_min condition is justified on the basis of an asymptotic analysis. The efficiency of the approach is tested with a supersonic flow over an aircraft. The evolution of a mesh adaptation/flow solution loop is shown, together with the influence of the parameters TOL and h_min. We verify numerically that the effect of varying h_min is concordant with the conclusions of the asymptotic analysis, giving hints on the selection of h_min with respect to TOL. Finally, we check that the results obtained with the a posteriori error estimator are at least as accurate as those obtained with anisotropic a priori error estimators. All the results presented can be obtained using a standard desktop computer, showing the efficiency of these adaptative methods. Copyright © 2008 John Wiley & Sons, Ltd.     

Adaptive mesh refinement and load balancing based on multi-level block-structured Cartesian mesh

We developed a framework for a distributed-memory parallel computer that enables dynamic data management for adaptive mesh refinement and load balancing. We employed simple data structure of the building cube method (BCM) where a computational domain is divided into multi-level cubic domains and each cube has the same number of grid points inside, realising a multi-level block-structured Cartesian mesh. Solution adaptive mesh refinement, which works efficiently with the help of the dynamic load balancing, was implemented by dividing cubes based on mesh refinement criteria. The framework was investigated with the Laplace equation in terms of adaptive mesh refinement, load balancing and the parallel efficiency. It was then applied to the incompressible Navier–Stokes equations to simulate a turbulent flow around a sphere. We considered wall-adaptive cube refinement where a non-dimensional wall distance y⁺ near the sphere is used for a criterion of mesh refinement. The result showed the load imbalance due to y⁺ adaptive mesh refinement was corrected by the present approach. To utilise the BCM framework more effectively, we also tested a cube-wise algorithm switching where an explicit and implicit time integration schemes are switched depending on the local Courant-Friedrichs-Lewy (CFL) condition in each cube.     

The use of dynamic grid adaptation algorithms for the modelling of flow around a circular cylinder in sub‐critical flow regime

In the present study a dynamic grid adaptation (DGA) algorithm is utilized for predicting flow around a circular cylinder in sub‐critical flow regime at a Reynolds number of 1.4×10⁵. The reason for adopting a DGA algorithm is the unsteadiness of the flow field which makes a conventional mesh inefficient. The concept being adopted is to concentrate mesh refinement in regions with high gradients and high turbulent viscosity, while in the region further downstream where the flow is fully developed a coarser mesh will develop and turbulence is modelled with the large eddy simulation (LES) turbulence model. The aim of the study is to present an appropriate variable for mesh refinement, which accomplishes a high rate of mesh refinement in the region with high gradients. The new variable is a product of the local mesh cell size and the rate of strain and includes two additional variables to allow control over the refinement behaviour. The results are compared with experimental data at the corresponding Reynolds number and also with numerical results obtained with conventional mesh. It is demonstrated that DGA algorithms can give results of a very high quality for a mesh that is significantly smaller than for a conventional mesh. Copyright © 2003 John Wiley & Sons, Ltd.     

Mesh adaptation for simulation of unsteady flow with moving immersed boundaries

In this work, an approach for performing mesh adaptation in the numerical simulation of two‐dimensional unsteady flow with moving immersed boundaries is presented. In each adaptation period, the mesh is refined in the regions where the solution evolves or the moving bodies pass and is unrefined in the regions where the phenomena or the bodies deviate. The flow field and the fluid–solid interface are recomputed on the adapted mesh. The adaptation indicator is defined according to the magnitude of the vorticity in the flow field. There is no lag between the adapted mesh and the computed solution, and the adaptation frequency can be controlled to reduce the errors due to the solution transferring between the old mesh and the new one. The preservation of conservation property is mandatory in long‐time scale simulations, so a P1‐conservative interpolation is used in the solution transferring. A nonboundary‐conforming method is employed to solve the flow equations. Therefore, the moving‐boundary flows can be simulated on a fixed mesh, and there is no need to update the mesh at each time step to follow the motion or the deformation of the solid boundary. To validate the present mesh adaptation method, we have simulated several unsteady flows over a circular cylinder stationary or with forced oscillation, a single self‐propelled swimming fish, and two fish swimming in the same or different directions. Copyright © 2012 John Wiley & Sons, Ltd.     

Interface‐fitted moving mesh method for axisymmetric two‐phase flow in microchannels

A boundary‐fitted moving mesh scheme is presented for the simulation of two‐phase flow in two‐dimensional and axisymmetric geometries. The incompressible Navier‐Stokes equations are solved using the finite element method, and the mini element is used to satisfy the inf‐sup condition. The interface between the phases is represented explicitly by an interface adapted mesh, thus allowing a sharp transition of the fluid properties. Surface tension is modelled as a volume force and is discretized in a consistent manner, thus allowing to obtain exact equilibrium (up to rounding errors) with the pressure gradient. This is demonstrated for a spherical droplet moving in a constant flow field. The curvature of the interface, required for the surface tension term, is efficiently computed with simple but very accurate geometric formulas. An adaptive moving mesh technique, where smoothing mesh velocities and remeshing are used to preserve the mesh quality, is developed and presented. Mesh refinement strategies, allowing tailoring of the refinement of the computational mesh, are also discussed. Accuracy and robustness of the present method are demonstrated on several validation test cases. The method is developed with the prospect of being applied to microfluidic flows and the simulation of microchannel evaporators used for electronics cooling. Therefore, the simulation results for the flow of a bubble in a microchannel are presented and compared to experimental data.     

Large-eddy simulation study of upstream boundary conditions influence upon a backward-facing step flowÉtude par simulation des grandes échelles de l'influence des conditions aux limites amont sur l'écoulement derrière une marche descendante

We use Large Eddy Simulation to investigate the influence of upstream boundary conditions on the development of a backward facing step flow. The first inlet condition consists of a mean turbulent boundary layer velocity profile perturbed by a white noise. The second relies upon a precursor calculation where the development of a quasi-temporal turbulent boundary layer is simulated. In this case, the quasi-longitudinal vortices in the upstream turbulent boundary-layer trigger the destabilization of the shear layer just behind the step, resulting in a shortening of the recirculation length and an increase of the characteristic frequency associated to the Kelvin–Helmholtz vortices. The mean flow and the characteristic frequencies of pressure fluctuations are strongly dependent of the upstream flow. It demonstrates the importance of realistic boundary conditions for the simulation of complex 3D flows or for flow control simulations. To cite this article: J.-L. Aider, A. Danet, C. R. Mecanique 334 (2006).     

Evaluation of adaptive mesh refinement and coarsening for the computation of compressible flows on unstructured meshes

The compressible gas flows of interest to aerospace applications often involve situations where shock and expansion waves are present. Decreasing the characteristic dimension of the computational cells in the vicinity of shock waves improves the quality of the computed flows. This reduction in size may be accomplished by the use of mesh adaption procedures. In this paper an analysis is presented of an adaptive mesh scheme developed for an unstructured mesh finite volume upwind computer code. This scheme is tailored to refine or coarsen the computational mesh where gradients of the flow properties are respectively high or low. The refinement and coarsening procedures are applied to the classical gas dynamic problems of the stabilization of shock waves by solid bodies. In particular, situations where oblique shock waves interact with an expansion fan and where bow shocks arise around solid bodies are considered. The effectiveness of the scheme in reducing the computational time, while increasing the solution accuracy, is assessed. It is shown that the refinement procedure alone leads to a number of computational cells which is 20% larger than when alternate passes of refinement and coarsening are used. Accordingly, a reduction of computational time of the same order of magnitude is obtained. Copyright © 2005 John Wiley & Sons, Ltd.     

Modeling particle sedimentation in viscous fluids with a coupled immersed boundary method and discrete element method

Numerical techniques have increasingly been used to model fluid–particle two-phase flows. Coupling the immersed boundary method (IBM) and discrete element method (DEM) is one promising approach for modeling particulate flows. In this study, IBM was coupled with DEM to improve the reliability and accuracy of IBM for determining the positions of particles during the sedimentation process within viscous fluids. The required ratio of the particle diameter to the grid size (D/dx) was determined by comparing the simulation results with the analytical solution and experimental data. A dynamic mesh refinement model was utilised in the IBM model to refine the computational fluid dynamics grid near the particles. In addition, an optimum coupling interval between the IBM and DEM models was determined based on the experimental results of a single particle sedimentation within silicon oil at a Reynolds number of 1.5. The experimental results and the analytical solution were then utilised to validate the IBM–DEM model at Reynolds numbers of 4.1, 11.6, and 31.9. Finally, the validated model was utilised to investigate the sedimentation process for more than one particle by modeling the drafting-kissing-tumbling process and the Boycott phenomenon. Benchmark tests showed that the IBM–DEM technique preserves the advantages of DEM for tracking a group of particles, while the IBM provides a reliable and accurate approach for modeling the particle–fluid interaction.     

The direct simulation Monte Carlo method using unstructured adaptive mesh and its application

The implementation of an adaptive mesh‐embedding (h‐refinement) scheme using unstructured grid in two‐dimensional direct simulation Monte Carlo (DSMC) method is reported. In this technique, local isotropic refinement is used to introduce new mesh where the local cell Knudsen number is less than some preset value. This simple scheme, however, has several severe consequences affecting the performance of the DSMC method. Thus, we have applied a technique to remove the hanging node, by introducing the an‐isotropic refinement in the interfacial cells between refined and non‐refined cells. Not only does this remedy increase a negligible amount of work, but it also removes all the difficulties presented in the originals scheme. We have tested the proposed scheme for argon gas in a high‐speed driven cavity flow. The results show an improved flow resolution as compared with that of un‐adaptive mesh. Finally, we have used triangular adaptive mesh to compute a near‐continuum gas flow, a hypersonic flow over a cylinder. The results show fairly good agreement with previous studies. In summary, the proposed simple mesh adaptation is very useful in computing rarefied gas flows, which involve both complicated geometry and highly non‐uniform density variations throughout the flow field. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical approximation of generalized Newtonian fluids using Powell–Sabin–Heindl elements: I. theoretical estimates

In this paper we consider the numerical approximation of steady and unsteady generalized Newtonian fluid flows using divergence free finite elements generated by the Powell–Sabin–Heindl elements. We derive a priori and a posteriori finite element error estimates and prove convergence of the method of successive approximations for the steady flow case. A priori error estimates of unsteady flows are also considered. These results provide a theoretical foundation and supporting numerical studies are to be provided in Part II. Copyright © 2003 John Wiley & Sons, Ltd.     

Application of adaptive mesh refinement method to complex flows in clean rooms

The adaptive mesh refinement (AMR) method is developed for three-dimensional turbulent complex flows in clean rooms using the finite volume method with a collocated grid arrangement. Clean rooms have many interesting and complex flow characteristics especially the secondary flows and the recirculation regions. The accurate numerical solution of the flows is important for the efficient design of clean rooms. The use of the conventional uniform grid requires such a high computational time and data storage capacity that they make computational fluid dynamics (CFD) less attractive for the design optimization. The AMR method is, therefore, applied by using the fine grid only in the required regions and using the coarse grid in the other regions. The velocity is chosen as the main parameter for the grid refinement because it is the most influential parameter in clean rooms. The results show that the present AMR method can reduce the computational time by eight times and the data storage requirement is only 37% of that using the conventional method, while the same order of accuracy can be maintained. The present AMR method is, therefore, proved to be a promising technique for solving three-dimensional turbulent complex flows in clean rooms.     

A Priori Assessment of Prediction Confidence for Data-Driven Turbulence Modeling

Although Reynolds-Averaged Navier–Stokes (RANS) equations are still the dominant tool for engineering design and analysis applications involving turbulent flows, standard RANS models are known to be unreliable in many flows of engineering relevance, including flows with separation, strong pressure gradients or mean flow curvature. With increasing amounts of 3-dimensional experimental data and high fidelity simulation data from Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS), data-driven turbulence modeling has become a promising approach to increase the predictive capability of RANS simulations. However, the prediction performance of data-driven models inevitably depends on the choices of training flows. This work aims to identify a quantitative measure for a priori estimation of prediction confidence in data-driven turbulence modeling. This measure represents the distance in feature space between the training flows and the flow to be predicted. Specifically, the Mahalanobis distance and the kernel density estimation (KDE) technique are used as metrics to quantify the distance between flow data sets in feature space. To examine the relationship between these two extrapolation metrics and the machine learning model prediction performance, the flow over periodic hills at Re = 10595 is used as test set and seven flows with different configurations are individually used as training sets. The results show that the prediction error of the Reynolds stress anisotropy is positively correlated with Mahalanobis distance and KDE distance, demonstrating that both extrapolation metrics can be used to estimate the prediction confidence a priori. A quantitative comparison using correlation coefficients shows that the Mahalanobis distance is less accurate in estimating the prediction confidence than KDE distance. The extrapolation metrics introduced in this work and the corresponding analysis provide an approach to aid in the choice of data source and to assess the prediction performance for data-driven turbulence modeling.     

The prediction of spatially periodic flows using a finite‐volume model

A periodic boundary condition has been developed that can be used in conjunction with a specified flow rate to produce accurate results in spatially periodic geometries. This condition is useful in situations where the flow rate is known, or more importantly, in cases where the pressure gradient is not known a priori, such as in countercurrent flows. Using the present condition, the flow rate is imposed at the inlet in terms of a bulk velocity, but the velocity field evolves as part of the solution. The condition is formulated to be suitable for both fixed and moving periodic domains. For the case of a moving domain, a correction is introduced to account for changes in the instantaneous velocity through the periodic edges. Under periodic conditions, these corrections integrate to zero over a complete (temporal) period. The new periodic condition is shown to produce accurate results for flat and wavy‐walled channels under both induced flow and countercurrent conditions. Copyright © 2003 John Wiley & Sons, Ltd.     

Stability of elliptical vortices from “Imperfect–Velocity–Impulse” diagrams

In 1875, Lord Kelvin proposed an energy-based argument for determining the stability of vortical flows. While the ideas underlying Kelvin’s argument are well established, their practical use has been the subject of extensive debate. In a forthcoming paper, the authors present a methodology, based on the construction of “Imperfect–Velocity–Impulse” (IVI) diagrams, which represents a rigorous and practical implementation of Kelvin’s argument for determining the stability of inviscid flows. In this work, we describe in detail the use of the theory by considering an example involving a well-studied classical flow, namely the family of elliptical vortices discovered by Kirchhoff. By constructing the IVI diagram for this family of vortices, we detect the first three bifurcations (which are found to be associated with perturbations of azimuthal wavenumber m = 3, 4 and 5). Examination of the IVI diagram indicates that each of these bifurcations contributes an additional unstable mode to the original family; the stability properties of the bifurcated branches are also determined. By using a novel numerical approach, we proceed to explore each of the bifurcated branches in its entirety. While the locations of the changes of stability obtained from the IVI diagram approach turn out to match precisely classical results from linear analysis, the stability properties of the bifurcated branches are presented here for the first time. In addition, it appears that the m = 3, 5 branches had not been computed in their entirety before. In summary, the work presented here outlines a new approach representing a rigorous implementation of Kelvin’s argument. With reference to the Kirchhoff elliptical vortices, this method is shown to be effective and reliable.     

An adaptive moving finite element method for steady low Mach number compressible combustion problems

This work surveys an r-adaptive moving mesh finite element method for the numerical solution of premixed laminar flame problems. Since the model of chemically reacting flow involves many different modes with diverse length scales, the computation of such a problem is often extremely time-consuming. Importantly, to capture the significant characteristics of the flame structure when using detailed chemistry, a much more stringent requirement on the spatial resolution of the interior layers of some intermediate species is necessary. Here, we propose a moving mesh method in which the mesh is obtained from the solution of so-called moving mesh partial differential equations. Such equations result from the variational formulation of a minimization problem for a given target functional that characterizes the inherent difficulty in the numerical approximation of the underlying physical equations. Adaptive mesh movement has emerged as an area of intense research in mesh adaptation in the last decade. With this approach, points are only allowed to be shifted in space leaving the topology of the grid unchanged. In contrast to methods with local refinement, data structure hence is unchanged and load balancing is not an issue as grid points remain on the processor where they are. We will demonstrate the high potential of moving mesh methods for effectively optimizing the distribution of grid points to reach the required resolution for chemically reacting flows with extremely thin boundary layers.     

Validation of adaptive unstructured hexahedral mesh computations of flow around a wind turbine airfoil

In this paper we investigate local adaptive refinement of unstructured hexahedral meshes for computations of the flow around the DU91 wind turbine airfoil. This is a 25% thick airfoil, found at the mid‐span section of a wind turbine blade. Wind turbine applications typically involve unsteady flows due to changes in the angle of attack and to unsteady flow separation at high angles of attack. In order to obtain reasonably accurate results for all these conditions one should use a mesh which is refined in many regions, which is not computationally efficient. Our solution is to apply an automated mesh adaptation technique. In this paper we test an adaptive refinement strategy developed for unstructured hexahedral meshes for steady flow conditions. The automated mesh adaptation is based on local flow sensors for pressure, velocity, density or a combination of these flow variables. This way the mesh is refined only in those regions necessary for high accuracy, retaining computational efficiency. A validation study is performed for two cases: attached flow at an angle of 6° and separated flow at 12°. The results obtained using our adaptive mesh strategy are compared with experimental data and with results obtained with an equally sized non‐adapted mesh. From these computations it can be concluded that for a given computing time, adapted meshes result in solutions closer to the experimental data compared to non‐adapted meshes for attached flow. Finally, we show results for unsteady computations. Copyright © 2005 John Wiley & Sons, Ltd.     

Compound wall treatment with low‐Re turbulence model

A generalized treatment for the wall boundary conditions relating to turbulent flows is developed that blends the integration to a solid wall with wall functions. The blending function ensures a smooth transition between the viscous and turbulent regions. An improved low Reynolds number k?ε model is coupled with the proposed compound wall treatment to determine the turbulence field. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. Computations with fine and coarse meshes of a few flow cases yield appreciably good agreement with the direct numerical simulation and experimental data. The method is recommended for computing the complex flows where computational grids cannot satisfy a priori the prerequisites of viscous/turbulence regions. Copyright © 2011 John Wiley & Sons, Ltd.     

A Mesh Adaptation Strategy to Predict Pressure Losses in LES of Swirled Flows

Large-Eddy Simulation (LES) has become a potent tool to investigate instabilities in swirl flows even for complex, industrial geometries. However, the accurate prediction of pressure losses on these complex flows remains difficult. The paper identifies localised near-wall resolution issues as an important factor to improve accuracy and proposes a solution with an adaptive mesh h-refinement strategy relying on the tetrahedral fully automatic MMG3D library of Dapogny et al. (J. Comput. Phys. 262, 358-378, 2014) using a novel sensor based on the dissipation of kinetic energy. Using a joint experimental and numerical LES study, the methodology is first validated on a simple diaphragm flow before to be applied on a swirler with two counter-rotating passages. The results demonstrate that the new sensor and adaptation approach can effectively produce the desired local mesh refinement to match the target losses, measured experimentally. Results shows that the accuracy of pressure losses prediction is mainly controlled by the mesh quality and density in the swirler passages. The refinement also improves the computed velocity and turbulence profiles at the swirler outlet, compared to PIV results. The significant improvement of results confirms that the sensor is able to identify the relevant physics of turbulent flows that is essential for the overall accuracy of LES. Finally, in the appendix, an additional comparison of the sensor fields on tetrahedral and hexahedral meshes demonstrates that the methodology is broadly applicable to all mesh types.     

Velocity field and parametric analysis of a subsonic,medium-Reynolds number cavity flow

Cavity flows are a class of flows bounded by material structures, where a recirculation region is present, and they are found in many practical applications. In the present study, the interaction between a boundary layer and an open parallelepipedic cavity develops a Kelvin–Helmholtz-like instability coupled with the cavity recirculation. PIV measurements of the flow are carried out in two orthogonal planes inside the cavity, for different aspect ratios, incompressible flow conditions, and Reynolds numbers in the range 1,900–12,000. Mean velocity and second-order moments of velocity fluctuations reveal the flow morphology. For particular conditions, centrifugal instabilities appear that are induced by flow curvature due to wall confinement. The use of an identification criterion indicates the presence of pairs of counter-rotating vortices winded around the recirculation. A parametric analysis is conducted, and the inviscid Rayleigh discriminant provides the potentially unstable flow regions inside the cavity. Finally, a stability parameter considering the ratio between centrifugal destabilizing effects and stabilizing viscous effects is carried out and gives thresholds for the emergence of the centrifugal instability. The study draws to an end with a comparison with a well-documented lid-driven cavity flow.