High‐order detached eddy simulation,zonal LES and URANS of cavity and labyrinth seal flows

Hybrid numerical large eddy simulation (NLES), detached eddy simulation (DES) and URANS methods are assessed on a cavity and a labyrinth seal geometry. A high sixth‐order discretization scheme is used and is validated using the test case of a two‐dimensional vortex. The hybrid approach adopts a new blending function. For the URANS simulations, the flow within the cavity remains steady, and the results show significant variation between models. Surprisingly, low levels of resolved turbulence are observed in the cavity for the DES simulation, and the cavity shear layer remains two dimensional. The hybrid RANS–NLES approach does not suffer from this trait. For the labyrinth seal, both the URANS and DES approaches give low levels of resolved turbulence. The zonal Hamilton–Jacobi approach on the other had given significantly more resolved content. Both DES and hybrid RANS–NLES give good agreement with the experimentally measured velocity profiles. Again, there is significant variation between the URANS models, and swirl velocities are overpredicted. Copyright © 2013 John Wiley & Sons, Ltd.     

Full coverage film cooling using compound angle

An experimental and numerical study is carried out on a cooling film issuing from a multiperforated wall of a simplified combustor. The objectives of this work are to achieve a better understanding of the dynamics of the film and to construct an experimental database on a simplified geometry in order to test numerical models. A parametric study of film cooling efficiency based on the direction of the cooling air injection is presented and shows that a swirling injection greatly enhances the cooling efficiency. As accounting for multiperforated walls in numerical simulations cannot be done at the jets scale because of computing resources, in this article are presented RANS computations performed using a uniform boundary condition to provide the injection of coolant. Two injection models are applied on this boundary and numerical results are compared to experimental data in the recovery region. The standard model is shown to be totally inappropriate while the multiperforation model delivers promising results although some weaknesses appear very close to the wall. To cite this article: B. Michel et al., C. R. Mecanique 337 (2009).     

Finite‐difference 4th‐order compact scheme for the direct numerical simulation of instabilities of shear layers

An algorithm based on the 4th‐order finite‐difference compact scheme is developed and applied in the direct numerical simulations of instabilities of channel flow. The algorithm is illustrated in the context of stream function formulation that leads to field equation involving 4th‐order spatial derivatives. The finite‐difference discretization in the wall‐normal direction uses five arbitrarily spaced points. The discretization coefficients are determined numerically, providing a large degree of flexibility for grid selection. The Fourier expansions are used in the streamwise direction. A hybrid Runge–Kutta/Crank–Nicholson low‐storage scheme is applied for the time discretization. Accuracy tests demonstrate that the algorithm does deliver the 4th‐order accuracy. The algorithm has been used to simulate the natural instability processes in channel flow as well as processes occurring when the flow is spatially modulated using wall transpiration. Extensions to three‐dimensional situations are suggested. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical prediction of turbulent flows using Reynolds‐averaged Navier–Stokes and large‐eddy simulation with uncertain inflow conditions

Numerous comparisons between Reynolds‐averaged Navier–Stokes (RANS) and large‐eddy simulation (LES) modeling have already been performed for a large variety of turbulent flows in the context of fully deterministic flows, that is, with fixed flow and model parameters. More recently, RANS and LES have been separately assessed in conjunction with stochastic flow and/or model parameters. The present paper performs a comparison of the RANS k ? ε model and the LES dynamic Smagorinsky model for turbulent flow in a pipe geometry subject to uncertain inflow conditions. The influence of the experimental uncertainties on the computed flow is analyzed using a non‐intrusive polynomial chaos approach for two flow configurations (with or without swirl). Measured quantities including an estimation of the measurement error are then compared with the statistical representation (mean value and variance) of their RANS and LES numerical approximations in order to check whether experiment/simulation discrepancies can be explained within the uncertainty inherent to the studied configuration. The statistics of the RANS prediction are found in poor agreement with experimental results when the flow is characterized by a strong swirl, whereas the computationally more expensive LES prediction remains statistically well inside the measurement intervals for the key flow quantities.Copyright © 2012 John Wiley & Sons, Ltd.     

