Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

Coupling of a Jahn-Teller pseudorotation with a hindered internal rotation in an isolated molecule: 9-hydroxytriptycene

The irregular vibronic structure in the S-1<--S-0 resonant two-photon ionization (R2PI) spectrum of supersonically cooled triptycene is a result of a classic Exe Jahn-Teller effect [A. Furlan et al., J. Chem. Phys. 96, 7306 (1992)]. This is well characterized and can be used as an effective probe of intramolecular perturbations. Here we examine the S-1<--S-0 R2PI spectrum of 9-hydroxytriptycene and the fluorescence from various excited state vibronic levels. In this system the pseudorotation of the Jahn-Teller vibration is strongly coupled to the torsional motion of the bridgehead hydroxy group. This torsional motion results in a tunneling splitting in both the ground and excited states. The population of the upper level in the ground electronic state results in additional vibronic transitions becoming symmetry allowed in the R2PI spectrum that are forbidden in the bare triptycene molecule. The assignment of the R2PI and fluorescence spectra allows the potential energy surfaces of these vibrational modes to be accurately quantified. The full C-3v vibronic point group must be used to interpret the spectra. The time scale of the internal rotation of the-OH group and the butterfly flapping of the Jahn-Teller pseudorotation are of similar magnitude. The tunneling between the nine minima on the three dimensional potential energy surface is such that the Jahn-Teller pseudorotation occurs in concert with the-OH internal rotation. The Berry phase that is acquired during this motion is discussed. The simple physical picture emerges of the angle between two of the three benzene moieties opening in three equivalent ways in the S-1 electronic state. This geometry follows the position of the hydroxy group, which preferentially orients itself to point between these two rings. (C) 1998 American Institute of Physics. [S0021-9606(98)02348-4].

Evolution of three-dimensional folds for a non-Newtonian plate in a viscous medium

We derive analytical solutions for the three-dimensional time-dependent buckling of a non-Newtonian viscous plate in a less viscous medium. For the plate we assume a power-law rheology. The principal, axes of the stretching D-ij in the homogeneously deformed ground state are parallel and orthogonal to the bounding surfaces of the plate in the flat state. In the model formulation the action of the less viscous medium is replaced by equivalent reaction forces. The reaction forces are assumed to be parallel to the normal vector of the deformed plate surfaces. As a consequence, the buckling process is driven by the differences between the in-plane stresses and out of plane stress, and not by the in-plane stresses alone as assumed in previous models. The governing differential equation is essentially an orthotropic plate equation for rate dependent material, under biaxial pre-stress, supported by a viscous medium. The differential problem is solved by means of Fourier transformation and largest growth coefficients and corresponding wavenumbers are evaluated. We discuss in detail fold evolutions for isotropic in-plane stretching (D-11 = D-22), uniaxial plane straining (D-22 = 0) and in-plane flattening (D-11 = -2D(22)). Three-dimensional plots illustrate the stages of fold evolution for random initial perturbations or initial embryonic folds with axes non-parallel to the maximum compression axis. For all situations, one dominant set of folds develops normal to D-11, although the dominant wavelength differs from the Biot dominant wavelength except when the plate has a purely Newtonian viscosity. However, in the direction parallel to D-22, there exist infinitely many modes in the vicinity of the dominant wavelength which grow only marginally slower than the one corresponding to the dominant wavelength. This means that, except for very special initial conditions, the appearance of a three-dimensional fold will always be governed by at least two wavelengths. The wavelength in the direction parallel to D-11 is the dominant wavelength, and the wavelength(s) in the direction parallel to D-22 is determined essentially by the statistics of the initial state. A comparable sensitivity to the initial geometry does not exist in the classic two-dimensional folding models. In conformity with tradition we have applied Kirchhoff's hypothesis to constrain the cross-sectional rotations of the plate. We investigate the validity of this hypothesis within the framework of Reissner's plate theory. We also include a discussion of the effects of adding elasticity into the constitutive relations and show that there exist critical ratios of the relaxation times of the plate and the embedding medium for which two dominant wavelengths develop, one at ca. 2.5 of the classical Biot dominant wavelength and the other at ca. 0.45 of this wavelength. We propose that herein lies the origin of parasitic folds well known in natural examples.

Simulation and optimization of heated, inviscid flows in scramjet ducts

