Example of an exact solution of the stationary problem of two-layer flows with evaporation at the interface

Stationary two-layer liquid and gas flows with fluid evaporation at the interface are studied. On the solid impermeable boundaries of the channel, no-slip conditions are satisfied and a linear temperature distribution along the longitudinal coordinate and a condition for the vapor concentration at the upper boundary are specified. On the thermocapillary interface, remaining undeformed, the following conditions are specified: kinematic and dynamic conditions, a condition for thermal flows with mass transfer, continuity conditions for the velocity, temperature, and mass balance, and a relation for the saturated vapor concentration. An exact solution of the stationary problem for a given gas flow rate is obtained. Examples of velocity profiles are given for stationary flows of the ethanol-nitrogen system under normal and reduced gravity are given. The effect of longitudinal temperature gradients specified at the boundaries of the channel on the flow pattern is investigated.     

Stability of two-layer fluid flows with evaporation at the interface

The problem of stability of two-layer (fluid-gas) flows with account of evaporation at the thermocapillary interface is studied under the condition of a fixed gas flow rate. In the upper gas-vapor layer, the Dufour effect is taken into account. A novel exact solution of the Navier–Stokes equations in the Boussinesq approximation is constructed. The effects of longitudinal temperature gradients, gravity, thicknesses of the gas and fluid layers, and the gas flow rate on the flow structure, the onset of recirculated flows near the interface, the evaporation rate, and the properties of characteristic disturbances are investigated.     

Convective flows in a cylindrical fluid zone in a high-frequency vibration field

Convective flows of a nonuniformly heated fluid in a cylindrical fluid zone in a high-frequency longitudinal vibration field are studied. Vibration frequencies which are high as compared with dissipative decrements and capillary frequencies, but small as compared with acoustic frequencies are considered. The general method formulated earlier for describing the behavior of inhomogeneous fluids under the influence of high-frequency vibrations is used. The interaction between the vibrational flow mechanisms and thermocapillary effects on a free surface is analyzed.Perm'. Marseilles (France). Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, pp. 53–61, September–October, 1994.     

Velocity-profile deformation in EHD flow in a two-dimensional channel

The study of unipolar-charged fluids in the presence of external and induced electric fields has recently taken on great importance. The characteristics of one-dimensional EGD flows [1, 2] and developed laminar flows of a viscous fluid [3] have been clarified in several studies made in this field. However, the study of three-dimensional flows of such media is actually just beginning. Here, along with the analysis of three-dimensional boundary layers and jets [4], there is considerable interest in the study of spatial (two-dimensional and three-dimensional) EHD flows of an inviscid fluid, since in many engineering devices the zone of interaction of the flow with the electric fields does not exceed a few channel diameters, which makes it possible to neglect viscous effects.In this paper we examine some aspects of two-dimensional EHD flows of a viscous incompressible medium for infinitely large electric Reynolds numbers. The perturbations of the hydrodynamic parameters of the flow downstream from the zone of action of the electrostatic forces are determined. It is shown that in many cases the flow parameters outside this zone may be determined without solving the complete system of EHD partial differential equations.     

Flow in rotating radial channels at small Rossby and Ekman numbers

The problem of flow of a viscous fluid in a rotating channel is considered in the region of very small Rossby and Ekman numbers and moderately large Reynolds numbers. Asymptotic expressions with respect to the Ekman number are found for the velocity components and the longitudinal pressure gradient by solving a system of linear differential equations using Fourier series. The stability limits of such flow are predicted. Attention is drawn to a similarity between the velocity profiles of these flows and flows of a magnetic fluid and a fluid executing longitudinal oscillations in a fixed channel.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 11–15, January–February, 1984.     

A numerical study of the longitudinal thermoconvective rolls in a mixed convection flow in a horizontal channel with a free surface

This paper presents a numerical study of three-dimensional laminar mixed convection within a liquid flowing on a horizontal channel heated uniformly from below. The upper surface is free and assumed to be flat. The coupled Navier–Stokes and energy equations are solved numerically by the finite volume method taking into account the thermocapillary effects (Marangoni effect). When the strength of the buoyancy, thermocapillary effects and forced convective currents are comparable (Ri ⋍ O(1) and Bd = Ra/Ma ⋍ O(1)), the results show that the development of instabilities in the form of steady longitudinal convective rolls is similar to those encountered in the Poiseuille–Rayleigh–Bénard flow. The number and spatial distribution of these rolls along the channel depend on the flow conditions. The objective of this work is to study the influence of parameters, such as the Reynolds, Rayleigh and Biot numbers, on the flow patterns and heat transfer characteristics. The effects of variations in the surface tension with temperature gradients (Marangoni effect) are also considered.     

