Transfer matrix eigenvalues of the anisotropic multiparametric U model

A multiparametric extension of the anisotropic U model is discussed which maintains integrability. The R-matrix solving the Yang-Baxter equation is obtained through a twisting construction applied to the underlying U-q(sl (2/1)) superalgebraic structure which introduces the additional free parameters that arise in the model. Three forms of Bethe ansatz solution for the transfer matrix eigenvalues are given which we show to be equivalent.

Similitude applied to centrifugal scaling of unsaturated flow

Centrifuge experiments modeling single-phase flow in prototype porous media typically use the same porous medium and permeant. Then, well-known scaling laws are used to transfer the results to the prototype. More general scaling laws that relax these restrictions are presented. For permeants that are immiscible with an accompanying gas phase, model-prototype (i.e., centrifuge model experiment-target system) scaling is demonstrated. Scaling is shown to be feasible for Miller-similar (or geometrically similar) media. Scalings are presented for a more, general class, Lisle-similar media, based on the equivalence mapping of Richards' equation onto itself. Whereas model-prototype scaling of Miller-similar media can be realized easily for arbitrary boundary conditions, Lisle-similarity in a finite length medium generally, but not always, involves a mapping to a moving boundary problem. An exception occurs for redistribution in Lisle-similar porous media, which is shown to map to spatially fixed boundary conditions. Complete model-prototype scalings for this example are derived.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Stabilization of native protein fold by intein-mediated covalent cyclization

A mutant version of the N-terminal domain of Escherichia coli DnaB helicase was used as a model system to assess the stabilization against unfolding gained by covalent cyclization. Cyclization was achieved in vivo by formation of an amide bond between the N and C termini with the help of a split mini-intein. Linear and circular proteins were constructed to be identical in amino acid sequence. Mutagenesis of Phe102 to Glu rendered the protein monomeric even at high concentration. A difference in free energy of unfolding, DeltaDeltaG, between circular and linear protein of 2.3(+/-0.5) kcal mol(-1) was measured at 10degreesC by circular dichroism. A theoretical estimate of the difference in conformational entropy of linear and circular random chains in a three-dimensional cubic lattice model predicted DeltaDeltaG = 2.3 kcal mol(-1), suggesting that stabilization by protein cyclization is driven by the reduced conformational entropy of the unfolded state. Amide-proton exchange rates measured by NMR spectroscopy and mass spectrometry showed a uniform, approximately tenfold decrease of the exchange rates of the most slowly exchanging amide protons, demonstrating that cyclization globally decreases the unfolding rate of the protein. The amide proton exchange was found to follow EX1 kinetics at near-neutral pH, in agreement with an unusually slow refolding I measured by stopped-flow circular dichroism. rate of less than 4 min(-1) The linear and circular proteins differed more in their unfolding than in their folding rates. Global unfolding of the N-terminal domain of E. coli DnaB is thus promoted strongly by spatial separation of the N and C termini, whereas their proximity is much less important for folding. (C) 2005 Elsevier Ltd. All rights reserved.

Developmental vitamin D deficiency alters adult behaviour: An informative animal model of schizophrenia

Background: There is growing evidence that vitamin D is active in the brain but until recently there was a lack of evidence about its role during brain development. Guided by certain features of the epidemiology of schizophrenia, we have explored the role of vitamin D in the developing brain and behaviour using whole animal models. Methods: Sprague-Dawley rats were fed a vitamin D deficient diet (DVD) or control diet 6 weeks prior to mating and housed under UVB-free lighting conditions. On the day of birth all rats were fed a control diet for the remainder of the study. We observed behaviour at two timepoints; on the day of birth to study maternal behaviour, and at 10 weeks of age to study offspring behaviour in adulthood, under baseline and drug induced conditions (MK-801, haloperidol, amphetamine). Results: Prenatal vitamin D deficiency results in subtle alterations in maternal behaviour as well as long lasting effects on the adult offspring, despite a return to normal vitamin D levels during postnatal life. These affects were specific to transient prenatal vitamin D depletion as adult vitamin D depletion, combined prenatal and chronic postnatal vitamin D depletion, or ablation of the vitamin D receptor in mice led to markedly different outcomes. Conclusions: The developmental vitamin D (DVD) model now draws strength from epidemiological evidence of schizophrenia and animal experiments. Although the DVD model does not replicate every aspect of schizophrenia, it has several attractive features: (1) the exposure is based on clues from epidemiology; (2) it reproduces the increase in lateral ventricles; (3) it reproduces well-regarded behavioural phenotypes associated with schizophrenia (e.g. MK- 801 induced hyperlocomotion); and (4) it implicates a disturbance in dopamine signaling. In summary, low prenatal levels of vitamin D can influence critical components of orderly brain development and that this has a long lasting effect on behaviour.

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.

