Marangoni effect in the Couette flow

Viscous heating in an axisymmetric creeping flow of a second-order fluid with free surface between two coaxially mounted cylinders produces a radial temperature gradient in the fluid. The dependence of the surface tension upon temperature is the cause for a secondary flow in the meridional plane of the flow field. This secondary flow (Marangoni effect), and its influence upon the shape of the free surface are studied. The deformation of the free surface caused by the Marangoni effect is compared with the deformation caused by inertia and normal stress differences.     

Emergence of Fluid Rotation in the Marangoni Boundary Layer in the Region of Local Cooling of the Free Boundary

Rotation bifurcation in a steady axisymmetric thermocapillary flow of an incompressible fluid filling a semi-infinite space bounded by the free surface with a nonuniform distribution of temperature is studied. The fluid flow is calculated on the basis of Navier–Stokes equations under the assumption of small diffusion coefficients. It is shown that the bifurcation triggers rotational motion in a thin Marangoni boundary layer in the case of local cooling of the free boundary near the axis of symmetry and in the presence of an external flow; there is no rotation outside this layer. In the case of local heating of the free boundary, rotation is not observed.     

Asymptotic analysis of the finite cone-and-plate flow of a non-newtonian fluid

Consider the shearing flow of a viscoelastic fluid trapped by surface tension between a cone and a plate. An asymptotic analysis of this problem in the limit of small gap angle has been done. This limit is realized in many practical situations. It is assumed that the Deborah number De, the Reynolds number Re, and the retardation parameter β are all order unity and that the shape of the free surface is very nearly spherical. Closed form analytic expressions are obtained for the leading terms of the primary and weak secondary motion of the fluid as well as the meniscus shape. It is found that the velocity field is bounded and continuous if and only if . There is a family of curves in the De-β plane on which the velocity field has a removable singularity at the origin. The secondary flow is made up of either one or two toroidal vortices. The meniscus has a bulge near the rotating cone and a trough near the stationary plate.     

Marangoni convection in weld pools with a free surface

A steady-state two-dimensional model of heat transfer and fluid flow was developed to describe Marangoni convection in the weld pool. Both the pool surface and the fusion boundary were calculated. The validity of the model was verified against an asymptotic solution for Marangoni-convection-induced free surface geometry. Two different cases were studied, i.e. a negative surface tension temperature coefficient ?γ/?T a positive one, and the resultant shapes of the weld pool surface were compared.     

Natural convection with evaporation in a vertical cylindrical cavity under the effect of temperature-dependent surface tension

The effect of surface tension on laminar natural convection in a vertical cylindrical cavity filled with a weak evaporating liquid has been analyzed numerically. The cylindrical enclosure is insulated at the bottom, heated by a constant heat flux from the side, and cooled by a non-uniform evaporative heat flux from the top free surface having temperature-dependent surface tension. Governing equations with corresponding boundary conditions formulated in dimensionless stream function, vorticity, and temperature have been solved by finite difference method of the second-order accuracy. The influence of Rayleigh number, Marangoni number, and aspect ratio on the liquid flow and heat transfer has been studied. Obtained results have revealed that the heat transfer rate at free surface decreases with Marangoni number and increases with Rayleigh number, while the average temperature inside the cavity has an opposite behavior; namely, it growths with Marangoni number and reduces with Rayleigh number.     

Effect of magnetic field on the buoyancy and thermocapillary driven convection of an electrically conducting fluid in an annular enclosure

The main objective of this article is to study the effect of magnetic field on the combined buoyancy and surface tension driven convection in a cylindrical annular enclosure. In this study, the top surface of the annulus is assumed to be free, and the bottom wall is insulated, whereas the inner and the outer cylindrical walls are kept at hot and cold temperatures respectively. The governing equations of the flow system are numerically solved using an implicit finite difference technique. The numerical results for various governing parameters of the problem are discussed in terms of the streamlines, isotherms, Nusselt number and velocity profiles in the annuli. Our results reveal that, in tall cavities, the axial magnetic field suppresses the surface tension flow more effectively than the radial magnetic field, whereas, the radial magnetic field is found to be better for suppressing the buoyancy driven flow compared to axial magnetic field. However, the axial magnetic field is found to be effective in suppressing both the flows in shallow cavities. From the results, we also found that the surface tension effect is predominant in shallow cavities compared to the square and tall annulus. Further, the heat transfer rate increases with radii ratio, but decreases with the Hartmann number.     