Control of Thermo- and Solutocapillary Flows in FZ Crystal Growth by High-Frequency Vibrations

The paper deals with the numerical investigation of the possibilities to control convective flows in the liquid bridge in zero gravity conditions applying axial vibrations. The surface tension is assumed to be dependent both on the temperature and on the solute concentration. The free surface deformations and the curvature of the phase change surfaces are neglected but pulsational deformations of the free surface are accounted for. The first part of the paper concerns axisymmetric steady flows. The calculations show that the evolution of convective flow with the variation of thermal Marangoni number at a fixed value of the solutal Marangoni number is accompanied by the hysteresis phenomenon, which is related to the existence of two stable steady regimes in a certain parameter range. One of these regimes is thermocapillary dominated, it corresponds to the two-vortex flow, and the other is solutocapillary dominated, it corresponds to the single-vortex flow. Under vibrations, the range of the Marangoni numbers where hysteresis is observed becomes narrower and is shifted to the area of larger values. The second part of the paper concerns the stability of axisymmetric thermo-and solutocapillary flows and the transition to three-dimensional regimes. Significant mutual influence of flows generated by each process on the stability of the other is discovered. Stability maps in the parametric plane for the thermal Marangoni number, the solutal Marangoni number, are obtained for different values of vibration parameters. It is shown, that vibrations exert a stabilizing effect, increasing critical Marangoni numbers for all modes of instability. However, this effect is different for different modes and at high vibration intensity destabilization is possible. Consequently, vibrations can modify the scenario of the transition to the three-dimensional mode.     

A finite element solutionof three-dimensional mixed convection gas flows in horizontal chnnels using preconditioned iterative metrix methods

An algorithm is presented for the finite element solution of three-dimensional mixed convection gas flows in channels heated from below. The algorithm uses Newton's method and iterative matrix methods. Two iterative solution algorithms, conjugate gradient squared (CGS) and generalized minimal residual (GMERS), are used in conjunction with a preconditioning technique that is simple to implement. The preconditioner is a subset of the full Jacobian matrix centered around the main diagonal but retaining the most fundamental axial coupling of the residual equations. A domain-renumbering scheme that enhances the overall algorithm performance is proposed on the basis of and analysis of the preconditioner. Comparison with the frontal elimination method demonstrates that the iterative method will be faster when the front width exceeds approximately 500. Techniques for the direct assembly f the problem into a compressed sparse row storage format are demonstrated. Elimination of fixed boundary conditions is shown to decrease the size of the matrix problem by up to 30%. Finally, fluid flow solutions obtained with the numerical technique are presented. These solutions reveal complex three-dimensional mixed convection fluid flow phenomena at low Reynolds numbers, including the reversal of the direction of longitudinal rolls in the presence of a strong recirculation in the entrance region of the channel.     

Selection of convective rolls in a thin evaporating-fluid layer in an air flow

The process of selection of longitudinal convective rolls in a thin layer of evaporating fluid immersed in an air turbulent boundary layer flow is studied numerically. The dependence of the two-dimensional flow patterns on the Rayleigh number and boundary conditions is analyzed. Calculations with account for the thermocapillary effect are carried out. The numerical results are compared with experimental data.     

Flow of a liquid with an anisotropic viscosity tensor

This paper presents a simple model of an anisotropic incompressible viscous fluid, whose equations of state involve one anisotropic physical constant tensor (in the sense of Oldroyd [1]). Attention is restricted to the case of a fluid that is everywhere transversely isotropic at some given instant, so that the model is essentially one of a liquid with initially just one privileged direction at each point. Transverse isotropy does not persist, however, in some flow situations.Predictions are made for simple shearing flows and for channel and pipe flows with different initial directions of orientation. In some cases, the volume rate of flow under constant pressure-gradient decreases steadily and tends to zero after a long time.     

Linear stability of two-dimensional steady flow in wavy-walled channels

