An extended computational model of the circulatory system for designing ventricular assist devices.

An extended computational model of the circulatory system has been developed to predict blood flow in the presence of ventricular assist devices (VADs). A novel VAD, placed in the descending aorta, intended to offload the left ventricle (LV) and augment renal perfusion is being studied. For this application, a better understanding of the global hemodynamic response of the VAD, in essence an electrically driven pump, and the cardiovascular system is necessary. To meet this need, a model has been established as a nonlinear, lumped-parameter electrical analog, and simulated results under different states [healthy, congestive heart failure (CHF), and postinsertion of VAD] are presented. The systemic circulation is separated into five compartments and the descending aorta is composed of three components to accurately yield the system response of each section before and after the insertion of the VAD. Delays in valve closing time and blood inertia in the aorta were introduced to deliver a more realistic model. Pump governing equations and optimization are based on fundamental theories of turbomachines and can serve as a practical initial design point for rotary blood pumps. The model's results closely mimic established parameters for the circulatory system and confirm the feasibility of the intra-aortic VAD concept. This computational model can be linked with models of the pump motor to provide a valuable tool for innovative VAD design.

Material model of lung parenchyma based on living precision-cut lung slice testing.

We describe a novel constitutive model of lung parenchyma, which can be used for continuum mechanics based predictive simulations. To develop this model, we experimentally determined the nonlinear material behavior of rat lung parenchyma. This was achieved via uni-axial tension tests on living precision-cut rat lung slices. The resulting force-displacement curves were then used as inputs for an inverse analysis. The Levenberg-Marquardt algorithm was utilized to optimize the material parameters of combinations and recombinations of established strain-energy density functions (SEFs). Comparing the best-fits of the tested SEFs we found Wpar = 4.1 kPa(I1-3)2 + 20.7 kPa(I1 - 3)3 + 4.1 kPa(-2 ln J + J2 - 1) to be the optimal constitutive model. This SEF consists of three summands: the first can be interpreted as the contribution of the elastin fibers and the ground substance, the second as the contribution of the collagen fibers while the third controls the volumetric change. The presented approach will help to model the behavior of the pulmonary parenchyma and to quantify the strains and stresses during ventilation.

A stochastic model of temporal variations in monthly temperature, precipitation, snowfall, and resulting snowpack

We report a Monte Carlo representation of the long-term inter-annual variability of monthly snowfall on a detailed (1 km) grid of points throughout the southwest. An extension of the local climate model of the southwestern United States (Stamm and Craig 1992) provides spatially based estimates of mean and variance of monthly temperature and precipitation. The mean is the expected value from a canonical regression using independent variables that represent controls on climate in this area, including orography. Variance is computed as the standard error of the prediction and provides site-specific measures of (1) natural sources of variation and (2) errors due to limitations of the data and poor distribution of climate stations. Simulation of monthly temperature and precipitation over a sequence of years is achieved by drawing from a bivariate normal distribution. The conditional expectation of precipitation. given temperature in each month, is the basis of a numerical integration of the normal probability distribution of log precipitation below a threshold temperature (3°C) to determine snowfall as a percent of total precipitation. Snowfall predictions are tested at stations for which long-term records are available. At Donner Memorial State Park (elevation 1811 meters) a 34-year simulation - matching the length of instrumental record - is within 15 percent of observed for mean annual snowfall. We also compute resulting snowpack using a variation of the model of Martinec et al. (1983). This allows additional tests by examining spatial patterns of predicted snowfall and snowpack and their hydrologic implications.