Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Satellite observation and climate system model simulation of the St. Lawrence Island polynya

The St. Lawrence Island polynya (SLIP) is a commonly occurring winter phenomenon in the Bering Sea, in which dense saline water produced during new ice formation is thought to flow northward through the Bering Strait to help maintain the Arctic Ocean halocline. Winter darkness and inclement weather conditions have made continuous in situ and remote observation of this polynya difficult. However, imagery acquired from the European Space Agency ERS-1 Synthetic Aperture Radar (SAR) has allowed observation of the St. Lawrence Island polynya using both the imagery and derived ice displacement products. With the development of ARCSyM, a high resolution regional model of the Arctic atmosphere/sea ice system, simulation of the SLIP in a climate model is now possible. Intercomparisons between remotely sensed products and simulations can lead to additional insight into the SLIP formation process. Low resolution SAR, SSM/I and AVHRR infrared imagery for the St. Lawrence Island region are compared with the results of a model simulation for the period of 24-27 February 1992. The imagery illustrates a polynya event (polynya opening). With the northerly winds strong and consistent over several days, the coupled model captures the SLIP event with moderate accuracy. However, the introduction of a stability dependent atmosphere-ice drag coefficient, which allows feedbacks between atmospheric stability, open water, and air-ice drag, produces a more accurate simulation of the SLIP in comparison to satellite imagery. Model experiments show that the polynya event is forced primarily by changes in atmospheric circulation followed by persistent favorable conditions: ocean surface currents are found to have a small but positive impact on the simulation which is enhanced when wind forcing is weak or variable.

Simulation of mixing in unsteady flow through a periodically obstructed channel

Fluid mixing in steady and unsteady Bow through a channel containing periodic square obstructions has been studied using a finite-difference simulation to determine fluid velocities, followed by the use of passive marker particle advection to look at fluid transport out of the cavities formed between each of the obstructions. The geometry and Bow conditions were chosen from the work by Perkins (1989, M.S. Thesis, Lehigh University; 1992, Ph.D. Thesis, Lehigh University); who investigated heat transfer enhancement due to unsteady flow through such an obstructed channel. Particle advection shows that Bow regimes which are predicted to give good mixing based on snapshots of instantaneous streamline contour plots were not necessarily able to efficiently mix fluid which started in the cavity regions throughout the channel. The use of Poincare sections shows regular regions existing under these conditions which inhibit efficient fluid transport. These regular regions are found to disappear when the unsteady Bow velocity is increased. (C) 1997 Elsevier Science Ltd.

Gas-surface interactions in the thermal and sub-thermal regime: A molecular dynamics study

Molecular dynamics simulations are used to study the interaction of low-energy Ar atoms with the Ni(001) surface, Angular scattering distributions, in and out of the plane of incidence, are investigated as a function of incident energy, angles of incidence, crystallographic orientation of the incident beam and surface temperature. The results show a clear transition to the structure scattering regime at around 2 eV. However, at lower energies, two sub-regimes are revealed by the simulations, Far energies up to 250 meV, scattering is mainly diffuse, and significant trapping on the surface is observed, At energies above this level, lobular patterns start to form and trapping decreases with the increase in energy, Generally, there is a weak temperature dependence, but variations in the angle of incidence and/or changes in the crystallographic direction, generate significant changes in the scattering patterns.

Application of a wrist dosimeter prototype for radiation monitoring ((153)Sm) during a therapeutic procedure simulation

Gamma and beta radiation emitting radiopharmaceuticals are handled in nuclear medicine services, and in many cases there is only individual monitoring of gamma radiation. In this paper, the results obtained using a wrist dosimeter prototype (CaSO(4):Dy + Teflon pellets) show that the doses for workers occupationally exposed to beta radiation from (153)Sm are not negligible. It is important that this dose is evaluated, and it has to be taken into consideration in the individual monitoring system.

Virtual plants - integrating architectural and physiological models

Recent advances in computer technology have made it possible to create virtual plants by simulating the details of structural development of individual plants. Software has been developed that processes plant models expressed in a special purpose mini-language based on the Lindenmayer system formalism. These models can be extended from their architectural basis to capture plant physiology by integrating them with crop models, which estimate biomass production as a consequence of environmental inputs. Through this process, virtual plants will gain the ability to react to broad environmental conditions, while crop models will gain a visualisation component. This integration requires the resolution of the fundamentally different time scales underlying the approaches. Architectural models are usually based on physiological time; each time step encompasses the same amount of development in the plant, without regard to the passage of real time. In contrast, physiological models are based in real time; the amount of development in a time step is dependent on environmental conditions during the period. This paper provides a background on the plant modelling language, then describes how widely-used concepts of thermal time can be implemented to resolve these time scale differences. The process is illustrated using a case study. (C) 1997 Elsevier Science Ltd.

Non-equilibrium energy and momentum accommodation coefficients of Ar atoms scattered from Ni(001) in the thermal regime: A molecular dynamics study