A DES method applied to a Backward Facing Step reactive flow

Hybrid Reynolds Averaged Navier Stokes–Large Eddy Simulation is a trend which is becoming of common use in aerodynamics but has seldom been employed to simulate reactive flows. Such methods, like the Delayed Detached Eddy Simulation (DDES) presented in this article, have been created to treat near wall flows with a RANS approach while switching to LES in the separated flow region. It is indeed an affordable solution to simulate complex and unsteady compressible flows and to have access to accurate skin friction and wall thermal fluxes. In order to validate this technique in combustion, we chose a simple and well documented Backward Facing Step combustor. To account for turbulent combustion a Dynamic Thickened Flame was used. The results obtained on this case show a good agreement with the experimental database and are of the same quality as LES in the separated region for both inert and reactive flows. To cite this article: B. Sainte-Rose et al., C. R. Mecanique 337 (2009).     

Efficient retrieval of the thermo‐acoustic flame transfer function from a linearized CFD simulation of a turbulent flame

The stability of thermo‐acoustic pressure oscillations in a lean premixed methane‐fired generic gas turbine combustor is investigated. A key element in predicting the acoustically unstable operating conditions of the combustor is the flame transfer function. This function represents the dynamic relationship between a fluctuation in the combustor inlet conditions and the flame's acoustic response. A transient numerical experiment involving spectral analysis in computational fluid dynamics (CFD) is usually conducted to predict the flame transfer function. An important drawback of this spectral method application to numerical simulations is the required computational effort. A much faster and more accurate method to calculate the transfer function is derived in this paper by using a most important basic assumption: the fluctuations must be small enough for the system to behave linear. This alternative method, which is called the linear coefficient method, uses a linear representation of the unsteady equations describing the CFD problem. This linearization is applied around a steady‐state solution of the problem, where it can consequently describe the dynamics of the system. Finally, the flame transfer function can be calculated from this linear representation. The advantage of this approach is that one only needs a steady‐state solution and linearization of the unsteady equations for calculating a dynamic transfer function, i.e. no time‐consuming transient simulations are necessary anymore. Nevertheless, as a consequence of the large number of degrees of freedom in a CFD problem, an extra order reduction step needs to be performed prior to calculating the transfer function from the linear representation. Still, the linear coefficient method shows a significant gain in both speed and accuracy when calculating the transfer function from the linear representation as compared to a spectral analysis‐based calculation. Hence, this method gives a major improvement to the application of the flame transfer function as a thermo‐acoustic design tool. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical solutions of systems with (p,δ)‐structure using local discontinuous Galerkin finite element methods

In this paper, we present LDG methods for systems with (p,δ)‐structure. The unknown gradient and the nonlinear diffusivity function are introduced as auxiliary variables and the original (p,δ) system is decomposed into a first‐order system. Every equation of the produced first‐order system is discretized in the discontinuous Galerkin framework, where two different nonlinear viscous numerical fluxes are implemented. An a priori bound for a simplified problem is derived. The ODE system resulting from the LDG discretization is solved by diagonal implicit Runge–Kutta methods. The nonlinear system of algebraic equations with unknowns the intermediate solutions of the Runge–Kutta cycle is solved using Newton and Picard iterative methodology. The performance of the two nonlinear solvers is compared with simple test problems. Numerical tests concerning problems with exact solutions are performed in order to validate the theoretical spatial accuracy of the proposed method. Further, more realistic numerical examples are solved in domains with non‐smooth boundary to test the efficiency of the method. Copyright © 2014 John Wiley & Sons, Ltd.     

Computational and experimental investigation of flow over a transient pitching hydrofoil

