A Compressible Turbulence Model for Pressure—Strain

This work focuses on the performance and validation of compressible turbulence models for the pressure-strain correlation. Considering the Launder Reece and Rodi (LRR) incompressible model for the pressure-strain correlation, Adumitroaie et al., Huang et al., and Marzougui et al., used different modeling approaches to develop turbulence models, taking into account compressibility effects for this term. Two numerical coefficients are dependent on the turbulent Mach number, and all of the remaining coefficients conserve the same values as in the original LRR model. The models do not correctly predict the compressible turbulence at a high-speed shear flow. So, the revision of these models is the major aim of this study. In the present work, the compressible model for the pressure-strain correlation developed by Khlifi−Lili, involving the turbulent Mach number, the gradient, and the convective Mach numbers, is used to modify the linear mean shear strain and the slow terms of the previous models. The models are tested in two compressible turbulent flows: homogeneous shear flow and the newly developed plane mixing layers. The predicted results of the proposed modifications of the Adumitroaie et al., Huang et al., and Marzougui et al., models and of its universal versions are compared with direct numerical simulation (DNS) and experiment data. The results show that the important parameters of compressibility in homogeneous shear flow and in the mixing layers are well predicted by the proposal models.

Chemical production on a deforming substrate

The interplay between chemical reaction and substrate deformation is discussed by adapting Ranz's formulation for scalar mixing to the case of a reactive mixture between segregated reactants, initially separated by an interface whose thickness may not be vanishingly small. Experiments in a simple shear flow demonstrate the existence of three regimes depending on the Damköhler number $Da=t_s/t_c$ where $t_s$ is the mixing time of the interface width and $t_c$ is the chemical time. Instead of treating explicitly the chemical cross-term, we rationalize these different regimes by globalizing it as a production term involving a flux which depends on the rate at which the reaction zone is fed by the reactants, a formulation valid for $Da>1$ . For $Da<1$ , the reactants interpenetrate before they react, giving rise to a ‘diffusio-chemical’ regime where chemical production occurs within a substrate whose width is controlled by molecular diffusion.

Thermally Fluctuating, Semiflexible Sheets in Simple Shear Flow

We perform Brownian dynamics simulations of semiflexible colloidal sheets with hydrodynamic interactions and thermal fluctuations in shear flow. As a function of the ratio of bending rigidity to shear energy...

Oblique instability of a stratified oscillatory boundary layer

The instability and dynamics of a vertical oscillatory boundary layer in a container filled with a stratified fluid are addressed. Past experiments have shown that when the boundary oscillation frequency is of the same order as the buoyancy frequency, the system is unstable to a herringbone pattern of oblique waves. Prior studies assuming the basic state to be a unidirectional oscillatory shear flow were unable to account for the oblique waves. By accounting for confinement effects present in the experiments, and the ensuing three-dimensional structure of the basic state, we are able to numerically reproduce the experimental observations, opening the door to fully analysing the impacts of stratification on such boundary layers.

Buckling of elastic fibers in a shear flow

Three-dimensional dynamics of flexible fibers in shear flow are studied numerically, with a qualitative comparison to experiments. Initially, the fibers are straight, with different orientations with respect to the flow. By changing the rotation speed of a shear rheometer, we change the ratio A of bending to shear forces. We observe fibers in the flow-vorticity plane, which gives insight into the motion out of the shear plane. The numerical simulations of moderately flexible fibers show that they rotate along effective Jeffery orbits, and therefore the fiber orientation rapidly becomes very close to the flow-vorticity plane, on average close to the flow direction, and the fiber remains in an almost straight configuration for a long time. This ``ordering'' of fibers is temporary since they alternately bend and straighten out while tumbling. We observe numerically and experimentally that if the fibers are initially in the compressional region of the shear flow, they can undergo a compressional buckling, with a pronounced deformation of shape along their whole length during a short time, which is in contrast to the typical local bending that originates over a long time from the fiber ends. We identify differences between local and compressional bending and discuss their competition, which depends on the initial orientation of the fiber and the bending stiffness ratio A. There are two main finding. First, the compressional buckling is limited to a certain small range of the initial orientations, excluding those from the flow-vorticity plane. Second, since fibers straighten out in the flow-vorticity plane while tumbling, the compressional buckling is transient - it does not appear for times longer than 1/4 of the Jeffery period. For larger times, bending of fibers is always driven by their ends.

Some aspects in Kelvin-Helmholtz instability with and without Boussinesq approximation

The topic of this paper is the Kelvin-Helmholtz instability, a phenomenon which occurs on the interface of a stratified fluid, in the presence of a parallel shear flow, when there is a velocity and density difference across the interface of two adjacent layers. This paper focuses on a numerical simulation modelled by the Taylor-Goldstein equation, which represents a more realistic case compared to the basic Kelvin-Helmholtz shear flow. The Euler system is solved with new modelled smooth velocity and density profiles at the interface. The flux at cell boundaries is reconstructed by implementing a third order WENO (Weighted Essentially Non-Oscillatory) method. Next, a Riemann solver builds the fluxes at cell interfaces. The use of both Rusanov and HLLC solvers is investigated. Temporal discretization is done by applying the second order TVD (total variation diminishing) Runge-Kutta method on a uniform grid. Numerical simulations are performed comparatively for both Kelvin-Helmholtz and Taylor-Goldstein instabilities, on the same simulation domains. We find that increasing the number of grid points leads to a better accuracy in shear layer vortices visualization. Thus, we can conclude that applying the Taylor-Goldstein equation improves the realism in the general fluid instability modelling.