A space-marching code for the simulation and optimization of inviscid supersonic flow in three dimensions is described. The now in a scramjet module with a relatively complex three-dimensional geometry is examined and wall-pressure estimates are compared with experimental data. Given that viscous effects are not presently included, the comparison is reasonable. The thermodynamic compromise of adding heat in a diverging combustor is also examined. The code is then used to optimize the shape of a thrust surface for a simpler (box-section) scramjet module in the presence of uniform and nonuniform heat distributions. The optimum two-dimensional profiles for the thrust surface are obtained via a perturbation procedure that requires about 30-50 now solutions. It is found that the final shapes are fairly insensitive to the details of the heat distribution.

Selectively excited luminescence and magnetic circular dichroism of Cr4+-doped YAG and YGG

Site selective luminescence and magnetic circular dichroism experiments on Cr4+-doped yttrium aluminum garnet and yttrium gallium garnet have been made at low temperature. The spectral assignments for these near-IR lasing materials have been made using experimental data and ligand field calculations guided by the known geometry of the lattices. [S0163-1829(99)07003-4].

Solution structure of a defensin-like peptide from platypus venom

Three defensin-like peptides (DLPs) were isolated from platypus venom and sequenced. One of these peptides, DLP-1, was synthesized chemically and its three-dimensional structure was determined using NMR spectroscopy. The main structural elements of this 42-residue peptide were an anti-parallel beta-sheet comprising residues 15-18 and 37-40 and a small 3(10) helix spanning residues 10-12. The overall three-dimensional fold is similar to that of beta-defensin-12, and similar to the sodium-channel neurotoxin ShI (Stichodactyla helianthus neurotoxin I). However, the side chains known to be functionally important in beta-defensin-12 and ShI are not conserved in DLP-1, suggesting that it has a different biological function. Consistent with this contention, we showed that DLP-1 possesses no anti-microbial properties and has no observable activity on rat dorsal-root-ganglion sodium-channel currents.

Using Lindenmayer systems to model morphogenesis in a tropical pasture legume Stylosanthes scabra

This paper summarizes the processes involved in designing a mathematical model of a growing pasture plant, Stylosanthes scabra Vog. cv. Fitzroy. The model is based on the mathematical formalism of Lindenmayer systems and yields realistic computer-generated images of progressive plant geometry through time. The processes involved in attaining growth data, retrieving useful growth rules, and constructing a virtual plant model are outlined. Progressive output morphological data proved useful for predicting total leaf area and allowed for easier quantification of plant canopy size in terms of biomass and total leaf area.

Determining forest structural attributes using an inverted geometric-optical model in mixed eucalypt forests, Southeast Queensland, Australia

The Montreal Process indicators are intended to provide a common framework for assessing and reviewing progress toward sustainable forest management. The potential of a combined geometrical-optical/spectral mixture analysis model was assessed for mapping the Montreal Process age class and successional age indicators at a regional scale using Landsat Thematic data. The project location is an area of eucalyptus forest in Emu Creek State Forest, Southeast Queensland, Australia. A quantitative model relating the spectral reflectance of a forest to the illumination geometry, slope, and aspect of the terrain surface and the size, shape, and density, and canopy size. Inversion of this model necessitated the use of spectral mixture analysis to recover subpixel information on the fractional extent of ground scene elements (such as sunlit canopy, shaded canopy, sunlit background, and shaded background). Results obtained fron a sensitivity analysis allowed improved allocation of resources to maximize the predictive accuracy of the model. It was found that modeled estimates of crown cover projection, canopy size, and tree densities had significant agreement with field and air photo-interpreted estimates. However, the accuracy of the successional stage classification was limited. The results obtained highlight the potential for future integration of high and moderate spatial resolution-imaging sensors for monitoring forest structure and condition. (C) Elsevier Science Inc., 2000.

Discovery and characterization of a family of insecticidal neurotoxins with a rare vicinal disulfide bridge

We have isolated a family of insect-selective neurotoxins from the venom of the Australian funnel-web spider that appear to be good candidates for biopesticide engineering. These peptides, which we have named the Janus-faced atracotoxins (J-ACTXs), each contain 36 or 37 residues, with four disulfide bridges, and they show no homology to any sequences in the protein/DNA databases. The three-dimensional structure of one of these toxins reveals an extremely rare vicinal disulfide bridge that we demonstrate to be critical for insecticidal activity. We propose that J-ACTX comprises an ancestral protein fold that we refer to as the disulfide-directed beta-hairpin.

Numerical modelling of single-layer folding: clarification of an issue regarding the possible effect of computer codes and the influence of initial irregularities

