Effect of intermediate fluid layer on Kelvin–Helmholtz instability

2018 ◽  
Vol 96 (10) ◽  
pp. 1145-1154 ◽  
Author(s):  
Ying Zhang ◽  
Wenqiang Shang ◽  
Mengjun Yao ◽  
Boheng Dong ◽  
Peisheng Li
Keyword(s):  
Growth Rate ◽  
Stokes Equations ◽  
Fluid Layer ◽  
Immiscible Fluids ◽  
Navier Stokes ◽  
Helmholtz Instability ◽  
Flow Configuration ◽  
Front Tracking Method ◽  
Internal Disturbance ◽  
Over Time

Two-dimensional K-H (Kelvin–Helmholtz) instability of the three-component immiscible fluids with an intermediate fluid layer is numerically studied using the front-tracking method (FTM). The instability is governed by the Navier–Stokes equations and the conservation of mass equation for incompressible flow. A finite difference method is used to discretize the governing system. This study focuses on the influence of flow configuration, the thickness of intermediate fluid layer and the distribution of intermediate fluid layer on K-H instability. It is shown that the larger the initial horizontal velocity difference is, the faster the internal disturbance increases, and the characteristic form of K-H instability becomes more obvious for different flow configuration. It is also observed that the thickness of the intermediate fluid layer is negatively correlated to the billow height and the numerical growth rate. In addition, when the intermediate fluid layer is thicker than 0.4 times the disturbance wavelength, the billow height and the numerical growth rate for the K-H instability of the upper and lower interfaces change over time synchronously. The higher the initial height of the lower interface is, the greater the growth rate and billow height of the upper interface are. Besides, the upper and lower interfaces are rolled up synchronously over time when the intermediate fluid layer is symmetrically distributed with y = 0.5 in the fluid system.

Pontryagin maximum principle and second order optimality conditions for optimal control problems governed by 2D nonlocal Cahn–Hilliard–Navier–Stokes equations

Analysis ◽  
2020 ◽  
Vol 40 (3) ◽  
pp. 127-150
Author(s):  
Tania Biswas ◽  
Sheetal Dharmatti ◽  
Manil T. Mohan
Keyword(s):  
Optimal Control ◽  
Optimal Control Problem ◽  
Control Problem ◽  
Stokes Equations ◽  
Minimum Principle ◽  
Immiscible Fluids ◽  
Second Order ◽  
Navier Stokes ◽  
Distributed Optimal Control ◽  
Navier Stokes Equations

AbstractIn this paper, we formulate a distributed optimal control problem related to the evolution of two isothermal, incompressible, immiscible fluids in a two-dimensional bounded domain. The distributed optimal control problem is framed as the minimization of a suitable cost functional subject to the controlled nonlocal Cahn–Hilliard–Navier–Stokes equations. We describe the first order necessary conditions of optimality via the Pontryagin minimum principle and prove second order necessary and sufficient conditions of optimality for the problem.

scholarly journals Numerical simulation of Faraday waves

2009 ◽  
Vol 635 ◽  
pp. 1-26 ◽  
Author(s):  
NICOLAS PÉRINET ◽  
DAMIR JURIC ◽  
LAURETTE S. TUCKERMAN
Keyword(s):  
Stokes Equations ◽  
Three Dimensional ◽  
Oscillation Period ◽  
Finite Amplitude ◽  
Three Dimensions ◽  
Navier Stokes ◽  
Faraday Waves ◽  
Temporal Behaviour ◽  
Two Fluids ◽  
Front Tracking Method

We simulate numerically the full dynamics of Faraday waves in three dimensions for two incompressible and immiscible viscous fluids. The Navier–Stokes equations are solved using a finite-difference projection method coupled with a front-tracking method for the interface between the two fluids. The critical accelerations and wavenumbers, as well as the temporal behaviour at onset are compared with the results of the linear Floquet analysis of Kumar & Tuckerman (J. Fluid Mech., vol. 279, 1994, p. 49). The finite-amplitude results are compared with the experiments of Kityk et al (Phys. Rev. E, vol. 72, 2005, p. 036209). In particular, we reproduce the detailed spatio-temporal spectrum of both square and hexagonal patterns within experimental uncertainty. We present the first calculations of a three-dimensional velocity field arising from the Faraday instability for a hexagonal pattern as it varies over its oscillation period.

scholarly journals Numerical stability analysis of a vortex ring with swirl

2019 ◽  
Vol 878 ◽  
pp. 5-36 ◽  
Author(s):  
Yuji Hattori ◽  
Francisco J. Blanco-Rodríguez ◽  
Stéphane Le Dizès
Keyword(s):  
Numerical Simulation ◽  
Growth Rate ◽  
Direct Numerical Simulation ◽  
Vortex Ring ◽  
Stokes Equations ◽  
Base Flow ◽  
Critical Layer ◽  
Navier Stokes ◽  
Instability Mode ◽  
Navier Stokes Equations

