An assumption-driven case-specific model editor

A case sensitive intelligent model editor has been developed for constructing consistent lumped dynamic process models and for simplifying them using modelling assumptions. The approach is based on a systematic assumption-driven modelling procedure and on the syntax and semantics of process,models and the simplifying assumptions.

Tidal fluctuations in a leaky confined aquifer: Dynamic effects of an overlying phreatic aquifer

Tidal fluctuations in a leaky confined coastal aquifer are damped significantly due to leakage into an overlying phreatic aquifer. Jiao and Tang [1999] presented an analytical solution to a simple model describing this phenomenon. Their solution assumes that the tidal fluctuations in the overlying phreatic aquifer are negligible (i.e,, a static phreatic aquifer), Here we examine dynamic effects of the overlying aquifer based on a new approximate analytical solution. The numerical results indicate that the dynamic effects can be significant for a relatively large leakage and a high transmissivity of the phreatic aquifer.

Using stochastic dynamic programming to determine optimal fire management for Banksia ornata

1. A model of the population dynamics of Banksia ornata was developed, using stochastic dynamic programming (a state-dependent decision-making tool), to determine optimal fire management strategies that incorporate trade-offs between biodiversity conservation and fuel reduction. 2. The modelled population of B. ornata was described by its age and density, and was exposed to the risk of unplanned fires and stochastic variation in germination success. 3. For a given population in each year, three management strategies were considered: (i) lighting a prescribed fire; (ii) controlling the incidence of unplanned fire; (iii) doing nothing. 4. The optimal management strategy depended on the state of the B. ornata population, with the time since the last fire (age of the population) being the most important variable. Lighting a prescribed fire at an age of less than 30 years was only optimal when the density of seedlings after a fire was low (< 100 plants ha(-1)) or when there were benefits of maintaining a low fuel load by using more frequent fire. 5. Because the cost of management was assumed to be negligible (relative to the value of the persistence of the population), the do-nothing option was never the optimal strategy, although lighting prescribed fires had only marginal benefits when the mean interval between unplanned fires was less than 20-30 years.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.

An analysis of US greenhouse gas cap-and-trade proposals using a forward-looking economic model

We develop a forward-looking version of the recursive dynamic MIT Emissions Prediction and Policy Analysis (EPPA) model, and apply it to examine the economic implications of proposals in the US Congress to limit greenhouse gas (GHG) emissions. We find that shocks in the consumption path are smoothed out in the forward-looking model and that the lifetime welfare cost of GHG policy is lower than in the recursive model, since the forward-looking model can fully optimize over time. The forward-looking model allows us to explore issues for which it is uniquely well suited, including revenue-recycling and early action crediting. We find capital tax recycling to be more welfare-cost reducing than labor tax recycling because of its long-term effect on economic growth. Also, there are substantial incentives for early action credits; however, when spread over the full horizon of the policy they do not have a substantial effect on lifetime welfare costs.

Forward-looking versus recursive-dynamic modeling in climate policy analysis: A comparison

This paper develops a multi-regional general equilibrium model for climate policy analysis based on the latest version of the MIT Emissions Prediction and Policy Analysis (EPPA) model. We develop two versions so that we can solve the model either as a fully inter-temporal optimization problem (forward-looking, perfect foresight) or recursively. The standard EPPA model on which these models are based is solved recursively, and it is necessary to simplify some aspects of it to make inter-temporal solution possible. The forward-looking capability allows one to better address economic and policy issues such as borrowing and banking of GHG allowances, efficiency implications of environmental tax recycling, endogenous depletion of fossil resources, international capital flows, and optimal emissions abatement paths among others. To evaluate the solution approaches, we benchmark each version to the same macroeconomic path, and then compare the behavior of the two versions under a climate policy that restricts greenhouse gas emissions. We find that the energy sector and CO(2) price behavior are similar in both versions (in the recursive version of the model we force the inter-temporal theoretical efficiency result that abatement through time should be allocated such that the CO(2) price rises at the interest rate.) The main difference that arises is that the macroeconomic costs are substantially lower in the forward-looking version of the model, since it allows consumption shifting as an additional avenue of adjustment to the policy. On the other hand, the simplifications required for solving the model as an optimization problem, such as dropping the full vintaging of the capital stock and fewer explicit technological options, likely have effects on the results. Moreover, inter-temporal optimization with perfect foresight poorly represents the real economy where agents face high levels of uncertainty that likely lead to higher costs than if they knew the future with certainty. We conclude that while the forward-looking model has value for some problems, the recursive model produces similar behavior in the energy sector and provides greater flexibility in the details of the system that can be represented. (C) 2009 Elsevier B.V. All rights reserved.

