Finite element modelling of dissipative structures for nonequilibrium chemical reactions in fluid-saturated porous media

We use the finite element method to model and predict the dissipative structures of chemical species for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. In particular, we explore the conditions under which dissipative structures of the species may exist in the Brusselator type of nonequilibrium chemical reaction. Since this is the first time the finite element method and related strategies have been used to study the chemical instability problems in a fluid-saturated porous medium, it is essential to validate the method and strategies before they are put into application. For this purpose, we have rigorously derived the analytical solutions for dissipative structures of chemical species in a benchmark problem, which geometrically is a square. Comparison of the numerical solutions with the analytical ones demonstrates that the proposed numerical method and strategy are robust enough to solve chemical instability problems in a fluid-saturated porous medium. Finally, the related numerical results from two application examples indicate that both the regime and the magnitude of pore-fluid flow have significant effects on the nature of the dissipative structures that developed for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. The motivation for this study is that self-organization under conditions of pore-fluid flow in a porous medium is a potential mechanism of the orebody formation and mineralization in the upper crust of the Earth. (C) 2000 Elsevier Science S.A. All rights reserved.

Shock standoff on hypersonic blunt bodies in nonequilibrium gas flows

This paper presents a new theory of hypersonic blunt-nose shock standoff, based on a compressibility coordinate transformation for inviscid flow. It embraces a wide range of nonequilibrium shock-layer chemistry and gas mixtures including ionization and freestream dissociation. An extended binary scaling property of the analysis is also demonstrated. Specific application is made here to the family of arbitrarily diluted dissociating diatomic gases, with parametric study results presented for the scaled shock standoff distance as a function of an appropriate blunt-nose region Damkohler number. Comparisons with other theories and data in the case of nitrogen are also given and discussed.

Oceanic interchange and nonequilibrium population structure in the estuarine dependent Indo-Pacific tasselfish, Polynemus sheridani

