Explicit/implicit partitioning and a new explicit form of the generalized alpha method

Most finite element packages use the Newmark algorithm for time integration of structural dynamics. Various algorithms have been proposed to better optimize the high frequency dissipation of this algorithm. Hulbert and Chung proposed both implicit and explicit forms of the generalized alpha method. The algorithms optimize high frequency dissipation effectively, and despite recent work on algorithms that possess momentum conserving/energy dissipative properties in a non-linear context, the generalized alpha method remains an efficient way to solve many problems, especially with adaptive timestep control. However, the implicit and explicit algorithms use incompatible parameter sets and cannot be used together in a spatial partition, whereas this can be done for the Newmark algorithm, as Hughes and Liu demonstrated, and for the HHT-alpha algorithm developed from it. The present paper shows that the explicit generalized alpha method can be rewritten so that it becomes compatible with the implicit form. All four algorithmic parameters can be matched between the explicit and implicit forms. An element interface between implicit and explicit partitions can then be used, analogous to that devised by Hughes and Liu to extend the Newmark method. The stability of the explicit/implicit algorithm is examined in a linear context and found to exceed that of the explicit partition. The element partition is significantly less dissipative of intermediate frequencies than one using the HHT-alpha method. The explicit algorithm can also be rewritten so that the discrete equation of motion evaluates forces from displacements and velocities found at the predicted mid-point of a cycle. Copyright (C) 2003 John Wiley Sons, Ltd.

Adaptations of strangler figs to life in the rainforest canopy

Figs are rainforest keystone species. Non-strangler figs establish on the forest floor; strangler figs establish epiphytically, followed by a dramatic transition from epiphyte to free-standing tree that kills its hosts. Free-standing figs display vigorous growth and resource demand suggesting that epiphytic strangler figs require special adaptations to deal with resource limitations imposed by the epiphytic environment. We studied epiphytic and free-standing strangler figs, and non-strangler figs in tropical rainforest and in cultivation, as well as strangler figs in controlled conditions. We investigated whether the transition from epiphyte to free-standing tree is characterised by morphological and physiological plasticity. Epiphyte substrate had higher levels of plant-available ammonium and phosphate, and similar levels of nitrate compared with rainforest soil, suggesting that N and P are initially not limiting resources. A relationship was found between taxonomic groups and plant N physiology; strangler figs, all members of subgenus Urostigma, had mostly low foliar nitrate assimilation rates whereas non-strangler figs, in subgenera Pharmacocycea, Sycidium, Sycomorus or Synoecia, had moderate to high rates. Nitrate is an energetically expensive N source, and low nitrate use may be an adaptation of strangler figs for conserving energy during epiphytic growth. Interestingly, significant amounts of nitrate were stored in fleshy taproot tubers of epiphytic stranglers. Supporting the concept of plasticity, leaves of epiphytic Ficus benjamina L. had lower N and C content per unit leaf area, lower stomatal density and 80% greater specific leaf area than leaves of conspecific free-standing trees. Similarly, glasshouse-grown stranglers strongly increased biomass allocation to roots under water limitation. Epiphytic and free-standing F. benjamina had similar average foliar delta C-13, but epiphytes had more extreme values; this indicates that both groups of plants use the C-3 pathway of CO2 fixation but that water availability is highly variable for epiphytes. We hypothesise that epiphytic figs use fleshy stem tubers to avoid water stress, and that nitrate acts as an osmotic compound in tubers. We conclude that strangler figs are a unique experimental system for studying the transition from rainforest epiphyte to tree, and the genetic and environmental triggers involved.

Kinetic energy conserving integrators for Gaussian thermostatted SLLOD

A new integration scheme is developed for nonequilibrium molecular dynamics simulations where the temperature is constrained by a Gaussian thermostat. The utility of the scheme is demonstrated by its application to the SLLOD algorithm which is the standard nonequilibrium molecular dynamics algorithm for studying shear flow. Unlike conventional integrators, the new integrators are constructed using operator-splitting techniques to ensure stability and that little or no drift in the kinetic energy occurs. Moreover, they require minimum computer memory and are straightforward to program. Numerical experiments show that the efficiency and stability of the new integrators compare favorably with conventional integrators such as the Runge-Kutta and Gear predictor-corrector methods. (C) 1999 American Institute of Physics. [S0021-9606(99)50125-6].

Correcting physical activity energy expenditure for body size in young children

