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.

Honey flows through fertile valleys : The cognitive and evolutionary foundations of paradise representations

The aim of the study is to explain how paradise beliefs are born from the viewpoint of mental functions of the human mind. The focus is on the observation that paradise beliefs across the world are mutually more similar than dissimilar. By using recent theories and results from the cognitive and evolutionary study of religion as well as from studies of environmental preferences, I suggest that this is because pan-human unconscious motivations, the architecture of mind, and the way the human mind processes information constrain the possible repertoire of paradise beliefs. The study is divided into two parts, theoretical and empirical. The arguments in the theoretical part are tested with data in the empirical part with two data sets. The first data set was collected using an Internet survey. The second data set was derived from literary sources. The first data test the assumption that intuitive conceptions of an environment of dreams generally follow the outlines set by evolved environmental preferences, but that they can be tweaked by modifying the presence of desirable elements. The second data test the assumption that familiarity is a dominant factor determining the content of paradise beliefs. The results of the study show that in addition to the widely studied belief in supernatural agents, belief in supernatural environments wells from the natural functioning of the human mind attesting the view that religious thinking and ideas are natural for human species and are produced by the same mental mechanisms as other cultural information. The results also help us to understand that the mental structures behind the belief in the supernatural have a wider scope than has been previously acknowledged.

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.

Structure and dynamics of two-dimensional sheared granular flows

The structure and dynamics of the two-dimensional linear shear flow of inelastic disks at high area fractions are analyzed. The event-driven simulation technique is used in the hard-particle limit, where the particles interact through instantaneous collisions. The structure (relative arrangement of particles) is analyzed using the bond-orientational order parameter. It is found that the shear flow reduces the order in the system, and the order parameter in a shear flow is lower than that in a collection of elastic hard disks at equilibrium. The distribution of relative velocities between colliding particles is analyzed. The relative velocity distribution undergoes a transition from a Gaussian distribution for nearly elastic particles, to an exponential distribution at low coefficients of restitution. However, the single-particle distribution function is close to a Gaussian in the dense limit, indicating that correlations between colliding particles have a strong influence on the relative velocity distribution. This results in a much lower dissipation rate than that predicted using the molecular chaos assumption, where the velocities of colliding particles are considered to be uncorrelated.

Linking the populations of barramundi and king threadfin to environmental flows in four rivers of tropical Australia

Linking the populations of barramundi and king threadfin to environmental flows in four rivers of tropical Australia

Invariance group properties and exact solutions of equations describing time-dependent free surface flows under gravity

Using the method of infinitesimal transformations, a 6-parameter family of exact solutions describing nonlinear sheared flows with a free surface are found. These solutions are a hybrid between the earlier self-propagating simple wave solutions of Freeman, and decaying solutions of Sachdev. Simple wave solutions are also derived via the method of infinitesimal transformations. Incomplete beta functions seem to characterize these (nonlinear) sheared flows in the absence of critical levels.

Dynamics of dense sheared granular flows. Part 1. Structure and diffusion

Shear flows of inelastic spheres in three dimensions in the Volume fraction range 0.4-0.64 are analysed using event-driven simulations.Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to the line Joining the centres). Here, we have considered both e(t) = +1 and e(t) = e(n) (rough particles) and e(t) =-1 (smooth particles), and the normal coefficient of restitution e(n) was varied in the range 0.6-0.98. Care was taken to avoid inelastic collapse and ensure there are no particle overlaps during the simulation. First, we studied the ordering in the system by examining the icosahedral order parameter Q(6) in three dimensions and the planar order parameter q(6) in the plane perpendicular to the gradient direction. It was found that for shear flows of sufficiently large size, the system Continues to be in the random state, with Q(6) and q(6) close to 0, even for volume fractions between phi = 0.5 and phi = 0.6; in contrast, for a system of elastic particles in the absence of shear, the system orders (crystallizes) at phi = 0.49. This indicates that the shear flow prevents ordering in a system of sufficiently large size. In a shear flow of inelastic particles, the strain rate and the temperature are related through the energy balance equation, and all time scales can be non-dimensionalized by the inverse of the strain rate. Therefore, the dynamics of the system are determined only by the volume fraction and the coefficients of restitution. The variation of the collision frequency with volume fraction and coefficient of estitution was examined. It was found, by plotting the inverse of the collision frequency as a function of volume fraction, that the collision frequency at constant strain rate diverges at a volume fraction phi(ad) (volume fraction for arrested dynamics) which is lower than the random close-packing Volume fraction 0.64 in the absence of shear. The volume fraction phi(ad) decreases as the coefficient of restitution is decreased from e(n) = 1; phi(ad) has a minimum of about 0.585 for coefficient of restitution e(n) in the range 0.6-0.8 for rough particles and is slightly larger for smooth particles. It is found that the dissipation rate and all components of the stress diverge proportional to the collision frequency in the close-packing limit. The qualitative behaviour of the increase in the stress and dissipation rate are well Captured by results derived from kinetic theory, but the quantitative agreement is lacking even if the collision frequency obtained from simulations is used to calculate the pair correlation function used In the theory.