The present study is developed within the framework of marine structure design operating in transient regimes. It deals with an experimental and numerical investigation of the time–space distribution of the wall-pressure field on a NACA66 hydrofoil undergoing a transient up-and-down pitching motion from 0° to 15° at four pitching velocities and a Reynolds number Re = 0.75 × 10⁶. The experimental investigation is performed using an array of wall-pressure transducers located on the suction side and by means of time–frequency analysis and Empirical Modal Decomposition method. The numerical study is conducted for the same flow conditions. It is based on a 2D RANS code including mesh reconstruction and an ALE formulation in order to take into account the foil rotation and the tunnel walls. Due to the moderate Reynolds number, a laminar to turbulent transition model was also activated. For the operating flow conditions of the study, experimental and numerical flow analysis revealed that the flow experiences complex boundary layer events as leading-edge laminar separation bubble, laminar to turbulent transition, trailing-edge separation and flow detachment at stall. Although the flow is relatively complex, the calculated wall pressure shows a quite good agreement with the experiment provided that the mesh resolution and the temporal discretization are carefully selected depending on the pitching velocity. It is particularly shown that the general trend of the wall pressure (low frequency) is rather well predicted for the four pitching velocities with for instance a net inflection of the wall pressure when transition occurs. The inflection zone is reduced as the pitching velocity increases and tends to disappear for the highest pitching velocity. Conversely, high frequency wall-pressure fluctuations observed experimentally are not captured by the RANS model. Based on the good agreement with experiment, the model is then used to investigate the effects of the pitching velocity on boundary layer events and on hydrodynamic loadings. It is shown that increasing the pitching velocity tends to delay the laminar-to-turbulence transition and even to suppress it for the highest pitching velocity during the pitch-up motion. It induces also an increase of the stall angle (compared to quasi-static one) and an increase of the hysteresis effect during pitch-down motion resulting to a significant increase of the hydrodynamic loading.     

Surface‐sampled simulations of turbulent flow at high Reynolds number

A new approach to turbulence simulation, based on a combination of large eddy simulation (LES) for the whole flow and an array of non–space‐filling quasi‐direct numerical simulations (QDNS), which sample the response of near‐wall turbulence to large‐scale forcing, is proposed and evaluated. The technique overcomes some of the cost limitations of turbulence simulation, since the main flow is treated with a coarse‐grid LES, with the equivalent of wall functions supplied by the near‐wall sampled QDNS. Two cases are tested, at friction Reynolds number Re_τ=4200 and 20000. The total grid point count for the first case is less than half a million and less than 2 million for the second case, with the calculations only requiring a desktop computer. A good agreement with published direct numerical simulation (DNS) is found at Re_τ=4200, both in the mean velocity profile and the streamwise velocity fluctuation statistics, which correctly show a substantial increase in near‐wall turbulence levels due to a modulation of near‐wall streaks by large‐scale structures. The trend continues at Re_τ=20000, in agreement with experiment, which represents one of the major achievements of the new approach. A number of detailed aspects of the model, including numerical resolution, LES‐QDNS coupling strategy and subgrid model are explored. A low level of grid sensitivity is demonstrated for both the QDNS and LES aspects. Since the method does not assume a law of the wall, it can in principle be applied to flows that are out of equilibrium.     

A local mesh refinement approach for large‐eddy simulations of turbulent flows

In this paper, a local mesh refinement (LMR) scheme on Cartesian grids for large‐eddy simulations is presented. The approach improves the calculation of ghost cell pressures and velocities and combines LMR with high‐order interpolation schemes at the LMR interface and throughout the rest of the computational domain to ensure smooth and accurate transition of variables between grids of different resolution. The approach is validated for turbulent channel flow and flow over a matrix of wall‐mounted cubes for which reliable numerical and experimental data are available. Comparisons of predicted first‐order and second‐order turbulence statistics with the validation data demonstrated a convincing agreement. Importantly, it is shown that mean streamwise velocities and fluctuating turbulence quantities transition smoothly across coarse‐to‐fine and fine‐to‐coarse interfaces. © 2016 The Authors International Journal for Numerical Methods in Fluids Published by John Wiley & Sons Ltd     

Investigation of the sensitivity of turbulent closures and coupling of hybrid RANS‐LES models for predicting flow fields with separation and reattachment

