Impulse-response function of splanchnic circulation with model-independent constraints: theory and experimental validation

Modeling physiological processes using tracer kinetic methods requires knowledge of the time course of the tracer concentration in blood supplying the organ. For liver studies, however, inaccessibility of the portal vein makes direct measurement of the hepatic dual-input function impossible in humans. We want to develop a method to predict the portal venous time-activity curve from measurements of an arterial time-activity curve. An impulse-response function based on a continuous distribution of washout constants is developed and validated for the gut. Experiments with simultaneous blood sampling in aorta and portal vein were made in 13 anesthetized pigs following inhalation of intravascular [O-15] CO or injections of diffusible 3-O[ C-11] methylglucose (MG). The parameters of the impulse-response function have a physiological interpretation in terms of the distribution of washout constants and are mathematically equivalent to the mean transit time ( T) and standard deviation of transit times. The results include estimates of mean transit times from the aorta to the portal vein in pigs: (T) over bar = 0.35 +/- 0.05 min for CO and 1.7 +/- 0.1 min for MG. The prediction of the portal venous time-activity curve benefits from constraining the regression fits by parameters estimated independently. This is strong evidence for the physiological relevance of the impulse-response function, which includes asymptotically, and thereby justifies kinetically, a useful and simple power law. Similarity between our parameter estimates in pigs and parameter estimates in normal humans suggests that the proposed model can be adapted for use in humans.

A unified model of flames as gasdynamic discontinuities

Viewed on a hydrodynamic scale, flames in experiments are often thin so that they may be described as gasdynamic discontinuities separating the dense cold fresh mixture from the light hot burned products. The original model of a flame as a gasdynamic discontinuity was due to Darrieus and to Landau. In addition to the fluid dynamical equations, the model consists of a flame speed relation describing the evolution of the discontinuity surface, and jump conditions across the surface which relate the fluid variables on the two sides of the surface. The Darrieus-Landau model predicts, in contrast to observations, that a uniformly propagating planar flame is absolutely unstable and that the strength of the instability grows with increasing perturbation wavenumber so that there is no high-wavenumber cutoff of the instability. The model was modified by Markstein to exhibit a high-wavenumber cutoff if a phenomenological constant in the model has an appropriate sign. Both models are postulated, rather than derived from first principles, and both ignore the flame structure, which depends on chemical kinetics and transport processes within the flame. At present, there are two models which have been derived, rather than postulated, and which are valid in two non-overlapping regions of parameter space. Sivashinsky derived a generalization of the Darrieus-Landau model which is valid for Lewis numbers (ratio of thermal diffusivity to mass diffusivity of the deficient reaction component) bounded away from unity. Matalon & Matkowsky derived a model valid for Lewis numbers close to unity. Each model has its own advantages and disadvantages. Under appropriate conditions the Matalon-Matkowsky model exhibits a high-wavenumber cutoff of the Darrieus-Landau instability. However, since the Lewis numbers considered lie too close to unity, the Matalon-Matkowsky model does not capture the pulsating instability. The Sivashinsky model does capture the pulsating instability, but does not exhibit its high-wavenumber cutoff. In this paper, we derive a model consisting of a new flame speed relation and new jump conditions, which is valid for arbitrary Lewis numbers. It captures the pulsating instability and exhibits the high-wavenumber cutoff of all instabilities. The flame speed relation includes the effect of short wavelengths, not previously considered, which leads to stabilizing transverse surface diffusion terms.

A modified Weibull distribution

A new lifetime distribution capable of modeling a bathtub-shaped hazard-rate function is proposed. The proposed model is derived as a limiting case of the Beta Integrated Model and has both the Weibull distribution and Type I extreme value distribution as special cases. The model can be considered as another useful 3-parameter generalization of the Weibull distribution. An advantage of the model is that the model parameters can be estimated easily based on a Weibull probability paper (WPP) plot that serves as a tool for model identification. Model characterization based on the WPP plot is studied. A numerical example is provided and comparison with another Weibull extension, the exponentiated Weibull, is also discussed. The proposed model compares well with other competing models to fit data that exhibits a bathtub-shaped hazard-rate function.

