Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Water wave-driven seepage in marine sediments

The action of water waves moving over a porous seabed drives a seepage flux into and out of the marine sediments. The volume of fluid exchange per wave cycle may affect the rate of contaminant transport in the sediments. In this paper, the dynamic response of the seabed to ocean waves is treated analytically on the basis of pore-elastic theory applied to a porous seabed. The seabed is modelled as a semi-infinite, isotropic, homogeneous material. Most previous investigations on the wave-seabed interaction problem have assumed quasi-static conditions within the seabed, although dynamic behaviour often occurs in natural environments. Furthermore, wave pressures used in the previous approaches were obtained from conventional ocean wave theories: which are based on the assumption of an impermeable rigid seabed. By introducing a complex wave number, we derive a new wave dispersion equation, which includes the seabed characteristics (such as soil permeability, shear modulus, etc.). Based on the new closed-form analytical solution, the relative differences of the wave-induced seabed response under dynamic and quasi-static conditions are examined. The effects of wave and soil parameters on the seepage flux per wave cycle are also discussed in detail. (C) 2000 Elsevier Science Ltd. All rights reserved.

Impact on perinatal mortality of missed opportunities to treat maternal syphilis in rural South Africa: baseline results from a clinic randomized controlled trial

OBJECTIVE To demonstrate the impact on perinatal mortality of inadequate treatment for maternal syphilis despite adequate screening. METHOD In 12 clinics providing antenatal care in Hlabisa, South Africa 1783 pregnant women were screened for syphilis at their first antenatal visit between June and October 1998. Pregnancy outcome was determined among those with syphilis. RESULTS A total of 158 women were diagnosed with syphilis: prevalence 9% (95% CI 8-10%). Mean gestation at first antenatal visit was 24 weeks. Thirty women (19%) received no treatment and 96 (61%) received all three recommended doses of penicillin. Among those receiving at least one dose, mean delay to the first dose was 20 days. Among those fully treated mean delay to treatment completion was 34 days. Pregnancy outcome was known for 142 women (90%) and there were 17 perinatal deaths among 15 women (11%). Eleven of 43 women (26%) who received one or fewer doses of penicillin experienced ii perinatal death whilst only four of 99 women (4%) who received two or more doses of penicillin did so (P = 0.0001). Protection from perinatal death increased with the number of doses of penicillin: linear modelling suggests that one dose reduced the risk by 41%, two doses by 65% and three doses by 79%, compared with no doses. A dose-specific, categorical model confirmed reduction in risk by 79% for all three doses. CONCLUSION Despite effective screening, many pregnant women with syphilis remain inadequately treated, resulting in avoidable perinatal mortality. Delays in starting and finishing treatment, as well as incomplete treatment occur. Near-patient syphilis testing in the antenatal clinic with early treatment could improve treatment of syphilis and reduce perinatal mortality, and a randomized trial to test this is underway.

Free energy approximations in simple lattice proteins

This work addresses the question of whether it is possible to define simple pairwise interaction terms to approximate free energies of proteins or polymers. Rather than ask how reliable a potential of mean force is, one can ask how reliable it could possibly be. In a two-dimensional, infinite lattice model system one can calculate exact free energies by exhaustive enumeration. A series of approximations were fitted to exact results to assess the feasibility and utility of pairwise free energy terms. Approximating the true free energy with pairwise interactions gives a poor fit with little transferability between systems of different size. Adding extra artificial terms to the approximation yields better fits, but does not improve the ability to generalize from one system size to another. Furthermore, one cannot distinguish folding from nonfolding sequences via the approximated free energies. Most usefully, the methodology shows how one can assess the utility of various terms in lattice protein/polymer models. (C) 2001 American Institute of Physics.

Systemic administration of antisense p75(NTR) oligodeoxynucleotides rescues axotomised spinal motor neurons

The 75 kD low-affinity neurotrophin receptor (p75(NTR)) is expressed in developing and axotomised spinal motor neurons. There is now convincing evidence that p75NTR can, under some circumstances, become cytotoxic and promote neuronal cell death. We report here that a single application of antisense p75(NTR) oligodeoxynucleotides to the proximal nerve stumps of neonatal rats significantly reduces the loss of axotomised motor neurons compared to controls treated with nonsense oligodeoxynucleotides or phosphate-buffered saline. Our investigations also show that daily systemic intraperitoneal injections of antisense p75(NTR) oligodeoxynucleotides for 14 days significantly reduce the loss of axotomised motor neurons compared to controls. Furthermore, we found that systemic delivery over a similar period continues to be effective following axotomy when intraperitoneal injections were 1) administered after a delay of 24 hr, 2) limited to the first 7 days, or 3) administered every third day. In addition, p75(NTR) protein levels were reduced in spinal motor neurons following treatment with antisense p75(NTR) oligodeoxynucleotides. There were also no obvious side effects associated with antisense p75(NTR) oligodeoxynucleotide treatments as determined by behavioural observations and postnatal weight gain. Our findings indicate that antisense-based strategies could be a novel approach for the prevention of motor neuron degeneration associated with injuries or disease. (C) 2001 Wiley-Liss, Inc.