Hybrid models have found widespread applications for simulation of wall‐bounded flows at high Reynolds numbers. Typically, these models employ Reynolds‐averaged Navier–Stokes (RANS) and large eddy simulation (LES) in the near‐body and off‐body regions, respectively. A number of coupling strategies between the RANS and LES regions have been proposed, tested, and applied in the literature with varying degree of success. Linear eddy‐viscosity models (LEVM) are often used for the closure of turbulent stress tensor in RANS and LES regions. LEVM incorrectly predicts the anisotropy of Reynolds normal stress at the RANS‐LES interface region. To overcome this issue, use of non‐linear eddy‐viscosity models (NLEVM) have started receiving attention. In this study, a generic non‐linear blended modeling framework for performing hybrid simulations is proposed. Flow over the periodic hills is used as the test case for model evaluation. This case is chosen due to complex flow physics with simplified geometry. Analysis of the simulations suggests that the non‐linear hybrid models show a better performance than linear hybrid models. It is also observed that the non‐linear closures are less sensitive to the RANS‐LES coupling and grid resolution. Copyright © 2016 John Wiley & Sons, Ltd.     

A local domain‐free discretization method to simulate three‐dimensional compressible inviscid flows

In this paper, the domain‐free discretization method (DFD) is extended to simulate the three‐dimensional compressible inviscid flows governed by Euler equations. The discretization strategy of DFD is that the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior‐dependent points are updated at each time step by extrapolation along the wall normal direction in conjunction with the wall boundary conditions and the simplified momentum equation in the vicinity of the wall. Spatial discretization is achieved with the help of the finite element Galerkin approximation. The concept of ‘osculating plane’ is adopted, with which the local DFD can be easily implemented for the three‐dimensional case. Geometry‐adaptive tetrahedral mesh is employed for three‐dimensional calculations. Finally, we validate the DFD method for three‐dimensional compressible inviscid flow simulations by computing transonic flows over the ONERA M6 wing. Comparison with the reference experimental data and numerical results on boundary‐conforming grid was displayed and the results show that the present DFD results compare very well with the reference data. Copyright © 2009 John Wiley & Sons, Ltd.     

Injection coupling with high amplitude transverse modes: Experimentation and simulation

High frequency combustion instabilities have technical importance in the design of liquid rocket engines. These phenomena involve a strong coupling between transverse acoustic modes and combustion. They are currently being investigated by combining experimentation and numerical simulations. On the experimental level, the coupling is examined in a model scale system featuring a multiple injector combustor (MIC) comprising five coaxial injectors fed with liquid oxygen and gaseous methane. This system is equipped with a novel VHAM actuator (Very High Amplitude Modulator) which comprises two nozzles and a rotating toothed wheel blocking the nozzles in an alternate fashion. This device was designed to obtain the highest possible levels of transverse oscillation in the MIC. After a brief review of the VHAM, this article reports cold flow experiments using this modulator. Velocity maps obtained under resonant conditions using the VHAM are examined at different instants during a cycle of oscillation. Experimental data are compared with numerical pressure and velocity fields obtained from an acoustic solver. The good agreement observed in the nozzle vicinity indicates that numerical simulations can be used to analyze the complex flow field generated by the VHAM. To cite this article: Y. Mery et al., C. R. Mecanique 337 (2009).     

Numerical simulation of transpiration cooling through porous material

Transpiration cooling using ceramic matrix composite materials is an innovative concept for cooling rocket thrust chambers. The coolant (air) is driven through the porous material by a pressure difference between the coolant reservoir and the turbulent hot gas flow. The effectiveness of such cooling strategies relies on a proper choice of the involved process parameters such as injection pressure, blowing ratios, and material structure parameters, to name only a few. In view of the limited experimental access to the subtle processes occurring at the interface between hot gas flow and porous medium, reliable and accurate simulations become an increasingly important design tool. In order to facilitate such numerical simulations for a carbon/carbon material mounted in the side wall of a hot gas channel that are able to capture a spatially varying interplay between the hot gas flow and the coolant at the interface, we formulate a model for the porous medium flow of Darcy–Forchheimer type. A finite‐element solver for the corresponding porous medium flow is presented and coupled with a finite‐volume solver for the compressible Reynolds‐averaged Navier–Stokes equations. The two‐dimensional and three‐dimensional results at Mach number Ma = 0.5 and hot gas temperature T_HG=540 K for different blowing ratios are compared with experimental data. Copyright © 2014 John Wiley & Sons, Ltd.     