Dynamic susceptibility investigation of maghemite nanoparticles incorporated in bovine serum albumin template

Room-temperature measurements of the magnetic susceptibility of Bovine Serum Albumin-based nanocapsules (50 to 300 nm in size) loaded with different amounts of maghemite nanoparticles (7.6 nm average diameter) have been carried out in this study The field (H) dependence of the imaginary peak susceptibility (f(P)) of the nanocomposite samples was investigated in the range of 0 to 4 kOe. From the analysis of the f(P) x H curves the concentration (N) dependence of the effective maghemite magnetocrystalline energy barrier (E) was obtained. Analysis of the E x N data was performed using a modified Morup-Tronc [Phys. Rev. Lett. 72, 3278 (1994)] model, from which a huge contribution from the magnetocrystalline surface anisotropy was observed.

Hyaluronidase Behavior at the Air/Liquid and Air/Lipid Interfaces and Improved Enzymatic Activity by Its Immobilization in a Biomembrane Model

Bovine testicular hyalurphidase (BT-HAase), a tetrameric enzyme responsible for randomly hyaluronic acid, catalytic hydrolysis, was successfully immobilized on Langmuir- Blodgett films prepared with the sodium salt of dihexadacylphosphoric acid, (DHP-Zn(II)) ending with dipalmitoylphosphatidylcholine, DPPC. Data of protein, adsorption at the air-liquid interface by means of pendant drop shipe analysis and interaction of the protein with Langmuir monolayers of DPPC, using a Langmuir trough, have provided information. about the conditions to be used in the protein immobilization. The dynamic surface pressure curves obtained from pendant drop experiments for the enzyme in buffer solutions indicate that, within the range of concentration investigated in this study, the enzyme exhibits the largest induction time at 5 mu g L(-1) attributed to diffusion processes. Nevertheless, it seems that, at this concentration, the most probable conformation should be the one which occupies the smallest area at pi -> 0. The surface pressure (pi) area curves obtained for BT-HAase and mixed DPPC- BT-HAase monolayers reveal the presence of the enzyme at the air-lipid interface up to 45 mN m(-1). Tests of enzymatic activity, using hyaluronic acid, HA, as the substrate, showed an increase of activity compared to the homogeneous medium. A simplified model of protein insertion into the lipid matrix is used to explain the obtained results.

A Comparison Between Pulse Pressure Variation and Right End Diastolic Volume Index as Guides to Resuscitation in a Model of Hemorrhagic Shock in Pigs

Background: Different hemodynamic parameters including static indicators of cardiac preload as right ventricular end-diastolic volume index (RVEDVI) and dynamic parameters as pulse pressure variation (PPV) have been used in the decision-making process regarding volume expansion in critically ill patients. The objective of this study was to compare fluid resuscitation guided by either PPV or RVEDVI after experimentally induced hemorrhagic shock. Methods: Twenty-six anesthetized and mechanically ventilated pigs were allocated into control (group I), PPV (group II), or RVEDVI (group III) group. Hemorrhagic shock was induced by blood withdrawal to target mean arterial pressure of 40 mm Hg, maintained for 60 minutes. Parameters were measured at baseline, time of shock, 60 minutes after shock, immediately after resuscitation with hydroxyethyl starch 6% (130/0.4), 1 hour and 2 hours thereafter. The endpoint of fluid resuscitation was determined as the baseline values of PPV and RVEDVI. Statistical analysis of data was based on analysis of variance for repeated measures followed by the Bonferroni test (p < 0.05). Results: Volume and time to resuscitation were higher in group III than in group II (group III = 1,305 +/- 331 mL and group II = 965 +/- 245 mL, p < 0.05; and group III = 24.8 +/- 4.7 minutes and group II = 8.8 +/- 1.3 minutes, p < 0.05, respectively). All static and dynamic parameters and biomarkers of tissue oxygenation were affected by hemorrhagic shock and nearly all parameters were restored after resuscitation in both groups. Conclusion: In the proposed model of hemorrhagic shock, resuscitation to the established endpoints was achieved within a smaller amount of time and with less volume when guided by PPV than when guided by pulmonary artery catheter-derived RVEDVI.