Explicit/implicit partitioning and a new explicit form of the generalized alpha method

Most finite element packages use the Newmark algorithm for time integration of structural dynamics. Various algorithms have been proposed to better optimize the high frequency dissipation of this algorithm. Hulbert and Chung proposed both implicit and explicit forms of the generalized alpha method. The algorithms optimize high frequency dissipation effectively, and despite recent work on algorithms that possess momentum conserving/energy dissipative properties in a non-linear context, the generalized alpha method remains an efficient way to solve many problems, especially with adaptive timestep control. However, the implicit and explicit algorithms use incompatible parameter sets and cannot be used together in a spatial partition, whereas this can be done for the Newmark algorithm, as Hughes and Liu demonstrated, and for the HHT-alpha algorithm developed from it. The present paper shows that the explicit generalized alpha method can be rewritten so that it becomes compatible with the implicit form. All four algorithmic parameters can be matched between the explicit and implicit forms. An element interface between implicit and explicit partitions can then be used, analogous to that devised by Hughes and Liu to extend the Newmark method. The stability of the explicit/implicit algorithm is examined in a linear context and found to exceed that of the explicit partition. The element partition is significantly less dissipative of intermediate frequencies than one using the HHT-alpha method. The explicit algorithm can also be rewritten so that the discrete equation of motion evaluates forces from displacements and velocities found at the predicted mid-point of a cycle. Copyright (C) 2003 John Wiley Sons, Ltd.

The tree cut and merge algorithm for estimation of network reliability

This article presents Monte Carlo techniques for estimating network reliability. For highly reliable networks, techniques based on graph evolution models provide very good performance. However, they are known to have significant simulation cost. An existing hybrid scheme (based on partitioning the time space) is available to speed up the simulations; however, there are difficulties with optimizing the important parameter associated with this scheme. To overcome these difficulties, a new hybrid scheme (based on partitioning the edge set) is proposed in this article. The proposed scheme shows orders of magnitude improvement of performance over the existing techniques in certain classes of network. It also provides reliability bounds with little overhead.

Association of severe coronary stenosis with subclinical left ventricular dysfunction in the absence of infarction