Algebraic Bethe ansatz for the anisotropic supersymmetric U model

We present an algebraic Bethe ansatz for the anisotropic supersymmetric U model for correlated electrons on the unrestricted 4(L)-dimensional electronic Hilbert space x(n=l)(L)C(4)(where L is the lattice length). The supersymmetry algebra of the local Hamiltonian is the quantum superalgebra U-q[gl(2\1)] and the model contains two symmetry-preserving free real parameters; the quantization parameter q and the Hubbard interaction parameter U. The parameter U arises from the one-parameter family of inequivalent typical four-dimensional irreps of U-q[gl(2\1)]. Eigenstates of the model are determined by the algebraic Bethe ansatz on a one-dimensional periodic lattice.

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.

Numerical modelling of tide-induced beach water table fluctuations

Field studies have shown that the elevation of the beach groundwater table varies with the tide and such variations affect significantly beach erosion or accretion. In this paper, we present a BEM (Boundary Element Method) model for simulating the tidal fluctuation of the beach groundwater table. The model solves the two-dimensional flow equation subject to free and moving boundary conditions, including the seepage dynamics at the beach face. The simulated seepage faces were found to agree with the predictions of a simple model (Turner, 1993). The advantage of the present model is, however, that it can be used with little modification to simulate more complicated cases, e.g., surface recharge from rainfall and drainage in the aquifer may be included (the latter is related to beach dewatering technique). The model also simulated well the field data of Nielsen (1990). In particular, the model replicated three distinct features of local water table fluctuations: steep rising phase versus flat falling phase, amplitude attenuation and phase lagging.

Isolated limb perfusion with melphalan for human melanoma xenografts in the hindlimb of nude rats: A surviving animal model

We have established a surviving model of isolated limb perfusion using xenografts of the human melanoma cell line MM 96L injected subcutaneously into the hindlimb of a nude rat, The femoral artery and vein were cannulated via the left renal artery and vein and the hind limb was isolated using tourniquets. The limb was perfused with Krebs Heinseleit buffer at 37 degrees C containing 4.7% bovine serum albumin at a constant flow rate of 4 mi per min for 30-60 min with 100% survival of the animals, Tumour vascularization and blood flow were demonstrated using vascular casts and [Cr-51]-microspheres. Following the addition of melphalan (15 or 100 mu g/ml), drug concentrations in the perfusate, tissues and systemic circulation were determined using high pressure liquid chromatography (HPLC), Systemic leakage, assessed using [I-125]albumin and melphalan and detected by a gamma-counter and HPLC respectively, was <0.5%. The melphalan concentration and tissue flow rate in the tumour deposits were 40 and 30% respectively, when compared with the surrounding subcutaneous tissue, At a dose of 15 mu g/ml, melphalan caused a reduction in tumour growth after 60 min perfusion, and a significant reduction in tumour size was seen when the melphalan dose was 100 mu g/ml. The surviving nude rat model of isolated limb perfusion for recurrent melanoma will allow examination of optimal perfusion conditions, along with the pharmacokinetics, pharmacodynamics and efficacy of melphalan and other drugs.

Hydraulic tracer studies in a pilot scale subsurface flow constructed wetland

