No Arabic abstract
Buoyant shear layers are encountered in many engineering and environmental applications and have been studied by researchers in the context of experiments and modeling for decades. Often, these flows have high Reynolds and Richardson numbers, and this leads to significant/intractable space-time resolution requirements for DNS or LES modeling. On the other hand, many of the important physical mechanisms in these systems, such as stress anisotropy, wake stabilization, and regime transition, inherently render eddy viscosity-based RANS modeling inappropriate. Accordingly, we pursue second-moment closure (SMC), i.e., full Reynolds stress/flux/variance modeling, for moderate Reynolds number non-stratified and stratified shear layers for which DNS is possible. A range of sub-model complexity is pursued for the diffusion of stresses, density fluxes and variance, pressure strain and scrambling, and dissipation. These sub-models are evaluated in terms of how well they are represented by DNS in comparison to the exact Reynolds averaged terms, and how well they impact the accuracy of the full RANS closure. For the non-stratified case, the SMC model predicts the shear layer growth rate and Reynolds shear stress profiles accurately. Stress anisotropy and budgets are captured only qualitatively. Comparing DNS of exact and modeled terms, inconsistencies in model performance and assumptions are observed, including inaccurate prediction of individual statistics, non-negligible pressure diffusion, and dissipation anisotropy. For the stratified case, shear layer and gradient Richardson number growth rates, and stress, flux, and variance decay rates, are captured with less accuracy than corresponding flow parameters in the non-stratified case. These studies lead to several recommendations for model improvement.
Instability of stratified multi-phase flow in a rotating platform becomes important because of a potential role in micro-mixing and micro-machines. Centrifugal actuation can play an important role in driving the flow and Coriolis force can enhance the mixing in a short span by destabilizing the flow. In this study, we focus on the impact of the Coriolis force on a rotating viscosity-stratified flow with a thin diffusive mixed layer between two fluid layers. Modal stability analysis is used to estimate the critical parameters, namely Rotation number, Reynolds number, and wave number, which are responsible to modulate the instability mechanism for different viscosity contrasts. Present study explores competing influences of rotational forces against the viscous and inertial forces. Correspondingly, rotational direction (clockwise/anticlockwise) shows a significant effect on the spatio-temporal instability mechanism and anticlockwise rotation promotes the instability. Usually, miscible viscosity stratified flow with respect to streamwise disturbance becomes more unstable for a thinner mixed layer. On the contrary, our numerical computation confirms a completely contrasting scenario, considering Coriolis force driven instability of a miscible system on account of spanwise disturbances. Possible physical mechanisms for the same are discussed in terms of base flow and energy fluctuation among perturbed and base flow. Comparison of two and three-dimensional instability is done to give a clear-cut idea about the linear instability of the flow system considered herein. Velocity and viscosity perturbation distributions display a critical bonding between the vortices near and away from mixed layer, which may be responsible for the variation of instability with respect to viscosity ratio and rotational direction.
The Lagrangian (LA) and Eulerian Acceleration (EA) properties of fluid particles in homogeneous turbulence with uniform shear and uniform stable stratification are studied using direct numerical simulations. The Richardson number is varied from $Ri=0$, corresponding to unstratified shear flow, to $Ri=1$, corresponding to strongly stratified shear flow. The probability density functions (pdfs) of both LA and EA have a stretched-exponential shape and they show a strong and similar influence on the Richardson number. The extreme values of the EA are stronger than those observed for the LA. Geometrical statistics explain that the magnitude of the EA is larger than its Lagrangian counterpart due to the mutual cancellation of the Eulerian and convective acceleration, as both vectors statistically show an anti-parallel preference. A wavelet-based scale-dependent decomposition of the LA and EA is performed. The tails of the acceleration pdfs grow heavier for smaller scales of turbulent motion. Hence the flatness increases with decreasing scale, indicating stronger intermittency at smaller scales. The joint pdfs of the LA and EA indicate a trend to stronger correlations with increasing Richardson number and at larger scales of the turbulent motion. A consideration of the terms in the Navier--Stokes equation shows that the LA is mainly determined by the pressure-gradient term, while the EA is dominated by the nonlinear convection term.
Exact solutions for laminar stratified flows of Newtonian/non-Newtonian shear-thinning fluids in horizontal and inclined channels are presented. An iterative algorithm is proposed to compute the laminar solution for the general case of a Carreau non-Newtonian fluid. The exact solution is used to study the effect of the rheology of the shear-thinning liquid on two-phase flow characteristics considering both gas/liquid and liquid/liquid systems. Concurrent and counter-current inclined systems are investigated, including the mapping of multiple solution boundaries. Aspects relevant to practical applications are discussed, such as the insitu hold-up, or lubrication effects achieved by adding a less viscous phase. A characteristic of this family of systems is that, even if the liquid has a complex rheology (Carreau fluid), the two-phase stratified flow can behave like the liquid is Newtonian for a wide range of operational conditions. The capability of the two-fluid model to yield satisfactory predictions in the presence of shear-thinning liquids is tested, and an algorithm is proposed to a priori predict if the Newtonian (zero shear rate viscosity) behaviour arises for a given operational conditions in order to avoid large errors in the predictions of flow characteristics when the power-law is considered for modelling the shear-thinning behaviour. Two-fluid model closures implied by the exact solution and the effect of a turbulent gas layer are also addressed.
Mixing layers can grow in time or space by vortex pairings that succeed each other in a nearly self-similar way. We use a point vortex model to study how confinement eventually limits this growth process, leading us to propose a wavelength selection mechanism for free shear layers with counterflow.
Linear stability of horizontal and inclined stratified channel flows of Newtonian/non-Newtonian shear-thinning fluids is investigated with respect to all wavelength perturbations. The Carreau model has been chosen for the modeling of the rheology of a shear-thinning fluid, owing to its capability to describe properly the constant viscosity limits (Newtonian behavior) at low and high shear rates. The results are presented in the form of stability boundaries on flow pattern maps (with the phases superficial velocities as coordinates) for several practically important gas-liquid and liquid-liquid systems. The stability maps are accompanied by spatial profiles of the critical perturbations, along with the distributions of the effective and tangent viscosities in the non-Newtonian layer, to show the influence of the complex rheological behavior of shear-thinning liquids on the mechanisms responsible for triggering instability. Due to the complexity of the considered problem, a working methodology is proposed to alleviate the search for the stability boundary. Implementation of the proposed methodology helps to reveal that in many cases the investigation of the simpler Newtonian problem is sufficient for the prediction of the exact (non-Newtonian) stability boundary of smooth stratified flow (i.e., in case of horizontal gas-liquid flow). Therefore, the knowledge gained from the stability analysis of Newtonian fluids is applicable to those (usually highly viscous) non-Newtonian systems. Since the stability of stratified flow involving highly viscous Newtonian liquids has not been researched in the literature, interesting findings on the viscosity effects are also obtained.