Contact with friction using the augmented Lagrangian Method: a conditional constrained minimization problem

This work presents a formulation of the contact with friction between elastic bodies. This is a non linear problem due to unilateral constraints (inter-penetration of bodies) and friction. The solution of this problem can be found using optimization concepts, modelling the problem as a constrained minimization problem. The Finite Element Method is used to construct approximation spaces. The minimization problem has the total potential energy of the elastic bodies as the objective function, the non-inter-penetration conditions are represented by inequality constraints, and equality constraints are used to deal with the friction. Due to the presence of two friction conditions (stick and slip), specific equality constraints are present or not according to the current condition. Since the Coulomb friction condition depends on the normal and tangential contact stresses related to the constraints of the problem, it is devised a conditional dependent constrained minimization problem. An Augmented Lagrangian Method for constrained minimization is employed to solve this problem. This method, when applied to a contact problem, presents Lagrange Multipliers which have the physical meaning of contact forces. This fact allows to check the friction condition at each iteration. These concepts make possible to devise a computational scheme which lead to good numerical results.

A new approach to evaluate imperfection sensitivity in asymmetric bifurcation buckling analysis

A direct procedure for the evaluation of imperfection sensitivity in bifurcation problems is presented. The problems arise in the context of the general theory of elastic stability for discrete structural systems, in which the energy criterion of stability of structures and the total potential energy formulation are employed. In cases of bifurcation buckling the sensitivity of the critical load with respect to an imperfection parameter e is singular at the state given by epsilon =0, so that, a regular perturbation expansion of the solution is not possible. In this work we describe a direct procedure to obtain the relations between the critical loads, the generalized coordinates at the critical state, the eigenvector, and the amplitude of the imperfection, using singular perturbation analysis. The expansions are assumed in terms of arbitrary powers of the imperfection parameter, so that both exponents and coefficients of the expansion are unknown. The solution of the series exponents is obtained by searching the least degenerate solution. The formulation is here applied to asymmetric bifurcations, for which explicit expressions of the coefficients are obtained. The use of the method is illustrated by a simple example, which allows consideration of the main features of the formulation.

Hydrologisen tilanteen ja tuulivoimatuotannon muutosten vaikutus vesiarvoihin pohjoismaisilla sähkön spot-markkinoilla

Vesivoimalla on merkittävä rooli pohjoismaisessa sähköntuotantojärjestelmässä. Spot-markkinoille tarjottavan vesivoiman hinta riippuu vaihtoehtoisen tuotannon hinnasta ja odotetusta vesivoimantuottajien käytettävissä olevasta veteen sitoutuneen potentiaalienergian määrästä. Hydrologisella tilanteella tarkoitetaan tässä tämän potentiaalienergian poikkeamaa normaalitasostaan. Viime vuosina tuulivoimatuotanto pohjoismaisella sähkömarkkina-alueella on kasvanut voimakkaasti, ja on tullut aiheelliseksi tarkastella, millaisia vaikutuksia tällä on vesivoimantuottajien toimintaan. Työssä määritellään vesivoimalle vaihtoehtoisen sähköntuotannon tuotantokustannus, joka pitkällä aikavälillä toimii vertailutasona, jonka perusteella vesivoimantuottajat määrittävät markkinoilla tarjontahinnan tuotannolleen. Tarkastellaan, kuinka hydrologisen tilanteen ja vaihtoehtoisen tuotannon tuotantokustannusten muutokset vaikuttavat vesiarvoon, joka on hinta, jolla hintariippuvaista eli säätyvää vesivoimaa tarjotaan spot-markkinoille. Todetaan, että hydrologisen tilanteen vahvistuminen ja vaihtoehtoisen tuotantokustannuksen aleneminen alentavat vesiarvoja. Todetaan lisäksi, että tuulivoima vaikuttaa sähkön hinnanmuodostukseen markkinoilla samankaltaisesti kuin hintariippumaton vesivoimatuotanto. Esitetään aikasarjamalli vesivoimatuotannon hintariippuvuuden mallintamiseksi. Vertaillaan vesivoimatuottajien toimintaa kahdella vesivoimatuotantoa sisältävällä hinta-alueella, joista toisella tuulivoimatuotanto on kasvanut voimakkaammin kuin toisella. Havaitaan, että molemmilla hinta-alueilla hydrologisen tilanteen vahvistuminen on alentanut ja heikkeneminen nostanut vesiarvoja. Lisäksi havaitaan, että alueella, jonka tuulivoimatuotanto on kasvanut enemmän, vesiarvot ovat laskeneet suhteessa alueen, jolla tuulivoimatuotanto on kasvanut vähemmän, vesiarvoihin. Tuulivoiman voidaan todeta syrjäyttäneen markkinoilta tuotantokustannuksiltaan kalliimpaa tuotantoa.

The vibrational spectra and the normal coordinate analysis of thiocarbonye halides /|nby C. R. Subramaniam. -- 260 St. Catharines, Ont. : [s. n.],

