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.

Effect of a magnetic field on the flow and blood oxygenation in channels of variable cross-section

The effect of a magnetic field on the flow and oxygenation of an incompressible Newtonian conducting fluid in channels with irregular boundaries has been investigated. The geometric parameter δ, which is a ratio of the mean half width of the channel d to the characteristic length λ along the channel over which the significant changes in the flow quantities occur, has been used for perturbing the governing equations. Closed form solutions of the various order equations are presented for the stream function. The equations for oxygen partial pressure remain nonlinear even after perturbation, therefore a numerical solution is presented. The expressions for shear stress at a wall and pressure distributions are derived. Here the separation in the flow occurs at a higher Reynolds number than the corresponding non-magnetic case. It is found that the magnetic field has an effect on local oxygen concentration but has a little effect on the saturation length.

Nonsimilar incompressible laminar boundary-layer flows in micropolar fluids

The solution of the steady laminar incompressible nonsimilar boundary-layer problem for micropolar fluids over two-dimensional and axisymmetric bodies has been presented. The partial differential equations governing the flow have been transformed into new co-ordinates having finite range. The resulting equations have been solved numerically using implicit finite-difference scheme. The computations have been carried out for a cylinder and a sphere. The results indicate that the separation in micropolar fluids occurs at earlier streamwise locations as compared to Newtonian fluids. The skin friction and velocity profiles depend on the shape of the body and are almost insensitive to microrotation or coupling parameter, provided the coupling parameter is small. On the other hand, the microrotation profiles and microrotation gradient depend on the microrotation parameter and they are insensitive to the coupling parameter.

The possible role of meridional flows in suppressing magnetic buoyancy

The equation of motion for a toroidal flux ring in a stellar convective envelope is derived, and the equilibrium of such a ring is considered. Necessary conditions for the stability of toroidal flux rings are derived, and results of stability calculations for a particular model of the meridional flow are presented. The motions of the flux rings when the rings are far from their equilibrium position or when equilibrium does not exist are considered. The results confirm the linear stability analysis, and show that in the absence of stable equilibrium, the rings move toward the solar surface along a trajectory which is parallel to the rotation axis. It is expected that viscosity will tend to reduce the rotational velocity difference between the flux ring and its surroundings, thus reducing the Coriolis force and altering the equilibrium. The storage time of toroidal flux rings is estimated, and some implications for the sun are discussed.

A class of one-dimensional steady-state flows in radiation-gas-dynamics

The study of steady-state flows in radiation-gas-dynamics, when radiation pressure is negligible in comparison with gas pressure, can be reduced to the study of a single first-order ordinary differential equation in particle velocity and radiation pressure. The class of steady flows, determined by the fact that the velocities in two uniform states are real, i.e. the Rankine-Hugoniot points are real, has been discussed in detail in a previous paper by one of us, when the Mach number M of the flow in one of the uniform states (at x=+∞) is greater than one and the flow direction is in the negative direction of the x-axis. In this paper we have discussed the case when M is less than or equal to one and the flow direction is still in the negative direction of the x-axis. We have drawn the various phase planes and the integral curves in each phase plane give various steady flows. We have also discussed the appearance of discontinuities in these flows.

What mainly controls recession flows in river basins?

The ubiquity of the power law relationship between dQ/dt and Q for recession periods (-dQ/dt kQ(alpha); Q being discharge at the basin outlet at time t) clearly hints at the existence of a dominant recession flow process that is common to all real basins. It is commonly assumed that a basin, during recession events, functions as a single phreatic aquifer resting on a impermeable horizontal bed or the Dupuit-Boussinesq (DB) aquifer, and with time different aquifer geometric conditions arise that give different values of alpha and k. The recently proposed alternative model, geomorphological recession flow model, however, suggests that recession flows are controlled primarily by the dynamics of the active drainage network (ADN). In this study we use data for several basins and compare the above two contrasting recession flow models in order to understand which of the above two factors dominates during recession periods in steep basins. Particularly, we do the comparison by selecting three key recession flow properties: (1) power law exponent alpha, (2) dynamic dQ/dt-Q relationship (characterized by k) and (3) recession timescale (time period for which a recession event lasts). Our observations suggest that neither drainage from phreatic aquifers nor evapotranspiration significantly controls recession flows. Results show that the value of a and recession timescale are not modeled well by DB aquifer model. However, the above mentioned three recession curve properties can be captured satisfactorily by considering the dynamics of the ADN as described by geomorphological recession flow model, possibly indicating that the ADN represents not just phreatic aquifers but the organization of various sub-surface storage systems within the basin. (C) 2014 Elsevier Ltd. All rights reserved.

