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.

On the relationship between finger width, velocity, and fluxes in thermohaline convection

Double-diffusive finger convection occurs in many natural processes.The theories for double-diffusive phenomena that exist at present consider systems with linear stratification in temperature and salinity. The double-diffusive systems with step change in salinity and temperature are, however, not amenable to simple stability analysis. Hence factors that control the width of the finger, velocity, and fluxes in systems that have step change in temperature and salinity have not been understood so far. In this paper we provide new physical insight regarding factors that influence finger convection in two-layer double-diffusive system through two-dimensional numerical simulations. Simulations have been carried out for density stability ratios (R-rho) from 1.5 to 10. For each density stability ratio, the thermal Rayleigh number (Ra-T) has been systematically varied from 7x10(3) to 7x10(8). Results from these simulations show how finger width, velocity, and flux ratios in finger convection are interrelated and the influence of governing parameters such as density stability ratio and the thermal Rayleigh number. The width of the incipient fingers at the time of onset of instability has been shown to vary as Ra-T-1/3. Velocity in the finger varies as Ra(T)1/3/R-rho. Results from simulation agree with the scale analysis presented in the paper. Our results demonstrate that wide fingers have lower velocities and flux ratios compared to those in narrow fingers. This result contradicts present notions about the relation between finger width and flux ratio. A counterflow heat-exchanger analogy is used in understanding the dependence of flux ratio on finger width and velocity.

Axially Loaded Timoshenko Rotors From A Prestressed Continuum Approach

We consider an axially loaded Timoshenko rotor rotating at a constant speed and derive its governing equations from a continuum viewpoint. The primary aim of this paper is to understand the source and role of gyroscopic terms, when the rotor is viewed not as a Timoshenko beam but as a genuine 3D continuum. We offer the primary insight that macroscopically observed gyroscopic terms may also, quite equivalently, be viewed as external manifestations of internally existing spin-induced prestresses at the continuum level. To demonstrate this idea with an analytical example (the Timoshenko rotor), we have studied the reliable equations of Choi et al. (Journal of Vibration and Acoustics, 114, 1992, 249-259). Using a straightforward application of our insight in the framework of nonlinear elasticity, we obtain equations that exactly match Choi et al. for the case with no axial load. For the case of axial preload, our straightforward formulation leads to a slightly different set of equations that have negligible numerical consequence for solid rotors. However, we offer a macroscopic, intuitive, justification for modifying our formulation so as to obtain the exact equations of Choi et al. with the axial load included.

Extended Kalman filters using explicit and derivative-free local linearizations

We propose three variants of the extended Kalman filter (EKF) especially suited for parameter estimations in mechanical oscillators under Gaussian white noises. These filters are based on three versions of explicit and derivative-free local linearizations (DLL) of the non-linear drift terms in the governing stochastic differential equations (SDE-s). Besides a basic linearization of the non-linear drift functions via one-term replacements, linearizations using replacements through explicit Euler and Newmark expansions are also attempted in order to ensure higher closeness of true solutions with the linearized ones. Thus, unlike the conventional EKF, the proposed filters do not need computing derivatives (tangent matrices) at any stage. The measurements are synthetically generated by corrupting with noise the numerical solutions of the SDE-s through implicit versions of these linearizations. In order to demonstrate the effectiveness and accuracy of the proposed methods vis-à-vis the conventional EKF, numerical illustrations are provided for a few single degree-of-freedom (DOF) oscillators and a three-DOF shear frame with constant parameters.

Oxygen transfer and energy dissipation rate in surface aerator

The dissipation rate of turbulent kinetic energy (epsilon) is a key parameter for mixing in surface aerators. In particular, determination epsilon across the impeller stream, where the most intensive mixing takes place, is essential to ascertain that an appropriate degree of mixing is achieved. Present work by using commercial software VisiMix (R) calculates the energy dissipation rate in geometrically similar unbaffled surface aeration systems in order to scale-up the oxygen transfer process. It is found that in geometrically similar system, oxygen transfer rate is uniquely correlated with dissipation rate of energy. Simulation or scale-up equation governing oxygen transfer rate and dissipation rate of energy has been developed in the present work.

Oxygen transfer and energy dissipation rate in surface aerator

The dissipation rate of turbulent kinetic energy(e)is a key parameter for mixing in surface aerators. In particular, determination e across the impeller stream, where the most intensive mixing takes place, is essential to ascertain that an appropriate degree of mixing is achieved. Present work by using commercial software VisiMix calculates the energy dissipation rate in geometrically similar unbaffled surface aeration systems in order to scale-up the oxygen transfer process. It is found that in geometrically similar system,oxygen transfer rate is uniquely correlated with dissipation rate of energy. Simulation or scale-up equation governing oxygen transfer rate and dissipation rate of energy has been developed in the present work.

Unsteady Natural Convection Flow in a Square Cavity Filled with a Porous Medium Due to Impulsive Change in Wall Temperature

Unsteady natural convection flow in a two- dimensional square cavity filled with a porous material has been studied. The flow is initially steady where the left- hand vertical wall has temperature T-h and the right- hand vertical wall is maintained at temperature T-c ( T-h > T-c) and the horizontal walls are insulated. At time t > 0, the left- hand vertical wall temperature is suddenly raised to (T-h) over bar ((T-h) over bar > T-h) which introduces unsteadiness in the flow field. The partial differential equations governing the unsteady natural convection flow have been solved numerically using a finite control volume method. The computation has been carried out until the final steady state is reached. It is found that the average Nusselt number attains a minimum during the transient period and that the time required to reach the final steady state is longer for low Rayleigh number and shorter for high Rayleigh number.

