Computer simulations of coupled problems in geological and geochemical systems

The finite element method is used to simulate coupled problems, which describe the related physical and chemical processes of ore body formation and mineralization, in geological and geochemical systems. The main purpose of this paper is to illustrate some simulation results for different types of modelling problems in pore-fluid saturated rock masses. The aims of the simulation results presented in this paper are: (1) getting a better understanding of the processes and mechanisms of ore body formation and mineralization in the upper crust of the Earth; (2) demonstrating the usefulness and applicability of the finite element method in dealing with a wide range of coupled problems in geological and geochemical systems; (3) qualitatively establishing a set of showcase problems, against which any numerical method and computer package can be reasonably validated. (C) 2002 Published by Elsevier Science B.V.

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.

Study of the antimicrobial peptide indolicidin and mutants in eukaryotic modelled membrane by molecular dynamics simulations

In this work the interaction of the antimicrobial peptide indolicidin (IND) and its mutants CP10A and CP11 with a eukaryotic membrane model was examined by molecular dynamics simulations. The aim was to analyse the behaviour of these antimicrobial peptides when they interact with a eukaryotic modelled membrane, thereby obtaining atomic detailed observations that are not experimentally available. In the simulations, the widely studied dipalmitoylphosphatidylcholine hydrated bilayer was used as a eukaryotic membrane model. In agreement with experimental observations, the peptides IND, CP10A, and CP11 insert into the bilayer differently; the peptides that insert more deeply present the major hemolytic activities. The hydrophobic residues are responsible for the insertion, but some Trp residues of the peptides remain at the bilayer/water interface because they interact with the bilayer choline groups by cation-pi interactions that should be important for recognition of eukaryotic membrane by the three studied peptides.

Stochastic gauges in quantum dynamics for many-body simulations

Quantum dynamics simulations can be improved using novel quasiprobability distributions based on non-orthogonal Hermitian kernel operators. This introduces arbitrary functions (gauges) into the stochastic equations. which can be used to tailor them for improved calculations. A possible application to full quantum dynamic simulations of BEC's is presented. (C) 2001 Elsevier Science B.V. All rights reserved.

Simulations of thermal Bose fields in the classical limit

We demonstrate that the time-dependent projected Gross-Pitaevskii equation (GPE) derived earlier [M. J. Davis, R. J. Ballagh, and K. Burnett, J. Phys. B 34, 4487 (2001)] can represent the highly occupied modes of a homogeneous, partially-condensed Bose gas. Contrary to the often held belief that the GPE is valid only at zero temperature, we find that this equation will evolve randomized initial wave functions to a state describing thermal equilibrium. In the case of small interaction strengths or low temperatures, our numerical results can be compared to the predictions of Bogoliubov theory and its perturbative extensions. This demonstrates the validity of the GPE in these limits and allows us to assign a temperature to the simulations unambiguously. However, the GPE method is nonperturbative, and we believe it can be used to describe the thermal properties of a Bose gas even when Bogoliubov theory fails. We suggest a different technique to measure the temperature of our simulations in these circumstances. Using this approach we determine the dependence of the condensate fraction and specific heat on temperature for several interaction strengths, and observe the appearance of vortex networks. Interesting behavior near the critical point is observed and discussed.

Exposure to automotive pollution increases plasma susceptibility to oxidation

Low-density lipoprotein oxidation is implicated in the development of atherosclerosis. Plasma susceptibility to oxidation may be used as a marker of low-density lipoprotein oxidation and thus predict atherosclerotic risk. In this study the authors investigated the relationship between plasma susceptibility to oxidation and exposure to automotive pollution in a group of automobile mechanics (n = 16) exposed to high levels of automotive pollution, vs. matched controls (n = 13). The authors induced plasma oxidation by a free radical initiator and they determined susceptibility to oxidation by (1) change in absorbance at 234 nm, (2) lag time to conjugated diene formation, and (3) linear slope of the oxidation curve. Mechanics had significantly higher values (mean standard error) for change in absorbance (1.60 +/- 0.05 vs. 1.36 +/- 0.05; p < .002), and slope (1.6 x 10(-3) +/- 0.1 x 10(-3) vs. 1.3 x 10(-3) +/- 0.1 x 10(-3); p < .001), compared with controls. These results indicate that regular exposure to automotive pollutants increases plasma susceptibility to oxidation and may, in the long term, increase the risk of developing atherosclerosis.

