Application of a Reynolds stress turbulence model to the compressible shear layer

Theoretically based turbulence models have had success in predicting many features of incompressible, free shear layers. However, attempts to extend these models to the high-speed, compressible shear layer have been less effective. In the present work, the compressible shear layer was studied with a second-order turbulence closure, which initially used only variable density extensions of incompressible models for the Reynolds stress transport equation and the dissipation rate transport equation. The quasi-incompressible closure was unsuccessful; the predicted effect of the convective Mach number on the shear layer growth rate was significantly smaller than that observed in experiments. Having thus confirmed that compressibility effects have to be explicitly considered, a new model for the compressible dissipation was introduced into the closure. This model is based on a low Mach number, asymptotic analysis of the Navier-Stokes equations, and on direct numerical simulation of compressible, isotropic turbulence. The use of the new model for the compressible dissipation led to good agreement of the computed growth rates with the experimental data. Both the computations and the experiments indicate a dramatic reduction in the growth rate when the convective Mach number is increased. Experimental data on the normalized maximum turbulence intensities and shear stress also show a reduction with increasing Mach number

Numerical study of three-dimensional separation and flow control at a wing/body junction

The problem of three-dimensional separation and flow control at a wing/body junction has been investigated numerically using a three-dimensional Navier-Stokes code. The numerical code employs an algebraic grid generation technique for generating the grid for unmodified junction and an elliptic grid generation technique for filleted fin junction. The results for laminar flow past a blunt fin/flat plate junction demonstrate that after grid refinement, the computations agree with experiment and reveal a strong dependency of the number of vortices at the junction on Mach number and Reynolds number. The numerical results for pressure distribution, particle paths and limiting streamlines for turbulent flow past a swept fin show a decrease in the peak pressure and in the extent of the separated flow region compared to the laminar case. The results for a filleted juncture indicate that the streamline patterns lose much of their vortical character with proper filleting. Fillets with a radius of three and one-half times the fin leading edge diameter or two times the incoming boundary layer thickness, significantly weaken the usual necklace interaction vortex for the Mach number and Reynolds number considered in the present study

Fluctuating hydrodynamics of multi-species, non-reactive mixtures

In this paper we discuss the formulation of the fuctuating Navier-Stokes (FNS) equations for multi-species, non-reactive fluids. In particular, we establish a form suitable for numerical solution of the resulting stochastic partial differential equations. An accurate and efficient numerical scheme, based on our previous methods for single species and binary mixtures, is presented and tested at equilibrium as well as for a variety of non-equilibrium problems. These include the study of giant nonequilibrium concentration fluctuations in a ternary mixture in the presence of a diffusion barrier, the triggering of a Rayleigh-Taylor instability by diffusion in a four-species mixture, as well as reverse diffusion in a ternary mixture. Good agreement with theory and experiment demonstrates that the formulation is robust and can serve as a useful tool in the study of thermal fluctuations for multi-species fluids. The extension to include chemical reactions will be treated in a sequel paper

Particle-Hole Asymmetry and Brightening of Solitons in A Strongly Repulsive BEC

We study solitary wave propagation in the condensate of a system of hard-core bosons with nearest-neighbor interactions. For this strongly repulsive system, the evolution equation for the condensate order parameter of the system, obtained using spin coherent state averages is different from the usual Gross-Pitaevskii equation (GPE). The system is found to support two kinds of solitons when there is a particle-hole imbalance: a dark soliton that dies out as the velocity approaches the sound velocity, and a new type of soliton which brightens and persists all the way up to the sound velocity, transforming into a periodic wave train at supersonic speed. Analogous to the GPE soliton, the energy-momentum dispersion for both solitons is characterized by Lieb II modes.Comment: Accepted for publication in PRL, Nov 12, 200

Velocity Correlations, Diffusion and Stochasticity in a One-Dimensional System

We consider the motion of a test particle in a one-dimensional system of equal-mass point particles. The test particle plays the role of a microscopic "piston" that separates two hard-point gases with different concentrations and arbitrary initial velocity distributions. In the homogeneous case when the gases on either side of the piston are in the same macroscopic state, we compute and analyze the stationary velocity autocorrelation function C(t). Explicit expressions are obtained for certain typical velocity distributions, serving to elucidate in particular the asymptotic behavior of C(t). It is shown that the occurrence of a non-vanishing probability mass at zero velocity is necessary for the occurrence of a long-time tail in C(t). The conditions under which this is a $t^{-3}$ tail are determined. Turning to the inhomogeneous system with different macroscopic states on either side of the piston, we determine its effective diffusion coefficient from the asymptotic behavior of the variance of its position, as well as the leading behavior of the other moments about the mean. Finally, we present an interpretation of the effective noise arising from the dynamics of the two gases, and thence that of the stochastic process to which the position of any particle in the system reduces in the thermodynamic limit.Comment: 22 files, 2 eps figures. Submitted to PR

Impact of the various spin and orbital ordering processes on multiferroic properties of orthovanadate DyVO3

The orthovanadate DyVO3 crystal, known to exhibit multiple structural, spin and orbital ordering transitions, is presently investigated on the basis of magnetization, heat capacity, resistivity, dielectric and polarization measurements. Our main result is experimental evidence for the existence of multiferroicity below a high TC of 108 K over a wide temperature range including different spin-orbital ordered states. The onset of ferroelectricity is found to coincide with the antiferromagnetic C-type spin ordering transition taking place at 108 K, which indicates that DyVO3 belongs to type II multiferroics exhibiting a coupling between magnetism and ferroelectricity. Some anomalies detected on the temperature dependence of electric polarization are discussed with respect to the nature of the spin-orbital ordered states of the V sublattice and the degree of spin alignment in the Dy sublattice. The orthovanadates RVO3 (R = rare earth or Y) form an important new category for searching for high-TC multiferroics.Comment: 25 pages, 7 figures, 68 references, one supplementary material, Physical Review B, Published 23 July 201

Hybrid High Performance Computing (HPC) + Cloud for Scientific Computing

The HPC+Cloud framework has been built to enable on-premise HPC jobs to use resources from cloud computing nodes. As part of designing the software framework, public cloud providers, namely Amazon AWS, Microsoft Azure and NeCTAR were benchmarked against one another, and Microsoft Azure was determined to be the most suitable cloud component in the proposed HPC+Cloud software framework. Finally, an HPC+Cloud cluster was built using the HPC+Cloud software framework and then was validated by conducting HPC processing benchmarks