Being able to compare the energy cost of physical activity across and between populations is important. However, energy expenditure is related to body size, so it is necessary to appropriately adjust for differences in body size when comparisons are made. This study examined the relationship between the daily energy cost of activity and body weight in 47 children aged 6-10 years. Log-log regression showed weight(1.0) to be an inappropriate adjustment for activity energy expenditure in children, with a more valid adjustment being weight(0.3). Clearly, both weight dependent and non-weight dependent activities are part of everyday living in children. This balance influences how energy expenditure is correctly adjusted for body size. Investigators interpreting data of energy expenditure in children from children of different body sizes need to take this into consideration.

Hydraulic Design of Stepped Spillways and Downstream Energy Dissipators for Embankment Dams

In recent years, the design flows of many dams were re-evaluated, often resulting in discharges larger than the original design. In many cases, the occurrence of the revised flows could result in dam overtopping because of insufficient storage and spillway capacity. An experimental study was conducted herein to gain a better understanding of the flow properties in stepped chutes with slopes typical of embankment dams. The work was based upon a Froude similitude in large-size experimental facilities. A total of 10 configurations were tested including smooth steps, steps equipped with devices to enhance energy dissipation and rough steps. The present results yield a new design procedure. The design method includes some key issues not foreseen in prior studies : e.g., gradually varied flow, type of flow regime, flow resistance. It is believed that the outcomes are valid for a wide range of chute geometry and flow conditions typical of embankment chutes.

Least energy solutions of a critical Neumann problem with a weight

We investigate the effect of the coefficient of the critical nonlinearity for the Neumann problem on the existence of least energy solutions. As a by-product we establish a Sobolev inequality with interior norm.

Algebraic Bethe ansatz method for the exact calculation of energy spectra and form factors: applications to models of Bose-Einstein condensates and metallic nanograins

In this review we demonstrate how the algebraic Bethe ansatz is used for the calculation of the-energy spectra and form factors (operator matrix elements in the basis of Hamiltonian eigenstates) in exactly solvable quantum systems. As examples we apply the theory to several models of current interest in the study of Bose-Einstein condensates, which have been successfully created using ultracold dilute atomic gases. The first model we introduce describes Josephson tunnelling between two coupled Bose-Einstein condensates. It can be used not only for the study of tunnelling between condensates of atomic gases, but for solid state Josephson junctions and coupled Cooper pair boxes. The theory is also applicable to models of atomic-molecular Bose-Einstein condensates, with two examples given and analysed. Additionally, these same two models are relevant to studies in quantum optics; Finally, we discuss the model of Bardeen, Cooper and Schrieffer in this framework, which is appropriate for systems of ultracold fermionic atomic gases, as well as being applicable for the description of superconducting correlations in metallic grains with nanoscale dimensions.; In applying all the above models to. physical situations, the need for an exact analysis of small-scale systems is established due to large quantum fluctuations which render mean-field approaches inaccurate.

Comparative public support for conserving reptile species is high: Australian evidence and its implications

This paper investigates factors influencing the public’s support for conservation of tropical reptile species in a focal group drawing on Australian data and an experiment involving a sample of the Australian public. The influences of the likeability of the species, their degree of endangerment, ethical considerations as well as knowledge are examined and found to be important. Likeability is found to be much less important than the existing literature suggests. This is highlighted by comparing the likeability of the focal group of reptiles with that for a group of birds and a group of mammals with differences in willingness to pay for their conservation.

Net uptake of dissolved free amino acids by the giant clam, Tridacna maxima: alternative sources of energy and nitrogen?

The role of dissolved free amino acids (DFAA) in nitrogen and energy budgets was investigated for the giant clam, Tridacna maxima, growing under field conditions at One Tree Island, at the southern end of the Great Barrier Reef, Australia. Giant clams (121.5-143.7 mm in shell length) took up neutral, acidic and basic amino acids. The rates of net uptake of DFAA did not differ between light and dark, nor for clams growing under normal or slightly enriched ammonium concentrations. Calculations based on the net uptake concentrations typical of the maximum concentrations of DFAA found in coral reef waters (similar to 0.1 mu M)revealed that DFAA could only contribute 0.1% and 1% of the energy and nitrogen demands of giant clams, respectively. These results suggest that DFAA does not supply significant amounts of energy or nitrogen for giant clams or their symbionts.

Assessment of limb muscle and adipose tissue by dual-energy X-ray absorptiometry using magnetic resonance imaging for comparison

