The term “immersed boundary method” was first used in reference to a method developed by Peskin (1972) to simulate cardiac mechanics and associated blood flow. The distinguishing feature of this method was that the entire simulation was carried out on a Cartesian grid, which did not conform to the geometry of the heart, and a novel procedure was formulated for imposing the effect of the immersed boundary (IB) on the flow. Since Peskin introduced this method, numerous modifications and refinements have been proposed and a number of variants of this approach now exist. In addition, there is another class of methods, usually referred to as “Cartesian grid methods,” which were originally developed for simulating inviscid flows with complex embedded solid boundaries on Cartesian grids (Berger & Aftosmis 1998, Clarke et al. 1986, Zeeuw & Powell 1991). These methods have been extended to simulate unsteady viscous flows (Udaykumar et al. 1996, Ye et al. 1999) and thus have capabilities similar to those of IB methods. In this review, we use the term immersed boundary (IB) method to encompass all such methods that simulate viscous flows with immersed (or embedded) boundaries on grids that do not conform to the shape of these boundaries. Furthermore, this review focuses mainly on IB methods for flows with immersed solid boundaries. Application of these and related methods to problems with liquid-liquid and liquid-gas boundaries was covered in previous reviews by Anderson et al. (1998) and Scardovelli & Zaleski (1999). Consider the simulation of flow past a solid body shown in Figure 1a. The conventional approach to this would employ structured or unstructured grids that conform to the body. Generating these grids proceeds in two sequential steps. First, a surface grid covering the boundaries b is generated. This is then used as a boundary condition to generate a grid in the volume f occupied by the fluid. If a finite-difference method is employed on a structured grid, then the differential form of the governing equations is transformed to a curvilinear coordinate system aligned with the grid lines (Ferziger & Peric 1996). Because the grid conforms to the surface of the body, the transformed equations can then be discretized in the

/pdf/immersed-boundary-methods-4id27kda4o.pdf

Immersed boundary methods

Detached-eddy simulation (DES) was first proposed in 1997 and first used in 1999, so its full history can be surveyed. A DES community has formed, with adepts and critics, as well as new branches. The initial motivation of high–Reynolds number, massively separated flows remains, for which DES is convincingly more capable presently than either unsteady Reynolds-averaged Navier-Stokes (RANS) or large-eddy simulation (LES). This review discusses compelling examples, noting the visual and quantitative success of DES. Its principal weakness is its response to ambiguous grids, in which the wall-parallel grid spacing is of the order of the boundary-layer thickness. In some situations, DES on a given grid is then less accurate than RANS on the same grid or DES on a coarser grid. Partial remedies have been found, yet dealing with thickening boundary layers and shallow separation bubbles is a central challenge. The nonmonotonic response of DES to grid refinement is disturbing to most observers, as is the absence of...

Detached-Eddy Simulation

An immersed-boundary finite-volume method for simulations of flow in complex geometries

Dynamics and stability of lean-premixed swirl-stabilized combustion

Large-eddy simulation (LES) of turbulent combustion is a relatively new research field. Much research has been carried out over the past years, but to realize the full predictive potential of combustion LES, many fundamental questions still have to be addressed, and common practices of LES of nonreacting flows revisited. The focus of the present review is to highlight the fundamental differences between Reynolds-averaged Navier-Stokes (RANS) and LES combustion models for nonpremixed and premixed turbulent combustion, to identify some of the open questions and modeling issues for LES, and to provide future perspectives.

/pdf/large-eddy-simulation-of-turbulent-combustion-39umcykkan.pdf

Large-eddy simulation of turbulent combustion

A numerical method for large-eddy simulation in complex geometries

The flow field around a sphere in an uniform flow has been analyzed numerically for conditions corresponding to the subcritical (laminar separation) and supercritical (turbulent separation) regimes spanning a wide range of Reynolds numbers (104–106). Particular attention has been devoted to assessing predictions of the pressure distribution, skin friction, and drag as well as to understanding the changes in the wake organization and vortex dynamics with the Reynolds number. The unsteady turbulent flow is computed using detached-eddy simulation, a hybrid approach that has Reynolds-averaged Navier–Stokes behavior near the wall and becomes a large eddy simulation in the regions away from solid surfaces. For both the subcritical and supercritical solutions, the agreement with experimental measurements for the mean drag and pressure distribution over the sphere is adequate; differences in skin friction exist due to the simplistic treatment of the attached boundary layers in the computations. Improved agreement in the skin-friction distribution is obtained for the supercritical flows in which boundary layer transition is fixed at the position observed in experiments conducted at the same Reynolds numbers. For the subcritical flows the Strouhal number, St, associated with the large-scale shedding is predicted at St∼0.195 along with a higher frequency component associated with the development of the Kelvin–Helmholtz instabilities in the detached shear layers. If in the subcritical regime the wake assumes a helical-like form due to the shedding of hairpin-like vortices at different azimuthal angles, in the supercritical regime the wake structure is characterized by “regular” shedding of hairpin-like vortices at approximately the same azimuthal angle and at a much higher frequency (St∼1.3) that is practically independent of the Reynolds number and not sensitive to the position of laminar-to-turbulent transition.

