Generation of representation models for complex systems using Lagrangian functions

In this article, a new methodology is presented to obtain representation models for a priori relation z = u(x1, x2, . . . ,xn) (1), with a known an experimental dataset zi; x1i ; x2i ; x3i ; . . . ; xni i=1;2;...;p· In this methodology, a potential energy is initially defined over each possible model for the relationship (1), what allows the application of the Lagrangian mechanics to the derived system. The solution of the Euler–Lagrange in this system allows obtaining the optimal solution according to the minimal action principle. The defined Lagrangian, corresponds to a continuous medium, where a n-dimensional finite elements model has been applied, so it is possible to get a solution for the problem solving a compatible and determined linear symmetric equation system. The computational implementation of the methodology has resulted in an improvement in the process of get representation models obtained and published previously by the authors.

Delamination Arrestment in Bonded-Bolted Composite Structures by Fasteners

Thesis (Ph.D.)--University of Washington, 2016-06

Hydrogen economy options for Australia

Global concerns over the effects of current carbon dioxide (CO2) emissions have lead to extensive research on the use of hydrogen as a potential energy carrier for a lower emissions society. Hydrogen can be produced from both fossil and renewable energy sources. The hydrogen economy, in which hydrogen will be a carrier of energy from renewable sources, is a long-term development and any increasing demand for hydrogen will probably be covered initially from fossil sources. Technologies for hydrogen generation from renewable energies are being explored, whereas technologies for hydrogen production from fossil fuels have to a certain extent reached maturity. This paper addresses the major hydrogen generation processes and utilisation technology (fuel cells) currently available for the move from a fossil fuelsbased economy to a hydrogen economy. In particular, it illustrates the applicability of different hydrogen sources using Australia as an example.

The rearrangements of naphthyinitrenes: UV/Vis and IR spectra of azirines, cyclic ketenimines, and cyclic nitrile ylides

Ar matrix photolysis of 1- and 2-naphthyl azides 3 and 4 at 313 nm initially affords the singlet naphthyl nitrenes, (1)1 and (1)2. Relaxation to the corresponding lower energy, persistent triplet nitrenes (3)1 and (3)2 competes with cyclization to the azirines 15 and 18, which can also be formed photochemically from the triplet nitrenes. On prolonged irradiation, the azirines can be converted to the seven-membered cyclic ketenimines 10 and 13, respectively, as described earlier by Dunkin and Thomson. However, instead of the o-quinoid ketenimines 16 and 19, which are the expected primary ring-opening products of azirines 15 and 18, respectively, we observed their novel bond-shift isomers 17 and 20, which may be formally regarded as cyclic nitrile ylides. The existence of such ylidic heterocumulenes has been predicted previously, but this work provides the first experimental observation of such species. The factors which are responsible for the special stability of the ylidic species 17 and 20 are discussed.

Transition-state-theory calculations for reactions of O(3P) with halogenated olefins

The potential energy surfaces for the reactions of atomic oxygen in its ground electronic state, O(P-3), with the olefins: CF2=CCl2 and CF2=CF - CF3, have been characterized using ab initio molecular orbital calculations. Geometry optimization and vibrational frequency calculations were performed for reactants, transition states and products at the MP2 and QCISD levels of theory using the 6-31G(d) basis set. This database was then used to calculate the rate constants by means of Transition-State-Theory. To obtain a better reference and to test the reliability of the activation barriers we have also carried out computations using the CCSD(T)(fc)/6-311Gdagger, MP4(SDQ)(fc)/CBSB4 and MP2(fc)/CBSB3 single point energy calculations at both of the above levels of theory, as well as with the composite CBS-RAD procedure ( P. M. Mayer, C. J. Parkinson, D. M. Smith and L. Radom, J. Chem. Phys., 1998, 108, 604) and a modi. cation of this approach, called: CBS-RAD( MP2, MP2). It was found that the kinetic parameters obtained in this work particularly with the CBS-RAD ( MP2, MP2) procedure are in reasonable agreement with the experimental values. For both reactions it is found that the channels leading to the olefin double-bond addition predominates with respect to any other reaction pathway. However, on account of the different substituents in the alkenes we have located, at all levels of theory, two transition states for each reaction. Moreover, we have found that, for the reactions studied, a correlation exists between the activation energies and the electronic structure of the transition states which can explain the influence of the substituent effect on the reactivity of the halo-olefins.

