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.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter = 9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 10(8). Under these conditions the convection is turbulent, and the time-averaged velocity at any point is `zero'. The Reynolds number based on the Taylor microscale, Re-lambda, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6-7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as (RaSc1/2)-Sc-1/2, and the Reynolds number would scale as (RaSc-1/2)-Sc-1/2. The velocity and the flux measurements appear to be consistent with the Ra-1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are Much higher compared to what would be obtained in Rayleigh-Benard (R-B) convection for similar density differences.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter=9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 108. Under these conditions the convection is turbulent, and the time-averaged velocity at any point is ‘zero’. The Reynolds number based on the Taylor microscale, Reλ, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6–7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as Ra1/2Sc1/2, and the Reynolds number would scale as Ra1/2Sc−1/2. The velocity and the flux measurements appear to be consistent with the Ra1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are much higher compared to what would be obtained in Rayleigh–Bénard (R–B) convection for similar density differences.

Quasiuniversal Transient Behavior of a Nonequilibrium Mott Insulator Driven by an Electric Field

We use a self-consistent strong-coupling expansion for the self-energy (perturbation theory in the hopping) to describe the nonequilibrium dynamics of strongly correlated lattice fermions. We study the three-dimensional homogeneous Fermi-Hubbard model driven by an external electric field showing that the damping of the ensuing Bloch oscillations depends on the direction of the field and that for a broad range of field strengths a long-lived transient prethermalized state emerges. This long-lived transient regime implies that thermal equilibrium may be out of reach of the time scales accessible in present cold atom experiments but shows that an interesting new quasiuniversal transient state exists in nonequilibrium governed by a thermalized kinetic energy but not a thermalized potential energy. In addition, when the field strength is equal in magnitude to the interaction between atoms, the system undergoes a rapid thermalization, characterized by a different quasiuniversal behavior of the current and spectral function for different values of the hopping. DOI: 10.1103/PhysRevLett.109.260402

Manually operated piston-driven shock tube

A simple hand-operated shock tube capable of producing Mach 2 shock waves is described. Performance of this miniature shock tube using compressed high pressure air created by a manually operated piston in the driver section of the shock tube as driver gas with air at 1 atm pressure as the test gas in the driven tube is presented. The performance of the shock tube is found to match well with the theoretically estimated values using normal shock relations. Applications of this shock tube named Reddy tube, include study of blast-induced traumatic brain injuries and high temperature chemical kinetics.

Dynamics of a driven magneto-martensitic ribbon

We investigate the dynamics of a sinusoidally driven ferromagnetic martensitic ribbon by adopting a recently introduced model that involves strain and magnetization as order parameters. Retaining only the dominant mode of excitation we reduce the coupled set of partial differential equations for strain and magnetization to a set of coupled ordinary nonlinear equations for the strain and magnetization amplitudes. The equation for the strain amplitude takes the form of parametrically driven oscillator. Finite strain amplitude can only be induced beyond a critical value of the strength of the magnetic field. Chaotic response is seen for a range of values of all the physically interesting parameters. The nature of the bifurcations depends on the choice of temperature relative to the ordering of the Curie and the martensite transformation temperatures. We have studied the nature of response as a function of the strength and frequency of the magnetic field, and magneto-elastic coupling. In general, the bifurcation diagrams with respect to these parameters do not follow any standard route. The rich dynamics exhibited by the model is further illustrated by the presence of mixed mode oscillations seen for low frequencies. The geometric structure of the mixed mode oscillations in the phase space has an unusual deep crater structure with an outer and inner cone on which the orbits circulate. We suggest that these features should be seen in experiments on driven magneto-martensitic ribbons. (C) 2014 Elsevier B. V. All rights reserved.

Void Geometry Driven Spin Crossover in Zeolite-Encapsulated Cobalt Tris(bipyridyl) Complex Ion

The cobalt(II) tris(bipyridyl) complex ion encapsulated in zeolite-Y supercages exhibits a thermally driven interconversion between a low-spin and a high-spin state-a phenomenon not observed for this ion either in solid state or in solution. From a comparative study of the magnetism and optical spectroscopy of the encapsulated and unencapsulated complex ion, supported by molecular modeling, such spin behavior is shown to be intramolecular in origin. In the unencapsulated or free state, the [Co(bipy)(3)](2+) ion exhibits a marked trigonal prismatic distortion, but on encapsulation, the topology of the supercage forces it to adopt a near-octahedral geometry. An analysis using the angular overlap ligand field model with spectroscopically derived parameters shows that the geometry does indeed give rise to a low-spin ground state, and suggests a possible scenario for the spin state interconversion.

Bounds on isoperimetric values of trees

Let G = (V, E) be a finite, simple and undirected graph. For S subset of V, let delta(S, G) = {(u, v) is an element of E : u is an element of S and v is an element of V - S} be the edge boundary of S. Given an integer i, 1 <= i <= vertical bar V vertical bar, let the edge isoperimetric value of G at i be defined as b(e)(i, G) = min(S subset of V:vertical bar S vertical bar=i)vertical bar delta(S, G)vertical bar. The edge isoperimetric peak of G is defined as b(e)(G) = max(1 <= j <=vertical bar V vertical bar)b(e)(j, G). Let b(v)(G) denote the vertex isoperimetric peak defined in a corresponding way. The problem of determining a lower bound for the vertex isoperimetric peak in complete t-ary trees was recently considered in [Y. Otachi, K. Yamazaki, A lower bound for the vertex boundary-width of complete k-ary trees, Discrete Mathematics, in press (doi: 10.1016/j.disc.2007.05.014)]. In this paper we provide bounds which improve those in the above cited paper. Our results can be generalized to arbitrary (rooted) trees. The depth d of a tree is the number of nodes on the longest path starting from the root and ending at a leaf. In this paper we show that for a complete binary tree of depth d (denoted as T-d(2)), c(1)d <= b(e) (T-d(2)) <= d and c(2)d <= b(v)(T-d(2)) <= d where c(1), c(2) are constants. For a complete t-ary tree of depth d (denoted as T-d(t)) and d >= c log t where c is a constant, we show that c(1)root td <= b(e)(T-d(t)) <= td and c(2)d/root t <= b(v) (T-d(t)) <= d where c(1), c(2) are constants. At the heart of our proof we have the following theorem which works for an arbitrary rooted tree and not just for a complete t-ary tree. Let T = (V, E, r) be a finite, connected and rooted tree - the root being the vertex r. Define a weight function w : V -> N where the weight w(u) of a vertex u is the number of its successors (including itself) and let the weight index eta(T) be defined as the number of distinct weights in the tree, i.e eta(T) vertical bar{w(u) : u is an element of V}vertical bar. For a positive integer k, let l(k) = vertical bar{i is an element of N : 1 <= i <= vertical bar V vertical bar, b(e)(i, G) <= k}vertical bar. We show that l(k) <= 2(2 eta+k k)

Local structure of Sr2FeMoxW1-xO6 double perovskites across the composition-driven metal to insulator transition

Sr2FeMoO6 oxides exhibit a half-metallic ferromagnetic (HM-FM) ground state and peculiar magnetic and magnetotransport properties, which are interesting for applications in the emerging field of spintronics and attractive for fundamental research in the field of heavily correlated electron systems. Sr2FeWO6 is an insulator with an antiferromagnetic (I-AFM) ground state. The solid solutions Sr2FeMoxW1-xO6 also have peculiar properties-W doping enhances chemical order which allows stabilization of the HM-FM state; as the W content exceeds a certain value a metal to insulator transition (MIT) occurs. The role of W in determining the physical properties of Sr2FeMoxW1-xO6 systems has been a matter of intense investigation. This work deals with the problem of the structural and electronic changes related to the MIT from a local perspective by means of x-ray absorption spectroscopy (XAS). This technique allows one to probe in detail the local structure and electronic modifications around selected absorber ions (W, Mo, Fe and Sr in our case). The results of XAS analysis in the whole composition range (0 <= x <= 1), in the near edge (XANES) and extended (EXAFS) regions, demonstrate an abrupt change of the local structure around the Fe and Mo sites at the critical composition, x(c). This change represents the microstructural counterpart associated with the MIT. Conversely, the local structure and electronic configuration of W ions remain unaltered in the whole composition range, suggesting indirect participation of W in the MIT.

Competing A-site and B-site driven ferroelectric instabilities in the (1-x)PbTiO3-(x)BiAlO3 system

The system (1-x)PbTiO3-(x)BiAlO3 has been investigated with regard to its solid solubility, crystal structure, microstructure, and ferroelectric transition. The unit cell volume and the tetragonality exhibit anomalous behavior near x=0.10. The Curie point (T-C) of PbTiO3 was however found to be nearly unchanged. The study seems to suggest that the decrease in the stability of the ferroelectric state due to dilution of the Ti-sublattice by smaller sized Al+3 ions is compensated by the increase in the ferroelectric stability by the Bi+3 ions.

Design of high-order filters with restricted component values

The feasibility of realising a high-order LC filter with a small set of different capacitor values, without sacrificing the frequency response specifications, is indicated. This idea could be conveniently adopted in other filter structures also—for example the FDNR transformed filter realisations.

Calorimetric values of composite solid propellants

The calorimetric values of composite solid propellant based on polystyrene, polyphenolformaldehyde, poly(vinyl chloride) and carboxy-terminated polybutadiene were determined using combustion calorimetry in order to assess the uncertainities in their measurements. The dependence of the calorimetric values on various propellant composition was obtained. The stoichiometry of oxidizer and fuel in the propellant for complete combustion obtained experimentally were compared with the theoretical stoichiometry calculated based on the oxidizer decomposition.

Precise Values of the "C Component Vicinal Coupling Constants for Calculating Side-Chain Conformations in Amino Acids

From consideration of 'H-lH vicinal coupling constants and '"G'H long-range coupling constants in a series of amino acid derivatives, the precise values of uC component vicinal coupling constants have been calculated for the three minimum energy staggered rotamers for the C(or)H-C(P)H, side-chains of amino acids.

New method of calculating calorific values from elemental compositions of fossil fuels

A new method of calculating the calorific values of fossil fuels from their chemical composition has been developed, based on the concept that heats of reaction of stoichiometric fuel-oxidizer systems are rectilinearly related with the total oxidizing or reducing valancies of the mixture. The calorific value of fossil fuels has been shown to be directly related to the net reducing valencies of the fuel. The proposed method is simple and compares favourably with the other prominent methods reported in the literature.

Stress and strain-driven algorithmic formulations for finite strain viscoplasticity for hybrid and standard finite elements

This work deals with the formulation and implementation of finite deformation viscoplasticity within the framework of stress-based hybrid finite element methods. Hybrid elements, which are based on a two-field variational formulation, are much less susceptible to locking than conventional displacement-based elements. The conventional return-mapping scheme cannot be used in the context of hybrid stress methods since the stress is known, and the strain and the internal plastic variables have to be recovered using this known stress field.We discuss the formulation and implementation of the consistent tangent tensor, and the return-mapping algorithm within the context of the hybrid method. We demonstrate the efficacy of the algorithm on a wide range of problems.