Direct Numerical Simulation of Passive Scalar in Decaying Compressible Turbulence

The passive scalars in the decaying compressible turbulence with the initial Reynolds number (defined by Taylor scale and RMS velocity) Re=72, the initial turbulent Mach numbers (defined by RMS velocity and mean sound speed) Mt=0.2-0.9, and the Schmidt numbers of passive scalar Sc=2-10 are numerically simulated by using a 7th order upwind difference scheme and 8th order group velocity control scheme. The computed results are validated with different numerical methods and different mesh sizes. The Batchelor scaling with k(-1) range is found in scalar spectra. The passive scalar spectra decay faster with the increasing turbulent Mach number. The extended self-similarity (ESS) is found in the passive scalar of compressible turbulence.

Turbulent passive scalar field of a small prandtl number

The statistical-mechanics theory of the passive scalar field convected by turbulence, developed in an earlier paper [Phys. Fluids 28, 1299 (1985)], is extended to the case of a small molecular Prandtl number. The set of governing integral equations is solved by the equation-error method. The resultant scalar-variance spectrum for the inertial range is F(k)~x−5/3/[1+1.21x1.67(1+0.353x2.32)], where x is the wavenumber scaled by Corrsin's dissipation wavenumber. This result reduces to the − (5)/(3) law in the inertial-convective range. It also approximately reduces to the − (17)/(3) law in the inertial-diffusive range, but the proportionality constant differs from Batchelor's by a factor of 3.6.

A passive scalar field convected by turbulence

Classical statistical mechanics is applied to the study of a passive scalar field convected by isotropic turbulence. A complete set of independent real parameters and dynamic equations are worked out to describe the dynamic state of the passive scalar field. The corresponding Liouville equation is solved by a perturbation method based upon a Langevin–Fokker–Planck model. The closure problem is treated by a variational approach reported in earlier papers. Two integral equations are obtained for two unknown functions: the scalar variance spectrum F(k) and the effective damping coefficient (k). The appearance of the energy spectrum of the velocity field in the two integral equations represents the coupling of the scalar field with the velocity field. As an application of the theory, the two integral equations are solved to derive the inertial-convective-range spectrum, obtaining F(k)=0.61 −1/3 k−5/3. Here is the dissipation rate of the scalar variance and is the dissipation rate of the energy of the velocity field. This theoretical value of the scalar Kolmogorov constant, 0.61, is in good agreement with experiments.

DIFFERENTIAL GEOMETRICAL METHODS IN THE STUDY OF OPTICAL-TRANSMISSION (SCALAR THEORY) .1. STATIC TRANSMISSION CASE

Static optical transmission is restudied by postulation of the optical path as the proper element in a three-dimensional Riemannian manifold (no torsion); this postulation can be applied to describe the light-medium interactive system. On the basis of the postulation, the behaviors of light transmitting through the medium with refractive index n are investigated, the investigation covering the realms of both geometrical optics and wave optics. The wave equation of light in static transmission is studied modally, the postulation being employed to derive the exact form of the optical field equation in a medium (in which the light is viewed as a single-component field). Correspondingly, the relationships concerning the conservation of optical fluid and the dynamic properties are given, and some simple applications of the theories mentioned are presented.

DIFFERENTIAL GEOMETRICAL METHODS IN THE STUDY OF OPTICAL-TRANSMISSION (SCALAR THEORY) .2. TIME-DEPENDENT TRANSMISSION THEORY

By generalization of the methods presented in Part I of the study [J. Opt. Soc. Am. A 12, 600 (1994)] to the four-dimensional (4D) Riemannian manifold case, the time-dependent behavior of light transmitting in a medium is investigated theoretically by the geodesic equation and curvature in a 4D manifold. In addition, the field equation is restudied, and the 4D conserved current of the optical fluid and its conservation equation are derived and applied to deduce the time-dependent general refractive index. On this basis the forces acting on the fluid are dynamically analyzed and the self-consistency analysis is given.

The effect of the scalar-isovector meson field on hyperon-rich neutron star matter

We investigate the effect of the calar-isovector delta-meson field on the equation of state (EOS) and composition of hyperonic neutron star matter, and the properties of hyperonic neutron stars within the frame work of the relativistic mean field theory. The influence of the delta-field turns out to be quite different and generally weaker for hyperonic neutron star matter as compared to that for npe mu neutron star matter. We find that inclusion of the delta-field enhances the strangeness content slightly and consequently moderately softens the EOS of neutron star matter in its hyperonic phase. As for the composition of hyperonic star matter, the effect of the delta-field is shown to shift the onset of the negatively-charged (positively-charged) hyperons to slightly lower (higher) densities and to enhance (reduce) their abundances. The influence of the delta-field on the maximum mass of hyperonic neutron stars is found to be fairly weak, where as inclusion of the delta-field turns out to enhance sizably both the radii and the moments of inertia of neutron stars with given masses. It is also shown that the effects of the delta-field on the properties of hyperonic neutron stars remain similar in the case of switching off the Sigma hyperons.

Electro-Osmotic Flow and Mixing in Heterogeneous Microchannels

Analytical and numerical studies of secondary electro-osmotic flow EOF and its mixing in microchannels with heterogeneous zeta potentials are carried out in the present work. The secondary EOFs are analyzed by solving the Stokes equation with heterogeneous slip velocity boundary conditions. The analytical results obtained are compared with the direct numerical simulation of the Navier-Stokes equations. The secondary EOFs could transport scalar in larger areas and increase the scalar gradients, which significantly improve the mixing rate of scalars. It is shown that the heterogeneous zeta potentials could generate complex flow patterns and be used to enhance scalar mixing.

Hierarchy of structure functions for passive scalars advected by turbulent flows

A hierarchical model is proposed for the joint moments of the passive scalar dissipation and the velocity dissipation in fluid turbulence. This model predicts that the joint probability density function (PDF) of the dissipations is a bivariate log-Poisson. An analytical calculation of the scaling exponents of structure functions of the passive scalar is carried out for this hierarchical model, showing a good agreement with the results of direct numerical simulations and experiments.

Two-Point Closure Strategy in the Mapping Closure Approximation Approach

A two-point closure strategy in mapping closure approximation (MCA) approach is developed for the evolution of the probability density function (PDF) of a scalar advected by stochastic velocity fields. The MCA approach is based on multipoint statistics. We formulate a MCA modeled system using the one-point PDFs and two-point correlations. The MCA models can describe both the evolution of the PDF shape and the rate at which the PDF evolves.

The cellular structure of a two-dimensional H2/O2/Ar detonation wave

In this paper, the cellular structure of a two-dimensional detonation wave in a low pressure H2/O2/Ar mixture calculated with a detailed chemical reaction model, high order scheme and high resolution grids is investigated. The regular cellular structure is produced about 1 ms after introducing perturbations in the reaction zone of a steady one-dimensional detonation wave. It is found from the present resolution study that the discrepancies concerning the structure type arising from the coarser grid employed can be resolved using a sufficiently fine grid size of 0.05 mm and below and shows a double-Mach-like strong-type configuration. During the structure evolution process, the structure configuration does not change much in the periods before and after the triple point collision. Through the triple point collision, three regular collision processes are observed and are followed by a quick change to the double-Mach-like configuration. The simulated structure tracks show that there are three different tracks associated with different triple points or the kink on the transverse wave. Comparisons with previous work and experiments indicate the presence of a strong structure for an ordinary detonation.