No Arabic abstract
We develop a framework for analyzing the momentum balance of laminar particle-laden flows based on immersed boundary methods, which solve the Navier-Stokes equations and resolve the particle surfaces. This framework differs from previous studies by explicitly accounting for the fluid inside the particles, which is a by-product of the immersed boundary method, allowing us to close the momentum balance for the flow around a single rolling sphere. We then compute a momentum balance of a laminar Poiseuille flow over a dense bed of particles, finding that the stresses remain in equilibrium even during unsteady flow conditions. While previous studies have focused on stresses for the streamwise momentum balance, the present approach also allows us to evaluate stress balances in the vertical direction, which are necessary to understand the role that collisions and hydrodynamic drag play during dilation and contraction of particle beds. While our analysis accounts for the fluid and particle phases separately, we attempt to establish a momentum balance for the fluid/particle mixture, but find that it does not completely close locally due to collision stresses not being resolved across the particle diameter. However, we find a correlation between the local shear rate and the gap in the mixture balance, which can potentially be used to close the balance for the mixture.
Rayleigh--Taylor fluid turbulence through a bed of rigid, finite-size, spheres is investigated by means of high-resolution Direct Numerical Simulations (DNS), fully coupling the fluid and the solid phase via a state-of-the art Immersed Boundary Method (IBM). The porous character of the medium reveals a totally different physics for the mixing process when compared to the well-known phenomenology of classical RT mixing. For sufficiently small porosity, the growth-rate of the mixing layer is linear in time (instead of quadratical) and the velocity fluctuations tend to saturate to a constant value (instead of linearly growing). We propose an effective continuum model to fully explain these results where porosity originated by the finite-size spheres is parameterized by a friction coefficient.
Unsteady laminar vortex shedding over a circular cylinder is predicted using a deep learning technique, a generative adversarial network (GAN), with a particular emphasis on elucidating the potential of learning the solution of the Navier-Stokes equations. Numerical simulations at two different Reynolds numbers with different time-step sizes are conducted to produce training datasets of flow field variables. Unsteady flow fields in the future at a Reynolds number which is not in the training datasets are predicted using a GAN. Predicted flow fields are found to qualitatively and quantitatively agree well with flow fields calculated by numerical simulations. The present study suggests that a deep learning technique can be utilized for prediction of laminar wake flow in lieu of solving the Navier-Stokes equations.
The impact of wall roughness on fully developed laminar pipe flow is investigated numerically. The roughness is comprised of square bars of varying size and pitch. Results show that the inverse relation between the friction factor and the Reynolds number in smooth pipes still persists in rough pipes, regardless of the rib height and pitch. At a given Reynolds number, the friction factor varies quadratically with roughness height and linearly with roughness pitch. We propose a single correlation for the friction factor that successfully collapses the data.
The mechanism of hydrodynamics-induced pairing of soft particles, namely closed bilayer membranes (vesicles, a model system for red blood cells) and drops, is studied numerically with a special attention paid to the role of the confinement (the particles are within two rigid walls). This study unveils the complexity of the pairing mechanism due to hydrodynamic interactions. We find both for vesicles and for drops that two particles attract each other and form a stable pair at weak confinement if their initial separation is below a certain value. If the initial separation is beyond that distance, the particles repel each other and adopt a longer stable interdistance. This means that for the same confinement we have (at least) two stable branches. To which branch a pair of particles relaxes with time depends only on the initial configuration. An unstable branch is found between these two stable branches. At a critical confinement the stable branch corresponding to the shortest interdistance merges with the unstable branch in the form of a saddle-node bifurcation. At this critical confinement we have a finite jump from a solution corresponding to the continuation of the unbounded case to a solution which is induced by the presence of walls. The results are summarized in a phase diagram, which proves to be of a complex nature. The fact that both vesicles and drops have the same qualitative phase diagram points to the existence of a universal behavior, highlighting the fact that with regard to pairing the details of mechanical properties of the deformable particles are unimportant. This offers an interesting perspective for simple analytical modeling.
We study theoretically and experimentally how a thin layer of liquid flows along a flexible beam. The flow is modelled using lubrication theory and the substrate is modelled as an elastica which deforms according to the Euler-Bernoulli equation. A constant flux of liquid is supplied at one end of the beam, which is clamped horizontally, while the other end of the beam is free. As the liquid film spreads, its weight causes the beam deflection to increase, which in turn enhances the spreading rate of the liquid. This feedback mechanism causes the front position ${sigma}$(t) and the deflection angle at the front ${phi}$(t) to go through a number of different power-law behaviours. For early times, the liquid spreads like a horizontal gravity current, with ${sigma}$(t) = $t^{4/5}$ and ${phi}$(t) = $t^{13/5}$. For intermediate times, the deflection of the beam leads to rapid acceleration of the liquid layer, with ${sigma}$(t) = $t^4$ and ${phi}$(t) = $t^9$. Finally, when the beam has sagged to become almost vertical, the liquid film flows downward with ${sigma}$(t) = $t$ and ${phi}$(t) ~ ${pi}$/2. We demonstrate good agreement between these theoretical predictions and experimental results.