The infrared and the Raman spectra of eSelF has been obtained for the first time and has been analysed to give the in-plane normal vibrational frequencies of the molecule, in the ground state. A normal co-ordinate analysis has been carried out for the molecules CSF2, CSClF and eSel 2 using a Urey-Bradley type of potential function and the elements of the [L] matrix elements, the distribution of the potential energy in Urey-Bradley space, and the displacement vector diagrams for the normal modes of vibration for these molecules, have been obtained. The bond for~e constants obtained through the normal co-ordinate analysis, have given some interesting results. The stretching force constant, Kes ' varies markedly with halogen substitution and the force constants KeF and Keel also vary with substitution.

Assessment of Convective Activity Using Stability Indices as Inferred from Radiosonde and MODIS Data

The combined use of both radiosonde data and three-dimensional satellite derived data over ocean and land is useful for a better understanding of atmospheric thermodynamics. Here, an attempt is made to study the ther-modynamic structure of convective atmosphere during pre-monsoon season over southwest peninsular India utilizing satellite derived data and radiosonde data. The stability indices were computed for the selected stations over southwest peninsular India viz: Thiruvananthapuram and Cochin, using the radiosonde data for five pre- monsoon seasons. The stability indices studied for the region are Showalter Index (SI), K Index (KI), Lifted In-dex (LI), Total Totals Index (TTI), Humidity Index (HI), Deep Convective Index (DCI) and thermodynamic pa-rameters such as Convective Available Potential Energy (CAPE) and Convective Inhibition Energy (CINE). The traditional Showalter Index has been modified to incorporate the thermodynamics over tropical region. MODIS data over South Peninsular India is also used for the study. When there is a convective system over south penin-sular India, the value of LI over the region is less than −4. On the other hand, the region where LI is more than 2 is comparatively stable without any convection. Similarly, when KI values are in the range 35 to 40, there is a possibility for convection. The threshold value for TTI is found to be between 50 and 55. Further, we found that prior to convection, dry bulb temperature at 1000, 850, 700 and 500 hPa is minimum and the dew point tem-perature is a maximum, which leads to increase in relative humidity. The total column water vapor is maximum in the convective region and minimum in the stable region. The threshold values for the different stability indices are found to be agreeing with that reported in literature.

An ab-initio calculation of the Coulomb explosion of N_2 after heavy-ion bombardement

Self-consistent-field calculations for the total potential energy of highly ionized N_2 molecules are presented. We compare these calculations to the experimentally observed energy released in the Coulomb explosion of ionized N_2 molecules created after collision with fast heavy ions. The most important electronic states of the fragment ions are determined.

Total differential scattering cross section of {Ar^+}-Ar at 15 to 400 keV

We report on the measurement of the total differential scattering cross section of {Ar^+}-Ar at laboratory energies between 15 and 400 keV. Using an ab initio relativistic molecular program which calculates the interatomic potential energy curve with high accuracy, we are able to reproduce the detailed structure found in the experiment.

Theoretical evidence for quasi-molecular structure at small internuclear distances in elastic ion-atom scattering

The potential energy curve of the system Ne-Ne is calculated for small internuclear distances from 0.005 to 3.0 au using a newly developed relativistic molecular Dirac-Fock-Slater code. A significant structure in the potential energy curve is found which leads to a nearly complete agreement with experimental differential elastic scattering cross sections. This demonstrates the presence of quasi-molecular effects in elastic ion-atom collisions at keV energies.

A systematic and feasible method for computing nuclear contributions to electrical properties of polyatomic molecules

An analytic method to evaluate nuclear contributions to electrical properties of polyatomic molecules is presented. Such contributions control changes induced by an electric field on equilibrium geometry (nuclear relaxation contribution) and vibrational motion (vibrational contribution) of a molecular system. Expressions to compute the nuclear contributions have been derived from a power series expansion of the potential energy. These contributions to the electrical properties are given in terms of energy derivatives with respect to normal coordinates, electric field intensity or both. Only one calculation of such derivatives at the field-free equilibrium geometry is required. To show the useful efficiency of the analytical evaluation of electrical properties (the so-called AEEP method), results for calculations on water and pyridine at the SCF/TZ2P and the MP2/TZ2P levels of theory are reported. The results obtained are compared with previous theoretical calculations and with experimental values

Variational calculation of static and dynamic vibrational nonlinear optical properties

The vibrational configuration interaction method used to obtain static vibrational (hyper)polarizabilities is extended to dynamic nonlinear optical properties in the infinite optical frequency approximation. Illustrative calculations are carried out on H2 O and N H3. The former molecule is weakly anharmonic while the latter contains a strongly anharmonic umbrella mode. The effect on vibrational (hyper)polarizabilities due to various truncations of the potential energy and property surfaces involved in the calculation are examined

Pressure-jump-induced kinetics reveals a hydration dependent folding/unfolding mechanism of ribonuclease A

