Surface tension effects on interaction between two fluids near a wall

Interaction between two fluids near a fixed solid surface ismodelled, with surface tension acting as an important influenceon the assumed planar motion. The two fluids are immiscible,incompressible and have small density and viscosity ratios;the heavier more viscous body of fluid is approaching the solidsurface and the other fluid is lying as a thin layer in between.In the so-called supercritical range where, for both fluids,inviscid forces dominate over viscous ones, a pair of pressure–shaperelations is found which leads to a nonlinear integro-differentialequation for the unknown interface shape. Analysis, computationand comparisons are applied to the equation. Travelling-statesolutions are found of periodic and non-periodic form, includinginteresting cases which exhibit parabolic growth of the layerthickness in the far field.     

Analytic solutions for unsteady flows over a superhydrophobic surface

A class of unsteady analytic solutions for the evolution of viscous Newtonian fluids over a superhydrophobic surface is derived. The surface is represented by regular rectilinear riblets and the fluid is assumed to flow along the grooves’ direction. A mixed boundary condition is imposed on the surface, consisting in homogeneous Neumann conditions over the riblet voids and homogeneous Dirichlet conditions on the wall intervals. The transition between the above conditions is modelled through a Robin condition with non-constant smooth coefficients in a completely general manner. A global solution is derived and relevant examples, that can be fruitfully adopted as benchmark solutions for testing numerical solvers, are discussed.     

The Immersed Interface Method for Simulating Two-Fluid Flows

We develop the immersed interface method (IIM) to simulate a two-fluid flow of two immiscible fluids with different density and viscosity. Due to the surface tension and the discontinuous fluid properties, the two-fluid flow has nonsmooth velocity and discontinuous pressure across the moving sharp interface separating the two fluids. The IIM computes the flow on a fixed Cartesian grid by incorporating into numerical schemes the necessary jump conditions induced by the interface. We present how to compute these necessary jump conditions from the analytical principal jump conditions derived in [Xu, DCDS, Supplement 2009, pp. 838-845]. We test our method on some canonical two-fluid flows. The results demonstrate that the method can handle large density and viscosity ratios, is second-order accurate in the infinity norm, and conserves mass inside a closed interface.     

A meshless numerical approach based on Integrated Radial Basis Functions and level set method for interfacial flows

This paper reports a new meshless Integrated Radial Basis Function Network (IRBFN) approach to the numerical simulation of interfacial flows in which the two-way interaction between a moving interface and the ambient viscous flow is fully investigated. When an interface between two immiscible fluids moves, not only its position and shape but also the flow variables (i.e. velocity field and pressure) change due to the presence of surface tension along the moving interface. The velocity field of the ambient flow, on the other hand, causes the interface to move and deform as a result of momentum transport between the two immiscible fluids on both sides of the interface. Numerical investigations of such a two-way interaction is reported in this paper where the level set method is used in combination with high-order projection schemes in the meshless framework of the IRBFN method. Numerical investigations on the meshless projection schemes are performed with typical benchmark incompressible viscous flow problems for verification purposes. The approach is then demonstrated with the numerical simulation of two bubbles moving, stretching and merging in an incompressible ambient fluid under the action of buoyancy force.     

Generalized viscous shock layer equations with slip conditions on a surface in a flow field and on a head shock

The plane and axisymmetric problems of super- and hypersonic flow of a homogeneous viscous heat-conducting perfect gas over a blunt body are considered. Generalized viscous shock layer equations that take into account all the second-order effects of boundary-layer theory, i.e., the terms O(Re^?1/2), are derived from the Navier–Stokes equations by the asymptotic method, and all the out-of-order third-order terms O(Re^?1) and higher-order terms are also retained, except terms with second derivations in the marching coordinate (Re is Reynolds number, determined from the free-stream density and velocity the linear dimension, which is equal to the nose radius of the blunt Body, and the free-stream shear viscosity at the stagnation temperature). Thus, only the presence of terms with second derivatives in the marching coordinate, which specify the elliptical properties of the complete system of Navier–Stokes equations, distinguish it from the generalized viscous shock layer equations, which do not contain these terms. Slip and a temperature jump conditions on a body surface are presented with the same degree of accuracy, and generalized Rankine–Hugoniot conditions on a head shock, which take into account the effects of the viscosity and heat conduction, including their influence on the determination of the pressure, are derived. The incorrect and unfounded approximations used in preceding studies and the efficiency of iterative marching techniques for solving the generalized viscous shock layer equations, as well as the ability of the latter to provide a correct solution for the drag and heat-transfer coefficients in the transitional flow regime if the solution is constructed taking the slip and temperature jump on a surface and on a head shock into account, are noted.     