AN IMPLICIT SURFACE TENSION ALGORITHM FOR PICARD SOLVERS OF SURFACE–TENSION–DOMINATED FREE AND MOVING BOUNDARY PROBLEMS

One of the methods for solving a free or moving boundary problem is the use of Picard solvers which solve the geometry and the velocity field successively. When, however, the kinematic condition is used for updating the geometry in this technique, numerical stability problems occur for surface-tension-dominated flow. These problems are shown here to originate from the unstable integration of the local smoothing of the surface by surface tension. By an extension of the surface tension contribution to the flow field an implicit treatment of surface tension is obtained which overcomes these stability problems. The algorithm is applicable to both free and moving boundary problems, as will be shown by examples in this paper.     

双自由面溶质?热毛细液层的不稳定性

溶质?热毛细对流是流体界面的浓度和温度分布不均导致的表面张力梯度驱动的流动, 它主要存在于空间微重力环境、小尺度流动等表面张力占主导的情况中, 例如晶体生长、微流控、合金浇筑凝固、有机薄液膜生长等. 对其流动进行稳定性分析具有重要意义. 本文采用线性稳定性理论研究了双自由面溶质?热毛细液层对流的不稳定性, 得到了两种负毛细力比(η)下的临界Marangoni数与Prandtl数(Pr)的函数关系, 并分析了临界模态的流场和能量机制. 研究发现: 溶质?热毛细对流和纯热毛细对流的临界模态有较大的差别, 前者是同向流向波、逆向流向波、展向稳态模态和逆向斜波, 后者是逆向斜波和逆向流向波. 在Pr较大时, Pr增加会降低流动稳定性; 在其他参数下, Pr增加会增强流动稳定性. 在中低Pr, 溶质毛细力使流动更加不稳定; 在大Pr时, 溶质毛细力的出现可能使流动更加稳定; 在其他参数下, 溶质毛细力会减弱流动稳定性. 流动稳定性不随η单调变化. 在多数情况下, 扰动浓度场与扰动温度场都是相似的. 能量分析表明: 扰动动能的主要能量来源是表面张力做功, 但其中溶质毛细力和热毛细力做功的正负性与参数有关.      

Forced Convection Heat Transfer from a Finite-Height Cylinder

This paper presents a large eddy simulation of forced convection heat transfer in the flow around a surface-mounted finite-height circular cylinder. The study was carried out for a cylinder with height-to-diameter ratio of 2.5, a Reynolds number based on the cylinder diameter of 44 000 and a Prandtl number of 1. Only the surface of the cylinder is heated while the bottom wall and the inflow are kept at a lower fixed temperature. The approach flow boundary layer had a thickness of about 10% of the cylinder height. Local and averaged heat transfer coefficients are presented. The heat transfer coefficient is strongly affected by the free-end of the cylinder. As a result of the flow over the top being downwashed behind the cylinder, a vortex-shedding process does not occur in the upper part, leading to a lower value of the local heat transfer coefficient in that region. In the lower region, vortex-shedding takes place leading to higher values of the local heat transfer coefficient. The circumferentially averaged heat transfer coefficient is 20 % higher near the ground than near the top of the cylinder. The spreading and dilution of the mean temperature field in the wake of the cylinder are also discussed.     

Unraveling surfactant transport on a thin liquid film

When a drop of insoluble surfactant is deposited on the surface of a thin liquid film, a radial flow is induced by the resulting surface tension gradient. It is difficult in practice to measure or visualize the evolution of the surfactant concentration and the corresponding surface tension field. In this contribution, we propose a numerical technique which allows, in theory, the reconstruction of the surfactant concentration and surface tension fields from the knowledge of the free surface velocity. The method also requires the knowledge of the equation of state relating the surfactant concentration to the surface tension. The proposed method is based on a reformulation of the lubrication approximation which then takes as an input the free surface velocity field. As a by-product, the film thickness is also reconstructed. We also show in this contribution, that the surface diffusion coefficient can also be estimated, in principle. The methodologies are successfully tested on ideal, synthetic data-sets but also on under-resolved, noisy, data-sets more representative of true experimental conditions. This contribution may help shed some light on the phenomena involved in surfactant transport.     

