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.

Distributed nonlinear optimization algorithms for general network problems

There are a number of large networks which occur in many problems dealing with the flow of power, communication signals, water, gas, transportable goods, etc. Both design and planning of these networks involve optimization problems. The first part of this paper introduces the common characteristics of a nonlinear network (the network may be linear, the objective function may be non linear, or both may be nonlinear). The second part develops a mathematical model trying to put together some important constraints based on the abstraction for a general network. The third part deals with solution procedures; it converts the network to a matrix based system of equations, gives the characteristics of the matrix and suggests two solution procedures, one of them being a new one. The fourth part handles spatially distributed networks and evolves a number of decomposition techniques so that we can solve the problem with the help of a distributed computer system. Algorithms for parallel processors and spatially distributed systems have been described.There are a number of common features that pertain to networks. A network consists of a set of nodes and arcs. In addition at every node, there is a possibility of an input (like power, water, message, goods etc) or an output or none. Normally, the network equations describe the flows amoungst nodes through the arcs. These network equations couple variables associated with nodes. Invariably, variables pertaining to arcs are constants; the result required will be flows through the arcs. To solve the normal base problem, we are given input flows at nodes, output flows at nodes and certain physical constraints on other variables at nodes and we should find out the flows through the network (variables at nodes will be referred to as across variables).The optimization problem involves in selecting inputs at nodes so as to optimise an objective function; the objective may be a cost function based on the inputs to be minimised or a loss function or an efficiency function. The above mathematical model can be solved using Lagrange Multiplier technique since the equalities are strong compared to inequalities. The Lagrange multiplier technique divides the solution procedure into two stages per iteration. Stage one calculates the problem variables % and stage two the multipliers lambda. It is shown that the Jacobian matrix used in stage one (for solving a nonlinear system of necessary conditions) occurs in the stage two also.A second solution procedure has also been imbedded into the first one. This is called total residue approach. It changes the equality constraints so that we can get faster convergence of the iterations.Both solution procedures are found to coverge in 3 to 7 iterations for a sample network.The availability of distributed computer systems — both LAN and WAN — suggest the need for algorithms to solve the optimization problems. Two types of algorithms have been proposed — one based on the physics of the network and the other on the property of the Jacobian matrix. Three algorithms have been deviced, one of them for the local area case. These algorithms are called as regional distributed algorithm, hierarchical regional distributed algorithm (both using the physics properties of the network), and locally distributed algorithm (a multiprocessor based approach with a local area network configuration). The approach used was to define an algorithm that is faster and uses minimum communications. These algorithms are found to converge at the same rate as the non distributed (unitary) case.

Evaluation of Precipitation Retrievals From Orbital Data Products of TRMM Over a Subtropical Basin in India

The spatial error structure of daily precipitation derived from the latest version 7 (v7) tropical rainfall measuring mission (TRMM) level 2 data products are studied through comparison with the Asian precipitation highly resolved observational data integration toward evaluation of the water resources (APHRODITE) data over a subtropical region of the Indian subcontinent for the seasonal rainfall over 6 years from June 2002 to September 2007. The data products examined include v7 data from the TRMM radiometer Microwave Imager (TMI) and radar precipitation radar (PR), namely, 2A12, 2A25, and 2B31 (combined data from PR and TMI). The spatial distribution of uncertainty from these data products were quantified based on performance metrics derived from the contingency table. For the seasonal daily precipitation over a subtropical basin in India, the data product of 2A12 showed greater skill in detecting and quantifying the volume of rainfall when compared with the 2A25 and 2B31 data products. Error characterization using various error models revealed that random errors from multiplicative error models were homoscedastic and that they better represented rainfall estimates from 2A12 algorithm. Error decomposition techniques performed to disentangle systematic and random errors verify that the multiplicative error model representing rainfall from 2A12 algorithm successfully estimated a greater percentage of systematic error than 2A25 or 2B31 algorithms. Results verify that although the radiometer derived 2A12 rainfall data is known to suffer from many sources of uncertainties, spatial analysis over the case study region of India testifies that the 2A12 rainfall estimates are in a very good agreement with the reference estimates for the data period considered.