The applicability of continuum models in the transitional regime of hypersonic flow over blunt bodies

Hypersonic rarefied gas flow over blunt bodies in the transitional flow regime (from continuum to free-molecule) is investigated. Asymptotically correct boundary conditions on the body surface are derived for the full and thin viscous shock layer models. The effect of taking into account the slip velocity and the temperature jump in the boundary condition along the surface on the extension of the limits of applicability of continuum models to high free-stream Knudsen numbers is investigated. Analytic relations are obtained, by an asymptotic method, for the heat transfer coefficient, the skin friction coefficient and the pressure as functions of the free-stream parameters and the geometry of the body in the flow field at low Reynolds number; the values of these coefficients approach their values in free-molecule flow (for unit accommodation coefficient) as the Reynolds number approaches zero. Numerical solutions of the thin viscous shock layer and full viscous shock layer equations, both with the no-slip boundary conditions and with boundary conditions taking into account the effects slip on the surface are obtained by the implicit finite-difference marching method of high accuracy of approximation. The asymptotic and numerical solutions are compared with the results of calculations by the Direct Simulation Monte Carlo method for flow over bodies of different shape and for the free-stream conditions corresponding to altitudes of 75–150 km of the trajectory of the Space Shuttle, and also with the known solutions for the free-molecule flow regine. The areas of applicability of the thin and full viscous shock layer models for calculating the pressure, skin friction and heat transfer on blunt bodies, in the hypersonic gas flow are estimated for various free-stream Knudsen numbers.     

On the Rayleigh-Taylor Instability for Two Uniform Viscous Incompressible Flows

The authors study the Rayleigh-Taylor instability for two incompressible immis- cible fluids with or without surface tension, evolving with a free interface in the presence of a uniform gravitational field in Eulerian coordinates. To deal with the free surface, instead of using the transformation to Lagrangian coordinates, the perturbed equations in Eule- rian coordinates are transformed to an integral form and the two-fluid flow is formulated as a single-fluid flow in a fixed domain, thus offering an alternative approach to deal with the jump conditions at the free interface. First, the linearized problem around the steady state which describes a denser immiscible fluid lying above a light one with a free interface separating the two fluids, both fluids being in （unstable） equilibrium is analyzed. By a general method of studying a family of modes, the smooth （when restricted to each fluid domain） solutions to the linearized problem that grow exponentially fast in time in Sobolev spaces are constructed, thus leading to a global instability result for the linearized problem. Then, by using these pathological solutions, the global instability for the corresponding nonlinear problem in an appropriate sense is demonstrated.     

Hydrodynamic and thermal slip flow boundary layers over a flat plate with constant heat flux boundary condition

In this paper the boundary layer flow over a flat plat with slip flow and constant heat flux surface condition is studied. Because the plate surface temperature varies along the x direction, the momentum and energy equations are coupled due to the presence of the temperature gradient along the plate surface. This coupling, which is due to the presence of the thermal jump term in Maxwell slip condition, renders the momentum and energy equations non-similar. As a preliminary study, this paper ignores this coupling due to thermal jump condition so that the self-similar nature of the equations is preserved. Even this fundamental problem for the case of a constant heat flux boundary condition has remained unexplored in the literature. It was therefore chosen for study in this paper. For the hydrodynamic boundary layer, velocity and shear stress distributions are presented for a range of values of the parameter characterizing the slip flow. This slip parameter is a function of the local Reynolds number, the local Knudsen number, and the tangential momentum accommodation coefficient representing the fraction of the molecules reflected diffusively at the surface. As the slip parameter increases, the slip velocity increases and the wall shear stress decreases. These results confirm the conclusions reached in other recent studies. The energy equation is solved to determine the temperature distribution in the thermal boundary layer for a range of values for both the slip parameter as well as the fluid Prandtl number. The increase in Prandtl number and/or the slip parameter reduces the dimensionless surface temperature. The actual surface temperature at any location of x is a function of the local Knudsen number, the local Reynolds number, the momentum accommodation coefficient, Prandtl number, other flow properties, and the applied heat flux.     

