Direct Numerical Simulation of turbulent Taylor-Couette flow

The direct numerical simulation (DNS) of the Taylor--Couette flow in the fully turbulent regime is described. The numerical method extends the work by Quadrio & Luchini (Eur. J. Mech. B / Fluids, v.21, pp.413--427, 2002), and is based on a parallel computer code which uses mixed spatial discretization (spectral schemes in the homogeneous directions, and fourth-order, compact explicit finite-difference schemes in the radial direction). A DNS is carried out to simulate for the first time the turbulent Taylor--Couette flow in the turbulent regime. Statistical quantities are computed to complement the existing experimental information, with a view to compare it to planar, pressure-driven turbulent flow at the same value of the Reynolds number. The main source for differences in flow statistics between plane and curved-wall flows is attributed to the presence of large-scale rotating structures generated by curvature effects.Comment: To appear in European Journal of Mechanics B / Fluid

Embedded direct numerical simulation for aeronautical CFD

AbstractA method is proposed by which a direct numerical simulation of the compressible Navier-Stokes equations may be embedded within a more general aeronautical CFD code. The method may be applied to any code which solves the Euler equations or the Favre-averaged Navier-Stokes equations. A formal decomposition of the flowfield is used to derive modified equations for use with direct numerical simulation solvers. Some preliminary applications for model flows with transitional separation bubbles are given.</jats:p

Direct numerical simulation of multiphase flows with unstable interfaces

Published under licence in Journal of Physics: Conference Series by IOP Publishing Ltd. Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.This paper presents a numerical model that intends to simulate efficiently the surface instability that arise in multiphase flows, typically liquid-gas, both for laminar or turbulent regimes. The model is developed on the in-house computing platform TermoFluids , and operates the finite-volume, direct numerical simulation (DNS) of multiphase flows by means of a conservative level-set method for the interface-capturing. The mesh size is optimized by means of an adaptive mesh refinement (AMR) strategy, that allows the dynamic re-concentration of the mesh in the vicinity of the interfaces between fluids, in order to correctly represent the diverse structures (as ligaments and droplets) that may rise from unstable phenomena. In addition, special attention is given to the discretization of the various terms of the momentum equations, to ensure stability of the flow and correct representation of turbulent vortices. As shown, the method is capable of truthfully simulate the interface phenomena as the Kelvin-Helmholtz instability and the Plateau-Rayleigh instability, both in the case of 2-D and 3-D configurations. Therefore it is suitable for the simulation of complex phenomena such as simulation of air-blast atomization, with several important application in the field of automotive and aerospace engines. A prove is given by our preliminary study of the 3-D coaxial liquid-gas jet.Peer ReviewedPostprint (published version

Direct numerical simulation of dispersed particles in a compressible fluid

We present a direct numerical simulation method for investigating the dynamics of dispersed particles in a compressible solvent fluid. The validity of the simulation is examined by calculating the velocity relaxation of an impulsively forced spherical particle with a known analytical solution. The simulation also gives information about the fluid motion, which provides some insight into the particle motion. Fluctuations are also introduced by random stress, and the validity of this case is examined by comparing the calculation results with the fluctuation-dissipation theorem.Comment: 7 pages, 5 figure

Proper Orthogonal Decomposition Closure Models For Turbulent Flows: A Numerical Comparison

This paper puts forth two new closure models for the proper orthogonal decomposition reduced-order modeling of structurally dominated turbulent flows: the dynamic subgrid-scale model and the variational multiscale model. These models, which are considered state-of-the-art in large eddy simulation, together with the mixing length and the Smagorinsky closure models, are tested in the numerical simulation of a 3D turbulent flow around a circular cylinder at Re = 1,000. Two criteria are used in judging the performance of the proper orthogonal decomposition reduced-order models: the kinetic energy spectrum and the time evolution of the POD coefficients. All the numerical results are benchmarked against a direct numerical simulation. Based on these numerical results, we conclude that the dynamic subgrid-scale and the variational multiscale models perform best.Comment: 28 pages, 6 figure

Direct Numerical Simulation Tests of Eddy Viscosity in Two Dimensions

Two-parametric eddy viscosity (TPEV) and other spectral characteristics of two-dimensional (2D) turbulence in the energy transfer sub-range are calculated from direct numerical simulation (DNS) with 512$^2$ resolution. The DNS-based TPEV is compared with those calculated from the test field model (TFM) and from the renormalization group (RG) theory. Very good agreement between all three results is observed.Comment: 9 pages (RevTeX) and 5 figures, published in Phys. Fluids 6, 2548 (1994

A coupled approximate deconvolution and dynamic mixed scale model for large-eddy simulation

Large-eddy simulations of incompressible Newtonian fluid flows with approximate deconvolution models based on the van Cittert method are reported. The Legendre spectral element method is used for the spatial discretization to solve the filtered Navier--Stokes equations. A novel variant of approximate deconvolution models blended with a mixed scale model using a dynamic evaluation of the subgrid-viscosity constant is proposed. This model is validated by comparing the large-eddy simulation with the direct numerical simulation of the flow in a lid-driven cubical cavity, performed at a Reynolds number of 12'000. Subgrid modeling in the case of a flow with coexisting laminar, transitional and turbulent zones such as the lid-driven cubical cavity flow represents a challenging problem. Moreover, the coupling with the spectral element method having very low numerical dissipation and dispersion builds a well suited framework to analyze the efficiency of a subgrid model. First- and second-order statistics obtained using this new model are showing very good agreement with the direct numerical simulation. Filtering operations rely on an invertible filter applied in a modal basis and preserving the C0-continuity across elements. No clipping on dynamic parameters was needed to preserve numerical stability

Dissipative Pulsar Magnetosphere

Dissipative axisymmetric pulsar magnetosphere is calculated by a direct numerical simulation of the Strong-Field Electrodynamics equations. The magnetic separatrix disappears, it is replaced by a region of enhanced dissipation. With a better numerical scheme, one should be able to calculate the bolometric lightcurves for a given conductivity.Comment: 2 pages, 10 figures, minor changes for the journa