We assayed mtDNA haplotype [300 base pairs (bp) control region] geography and genealogy in the Indo-Pacific tasselfish, Polynemus sheridani from its contiguous estuarine distribution across northern Australia (n = 169). Eight estuaries were sampled from three oceanographic regions (Timor Sea, Gulf of Carpentaria and the Coral Sea) to assess the impact of Pleistocene sea level changes on the historical connectivity among P. sheridani populations. Specifically, we investigated the genetic consequences of disruption to Indian-Pacific Ocean connectivity brought about by the closure of the Torres Strait. Overall there was significant population subdivision among estuaries (F-ST = 0.161, (Phi(ST) = 0.187). Despite a linear distribution, P. sheridani did not show isolation by distance over the entire sampled range because of genetic similarity of estuaries greater than 3000 km apart. However, significant isolation by distance was detected between estuaries separated by less than 3000 km of coastline. Unlike many genetic studies of Indo-Pacific marine species, there was no evidence for an historical division between eastern and western populations. Instead, phylogeographical patterns were dominated by a starlike intraspecific phylogeny coupled with evidence for population expansion in both the Gulf of Carpentaria and the Coral Sea but not the Timor Sea. This was interpreted as evidence for recent west to east recolonization across of northern Australia following the last postglacial marine advance. We argue that although sufficient time has elapsed postcolonization for populations to approach gene flow-drift equilibrium over smaller spatial scales (< 3000 km), the signal of historical colonization persists to obscure the expected equilibrium pattern of isolation by distance over large spatial scales (> 3000 km).

Expansion tube operating conditions for studying nonequilibrium radiation relevant to Titan Aerocapture

The predictions of nonequilibrium radiation in the shock layer for a Titan aerocapture aeroshell vary significantly amongst Computational Fluid Dynamics (CFD) analyses and are limited by the physical models of the nonequilibrium flow processes. Of particular interest are nonequilibrium processes associated with the CN molecule which is a strong radiator. It is necessary to have experimental data for these radiating shock layers which will allow for validation of the CFD models. This paper describes the development of a test flow condition for subscale aeroshell models in a superorbital expansion tunnel. We discuss the need for a Titan gas condition that closely simulates the atmospheric composition and present experimental data of the free stream test flow conditions. Furthermore, we present finite-rate CFD calculations of the facility to estimate the remaining free stream conditions, which cannot be directly measured during experiments.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Quantifying the influence of intra-aggregate concentration gradients on solute transport

The physical nonequilibrium of solute concentration resulting from preferential now of soil water has often led to models where the soil is partitioned into two regions: preferential flow paths, where solute transport occurs mainly by advection, and the remaining region, where significant solute transport occurs through diffusive exchange with the flow paths. These two-region models commonly ignore concentration gradients within the regions. Our objective was to develop a simple model to assess the influence of concentration gradients on solute transport and to compare model results with experiments conducted on structured materials. The model calculates the distribution of solutes in a single spherical aggregate surrounded by preferential now paths and subjected to alternating boundary conditions representing either an exchange of solutes between the two regions (a wet period) or no exchange but redistribution of solutes within the aggregate (a dry period). The key parameter in the model is the aggregate radius, which defines the diffusive time scales. We conducted intermittent leaching experiments on a column of packed porous spheres and on a large (300 mm long by 216 mm diameter) undisturbed field soil core to test the validity of the model and its application to field soils. Alternating wet and dry periods enhanced leaching by up to 20% for this soil, which was consistent with the model's prediction, given a fitted equivalent aggregate radius of 1.8 cm, If similar results are obtained for other soils, use of alternating wet and dry periods could improve management of solutes, for example in salinity control and in soil remediation.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

Quantum criticality

We investigate the theory of quantum fluctuations in non-equilibrium systems having large critical fluctuations. This allows us to treat the limits imposed by nonlinearities to quantum squeezing and noise reduction, and also to envisage future tests of quantum theory in regions of macroscopic quantum fluctuations. A long-term objective of this research is to identify suitable physical systems in which macroscopic 'Schrodinger cat'-like behaviour may be observed. We investigate two systems in particular of much current experimental interest, namely the degenerate parametric oscillator near threshold, and the evaporatively cooled (BEC). We compare the results obtained in the positive-P representation, as a fully quantum mechanical calculation, with the truncated Wigner phase space equation, also known as semi-classical theory. We show when these results agree and differ in calculations taken beyond the linearized approximation. In the region where the largest quantum fluctuations and Schrodinger cat-like behaviour might be expected, we find that the quantum predictions correspond very closely to the semi-classical theory. Nature abhors observing a Schrodinger car.

Inferring population history from microsatellite and enzyme data in serially introduced cane toads, Bufo marinus

Much progress has been made on inferring population history from molecular data. However, complex demographic scenarios have been considered rarely or have proved intractable. The serial introduction of the South-Central American cane Load Bufo marinas in various Caribbean and Pacific islands involves four major phases: a possible genetic admixture during the first introduction, a bottleneck associated with founding, a transitory, population boom, and finally, a demographic stabilization. A large amount of historical and demographic information is available for those introductions and can be combined profitably with molecular data. We used a Bayesian approach to combine this information With microsatellite (10 loci) and enzyme (22 loci) data and used a rejection algorithm to simultaneously estimate the demographic parameters describing the four major phases of the introduction history,. The general historical trends supported by microsatellites and enzymes were similar. However, there was a stronger support for a larger bottleneck at introductions for microsatellites than enzymes and for a more balanced genetic admixture for enzymes than for microsatellites. Verb, little information was obtained from either marker about the transitory population boom observed after each introduction. Possible explanations for differences in resolution of demographic events and discrepancies between results obtained with microsatellites and enzymes were explored. Limits Of Our model and method for the analysis of nonequilibrium populations were discussed.

Critical quantum fluctuations in the degenerate parametric oscillator

We develop a systematic theory of critical quantum fluctuations in the driven parametric oscillator. Our analytic results agree well with stochastic numerical simulations. We also compare the results obtained in the positive-P representation, as a fully quantum-mechanical calculation, with the truncated Wigner phase-space equation, also known as the semiclassical theory. We show when these results agree and differ in calculations taken beyond the linearized approximation. We find that the optimal broadband noise reduction occurs just above threshold. In this region where there are large quantum fluctuations in the conjugate variance and macroscopic quantum superposition states might be expected, we find that the quantum predictions correspond very closely to the semiclassical theory.

Limits to squeezing in the degenerate optical parametric oscillator

We develop a systematic theory of quantum fluctuations in the driven optical parametric oscillator, including the region near threshold. This allows us to treat the limits imposed by nonlinearities to quantum squeezing and noise reduction in this nonequilibrium quantum phase transition. In particular, we compute the squeezing spectrum near threshold and calculate the optimum value. We find that the optimal noise reduction occurs at different driving fields, depending on the ratio of damping rates. The largest spectral noise reductions are predicted to occur with a very high-Q second-harmonic cavity. Our analytic results agree well with stochastic numerical simulations. We also compare the results obtained in the positive-P representation, as a fully quantum-mechanical calculation, with the truncated Wigner phase-space equation, also known as the semiclassical theory.

The effect of volition on the stability of bimanual coordination

The authors investigated how the intention to passively perform a behavior and the intention to persist with a behavior impact upon the spatial and temporal properties of bimanual coordination. Participants (N = 30) were asked to perform a bimanual coordination task that demanded the continuous rhythmic extension-flexion of the wrists. The frequency of movement was scaled by an auditory metronome beat from 1.5 Hz, increasing to 3.25 Hz in .25-Hz increments. The task was further defined by the requirement that the movements be performed initially in a prescribed pattern of coordination (in-phase or antiphase) while the participants assumed one of two different intentional states: stay with the prescribed pattern should it become unstable or do not intervene should the pattern begin to change. Transitions away from the initially prescribed pattern were observed only in trials conducted in the antiphase mode of coordination. The time at which the antiphase pattern of coordination became unstable was not found to be influenced by the intentional state. In addition, the do-not-intervene set led to a switch to an in-phase pattern of coordination whereas the stay set led to phase wandering. Those findings are discussed within the framework of a dynamic account of bimanual coordination.