NONLINEAR SURFACE-WAVE EXCITATIONS IN THE BENARD-MARANGONI SYSTEM

We examine the appearance of surface waves governed by Burgers and Korteweg-de Vries equations in a shallow viscous heated fluid. We consider waves triggered by a surface-tension variation induced by both temperature and concentration gradients. We also establish the range of parameters for which the above-mentioned equations appear.

Modulational instability analysis of surface-waves in the Bénard-Marangoni phenomenon

By using the long-wave approximation, a system of coupled evolutions equations for the bulk velocity and the surface perturbations of a Bénard-Marangoni system is obtained. It includes nonlinearity, dispersion and dissipation, and it is interpreted as a dissipative generalization of the usual Boussinesq system of equations. Then, by considering that the Marangoni number is near the critical value M = -12, we show that the modulation of the Boussinesq waves is described by a perturbed Nonlinear Schrödinger Equation, and we study the conditions under which a Benjamin-Feir instability could eventually set in. The results give sufficient conditions for stability, but are inconclusive about the existence or not of a Benjamin-Feir instability in the long-wave limit. © 1995.

In situ studies of phase separation and crystallization directed by Marangoni instabilities during spin-coating

Results of a pioneering study are presented in which for the first time, crystallization, phase separation and Marangoni instabilities occurring during the spin-coating of polymer blends are directly visualized, in real-space and real-time. The results provide exciting new insights into the process of self-assembly, taking place during spin-coating, paving the way for the rational design of processing conditions, to allow desired morphologies to be obtained. © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Finite Element Solution Of Surface-Tension Driven Flows In Laser Surface-Melting

Lasers are very efficient in heating localized regions and hence they find a wide application in surface treatment processes. The surface of a material can be selectively modified to give superior wear and corrosion resistance. In laser surface-melting and welding problems, the high temperature gradient prevailing in the free surface induces a surface-tension gradient which is the dominant driving force for convection (known as thermo-capillary or Marangoni convection). It has been reported that the surface-tension driven convection plays a dominant role in determining the melt pool shape. In most of the earlier works on laser-melting and related problems, the finite difference method (FDM) has been used to solve the Navier Stokes equations [1]. Since the Reynolds number is quite high in these cases, upwinding has been used. Though upwinding gives physically realistic solutions even on a coarse grid, the results are inaccurate. McLay and Carey have solved the thermo-capillary flow in welding problems by an implicit finite element method [2]. They used the conventional Galerkin finite element method (FEM) which requires that the pressure be interpolated by one order lower than velocity (mixed interpolation). This restricts the choice of elements to certain higher order elements which need numerical integration for evaluation of element matrices. The implicit algorithm yields a system of nonlinear, unsymmetric equations which are not positive definite. Computations would be possible only with large mainframe computers.Sluzalec [3] has modeled the pulsed laser-melting problem by an explicit method (FEM). He has used the six-node triangular element with mixed interpolation. Since he has considered the buoyancy induced flow only, the velocity values are small. In the present work, an equal order explicit FEM is used to compute the thermo-capillary flow in the laser surface-melting problem. As this method permits equal order interpolation, there is no restriction in the choice of elements. Even linear elements such as the three-node triangular elements can be used. As the governing equations are solved in a sequential manner, the computer memory requirement is less. The finite element formulation is discussed in this paper along with typical numerical results.

Thermal instability in a three-dimensional rigid container with prescribed heat flux at lower boundary

The Bénard–Marangoni convection is studied in a three-dimensional container with thermally insulated lateral walls and prescribed heat flux at lower boundary. The upper surface of the incompressible, viscous fluid is assumed to be flat with temperature dependent surface tension. A Galerkin–Tau method with odd and even trial functions satisfying all the essential boundary conditions except the natural boundary conditions at the free surface has been used to solve the problem. The critical Marangoni and Rayleigh numbers are determined for the onset of steady convection as a function of aspect ratios x0 and y0 for the cases of Bénard–Marangoni, pure Marangoni and pure Bénard convections. It is observed that critical parameters are decreasing with an increase in aspect ratios. The flow structures corresponding to the values of the critical parameters are presented in all the cases. It is observed that the critical parameters are higher for case with heat flux prescribed than those corresponding to the case with prescribed temperature. The critical Marangoni number for pure Marangoni convection is higher than critical Rayleigh number corresponding to pure Bénard convection for a given aspect ratio whereas the reverse was observed for two-dimensional infinite layer.

Dewetting on the Surface of Rutile

After annealing a continuous SiO2 film on the (001) surface of TiO2, the film dewets and then spreads to form a complex pattern. The final droplet morphology displays a densely branching morphology similar to those seen in computer-simulated models. It is proposed that Bénard-Marangoni convection cells form within the film before dewetting occurs. The formation of Bénard-Marangoni convection cells prior to dewetting results in the uniform size and spacing of the droplets on the surface. These convection cells form at temperature when the TiO2 substrate dissolves into the SiO2 thin film. The change in composition results in regions of differing surface tensions and therefore leads to the formation of the convection cells.

