Numerical investigation of the underexpanded sonic coaxial jets is carried out using large eddy simulation for three typical inner nozzle lip-thicknesses. Various fundamental mechanisms dictating the flow phenomena including shock structure, shear layer evolution and sound production are investigated. It is found that the inner nozzle lip induces a recirculation zone between inner and outer jets, which significantly influences the behaviors of shock structures and shear layers. The sound properties of the coaxial jets are further analyzed in detail. As the inner lip-thickness increases, the helical screech mode switches to an axisymmetric one and high-frequency screech also occurs with an oscillation frequency of recirculation zone. Based on the temporal Fourier transform and correlation analysis, the primary sources of low- and high-frequency screeches are associated with the downstream shock cells in the jet column and the secondary shock structures in the outer annular jet, respectively. The proper orthogonal decomposition analysis reveals that the dominant structures constructed by the most energetic modes shift from the downstream shock cells region to the upstream secondary shock region as the lip-thickness increases. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to the coherent structures and sound properties in sonic coaxial jets.

Homogeneous equilibrium model (HEM) has been widely used in cavitating flow simulations. The major feature of this model is that a single equation of state (EOS) is proposed to describe the thermal behavior of bubbly liquid, where both kinematic and thermal equilibrium is assumed between two phases. In this paper, the HEM was coupled with multi-relaxation-time lattice Boltzmann model (MRT-LBM) and the acoustic behavior was simulated. Two approaches were applied alternatively: adjusting speed of sound (Buick, J. Phys. A, 2006, 39:13807-13815) and setting real gas EOS. Both approaches result in high accuracy in acoustic speed predictions for different void (gas) volume of fractions. It is demonstrated that LBM could be successfully applied as a Navier-Stokes equation solver for industrial applications. However, further dissipation and dispersion analysis shows that Shan-Chen type approaches of LBM are deficient, especially in large wave-number region.

In this paper, a conservative fifth-order upwind compact scheme using centered stencil is introduced. This scheme uses asymmetric coefficients to achieve the upwind property since the stencil is symmetric. Theoretical analysis shows that the proposed scheme is low-dissipative and has a relatively large stability range. To maintain the convergence rate of the whole spatial discretization, a proper non-periodic boundary scheme is also proposed. A detailed analysis shows that the spatial discretization implemented with the boundary scheme proposed by Pirozzoli [J. Comput. Phys., 178 (2001), pp. 81–117] is approximately fourth-order. Furthermore, a hybrid methodology, coupling the compact scheme with WENO scheme, is adopted for problems with discontinuities. Numerical results demonstrate the effectiveness of the proposed scheme.

As the numerical resolution is increased and the discretisation error decreases, the lattice Boltzmann method tends towards the discrete-velocity Boltzmann equation (DVBE). An expression for the propagation properties of plane sound waves is found for this equation. This expression is compared to similar ones from the Navier-Stokes and Burnett models, and is found to be closest to the latter. The anisotropy of sound propagation with the DVBE is examined using a two-dimensional velocity set. It is found that both the anisotropy and the deviation between the models is negligible if the Knudsen number is less than 1 by at least an order of magnitude.

A surface based lattice Boltzmann impedance boundary condition (BC) using Ozyoruk’s model [J. Comput. Phys., 146 (1998), pp. 29-57] is proposed and implemented in PowerFLOW. In Ozyoruk’s model, pressure fluctuation is directly linked to normal velocity on an impedance surface. In the present study, the relation between pressure and normal velocity is realized precisely by imposing a mass flux on the surface. This impedance BC is generalized and can handle complex geometry. Combined with the turbulence model in the lattice Boltzmann solver PowerFLOW, this BC can be used to model the effect of a liner in presence of a complex 3D turbulent flow. Preliminary simulations of the NASA Langley grazing flow tube and Kundt tube show satisfying agreement with experimental results.

When gas flows through corrugated pipes, pressure waves interacting with vortex shedding can produce distinct tonal noise and structural vibration. Based on established observations, a model is proposed which couples an acoustic pipe and self-excited oscillations with vortex shedding over the corrugation cavities. In the model, the acoustic response of the corrugated pipe is simulated by connecting the lossless medium moving with a constant velocity with a source based on a discrete distribution of van der Pol oscillators arranged along the pipe. Our time accurate solutions exhibit dynamic behavior consistent with that experimentally observed, including the lock-in frequency of vortex shedding, standing waves and the onset fluid velocity capable of generating the lock-in.

The paper offers a study of the inverse Laplace transforms of the functions ${{I}_{n}}\left( rs \right)\,{{\{sI_{n}^{\prime }\,\left( s \right)\}}^{-1}}$ where ${{I}_{n}}$ is the modified Bessel function of the first kind and $r$ is a parameter. The present study is a continuation of the author's previous work on the singular behavior of the special case of the functions in question, $r=1$. The general case of $r\,\in \,\left[ 0,\,1 \right]$ is addressed, and it is shown that the inverse Laplace transforms for such $r$ exhibit significantly more complex behavior than their predecessors, even though they still only have two different types of points of discontinuity: singularities and finite discontinuities. The functions studied originate from non-stationary fluid-structure interaction, and as such are of interest to researchers working in the area.