The influence of initial perturbation geometry and material propel-ties on final fold geometry has been investigated using finite-difference (FLAC) and finite-element (MARC) numerical models. Previous studies using these two different codes reported very different folding behaviour although the material properties, boundary conditions and initial perturbation geometries were similar. The current results establish that the discrepancy was not due to the different computer codes but due to the different strain rates employed in the two previous studies (i.e. 10(-6) s(-1) in the FLAC models and 10(-14) s(-1) in the MARC models). As a result, different parts of the elasto-viscous rheological field were bring investigated. For the same material properties, strain rate and boundary conditions, the present results using the two different codes are consistent. A transition in Folding behaviour, from a situation where the geometry of initial perturbation determines final fold shape to a situation where material properties control the final geometry, is produced using both models. This transition takes place with increasing strain rate, decreasing elastic moduli or increasing viscosity (reflecting in each case the increasing influence of the elastic component in the Maxwell elastoviscous rheology). The transition described here is mechanically feasible but is associated with very high stresses in the competent layer (on the order of GPa), which is improbable under natural conditions. (C) 2000 Elsevier Science Ltd. All rights reserved.

Heat flow for Yang-Mills-Higgs fields, part I

The Yang-Mills-Higgs field generalizes the Yang-Mills field. The authors establish the local existence and uniqueness of the weak solution to the heat flow for the Yang-Mills-Higgs field in a vector bundle over a compact Riemannian 4-manifold, and show that the weak solution is gauge-equivalent to a smooth solution and there are at most finite singularities at the maximum existing time.

Circular proteins in plants - Solution structure of a novel macrocyclic trypsin inhibitor from Momordica cochinchinensis

Much interest has been generated by recent reports on the discovery of circular (i.e. head-to-tail cyclized) proteins in plants. Here we report the three-dimensional structure of one of the newest such circular proteins, MCoTI-II, a novel trypsin inhibitor from Momordica cochinchinensis, a member of the Cucurbitaceae plant family. The structure consists of a small beta -sheet, several turns, and a cystine knot arrangement of the three disulfide bonds. Interestingly, the molecular topology is similar to that of the plant cyclotides (Craik, D. J., Daly, N. L., Bond, T., and Waine, C. (1999) J. Mol. Biol, 294, 1327-1336), which derive from the Rubiaceae and Violaceae plant families, have antimicrobial activities, and exemplify the cyclic cystine knot structural motif as part of their circular backbone. The sequence, biological activity, and plant family of MCoTI-II are all different from known cyclotides. However, given the structural similarity, cyclic backbone, and plant origin of MCoTI-II, we propose that MCoTI-II can be classified as a new member of the cyclotide class of proteins. The expansion of the cyclotides to include trypsin inhibitory activity and a new plant family highlights the importance and functional variability of circular proteins and the fact that they are more common than has previously been believed, Insights into the possible roles of backbone cyclization have been gained by a comparison of the structure of MCoTI-II with the homologous acyclic trypsin inhibitors CMTI-I and EETI-II from the Cucurbitaceae plant family.

Clarification of the Crystal Structure of [(HC(Ph2PO)3)2Cu](ClO4)2·2H2O

EXAFS spectra of [(HC(Ph2PO)(3))(2)Cu](ClO4)(2). 2H(2)O have been measured at room temperature. These show that the CuO6 unit is tetragonally elongated, rather than having the compressed tetragonal geometry previously inferred from the X-ray crystal structure determination. [GRAPHICS]

Model for gas-liquid flow through wet porous medium

A method involving bubbling of air through a fibrous filter immersed in water has recently been investigated (Agranovski et al. [1]). Experimental results showed that the removal efficiency for ultra-fine aerosols by such filters was greatly increased compared to dry filters. Nuclear Magnetic Resonance (NMR) imaging was used to examine the wet filter and to determine the nature of the gas flow inside the filter (Agranovski et al. [2]). It was found that tortuous preferential pathways (or flow tubes) develop within the filter through which the air flows and the distribution of air and water inside the porous medium has been investigated. The aim of this paper is to investigate the geometry of the pathways and to make estimates of the flow velocities and particle removal efficiency in such pathways. A mathematical model of the flow of air along the preferred pathways has been developed and verified experimentally. Even for the highest realistic gas velocity the flow field was essentially laminar (Re approximate to 250). We solved Laplace's equation for stream function to map trajectories of particles and gas molecules to investigate the possibility of their removal from the carrier.