Surface Waves for a Compressible Viscous Fluid

In this paper, we are concerned with free boundary problem for compressible viscous isotropic Newtonian fluid. Our problem is to find the three-dimensional domain occupied by the fluid which is bounded below by the fixed bottom and above by the free surface together with the density, the velocity vector field and the absolute temperature of the fluid satisfying the system of Navier-Stokes equations and the initial-boundary conditions. The Navier-Stokes equations consist of the conservations of mass, momentum under the gravitational field in a downward direction and energy. The effect of the surface tension on the free surface is taken into account. The purpose of this paper is to establish two existence theorems to the problem mentioned above: the first concerns with the temporary local solvability in anisotropic Sobolev-Slobodetskiĭ spaces and the second the global solvability near the equilibrium rest state. Here the equilibrium rest state (heat conductive state) means that the temperature distribution is a linear function with respect to a vertical direction and the density is determined by an ordinary differential equation which involves equation of state. For the proof, we rely on the methods due to Solonnikov in the case of incompressible fluid with some modifications, since our problem is hyperbolic-parabolic coupled system. Dedicated to Professors Takaaki Nishida and Masayasu Mimura on their sixtieth birthdays     

Numerical solutions of free convection boundary layer flow on a solid sphere with Newtonian heating in a micropolar fluid

In this paper, the problem of free convection boundary layer flow on a solid sphere in a micropolar fluid with Newtonian heating, in which the heat transfer from the surface is proportional to the local surface temperature, is considered. The transformed boundary layer equations in the form of partial differential equations are solved numerically using an implicit finite-difference scheme. Numerical solutions are obtained for the local wall temperature, the local skin friction coefficient, as well as the velocity, angular velocity and temperature profiles. The features of the flow and heat transfer characteristics for different values of the material or micropolar parameter K, the Prandtl number Pr and the conjugate parameter γ are analyzed and discussed.     

Similarity solutions for MHD thermosolutal Marangoni convection over a flat surface in the presence of heat generation or absorption effects

The problem of steady, laminar, thermosolutal Marangoni convection flow of an electrically-conducting fluid along a vertical permeable surface in the presence of a magnetic field, heat generation or absorption and a first-order chemical reaction effects is studied numerically. The general governing partial differential equations are converted into a set of self-similar equations using unique similarity transformations. Numerical solution of the similarity equations is performed using an implicit, iterative, tri-diagonal finite-difference method. Comparisons with previously published work is performed and the results are found to be in excellent agreement. Approximate analytical results for the temperature and concentration profiles as well as the local Nusselt and sherwood numbers are obtained for the conditions of small and large Prandtl and Schmidt numbers are obtained and favorably compared with the numerical solutions. The effects of Hartmann number, heat generation or absorption coefficient, the suction or injection parameter, the thermo-solutal surface tension ratio and the chemical reaction coefficient on the velocity, temperature and concentration profiles as well as quantitites related to the wall velocity, boundary-layer mass flow rate and the Nusselt and Sherwood numbers are presented in graphical and tabular form and discussed. It is found that a first-order chemical reaction increases all of the wall velocity, Nusselt and Sherwood numbers while it decreases the mass flow rate in the boundary layer. Also, as the thermo-solutal surface tension ratio is increased, all of the wall velocity, boundary-layer mass flow rate and the Nusselt and Sherwood numbers are predicted to increase. However, the exact opposite behavior is predicted as the magnetic field strength is increased.     

Bingham流体双自由面热毛细液层的稳定性分析

