2,044 research outputs found
Comparison of Subgrid-scale Viscosity Models and Selective Filtering Strategy for Large-eddy Simulations
Explicitly filtered large-eddy simulations (LES), combining high-accuracy schemes with the use of a selective filtering without adding an explicit subgrid-scales (SGS) model, are carried out for the Taylor-Green-vortex and the supersonic-boundary-layer cases. First, the present approach is validated against direct numerical simulation (DNS) results. Subsequently, several SGS models are implemented in order to investigate if they can improve the initial filter-based methodology. It is shown that the most accurate results are obtained when the filtering is used alone as an implicit model, and for a minimal cost. Moreover, the tests for the Taylor-Green vortex indicate that the discretization error from the numerical methods, notably the dissipation error from the high-order filtering, can have a greater influence than the SGS models
An Asymptotic Preserving Maxwell Solver Resulting in the Darwin Limit of Electrodynamics
In plasma simulations, where the speed of light divided by a characteristic
length is at a much higher frequency than other relevant parameters in the
underlying system, such as the plasma frequency, implicit methods begin to play
an important role in generating efficient solutions in these multi-scale
problems. Under conditions of scale separation, one can rescale Maxwell's
equations in such a way as to give a magneto static limit known as the Darwin
approximation of electromagnetics. In this work, we present a new approach to
solve Maxwell's equations based on a Method of Lines Transpose (MOL)
formulation, combined with a fast summation method with computational
complexity , where is the number of grid points (particles).
Under appropriate scaling, we show that the proposed schemes result in
asymptotic preserving methods that can recover the Darwin limit of
electrodynamics
Numerics and subgrid-scale modeling in large eddy simulations of stratocumulus clouds
Stratocumulus clouds are the most common type of boundary layer cloud; their radiative effects strongly modulate climate. Large eddy simulations (LES) of stratocumulus clouds often struggle to maintain fidelity to observations because of the sharp gradients occurring at the entrainment interfacial layer at the cloud top. The challenge posed to LES by stratocumulus clouds is evident in the wide range of solutions found in the LES intercomparison based on the DYCOMS-II field campaign, where simulated liquid water paths for identical initial and boundary conditions varied by a factor of nearly 12. Here we revisit the DYCOMS-II RF01 case and show that the wide range of previous LES results can be realized in a single LES code by varying only the numerical treatment of the equations of motion and the nature of subgrid-scale (SGS) closures. The simulations that maintain the greatest fidelity to DYCOMS-II observations are identified. The results show that using weighted essentially non-oscillatory (WENO) numerics for all resolved advective terms and no explicit SGS closure consistently produces the highest-fidelity simulations. This suggests that the numerical dissipation inherent in WENO schemes functions as a high-quality, implicit SGS closure for this stratocumulus case. Conversely, using oscillatory centered difference numerical schemes for momentum advection, WENO numerics for scalars, and explicitly modeled SGS fluxes consistently produces the lowest-fidelity simulations. We attribute this to the production of anomalously large SGS fluxes near the cloud tops through the interaction of numerical error in the momentum field with the scalar SGS model
Hermite regularization of the Lattice Boltzmann Method for open source computational aeroacoustics
The lattice Boltzmann method (LBM) is emerging as a powerful engineering tool
for aeroacoustic computations. However, the LBM has been shown to present
accuracy and stability issues in the medium-low Mach number range, that is of
interest for aeroacoustic applications. Several solutions have been proposed
but often are too computationally expensive, do not retain the simplicity and
the advantages typical of the LBM, or are not described well enough to be
usable by the community due to proprietary software policies. We propose to use
an original regularized collision operator, based on the expansion in Hermite
polynomials, that greatly improves the accuracy and stability of the LBM
without altering significantly its algorithm. The regularized LBM can be easily
coupled with both non-reflective boundary conditions and a multi-level grid
strategy, essential ingredients for aeroacoustic simulations. Excellent
agreement was found between our approach and both experimental and numerical
data on two different benchmarks: the laminar, unsteady flow past a 2D cylinder
and the 3D turbulent jet. Finally, most of the aeroacoustic computations with
LBM have been done with commercial softwares, while here the entire theoretical
framework is implemented on top of an open source library (Palabos).Comment: 34 pages, 12 figures, The Journal of the Acoustical Society of
America (in press
A space-time discretization procedure for wave propagation problems
Higher order compact algorithms are developed for the numerical simulation of wave propagation by using the concept of a discrete dispersion relation. The dispersion relation is the imprint of any linear operator in space-time. The discrete dispersion relation is derived from the continuous dispersion relation by examining the process by which locally plane waves propagate through a chosen grid. The exponential structure of the discrete dispersion relation suggests an efficient splitting of convective and diffusive terms for dissipative waves. Fourth- and eighth-order convection schemes are examined that involve only three or five spatial grid points. These algorithms are subject to the same restrictions that govern the use of dispersion relations in the constructions of asymptotic expansions to nonlinear evolution equations. A new eighth-order scheme is developed that is exact for Courant numbers of 1, 2, 3, and 4. Examples are given of a pulse and step wave with a small amount of physical diffusion
- …