Numerical investigations of flow over a sphere in the subcritical and supercritical regimes

[1] River confluences are complex hydrodynamic environments where convergence of incoming flows produces complicated patterns of fluid motion, including the development of large-scale turbulence structures. Accurately simulating confluence hydrodynamics represents a considerable challenge for numerical modeling of river flows. This study uses an eddy-resolving numerical model to simulate the mean flow and large-scale turbulence structure at an asymmetrical river confluence with a concordant bed when the momentum ratio between the two incoming streams is close to 1. Results of the simulation are compared with field data on mean flow and turbulence structure. The simulation shows that the mixing interface is populated by quasi-two-dimensional eddies. Successive eddies have opposing senses of rotation. The mixing layer structure resembles that of a wake behind a bluff body (wake mode). Strong streamwise-oriented vortical (SOV) cells form on both sides of the mixing layer, a finding consistent with patterns inferred from the field data. The predicted mean flow fields show that flow curvature has an important influence on streamwise variation of circulation within the cores of the two primary SOV cells. These SOV cells, along with vortices generated by flow over a submerged block of sediment at one of the banks, strongly influence distributions of the streamwise velocity and turbulent kinetic energy downstream of the junction. Comparison of the eddy-resolving simulation results with predictions from the steady Spalart-Allmaras RANS model shows that the latter fails to predict important features of the measured distributions of streamwise velocity and turbulent kinetic energy because the RANS model underpredicts the strength of the SOV cells. Analysis of instantaneous and mean bed shear stress distributions indicates that the SOV cells enhance bed shear stresses to a greater degree than the quasi-two-dimensional eddies in the mixing interface.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2010WR010018

Structure of turbulent flow at a river confluence with momentum and velocity ratios close to 1: Insight provided by an eddy‐resolving numerical simulation

Large Eddy Simulation (LES) using a dynamic Smagorinsky type subgridstress (SGS) model and Detached Eddy Simulation (DES) are applied toprediction and investigation of the flow around a sphere at a Reynoldsnumber of 104 in the subcritical regime. In this regime the boundarylayers at separation are laminar, and transition to turbulence occursfarther downstream in the separated shear layers via Kelvin–Helmholtz(K–H) instabilities. The dynamic eddy viscosity model of Germano et al.(Physics of Fluids3 (1991) 1760–1765) is used in the LES, while the current implementation of the DESemploys a formulation based on the Spalart–Allmaras (S–A) model. DES isa hybrid approach in which the closure is a modification to theproduction/destruction term of the original Reynolds-AveragedNavier–Stokes (RANS) model, reducing to RANS in the attached regions,and to LES away from the wall. In the present work where we simulate theflow over a sphere in the subcritical regime in which the boundarylayers at separation are laminar, DES can be viewed as LES with adifferent SGS model. Effects of the discretization scheme used toapproximate the convective terms are considered, along with sensitivityof predictions to changes in the additional model coefficient, CDES, in the DES formulation. DES and LES yield similar predictions of the wakestructure, large-scale vortex shedding and the Strouhal numberassociated with the low frequency mode in the wake. Predictions ofquantities such as the drag coefficient, wake frequencies, position oflaminar separation on the sphere, and the mean pressure andskin-friction distributions along the sphere are in good agreement withthe measurements of Achenbach (Journal of Fluid Mechanics54 (1972) 565–575). Predictions of the primaryReynolds shear stress, turbulent kinetic energy, eddy viscosity, andturbulent dissipation for the two models are also similar. In addition,both models successfully resolve the formation of the vortex tubes inthe detached shear layers along with the value of the Strouhal numberassociated with the high frequency instability mode, provided that thelevel of numerical dissipation introduced by the discretization schemeremains sufficiently low. Flow physics investigations are focused onunderstanding the wake structure in the subcritical regime.

LES and DES Investigations of Turbulent Flow over a Sphere at Re = 10,000

Publisher Summary
This chapter discusses the development of a numerical algorithm and solver capable of performing large-eddy simulation (LES) in geometries as complex as the combustor of a gas-turbine engine. The algorithm developed for unstructured grids is nondissipative yet robust at high Reynolds numbers on highly skewed grids. Results from validation in simple geometries are shown in the chapter along with the simulation results in the exceedingly complex geometry of a Pratt & Whitney gas-turbine combustor.

George Constantinescu

Papers

A numerical method for large-eddy simulation in complex geometries

Numerical investigations of flow over a sphere in the subcritical and supercritical regimes

Structure of turbulent flow at a river confluence with momentum and velocity ratios close to 1: Insight provided by an eddy‐resolving numerical simulation

LES and DES Investigations of Turbulent Flow over a Sphere at Re = 10,000

A numerical method for large-eddy simulation in complex geometries