A new parallel overlapped domain decomposition method for nonlinear dynamic finite element analysis

In this paper a new parallel algorithm for nonlinear transient dynamic analysis of large structures has been presented. An unconditionally stable Newmark-beta method (constant average acceleration technique) has been employed for time integration. The proposed parallel algorithm has been devised within the broad framework of domain decomposition techniques. However, unlike most of the existing parallel algorithms (devised for structural dynamic applications) which are basically derived using nonoverlapped domains, the proposed algorithm uses overlapped domains. The parallel overlapped domain decomposition algorithm proposed in this paper has been formulated by splitting the mass, damping and stiffness matrices arises out of finite element discretisation of a given structure. A predictor-corrector scheme has been formulated for iteratively improving the solution in each step. A computer program based on the proposed algorithm has been developed and implemented with message passing interface as software development environment. PARAM-10000 MIMD parallel computer has been used to evaluate the performances. Numerical experiments have been conducted to validate as well as to evaluate the performance of the proposed parallel algorithm. Comparisons have been made with the conventional nonoverlapped domain decomposition algorithms. Numerical studies indicate that the proposed algorithm is superior in performance to the conventional domain decomposition algorithms. (C) 2003 Elsevier Ltd. All rights reserved.

Decomposition temperatures of Cu2Ln2O5 (Ln = Tb, Dy, Ho, Er, Tm, Yb, and Lu) compounds

The equilibrium decomposition temperatures of Cu2Ln2O5 (Ln = Tb, Dy, Ho, Er, Tm, Yb, and Lu) compounds have been measured using a combined DTA-TGA apparatus under a flowing Ar + O2 gas mixture, in which the partial pressure of oxygen was controlled at 5.0 × 103 Pa. The Cu2Ln2O5 compounds yield Ln2O3 and Cu2O on decomposition. The decomposition temperature increases monotonically with the atomic number of the lanthanide element. This suggests that the stability of the Cu2Ln2O5 compounds with respect to the component binary oxides increases with decreasing radius of the Ln3+ ion.

Kinetic decomposition of Ni2SiO4 in oxygen potential gradients

Nickel orthosilicate (Ni2SiO4) has been found to decompose into its component binary oxides in oxygen potential gradients at 1373 K. Nickel oxide was formed at the high oxygen potential boundary, while silica was detected at the low oxygen potential side. Significant porosity and fissures were observed near the Ni2SiO4/SiO2 interface and the SiO2 layer. The critical oxygen partial pressure ratio required for decomposition varied from 1.63 to 2.15 as the oxygen pressures were altered from 1.01 ⊠ 105 to 2.7X 10−4 Pa, well above the dissociation pressure of Ni2SiO4. Platinum markers placed at the boundaries of the Ni2SiO4 sample indicated growth of NiO at the higher oxygen potential boundary, without any apparent transport of material to the low oxygen potential side. However, significant movement of the bulk Ni2SiO4 crystal with respect to the marker was not observed. The decomposition of the silicate occurs due to the unequal rates of transport of Ni and Si. The critical oxygen partial pressure ratio required for decomposition is related both to the thermodynamic stability of Ni2SiO4 with respect to component oxides and the ratio of diffusivities of nickel and silicon. Kinetic decomposition of multicomponent oxides, first discovered by Schmalzried, Laqua, and co-workers [H. Schmalzried, W. Laqua, and P. L. Lin, Z. Natur Forsch. Teil A 34, 192 (1979); H. Schmalzried and W. Laqua, Oxid. Met. 15, 339 (1981); W. Laqua and H. Schmalzried, Chemical Metallurgy—A Tribute to Carl Wagner (Metallurgical Society of the AIME, New York, 1981), p. 29] has important consequences for their use at high temperatures and in geochemistry.

An investigation of carbon nanotubes obtained from the decomposition of methane over reduced Mg1− xM xAl2O4 spinel catalysts

Carbon nanotubes produced by the treatment of Mg1−xMxAl2O4 (M = Fe, Co, or Ni; x = 0.1, 0.2, 0.3, or 0.4) spinels with an H2–CH4 mixture at 1070 °C have been investigated systematically. The grains of the oxide-metal composite particles are uniformly covered by a weblike network of carbon nanotube bundles, several tens of micrometers long, made up of single-wall nanotubes with a diameter close to 4 nm. Only the smallest metal particles (<5 nm) are involved in the formation of the nanotubes. A macroscopic characterization method involving surface area measurements and chemical analysis has been developed in order to compare the different nanotube specimens. An increase in the transition metal content of the catalyst yields more carbon nanotubes (up to a metal content of 10.0 wt% or x = 0.3), but causes a decrease in carbon quality. The best compromise is to use 6.7 wt% of metal (x = 0.2) in the catalyst. Co gives superior results with respect to both the quantity and quality of the nanotubes. In the case of Fe, the quality is notably hampered by the formation of Fe3C particles.

Thermal decomposition of haloethanols: single pulse shock tube and ab initio studies

This paper reports single pulse shock tube and ab initio studies on thermal decomposition of 2-fluoro and 2-chloroethanol at T=1000–1200 K. Both molecules have HX (X = F/Cl) and H2O molecular elimination channels. The CH3CHO formed by HX elimination is chemically active and undergoes secondary decomposition resulting in the formation of CH4, C2H6, and C2H4. A detailed kinetic simulation indicates that the formation of C2H4 could not be quantitatively explained as arising exclusively from secondary CH3CHO decomposition. Contributions from primary radical processes need to be considered to explain C2H4 quantitatively. Ab initio calculations on HX and H2O elimination reactions from the haloethanols at HF, MP2, and DFT levels with various basis sets up to 6/311++G**are reported. It is pointed out that due to strong correlations between A and Eα, comparison of these two parameters between experimental and theoretical results could be misleading.

Mechanochemical decomposition of Gd3Fe5O12 garnet phase

Gadolinium iron garnet was milled in a high energy ball mill to study its magnetic properties in the nanocrystalline regime. XRD reveals the decomposition of the garnet phase into Gd-orthoferrite and Gd2O3 on milling. The variation of saturation magnetization and coercivity with milling is attributed to a possible shift in the compensation temperature on grain size reduction and an increase in the orthoferrite content. The Mössbauer spectrum at 16 K is characteristic of the magnetically ordered state corresponding to GdIG, GdFeO3 and α-Fe2O3 whereas at room temperature it is a superparamagnetic doublet.

Effect of convection and viscoelasticity on spinodal decomposition in fluids

Current analytical work on the effect of convection and viscoelasticity on the early and late stages of spinodal decomposition is briefly described. In the early stages, the effect of viscoelastic stresses was analysed using a simple Maxwell model for the stress, which was incorporated in the Langevin equation for the momentum field. The viscoelastic stresses are found to enhance the rate of decomposition. In the late stages, the pattern formed depends on the relative composition of the two species. Droplet spinodal decomposition occurs when the concentration of one of the species is small. Convective transport does not have a significant effect on the growth of a single droplet, but it does result in an attractive interaction between non - Brownian droplets which could lead to coalescence. The effect of convective transport for the growth of random interfaces in a near symmetric quench was analysed using an 'area distribution function', which gives the distribution of surface area of the interface in curvature space. It was found that the curvature of the interface decreases proportional to t in the late stages of spinodal decomposition, and the surface area also decreases proportional to t.