An overview of diversion strategies for Australian drug-related offenders

This paper provides a descriptive overview of options for diversion of drug-related offenders from the criminal justice system. Drug-related offences include drug offences (for example, possession of a prohibited substance); offences that are directly linked to intoxication (for example, drink-driving or assault); and offences committed to support drug use (for example, theft). After an offence has been detected by police, multiple opportunities for diversion occur throughout the criminal justice process. (a) Pre-arrest: when an offence is first detected, prior to a charge being laid. This is known as police diversion and includes fines, warnings and cautions, sometimes with educational information or referral to assessment and treatment. (b) Pre-trial: when a charge is made but before the matter is heard at court. Examples are treatment as a condition of bail, conferencing and prosecutor discretion. (c) Pre-sentence: a delay of sentence while assessment and treatment are sought. (d) Post-sentence: as part of sentencing, for example suspended sentences, drug courts, noncustodial sentences and circle sentencing. (e) Pre-release: prior to release from a sentence, on parole. Issues for diversion programmes include net widening, the ethics of coercion to treatment, the needs of families and intersectoral collaboration. A framework for diversion is presented in which increasingly treatment-focused and coercive diversion strategies are used as offenders' criminal careers and drug problems increase.

The q-deformed supersymmetric t-J model with a boundary

The q-deformed supersymmetric t-J model on a semi-infinite lattice is diagonalized by using the level-one vertex operators of the quantum affine superalgebra U-q[sl(2\1)]. We. give the bosonization of the boundary states. We give an integral expression for the correlation functions of the boundary model, and derive the difference equations which they satisfy.

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.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.

Inter-subband redundancy prediction using neural network for video coding

High performance video codec is mandatory for multimedia applications such as video-on-demand and video conferencing. Recent research has proposed numerous video coding techniques to meet the requirement in bandwidth, delay, loss and Quality-of-Service (QoS). In this paper, we present our investigations on inter-subband self-similarity within the wavelet-decomposed video frames using neural networks, and study the performance of applying the spatial network model to all video frames over time. The goal of our proposed method is to restore the highest perceptual quality for video transmitted over a highly congested network. Our contributions in this paper are: (1) A new coding model with neural network based, inter-subband redundancy (ISR) prediction for video coding using wavelet (2) The performance of 1D and 2D ISR prediction, including multiple levels of wavelet decompositions. Our result shows a short-term quality enhancement may be obtained using both 1D and 2D ISR prediction.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Controlling chaos in a Lorenz-like system using feedback

We demonstrate that the dynamics of an autonomous chaotic laser can be controlled to a periodic or steady state under self-synchronization. In general, past the chaos threshold the dependence of the laser output on feedback applied to the pump is submerged in the Lorenz-like chaotic pulsation. However there exist specific feedback delays that stabilize the chaos to periodic behavior or even steady state. The range of control depends critically on the feedback delay time and amplitude. Our experimental results are compared with the complex Lorenz equations which show good agreement.

A model of seller holdout

We model a buyer who wishes to combine objects owned by two separate sellers in order to realize higher value. Sellers are able to avoid entering into negotiations with the buyer, so that the order in which they negotiate is endogenous. Holdout occurs if at least one of the sellers is not present in the first round of negotiations. We demonstrate that complementarity of the buyer's technology is a necessary condition for equilibrium holdout. Moreover, a rise in complementarity leads to an increased likelihood of holdout, and an increased efficiency loss. Applications include patents, the land assembly problem, and mergers.

The land assembly problem revisited

As in the standard land assembly problem, a developer wants to buy two adjacent blocks of land belonging to two different owners. The value of the two blocks of land to the developer is greater than the sum of the individual values of the blocks for each owner. Unlike the land assembly literature, however, our focus is on the incentive that each lot owner has to delay the start of negotiations, rather than on the public goods nature of the problem. An incentive for delay exists, for example, when owners perceive that being last to sell will allow them to capture a larger share of the joint surplus from the development. We show that competition at point of sale can cause equilibrium delay, and that cooperation at point of sale will eliminate delay. This suggests that strategic delay is another source for the inefficient allocation of land, in addition to the public-good type externality pointed out by Grossman and Hart [Bell Journal of Economics 11 (1980) 42] and O'Flaherty [Regional Science and Urban Economics 24 (1994) 287]. (C) 2004 Elsevier B.V. All rights reserved.