Mechanism of enhanced corrosion of alumina graphite refractories near slag/metal interface

All refractories show enhanced corrosion near the slag/metal interface due to Marangoni and convective flows. However, in the case of oxide refractories containing graphite flakes, corrosion is severe due to periodic oscillations in the contact angle at the slag/metal interface, resulting in cyclic dissolution of oxide and graphite into the slag and metal, respectively. Alumina--graphite (AG) refractories should be used only where they are not in simultaneous contact with slag (flux) and low carbon steel.

A k-epsilon model for turbulent weld pool convection in gas metal arc welding process

In this paper, we present a modified k - epsilon model capable of addressing turbulent weld-pool convection in a GMAW process, taking into account the morphology of the phase change interface during a Gas Metal Arc Welding (GMAW) process. A three-dimensional turbulence mathematical model has been developed to study the heat transfer and fluid flow within the weld pool by considering the combined effect of three driving forces, viz., buoyancy, Lorentz force and surface tension (Marangoni convection). Mass and energy transports by the droplets are considered through the thermal analysis of the electrode. The falling droplet's heat addition to the molten pool is considered to be a volumetric heat source distributed in an imaginary cylindrical cavity ("cavity model") within the weld pool. This nature of heat source distribution takes into account the momentum and the thermal, energy of the falling droplets. The numerically predicted weld pool dimensions both from turbulence and laminar models are then compared with the experimental post-weld results sectioned across the weld axis. The above comparison enables us to analyze the overall effects of turbulent convection on the nature of heat and fluid flow and hence on the weld pool shape/size during the arc welding processes.

Arbitrary Lagrangian-Eulerian finite-element method for computation of two-phase flows with soluble surfactants

A finite-element scheme based on a coupled arbitrary Lagrangian-Eulerian and Lagrangian approach is developed for the computation of interface flows with soluble surfactants. The numerical scheme is designed to solve the time-dependent Navier-Stokes equations and an evolution equation for the surfactant concentration in the bulk phase, and simultaneously, an evolution equation for the surfactant concentration on the interface. Second-order isoparametric finite elements on moving meshes and second-order isoparametric surface finite elements are used to solve these equations. The interface-resolved moving meshes allow the accurate incorporation of surface forces, Marangoni forces and jumps in the material parameters. The lower-dimensional finite-element meshes for solving the surface evolution equation are part of the interface-resolved moving meshes. The numerical scheme is validated for problems with known analytical solutions. A number of computations to study the influence of the surfactants in 3D-axisymmetric rising bubbles have been performed. The proposed scheme shows excellent conservation of fluid mass and of the total mass of the surfactant. (C) 2012 Elsevier Inc. All rights reserved.

Simulations of impinging droplets with surfactant-dependent dynamic contact angle

An arbitrary Lagrangian-Eulerian (ALE) finite element scheme for computations of soluble surfactant droplet impingement on a horizontal surface is presented. The numerical scheme solves the time-dependent Navier-Stokes equations for the fluid flow, scalar convection-diffusion equation for the surfactant transport in the bulk phase, and simultaneously, surface evolution equations for the surfactants on the free surface and on the liquid-solid interface. The effects of surfactants on the flow dynamics are included into the model through the surface tension and surfactant-dependent dynamic contact angle. In particular, the dynamic contact angle (theta(d)) of the droplet is defined as a function of the surfactant concentration at the contact line and the equilibrium contact angle (theta(0)(e)) of the clean surface using the nonlinear equation of state for surface tension. Further, the surface forces are included into the model as surface divergence of the surface stress tensor that allows to incorporate the Marangoni effects without calculating the surface gradient of the surfactant concentration on the free surface. In addition to a mesh convergence study and validation of the numerical results with experiments, the effects of adsorption and desorption surfactant coefficients on the flow dynamics in wetting, partially wetting and non-wetting droplets are studied in detail. It is observed that the effects of surfactants are more in wetting droplets than in the non-wetting droplets. Further, the presence of surfactants at the contact line reduces the equilibrium contact angle further when theta(0)(e) is less than 90 degrees, and increases it further when theta(0)(e) is greater than 90 degrees. Nevertheless, the presence of surfactants has no effect on the contact angle when theta(0)(e) = 90 degrees. The numerical study clearly demonstrates that the surfactant-dependent contact angle has to be considered, in addition to the Marangoni effect, in order to study the flow dynamics and the equilibrium states of surfactant droplet impingement accurately. The proposed numerical scheme guarantees the conservation of fluid mass and of the surfactant mass accurately. (C) 2015 Elsevier Inc. All rights reserved.

Modeling of the non-isothermal liquid droplet impact on a heated solid substrate with heterogeneous wettability

