A new code for astrophysical magnetohydrodynamics (MHD) is described. The code has been designed to be easily extensible for use with static and adaptive mesh refinement. It combines higher order Godunov methods with the constrained transport (CT) technique to enforce the divergence-free constraint on the magnetic field. Discretization is based on cell-centered volume averages for mass, momentum, and energy, and face-centered area averages for the magnetic field. Novel features of the algorithm include (1) a consistent framework for computing the time- and edge-averaged electric fields used by CT to evolve the magnetic field from the time- and area-averaged Godunov fluxes, (2) the extension to MHD of spatial reconstruction schemes that involve a dimensionally split time advance, and (3) the extension to MHD of two different dimensionally unsplit integration methods. Implementation of the algorithm in both C and FORTRAN95 is detailed, including strategies for parallelization using domain decomposition. Results from a test suite which includes problems in one-, two-, and three-dimensions for both hydrodynamics and MHD are given, not only to demonstrate the fidelity of the algorithms, but also to enable comparisons to other methods. The source code is freely available for download on the web.

/pdf/athena-a-new-code-for-astrophysical-mhd-1de7wtkkl8.pdf

Athena: a new code for astrophysical mhd

The Structure of Turbulent Shear Flow By Dr. A. A. Townsend. Pp. xii + 315. 8¾ in. × 5½ in. (Cambridge: At the University Press.) 40s.

/pdf/the-structure-of-turbulent-shear-flow-2yf3szhaqb.pdf

The Structure of Turbulent Shear Flow

Fluid dynamic turbulence is one of the most challenging computational physics problems because of the extremely wide range of time and space scales involved, the strong nonlinearity of the governing equations, and the many practical and important applications. While most linear fluid instabilities are well understood, the nonlinear interactions among them makes even the relatively simple limit of homogeneous isotropic turbulence difficult to treat physically, mathematically, and computationally. Turbulence is modeled computationally by a two-stage bootstrap process. The first stage, direct numerical simulation, attempts to resolve the relevant physical time and space scales but its application is limited to diffusive flows with a relatively small Reynolds number (Re). Using direct numerical simulation to provide a database, in turn, allows calibration of phenomenological turbulence models for engineering applications. Large eddy simulation incorporates a form of turbulence modeling applicable when the large-scale flows of interest are intrinsically time dependent, thus throwing common statistical models into question. A promising approach to large eddy simulation involves the use of high-resolution monotone computational fluid dynamics algorithms such as flux-corrected transport or the piecewise parabolic method which have intrinsic subgrid turbulence models coupled naturally to the resolved scales in the computed flow. The physical considerations underlying and evidence supporting this monotone integrated large eddy simulation approach are discussed.

New insights into large eddy simulation

After several years of planning, the 1st International Workshop on High-Order CFD Methods was successfully held in Nashville, Tennessee, on January 7-8, 2012, just before the 50th Aerospace Sciences Meeting. The American Institute of Aeronautics and Astronautics, the Air Force Office of Scientific Research, and the German Aerospace Center provided much needed support, financial and moral. Over 70 participants from all over the world across the research spectrum of academia, government labs, and private industry attended the workshop. Many exciting results were presented. In this review article, the main motivation and major findings from the workshop are described. Pacing items requiring further effort are presented. © 2013 John Wiley & Sons, Ltd.

/pdf/high-order-cfd-methods-current-status-and-perspective-1hm5p7dcff.pdf

High-order CFD methods: Current status and perspective

We review wall-bounded turbulent flows, particularly high–Reynolds number, zero–pressure gradient boundary layers, and fully developed pipe and channel flows. It is apparent that the approach to an asymptotically high–Reynolds number state is slow, but at a sufficiently high Reynolds number the log law remains a fundamental part of the mean flow description. With regard to the coherent motions, very-large-scale motions or superstructures exist at all Reynolds numbers, but they become increasingly important with Reynolds number in terms of their energy content and their interaction with the smaller scales near the wall. There is accumulating evidence that certain features are flow specific, such as the constants in the log law and the behavior of the very large scales and their interaction with the large scales (consisting of vortex packets). Moreover, the refined attached-eddy hypothesis continues to provide an important theoretical framework for the structure of wall-bounded turbulent flows.

/pdf/high-reynolds-number-wall-turbulence-1s7ld2qfbi.pdf

High–Reynolds Number Wall Turbulence

We present results from 1283 and 2563 direct numerical simulations (DNS) of decaying compressible, isotropic turbulence at fluctuation Mach numbers of Mt∼0.1–0.5 and at Taylor Reynolds numbers Reλ=O(50–100). The presence or absence of fluctuations of thermodynamic quantities as well as velocity divergence in the initial conditions are found to have a negligible effect on the decay of turbulent kinetic energy. The decay of the turbulent kinetic energy shows no significant effect of Mt and power laws fitted to the timewise decay exhibit exponents n=1.3–1.7 that are similar to those found for decaying incompressible turbulence. The main new phenomenon produced by compressibility is the appearance of random shocklets which form during the main part of the decay. An algorithm is developed to extract and quantify the shocklet statistics from the DNS fields. A model for the probability density function (PDF) of the shocklet strength Mn−1 (Mn is the normal shock Mach number) is derived based on combining weak-sho...

Direct numerical simulation of decaying compressible turbulence and shocklet statistics

Large-eddy simulations of the Richtmyer–Meshkov instability with reshock are pre- sented and the results are compared with experiments. Several configurations of shocks initially travelling from light (air) to heavy (sulfur hexafluoride, SF6) have been simulated to match previous experiments and good agreement is found in the growth rates of the turbulent mixing zone (TMZ). The stretched-vortex subgrid model used in this study allows for subgrid continuation modelling, where statistics of the unresolved scales of the flow are estimated. In particular, this multiscale modelling allows the anisotropy of the flow to be extended to the dissipation scale, eta, and estimates to be formed for the subgrid probability density function of the mixture fraction of air/SF6 based on the subgrid variance, including the effect of Schmidt number.

/pdf/large-eddy-simulation-and-multiscale-modelling-of-a-3wtuw7mj0u.pdf

Large-eddy simulation and multiscale modelling of a Richtmyer–Meshkov instability with reshock

https://scholar.harvard.edu/files/wimvanrees/files/rees2011.pdf

A comparison of vortex and pseudo-spectral methods for the simulation of periodic vortical flows at high Reynolds numbers

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_τ.

/pdf/large-eddy-simulation-and-wall-modelling-of-turbulent-41eow14qtf.pdf

Dale Pullin

Papers

Direct numerical simulation of decaying compressible turbulence and shocklet statistics

Large-eddy simulation and multiscale modelling of a Richtmyer–Meshkov instability with reshock

A comparison of vortex and pseudo-spectral methods for the simulation of periodic vortical flows at high Reynolds numbers

Large-eddy simulation and wall modelling of turbulent channel flow

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