Pressure-jump (p-jump)-induced relaxation kinetics was used to explore the energy landscape of protein folding/unfolding of Y115W, a fluorescent variant of ribonuclease A. Pressure-jumps of 40MPa amplitude (5ms dead-time) were conducted both to higher (unfolding) and to lower (folding) pressure, in the range from 100 to 500MPa, between 30 and 50°C. Significant deviations from the expected symmetrical protein relaxation kinetics were observed. Whereas downward p-jumps resulted always in single exponential kinetics, the kinetics induced by upward p-jumps were biphasic in the low pressure range and monophasic at higher pressures. The relative amplitude of the slow phase decreased as a function of both pressure and temperature. At 50°C, only the fast phase remained. These results can be interpreted within the framework of a two-dimensional energy surface containing a pressure- and temperature-dependent barrier between two unfolded states differing in the isomeric state of the Asn-113–Pro-114 bond. Analysis of the activation volume of the fast kinetic phase revealed a temperature-dependent shift of the unfolding transition state to a larger volume. The observed compensation of this effect by glycerol offers an explanation for its protein stabilizing effect

Understanding mixing efficiency in the oceans. Do the nonlinearities of the equation of state matter?

There exist two central measures of turbulent mixing in turbulent stratified fluids that are both caused by molecular diffusion: 1) the dissipation rate D(APE) of available potential energy APE; 2) the turbulent rate of change Wr, turbulent of background gravitational potential energy GPEr. So far, these two quantities have often been regarded as the same energy conversion, namely the irreversible conversion of APE into GPEr, owing to the well known exact equality D(APE)=Wr, turbulent for a Boussinesq fluid with a linear equation of state. Recently, however, Tailleux (2009) pointed out that the above equality no longer holds for a thermally-stratified compressible, with the ratio ξ=Wr, turbulent/D(APE) being generally lower than unity and sometimes even negative for water or seawater, and argued that D(APE) and Wr, turbulent actually represent two distinct types of energy conversion, respectively the dissipation of APE into one particular subcomponent of internal energy called the "dead" internal energy IE0, and the conversion between GPEr and a different subcomponent of internal energy called "exergy" IEexergy. In this paper, the behaviour of the ratio ξ is examined for different stratifications having all the same buoyancy frequency N vertical profile, but different vertical profiles of the parameter Υ=α P/(ρCp), where α is the thermal expansion coefficient, P the hydrostatic pressure, ρ the density, and Cp the specific heat capacity at constant pressure, the equation of state being that for seawater for different particular constant values of salinity. It is found that ξ and Wr, turbulent depend critically on the sign and magnitude of dΥ/dz, in contrast with D(APE), which appears largely unaffected by the latter. These results have important consequences for how the mixing efficiency should be defined and measured in practice, which are discussed.

Mesoscale simulations of organized convection: importance of convective equilibrium

The validity of convective parametrization breaks down at the resolution of mesoscale models, and the success of parametrized versus explicit treatments of convection is likely to depend on the large-scale environment. In this paper we examine the hypothesis that a key feature determining the sensitivity to the environment is whether the forcing of convection is sufficiently homogeneous and slowly varying that the convection can be considered to be in equilibrium. Two case studies of mesoscale convective systems over the UK, one where equilibrium conditions are expected and one where equilibrium is unlikely, are simulated using a mesoscale forecasting model. The time evolution of area-average convective available potential energy and the time evolution and magnitude of the timescale of convective adjustment are consistent with the hypothesis of equilibrium for case 1 and non-equilibrium for case 2. For each case, three experiments are performed with different partitionings between parametrized and explicit convection: fully parametrized convection, fully explicit convection and a simulation with significant amounts of both. In the equilibrium case, bulk properties of the convection such as area-integrated rain rates are insensitive to the treatment of convection. However, the detailed structure of the precipitation field changes; the simulation with parametrized convection behaves well and produces a smooth field that follows the forcing region, and the simulation with explicit convection has a small number of localized intense regions of precipitation that track with the mid-levelflow. For the non-equilibrium case, bulk properties of the convection such as area-integrated rain rates are sensitive to the treatment of convection. The simulation with explicit convection behaves similarly to the equilibrium case with a few localized precipitation regions. In contrast, the cumulus parametrization fails dramatically and develops intense propagating bows of precipitation that were not observed. The simulations with both parametrized and explicit convection follow the pattern seen in the other experiments, with a transition over the duration of the run from parametrized to explicit precipitation. The impact of convection on the large-scaleflow, as measured by upper-level wind and potential-vorticity perturbations, is very sensitive to the partitioning of convection for both cases. © Royal Meteorological Society, 2006. Contributions by P. A. Clark and M. E. B. Gray are Crown Copyright.

Normal coordinates in the methyl halides

The potential‐energy functions found by Chang for the methyl halides have been put into valence‐type form and revised to eliminate inconsistencies and to accord with the true (nontetrahedral) geometry and the normal frequencies (corrected for Fermi resonance and anharmonicity). The resulting valence‐type force constants and normal coordinates are given for light (CH3) and heavy (CD3) chlorides, bromides, and iodides.