D. I. Pullin

A low numerical dissipation patch-based adaptive mesh refinement method for large-eddy simulation of compressible flows

We present a methodology for the large-eddy simulation of compressible flows with a low-numerical dissipation scheme and structured adaptive mesh refinement (SAMR). A description of a conservative, flux-based hybrid numerical method that uses both centered finite-difference and a weighted essentially non-oscillatory (WENO) scheme is given, encompassing the cases of scheme alternation and internal mesh interfaces resulting from SAMR. In this method, the centered scheme is used in turbulent flow regions while WENO is employed to capture shocks. One-, two- and three-dimensional numerical experiments and example simulations are presented including homogeneous shock-free turbulence, a turbulent jet and the strongly shock-driven mixing of a Richtmyer-Meshkov instability.

Subgrid-scale modeling for large-eddy simulations of compressible turbulence

We present two phenomenological subgrid-scale (SGS) models for large-eddy simulations (LES) of compressible turbulent flows. A nonlinear model and a stretched-vortex model are tested in LES of compressible decaying isotropic turbulence. Results of LES at 323, 483, and 643 resolution are compared to corresponding 2563 direct numerical simulations (DNS) at a turbulent Mach number, Mt∼0.4. We use numerical schemes based on compact finite differences and study the effects of their order of accuracy on LES results. Both models give satisfactory agreement with DNS for the decay of the total turbulent kinetic energy. The probability densities (pdf) of energy transfer to subgrid scales obtained from filtered DNS and the SGS models are compared. Both models produce a narrower distribution of energy transfer than corresponding filtered DNS data, with less backscatter. The pdf of the alignment of components of the subgrid stress tensor and the eigenvectors of the rate-of-strain tensor obtained from the models reprod...

Unsteady forces on an accelerating plate and application to hovering insect flight

The aerodynamic forces on a flat plate accelerating from rest at fixed incidence in two-dimensional power-law flow are studied analytically and numerically. An inviscid approximation is made in which separation at the two plate edges is modelled by growing spiral vortex sheets, whose evolution is determined by the Birkhoff–Rott equation. A solution based on a similarity expansion is developed, valid when the scale of the separated vortex is much smaller than the plate dimension. The leading order is given by the well-known similarity growth of a vortex sheet from a semi-infinite flat plate, while equations at the second order describe the asymmetric sweeping effect of that component of the free-stream parallel to the plate. Owing to subtle cancellation, the unsteady vortex force exerted on the plate during the starting motion is independent of the sweeping effect and is determined by the similarity solution, to the order calculated. This gives a mechanism for dynamic stall based on a combination of unsteady vortex lift and pure added mass; the incidence angle for maximum vortex lift is

An adaptive high-order hybrid scheme for compressive, viscous flows with detailed chemistry

\arccos \sqrt{3/8}\,{\approx}\,52.2^\circ

Large-eddy simulation and wall modelling of turbulent channel flow

independent of the acceleration profile. Circulation on the flat plate makes no direct contribution. Both lift and drag force predictions from the unsteady inviscid theory are compared with those obtained from numerical solutions of the two-dimensional unsteady Navier–Stokes equations for an ellipse of high aspect ratio, and with predictions of Wagners classical theory. There is good agreement with numerical results at high incidence and moderate Reynolds number. The force per unit span predicted by the vortex theory is evaluated for parameters typical of insect wings and is found to be in reasonable agreement with numerical simulations. Estimates for the shed circulation and the size of the start-up vortices are also obtained. The significance of this flow as a mechanism for insect hovering flight is discussed.

On the non-local geometry of turbulence

A hybrid weighted essentially non-oscillatory (WENO)/centered-difference numerical method, with low numerical dissipation, high-order shock-capturing, and structured adaptive mesh refinement (SAMR), has been developed for the direct numerical simulation of the multicomponent, compressible, reactive Navier–Stokes equations. The method enables accurate resolution of diffusive processes within reaction zones. The approach combines time-split reactive source terms with a high-order, shock-capturing scheme specifically designed for diffusive flows. A description of the order-optimized, symmetric, finite difference, flux-based, hybrid WENO/centered-difference scheme is given, along with its implementation in a high-order SAMR framework. The implementation of new techniques for discontinuity flagging, scheme-switching, and high-order prolongation and restriction is described. In particular, the refined methodology does not require upwinded WENO at grid refinement interfaces for stability, allowing high-order prolongation and thereby eliminating a significant source of numerical diffusion within the overall code performance. A series of one-and two-dimensional test problems is used to verify the implementation, specifically the high-order accuracy of the diffusion terms. One-dimensional benchmarks include a viscous shock wave and a laminar flame. In two-space dimensions, a Lamb–Oseen vortex and an unstable diffusive detonation are considered, for which quantitative convergence is demonstrated. Further, a two-dimensional high-resolution simulation of a reactive Mach reflection phenomenon with diffusive multi-species mixing is presented.

On imploding cylindrical and spherical shock waves in a perfect gas