Silica from rice husk through thermal decomposition

The thermal decomposition characteristics of rice husk have been investigated by dynamic thermoanalytical techniques: DTA, TG, DTG and isothermal heating. The observed thermal behaviour is explained on the basis of a superposition of the decomposition of cellulose and lignin, which are the major organic constituents of rice husk. Morphological features of silica in husk as well as the ash are examined by scanning electron microscopy. Silica in the residual ash has been characterised by X-ray diffraction and infrared spectroscopy. Controlled thermal decomposition of rice husk has been shown to be a convenient method for the liberation of silica.

Thermal decomposition and combustion of ammine: AP pellets

Thermal decomposition and combustion of lithium perchlorate ammine:ammonium perchlorate (LPA:AP) and magnesium perchlorate ammine:ammonium perchlorate (MPA:AP) pellets have been studied using DTA, TG, and strand burner techniques. The DTA results of the ammine:AP pellets show that the addition of ammines lowers the ignition temperature of AP. However, isothermal TG of the ammine:AP pellets show that in the case of LPA:AP pellets the extent of decomposition increases with the increase in the concentration of LPA; whereas in the case of MPA:AP pellets the extent of decomposition decreases with the increase in the concentration of MPA. Similarly, LPA:AP pellets show higher burning rates compared to AP pellets. On the other hand, MPA:AP pellets show lower burning rates compared to AP pellets. Increasing the concentration of MPA in MPA:AP pellets completely suppresses the combustion. These results are explained on the basis of the thermal characteristics of the additives and their decomposition products.

Formation of Pb(Ti,Zr)O3, from thermal decomposition of oxalate precursors

Lead zir conyl oxalate hexahydrate (LZO) and lead titanyl zirconyl oxalate hydrate (LTZO) are prepared and characterized. Their thermal decompositions have been investigated by thermoanalytical and gas analysis techniques. The decomposition in air or oxygen has three steps — dehydration, decomposition of the oxalate to a carbonate and the decomposition of carbonate to PbZrO3. In non oxidising atmosphere, partial reduction of Pb(II) to Pb(0) takes place at the oxalate decomposition step. The formation of free metallic lead affects the stoichiometry of the intermediate carbonate and yields a mixture of Pb(Ti,Zr)O3 and ZrO2 as the final products. By maintaining oxidising atmosphere and low heating rate, direct preparation of stoichiometric, crystalline Pb(Ti,Zr)O3 at 550°C is possible from the corresponding oxalate precursor.

Thermal decomposition of rare-earth trioxalatocobaltates

The thermal decomposition of rare-earth trioxalatocobaltates LnCo(C2O4)3 · x H2O, where Ln = La, Pr, Nd, has been studied in flowing atmospheres of air/oxygen, argon/ nitrogen, carbon dioxide and a vacuum. The compounds decompose through three major steps, viz. dehydration, decomposition of the oxalate to an intermediate carbonate, which further decomposes to yield rare-earth cobaltite as the final product. The formation of the final product is influenced by the surrounding gas atmosphere. Studies on the thermal decomposition of photodecomposed lanthanum trioxalatocobaltate and a mechanical mixture of lanthanum oxalate and cobalt oxalate in 1 : 2 molar ratio reveal that the decomposition behaviour of the two samples is different. The drawbacks of the decomposition scheme proposed earlier have been pointed out, and logical schemes based on results obtained by TG, DTA, DTG, supplemented by various physico-chemical techniques such as gas and chemical analyses, IR and mass spectroscopy, surface area and magnetic susceptibility measurements and X-ray powder diffraction methods, have been proposed for the decomposition in air of rare-earth trioxalatocobaltates as well as for the photoreduced lanthanum salt and a mechanical mixture of lanthanum and cobalt oxalates.

Thermal Decomposition Studies on Ammonium Perchlorate—Polybutadiene Composite Propellant

The thermal decomposition of ammonium perchlorate based solid composite propellant using carboxyl terminated polybutadiene as binder has been studied employing thermogravimetry and differential thermal analysis techniques. The thermal decomposition characteristics of the propellant have been found to be quite similar to those of pure ammonium perchlorate with activation energy, 32 Kcal/mole and 60 Kcal/mole respectively in the low and high temperature regions. The effect of the sample size and shape on the thermal decomposition has also been evaluated.

