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.

Role of surface tension driven convection in pulsed - laser melting and solidification

Surface melting by a stationary, pulsed laser has been modelled by the finite element method. The role of the surface tension driven convection is investigated in detail. Numerical results are presented for a triangular laser pulse of durations 10, 50 and 200 ms. Though the magnitude of the velocity is high due to the surface tension forces, the present results indicate that a finite time is required for convection to affect the temperature distribution within the melt pool. The effect of convection is very significant for pulse durations longer than 10 ms.

A regular perturbation approach to surface tension driven flows

Surface tension induced convection in a liquid bridge held between two parallel, coaxial, solid disks is considered. The surface tension gradient is produced by a small temperature gradient parallel Co the undisturbed surface. The study is performed by using a mathematical regular perturbation approach based on a small parameter, e, which measures the deviation of the imposed temperature field from its mean value. The first order velocity field is given by a Stokes-type problem (viscous terms are dominant) with relatively simple boundary conditions. The first order temperature field is that imposed from the end disks on a liquid bridge immersed in a non-conductive fluid. Radiative effects are supposed to be negligible. The second order temperature field, which accounts for convective effects, is split into three components, one due to the bulk motion, and the other two to the distortion of the free surface. The relative importance of these components in terms of the heat transfer to or from the end disks is assessed

Surface tension and microgravity

The behaviour of confined liquids on board an orbiting spacecraft is mainly driven by surface tension phenomena, which cause an apparently anomalous response of the liquid when compared with the behaviour that can be observed on an Earth laboratory provided that the amount of liquid is high enough. The reason is that in an orbiting spacecraft the different inertial forces acting on the bulk of the liquid are almost zero, causing thus capillary forces to be the dominant ones. Of course, since gravity forces are proportional to the liquid volume, whereas surface tension forces are proportional to the liquid surface, there are situations on Earth where capillarity can be the dominant effect, as it happens when very small volume liquid samples are considered. However, work with small size samples may require the use of sophisticated optical devices. Leaving aside the neutral buoyancy technique, a way of handling large liquid interfaces is by using drop towers, where the sample falls subjected to the action of Earth's gravity. This approach is suitable when the characteristic time of the problem under consideration is much smaller than the drop time. In this work the transformation of an out-of-use chimney into a drop tower is presented. Because of the miniaturization, hardiness and low cost of current electronic devices, a drop tower can be used as an inexpensive tool for undergraduate students to experimentally analyse a large variety of surface tension driven phenomena.

The dynamics of a topological change in a system of soap films

Hutzler, S., Saadatfar, M., van der Net, A., Weaire, D. and Cox, S.J. (2007) The dynamics of a topological change in a system of soap films. Coll. Surf. A, 323:123-131. Sponsorship: This research was supported by the European Space Agency (contracts 14914/02/NL/SH and 14308/00/NL/SH), Science Foundation Ireland. (RFP 05/REP/PHY00/6), and the EU program COST P21 (The Physics of droplets). SJC acknowledges support from EPSRC (EP/D071127/1). MS is supported by the Irish Higher Education Authority (PRTLI-IITAC).

The effect of viscosity and surface tension on the Taylor instability of the surface of a cavitation bubble*

This paper is devoted to a consideration of the following problem: A spherical mass of fluid of density varrho1, viscosity μ1 and external radius R is surrounded by a fluid of density varrho2 and viscosity μ2.The fluids are immiscible and incompressible. The interface is accelerated radially by g1: to study the effect of viscosity and surface tension on the stability of the interface. By analyzing the problem in spherical harmonics the mathematical problem is reduced to one of solution of the characteristic determinant equation. The particular case of a cavity bubble, where the viscosity μ1 of the fluid inside the bubble is negligible in comparison with the viscosity μ2 of the fluid outside the bubble, is considered in some detail. It is shown that viscosity has a stabilizing role on the interface; and when g1 > T(n − 1) (n + 2)/R2(varrho2 − varrho1) the stabilizing role of both viscosity and surface tension is more pronounced than would result when either of them is taken individually.

Surface tension-induced global instability of planar jets and wakes