Dynamics of dense sheared granular flows. Part II. The relative velocity distributions

The distribution of relative velocities between colliding particles in shear flows of inelastic spheres is analysed in the Volume fraction range 0.4-0.64. Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to line joining the centres). The distribution or pre-collisional normal relative velocities (along the line Joining the centres of the particles) is Found to be an exponential distribution for particles with low normal coefficient of restitution in the range 0.6-0.7. This is in contrast to the Gaussian distribution for the normal relative velocity in all elastic fluid in the absence of shear. A composite distribution function, which consists of an exponential and a Gaussian component, is proposed to span the range of inelasticities considered here. In the case of roughd particles, the relative velocity tangential to the surfaces at contact is also evaluated, and it is found to be close to a Gaussian distribution even for highly inelastic particles.Empirical relations are formulated for the relative velocity distribution. These are used to calculate the collisional contributions to the pressure, shear stress and the energy dissipation rate in a shear flow. The results of the calculation were round to be in quantitative agreement with simulation results, even for low coefficients of restitution for which the predictions obtained using the Enskog approximation are in error by an order of magnitude. The results are also applied to the flow down an inclined plane, to predict the angle of repose and the variation of the volume fraction with angle of inclination. These results are also found to be in quantitative agreement with previous simulations.

Solutions of Boundary Layer Equations Governing Free Convective and Magnetohydrodynamic Flows through Parametric Differentiation

In der vorliegenden Arbeit wird die Methode der parametrischen Differentiation angewendet, um ein System nichtlinearer Gleichungen zu lösen, das zwei- und dreidimensionale freie, konvektive Grenzschichströmungen bzw. eine zweidimensionale magnetohydrodynamische Grenzschichtströmung beherrscht. Der Hauptvorteil dieser Methode besteht darin, daß die nichlinearen Gleichungen auf lineare reduziert werden und die Nichtlinearität auf ein System von Gleichungen erster Ordnung beschränkt wird, das, verglichen mit den ursprünglichen Nichtlinearen Gleichungen, viel leichter gelöst werden kann. Ein anderer Vorzug der Methode ist, daß sie es ermöglicht, die Lösung von einer bekannten, zu einem bestimmten Parameterwert gehörigen Lösung aus durch schrittweises Vorgehen die Lösung für den gesamten Parameterbereich zu erhalten. Die mit dieser Methode gewonnenen Ergebnisse stimmen gut mit den entsprechenden, mit anderen numerischen Verfahren erzielten überein.

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.

The effect of suction on boundary layer for rotating flows with or without magnetic field

The effect of suction on the steady laminar incompressible boundarylayer flow for a stationary infinite disc with or without magnetic field, when the fluid at a large distance from the surface of the disc undergoes a solid body rotation, has been studied. The governing coupled nonlinear equations have been solved numerically using the shooting method with least square convergence criterion. It has been found that suction tends to reduce the velocity overshoot and damp the oscillation.

Two temperature viscous accretion flows around rotating black holes: Description of under-fed systems to ultra-luminous X-ray sources

We discuss two temperature accretion disk flows around rotating black holes. As we know that to explain observed hard X-rays the choice of Keplerian angular momentum profile is not unique, we consider the sub-Keplerian regime of the disk. Without any strict knowledge of the magnetic field structure, we assume the cooling mechanism is dominated by bremsstrahlung process. We show that in a range of Shakura-Sunyaev viscosity parameter 0.2 greater than or similar to alpha greater than or similar to 0.0005, flow behavior varies widely, particularly by means of the size of disk, efficiency of cooling and corresponding temperatures of ions and electrons. We also show that the disk around a rotating black hole is hotter compared to that around a Schwarzschild black hole, rendering a larger difference between ion and electron temperatures in the former case. With all the theoretical solutions in hand, finally we reproduce the observed luminosities (L) of two extreme cases-the under-fed AGNs and quasars (e.g. Sgr A') with L greater than or similar to 10(33) erg/s to ultra-luminous X-ray sources with L similar to 10(41) erg/s, at different combinations of mass accretion rate, ratio of specific heats, Shakura-Sunyaev viscosity parameter and Kerr parameter, and conclude that Sgr A' may be an intermediate spinning black hole.

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.