Thermal decomposition of titanyl oxalates-I : Barium titanyl oxalate

Thermal decomposition of barium titanyl oxalate tetrahydrate (BTO) has been investigated employing TGA, DTG and DTA techniques and gas and chemical analysis. The decomposition proceeds through five steps and is not affected much by the surrounding gas atmosphere. The first step which is the dehydration of the tetrahydrate is followed by a low-temperature decomposition of the oxalate groups. In the temperature range 190–250°C half a mole of carbon monoxide is evolved with the formation of a transient intermediate containing both oxalate and carbonate groups. The oxalate groups are completely destroyed in the range 250–450°C, resulting in the formation of a carbonate which retains free carbon dioxide in the matrix. The trapped carbon dioxide is released in the temperature range of 460–600°C. The final decomposition of the carbonate takes place between 600–750°C and yields barium titanate. The i.r. spectra, surface area measurements and X-ray, powder diffraction data support entrapment of carbon dioxide in the matrix.

Thermal Decomposition and Combustion Studies on Potassium Perchlorat e/Polystyrene Propellant Systems

STUDIES on potassium perchlorate/polystyrene (KP/PS) propellant systems have been carried out by using such techniques as thermogravimetry (TG), differential thermal analysis (DTA), and mass spectrometry (MS). It has been found that the thermal decomposition (TD) behavior of the KP/PS propellant is similar to that of the AP/PS propellant studied earlier.! It has also been observed that the TD of KP in the melt has a correlation with the burning rate (r) of KP/PS propellant at atmospheric pressure.

Thermal decomposition studies on copolyurethanes based on hydroxyl terminated polybutadiene and poly(12-hydroxy stearic acid-co-TMP) ester polyol

Thermal degradation of copolyurethanes based on hydroxyl terminated polybutadiene (HTPB) and poly(12-hydroxy stearic acid-co-TMP) ester polyol (PEP) with varying compositions has been studied by thermo-gravimetric and pyrolysis-GC techniques. The copolyurethanes were found to decompose in multiple stages and the kinetic parameters were found to be dependent on the method of their evaluation. The activation energy for the initial stage of decomposition was found to increase, and for the main stage decreases with the increase in PEP content. The pyrolysis-GC studies on the ammonium perchlorate filled copolyurethanes (solid propellants) showed that the major products during the pyrolysis were C-2, C-3 hydrocarbons and butadiene. The amount of C-2 fraction in the pyrolyslate increased with solid loading, as well as with the HTPB content in the copolyurethanes. A linear relationship apparently exists between the amount of C-2 fraction and the burn rates of the solid propellants. (C) 2000 Elsevier Science Ltd. All rights reserved.

Thermal decomposition of zirconyl oxalates I. Barium zirconyl oxalate

Conditions for the preparation of stoichiometric barium zirconyl oxalate heptahydrate (BZO) have been standardized. The thermal decomposition of BZO has been investigated employing TG, DTG and DTA techniques and chemical and gas analysis. The decomposition proceeds through four steps and is not affected much by the surrounding gas atmosphere. Both dehydration and oxalate decomposition take place in two steps. The formation of a transient intermediate containing both oxalate and carbonate groups is inferred. The decomposition of oxalate groups results in a carbonate of composition Ba2Zr2OsCO3, which decomposes between 600 and 800 ~ and yields barium zirconate. Chemical analysis, IR spectra and X-ray powder diffraction data support the identity of the intermediate as a separate entity.

Thermal decomposition of zirconyl oxalates I. Barium zirconyl oxalate

