European Journal of Mechanics B-fluids最新文献

Wave breaking is a multifaceted physical phenomenon that is not fully understood and remains challenging to model. An effective method for investigating wave breaking involves utilising the two-phase Reynolds-averaged Navier–Stokes (RANS) equations to directly simulate breaking waves. In this study, we apply a RANS model with an adaptively refined mesh to simulate breaking waves in deep water using the stabilised RANS model proposed by Larsen and Fuhrman. This approach enables a more efficient simulation of the physics of breaking waves compared to Direct Numerical Simulations, as it places less stringent demands on grid resolution. Our findings demonstrate that the RANS model compares well with deep water wave breaking experiments in terms of surface elevation. We also give estimates of the breaking strength parameter of our RANS simulations and compared them with the literature.

We study theoretically the steady Newtonian flow in confined and hyperbolic long tubes (symmetric channels and axisymmetric pipes) considering slip along the walls. Using a stream function formulation, and the extended (or high-order) lubrication method in terms of the square of the aspect ratio of the tube, ε, the solution for the stream function is found analytically up to twentieth order in ε. At the classic lubrication limit, i.e. i.e. for a vanishing small aspect ratio, and for perfect slip conditions, the analysis predicts a plug-like velocity profile and a constant strain-rate on the midplane/axis of symmetry of the tube. A constant strain-rate is also predicted for the non-slip case. Furthermore, the high order asymptotic results for the stream function and fluid velocity are post-processed with an acceleration technique to investigate the convergence and accuracy of the solution. The results reveal the existence of a boundary layer at the inlet of the tube, the influence of which diminishes in a very short distance from the entrance. We discuss the effect of the contraction ratio of the tube and the dimensionless slip coefficient on the midplane/centerline and wall (slip) velocities, as well as on the average pressure-drop, required to maintain a constant flow-rate. The acceleration of converge technique on the solution for the pressure-drop revealed a remarkable convergence at a value slightly larger (∼1 %) than the value predicted by the classic lubrication theory. Finally, we comment on the common practice in the literature for approaching the velocity profile with the velocity profile at the classic lubrication limit, and we compare the high-order results for the strain rate at the midplane/centerline with the effective strain rate previously derived in the literature by Housiadas & Beris, J. Rheology, 68(3), 327–339, 2024.

This paper presents mathematical modelling and simulation of thin free-surface flows of viscoplastic fluids with a Herschel–Bulkley rheology over complex topographies with basal perturbations. Using the asymptotic expansion method, depth-averaged models (lubrication and shallow water type models) are derived for 3D (three-dimensional) multi-regime flows on non-flat inclined topographies with varying basal slipperiness. Starting from the Navier–Stokes equations, two flow regimes corresponding to different balances between shear and pressure forces are presented. Flow models corresponding to these regimes are calculated as perturbations of the zeroth-order solutions. The classical reference models in the literature are recovered by considering their respective cases on a flat-inclined surface. In the second regime case, a pressure term is non-negligible. Mathematically, it leads to a corrective term to the classical regime equations. Flow solutions of the two regimes are compared; the difference appears in particular in the vicinity of sharp changes of slopes. Nonetheless, both regime models are compared with experiments and are found to be in good agreement. Furthermore, numerical examples are shown to illustrate the robustness of the present shallow water models to simulate viscoplastic flows in 3D and over an inclined topography with local perturbations in basal elevation and basal slipperiness. The derived models are adequate for direct (engineering and geophysical) applications to real-world flow problems presenting Herschel–Bulkley rheology like lava and mud flows.

We compared the accuracy of volume of fluid (VOF) methods in unstructured solvers using the following five different methods: 1 - the algebraically compressive VOF method, 2 – simple coupled VOF method with Level Set (S-CLSVOF) method, 3 - interface-compressing VOF method incorporated with Laplacian filter (VOFL), 4 - isoAdvector method, and 5 - isoAdvector method incorporated with Laplacian filter (isoAdvectorL) by incorporating them into OpenFOAM®, an open-source software. To evaluate these methods under proper conditions, we compared the calculation accuracy using the optimized parameters, which are explored by Bayesian optimization. The test cases for advection accuracy of volume fraction and for imbalance of surface tension force in static multiphase fluid fields were considered. In this study, we found that the compression parameters and maximum Courant number should be adjusted to obtain high accuracy simulation according to the simulation condition in VOF and S-CLSVOF method. In VOFL and isoAdvectorL methods, the spurious current can be extremely reduced, which means that these methods are suitable for slow flow with higher Laplace number conditions.