Unimolecular Rovibrational Bound and Resonance States for Large Angular Momentum: J=20 Calculations for HO2

We explore the calculation of unimolecular bound states and resonances for deep-well species at large angular momentum using a Chebychev filter diagonalization scheme incorporating doubling of the autocorrelation function as presented recently by Neumaier and Mandelshtam [Phys. Rev. Lett. 86, 5031 (2001)]. The method has been employed to compute the challenging J=20 bound and resonance states for the HO2 system. The methodology has firstly been tested for J=2 in comparison with previous calculations, and then extended to J=20 using a parallel computing strategy. The quantum J-specific unimolecular dissociation rates for HO2-> H+O-2 in the energy range from 2.114 to 2.596 eV have been reported for the first time, and comparisons with the results of Troe and co-workers [J. Chem. Phys. 113, 11019 (2000) Phys. Chem. Chem. Phys. 2, 631 (2000)] from statistical adiabatic channel method/classical trajectory calculations have been made. For most of the energies, the reported statistical adiabatic channel method/classical trajectory rate constants agree well with the average of the fluctuating quantum-mechanical rates. Near the dissociation threshold, quantum rates fluctuate more severely, but their average is still in agreement with the statistical adiabatic channel method/classical trajectory results.

21-residue peptide's dynamics at and between elementary structural transitions

Elementary conformational changes of the backbone of a 21-residue peptide A5(A3RA)3A are studied using molecular dynamics simulations in explicit water. The processes of the conformational transitions and the regimes of stationary fluctuations between them are investigated using minimal perturbations of the system. The perturbations consist of a few degrees rotation of the velocity of one of the systems' atoms and keep the system on the same energy surface. It is found that (i) the system dynamics is insignificantly changed by the perturbations in the regimes between the transitions; (ii) it is very sensitive to the perturbations just before the transitions that prevents the peptide from making the transitions; and (iii) the perturbation of any atom of the system, including distant water molecules is equally effective in preventing the transition. The latter implies strongly collective dynamics of the peptide and water during the transitions.

Investigation of gas engine driven heat pump systems

Heat pumps are becoming increasingly popular, but poor electricity generating efficiency limits the potential energy savings of electrically powered units. Thus the work reported in this thesis concerns the development of a range of gas engine driven heat pumps for industrial and commercial heating applications, which recover heat from the prime mover, normally rejected to waste. Despite the convenience of using proprietary engine heat recovery packages, investigations have highlighted the necessity to ensure the engine and the heat recovery equipment are compatible. A problem common •to all air source heat pumps is the formation of frost on the evaporator, which must be removed periodically, with the expenditure of energy, to ensure the continued operation of the plant. An original fluidised bed defrosting mechanism is proposed, which prevents the build-up of this frost, and also improves system performance. Criticisms have been levelled against the rotary sliding vane compressor, in particular the effects of lubrication, which is essential. This thesis compares the rotary sliding vane compressor with other machines, and concludes that many of these criticisms are unfounded. A confidential market survey indicates an increasing demand for heat pumps up to and including 1990, and the technical support needed to penetrate this market is presented. Such support includes the development of a range of modular gas engine driven heat pumps, and a computer aided design for the selection of the optimum units. A case study of a gas engine driven heat pump for a swimming pool application which provided valuable experience is included.

The intermediate pyrolysis and catalytic steam reforming of brewers spent grain