MHD Flow of a 2nd Grade Fluid in an Orthogonal Rheometer

The magnetohydrodynamics (MHD) flow of a conducting, homogeneous incompressible Rivlin-Ericksen fluid of second grade contained between two infinite, parallel, insulated disks rotating with the same angular velocity about two noncoincident axes, under the application of a uniform transverse magnetic field, is investigated. This model represents the MHD flow of the fluid in the instrument called an orthogonal rheometer, except for the fact that in the rheometer the rotating plates are necessarily finite. An exact solution of the governing equations of motion is presented. The force components in the x and y directions on the disks are calculated. The effects of magnetic field and the viscoelastic parameter on the forces are discussed in detail.

Unsteady incompressible boundary layer flow of a micropolar fluid at a stagnation point

The flow, heat and mass transfer on the unsteady laminar incompressible boundary layer in micropolar fluid at the stagnation point of a 2-dimensional and an axisymmetric body have been studied when the free stream velocity and the wall temperature vary arbitrarily with time. The partial defferential equations governing the flow have been solved numerically using a quasilinear finite-difference scheme. The skin friction, microrotation gradient and heat transfer parameters are found to be strongly dependent on the coupling parameter, mass transfer and time, whereas the effect of the microrotation parameter on the skin friction and heat transfer is rather weak, but microrotation gradient is strongly affected by it. The Prandtl number and the variation of the wall temperature with time affect the heat-transfer very significantly but the skin friction and micrortation gradient are unaffected by them.

Unsteady flow with attenuation in a fluid filled elastic tube with a stenosis

An equation governing the excess pressure has been derived, for an axially tethered and stenosed elastic tube filled with viscous liquid, by introducing the elasticity of the tube through pressure-area relation. This equation is solved numerically for large Womersley parameter and the results are presented for different types of pressure-radius relations and geometries by prescribing an outgoing wave suffering attenuation at some axial point of the tube. For a locally constricted tube it is observed that the pressure oscillates more and generates sound on the down stream side of the constriction.

Compressible boundary-layer flow at a three-dimensional stagnation point with massive blowing

The effect of massive blowing rates on the steady laminar compressible boundary-layer flow with variable gas properties at a 3-dim. stagnation point (which includes both nodal and saddle points of attachment) has been studied. The equations governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique for nodal points of attachment but employing a parametric differentiation technique instead of quasilinearization for saddle points of attachment. It is found that the effect of massive blowing rates is to move the viscous layer away from the surface. The effect of the variation of the density- viscosity product across the boundary layer is found to be negligible for massive blowing rates but significant for moderate blowing rates. The velocity profiles in the transverse direction for saddle points of attachment in the presence of massive blowing show both the reverse flow as well as velocity overshoot.

Analytical studies of gain optimization in CO2-N2 gasdynamic lasers employing two-dimensional wedge nozzles

An analytical method has been proposed to optimise the small-signaloptical gain of CO2-N2 gasdynamic lasers (gdl) employing two-dimensional (2D) wedge nozzles. Following our earlier work the equations governing the steady, inviscid, quasi-one-dimensional flow in the wedge nozzle of thegdl are reduced to a universal form so that their solutions depend on a single unifying parameter. These equations are solved numerically to obtain similar solutions for the various flow quantities, which variables are subsequently used to optimize the small-signal-gain. The corresponding optimum values like reservoir pressure and temperature and 2D nozzle area ratio also have been predicted and graphed for a wide range of laser gas compositions, with either H2O or He as the catalyst. A large number of graphs are presented which may be used to obtain the optimum values of small signal gain for a wide range of laser compositions without further computations.

Unsteady Attachment Line Flow On A Flat-Plate With Attached Cylinder

The unsteady laminar incompressible boundary-layer attachment-line flow on a flat plate with attached cylinder with heat and mass transfer has been studied when the free stream velocity, mass transfer and surface wall temperature vary arbitrarily with time. The governing partial differential equations with three independent variables have been solved numerically using an implicit finite-difference scheme. The heat transfer was found to be strongly dependent on the Prandtl number, variation of wall temperature with time and dissipation parameter (for large times). However, the free stream velocity distribution and mass transfer affect both the heat transfer and skin friction.

Bending Of Circular Plate On Elastic-Foundation

The governing differential equation of linear, elastic, thin, circular plate of uniform thickness, subjected to uniformly distributed load and resting on Winkler-Pasternak type foundation is solved using ``Chebyshev Polynomials''. Analysis is carried out using Lenczos' technique, both for simply supported and clamped plates. Numerical results thus obtained by perturbing the differential equation for plates without foundation are compared and are found to be in good agreement with the available results. The effect of foundation on central deflection of the plate is shown in the form of graphs.

Scheme of squares: Part I. Systems formalism for potentiostatic studies

The so-called “Scheme of Squares”, displaying an interconnectivity of heterogeneous electron transfer and homogeneous (e.g., proton transfer) reactions, is analysed. Explicit expressions for the various partial currents under potentiostatic conditions are given. The formalism is applicable to several electrode geometries and models (e.g., semi-infinite linear diffusion, rotating disk electrodes, spherical or cylindrical systems) and the analysis is exact. The steady-state (t→∞) expressions for the current are directly given in terms of constant matrices whereas the transients are obtained as Laplace transforms that need to be inverted by approximation of numerical methods. The methodology employs a systems approach which replaces a system of partial differential equations (governing the concentrations of the several electroactive species) by an equivalent set of difference equations obeyed by the various partial currents.