In fluid dynamics, the Atwood number is a dimensionless parameter that quantifies the density difference between two fluids. It is calculated as $At=({\rho }_1-{\rho }_2)/({\rho }_1+{\rho }_2)$, where ${\rho }_1$ and ${\rho }_2$ represent the densities of the respective fluids. This research employs high-fidelity numerical simulations to examine the Atwood number impacts on Richtmyer–Meshkov (RM) flows triggered by a shocked forward-pentagonal bubble. Five distinct gases — ${\text{SF}}_6$, Kr, Ar, Ne, and He — are considered within the forward-pentagonal bubble, encompassed by ${\text{N}}_2$ gas. In these simulations, a third-order discontinuous Galerkin approach is applied to solve a two-dimensional set of compressible Navier–Stokes-Fourier (NSF) equations for two-component gas flows. To discretize space, hierarchical modal basis functions based on orthogonal-scaled Legendre polynomials are employed. This approach simplifies the NSF equations into a set of ordinary differential equations over time, which are solved using an explicit third-order SSP Runge–Kutta algorithm. The numerical results highlight the notable impact of the Atwood number on the evolution of RM flows in the shocked forward-pentagonal bubble, a phenomenon not previously reported in the literature. The Atwood number exerts a significant influence on the flow patterns, leading to intricate wave formations, shock focusing, jet generation, and interface distortion. Moreover, a comprehensive analysis of the these impact elucidates the mechanisms driving vorticity formation during the interaction process. Additionally, the study conducts a thorough quantitative examination of the Atwood number impacts on the flow fields based on integral quantities and interface features.

Diffraction is a fundamental phenomenon that occurs when blast or shock waves pass over sudden discontinuous surfaces. It generates a complex flow field consisting of diffracted waves, expansion waves, slipstream, contact surface, and an unstable shear layer, in addition to emitting acoustic waves. In this study, we investigated the diffraction of a blast wave passing over a series of grooved structures with different aspect ratios and geometrical shapes (rectangular, circular, and triangular) using high-speed shadowgraph images. The blast wave Mach number considered in our investigation is 1.34. The grooves feature leading-edge geometrical variations such as rectangular, circular arc, and wedge shapes positioned at various lateral locations from the exit of the shock tube. The aspect ratios of the rectangular grooves vary from 0.33, 0.5, and 0.67. The circular and triangular grooves have an aspect ratio of 0.33. The trajectories and velocities of the primary vortex are calculated by tracking the location of the vortex in the shadowgraph images. Our observations revealed that a large portion of the incident blast wave is abducted inside the groove as the aspect ratio increases in rectangular grooves, resulting in better attenuation of the blast wave. The grooves, which have circular shapes, produced weaker diffraction, which resulted in delayed and weak primary vortex. The triangular grooves produced the strongest primary vortex and have the highest attenuation characteristics among other grooves. The strength and trajectory of the primary vortex formed over the grooves strongly depend on the aspect ratio and the curvature of the leading edge for a given Mach number. Vortices generated from rectangular and triangular grooves exhibit considerable strength and longevity.

An airlift pump is a vertical tube that utilizes the buoyant effects of a gas to lift a liquid. Unlike a standard mechanical pump, the liquid flow rate through the airlift pump is not directly controlled; rather, it depends on the supplied gas flow rate, the tube length and diameter, and the relative height of the liquid supply free surface (submergence ratio). The present study uses the commercial CFD code ANSYS CFX to model the isothermal, 3D, transient flow in an airlift pump using water and air. The model applies pressure boundary conditions at both ends of the tube and specifies the mass flow rate of air through multiple openings in the side of the tube. The bottom of the tube is an inlet of water only and the outlet is a two-phase flow opening. A time-dependent, homogeneous, VOF two-phase RANS CFD modelling approach is used with the air treated as an ideal gas. This work found that a complete 3D domain was necessary for consistent prediction of the airlift performance and physically realistic two-phase flow structures. Statistical analysis of the two-phase flow structures was applied to characterize airlift pump instability and better understand the physics of the airlift pump.

Transport of particles in flows is often modeled in a combined Eulerian–Lagrangian framework. The flow is evaluated on an Eulerian grid, while particles are modeled as Lagrangian points whose positions and velocities are evolved in time, resulting in particle trajectories embedded in the time-dependent flow field. The method essentially resolves the flow field in complex geometries in detail but uses a closure model for the particle dynamics aimed at including the essential particle–fluid interactions at the cost of detailed small-scale physics. Rarefaction effects are usually included through the phenomenological Cunningham correction on the drag force experienced by the particles. In this Lagrangian point-particle approach, any explicit reference to the finite size and the shape of the particles, and their local orientation in the flow field, is typically ignored. In this work we aim to address this gap by deriving, from fully-resolved Direct Simulation Monte Carlo (DSMC) studies, heuristic or closure models for the drag force acting on prolate and oblate spheroidal particles with different aspect ratios, and a fixed orientation, in uniform ambient rarefied flows covering the transition regime between the continuum and free-molecular limits. These closure models predict the drag in the transition regime for all considered parameter settings (validated with DSMC data). The continuum limit is enforced a priori and we retrieve the free-molecular limit with reasonable accuracy (based on comparisons with literature data). We also include in the models the capability to predict effects related to basic gas-surface interactions via the tangential momentum accommodation coefficient. We furthermore assess the validity of the proposed closure model for particle dynamics in proximity to solid walls. This investigation extends our previous work, which focused on small aspect ratio spheroids with exclusively diffusive gas-surface interactions [see Livi et al. (2022)]. The derived models are obtained for isothermal, subsonic flows relevant for particle contamination control in semiconductor manufacturing.