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.

Numerical simulation of dendrite growth during solidification

A comprehensive probabilistic model for simulating dendrite morphology and investigating dendritic growth kinetics during solidification has been developed, based on a modified Cellular Automaton (mCA) for microscopic modeling of nucleation, growth of crystals and solute diffusion. The mCA model numerically calculated solute redistribution both in the solid and liquid phases, the curvature of dendrite tips and the growth anisotropy. This modeling takes account of thermal, curvature and solute diffusion effects. Therefore, it can simulate microstructure formation both on the scale of the dendrite tip length. This model was then applied for simulating dendritic solidification of an Al-7%Si alloy. Both directional and equiaxed dendritic growth has been performed to investigate the growth anisotropy and cooling rate on dendrite morphology. Furthermore, the competitive growth and selection of dendritic crystals have also investigated.

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.

Direct measurement of pulse distortion near the zero-dispersion wavelength in an optical fiber by frequency-resolved optical gating

The technique of frequency-resolved optical gating is used to characterize the intensity and the phase of picosecond pulses after propagation through 700 m of fiber at close to the zero-dispersion wavelength. Using the frequency-resolved optical gating technique, we directly measure the severe temporal distortion resulting from the interplay between self-phase modulation and higher-order dispersion in this regime. The measured intensity and phase of the pulses after propagation are found to be in good agreement with the predictions of numerical simulations with the nonlinear Schrodinger equation. (C) 1997 Optical Society of America.

Local dispersion of the Eucalyptus leaf beetle Chrysophtharta bimaculata (Coleoptera: Chrysomelidae), and implications for forest protection

1. Chrysophtharta bimaculata is a native chrysomelid species that can cause chronic defoliation of plantation and regrowth Eucalyptus forests in Tasmania, Australia. Knowledge of the dispersion pattern of C. bimaculata was needed in order to assess the efficiency of an integrated pest management (IPM) programme currently used for its control. 2. Using data from yellow flight traps, local populations of C. bimaculata adults were monitored over a season at spatial scales relevant to commercial forestry: within a 50-ha operational management unit (a forestry 'coupe') and between coupes. In addition, oviposition was monitored over a season at a subset of the between-coupe sites. 3. Dispersion indices (Taylor's Power Law and Iwao's Mean Crowding regression method) demonstrated that C. bimaculata adults were spatially aggregated within and between coupes, although the number of egg-batches laid at the between-coupe scale was uniform. Spatial autocorrelation analysis showed that trap-catches at the within-coupe level were similar (positively autocorrelated) to a radius distance of approximately 110 m, and then dissimilar (negatively autocorrelated) at approximately 250 m. At the between-coupe scale, no repeatable spatial autocorrelation patterns were observed. 4. For any individual site, rapid changes in beetle density were observed to be associated with loosely aggregated flights of beetles into and out of that site. Peak adult catches (> the weekly mean plus standard deviation trap-catch) for a site occurred for a period of 2.0 +/- 0.22 weeks at a time (n = 37), with normally only one or two peaks per site per season. Peak oviposition events for a site occurred on average 1.4 +/- 0.11 times per season and lasted 1.5 +/- 0.12 weeks. 5. Analysis of an extensive data set (n = 417) demonstrated that adult abundance at a site was positively correlated with egg density, but negatively correlated with tree damage (caused by conspecifics) and the presence of conspecific larvae. There was no relationship between adult abundance and a visual estimate of the amount of young foliage on trees. 6. Adults of C. bimaculata are show n to occur in relatively small, mobile aggregations. This means that pest surveys must be both regular (less than 2 weeks apart) and intensive (with sampling points no more than 150 m apart) if beetle populations are to be monitored with confidence. Further refinement of the current IPM strategy must recognize the problems posed by this temporal and spatial patchiness, particularly with regard to the use of biological insecticides, such as Bacillus thuringiensis, for which only a very short operational window exists.

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.

On the validity of the dispersion model of hepatic drug elimination when intravascular transit time densities are long-tailed