Computational fluid dynamics modelling of gas jets impinging onto liquid pools

Gas jets impinging onto a gas–liquid interface of a liquid pool are studied using computational fluid dynamics modelling, which aims to obtain a better understanding of the behaviour of the gas jets used metallurgical engineering industry. The gas and liquid flows are modelled using the volume of fluid technique. The governing equations are formulated using the density and viscosity of the “gas–liquid mixture”, which are described in terms of the phase volume fraction. Reynolds averaging is applied to yield a set of Reynolds-averaged conservation equations for the mass and momentum, and the k–ε turbulence model. The deformation of the gas–liquid interface is modelled by the pressure jump across the interface via the Young–Laplace equation. The governing equations in the axisymmetric cylindrical coordinates are solved using the commercial CFD code, FLUENT. The computed results are compared with experimental and theoretical data reported in the literature. The CFD modelling allows the simultaneous evaluation of the gas flow field, the free liquid surface and the bulk liquid flow, and provides useful insight to the highly complex, and industrially significant flows in the jetting system.     

Solvability in Hölder spaces of a model initial-boundary value problem generated by a problem on the motion of two fluids

One studies the initial-boundary value problem for the Stokes' system, arising at the investigation of the nonstationary motion of two viscous fluids, separated by a free surface. Junction conditions are prescribed in the plane ×₃=0}. The consideration of the surface tension leads to a noncoercive integral term in the condition for the jump of the normal stresses. The unique solvability and estimates of the solution in Hölder classes of functions of the given model problem are proved with the aid of a theorem on Fourier multipliers and a significant part of the paper is devoted to the proof of the required modifications of this theorem.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova Akademii Nauk SSSR, Vol. 188, pp. 5–44, 1991.     

Electrohydrodynamic instability of two superposed viscous miscible streaming fluids

The linear stability of two dielectric viscous fluids separated by a horizontal interface is investigated. The interface admits heat and mass transfer. The system is stressed by a normal periodic electric field producing surface charges at the interface. The effect of surface tension, small viscosity, velocity streaming and gravity on the critical surface charges density and on corresponding electric field are analyzed. The contribution of viscosity with the existence of surface charges and streaming are discussed. The investigation includes the stability analysis of the presence of the periodic electric field as well as the constant one. It is found that the presence of the surface charges made by the normal electric field play a dual role in the stability criterion, which shows some analogy with the nonlinear theory of stability. Some previous studies are compared using appropriate data. The marginal state of stability is also considered. It is found that the surface charges vanish under certain conditions. This study shows that the mass and heat transfer parameter has a destabilizing effect whether the electric field is static or periodic. Parametric excitation of the electrohydrodynamic (EHD) surface waves is analyzed in the case of Rayleigh–Taylor (R–T) instability. The transition curves are obtained by means of Whittaker's technique. The analytical results are numerically confirmed.     

A comparative analysis of approaches for investigating hypersonic flow over blunt bodies in a transitional regime

Hypersonic flows of a viscous perfect rarefied gas over blunt bodies in a transitional flow regime from continuum to free molecular, characteristic when spacecraft re-enter Earth's atmosphere at altitudes above 90-100 km, are considered. The two-dimensional problem of hypersonic flow is investigated over a wide range of free stream Knudsen numbers using both continuum and kinetic approaches: by numerical and analytical solutions of the continuum equations, by numerical solution of the Boltzmann kinetic equation with a model collision integral in the form of the S-model, and also by the direct simulation Monte Carlo method. The continuum approach is based on the use of asymptotically correct models of a thin viscous shock layer and a viscous shock layer. A refinement of the condition for a temperature jump on the body surface is proposed for the viscous shock layer model. The continuum and kinetic solutions, and also the solutions obtained by the Monte Carlo method are compared. The effectiveness, range of application, advantages and disadvantages of the different approaches are estimated.     

Unsteady convective boundary layer flow of a viscous fluid at a vertical surface with variable fluid properties