Experimental evidence of the kinematic Cosserat effect in dense granular flows

To capture shear localization in the flow of dense granular materials in a continuum description, it has earlier been proposed that granular materials be treated as Cosserat, or micropolar, continua. Here, we provide experimental verification of the kinematic Cosserat effect, or the deviation of the particle spin from the material spin induced by the velocity gradient. Contrary to earlier belief, we find this effect to be sizable even outside the shear layers. Remarkably, the particles and material elements spin in opposite directions in flow through a hopper.

Instabilities induced by variation of Brunt-Vaisala frequency in compressible stratified shear flows

The stability characteristics of a Helmholtz velocity profile in a stably stratified, compressible fluid in the presence of a lower rigid boundary are studied. A jump in the Brunt-Vaisala frequency at a level different from the shear zone is introduced and the variation of the Brunt-Vaisala frequency with respect to the vertical coordinate in the middle layer of the three-layered model is considered. An analytic solution in each of the layers is obtained, and the dispersion relation is solved numerically for parameters relevant to the model. The effect of shear in the lowermost layer of the three-layered model for a Boussinesq fluid is discussed. The results are compared with the earlier studies of Lindzen and Rosenthal, and Sachdev and Satya Narayanan. In the present model, new unstable modes with larger growth rates are obtained and the most unstable gravity wave modes are found to agree closely with the observed ones at various heights. Physics of Fluids is copyrighted by The American Institute of Physics.

Drag reduction by a forward facing aerospike for a large angle blunt cone in high enthalpy flows

Drag reduction studies are conducted using a flat disc tipped aerospike for a 120-degree apex angle blunt cone model in high enthalpy flows. Accelerometer based force balance is used for the drag force measurement in the newly established free piston driven shock tunnel, HST3. Drag reduction upto about 58 percent has been achieved for Mach 8 flow of 5 MJ/kg specific enthalpy at zero degree angle of attack.

Drag reduction by counterflow supersonic jet for a blunt cone in high enthalpy flows

Counterflow supersonic jet is used as a drag reduction device during the experiments in free piston driven shock tunnel, HST3. Accelerometer based force balance is employed to measure the drag force experienced by the 60-degree apex angle blunt cone model without and with the supersonic jet opposing the hypersonic flow. It is observed that the drag force decreases with increase in injection pressure ratio until the critical injection pressure is reached. Maximum reduction in drag force of 44 percent is recorded at the critical injection pressure ratio 22.36. Further increase in injection pressure ratio has reduced the percentage drag reduction. Change in nature of the flowfield around the model has also been observed across the critical injection pressure ratio.

Time-dependent flows of rotating and stratified fluids in geometries with non-uniform cross-sections

Unsteady rotating and stratified flows in geometries with non-uniform cross-sections are investigated under Oseen approximation using Laplace transform technique. The solutions are obtained in closed form and they reveal that the flow remains oscillatory even after infinitely large time. The existence of inertial waves propagating in both positive and negative directions of the flow is observed. When the Rossby or Froude number is close to a certain infinite set of critical values the blocking and back flow occur and the flow pattern becomes more and more complicated with increasing number of stagnant zones when each critical value is crossed. The analogy that is observed in the solutions for rotating and stratified flows is also discussed.

A new method of producing local enhancement of buoyancy in liquid flows

We describe here a novel method of generating large volumetric heating in a liquid. The method uses the principle of ohmic heating of the liquid, rendered electrically conducting by suitable additives if necessary. Electrolysis is prevented by the use of high frequency alternating voltage and chemically treated electrodes. The technique is demonstrated by producing substantial heating in an initially neutral jet of water. Simple flow visualisation studies, made by adding dye to the jet, show marked changes in the growth and development of the jet with heat addition.