热毛细对流是流体界面温度分布不均导致的表面张力梯度驱动的流动.它主要存在于空间等微重力环境或小尺度流动等表面张力占主导的情况中.在很多工业领域,如晶体生长、聚合物加工、喷墨打印、微流控,产品质量都与热毛细对流密切相关.空间3D打印是太空制造的重要技术,可以支持空间站的在轨长期有人照料的运行和维护,实现按需制造.本文以聚合物流体的空间3D打印为应用背景,采用线性稳定性理论研究了Bingham流体双自由面热毛细液层的稳定性,得到了在不同Bingham数(B)下的临界Marangoni数(Mac)与Prandtl数(Pr)的函数关系,分析了临界模态的流场和能量机制.研究发现:该流动的临界模态包括流向波和斜波模态,与B, Bi和两界面垂直方向上的温差(Q)相关. B和Bi的增加会增强热毛细对流的稳定性.当Q=0时,扰动温度分布分成对称和反对称两种情况.当Q> 0时, Pr的增加会减弱流动稳定性.在小Pr情况下,扰动温度分布在整个流场,在大Pr情况下,扰动温度在栓塞区为零.能量分析表明:扰动动能的主要能量来源是表面张力做功,但小Pr数下基本流也有一定贡献.     

微重力环境下Marangoni对流的有限元数值模拟

本文以悬浮区为背景研究液桥中气液交界面上由表面张力所驱动的流体对流。我们采用有限元方法对轴对称、定常运动方程,能量方程,自由面切向、法向应力平衡条件迭代求解,首次考虑了边界形状的影响,得到了自洽的流区位形和流场、温度场、表面压强分布。一般说来,流区自由面呈弯月形。结果表明,只要表面张力数 C_a<1,静态平衡的流区位形就是动态情况的良好近似。本文还分析了C_a数及G_r 数对流区位形的影响,得到了不同 M_a 数及散热条件下的温度场和流函数分布。     

自由面对潜艇尾流场流动特性影响研究

潜艇周围绕流场流动特性会影响潜艇的机动性能, 特别是近水面航行时, 自由面的存在会增大潜艇尾流场的复杂程度.为探究潜艇在近水面航行时自由液面对潜艇尾流场流动特性的影响机理,借助大型水下三维粒子图像测速技术开展潜艇尾流场流动特性研究.首先通过美国泰勒水池标准模型实验结果对试验方法准确性进行验证; 随后,用验证后的模型试验方法对潜艇尾流场进行测量,得到不同潜深工况、不同速度下的桨盘面轴向速度以及脉动速度,同时辅以数值模拟对试验无法测得的兴波波系及中纵剖面速度场加以补充,从兴波角度阐述了自由面对潜艇尾流场流动特性影响机理. 研究结果表明:潜艇在近水面航行时, 随着Fr增大,桨盘面处轴向速度云图中上方等值线整体趋于扁平化, 较4D潜深工况,1.5D潜深工况出现局部脉动速度极大值, 且脉动速度结构整体下移; 自由面存在时,艇体与自由面间流场速度明显增大, 特别在桨盘面区域, 流场速度明显提升.随着Fr增大, 桨盘面处的自由液面高度逐渐降低,这就导致了桨盘面位置出现更大的流体速度, 即造成了桨盘面伴流场挤压现象.      

Axisymmetric thermocapillary convection in a slowly rotating annular cylindrical container

In a slowly rotating annular cylindrical container the free liquid surface (liquid-gas interface) is subjected to a temperature gradient in radial direction. The temperature dependent surface tension creates a shear stress on the interface which is transmitting a thermocapillary convection in the bulk of the liquid. For constant temperature T ₁ of the inner and T ₂ of the outer wall a steady Marangoni convection takes place, exhibiting a double vortex ring of equal directional flow. For time-oscillatory temperatures of the walls a time-dependent thermocapillary convection appears, which will create on the free liquid surface various wave patterns. They shall, depending on the forcing frequency of the temperature, exhibit resonance peaks. The velocity distribution and the response magnitude inside the container has been determined. Received on 3 September 1997     

SEMI-IMPLICIT FINITE VOLUME SHALLOW-WATER FLOW AND SOLUTE TRANSPORT SOLVER WITH k–ε TURBULENCE MODEL