A fuzzy Reed-Frost model for epidemic spreading

In this paper, we present a fuzzy approach to the Reed-Frost model for epidemic spreading taking into account uncertainties in the diagnostic of the infection. The heterogeneities in the infected group is based on the clinical signals of the individuals (symptoms, laboratorial exams, medical findings, etc.), which are incorporated into the dynamic of the epidemic. The infectivity level is time-varying and the classification of the individuals is performed through fuzzy relations. Simulations considering a real problem with data of the viral epidemic in a children daycare are performed and the results are compared with a stochastic Reed-Frost generalization.

POLYETHYLENE GLYCOL ADDITION DOES NOT IMPROVE EXOGENOUS SURFACTANT FUNCTION IN AN EXPERIMENTAL MODEL OF MECONIUM ASPIRATION SYNDROME

Meconium (MEC) is a potent inactivator of pulmonary surfactant. The authors studied the effects of polyethylene glycol addition to the exogenous surfactant over the lung mechanics and volumes. Human meconium was administrated to newborn rabbits. Animals were ventilated for 20 minutes and dynamic compliance, ventilatory pressure, and tidal volume were recorded. Animals were randomized into 3 study groups: MEC group (without surfactant therapy); S100 group (100 mg/kg surfactant); and PEG group (100 mg/kg porcine surfactant plus 5% PEG). After ventilation, a pulmonary pressure-volume curve was built. Histological analysis was carried out to calculate the mean alveolar size (Lm) and the distortion index (DI). Both groups treated with surfactant showed higher values of dynamic pulmonary compliance and lower ventilatory pressure, compared with the MEC group (P .05). S100 group had a larger maximum lung volume, V30, compared with the MEC group (P .05). Lm and DI values were smaller in the groups treated with surfactant than in the MEC group (P .05). No differences were observed between the S100 and PEG groups. Animals treated with surfactant showed significant improvement in pulmonary function as compared to nontreated animals. PEG added to exogenous surfactant did not improve lung mechanics or volumes.

Comparative study of pressure- and volume-controlled ventilation on pulse pressure variation in a model of hypovolaemia in rabbits

Background and objective: Dynamic indices represented by systolic pressure variation and pulse pressure variation have been demonstrated to be more accurate than filling pressures in predicting fluid responsiveness. However, the literature is scarce concerning the impact of different ventilatory modes on these indices. We hypothesized that systolic pressure variation or pulse pressure variation could be affected differently by volume-controlled ventilation and pressure-controlled ventilation in an experimental model, during normovolaemia and hypovolaemia. Method: Thirty-two anaesthetized rabbits were randomly allocated into four groups according to ventilatory modality and volaemic status where G1-ConPCV was the pressure-controlled ventilation control group, G2-HemPCV was associated with haemorrhage, G3-ConVCV was the volume-controlled ventilation control group and G4-HemVCV was associated with haemorrhage. In the haemorrhage groups, blood was removed in two stages: 15% of the estimated blood volume withdrawal at M1, and, 30 min later, an additional 15% at M2. Data were submitted to analysis of variance for repeated measures; a value of P < 0.05 was considered to be statistically significant. Results: At MO (baseline), no significant differences were observed among groups. At M1, dynamic parameters differed significantly among the control and hypovolaemic groups (P < 0.05) but not between ventilation modes. However, when 30% of the estimated blood volume was removed (M2), dynamic parameters became significantly higher in animals under volume-controlled ventilation when compared with those under pressure-controlled ventilation. Conclusions: Under normovolaemia and moderate haemorrhage, dynamic parameters were not influenced by either ventilatory modalities. However, in the second stage of haemorrhage (30%), animals in volume-controlled ventilation presented higher values of systolic pressure variation and pulse pressure variation when compared with those submitted to pressure-controlled ventilation.