In this paper we present numerical solutions to the unsteady convective boundary layer flow of a viscous fluid at a vertical stretching surface with variable transport properties and thermal radiation. Both assisting and opposing buoyant flow situations are considered. Using a similarity transformation, the governing time-dependent partial differential equations are first transformed into coupled, non-linear ordinary differential equations with variable coefficients. Numerical solutions to these equations subject to appropriate boundary conditions are obtained by a second order finite difference scheme known as the Keller-Box method. The numerical results thus obtained are analyzed for the effects of the pertinent parameters namely, the unsteady parameter, the free convection parameter, the suction/injection parameter, the Prandtl number, the thermal conductivity parameter and the thermal radiation parameter on the flow and heat transfer characteristics. It is worth mentioning that the momentum and thermal boundary layer thicknesses decrease with an increase in the unsteady parameter.     

The Steady Viscous Flow of Two Differentially Rotating Immiscible Fluids

We have computed the steady, axisymmetric viscous boundary layers on either side of an interface between two immiscible, incompressible fluids that are in rigid body rotation far from the interface. The internal rotational Froude number is assumed small so that the interface may be considered horizontal. An application of our results to the spinup from rest of two immiscible slightly viscous fluids in a vertically mounted cylinder is discussed.     

The balance of material momentum at a shock wave propagating in a thermo-elastic material

At a surface of discontinuity, the mechanical balance laws are represented by jump conditions. It is shown how the balance of material momentum at an adiabatic shock propagating in a thermo-elastic material is obtained from the discontinuous balances of physical momentum and energy. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Global existence of two phase nonhomogeneous viscous incompbessible fluid flow

We shall discuss the temporarily global solution for the two phase free boundary problem. Both fluids are regarded as immiscible, nonhomogeneous, viscous and incompressible and subject to surface tention on the interface. The global solution is obtained near the equilibrium state under the sufficiently small initial data and external forces.     

Nonlinear analysis of capillary instability with heat and mass transfer

The nonlinear capillary instability of the cylindrical interface between the vapor and liquid phases of a fluid is studied when there is heat and mass transfer across the interface, using viscous potential flow theory. The fluids are considered to be viscous and incompressible with different kinematic viscosities. Both asymmetric and axisymmetric disturbances are considered. The analysis is based on the method of multiple scale perturbation and the nonlinear stability is governed by first-order nonlinear partial differential equation. The stability conditions are obtained and discussed theoretically as well as numerically. Regions of stability and instability have been shown graphically indicating the effect of various parameters. It has been observed that the heat and mass transfer has stabilizing effect on the stability of the system in the nonlinear analysis for both axisymmetric as well as asymmetric disturbances.     

Kinetic Equation for Soliton Gas: Integrable Reductions

We study the dynamic behaviour of two viscous fluid films confined between two concentric cylinders rotating at a small relative velocity. It is assumed that the fluids are immiscible and that the volume of the outer fluid film is large compared to the volume of the inner one. Moreover, while the outer fluid is considered to have constant viscosity, the rheological behaviour of the inner thin film is determined by a strain-dependent power-law. Starting from a Navier–Stokes system, we formally derive evolution equations for the interface separating the two fluids. Two competing effects drive the dynamics of the interface, namely the surface tension and the shear stresses induced by the rotation of the cylinders. When the two effects are comparable, the solutions behave, for large times, as in the Newtonian regime. We also study the regime in which the surface tension effects dominate the stresses induced by the rotation of the cylinders. In this case, we prove local existence of positive weak solutions both for shear-thinning and shear-thickening fluids. In the latter case, we show that interfaces which are initially close to a circle converge to a circle in finite time and keep that shape for later times.

     

Simulations of swelling media by weakly fulfilled Dirichlet boundary conditions

Charged hydrated porous media, which are found in biomechanics as well as in geomechanics, have the capability to change their volume under varying chemical conditions of the environment. In this contribution, these materials are modelled in the framework of the thermodynamically consistent Theory of Porous Media (TPM). The underlying model consists of four constituents, a charged solid and an aqueous solution composed of water and the ions of dissolved salt. The solid is modelled by a finite elasticity law accounting for the multiphasic micro structure, whereas the fluid is considered as a viscous Newton ian fluid. One finally ends up with four balance relations, the volume balance of the fluid, the concentration balance of the cations, the momentum balance and the balance of charges of the overall mixture. The resulting set of partial differential equations is solved within the framework of the FEM. Therefore, the weak forms are derived and the resulting set of equations for the primary variables pore pressure p, cation concentration c_m and solid displacement u _S is implemented into the FE tool PANDAS. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)