Brewers spent grain (BSG) is a widely available feedstock representing approximately 85% of the total by-products generated in the brewing industry. This is currently either disposed of to landfill or used as cattle feed due to its high protein content. BSG has received little or no attention as a potential energy resource, but increasing disposal costs and environmental constraints are now prompting the consideration of this. One possibility for the utilisation of BSG for energy is via intermediate pyrolysis to produce gases, vapours and chars. Intermediate pyrolysis is characterised by indirect heating in the absence of oxygen for short solids residence times of a few minutes, at temperatures of 350-450 °C. In the present work BSG has been characterised by chemical, proximate, ultimate and thermo-gravimetric analysis. Intermediate pyrolysis of BSG at 450 °C was carried out using a twin coaxial screw reactor known as Pyroformer to give yields of char 29%, 51% of bio-oil and 19% of permanent gases. The bio-oil liquid was found to separate in to an aqueous phase and organic phase. The organic phase contained viscous compounds that could age over time leading to solid tars that can present problems in CHP application. The quality of the pyrolysis vapour products before quenching can be upgraded to achieve much improved suitability as a fuel by downstream catalytic reforming. A Bench Scale batch pyrolysis reactor has then been used to pyrolyse small samples of BSG under a range of conditions of heating rate and temperature simulating the Pyroformer. A small catalytic reformer has been added downstream of the reactor in which the pyrolysis vapours can be further cracked and reformed. A commercial reforming nickel catalyst was used at 500, 750 and 850 °C at a space velocity about 10,000 L/h with and without the addition of steam. Results are presented for the properties of BSG, and the products of the pyrolysis process both with and without catalytic post-processing. Results indicate that catalytic reforming produced a significant increase in permanent gases mainly (H2 and CO) with H2 content exceeding 50 vol% at higher reforming temperatures. Bio-oil yield decreased significantly as reforming temperature increased with char remaining the same as pyrolysis condition remained unchanged. The process shows an increase in heating value for the product gas ranging between 10.8-25.2 MJ/m as reforming temperature increased. © 2012 Elsevier B.V. All rights reserved.

Formation of Polycyclic Aromatic Hydrocarbons and Nitrogen Containing Polycyclic Aromatic Compounds in Titan's Atmosphere, the Interstellar Medium and Combustion

Several different mechanisms leading to the formation of (substituted) naphthalene and azanaphthalenes were examined using theoretical quantum chemical calculations. As a result, a series of novel synthetic routes to Polycyclic Aromatic Hydrocarbons (PAHs) and Nitrogen Containing Polycyclic Aromatic Compounds (N-PACs) have been proposed. On Earth, these aromatic compounds originate from incomplete combustion and are released into our environment, where they are known to be major pollutants, often with carcinogenic properties. In the atmosphere of a Saturn's moon Titan, these PAH and N-PACs are believed to play a critical role in organic haze formation, as well as acting as chemical precursors to biologically relevant molecules. The theoretical calculations were performed by employing the ab initio G3(MP2,CC)/B3LYP/6-311G** method to effectively probe the Potential Energy Surfaces (PES) relevant to the PAH and N-PAC formation. Following the construction of the PES, Rice-Ramsperger-Kassel-Markus (RRKM) theory was used to evaluate all unimolecular rate constants as a function of collision energy under single-collision conditions. Branching ratios were then evaluated by solving phenomenological rate expressions for the various product concentrations. The most viable pathways to PAH and N-PAC formation were found to be those where the initial attack by the ethynyl (C2H) or cyano (CN) radical toward a unsaturated hydrocarbon molecule led to the formation of an intermediate which could not effectively lose a hydrogen atom. It is not until ring cyclization has occurred, that hydrogen elimination leads to a closed shell product. By quenching the possibility of the initial hydrogen atom elimination, one of the most competitive processes preventing the PAH or N-PAC formation was avoided, and the PAH or N-PAC formation was allowed to proceed. It is concluded that these considerations should be taken into account when attempting to explore any other potential routes towards aromatic compounds in cold environments, such as on Titan or in the interstellar medium.

Contribuição ao estudo de treliças não lineares considerando efeitos térmicos