A 3D semi-implicit finite volume scheme for shallow- water flow with the hydrostatic pressure assumption has been developed using the σ-co-ordinate system, incorporating a standard k–ε turbulence transport model and variable density solute transport with the Boussinesq approximation for the resulting horizontal pressure gradients. The mesh spacing in the vertical direction varies parabolically to give fine resolution near the bed and free surface to resolve high gradients of velocity, k and ε. In this study, wall functions are used at the bed (defined by the bed roughness) and wind stress at the surface is not considered. Surface elevation gradient terms and vertical diffusion terms are handled implicitly and horizontal diffusion and source terms explicitly, including the Boussinesq pressure gradient term due to the horizontal density gradient. The advection terms are handled in explicit (conservative) form using linear upwind interpolation giving second-order accuracy. A fully coupled solution for the flow field is obtained by substi- tuting for velocity in the depth-integrated continuity equation and solving for surface elevation using a conjugate gradient equation solver. Evaluation of horizontal gradients in the σ-co-ordinate system requires high-order derivatives which can cause spurious flows and this is avoided by obtaining these gradients in real space. In this paper the method is applied to parallel oscillatory (tidal) flow in deep and shallow water and compared with field measurements. It is then applied to current flow about a conical island of small side slope where vortex shedding occurs and velocities are compared with data from the laboratory. Computed concentration distributions are also compared with dye visualization and an example of the influence of temperature on plume dispersion is presented. © 1997 John Wiley & Sons, Ltd.     

Phase-field modeling of interfacial dynamics in emulsion flows: Nonequilibrium surface tension

A weakly nonlocal phase-field model is used to define the surface tension in liquid binary mixtures in terms of the composition gradient in the interfacial region so that, at equilibrium, it depends linearly on the characteristic length that defines the interfacial width. Contrary to previous works suggesting that the surface tension in a phase-field model is fixed, we define the surface tension for a curved interface and far-from-equilibrium conditions as the integral of the free energy excess (i.e., above the thermodynamic component of the free energy) across the interface profile in a direction parallel to the composition gradient. Consequently, the nonequilibrium surface tension can be widely different from its equilibrium value under dynamic conditions, while it reduces to its thermodynamic value for a flat interface at local equilibrium. In nonequilibrium conditions, the surface tension changes with time: during mixing, it decreases as the inverse square root of time, while in the linear regime of spinodal decomposition, it increases exponentially to its equilibrium value, as nonlinear effects saturate the exponential growth. In addition, since temperature gradients modify the steepness of the concentration profile in the interfacial region, they induce gradients in the nonequilibrium surface tension, leading to the Marangoni thermocapillary migration of an isolated drop. Similarly, Marangoni stresses are induced in a composition gradient, leading to diffusiophoresis. We also review results on the nonequilibrium surface tension for a wall-bound pendant drop near detachment, which help to explain a discrepancy between our numerically determined static contact angle dependence of the critical Bond number and its sharp-interface counterpart from a static stability analysis of equilibrium shapes after numerical integration of the Young-Laplace equation. Finally, we present new results from phase-field simulations of the motion of an isolated droplet down an incline in gravity, showing that dynamic contact angle hysteresis can be explained in terms of the nonequilibrium surface tension.     

Solidification of a binary mixture in a square cavity with a free surface

Solidification of an NH₄Cl-H₂O solution in a square cavity with a free surface has been studied numerically and experimentally. Solidification is characterized by the coexistence of solid, mushy (solid plus liquid), and liquid regions, and flow within the liquid and mushy regions is driven by surface tension gradients, as well as by opposing thermal/solutal buoyancy forces. Theoretical and experimental results indicate that (i) surface tension forces spawn a cellular flow at the top of the cavity and enhance its development with increasing time, (ii) water enrichment of the top cell induces remelting of adjoining solid, and (iii) final steady-state conditions are characterized by a liquidus front whose thickness increases with increasing depth. However, predicted and observed results differ with respect to the chronological pattern of multicellular flow development and the time required to achieve steady-state conditions. Differences are attributed to selected model assumptions and to uncertainties concerning mushy region structure.