Linear stability of fully developed two-dimensional periodic steady flows in sinusoidal wavy-walled channels is investigated numerically. Two types of channels are considered: the geometry of wavy walls is identical and the location of the crest of the lower and upper walls coincides (symmetric channel) or the crest of the lower wall corresponds to the furrow of the upper wall (sinuous channel). It is found that the critical Reynolds number is substantially lower than that for plane channel flow and that when the non-dimensionalized wall variation amplitude is smaller than a critical value (about 0.26 for symmetric channel, 0.28 for sinuous channel), critical modes are three-dimensional stationary and for larger , two-dimensional oscillatory instabilities set in. Critical Reynolds numbers of sinuous channel flows are smaller for three-dimensional disturbances and larger for two-dimensional disturbances than those of symmetric channel flows. The disturbance velocity distribution obtained by the linear stability analysis suggests that the three-dimensional stationary instability is mainly caused by local concavity of basic flows near the reattachment point, while the critical two-dimensional mode resembles closely the Tollmien–Schlichting wave for plane Poiseuille flow.     

Thermocapillary convection stability in a cylindrically-symmetric two-phase system

Thermocapillary flows in an infinitely long liquid cylinder surrounded by a coaxial gas layer with a controlled flow rate and the stability of such flows are investigated. In the layers a constant axial temperature gradient is maintained. An exact solution of the equations of motion describing the steady-state flow in this two-phase system is derived. Possible flow regimes and their stability in the linear approximation are studied. It is shown that in the liquid phase the thermocapillary flow can be completely stopped by the gas flow at the expense of the interaction between mechanical stresses at the interface. The results obtained indicate the possibility of controlling thermocapillary flows and their stability by means of gas flows.     

Influence of rotating magnetic field strength on three-dimensional thermocapillary flow in a floating half-zone model

The effects of rotating magnetic field (RMF) on the three-dimensional thermocapillary flow of semiconductor melt (Pr?=?0.01) in a floating half-zone model under microgravity are investigated numerically by the finite volume method. The results indicate that the thermocapillary flow without magnetic field is a steady three-dimensional convection for Ma?=?40 in a floating half-zone model with As?=?1, and the convection evolves to an oscillatory three-dimensional flow by applying 1–6?mT RMF with 50?Hz rotating frequency. Based on the fast Fourier transform spectrum, the convection is confirmed to be a periodically oscillating flow, the oscillatory main frequency, 1.59?×?10^?3?Hz for 1?mT RMF and 5.84?×?10^?2?Hz for 6?mT RMF, increases with the magnetic strength. However, with increasing the magnetic field strength up to 7?mT, the three-dimensional thermocapillary flow is effectively controlled and the convection turns into a steady axisymmetrical one.     

Validation of the S‐CLSVOF method with the density‐scaled balanced continuum surface force model in multiphase systems coupled with thermocapillary flows

Three numerical methods, namely, volume of fluid (VOF), simple coupled volume of fluid with level set (S‐CLSVOF), and S‐CLSVOF with the density‐scaled balanced continuum surface force (CSF) model, have been incorporated into OpenFOAM source code and were validated for their accuracy for three cases: (i) an isothermal static case, (ii) isothermal dynamic cases, and (iii) non‐isothermal dynamic cases with thermocapillary flow including dynamic interface deformation. Results have shown that the S‐CLSVOF method gives accurate results in the test cases with mild computation conditions, and the S‐CLSVOF technique with the density‐scaled balanced CSF model leads to accurate results in the cases of large interface deformations and large density and viscosity ratios. These show that these high accuracy methods would be appropriate to obtain accurate predictions in multiphase flow systems with thermocapillary flows. Copyright © 2016 John Wiley & Sons, Ltd.     

Fluid flow through heterogeneous porous media in a rotating square channel

An analytical three-dimensional solution to the fluid flow problem through heterogeneous porous media in a rotating square channel is presented. The permeability of the fluid saturated porous domain varies in the vertical direction, thus affecting the imposed main flow in the channel. As a result of Coriolis acceleration, secondary circulation in a plane perpendicular to the main flow direction is created. A particular example of a monotonic distribution of the permeability function is analyzed leading to a single vortex secondary circulation. Nevertheless, multiple vortex secondary flow solutions are possible depending on the particular variation of the permeability in the vertical direction. No secondary motion is expected for isothermal flows in homogeneous porous media.     

Turbulence manipulation in air–water flows on a stepped chute: An experimental study

In a stepped channel operating with large flow rates, the flow skims over the pseudo-bottom formed by the step edges as a coherent stream. Intense three-dimensional recirculation is maintained by shear stress transmission from the mainstream to the step cavities, while significant free-surface aeration takes place. The interactions between free-surface aeration and cavity recirculation are investigated herein with seven step cavity configurations. The experiments were conducted in a large stepped channel operating at large Reynolds numbers. For some experiments, triangular vanes, or longitudinal ribs, were placed across the step cavities to manipulate the flow turbulence to enhance the interactions between the mainstream flow and the cavity recirculation region. The results showed a strong influence of the vanes on the air–water flow properties in both free-stream and cavity flows. The findings demonstrate some passive turbulence manipulation in highly turbulent air–water flows.     