Database technologies for l-system simulations in virtual plant applications on bioinformatics

One of the most important advantages of database systems is that the underlying mathematics is rich enough to specify very complex operations with a small number of statements in the database language. This research covers an aspect of biological informatics that is the marriage of information technology and biology, involving the study of real-world phenomena using virtual plants derived from L-systems simulation. L-systems were introduced by Aristid Lindenmayer as a mathematical model of multicellular organisms. Not much consideration has been given to the problem of persistent storage for these simulations. Current procedures for querying data generated by L-systems for scientific experiments, simulations and measurements are also inadequate. To address these problems the research in this paper presents a generic process for data-modeling tools (L-DBM) between L-systems and database systems. This paper shows how L-system productions can be generically and automatically represented in database schemas and how a database can be populated from the L-system strings. This paper further describes the idea of pre-computing recursive structures in the data into derived attributes using compiler generation. A method to allow a correspondence between biologists' terms and compiler-generated terms in a biologist computing environment is supplied. Once the L-DBM gets any specific L-systems productions and its declarations, it can generate the specific schema for both simple correspondence terminology and also complex recursive structure data attributes and relationships.

Numerical simulations of wind and temperature structure within an Alpine lake basin, Lake Tekapo, New Zealand

High-resolution numerical model simulations have been used to study the local and mesoscale thermal circulations in an Alpine lake basin. The lake (87 km(2)) is situated in the Southern Alps, New Zealand and is located in a glacially excavated rock basin surrounded by mountain ranges that reach 3000 m in height. The mesoscale model used (RAMS) is a three-dimensional non-hydrostatic model with a level 2.5 turbulence closure scheme. The model demonstrates that thermal forcing at local (within the basin) and regional (coast-to-basin inflow) scales drive the observed boundary-layer airflow in the lake basin during clear anticyclonic summertime conditions. The results show that the lake can modify (perturb) both the local and regional wind systems. Following sunrise, local thermal circulations dominate, including a lake breeze component that becomes embedded within the background valley wind system. This results in a more divergent flow in the basin extending across the lake shoreline. However, a closed lake breeze circulation is neither observed nor modelled. Modelling results indicate that in the latter part of the day when the mesoscale (coast-to-basin) inflow occurs, the relatively cold pool of lake air in the basin can cause the intrusion to decouple from the surface. Measured data provide qualitative and quantitative support for the model results.

SUPPLY CHAIN RESILIENCE ANALYSIS: A BRAZILIAN AUTOMOTIVE CASE

Supply chain (SC) resilience and flexibility are important research topics receiving growing attention. However, the academic literature needs empirical studies on SC resilience capable of investigating the inter-organizational components of flexibility along different tiers. Therefore, this paper analyzes the main lack of flexibilities in three Brazilian automotive SCs that limit their resilience and therefore their capacity to better support and meet the demand changes in the marketplace. A multi-tier case study approach is adopted. Research findings identify lack of flexibilities in different tiers that inhibit the SC resilience as well as manufacturing and SC flexibilities that build SC resilience. The findings also highlight that the same SC may have the flexibility to be resilient for one of its products but not for another product, what sheds new lights on the academic literature. Finally, flexible SCs should be designed to increase SC resilience to cope with mishaps as significant demand changes.

Climate strategies of firms in the automotive and pulp & paper industries in Brazil: insights from an international perspective

This article focuses on the results of the final stage of research into the climate strategies of firms in the automotive and pulp-and-paper industries in Brazil, a country that is becoming increasingly important also in terms of climate change issues. In the first stage, the Climate Strategy Model (CSM) was developed to assess whether firms were adopting the necessary practices to assure the successful implementation of climate strategies. In the second, the CSM was applied to firms in the above mentioned industries that were chosen because of their important role in the domestic economy. In the final stage, interviews with executives of these firms were conducted to identify root causes of climate strategy implementation deficiencies and obtain new insights from an international perspective.