Current design procedures for Subsurface Flow (SSF) Wetlands are based on the simplifying assumptions of plug flow and first order decay of pollutants. These design procedures do yield functional wetlands but result in over-design and inadequate descriptions of the pollutant removal mechanisms which occur within them. Even though these deficiencies are often noted, few authors have attempted to improve modelling of either flow or pollutant removal in such systems. Consequently the Oxley Creek Wetland, a pilot scale SSF wetland designed to enable rigorous monitoring, has recently been constructed in Brisbane, Australia. Tracer studies have been carried out in order to determine the hydraulics of this wetland prior to commissioning it with sealed sewage. The tracer studies will continue during the wetland's commissioning and operational phases. These studies will improve our understanding of the hydraulics of newly built SSF wetlands and the changes brought on by operational factors such as biological films and wetland plant root structures. Results to date indicate that the flow through the gravel beds is not uniform and cannot be adequately modelled by a single parameter, plug flow with dispersion, model. We have developed a multiparameter model, incorporating four plug flow reactors, which provides a better approximation of our experimental data. With further development this model will allow improvements to current SSF wetland design procedures and operational strategies, and will underpin investigations into the pollutant removal mechanisms at the Oxley Creek Wetland. (C) 1997 IAWQ. Published by Elsevier Science Ltd.

Melphalan dosing regimens for management of recurrent melanoma by isolated limb perfusion: Application of a physiological pharmacokinetic model based on melphalan distribution in the isolated perfused rat hindlimb

The optimal dosing schedule for melphalan therapy of recurrent malignant melanoma in isolated limb perfusions has been examined using a physiological pharmacokinetic model with data from isolated rat hindlimb perfusions (IRHP), The study included a comparison of melphalan distribution in IRHP under hyperthermia and normothermia conditions. Rat hindlimbs were perfused with Krebs-Henseleit buffer containing 4.7% bovine serum albumin at 37 or 41.5 degrees C at a flow rate of 4 ml/min. Concentrations of melphalan in perfusate and tissues were determined by high performance liquid chromatography with fluorescence detection, The concentration of melphalan in perfusate and tissues was linearly related to the input concentration. The rate and amount of melphalan uptake into the different tissues was higher at 41.5 degrees C than at 37 degrees C. A physiological pharmacokinetic model was validated from the tissue and perfusate time course of melphalan after melphalan perfusion. Application of the model involved the amount of melphalan exposure in the muscle, skin and fat in a recirculation system was related to the method of melphalan administration: single bolus > divided bolus > infusion, The peak concentration of melphalan in the perfusate was also related to the method of administration in the same order, Infusing the total dose of melphalan over 20 min during a 60 min perfusion optimized the exposure of tissues to melphalan whilst minimizing the peak perfusate concentration of melphalan. It is suggested that this method of melphalan administration may be preferable to other methods in terms of optimizing the efficacy of melphalan whilst minimizing the limb toxicity associated with its use in isolated limb perfusion.

A compact shock-assisted free-piston driver for impulse facilities

A free-piston driver that employs entropy-raising shock processes with diaphragm rupture has been constructed, which promises significant theoretical advantages over isentropic compression. Results from a range of conditions with helium and argon driver gases are reported. Significant performance gains were achieved in some test cases. Heat losses are shown to have a strong effect on driver processes. Measurements compare well with predictions from a quasi-one-dimensional numerical code.

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.

Interconnected-tubes model of hepatic elimination

The distributed-tubes model of hepatic elimination is extended to include intermixing between sinusoids, resulting in the formulation of a new, interconnected-tubes model. The new model is analysed for the simple case of two interconnected tubes, where an exact solution is obtained. For the case of many strongly-interconnected tubes, it is shown that a zeroth-order approximation leads to the convection-dispersion model. As a consequence the dispersion number is expressed, for the first time, in terms of its main physiological determinants: heterogeneity of flow and density of interconnections between sinusoids. The analysis of multiple indicator dilution data from a perfused liver preparation using the simplest version of the model yields the estimate 10.3 for the average number of interconnections. The problem of boundary conditions for the dispersion model is considered from the viewpoint that the dispersion-convection equation is a zeroth-order approximation to the equations for the interconnected-tubes model. (C) 1997 Academic Press Limited.