Two dimensional computational fluid dynamic models for waste stabilisation ponds

Traditional waste stabilisation pond (WSP) models encounter problems predicting pond performance because they cannot account for the influence of pond features, such as inlet structure or pond geometry, on fluid hydrodynamics. In this study, two dimensional (2-D) computational fluid dynamics (CFD) models were compared to experimental residence time distributions (RTD) from literature. In one of the-three geometries simulated, the 2-D CFD model successfully predicted the experimental RTD. However, flow patterns in the other two geometries were not well described due to the difficulty of representing the three dimensional (3-D) experimental inlet in the 2-D CFD model, and the sensitivity of the model results to the assumptions used to characterise the inlet. Neither a velocity similarity nor geometric similarity approach to inlet representation in 2-D gave results correlating with experimental data. However. it was shown that 2-D CFD models were not affected by changes in values of model parameters which are difficult to predict, particularly the turbulent inlet conditions. This work suggests that 2-D CFD models cannot be used a priori to give an adequate description of the hydrodynamic patterns in WSP. (C) 1998 Elsevier Science Ltd. All rights reserved.

Metabolite mean transit times in the liver as predicted by various models of hepatic elimination

Predicted area under curve (AUC), mean transit time (MTT) and normalized variance (CV2) data have been compared for parent compound and generated metabolite following an impulse input into the liver, Models studied were the well-stirred (tank) model, tube model, a distributed tube model, dispersion model (Danckwerts and mixed boundary conditions) and tanks-in-series model. It is well known that discrimination between models for a parent solute is greatest when the parent solute is highly extracted by the liver. With the metabolite, greatest model differences for MTT and CV2 occur when parent solute is poorly extracted. In all cases the predictions of the distributed tube, dispersion, and tasks-in-series models are between the predictions of the rank and tube models. The dispersion model with mixed boundary conditions yields identical predictions to those for the distributed tube model (assuming an inverse gaussian distribution of tube transit times). The dispersion model with Danckwerts boundary conditions and the tanks-in series models give similar predictions to the dispersion (mixed boundary conditions) and the distributed tube. The normalized variance for parent compound is dependent upon hepatocyte permeability only within a distinct range of permeability values. This range is similar for each model but the order of magnitude predicted for normalized variance is model dependent. Only for a one-compartment system is the MIT for generated metabolite equal to the sum of MTTs for the parent compound and preformed metabolite administered as parent.

Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

Building better wildlife-habitat models

Wildlife-habitat models are an important tool in wildlife management toda?, and by far the majority of these predict aspects of species distribution (abundance or presence) as a proxy measure of habitat quality. Unfortunately, few are tested on independent data, and of those that are, few show useful predictive st;ill. We demonstrate that six critical assumptions underlie distribution based wildlife-habitat models, all of which must be valid for the model to predict habitat quality. We outline these assumptions in a mete-model, and discuss methods for their validation. Even where all sis assumptions show a high level of validity, there is still a strong likelihood that the model will not predict habitat quality. However, the meta-model does suggest habitat quality can be predicted more accurately if distributional data are ignored, and variables more indicative of habitat quality are modelled instead.

Phase diagram for a class of spin-1/2 Heisenberg models interpolating between the square-lattice, the triangular-lattice, and the linear-chain limits

We study the spin-1/2 Heisenberg models on an anisotropic two-dimensional lattice which interpolates between the square lattice at one end, a set of decoupled spin chains on the other end, and the triangular-lattice Heisenberg model in between. By series expansions around two different dimer ground states and around various commensurate and incommensurate magnetically ordered states, we establish the phase diagram for this model of a frustrated antiferromagnet. We find a particularly rich phase diagram due to the interplay of magnetic frustration, quantum fluctuations, and varying dimensionality. There is a large region of the usual two-sublattice Neel phase, a three-sublattice phase for the triangular-lattice model, a region of incommensurate magnetic order around the triangular-lattice model, and regions in parameter space where there is no magnetic order. We find that the incommensurate ordering wave vector is in general altered from its classical value by quantum fluctuations. The regime of weakly coupled chains is particularly interesting and appears to be nearly critical. [S0163-1829(99)10421-1].

Comparison of coherent and weakly incoherent transport models for the interlayer magnetoresistance of layered Fermi liquids

The interlayer magnetoresistance of layered metals in a tilted magnetic field is calculated for two distinct models for the interlayer transport. The first model involves coherent interlayer transport, and makes use of results of semiclassical or Bloch-Boltzmann transport theory. The second model involves weakly incoherent interlayer transport where the electron is scattered many times within a layer before tunneling into the next layer. The results are relevant to the interpretation of experiments on angular-dependent magnetoresistance oscillations (AMRO) in quasi-one- and quasi-two-dimensional organic metals. We find that the dependence of the magnetoresistance on the direction of the magnetic field is identical for both models except when the field is almost parallel to the layers. An important implication of this result is that a three-dimensional Fermi surface is not necessary for the observation of the Yamaji and Danner oscillations seen in quasi-two- and quasi-one-dimensional metals, respectively. A universal expression is given for the dependence of the resistance at AMRO maxima and minima on the magnetic field and scattering time (and thus the temperature). We point out three distinctive features of coherent interlayer transport: (i) a beat frequency in the magnetic oscillations of quasi-two-dimensional systems, (ii) a peak in the angular-dependent magnetoresistance when the field is sufficiently large and parallel to the layers, and (iii) a crossover from a linear to a quadratic field dependence for the magnetoresistance when the field is parallel to the layers. Properties (i) and (ii) are compared with published experimental data for a range of quasi-two-dimensional organic metals. [S0163-1829(99)02236-5].