Modeling of Internal Reacting Flows and External Static Stall Flows Using RANS and PRNS

The reacting flow in a research lean direct injection (LDI) hydrogen combustor and the static stall of NACA0012 airfoil were simulated using both Reynolds averaged Navier–Stokes (RANS) and partially resolved numerical simulation (PRNS) approaches. The concept and the main features of the PRNS approach are briefly described. The PRNS basic equations are grid independent or grid invariant; the subscale models are a dynamic equation system. We consider PRNS as an engineering tool for the very large eddy simulation of complex turbulent flows. Two CFD codes, NCC and Wind-US, with two different subscale models (i.e. two- and one-transport equation models, respectively) are used in the presented PRNS simulations. Based on the comparisons with available experimental data, the numerical results indicate that the PRNS subscale models seem to be able to capture important large scale turbulent structures and to improve the quality of numerical simulations while keeping a relatively low cost comparable to the unsteady RANS simulations.     

Convergence issues in using high‐resolution schemes and lower–upper symmetric Gauss–Seidel method for steady shock‐induced combustion problems

This paper reports numerical convergence study for simulations of steady shock‐induced combustion problems with high‐resolution shock‐capturing schemes. Five typical schemes are used: the Roe flux‐based monotone upstream‐centered scheme for conservation laws (MUSCL) and weighted essentially non‐oscillatory (WENO) schemes, the Lax–Friedrichs splitting‐based non‐oscillatory no‐free parameter dissipative (NND) and WENO schemes, and the Harten–Yee upwind total variation diminishing (TVD) scheme. These schemes are implemented with the finite volume discretization on structured quadrilateral meshes in dimension‐by‐dimension way and the lower–upper symmetric Gauss–Seidel (LU–SGS) relaxation method for solving the axisymmetric multispecies reactive Navier–Stokes equations. Comparison of iterative convergence between different schemes has been made using supersonic combustion flows around a spherical projectile with Mach numbers M = 3.55 and 6.46 and a ram accelerator with M = 6.7. These test cases were regarded as steady combustion problems in literature. Calculations on gradually refined meshes show that the second‐order NND, MUSCL, and TVD schemes can converge well to steady states from coarse through fine meshes for M = 3.55 case in which shock and combustion fronts are separate, whereas the (nominally) fifth‐order WENO schemes can only converge to some residual level. More interestingly, the numerical results show that all the schemes do not converge to steady‐state solutions for M = 6.46 in the spherical projectile and M = 6.7 in the ram accelerator cases on fine meshes although they all converge on coarser meshes or on fine meshes without chemical reactions. The result is based on the particular preconditioner of LU–SGS scheme. Possible reasons for the nonconvergence in reactive flow simulation are discussed.Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical simulation of a single bubble by compressible two‐phase fluids