This work presents the numerical analysis of nonlinear trusses summited to thermomechanical actions with Finite Element Method (FEM). The proposed formulation is so-called positional FEM and it is based on the minimum potential energy theorem written according to nodal positions, instead of displacements. The study herein presented considers the effects of geometric and material nonlinearities. Related to dynamic problems, a comparison between different time integration algorithms is performed. The formulation is extended to impact problems between trusses and rigid wall, where the nodal positions are constrained considering nullpenetration condition. In addition, it is presented a thermodynamically consistent formulation, based on the first and second law of thermodynamics and the Helmholtz free-energy for analyzing dynamic problems of truss structures with thermoelastic and thermoplastic behavior. The numerical results of the proposed formulation are compared with examples found in the literature.

Generalized atoms in molecules approach (GAIM approach)

This thesis involves two parts. The first is a new-proposed theoretical approach called generalized atoms in molecules (GAIM). The second is a computational study on the deamination reaction of adenine with OH⁻/nH₂O (n=0, 1, 2, 3) and 3H₂O. The GAIM approach aims to solve the energy of each atom variationally in the first step and then to build the energy of a molecule from each atom. Thus the energy of a diatomic molecule (A-B) is formulated as a sum of its atomic energies, EA and EB. Each of these atomic energies is expressed as, EA = Hᴬ + Vₑₑᴬᴬ + 1/2Vₑₑᴬ<>ᴮ EB = Hᴮ + Vₑₑᴮᴮ + 1/2Vₑₑᴬ<>ᴮ where; Hᴬ and Hᴮ are the kinetic and nuclear attraction energy of electrons of atoms A and B, respectively; Vₑₑᴬᴬ and Vₑₑᴮᴮ are the interaction energy between the electrons on atoms A and B, respectively; and Vₑₑᴬ<>ᴮ is the interaction energy between the electrons of atom A with the electrons of atom B. The energy of the molecule is then minimized subject to the following constraint, |ρA(r)dr + |ρB(r)dr = N where ρA(r) and ρB(r) are the electron densities of atoms A and B, respectively, and N is the number of electrons. The initial testing of the performance of GAIM was done through calculating dissociation curves for H₂, LiH, Li₂, BH, HF, HCl, N₂, F₂, and Cl₂. The numerical results show that GAIM performs very well with H₂, LiH, Li₂, BH, HF, and HCl. GAIM shows convergence problems with N₂, F₂, and Cl₂ due to difficulties in reordering the degenerate atomic orbitals Pₓ, Py, and Pz in N, F, and Cl atoms. Further work for the development of GAIM is required. Deamination of adenine results in one of several forms of premutagenic lesions occurring in DNA. In this thesis, mechanisms for the deamination reaction of adenine with OH⁻/nH₂O, (n = 0, 1, 2, 3) and 3H₂O were investigated. HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), and B3LYP/6-31+G(d) levels of theory were employed to optimize all the geometries. Energies were calculated at the G3MP2B3 and CBS-QB3 levels of theory. The effect of solvent (water) was computed using the polarizable continuum model (PCM). Intrinsic reaction coordinate (IRC) calculations were performed for all transition states. Five pathways were investigated for the deamination reaction of adenine with OH⁻/nH₂O and 3H₂O. The first four pathways (A-D) begin with by deprotonation at the amino group of adenine by OH⁻, while pathway E is initiated by tautomerization of adenine. For all pathways, the next two steps involve the formation of a tetrahedral intermediate followed by dissociation to yield products via a 1,3-hydrogen shift. Deamination with a single OH⁻ has a high activation barrier (190 kJ mol⁻¹ using G3MP2B3 level) for the rate-determining step. Addition of one water molecule reduces this barrier by 68 kJ mol⁻¹ calculated at G3MP2B3 level. Adding more water molecules decreases the overall activation energy of the reaction, but the effect becomes smaller with each additional water molecule. The most plausible mechanism is pathway E, the deamination reaction of adenine with 3H₂O, which has an overall G3MP2B3 activation energy of 139 and 137 kJ mol⁻¹ in the gas phase and PCM, respectively. This barrier is lower than that for the deamination with OH⁻/3H₂O by 6 and 2 kJ mol⁻¹ in the gas phase and PCM, respectively.

Age determination of sediment core NP90-52 (Table 1)