The dispersion model with mixed boundary conditions uses a single parameter, the dispersion number, to describe the hepatic elimination of xenobiotics and endogenous substances. An implicit a priori assumption of the model is that the transit time density of intravascular indicators is approximated by an inverse Gaussian distribution. This approximation is limited in that the model poorly describes the tail part of the hepatic outflow curves of vascular indicators. A sum of two inverse Gaussian functions is proposed as ail alternative, more flexible empirical model for transit time densities of vascular references. This model suggests that a more accurate description of the tail portion of vascular reference curves yields an elimination rate constant (or intrinsic clearance) which is 40% less than predicted by the dispersion model with mixed boundary conditions. The results emphasize the need to accurately describe outflow curves in using them as a basis for determining pharmacokinetic parameters using hepatic elimination models. (C) 1997 Society for Mathematical Biology.

Relative dispersion of intravascular transit times during isolated human limb perfusions for recurrent melanoma

Aims We have characterized the relative dispersion of vascular and extravascular markers in the limbs of three patients undergoing isolated limb perfusions with the cytotoxic melphalan for recurrent malignant melanoma both before and after melphalan dosing. Methods A bolus of injectate containing [Cr-51] labelled red blood cells, [C-14]-sucrose and [H-3]-water was injected into an iliac or femoral artery and outflow samples collected at 1 s intervals by a fraction collector. The radioactivity due to each isotype was analysed by either gamma [Cr-51] or beta [C-14 and H-3] counting. The moments of the outflow fraction-time profiles were estimated by a nonparametric (numerical integration) method and a parametric model (sum of two inverse Gaussian functions). Results The availability, mean transit time and normalised variance (CV2) obtained for labelled red blood cells, sucrose and water were similar before and after melphalan dosing and with the two methods of calculation but varied between the patients. Conclusions The vascular space is not well-stirred but characterized by a CV2 similar that reported previously for in situ rat hind limb and rat liver perfusions. A flow-limited blood-tissue exchange was observed for the permeating indicators. Administration of melphalan did not influence the distribution characteristics of the indicators.

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.

Effects of oxide promoters on metal dispersion and metal-support interactions in Ni catalysts supported on activated carbon

Various oxide-promoted Ni catalysts supported on activated carbon were prepared, and the effect of promoters on the surface structure and properties of Ni catalysts was studied. Physical adsorption (Na adsorption), thermogravimetric analysis (TGA), temperature-programmed desorption (TPD), X-ray diffraction (XRD), and X-ray photoelectron spectroscopy (XPS) were used to characterize the catalysts. It is found that nickel is fairly uniformly distributed in the pores of the carbon support. Addition of promoters produces a more homogeneous distribution of nickel ion in carbon. However, distributions of promoters in the pores are varying. Addition of promoters increases the dispersion of nickel in carbon. Promoters also change the interaction between the carbon and Ni, resulting in significantly different behaviors of catalysts under various environments. CaO and MgO promoters improve the reactivity of nickel catalysts with O-2 but retard the interaction between nickel oxide and carbon. La2O3 shows some inhibiting effect on the interactions between nickel oxide and oxygen as well as carbon.

Slaking, dispersion, and crust formation

Soil structure is generally defined as the arrangement, orientation, and organization of the primary particles of sand, silt, and clay into compound aggregates, which exhibit properties that are unequal to the properties of a mass of nonaggregated material with a similar texture.6 Therefore the nature of soil structure is that it conveys specific properties to the soil and any alteration, i.e., breakdown or structural development, to the soil structural units will affect the physical properties of the soil. The aggregation and organization of the soil particles tend to form a hierarchical order4, 5 where the lower orders tend to have higher densities and greater internal strength than the higher orders. A schematic diagram of the hierarchical nature of soil structural elements in a clay soil is given in Fig. 1.4 Clay particles tend to form domains (packets of parallel clay sheets, generally consisting of 5-7 sheets), in turn several domains form clusters, followed by several orders of clusters, micro- and macroaggregates. The hierarchical nature implies that the destruction of a lower order will result in the destruction of all higher hierarchical orders. An example is the dispersion of sodic clay domains which results in the destruction of all higher orders, resulting in a dense soil with low hydraulic conductivity. Hence the clay domains are the fundamental building blocks of the soil and its integrity may determine the soil's physical properties and behavior.