Some self-similar viscous gas flows in a channel

We consider the class of self-similar axisymmetric and two-dimensional laminar flows of a viscous gas in a long channel with smooth contour, in which the longitudinal component of the velocity and the gas temperature are functions of a single dimensionless transverse coordinate. Such flows correspond to exponential (axisymmetric flow) or linear (two-dimensional flow) increase of the radius or height of the channel and corresponding exponential or hyperbolic decrease of the static pressure along the channel.     

Conical and Quasi-Conical Incompressible Fluid Flows

The solution of the problem of fluid flow inside a cone with a small vertex angle is obtained in closed form. The conditions of occurrence of singular separation are considered within the framework of conical flow theory. A class of conical flows in which the vorticity is transported along streamlines of the potential velocity component is detected.Quasi-conical incompressible fluid flow, i.~e. a flow inside and outside an axisymmetric body with power-law generators is defined by analogy with supersonic compressible fluid flow. The conditions under which the effect of vorticity and swirling is significant are found as a result of an inspection analysis. An approximate solution of the problem of fluid flow inside a zero corner is found.A coordinate expansion representing a plane analog of conical flow is constructed in the neighborhood of the separation point of a creeping flow on a smooth surface.     

Computation of Stokes Flow in a Channel with a Collapsible Segment

The numerical solution of Stokes flow in two-dimensional channel in which a segment of one wall is formed by an elastic membrane under longitudinal tension and the remaining channel boundary is rigid is considered. This model problem is being used to gain an understanding of the complex interactions that occurs between the fluid flow and the wall mechanics when fluid flows through a collapsible tube, examples of which are widespread in physiology. Previous work by Pedley considered a similar system using lubrication theory in which the wall slopes are assumed small. The results showed that as the longitudinal wall tension is reduce, the downstream end of the collapsible segment becomes ever steeper, thus violating the assumptions. Here, lubrication theory is abandoned and a numerical solution of the full governing equations, including the complete expression for wall curvature, is sought using an iterative scheme. The effect of the variation in wall tension due to the fluid shear stresses at the compliant boundary is also included.Results are presented for a range of transmural (internal minus external) pressures and wall tensions. It is found, however, that as the wall tension is reduced, the iterative scheme considered fails to converge. This similar behaviour to that seen by Silliman & Scriven in viscous free-surface flows. Possible reasons for this breakdown together with alternative solution strategies are discussed.     

Gas-kinetic numerical method for solving mesoscopic velocity distribution function equation

A gas-kinetic numerical method for directly solving the mesoscopic velocity distribution function equation is presented and applied to the study of three-dimensional complex flows and micro-channel flows covering various flow regimes. The unified velocity distribution function equation describing gas transport phenomena from rarefied transition to continuum flow regimes can be presented on the basis of the kinetic Boltzmann–Shakhov model equation. The gas-kinetic finite-difference schemes for the velocity distribution function are constructed by developing a discrete velocity ordinate method of gas kinetic theory and an unsteady time-splitting technique from computational fluid dynamics. Gas-kinetic boundary conditions and numerical modeling can be established by directly manipulating on the mesoscopic velocity distribution function. A new Gauss-type discrete velocity numerical integration method can be developed and adopted to attack complex flows with different Mach numbers. HPF parallel strategy suitable for the gas-kinetic numerical method is investigated and adopted to solve three-dimensional complex problems. High Mach number flows around three-dimensional bodies are computed preliminarily with massive scale parallel. It is noteworthy and of practical importance that the HPF parallel algorithm for solving three-dimensional complex problems can be effectively developed to cover various flow regimes. On the other hand, the gas-kinetic numerical method is extended and used to study micro-channel gas flows including the classical Couette flow, the Poiseuille- channel flow and pressure-driven gas flows in two-dimensional short micro-channels. The numerical experience shows that the gas-kinetic algorithm may be a powerful tool in the numerical simulation of micro-scale gas flows occuring in the Micro-Electro-Mechanical System (MEMS). The project supported by the National Natural Science Foundation of China (90205009 and 10321002), and the National Parallel Computing Center in Beijing. The English text was polished by Yunming Chen.