The present work deals with the numerical investigation of a collapsing bubble in a liquid–gas fluid, which is modeled as a single compressible medium. The medium is characterized by the stiffened gas law using different material parameters for the two phases. For the discretization of the stiffened gas model, the approach of Saurel and Abgrall is employed where the flow equations, here the Euler equations, for the conserved quantities are approximated by a finite volume scheme, and an upwind discretization is used for the non‐conservative transport equations of the pressure law coefficients. The original first‐order discretization is extended to higher order applying second‐order ENO reconstruction to the primitive variables. The derivation of the non‐conservative upwind discretization for the phase indicator, here the gas fraction, is presented for arbitrary unstructured grids. The efficiency of the numerical scheme is significantly improved by employing local grid adaptation. For this purpose, multiscale‐based grid adaptation is used in combination with a multilevel time stepping strategy to avoid small time steps for coarse cells. The resulting numerical scheme is then applied to the numerical investigation of the 2‐D axisymmetric collapse of a gas bubble in a free flow field and near to a rigid wall. The numerical investigation predicts physical features such as bubble collapse, bubble splitting and the formation of a liquid jet that can be observed in experiments with laser‐induced cavitation bubbles. Opposite to the experiments, the computations reveal insight to the state inside the bubble clearly indicating that these features are caused by the acceleration of the gas due to shock wave focusing and reflection as well as wave interaction processes. While incompressible models have been used to provide useful predictions on the change of the bubble shape of a collapsing bubble near a solid boundary, we wish to study the effects of shock wave emissions into the ambient liquid on the bubble collapse, a phenomenon that may not be captured using an incompressible fluid model. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical simulation of free‐surface flow using the level‐set method with global mass correction

A new numerical method that couples the incompressible Navier–Stokes equations with the global mass correction level‐set method for simulating fluid problems with free surfaces and interfaces is presented in this paper. The finite volume method is used to discretize Navier–Stokes equations with the two‐step projection method on a staggered Cartesian grid. The free‐surface flow problem is solved on a fixed grid in which the free surface is captured by the zero level set. Mass conservation is improved significantly by applying a global mass correction scheme, in a novel combination with third‐order essentially non‐oscillatory schemes and a five stage Runge–Kutta method, to accomplish advection and re‐distancing of the level‐set function. The coupled solver is applied to simulate interface change and flow field in four benchmark test cases: (1) shear flow; (2) dam break; (3) travelling and reflection of solitary wave and (4) solitary wave over a submerged object. The computational results are in excellent agreement with theoretical predictions, experimental data and previous numerical simulations using a RANS‐VOF method. The simulations reveal some interesting free‐surface phenomena such as the free‐surface vortices, air entrapment and wave deformation over a submerged object. Copyright © 2009 John Wiley & Sons, Ltd.     

Calculation of turbulent fluid flow and heat transfer in ducts by a full Reynolds stress model

A computational method has been developed to predict the turbulent Reynolds stresses and turbulent heat fluxes in ducts by different turbulence models. The turbulent Reynolds stresses and other turbulent flow quantities are predicted with a full Reynolds stress model (RSM). The turbulent heat fluxes are modelled by a SED concept, the GGDH and the WET methods. Two wall functions are used, one for the velocity field and one for the temperature field. All the models are implemented for an arbitrary three‐dimensional channel. Fully developed condition is achieved by imposing cyclic boundary conditions in the main flow direction. The numerical approach is based on the finite volume technique with a non‐staggered grid arrangement. The pressure–velocity coupling is handled by using the SIMPLEC‐algorithm. The convective terms are treated by the van Leer scheme while the diffusive terms are handled by the central‐difference scheme. The hybrid scheme is used for solving the ε equation. The secondary flow generation using the RSM model is compared with a non‐linear k–ε model (non‐linear eddy viscosity model). The overall comparison between the models is presented in terms of the friction factor and Nusselt number. Copyright © 2003 John Wiley & Sons, Ltd.     

A FE‐based algorithm for the inverse natural convection problem

Several numerical algorithms for solving inverse natural convection problems are revisited and studied. Our aim is to identify the unknown strength of a time‐varying heat source via a set of coupled nonlinear partial differential equations obtained by the so‐called finite element consistent splitting scheme (CSS) in order to get a good approximation of the unknown heat source from both the measured data and model results, by minimizing a functional that measures discrepancies between model and measured data. Viewed as an optimization problem, the solutions are obtained by means of the conjugate gradient method. A second‐order CSS in time involving the direct problem, the adjoint problem, the sensitivity problem and a system of sensitivity functions is used in order to enhance the numerical accuracy obtained for the unknown heat source function. A spatial discretization of all field equations is implemented using equal‐order and mixed finite element methods. Numerical experiments validate the proposed optimization algorithms that are in good agreement with the existing results. Copyright © 2010 John Wiley & Sons, Ltd.