We report large-eddy simulation (LES) of turbulent channel flow. This LES neither resolves nor partially resolves the near-wall region. Instead, we develop a special near-wall subgrid-scale (SGS) model based on wall-parallel filtering and wall-normal averaging of the streamwise momentum equation, with an assumption of local inner scaling used to reduce the unsteady term. This gives an ordinary differential equation (ODE) for the wall shear stress at every wall location that is coupled with the LES. An extended form of the stretched-vortex SGS model, which incorporates the production of near-wall Reynolds shear stress due to the winding of streamwise momentum by near-wall attached SGS vortices, then provides a log relation for the streamwise velocity at the top boundary of the near-wall averaged domain. This allows calculation of an instantaneous slip velocity that is then used as a ‘virtual-wall’ boundary condition for the LES. A Karman-like constant is calculated dynamically as part of the LES. With this closure we perform LES of turbulent channel flow for Reynolds numbers Re_τ based on the friction velocity u_τ and the channel half-width δ in the range 2 × 10^3 to 2 × 10^7. Results, including SGS-extended longitudinal spectra, compare favourably with the direct numerical simulation (DNS) data of Hoyas & Jimenez (2006) at Re_τ = 2003 and maintain an O(1) grid dependence on Re_τ.

Direct numerical simulation and large-eddy simulation of stationary buoyancy-driven turbulence

A multi-scale methodology for the study of the non-local geometry of eddy structures in turbulence is developed. Starting from a given three-dimensional field, this consists of three main steps: extraction, characterization and classification of structures. The extraction step is done in two stages. First, a multi-scale decomposition based on the curvelet transform is applied to the full three-dimensional field, resulting in a finite set of component three-dimensional fields, one per scale. Second, by iso-contouring each component field at one or more iso-contour levels, a set of closed iso-surfaces is obtained that represents the structures at that scale. The characterization stage is based on the joint probability density function (p.d.f.), in terms of area coverage on each individual iso-surface, of two differential-geometry properties, the shape index and curvedness, plus the stretching parameter, a dimensionless global invariant of the surface. Taken together, this defines the geometrical signature of the iso-surface. The classification step is based on the construction of a finite set of parameters, obtained from algebraic functions of moments of the joint p.d.f. of each structure, that specify its location as a point in a multi-dimensional ‘feature space’. At each scale the set of points in feature space represents all structures at that scale, for the specified iso-contour value. This then allows the application, to the set, of clustering techniques that search for groups of structures with a common geometry. Results are presented of a first application of this technique to a passive scalar field obtained from 512 3 direct numerical simulation of scalar mixing by forced, isotropic turbulence ( Re λ = 265). These show transition, with decreasing scale, from blob-like structures in the larger scales to blob- and tube-like structures with small or moderate stretching in the inertial range of scales, and then toward tube and, predominantly, sheet-like structures with high level of stretching in the dissipation range of scales. Implications of these results for the dynamical behaviour of passive scalar stirring and mixing by turbulence are discussed.

An Eulerian hybrid WENO centered-difference solver for elastic-plastic solids

The problem of a cylindrically or spherically imploding and reflecting shock wave in a flow initially at rest is studied without the use of the strong-shock approximation. Dimensional arguments are first used to show that this flow admits a general solution where an infinitesimally weak shock from infinity strengthens as it converges towards the origin. For a perfect-gas equation of state, this solution depends only on the dimensionality of the flow and on the ratio of specific heats. The Guderley power-law result can then be interpreted as the leading-order, strong-shock approximation, valid near the origin at the implosion centre. We improve the Guderley solution by adding two further terms in the series expansion solution for both the incoming and the reflected shock waves. A series expansion, valid where the shock is still weak and very far from the origin, is also constructed. With an appropriate change of variables and using the exact shock-jump conditions, a numerical, characteristics-based solution is obtained describing the general shock motion from almost infinity to very close to the reflection point. Comparisons are made between the series expansions, the characteristics solution, and the results obtained using an Euler solver. These show that the addition of two terms to the Guderley solution significantly extends the range of validity of the strong-shock series expansion.

Eulerian adaptive finite-difference method for high-velocity impact and penetration problems

We report direct numerical simulation (DNS) and large-eddy simulation (LES) of statistically stationary buoyancy-driven turbulent mixing of an active scalar. We use an adaptation of the fringe-region technique, which continually supplies the flow with unmixed fluids at two opposite faces of a triply periodic domain in the presence of gravity, effectively maintaining an unstably stratified, but statistically stationary flow. We also develop a new method to solve the governing equations, based on the Helmholtz–Hodge decomposition, that guarantees discrete mass conservation regardless of iteration errors. Whilst some statistics were found to be sensitive to the computational box size, we show, from inner-scaled planar spectra, that the small scales exhibit similarity independent of Reynolds number, density ratio and aspect ratio. We also perform LES of the present flow using the stretched-vortex subgridscale (SGS) model. The utility of an SGS scalar flux closure for passive scalars is demonstrated in the present active-scalar, stably stratified flow setting. The multi-scale character of the stretched-vortex SGS model is shown to enable extension of some second-order statistics to subgrid scales. Comparisons with DNS velocity spectra and velocity-density cospectra show that both the resolved-scale and SGS-extended components of the LES spectra accurately capture important features of the DNS spectra, including small-scale anisotropy and the shape of the viscous roll-off.