OBJECTIVE: To use magnetic resonance imaging (MRI) to validate estimates of muscle and adipose tissue (AT) in lower limb sections obtained by dual-energy X-ray absorptiometry (DXA) modelling. DESIGN: MRI measurements were used as reference for validating limb muscle and AT estimates obtained by DXA models that assume fat-free soft tissue (FFST) comprised mainly muscle: model A accounted for bone hydration only; model B also applied constants for FFST in bone and skin and fat in muscle and AT; model C was as model B but allowing for variable fat in muscle and AT. SUBJECTS: Healthy men (n = 8) and women (n = 8), ages 41 - 62 y; mean (s.d.) body mass indices (BMIs) of 28.6 (5.4) kg/m(2) and 25.1 (5.4) kg/m2, respectively. MEASUREMENTS: MRI scans of the legs and whole body DXA scans were analysed for muscle and AT content of thigh (20 cm) and lower leg (10 cm) sections; 24 h creatinine excretion was measured. RESULTS: Model A overestimated thigh muscle volume (MRI mean, 2.3 l) substantially (bias 0.36 l), whereas model B underestimated it by only 2% (bias 0.045 l). Lower leg muscle (MRI mean, 0.6 l) was better predicted using model A (bias 0.04 l, 7% overestimate) than model B (bias 0.1 l, 17% underestimate). The 95% limits of agreement were high for these models (thigh,+/- 20%; lower leg,+/- 47%). Model C predictions were more discrepant than those of model B. There was generally less agreement between MRI and all DXA models for AT. Measurement variability was generally less for DXA measurements of FFST (coefficient of variation 0.7 - 1.8%) and fat (0.8 - 3.3%) than model B estimates of muscle (0.5-2.6%) and AT (3.3 - 6.8%), respectively. Despite strong relationships between them, muscle mass was overestimated by creatinine excretion with highly variable predictability. CONCLUSION: This study has shown the value of DXA models for assessment of muscle and AT in leg sections, but suggests the need to re-evaluate some of the assumptions upon which they are based.

Determination of activation energy distributions for chemisorption of oxygen on carbon: an improved approach

The present paper proposes an approach to obtaining the activation energy distribution for chemisorption of oxygen onto carbon surfaces, while simultaneously allowing for the activation energy dependence of the pre-exponential factor of the rate constant. Prior studies in this area have considered this factor to be uniform, thereby biasing estimated distributions. The results show that the derived activation energy distribution is not sensitive to the chemisorption mechanism because of the step function like property of the coverage. The activation energy distribution is essentially uniform for some carbons, and has two or possibly more discrete stages, suggestive of at least two types of sites, each with its own uniform distribution. The pre-exponential factors of the reactions are determined directly from the experimental data, and are found not to be constant as assumed in earlier work, but correlated with the activation energy. The latter results empirically follow an exponential function, supporting some earlier statistical and experimental work. The activation energy distribution obtained in the present paper permits improved correlation of chemisorption data in comparison to earlier studies. (C) 2000 Elsevier Science Ltd. All rights reserved.

Accelerating seismic energy release and evolution of event time and size statistics: Results from two heterogeneous cellular automaton models

The evolution of event time and size statistics in two heterogeneous cellular automaton models of earthquake behavior are studied and compared to the evolution of these quantities during observed periods of accelerating seismic energy release Drier to large earthquakes. The two automata have different nearest neighbor laws, one of which produces self-organized critical (SOC) behavior (PSD model) and the other which produces quasi-periodic large events (crack model). In the PSD model periods of accelerating energy release before large events are rare. In the crack model, many large events are preceded by periods of accelerating energy release. When compared to randomized event catalogs, accelerating energy release before large events occurs more often than random in the crack model but less often than random in the PSD model; it is easier to tell the crack and PSD model results apart from each other than to tell either model apart from a random catalog. The evolution of event sizes during the accelerating energy release sequences in all models is compared to that of observed sequences. The accelerating energy release sequences in the crack model consist of an increase in the rate of events of all sizes, consistent with observations from a small number of natural cases, however inconsistent with a larger number of cases in which there is an increase in the rate of only moderate-sized events. On average, no increase in the rate of events of any size is seen before large events in the PSD model.

Free energy approximations in simple lattice proteins

This work addresses the question of whether it is possible to define simple pairwise interaction terms to approximate free energies of proteins or polymers. Rather than ask how reliable a potential of mean force is, one can ask how reliable it could possibly be. In a two-dimensional, infinite lattice model system one can calculate exact free energies by exhaustive enumeration. A series of approximations were fitted to exact results to assess the feasibility and utility of pairwise free energy terms. Approximating the true free energy with pairwise interactions gives a poor fit with little transferability between systems of different size. Adding extra artificial terms to the approximation yields better fits, but does not improve the ability to generalize from one system size to another. Furthermore, one cannot distinguish folding from nonfolding sequences via the approximated free energies. Most usefully, the methodology shows how one can assess the utility of various terms in lattice protein/polymer models. (C) 2001 American Institute of Physics.