The linear instability of a vortex ring with swirl with Gaussian distributions of azimuthal vorticity and velocity in its core is studied by direct numerical simulation. The numerical study is carried out in two steps: first, an axisymmetric simulation of the Navier–Stokes equations is performed to obtain the quasi-steady state that forms a base flow; then, the equations are linearized around this base flow and integrated for a sufficiently long time to obtain the characteristics of the most unstable mode. It is shown that the vortex rings are subjected to curvature instability as predicted analytically by Blanco-Rodríguez & Le Dizès (J. Fluid Mech., vol. 814, 2017, pp. 397–415). Both the structure and the growth rate of the unstable modes obtained numerically are in good agreement with the analytical results. However, a small overestimation (e.g. 22 % for a curvature instability mode) by the theory of the numerical growth rate is found for some instability modes. This is most likely due to evaluation of the critical layer damping which is performed for the waves on axisymmetric line vortices in the analysis. The actual position of the critical layer is affected by deformation of the core due to the curvature effect; as a result, the damping rate changes since it is sensitive to the position of the critical layer. Competition between the curvature and elliptic instabilities is also investigated. Without swirl, only the elliptic instability is observed in agreement with previous numerical and experimental results. In the presence of swirl, sharp bands of both curvature and elliptic instabilities are obtained for $\unicode[STIX]{x1D700}=a/R=0.1$, where $a$ is the vortex core radius and $R$ the ring radius, while the elliptic instability dominates for $\unicode[STIX]{x1D700}=0.18$. New types of instability mode are also obtained: a special curvature mode composed of three waves is observed and spiral modes that do not seem to be related to any wave resonance. The curvature instability is also confirmed by direct numerical simulation of the full Navier–Stokes equations. Weakly nonlinear saturation and subsequent decay of the curvature instability are also observed.

scholarly journals A ternary Cahn–Hilliard–Navier–Stokes model for two-phase flow with precipitation and dissolution

2020 ◽  
pp. 1-35
Author(s):  
Christian Rohde ◽  
Lars von Wolff
Keyword(s):  
Phase Field ◽  
Stokes Equations ◽  
Solid Phase ◽  
Immiscible Fluids ◽  
Sharp Interface ◽  
Ion Concentration ◽  
Navier Stokes ◽  
Interface Model ◽  
Sharp Interface Model ◽  
Two Phase

We consider the incompressible flow of two immiscible fluids in the presence of a solid phase that undergoes changes in time due to precipitation and dissolution effects. Based on a seminal sharp interface model a phase-field approach is suggested that couples the Navier–Stokes equations and the solid’s ion concentration transport equation with the Cahn–Hilliard evolution for the phase fields. The model is shown to preserve the fundamental conservation constraints and to obey the second law of thermodynamics for a novel free energy formulation. An extended analysis for vanishing interfacial width reveals that in this limit the sharp interface model is recovered, including all relevant transmission conditions. Notably, the new phase-field model is able to realize Navier-slip conditions for solid–fluid interfaces in the limit.

scholarly journals Effects of surfactant on propagation and rupture of a liquid plug in a tube

2019 ◽  
Vol 872 ◽  
pp. 407-437 ◽  
Author(s):  
M. Muradoglu ◽  
F. Romanò ◽  
H. Fujioka ◽  
J. B. Grotberg
Keyword(s):  
Shear Stress ◽  
Evolution Equations ◽  
Stokes Equations ◽  
Thin Liquid Film ◽  
Front Tracking ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Liquid Plug ◽  
Front Tracking Method ◽  
Airway Walls

Surfactant-laden liquid plug propagation and rupture occurring in lower lung airways are studied computationally using a front-tracking method. The plug is driven by an applied constant pressure in a rigid axisymmetric tube whose inner surface is coated by a thin liquid film. The evolution equations of the interfacial and bulk surfactant concentrations coupled with the incompressible Navier–Stokes equations are solved in the front-tracking framework. The numerical method is first validated for a surfactant-free case and the results are found to be in good agreement with the earlier simulations of Fujioka et al. (Phys. Fluids, vol. 20, 2008, 062104) and Hassan et al. (Intl J. Numer. Meth. Fluids, vol. 67, 2011, pp. 1373–1392). Then extensive simulations are performed to investigate the effects of surfactant on the mechanical stresses that could be injurious to epithelial cells, such as pressure and shear stress. It is found that the liquid plug ruptures violently to induce large pressure and shear stress on airway walls and even a tiny amount of surfactant significantly reduces the pressure and shear stress and thus improves cell survivability. However, addition of surfactant also delays the plug rupture and thus airway reopening.

Dynamics of the interface between miscible liquids subjected to horizontal vibration

2015 ◽  
Vol 784 ◽  
pp. 342-372 ◽  
Author(s):  
Y. A. Gaponenko ◽  
M. Torregrosa ◽  
V. Yasnou ◽  
A. Mialdun ◽  
V. Shevtsova
Keyword(s):  
Time Scales ◽  
Stokes Equations ◽  
Immiscible Fluids ◽  
Interfacial Instability ◽  
Navier Stokes ◽  
Layer System ◽  
Navier Stokes Equations ◽  
Horizontal Vibration ◽  
Natural Time ◽  
Miscible Liquids