Background. Regional left ventricular (LV) dysfunction may occur in patients with coronary artery disease (CAD) in the absence of infarction, but the causes of this phenomenon are unclear. We sought to identify whether changes in regional LV function were related to stenosis severity, using sensitive new ultrasound markers of function. Methods: We studied 67 individuals with no history of infarction and with normal LV systolic function: 49 patients with CAD and 18 control subjects without CAD. All patients underwent color Doppler tissue imaging, integrated backscatter (IB), anatomic M-mode echocardiography, and strain rate imaging to detect changes in structure and function. Peak early and late diastolic myocardial velocity, cyclic variation of IB, wall thickness, and percent wall thickening were measured in each basal and mid segment. Strain rate and peak systolic strain were calculated in each wall. CAD was defined as greater than or equal to 50% diameter stenosis. Normokinetic segments (n = 354) subtended by CAD were divided according to stenosis severity into 3 groups: group 1 (subtended by 50%-69% stenosis); group 2 (subtended by 70%-98% stenosis); and group 3 (subtended by greater than or equal to99% stenosis). Each parameter in each group was compared with that in 216 segments from control subjects. Results: Segments subtended by significant CAD showed lower peak early and late diastolic myocardial velocity compared with control segments. Group 3 showed significantly lower myocardial velocities than group 2 for both peak early (4.8 +/- 1.8 vs 6.0 +/- 2.0 cm/s, P

Integrable generalized spin ladder models based on the SU(1|3) and SU(3|1) algebras

We present two integrable spin ladder models which possess a general free parameter besides the rung coupling J. The models are exactly solvable by means of the Bethe ansatz method and we present the Bethe ansatz equations. We analyze the elementary excitations of the models which reveal the existence of a gap for both models that depends on the free parameter. (C) 2003 American Institute of Physics.

Effect of vehicle pretreatment on the flux, retention, and diffusion of topically applied penetrants in vitro

Purpose. The flux of a topically applied drug depends on the activity in the skin and the interaction between the vehicle and skin. Permeation of vehicle into the skin can alter the activity of drug and the properties of the skin barrier. The aim of this in vitro study was to separate and quantify these effects. Methods. The flux of four radiolabeled permeants (water, phenol, diflunisal, and diazepam) with log K-oct/water values from 1.4 to 4.3 was measured over 4 h through heat-separated human epidermis pretreated for 30 min with vehicles having Hildebrand solubility parameters from 7.9 to 23.4 (cal/cm(3))(1/2). Results. Enhancement was greatest after pretreatment with the more lipophilic vehicles. A synergistic enhancement was observed using binary mixtures. The flux of diazepam was not enhanced to the same extent as the other permeants, possibly because its partitioning into the epidermis is close to optimal (log K-oct 2.96). Conclusion. An analysis of the permeant remaining in the epidermis revealed that the enhancement can be the result of either increased partitioning of permeant into the epidermis or an increasing diffusivity of permeants through the epidermis.

Is maximizing resilience compatible with established ecological goal functions?

Cropp and Gabric [Ecosystem adaptation: do ecosystems maximise resilience? Ecology. In press] used a simple phytoplanktonzooplankton-nutrient model and a genetic algorithm to determine the parameter values that would maximize the value of certain goal functions. These goal functions were to maximize biomass, maximize flux, maximize flux to biomass ratio, and maximize resilience. It was found that maximizing goal functions maximized resilience. The objective of this study was to investigate whether the Cropp and Gabric [Ecosystem adaptation: do ecosystems maximise resilience? Ecology. In press] result was indicative of a general ecosystem principle, or peculiar to the model and parameter ranges used. This study successfully replicated the Cropp and Gabric [Ecosystem adaptation: do ecosystems maximise resilience? Ecology. In press] experiment for a number of different model types, however, a different interpretation of the results is made. A new metric, concordance, was devised to describe the agreement between goal functions. It was found that resilience has the highest concordance of all goal functions trialled. for most model types. This implies that resilience offers a compromise between the established ecological goal functions. The parameter value range used is found to affect the parameter versus goal function relationships. Local maxima and minima affected the relationship between parameters and goal functions, and between goal functions. (C) 2003 Elsevier B.V. All rights reserved.

New fluorescence parameters for monitoring photosynthesis in plants - 1. The effect of illumination on the fluorescence parameters of the JIP-test

Chlorophyll fluorescence measurements have a wide range of applications from basic understanding of photosynthesis functioning to plant environmental stress responses and direct assessments of plant health. The measured signal is the fluorescence intensity (expressed in relative units) and the most meaningful data are derived from the time dependent increase in fluorescence intensity achieved upon application of continuous bright light to a previously dark adapted sample. The fluorescence response changes over time and is termed the Kautsky curve or chlorophyll fluorescence transient. Recently, Strasser and Strasser (1995) formulated a group of fluorescence parameters, called the JIP-test, that quantify the stepwise flow of energy through Photosystem II, using input data from the fluorescence transient. The purpose of this study was to establish relationships between the biochemical reactions occurring in PS II and specific JIP-test parameters. This was approached using isolated systems that facilitated the addition of modifying agents, a PS II electron transport inhibitor, an electron acceptor and an uncoupler, whose effects on PS II activity are well documented in the literature. The alteration to PS II activity caused by each of these compounds could then be monitored through the JIP-test parameters and compared and contrasted with the literature. The known alteration in PS II activity of Chenopodium album atrazine resistant and sensitive biotypes was also used to gauge the effectiveness and sensitivity of the JIP-test. The information gained from the in vitro study was successfully applied to an in situ study. This is the first in a series of four papers. It shows that the trapping parameters of the JIP-test were most affected by illumination and that the reduction in trapping had a run-on effect to inhibit electron transport. When irradiance exposure proceeded to photoinhibition, the electron transport probability parameter was greatly reduced and dissipation significantly increased. These results illustrate the advantage of monitoring a number of fluorescence parameters over the use of just one, which is often the case when the F-V/F-M ratio is used.

On the induced electric field gradients in the human body for magnetic stimulation by gradient coils in MRI

Prior theoretical studies indicate that the negative spatial derivative of the electric field induced by magnetic stimulation may he one of the main factors contributing to depolarization of the nerve fiber. This paper studies this parameter for peripheral nerve stimulation (PNS) induced by time.-varying gradient fields during MRI scans. The numerical calculations are based on an efficient, quasi-static, finite-difference scheme and an anatomically realistic human, full-body model. Whole-body cylindrical and planar gradient sets in MRI systems and various input signals have been explored. The spatial distributions of the induced electric field and their gradients are calculated and attempts are made to correlate these areas with reported experimental stimulation data. The induced electrical field pattern is similar for both the planar coils and cylindrical coils. This study provides some insight into the spatial characteristics of the induced field gradients for PNS in MRI, which may be used to further evaluate the sites where magnetic stimulation is likely to occur and to optimize gradient coil design.

Model-based analysis of anaerobic acetate uptake by a mixed culture of polyphosphate-accumulating and glycogen-accumulating organisms

An increasing number of studies shows that the glycogen-accumulating organisms (GAOs) can survive and may indeed proliferate under the alternating anaerobic/aerobic conditions found in EBPR systems, thus forming a strong competitor of the polyphosphate-accumulating organisms (PAOs). Understanding their behaviors in a mixed PAO and GAO culture under various operational conditions is essential for developing operating strategies that disadvantage the growth of this group of unwanted organisms. A model-based data analysis method is developed in this paper for the study of the anaerobic PAO and GAO activities in a mixed PAO and GAO culture. The method primarily makes use of the hydrogen ion production rate and the carbon dioxide transfer rate resulting from the acetate uptake processes by PAOs and GAOs, measured with a recently developed titration and off-gas analysis (TOGA) sensor. The method is demonstrated using the data from a laboratory-scale sequencing batch reactor (SBR) operated under alternating anaerobic and aerobic conditions. The data analysis using the proposed method strongly indicates a coexistence of PAOs and GAOs in the system, which was independently confirmed by fluorescent in situ hybridization (FISH) measurement. The model-based analysis also allowed the identification of the respective acetate uptake rates by PAOs and GAOs, along with a number of kinetic and stoichiometric parameters involved in the PAO and GAO models. The excellent fit between the model predictions and the experimental data not involved in parameter identification shows that the parameter values found are reliable and accurate. It also demonstrates that the current anaerobic PAO and GAO models are able to accurately characterize the PAO/GAO mixed culture obtained in this study. This is of major importance as no pure culture of either PAOs or GAOs has been reported to date, and hence the current PAO and GAO models were developed for the interpretation of experimental results of mixed cultures. The proposed method is readily applicable for detailed investigations of the competition between PAOs and GAOs in enriched cultures. However, the fermentation of organic substrates carried out by ordinary heterotrophs needs to be accounted for when the method is applied to the study of PAO and GAO competition in full-scale sludges. (C) 2003 Wiley Periodicals, Inc.

Effects of heat treatment on the nutritional value of raw soybean selected for low trypsin inhibitor activity

1. Two broiler experiments and a layer experiments were conducted on Kunitz trypsin inhibitor (Kti) soybeans (SB) of low trypsin inhibitor (TI) activity to determine their nutritive value when included as mash in least-cost poultry diets. 2. Experiment 1 compared chick performance on the Kti or raw SB using a commercial full-fat SB meal (FFSBM) and a solvent extracted SB meal (SBM) as controls during a 20 d experimental period. Broiler experiment 2 compared Kti and raw SB, non-steamed, or steam-pelleted with and without DL-methionine supplementation added to every treatment containing 170 g SB/kg. For each broiler experiment the levels of each SB were 70, 120 and 170 g/kg with the control birds fed only 170 g SB/kg. 3. The layer experiment, compared steam-pelleted Kti and raw SB against a non-steamed Kti and raw SB each fed at two levels (70 and 110 g/kg) x 30 replicates from 29 weeks of age for 19 weeks in a completely randomised design. Production parameters were measured when diets were formulated to contain minimum required specifications and calculated apparent metabolisable energy (AME). At the completion of each trial, 2 broiler birds from each cage and 5 layer birds per treatment were killed, weighed, and their liver and pancreas weighed. 4. Both broiler experiments indicated that production parameters on the Kti SB treatments were significantly lower (P < 0.05) than on the two commercial control SB treatments. However, the Kti treatments were superior to the raw SB treatments. 5. Pancreas weight increased with increasing inclusion of both raw and Kti SB, suggesting that a TI was causing the depression in performance. The AME of the Kti SB was similar to that of commercial FFSB meal. After steam conditioning, the raw SB meal AME value of 9.5 MJ/kg dry matter (DM) was improved to 14.1 MJ/kg DM by reduced TI activity, but this AME improvement with TI activity reduction, plus the supplementation with DL-methionine on birds fed the raw SB had no effect (P > 0.05) on any parameter evaluated in experiment 2. 6. The layer experiment showed that hens on the Kti SB treatments had significantly greater live weight gain (LWG), egg weight and daily egg mass than birds given raw SB. A reduced food intake (FI) was observed in the Kti treatments but egg mass was generally similar to that on the FFSB control diet, indicating that Kti SB supported excellent egg production at an inclusion of 110 g/kg. The depressed performance observed for broiler chicks suggest that younger birds are more susceptible to the effects of SB TI.

Stokes-operator-squeezed continuous-variable polarization states

We investigate nonclassical Stokes-operator variances in continuous-wave polarization-squeezed laser light generated from one and two optical parametric amplifiers. A general expression of how Stokes-operator variances decompose into two-mode quadrature operator variances is given. Stokes parameter variance spectra for four different polarization-squeezed states have been measured and compared with a coherent state. Our measurement results are visualized by three-dimensional Stokes-operator noise volumes mapped on the quantum Poincare sphere. We quantitatively compare the channel capacity of the different continuous-variable polarization states for communication protocols. It is shown that squeezed polarization states provide 33% higher channel capacities than the optimum coherent beam protocol.

Anomalous water absorption in porous materials

The absorption of fluid by unsaturated, rigid porous materials may be characterized by the sorptivity. This is a simple parameter to determine and is increasingly being used as a measure of a material's resistance to exposure to fluids (especially moisture and reactive solutes) in aggressive environments. The complete isothermal absorption process is described by a nonlinear diffusion equation, with the hydraulic diffusivity being a strongly nonlinear function of the degree of saturation of the material. This diffusivity can be estimated from the sorptivity test. In a typical test the cumulative absorption is proportional to the square root of time. However, a number of researchers have observed deviation from this behaviour when the infiltrating fluid is water and there is some potential for chemo-mechanical interaction with the material. In that case the current interpretation of the test and estimation of the hydraulic diffusivity is no longer appropriate. Kuntz and Lavallee (2001) discuss the anomalous behaviour and propose a non-Darcian model as a more appropriate physical description. We present an alternative Darcian explanation and theory that retrieves the earlier advantages of the simple sorptivity test in providing parametric information about the material's hydraulic properties and allowing simple predictive formulae for the wetting profile to be generated.