Molecular dynamics simulations are used to study energy and momentum transfer of low-energy Ar atoms scattered from the Ni(001) surface. The investigation concentrates on the dependence of these processes on incident energy, angles of incidence and surface temperature. Energy transfer exhibits a strong dependence on the surface temperature, at incident energies below 500 meV, and incident angles close to specular incidence. Above 500 meV, the surface temperature dependence vanishes, and a limiting value in the amount of energy transferred to the surface is attained. Momentum exchange is investigated in terms of tangential and normal components. Both components exhibit a weak surface temperature dependence, but they have opposite behaviours at all incidence angles. In each component, momentum can be lost or gained following the interaction with the surface. (C) 1997 Elsevier Science B.V.

Improving wheat simulation capabilities in Australia from a cropping systems perspective: water and nitrogen effects on spring wheat in a semi-arid environment

Systems approaches can help to evaluate and improve the agronomic and economic viability of nitrogen application in the frequently water-limited environments. This requires a sound understanding of crop physiological processes and well tested simulation models. Thus, this experiment on spring wheat aimed to better quantify water x nitrogen effects on wheat by deriving some key crop physiological parameters that have proven useful in simulating crop growth. For spring wheat grown in Northern Australia under four levels of nitrogen (0 to 360 kg N ha(-1)) and either entirely on stored soil moisture or under full irrigation, kernel yields ranged from 343 to 719 g m(-2). Yield increases were strongly associated with increases in kernel number (9150-19950 kernels m(-2)), indicating the sensitivity of this parameter to water and N availability. Total water extraction under a rain shelter was 240 mm with a maximum extraction depth of 1.5 m. A substantial amount of mineral nitrogen available deep in the profile (below 0.9 m) was taken up by the crop. This was the source of nitrogen uptake observed after anthesis. Under dry conditions this late uptake accounted for approximately 50% of total nitrogen uptake and resulted in high (>2%) kernel nitrogen percentages even when no nitrogen was applied,Anthesis LAI values under sub-optimal water supply were reduced by 63% and under sub-optimal nitrogen supply by 50%. Radiation use efficiency (RUE) based on total incident short-wave radiation was 1.34 g MJ(-1) and did not differ among treatments. The conservative nature of RUE was the result of the crop reducing leaf area rather than leaf nitrogen content (which would have affected photosynthetic activity) under these moderate levels of nitrogen limitation. The transpiration efficiency coefficient was also conservative and averaged 4.7 Pa in the dry treatments. Kernel nitrogen percentage varied from 2.08 to 2.42%. The study provides a data set and a basis to consider ways to improve simulation capabilities of water and nitrogen effects on spring wheat. (C) 1997 Elsevier Science B.V.

Influence of thermal and mechanical cycling on the flexural strength of ceramics with titanium or gold alloy frameworks

Objectives. The aim of this study was to evaluate the effect of thermal and mechanical cycling alone or in combination, on the flexural strength of ceramic and metallic frameworks cast in gold alloy or titanium. Methods. Metallic frameworks (25 mm x 3 mm x 0.5 mm) (N = 96) cast in gold alloy or commercial pure titanium (Ti cp) were obtained using acrylic templates. They were airborne particle-abraded with 150 mu m aluminum oxide at the central area of the frameworks (8 mm x 3 mm). Bonding agent and opaque were applied on the particle-abraded surfaces and the corresponding ceramic for each metal was fired onto them. The thickness of the ceramic layer was standardized by positioning the frameworks in a metallic template (height: I mm). The specimens from each ceramic-metal combination (N = 96, n = 12 per group) were randomly assigned into four experimental fatigue conditions, namely water storage at 37 degrees C for 24 h (control group), thermal cycling (3000 cycles, between 4 and 55 degrees C, dwell time: 10 s), mechanical cycling (20,000 cycles under 10 N load, immersion in distilled water at 37 degrees C) and, thermal and mechanical cycling. A flexural strength test was performed in a universal testing machine (crosshead speed: 1.5 mm/min). Data were statistically analyzed using two-way ANOVA and Tukey`s test (alpha = 0.05). Results. The mean flexural strength values for the ceramic-gold alloy combination (55 +/- 7.2MPa) were significantly higher than those of the ceramic-Ti cp combination (32 +/- 6.7 MPa) regardless of the fatigue conditions performed (p < 0.05). Mechanical and thermo-mechanical fatigue decreased the flexural strength results significantly for both ceramic-gold alloy (52 +/- 6.6 and 53 +/- 5.6 MPa, respectively) and ceramic-Ti cp combinations (29 +/- 6.8 and 29 +/- 6.8 MPa, respectively) compared to the control group (58 +/- 7.8 and 39 SA MPa, for gold and Ti cp, respectively) (p < 0.05) (Tukey`s test). While ceramic-Ti cp combinations failed adhesively at the metal-opaque interface, gold alloy frameworks exhibited a residue of ceramic material on the surface in all experimental groups. Significance. Mechanical and thermo-mechanical fatigue conditions decreased the flexural strength values for both ceramic-gold alloy and ceramic-Ti cp combinations with the results being significantly lower for the latter in all experimental conditions. (C) 2007 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.