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.

Ladder operators for integrable one-dimensional lattice models

A generalised ladder operator is used to construct the conserved operators for any one-dimensional lattice model derived from the Yang-Baxter equation. As an example, the low order conserved operators for the XYh model are calculated explicitly.

Signature of mott-insulator transition with ultracold fermions in a one-dimensional optical lattice

Using the exact Bethe ansatz solution of the Hubbard model and Luttinger liquid theory, we investigate the density profiles and collective modes of one-dimensional ultracold fermions confined in an optical lattice with a harmonic trapping potential. We determine a generic phase diagram in terms of a characteristic filling factor and a dimensionless coupling constant. The collective oscillations of the atomic mass density, a technique that is commonly used in experiments, provide a signature of the quantum phase transition from the metallic phase to the Mott-insulator phase. A detailed experimental implementation is proposed.

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.

Ferromagnetism, paramagnetism, and a Curie-Weiss metal in an electron-doped Hubbard model on a triangular lattice

Motivated by the unconventional properties and rich phase diagram of NaxCoO2 we consider the electronic and magnetic properties of a two-dimensional Hubbard model on an isotropic triangular lattice doped with electrons away from half-filling. Dynamical mean-field theory (DMFT) calculations predict that for negative intersite hopping amplitudes (t < 0) and an on-site Coulomb repulsion, U, comparable to the bandwidth, the system displays properties typical of a weakly correlated metal. In contrast, for t > 0 a large enhancement of the effective mass, itinerant ferromagnetism, and a metallic phase with a Curie-Weiss magnetic susceptibility are found in a broad electron doping range. The different behavior encountered is a consequence of the larger noninteracting density of states (DOS) at the Fermi level for t > 0 than for t < 0, which effectively enhances the mass and the scattering amplitude of the quasiparticles. The shape of the DOS is crucial for the occurrence of ferromagnetism as for t > 0 the energy cost of polarizing the system is much smaller than for t < 0. Our observation of Nagaoka ferromagnetism is consistent with the A-type antiferromagnetism (i.e., ferromagnetic layers stacked antiferromagnetically) observed in neutron scattering experiments on NaxCoO2. The transport and magnetic properties measured in NaxCoO2 are consistent with DMFT predictions of a metal close to the Mott insulator and we discuss the role of Na ordering in driving the system towards the Mott transition. We propose that the Curie-Weiss metal phase observed in NaxCoO2 is a consequence of the crossover from a bad metal with incoherent quasiparticles at temperatures T > T-* and Fermi liquid behavior with enhanced parameters below T-*, where T-* is a low energy coherence scale induced by strong local Coulomb electron correlations. Our analysis also shows that the one band Hubbard model on a triangular lattice is not enough to describe the unusual properties of NaxCoO2 and is used to identify the simplest relevant model that captures the essential physics in NaxCoO2. We propose a model which allows for the Na ordering phenomena observed in the system which, we propose, drives the system close to the Mott insulating phase even at large dopings.

Pentoxifylline modifies three-dimensional collagen lattice model contraction and expression of collagen types I and III by human fibroblasts derived from post-burn hypertrophic scars and from normal skin

Fibroblasts are thought to be partially responsible for the persisting contractile forces that result in burn contractures. Using a monolayer cell culture and fibroblast populated collagen lattice (FPCL) three-dimensional model we subjected hypertrophic scar and non-cicatricial fibroblasts to the antifibrogenic agent pentoxifylline (PTF - 1 mg/mL) in order to reduce proliferation, collagen types I and III synthesis and model contraction. Fibroblasts were isolated from post-burn hypertrophic scars (HSHF) and non-scarred skin (NHF). Cells were grown in monolayers or incorporated into FPCL`s and exposed to PTF. In monolayer, cell number proliferation was reduced (46.35% in HSHF group and 37.73% in NHF group, p < 0.0001). PTF selectively inhibited collagen III synthesis in the HSHF group while inhibition was more evident to type I collagen synthesis in the NHF group. PTF also reduced contraction in both (HSHF and NHF) FPCL. (C) 2009 Elsevier Ltd and ISBI. 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.

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.

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.