Data on glacial erosion have been compiled and synthesised using a wide range of sediment budget and sediment yield studies from the Svalbard-Barents Sea region. The data include studies ranging in timescale from 1 to 10**6 yr, and in size of drainage basin from 101 to 105 km**2. They show a clear dependence of sediment yield on the mode of glacierization. Polar glaciers erode at rates comparable to those found in Arctic fluvial basins, or about 40 t/km**-2/ yr or 0.02 mm/yr. In contrast, rates of erosion by polythermal glaciers are 800-1000 t/km**2/ yr (or ca 0.3-0.4 mm/yr), while rates from fast-flowing glaciers are slightly more than twice this: 2100 t/km**2/yr (or 1 mm/yr). Similar rates are also found for large glacierized basins like those in the southwestern parts of the Barents Sea. In contrast to the situation in fluvial basins, in which sediment yield typically decreases with increasing basin size, the tendency in glacierized basins is for erosion to be independent of basin size. In studies of sediment yield from glaciers it is sometimes difficult to distinguish between material actually dislodged from the bedrock by glaciers and material dislodged by other processes in interglacial times and simply transported to a depocenter by a glacier. Our data suggest that pulses of sediment resulting from advance of a glacier over previously-dislodged material last on the order of 10**3 yr, and result in inferred erosion rates that are approximately 25% higher than long-term average rates of glacial erosion. The maximum sediment fluxes from the large Storfjorden and Bear Island drainage basins occurred in mid-Pleistocene. The onset of this period of high sediment yield coincided with the onset of the 100 kyr glacial cycle. We presume that this was the beginning of a period of increased glacial activity, but one in which glaciers still advanced and retreated frequently. During the last two to four 100 kyr cycles, however, sediment yields appear to have decreased. This decrease may be the result of the submergence of the Barents Sea. Glacier erosion would be much higher for a subaerial Barents Sea setting than it would be for a present day subsea Barents Sea. A classical question in Quaternary Geology is whether glaciers are more erosive than rivers. We surmise that if factors such as the lithology and the available potential energy (mgh) of the precipitation falling at a given altitude, whether in liquid or solid form, are held constant, then glaciers are vastly more effective agents of erosion than rivers.

Nonholonomically Constrained Dynamics and Optimization of Rolling Isolation Systems

Rolling Isolation Systems provide a simple and effective means for protecting components from horizontal floor vibrations. In these systems a platform rolls on four steel balls which, in turn, rest within shallow bowls. The trajectories of the balls is uniquely determined by the horizontal and rotational velocity components of the rolling platform, and thus provides nonholonomic constraints. In general, the bowls are not parabolic, so the potential energy function of this system is not quadratic. This thesis presents the application of Gauss's Principle of Least Constraint to the modeling of rolling isolation platforms. The equations of motion are described in terms of a redundant set of constrained coordinates. Coordinate accelerations are uniquely determined at any point in time via Gauss's Principle by solving a linearly constrained quadratic minimization. In the absence of any modeled damping, the equations of motion conserve energy. This mathematical model is then used to find the bowl profile that minimizes response acceleration subject to displacement constraint.

Dynamics of the Disk-Pendulum Coupled System With Vertical Excitation

This paper investigates the static and dynamic characteristics of the semi-elliptical rocking disk on which a pendulum pinned. This coupled system’s response is also analyzed analytically and numerically when a vertical harmonic excitation is applied to the bottom of the rocking disk. Lagrange’s Equation is used to derive the motion equations of the disk-pendulum coupled system. The second derivative test for the system’s potential energy shows how the location of the pendulum’s pivotal point affects the number and stability of equilibria, and the change of location presents different bifurcation diagrams for different geometries of the rocking disk. For both vertically excited and unforced cases, the coupled system shows chaos easily, but the proper chosen parameters can still help the system reach and keep the steady state. For the steady state of the vertically excited rocking disk without a pendulum, the variation of the excitation’s amplitude and frequency result in the hysteresis for the amplitude of the response. When a pendulum is pinned on the rocking disk, three major categories of steady states are presently in the numerical way.