Combining two potential causes of metalloproteinase secretion causes abdominal aortic aneurysms in rats: a new experimental model

P>Progress in understanding the pathophysiology of abdominal aortic aneurysms (AAA) is dependent in part on the development and application of effective animal models that recapitulate key aspects of the disease. The objective was to produce an experimental model of AAA in rats by combining two potential causes of metalloproteinase (MMP) secretion: inflammation and turbulent blood flow. Male Wistar rats were randomly divided in four groups: Injury, Stenosis, Aneurysm and Control (40/group). The Injury group received a traumatic injury to the external aortic wall. The Stenosis group received an extrinsic stenosis at a corresponding location. The Aneurysm group received both the injury and stenosis simultaneously, and the Control group received a sham operation. Animals were euthanized at days 1, 3, 7 and 15. Aorta and/or aneurysms were collected and the fragments were fixed for morphologic, immunohistochemistry and morphometric analyses or frozen for MMP assays. AAAs had developed by day 3 in 60-70% of the animals, reaching an aortic dilatation ratio of more than 300%, exhibiting intense wall remodelling initiated at the adventitia and characterized by an obvious inflammatory infiltrate, mesenchymal proliferation, neoangiogenesis, elastin degradation and collagen deposition. Immunohistochemistry and zymography studies displayed significantly increased expressions of MMP-2 and MMP-9 in aneurysm walls compared to other groups. The haemo-dynamic alterations caused by the stenosis may have provided additional contribution to the MMPs liberation. This new model illustrated that AAA can be multifactorial and confirmed the key roles of MMP-2 and MMP-9 in this dynamic remodelling process.

Izergin-Korepin model with a boundary

The Izergin-Korepin model on a semi-infinite lattice is diagonalized by using the level-one vertex operators of the twisted quantum affine algebra U-q[((2))(2)]. We give the bosonization of the vacuum state with zero particle content. Excitation states are given by the action of the vertex operators on the vacuum state. We derive the boundary S-matrix. We give an integral expression of the correlation functions of the boundary model, and derive the difference equations which they satisfy. (C) 2001 Elsevier Science B.V. All rights reserved.

The acquisition of movement skills: Practice enhances the dynamic stability of bimanual coordination

During bimanual movements, two relatively stable inherent patterns of coordination (in-phase and anti-phase) are displayed (e.g., Kelso, Am. J. Physiol. 246 (1984) R1000). Recent research has shown that new patterns of coordination can be learned. For example, following practice a 90 degrees out-of-phase pattern can emerge as an additional, relatively stable, state (e.g., Zanone & Kelso, J. Exp. Psychol.: Human Performance and Perception 18 (1992) 403). On this basis, it has been concluded that practice leads to the evolution and stabilisation of the newly learned pattern and that this process of learning changes the entire attractor layout of the dynamic system. A general feature of such research has been to observe the changes of the targeted pattern's stability characteristics during training at a single movement frequency. The present study was designed to examine how practice affects the maintenance of a coordinated pattern as the movement frequency is scaled. Eleven volunteers were asked to perform a bimanual forearm pronation-supination task. Time to transition onset was used as an index of the subjects' ability to maintain two symmetrically opposite coordinated patterns (target task - 90 degrees out-of-phase - transfer task - 270 degrees out-of-phase). Their ability to maintain the target task and the transfer task were examined again after five practice sessions each consisting of 15 trials of only the 90 degrees out-of-phase pattern. Concurrent performance feedback (a Lissajous figure) was available to the participants during each practice trial. A comparison of the time to transition onset showed that the target task was more stable after practice (p = 0.025). These changes were still observed one week (p = 0.05) and two months (p = 0.075) after the practice period. Changes in the stability of the transfer task were not observed until two months after practice (p = 0.025). Notably, following practice, transitions from the 90 degrees pattern were generally to the anti-phase (180 degrees) pattern, whereas, transitions from the 270 degrees pattern were to the 90 degrees pattern. These results suggest that practice does improve the stability of a 90 degrees pattern, and that such improvements are transferable to the performance of the unpractised 270 degrees pattern. In addition, the anti-phase pattern remained more stable than the practised 90 degrees pattern throughout. (C) 2001 Elsevier Science B.V. All rights reserved.