We present experimental evidence of the existence of an interfacial instability between two miscible liquids of similar (but non-identical) viscosities and densities under horizontal vibration. A stably stratified two-layer system is composed of the same binary mixture with different concentrations placed in a confined cell (with length twice as large as the height). Unlike the case of immiscible fluids, here, the interface is a transient layer of small but non-zero thickness. In the experiments, the frequency and amplitude were varied within the ranges 2–24 Hz and 1.5–16 mm, respectively. When the value of the oscillatory forcing increases, the amplitudes of the interface perturbations grow continuously, forming a saw-tooth frozen structure. This evolution is also examined numerically. In addition to the solutions of full 3-D Navier–Stokes equations, an averaging approach with separation of time scales is used for situations in which the forcing period is very small compared to the natural time scales of the problem. The simulation of averaged equations provides the explanation of the instability development, the calculations of the full nonlinear equations shed light on the decay of a wavy pattern. The results of numerical modelling perfectly support the experimental observations.

Normal modes for a stratified viscous fluid layer

2002 ◽  
Vol 132 (3) ◽  
pp. 611-625 ◽  
Author(s):  
K. F. Gurski ◽  
R. L. Pego
Keyword(s):  
Stokes Equations ◽  
Fluid Layer ◽  
Normal Modes ◽  
Internal Gravity Waves ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Damped Oscillations ◽  
Small Constant ◽  
Constant Viscosity ◽  
Stratified Viscous Fluid

We consider internal gravity waves in a stratified fluid layer with rigid horizontal boundaries and periodic boundary conditions on the sides at constant temperature with a small constant viscosity, modelled using the incompressible Navier-Stokes equations. Using operator-theoretic methods to study the damping rates of internal waves we prove there are non-oscillatory wave modes with arbitrarily small damping rates. We provide an asymptotic approximation for these non-oscillatory modes. Additionally, we find that the eigenvalues for damped oscillations are in an explicitly describable half-ring.

scholarly journals Dynamics of a contracting fluid compound filament with a variable density ratio

2021 ◽  
Vol 24 (2) ◽  
pp. first
Author(s):  
Truong V. Vu ◽  
Vinh T. Nguyen ◽  
Phan H. Nguyen ◽  
Nang X. Ho ◽  
Binh D. Pham ◽  
...  
Keyword(s):  
Interfacial Tension ◽  
Stokes Equations ◽  
Weighting Function ◽  
Density Ratio ◽  
Variable Density ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Front Tracking Method ◽  
Fixed Grid ◽  
Eulerian Grid

Introduction: Compound fluid filaments appear in many applications, e.g., drug delivery and processing or microfluidic systems. This paper focuses on the numerical simulation of an incompressible, immiscible, and Newtonian fluid for the contraction process of a fluid compound filament by solving the Navier-Stokes equations. The front-tracking method is used to solve this problem, which uses connected segments (Lagrangian grid) that move on a fixed grid (Eulerian grid) to represent the interface between the liquids. Methods: The interface points are advected by the velocity interpolated from those of the fixed grid using the area weighting function. The coordinates of the interface points are used to construct the indicators specifying the different fluids and compute the interfacial tension force. Results: The simulation results show that under the effects of the interfacial tension, the capsuleshaped filament can transform into a spherical compound droplet (i.e., non-breakup) or can break up into smaller spherical compound and simple droplets (i.e., breakup). When the density ratio of the outer to middle fluids increases, the filament changes from non-breakup to breakup upon contraction. Conclusion: Increasing the density ratio enhances the breakup of the compound filament during contraction. The breakup is also promoted by increasing the initial length of the filament.

scholarly journals Direct numerical simulation of solidification with effects of density difference

2016 ◽  
Vol 38 (3) ◽  
pp. 193-204
Author(s):  
Vu Van Truong
Keyword(s):  
Stokes Equations ◽  
Solid Phase ◽  
Density Difference ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Front Tracking Method ◽  
Liquid Phases ◽  
Solidification Processes ◽  
Solid Liquid Interface ◽  
Effects Of Density

In this paper, direct numerical simulations are presented for solidification with the effects of density difference between the solid and liquid phases. A front-tracking method is used. The solidification front, i.e. the solid-liquid interface separating solid and liquid, is represented by connected elements that move on a rectangular and stationary grid. The Navier-Stokes equations are solved by a projection method on the entire domain including the solid phase. An indicator function reconstructed from the front information is used to set the velocities in the solid phase to zero, and thus to enforce the no-slip condition at the interface. The method is validated through comparisons with exact solutions for one- and two-dimensional problems. The method is then used to simulate the solidification processes with the effects of volume change due to density difference

Sign in / Sign up

Export Citation Format

Share Document