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.

Mesoionic pyridopyrimidinylium and pyridooxazinylium olates and non-mesoionic pyridopyrimidinones. Structures in the solid state, solution, and matrices

X-Ray crystal structures, C-13 NMR spectra and theoretical calculations (B3LYP/6-31G*) are reported for the mesoionic (zwitterionic) pyridopyrimidinylium- and pyridooxazinyliumolates 2a, 3a and 5a,b as well as the enol ether 11b and the enamine 11c. The 1-NH compounds like 1a, 2a and 3a exist in the mesoionic form in the crystal and in solution, but the OH tautomers such as 1b and 2b dominate in the gas phase as revealed by the Ar matrix IR spectra in conjunction with DFT calculations. All data indicate that the mesoionic compounds can be regarded as intramolecular pyridine-ketene zwitterions (cf. 16 --> 17) with a high degree of positive charge on the pyridinium nitrogen, a long pyridinium N-CO bond (ca. 1.44-1.49 Angstrom), and normal C=O double bonds (ca. 1.22 Angstrom). All mesoionic compounds exhibit a pronounced tilting of the olate C=O groups (the C=O groups formally derived from a ketene) towards the pyridinium nitrogen, giving NCO angles of 110-118 degrees. Calculations reveal a hydrogen bond with 6-CH, analogous to what is found in ketene-pyridine zwitterions and the C3O2-pyridine complex. The 2-OH tautomers of type 1b, 2b, and 11 also show a high degree of zwitterionic character as indicated by the canonical structures 11 12.

Commentary: Using the convection-dispersion model and transit time density functions in the analysis of organ distribution kinetics

The convection-dispersion model and its extended form have been used to describe solute disposition in organs and to predict hepatic availabilities. A range of empirical transit-time density functions has also been used for a similar purpose. The use of the dispersion model with mixed boundary conditions and transit-time density functions has been queried recently by Hisaka and Sugiyanaa in this journal. We suggest that, consistent with soil science and chemical engineering literature, the mixed boundary conditions are appropriate providing concentrations are defined in terms of flux to ensure continuity at the boundaries and mass balance. It is suggested that the use of the inverse Gaussian or other functions as empirical transit-time densities is independent of any boundary condition consideration. The mixed boundary condition solutions of the convection-dispersion model are the easiest to use when linear kinetics applies. In contrast, the closed conditions are easier to apply in a numerical analysis of nonlinear disposition of solutes in organs. We therefore argue that the use of hepatic elimination models should be based on pragmatic considerations, giving emphasis to using the simplest or easiest solution that will give a sufficiently accurate prediction of hepatic pharmacokinetics for a particular application. (C) 2000 Wiley-Liss Inc. and the American Pharmaceutical Association J Pharm Sci 89:1579-1586, 2000.

Tidal fluctuations in a leaky confined aquifer: Dynamic effects of an overlying phreatic aquifer

Tidal fluctuations in a leaky confined coastal aquifer are damped significantly due to leakage into an overlying phreatic aquifer. Jiao and Tang [1999] presented an analytical solution to a simple model describing this phenomenon. Their solution assumes that the tidal fluctuations in the overlying phreatic aquifer are negligible (i.e,, a static phreatic aquifer), Here we examine dynamic effects of the overlying aquifer based on a new approximate analytical solution. The numerical results indicate that the dynamic effects can be significant for a relatively large leakage and a high transmissivity of the phreatic aquifer.

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.

Comments on "Exact solution for horizontal redistribution by general similarity"

No abstract

Sudden drawdown and drainage of a horizontal aquifer

Drainage of a saturated horizontal aquifer following a sudden drawdown is reanalyzed using the Boussinesq equation. The effect of the finite length of the aquifer is considered in detail. An analytical approximation based on a superposition principle yields a very good estimate of the outflow when compared to accurate numerical solutions. An illustration of the new analytical approach to analyze basin-scale field data is used to demonstrate possible field applications of the new solution.

Conformations of platypus venom C-type natriuretic peptide in aqueous solution and sodium dodecyl sulfate micelles

Nuclear magnetic resonance spectroscopy was used to investigate the conformations of the platypus venom C-type natriuretic peptide A (OvCNPa) in aqueous solutions and in solutions containing sodium dodecyl sulfate (SDS) micelles. The chemically synthesized OvCNPa showed a substantial decrease in flexibility in aqueous solution at 10 degreesC, allowing the observation of medium- and long-range nuclear Overhauser enhancement (NOE) connectivities. Three-dimensional structures calculated using these data showed flexible and reasonably well-defined regions, the locations of which were similar in the two solvents. In aqueous solution, the linear part that spans residues 3-14 was basically an extended conformation while the cyclic portion, defined by residues 23-39, contained a series of beta-turns. The overall shape of the cyclic portion was similar to that observed for an atrial natriuretic peptide (ANP) variant in aqueous solution. OvCNPa adopted a different conformation in SDS micelles wherein the N-terminal region, defined by residues 2-10, was more compact, characterised by turns and a helix, while the cyclic region had turns and an overall shape that was fundamentally different from those structures observed in aqueous solution. The hydrophobic cluster, situated at the centre of the ring of the structure in aqueous solution, was absent in the structure in the presence of SDS micelles. Thus, OvCNPa interacts with SDS micelles and can possibly form ion-channels in cell membranes. (C) 2002 Elsevier Science Ltd. All rights reserved.

Analytical solution for tidal propagation in a coupled semi-confined/phreatic coastal aquifer

An analytical solution is derived for tidal fluctuations in a coupled coastal aquifer system consisting of a semi-confined aquifer, a thin semi-permeable layer and a phreatic aquifer. Based on the solution, we study the interactions (via leakage) between the confined and unconfined aquifers in response to tides. The results show that, under certain conditions, leakage from the confined aquifer can affect considerably the tidal water table fluctuation in the phreatic aquifer and vice versa. Ignoring these effects could lead to errors in estimating aquifer properties based on tidal signals. (C) 2002 Elsevier Science Ltd. All rights reserved.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Mantle convection modeling with viscoelastic/brittle lithosphere: Numerical methodology and plate tectonic modeling

The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.