Subgroup relations: A comparison of mutual intergroup differentiation and common ingroup identity models of prejudice reduction

Two studies examined relations between groups (humanities and math-science students) that implicitly or explicitly share a common superordinate category (university student). In Experiment 1, 178 participants performed a noninteractive decision-making task during which category salience was manipulated in a 2 (superordinate category salience) x 2 (subordinate category salience) between-groups design. Consistent with the mutual intergroup differentiation model, participants for whom both categories were salient exhibited the lowest levels of bias, whereas bias was strongest when the superordinate category alone was made salient. This pattern of results was replicated in Experiment 2 (N = 135). In addition, Experiment 2 demonstrated that members of subgroups that are nested within a superordinate category are more sensitive to how the superordinate category is represented than are members of subgroups that extend beyond the boundaries of the superordinate category.

Incorporating phase retardation effects into radiofrequency resonator models: application to magnetic resonance microscopy

A method is presented for including path propagation effects into models of radiofrequency resonators for use in magnetic resonance imaging. The method is based on the use of Helmholtz retarded potentials and extends our previous work on current density models of resonators based on novel inverse finite Hilbert transform solutions to the requisite integral equations. Radiofrequency phase retardation effects are most pronounced at high field strengths (frequencies) as are static field perturbations due to the magnetic materials in the resonators themselves. Both of these effects are investigated and a novel resonator structure presented for use in magnetic resonance microscopy.

Empirical catchment-wide rainfall erosivity models for two rivers in the humid tropics of Australia

The concept of rainfall erosivity is extended to the estimation of catchment sediment yield and its variation over time. Five different formulations of rainfall erosivity indices, using annual, monthly and daily rainfall data, are proposed and tested on two catchments in the humid tropics of Australia. Rainfall erosivity indices, using simple power functions of annual and daily rainfall amounts, were found to be adequate in describing the interannual and seasonal variation of catchment sediment yield. The parameter values of these rainfall erosivity indices for catchment sediment yield are broadly similar to those for rainfall erosivity models in relation to the R-factor in the Universal Soil Loss Equation.

Integrable Kondo impurities in one-dimensional extended Hubbard models

Three kinds of integrable Kondo problems in one-dimensional extended Hubbard models are studied by means of the boundary graded quantum inverse scattering method. The boundary K matrices depending on the local moments of the impurities are presented as a nontrivial realization of the graded reflection equation algebras acting in a (2s alpha + 1)-dimensional impurity Hilbert space. Furthermore, these models are solved using the algebraic Bethe ansatz method, and the Bethe ansatz equations are obtained.

Swapping space for time and unfair tests of ecological models

Testing ecological models for management is an increasingly important part of the maturation of ecology as an applied science. Consequently, we need to work at applying fair tests of models with adequate data. We demonstrate that a recent test of a discrete time, stochastic model was biased towards falsifying the predictions. If the model was a perfect description of reality, the test falsified the predictions 84% of the time. We introduce an alternative testing procedure for stochastic models, and show that it falsifies the predictions only 5% of the time when the model is a perfect description of reality. The example is used as a point of departure to discuss some of the philosophical aspects of model testing.

Evolution of stress deficit and changing rates of seismicity in cellular automaton models of earthquake faults

We investigate the internal dynamics of two cellular automaton models with heterogeneous strength fields and differing nearest neighbour laws. One model is a crack-like automaton, transferring ail stress from a rupture zone to the surroundings. The other automaton is a partial stress drop automaton, transferring only a fraction of the stress within a rupture zone to the surroundings. To study evolution of stress, the mean spectral density. f(k(r)) of a stress deficit held is: examined prior to, and immediately following ruptures in both models. Both models display a power-law relationship between f(k(r)) and spatial wavenumber (k(r)) of the form f(k(r)) similar tok(r)(-beta). In the crack model, the evolution of stress deficit is consistent with cyclic approach to, and retreat from a critical state in which large events occur. The approach to criticality is driven by tectonic loading. Short-range stress transfer in the model does not affect the approach to criticality of broad regions in the model. The evolution of stress deficit in the partial stress drop model is consistent with small fluctuations about a mean state of high stress, behaviour indicative of a self-organised critical system. Despite statistics similar to natural earthquakes these simplified models lack a physical basis. physically motivated models of earthquakes also display dynamical complexity similar to that of a critical point system. Studies of dynamical complexity in physical models of earthquakes may lead to advancement towards a physical theory for earthquakes.

Accelerating seismic energy release and evolution of event time and size statistics: Results from two heterogeneous cellular automaton models

The evolution of event time and size statistics in two heterogeneous cellular automaton models of earthquake behavior are studied and compared to the evolution of these quantities during observed periods of accelerating seismic energy release Drier to large earthquakes. The two automata have different nearest neighbor laws, one of which produces self-organized critical (SOC) behavior (PSD model) and the other which produces quasi-periodic large events (crack model). In the PSD model periods of accelerating energy release before large events are rare. In the crack model, many large events are preceded by periods of accelerating energy release. When compared to randomized event catalogs, accelerating energy release before large events occurs more often than random in the crack model but less often than random in the PSD model; it is easier to tell the crack and PSD model results apart from each other than to tell either model apart from a random catalog. The evolution of event sizes during the accelerating energy release sequences in all models is compared to that of observed sequences. The accelerating energy release sequences in the crack model consist of an increase in the rate of events of all sizes, consistent with observations from a small number of natural cases, however inconsistent with a larger number of cases in which there is an increase in the rate of only moderate-sized events. On average, no increase in the rate of events of any size is seen before large events in the PSD model.