A comprehensive numerical investigation on the impingement and spreading of a non-isothermal liquid droplet on a solid substrate with heterogeneous wettability is presented in this work. The time-dependent incompressible Navier-Stokes equations are used to describe the fluid flow in the liquid droplet, whereas the heat transfer in the moving droplet and in the solid substrate is described by the energy equation. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the time-dependent incompressible Navier-Stokes equation and the energy equation in the time-dependent moving domain. Moreover, the Marangoni convection is included in the variational form of the Navier-Stokes equations without calculating the partial derivatives of the temperature on the free surface. The heterogeneous wettability is incorporated into the numerical model by defining a space-dependent contact angle. An array of simulations for droplet impingement on a heated solid substrate with circular patterned heterogeneous wettability are presented. The numerical study includes the influence of wettability contrast, pattern diameter, Reynolds number and Weber number on the confinement of the spreading droplet within the inner region, which is more wettable than the outer region. Also, the influence of these parameters on the total heat transfer from the solid substrate to the liquid droplet is examined. We observe that the equilibrium position depends on the wettability contrast and the diameter of the inner surface. Consequently. the heat transfer is more when the wettability contrast is small and/or the diameter of inner region is large. The influence of the Weber number on the total heat transfer is more compared to the Reynolds number, and the total heat transfer increases when the Weber number increases. (C) 2015 Elsevier Ltd. All rights reserved.

Experimental Investigations on Interaction of Two Drops by Thermocapillary-Buoyancy Migration

Experiments were performed, in a terrestrial environment, to study the migration and interaction of two drops with different diameters in matrix liquid under temperature gradient field. Pure soybean oil and silicon oil were used as matrix liquid and the drop liquid, respectively. The information on the motions of two drops was recorded by CCD camera system in the experiments to analyze the trajectories and velocities of the drops. Our experiments showed that, upon two drops approaching each other, the influence of the larger drop on the motion of the smaller one became significant. Meanwhile the smaller drop had a little influence on the larger one all the time. The oscillation of migration velocities of both drops was observed as they were approaching. For a short period the smaller drop even moved backward when it became side by side with the larger one during the migration. Although our experimental results on the behavior of two drops are basically consistent with the theoretical predictions, there are also apparent differences. 2006 Elsevier Ltd. All rights reserved. Keywords: Thermocapillary migration; Drop; Interaction; Oscillation 1. Introduction A bubble or drop will move when placed in another fluid with temperature gradient. This motion happens as a consequence of the variation of interfacial tension with temperature. Such a phenomenon is already known as Marangoni migration problem. With the development of microgravity science, bubble dynamics and droplet dynamics became a hot point problem of research because this investigation is very important for basic research as well as for applications in reduced gravity environment, such as space material science, chemical engineering and so on. Young et al. first investigated the thermocapillary migration of

Experimental investigation of thermocapillary migration of isolated drops

Experimental hardware has been developed to perform experiments on the Marangoni migration of drops in the case of intermediate Reynolds numbers in a microgravity environment. The experiments were conducted using the drop shaft free fall facility with a 4.5 second microgravity period in the Microgravity Laboratory of Japan. In this experiment, the thermocapillary velocity of drop migration was measured for drops of different sizes in a series of temperature gradients.

两层流体热毛细对流空间实验研究=Space Experiments of Thermo Capillary Convection in Two-liquid Layers

1999年,在我国实践5号卫星上完成了两层流体空间实验,实验研究两层不相混合流体的纯Marangoni对流（温度梯度与界面垂直）与热毛细对流（温度梯度方向与流体界面平行）.前者存在发生Marangoni对流的最小临界温差值△Tc,低于该值流体系统处于静止状态;后者中只要存在沿界面的温度梯度便会产生热毛细对流.空间实验采用石蜡和氟化液两层流体新体系,实现了平整的液-液交界面,并从卫星上传回上万幅数字图像.通过多幅图像叠加处理得到了定量的流速场.数值模拟计算分析了相应工况时对流流动的速度场,两者的流场结构和速度大小基本一致,实验验证了理论模型.

Effect of liquid bridge volume on the instability in small-Prandtl-number half zones

A perturbation method is used to examine the linear instability of thermocapillary convection in a liquid bridge of floating half-zone filled with a small Prandtl number fluid. The influence of liquid bridge volume on critical Marangoni number and flow features is analyzed. The neutral modes show that the instability is mainly caused by the bulk flow that is driven by the nonuniform thermocapillary forces acting on the free surface. The hydrodynamic instability is dominant in the case of small Prandtl number fluid and the first instability mode is a stationary bifurcation. The azimuthal wave number for the most dangerous mode depends on the liquid bridge volume, and is not always two as in the case of a cylindrical liquid bridge with aspect ratio near 0.6. Its value may be equal to unity when the liquid bridge is relatively slender.