The effect of surface tension on global stability of co-flow jets and wakes at a moderate Reynolds number is studied. The linear temporal two-dimensional global modes are computed without approximations. All but one of the flow cases under study are globally stable without surface tension. It is found that surface tension can cause the flow to be globally unstable if the inlet shear (or equivalently, the inlet velocity ratio) is strong enough. For even stronger surface tension, the flow is re-stabilized. As long as there is no change of the most unstable mode, increasing surface tension decreases the oscillation frequency. Short waves appear in the high-shear region close to the nozzle, and their wavelength increases with increasing surface tension. The critical shear (the weakest inlet shear at which a global instability is found) gives rise to antisymmetric disturbances for the wakes and symmetric disturbances for the jets. However, at stronger shear, the opposite symmetry can be the most unstable one, in particular for wakes at high surface tension. The results show strong effects of surface tension that should be possible to reproduce experimentally as well as numerically.

Bubble extinction in Hele-Shaw flow with surface tension and kinetic undercooling regularisation

We perform an analytic and numerical study of an inviscid contracting bubble in a two-dimensional Hele-Shaw cell, where the effects of both surface tension and kinetic undercooling on the moving bubble boundary are not neglected. In contrast to expanding bubbles, in which both boundary effects regularise the ill-posedness arising from the viscous (Saffman-Taylor) instability, we show that in contracting bubbles the two boundary effects are in competition, with surface tension stabilising the boundary, and kinetic undercooling destabilising it. This competition leads to interesting bifurcation behaviour in the asymptotic shape of the bubble in the limit it approaches extinction. In this limit, the boundary may tend to become either circular, or approach a line or "slit" of zero thickness, depending on the initial condition and the value of a nondimensional surface tension parameter. We show that over a critical range of surface tension values, both these asymptotic shapes are stable. In this regime there exists a third, unstable branch of limiting self-similar bubble shapes, with an asymptotic aspect ratio (dependent on the surface tension) between zero and one. We support our asymptotic analysis with a numerical scheme that utilises the applicability of complex variable theory to Hele-Shaw flow.

A hybrid surface tension seal for pneumatic and hydraulic microactuators

Recent research revealed that microactuators driven by pressurized fluids are able to generate high power and force densities at microscale. One of the main technological barriers in the development of these actuators is the fabrication low friction seals. This paper presents a novel scalable seal technology, which resists the actuation pressure relying on a combination of a clearance seal and a surface tension seal. This approach allows to seal pressures of more than 800 kPa without leakage. The seal is tested on an actuator with a bore of 0.8 mm2 and a length of 13 mm, which was able to generate forces up to 0.32 N. © 2008 Springer-Verlag.

Three-Dimensional Linear Instability Analysis of Thermocapillary Return Flow on a Porous Plane

A three-dimensional linear instability analysis of thermocapillary convection in a fluid-porous double layer system, imposed by a horizontal temperature gradient, is performed. The basic motion of fluid is the surface-tension-driven return flow, and the movement of fluid in the porous layer is governed by Darcy's law. The slippery effect of velocity at the fluid-porous interface has been taken into account, and the influence of this velocity slippage on the instability characteristic of the system is emphasized. The new behavior of the thermocapillary convection instability has been found and discussed through the figures of the spectrum.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

The effect of surface tension and kinetic undercooling on a radially-symmetric melting problem

The addition of surface tension to the classical Stefan problem for melting a sphere causes the solution to blow up at a finite time before complete melting takes place. This singular behaviour is characterised by the speed of the solid-melt interface and the flux of heat at the interface both becoming unbounded in the blow-up limit. In this paper, we use numerical simulation for a particular energy-conserving one-phase version of the problem to show that kinetic undercooling regularises this blow-up, so that the model with both surface tension and kinetic undercooling has solutions that are regular right up to complete melting. By examining the regime in which the dimensionless kinetic undercooling parameter is small, our results demonstrate how physically realistic solutions to this Stefan problem are consistent with observations of abrupt melting of nanoscaled particles.

The effect of surface tension in porous wave maker problems

Using a mixed-type Fourier transform of a general form in the case of water of infinite depth and the method of eigenfunction expansion in the case of water of finite depth, several boundary-value problems involving the propagation and scattering of time harmonic surface water waves by vertical porous walls have been fully investigated, taking into account the effect of surface tension also. Known results are recovered either directly or as particular cases of the general problems under consideration.