No Arabic abstract
In this investigation we revisit the concept of effective free surfaces arising in the solution of the time-averaged fluid dynamics equations in the presence of free boundaries. This work is motivated by applications of the optimization and optimal control theory to problems involving free surfaces, where the time-dependent formulations lead to many technical difficulties which are however alleviated when steady governing equations are used instead. By introducing a number of precisely stated assumptions, we develop and validate an approach in which the interface between the different phases, understood in the time-averaged sense, is sharp. In the proposed formulation the terms representing the fluctuations of the free boundaries and of the hydrodynamic quantities appear as boundary conditions on the effective surface and require suitable closure models. As a simple model problem we consider impingement of free-falling droplets onto a fluid in a pool with a free surface, and a simple algebraic closure model is proposed for this system. The resulting averaged equations are of the free-boundary type and an efficient computational approach based on shape optimization formulation is developed for their solution. The computed effective surfaces exhibit consistent dependence on the problem parameters and compare favorably with the results obtained when the data from the actual time-dependent problem is used in lieu of the closure model.
In this paper, the instability of layered two-phase flows caused by the presence of a soluble surfactant (or a surface active solute) is studied. The fluids have different viscosities, but are density matched to focus on Marangoni effects. The fluids flow between two flat plates, which are maintained at different solute concentrations. This establishes a constant flux of solute from one fluid to the other in the base state. A linear stability analysis is performed, using a combination of asymptotic and numerical methods. In the creeping flow regime, Marangoni stresses destabilize the flow, provided a concentration gradient is maintained across the fluids. One long wave and two short wave Marangoni instability modes arise, in different regions of parameter space. A well-defined condition for the long wave instability is determined in terms of the viscosity and thickness ratios of the fluids, and the direction of mass transfer. Energy budget calculations show that the Marangoni stresses that drive long and short wave instabilities have distinct origins. The former is caused by interface deformation while the latter is associated with convection by the disturbance flow. Consequently, even when the interface is non-deforming (in the large interfacial tension limit), the flow can become unstable to short wave disturbances. On increasing $Re$, the viscosity-induced interfacial instability comes into play. This mode is shown to either suppress or enhance the Marangoni instability, depending on the viscosity and thickness ratios. This analysis is relevant to applications such as solvent extraction in microchannels, in which a surface-active solute is transferred between fluids in parallel stratified flow. It is also applicable to the thermocapillary problem of layered flow between heated plates.
We study numerically joint mixing of salt and colloids by a chaotic velocity field $mathbf{V}$, and how salt inhomogeneities accelerate or delay colloid mixing by inducing a velocity drift $mathbf{V}_{rm dp}$ between colloids and fluid particles as proposed in recent experiments cite{Deseigne2013}. We demonstrate that because the drift velocity is no longer divergence free, small variations to the total velocity field drastically affect the evolution of colloid variance $sigma^2=langle C^2 rangle - langle C rangle^2$. A consequence is that mixing strongly depends on the mutual coherence between colloid and salt concentration fields, the short time evolution of scalar variance being governed by a new variance production term $P=- langle C^2 abla cdot mathbf{V}_{rm dp} rangle/2$ when scalar gradients are not developed yet so that dissipation is weak. Depending on initial conditions, mixing is then delayed or enhanced, and it is possible to find examples for which the two regimes (fast mixing followed by slow mixing) are observed consecutively when the variance source term reverses its sign. This is indeed the case for localized patches modeled as gaussian concentration profiles.
We investigate the effective rheology of two-phase flow in a bundle of parallel capillary tubes carrying two immiscible fluids under an external pressure drop. The diameter of each tube varies along its length and the corresponding capillary threshold pressures are considered to be distributed randomly according to a uniform probability distribution. We demonstrate through analytical calculations that a transition from a linear Darcy regime to a non-linear behavior occurs while decreasing the pressure drop $Delta P$, where the total flow rate $langle Q rangle$ varies with $Delta P$ with an exponent $2$. This exponent for the non-linear regime changes when a lower cut-off $P_m$ is introduced in the threshold distribution. We demonstrate analytically that, in the limit where $Delta P$ approaches $P_m$, the flow rate scales as $langle Q rangle sim (|Delta P|-P_m)^{3/2}$. We have also provided some numerical results in support to our analytical findings.
Most fluid flow problems that are vital in engineering applications involve at least one of the following features: turbulence, shocks, and/or material interfaces. While seemingly different phenomena, these flows all share continuous generation of high wavenumber modes, which we term the $k_infty$ irregularity. In this work, an inviscid regularization technique called observable regularization is proposed for the simulation of two-phase compressible flows. The proposed approach regularizes the equations at the level of the partial differential equation and as a result, any numerical method can be used to solve the system of equations. The regularization is accomplished by introducing an observability limit that represents the length scale below which one cannot properly model or continue to resolve flow structures. An observable volume fraction equation is derived for capturing the material interface, which satisfies the pressure equilibrium at the interface. The efficacy of the observable regularization method is demonstrated using several test cases, including a one-dimensional material interface tracking, one-dimensional shock-tube and shock-bubble problems, and two-dimensional simulations of a shock interacting with a cylindrical bubble. The results show favorable agreement, both qualitatively and quantitatively, with available exact solutions or numerical and experimental data from the literature. The computational saving by using the current method is estimated to be about one order of magnitude in two-dimensional computations and significantly higher in three-dimensional computations. Lastly, the effect of the observability limit and best practices to choose its value are discussed.
We develop an adversarial-reinforcement learning scheme for microswimmers in statistically homogeneous and isotropic turbulent fluid flows, in both two (2D) and three dimensions (3D). We show that this scheme allows microswimmers to find non-trivial paths, which enable them to reach a target on average in less time than a naive microswimmer, which tries, at any instant of time and at a given position in space, to swim in the direction of the target. We use pseudospectral direct numerical simulations (DNSs) of the 2D and 3D (incompressible) Navier-Stokes equations to obtain the turbulent flows. We then introduce passive microswimmers that try to swim along a given direction in these flows; the microswimmers do not affect the flow, but they are advected by it.