Conditions for the preparation of stoichiometric barium zirconyl oxalate heptahydrate (BZO) have been standardized. The thermal decomposition of BZO has been investigated employing TG, DTG and DTA techniques and chemical and gas analysis. The decomposition proceeds through four steps and is not affected much by the surrounding gas atmosphere. Both dehydration and oxalate decomposition take place in two steps. The formation of a transient intermediate containing both oxalate and carbonate groups is inferred. The decomposition of oxalate groups results in a carbonate of composition Ba2Zr2O5CO3, which decomposes between 600 and 800° and yields barium zirconate. Chemical analysis, IR spectra and X-ray powder diffraction data support the identity of the intermediate as a separate entity.Die Bedingungen für die Herstellung von stöchiometrischem Barium-zirconyl-oxalat Heptahydrat (BZO) wurden standardisiert. Die thermische Zersetzung von BZO wurde unter Einsatz der TG-, DTG- und DTA, sowie der chemischen und Gasanalyse untersucht. Die Zersetzung verläuft über vier Stufen und wird von der umgebenden Gasathmosphäre nicht besonders beeinflusst. Sowohl die Dehydratisierung als auch die Oxalatzersetzung erfolgt in zwei Stufen. Die Bildung einer intermediären Übergangsverbindung mit sowohl Oxalat- als auch Carbonatgruppen wirken hierbei mit. Die Zersetzung der Oxalatgruppen ergibt ein Carbonat der Zusammensetzung Ba2Zr2O5CO3, das zwischen 600 und 800° zersetzt wird und Bariumzirconat ergibt. Die Angaben der chemischen Analyse, der IR-Spekren und der Röntgen-Pulver-Diffraktion unterstützen die Identität der Intermediärverbindung als eine separate Einheit.On a standardisé les conditions de préparation de l'oxalate heptahydraté de zirconyle et de baryum (BZO) stoechiométrique. On a étudié la décomposition thermique de BZO par TG, TGD et ATD ainsi que par analyses chimiques et analyses des gaz. La décomposition a lieu en quatre étapes et n'est pas trop influencée par l'atmosphère ambiante. La déshydratation et la décomposition de l'oxalate ont lieu en deux étapes. Il se forme un composé intermédiaire de transition contenant à la fois les groupes oxalate et carbonate. La décomposition des groupes oxalate fournit un carbonate de composition Ba2Zr2O5CO3 qui se décompose entre 600 et 800° pour fournir du zirconate de baryum. L'analyse chimique, les spectres IR et la diffraction des rayons X sur poudre, apportent les preuves de l'existence d'un composé intermédiaire comme entité séparée.

Reduced Order Suboptimal Control Design for a Class of Nonlinear Distributed Parameter Systems Using POD and /spl thetas/-D Techniques

A new computational tool is presented in this paper for suboptimal control design of a class of nonlinear distributed parameter systems. First proper orthogonal decomposition based problem-oriented basis functions are designed, which are then used in a Galerkin projection to come up with a low-order lumped parameter approximation. Next, a suboptimal controller is designed using the emerging /spl thetas/-D technique for lumped parameter systems. This time domain sub-optimal control solution is then mapped back to the distributed domain using the same basis functions, which essentially leads to a closed form solution for the controller in a state feedback form. Numerical results for a real-life nonlinear temperature control problem indicate that the proposed method holds promise as a good suboptimal control design technique for distributed parameter systems.

Coherent-mode decomposition of general anisotropic Gaussian Schell-model beams

We present a complete solution to the problem of coherent-mode decomposition of the most general anisotropic Gaussian Schell-model (AGSM) beams, which constitute a ten-parameter family. Our approach is based on symmetry considerations. Concepts and techniques familiar from the context of quantum mechanics in the two-dimensional plane are used to exploit the Sp(4, R) dynamical symmetry underlying the AGSM problem. We take advantage of the fact that the symplectic group of first-order optical system acts unitarily through the metaplectic operators on the Hilbert space of wave amplitudes over the transverse plane, and, using the Iwasawa decomposition for the metaplectic operator and the classic theorem of Williamson on the normal forms of positive definite symmetric matrices under linear canonical transformations, we demonstrate the unitary equivalence of the AGSM problem to a separable problem earlier studied by Li and Wolf [Opt. Lett. 7, 256 (1982)] and Gori and Guattari [Opt. Commun. 48, 7 (1983)]. This conn ction enables one to write down, almost by inspection, the coherent-mode decomposition of the general AGSM beam. A universal feature of the eigenvalue spectrum of the AGSM family is noted.