Relative impact of respiratory muscle activity on tidal flow and end expiratory volume in healthy neonates

INTRODUCTION: It has been suggested that infants dynamically regulate their tidal flow and end-expiratory volume level. The interaction between muscle activity, flow and lung volume in spontaneously sleeping neonates is poorly studied, since it requires the assessment of transcutaneous electromyography of respiratory muscles (rEMG) in matched comparison to lung function measurements. METHODS: After determining feasibility and repeatability of rEMG in 20 spontaneously sleeping healthy neonates, we measured the relative impact of intercostal and diaphragmatic EMG activity in direct comparison to the resulting tidal flow and FRC. RESULTS: We found good feasibility, repeatability and correlation of timing indices between rEMG activity and flow. The rEMG amplitude was significantly dependent on the resistive load of the face mask. Diaphragm and intercostal muscle activity commenced prior to the onset of flow and remained active during the expiratory cycle. The relative contribution of intercostal and diaphragmatic activity to flow was variable and changed dynamically. CONCLUSION: Using matched rEMG, air flow and lung volume measurements, we have found good feasibility and repeatability of intercostal and diaphragm rEMG measurements and provide the first quantitative measures of the temporal relationship between muscle activity and flow in spontaneously sleeping healthy neonates. Lung mechanical function is dynamically regulated and adapts on a breath to breath basis. So, non-invasive rEMG measurements alone or in combination with lung function might provide a more comprehensive picture of pulmonary mechanics in future studies. The data describing the timing of EMG and flow may be important for future studies of EMG triggered mechanical ventilation.

A patient-specific finite element methodology to predict damage accumulation in vertebral bodies under axial compression, sagittal flexion and combined loads

Due to the inherent limitations of DXA, assessment of the biomechanical properties of vertebral bodies relies increasingly on CT-based finite element (FE) models, but these often use simplistic material behaviour and/or single loading cases. In this study, we applied a novel constitutive law for bone elasticity, plasticity and damage to FE models created from coarsened pQCT images of human vertebrae, and compared vertebral stiffness, strength and damage accumulation for axial compression, anterior flexion and a combination of these two cases. FE axial stiffness and strength correlated with experiments and were linearly related to flexion properties. In all loading modes, damage localised preferentially in the trabecular compartment. Damage for the combined loading was higher than cumulated damage produced by individual compression and flexion. In conclusion, this FE method predicts stiffness and strength of vertebral bodies from CT images with clinical resolution and provides insight into damage accumulation in various loading modes.

Finite element 3D reconstruction of the pulmonary acinus imaged by synchrotron X-ray tomography

The alveolated structure of the pulmonary acinus plays a vital role in gas exchange function. Three-dimensional (3D) analysis of the parenchymal region is fundamental to understanding this structure-function relationship, but only a limited number of attempts have been conducted in the past because of technical limitations. In this study, we developed a new image processing methodology based on finite element (FE) analysis for accurate 3D structural reconstruction of the gas exchange regions of the lung. Stereologically well characterized rat lung samples (Pediatr Res 53: 72-80, 2003) were imaged using high-resolution synchrotron radiation-based X-ray tomographic microscopy. A stack of 1,024 images (each slice: 1024 x 1024 pixels) with resolution of 1.4 mum(3) per voxel were generated. For the development of FE algorithm, regions of interest (ROI), containing approximately 7.5 million voxels, were further extracted as a working subunit. 3D FEs were created overlaying the voxel map using a grid-based hexahedral algorithm. A proper threshold value for appropriate segmentation was iteratively determined to match the calculated volume density of tissue to the stereologically determined value (Pediatr Res 53: 72-80, 2003). The resulting 3D FEs are ready to be used for 3D structural analysis as well as for subsequent FE computational analyses like fluid dynamics and skeletonization.

Three dimensional finite element simulation of polymer melting and flow in a single-screw extruder : optimization of screw channel geometry

Single-screw extrusion is one of the widely used processing methods in plastics industry, which was the third largest manufacturing industry in the United States in 2007 [5]. In order to optimize the single-screw extrusion process, tremendous efforts have been devoted for development of accurate models in the last fifty years, especially for polymer melting in screw extruders. This has led to a good qualitative understanding of the melting process; however, quantitative predictions of melting from various models often have a large error in comparison to the experimental data. Thus, even nowadays, process parameters and the geometry of the extruder channel for the single-screw extrusion are determined by trial and error. Since new polymers are developed frequently, finding the optimum parameters to extrude these polymers by trial and error is costly and time consuming. In order to reduce the time and experimental work required for optimizing the process parameters and the geometry of the extruder channel for a given polymer, the main goal of this research was to perform a coordinated experimental and numerical investigation of melting in screw extrusion. In this work, a full three-dimensional finite element simulation of the two-phase flow in the melting and metering zones of a single-screw extruder was performed by solving the conservation equations for mass, momentum, and energy. The only attempt for such a three-dimensional simulation of melting in screw extruder was more than twenty years back. However, that work had only a limited success because of the capability of computers and mathematical algorithms available at that time. The dramatic improvement of computational power and mathematical knowledge now make it possible to run full 3-D simulations of two-phase flow in single-screw extruders on a desktop PC. In order to verify the numerical predictions from the full 3-D simulations of two-phase flow in single-screw extruders, a detailed experimental study was performed. This experimental study included Maddock screw-freezing experiments, Screw Simulator experiments and material characterization experiments. Maddock screw-freezing experiments were performed in order to visualize the melting profile along the single-screw extruder channel with different screw geometry configurations. These melting profiles were compared with the simulation results. Screw Simulator experiments were performed to collect the shear stress and melting flux data for various polymers. Cone and plate viscometer experiments were performed to obtain the shear viscosity data which is needed in the simulations. An optimization code was developed to optimize two screw geometry parameters, namely, screw lead (pitch) and depth in the metering section of a single-screw extruder, such that the output rate of the extruder was maximized without exceeding the maximum temperature value specified at the exit of the extruder. This optimization code used a mesh partitioning technique in order to obtain the flow domain. The simulations in this flow domain was performed using the code developed to simulate the two-phase flow in single-screw extruders.

Automated optimization of polymer extrusion dies using finite element analysis

An extrusion die is used to continuously produce parts with a constant cross section; such as sheets, pipes, tire components and more complex shapes such as window seals. The die is fed by a screw extruder when polymers are used. The extruder melts, mixes and pressures the material by the rotation of either a single or double screw. The polymer can then be continuously forced through the die producing a long part in the shape of the die outlet. The extruded section is then cut to the desired length. Generally, the primary target of a well designed die is to produce a uniform outlet velocity without excessively raising the pressure required to extrude the polymer through the die. Other properties such as temperature uniformity and residence time are also important but are not directly considered in this work. Designing dies for optimal outlet velocity variation using simple analytical equations are feasible for basic die geometries or simple channels. Due to the complexity of die geometry and of polymer material properties design of complex dies by analytical methods is difficult. For complex dies iterative methods must be used to optimize dies. An automated iterative method is desired for die optimization. To automate the design and optimization of an extrusion die two issues must be dealt with. The first is how to generate a new mesh for each iteration. In this work, this is approached by modifying a Parasolid file that describes a CAD part. This file is then used in a commercial meshing software. Skewing the initial mesh to produce a new geometry was also employed as a second option. The second issue is an optimization problem with the presence of noise stemming from variations in the mesh and cumulative truncation errors. In this work a simplex method and a modified trust region method were employed for automated optimization of die geometries. For the trust region a discreet derivative and a BFGS Hessian approximation were used. To deal with the noise in the function the trust region method was modified to automatically adjust the discreet derivative step size and the trust region based on changes in noise and function contour. Generally uniformity of velocity at exit of the extrusion die can be improved by increasing resistance across the die but this is limited by the pressure capabilities of the extruder. In optimization, a penalty factor that increases exponentially from the pressure limit is applied. This penalty can be applied in two different ways; the first only to the designs which exceed the pressure limit, the second to both designs above and below the pressure limit. Both of these methods were tested and compared in this work.

Applications of finite geometries to designs and codes

This dissertation concerns the intersection of three areas of discrete mathematics: finite geometries, design theory, and coding theory. The central theme is the power of finite geometry designs, which are constructed from the points and t-dimensional subspaces of a projective or affine geometry. We use these designs to construct and analyze combinatorial objects which inherit their best properties from these geometric structures. A central question in the study of finite geometry designs is Hamada’s conjecture, which proposes that finite geometry designs are the unique designs with minimum p-rank among all designs with the same parameters. In this dissertation, we will examine several questions related to Hamada’s conjecture, including the existence of counterexamples. We will also study the applicability of certain decoding methods to known counterexamples. We begin by constructing an infinite family of counterexamples to Hamada’s conjecture. These designs are the first infinite class of counterexamples for the affine case of Hamada’s conjecture. We further demonstrate how these designs, along with the projective polarity designs of Jungnickel and Tonchev, admit majority-logic decoding schemes. The codes obtained from these polarity designs attain error-correcting performance which is, in certain cases, equal to that of the finite geometry designs from which they are derived. This further demonstrates the highly geometric structure maintained by these designs. Finite geometries also help us construct several types of quantum error-correcting codes. We use relatives of finite geometry designs to construct infinite families of q-ary quantum stabilizer codes. We also construct entanglement-assisted quantum error-correcting codes (EAQECCs) which admit a particularly efficient and effective error-correcting scheme, while also providing the first general method for constructing these quantum codes with known parameters and desirable properties. Finite geometry designs are used to give exceptional examples of these codes.

Muffler characterization with implementation of the finite element method and experimental techniques

Determining how an exhaust system will perform acoustically before a prototype muffler is built can save the designer both a substantial amount of time and resources. In order to effectively use the simulation tools available it is important to understand what is the most effective tool for the intended purpose of analysis as well as how typical elements in an exhaust system affect muffler performance. An in-depth look at the available tools and their most beneficial uses are presented in this thesis. A full parametric study was conducted using the FEM method for typical muffler elements which was also correlated to experimental results. This thesis lays out the overall ground work on how to accurately predict sound pressure levels in the free field for an exhaust system with the engine properties included. The accuracy of the model is heavily dependent on the correct temperature profile of the model in addition to the accuracy of the source properties. These factors will be discussed in detail and methods for determining them will be presented. The secondary effects of mean flow, which affects both the acoustical wave propagation and the flow noise generation, will be discussed. Effective ways for predicting these secondary effects will be described. Experimental models will be tested on a flow rig that showcases these phenomena.

Erupted magma volume estimates at Santiaguito and Pacaya Volcanoes, Guatemala using digital elevation models

High resolution digital elevation models (DEMs) of Santiaguito and Pacaya volcanoes, Guatemala, were used to estimate volume changes and eruption rates between 1954 and 2001. The DEMs were generated from contour maps and aerial photography, which were analyzed in ArcGIS 9.0®. Because both volcanoes were growing substantially over the five decade period, they provide a good data set for exploring effective methodology for estimating volume changes. The analysis shows that the Santiaguito dome complex grew by 0.78 ± 0.07 km3 (0.52 ± 0.05 m3 s-1) over the 1954-2001 period with nearly all the growth occurring on the El Brujo (1958-75) and Caliente domes (1971-2001). Adding information from field data prior to 1954, the total volume extruded from Santiaguito since 1922 is estimated at 1.48 ± 0.19 km3. Santiaguito’s growth rate is lower than most other volcanic domes, but it has been sustained over a much longer period and has undergone a change toward more exogenous and progressively slower extrusion with time. At Santiaguito some of the material being added at the dome is subsequently transported downstream by block and ash flows, mudflows and floods, creating channel shifting and areas of aggradation and erosion. At Pacaya volcano a total volume of 0.21 ± 0.05 km3 was erupted between 1961 and 2001 for an average extrusion rate of 0.17 ± 0.04 m3 s-1. Both the Santiaguito and Pacaya eruption rate estimates reported here are minima, because they do not include estimates of materials which are transported downslope after eruption and data on ashfall which may result in significant volumes of material spread over broad areas. Regular analysis of high resolution DEMs using the methods outlined here, would help quantify the effects of fluvial changes to downstream populated areas, as well as assist in tracking hazards related to dome collapse and eruption.

Lung volume reduction surgery for severe emphysema

BACKGROUND: We report mid-term results after 25 consecutive lung volume reduction operations (LVRS) for the treatment of severe dyspnea due to advanced emphysema. METHODS: Study design: patients were studied prospectively up to 12 months after surgery. Setting: preoperative evaluation, surgery and postoperative care took place in our university hospital. Patients: patient selection was based on severe dyspnea and airway obstruction despite optimal medical treatment, lung overinflation and completed rehabilitation programme. Patients with severe hypercarbia (PCO2>50 mmHg) were excluded. Nineteen rehabilitated patients who fulfilled our inclusion criteria but postponed or denied LVRS were followed up clinically. Interventions: LVRS was performed bilaterally in 22 patients (median sternotomy) and unilaterally in 3 patients (limited thoracotomy). Measures: Outcome was measured by dyspnea evaluation, 6-minute-walking distance and pulmonary function tests. RESULTS: Twelve months postoperatively dyspnea and mobility improved significantly (MRC score from 3.3+/-0.7 to 2.12+/-0.8, 6-min-walk from 251+/-190 to 477+/-189 m). These results were superior compared to the results of the conservatively treated patients. Significant improvement could also be documented in airway obstruction (FEV1 from 960+/-369 to 1438+/-610 ml) and overinflation (TLC from 133+/-14 to 118+/-21% predicted and RV from 280+/-56 to 186+/-59% predicted). CONCLUSIONS: LVRS is an effective and promising treatment option for selected patients with end-stage emphysema and could be offered as an alternative and / or bridge to lung transplantation.

[Bilateral surgical lung volume reduction in severe emphysema]

OBJECTIVE: To report preliminary results with a new surgical method of treating terminal emphysema by bilateral reduction of lung volume. PATIENTS AND METHODS: In a prospective study, the results obtained with the first 20 consecutive patients (mean FEV1: 590 +/- 180 ml) who underwent operative reduction of lung volume were recorded. 19 of the 20 patients had required continuous oxygen supply. RESULTS: The patients were extubated 8.5 +/- 6 h postoperatively; thoracic drainage was removed after 9 +/- 6 days. The degree of dyspnoea was decreased in all patients (3.5 +/- 0.5 vs 0.5 +/- 0.1). Significant reduction of overinflation occurred soon after the operation (residual volume 273 +/- 125 to 201 +/- 107% of normal; total capacity from 142 +/- 18 to 109 +/- 22% of normal), as well as reduction in the degree of obstruction (FEV1 from 18 +/- 6 to 24 +/- 7% of normal; for each, P < 0.05). One patient died 3 weeks post-operatively of Candida infection. CONCLUSION: The method looks promising for the treatment of selected patients and may thus provide an alternative to lung transplantation.

Sirolimus-eluting stent treatment at high-volume centers confers lower mortality at 6-month follow-up: results from the prospective multicenter German Cypher Registry

BACKGROUND: Studies continue to identify percutaneous coronary intervention procedural volume both at the institutional level and at the operator level as being strongly correlated with outcome. High-volume centers have been defined as those that perform >400 percutaneous coronary intervention procedures per year. The relationship between drug-eluting stent procedural volume and outcome is unknown. We investigated this relationship in the German Cypher Registry. METHODS AND RESULTS: The present analysis included 8201 patients treated with sirolimus-eluting stents between April 2002 and September 2005 in 51 centers. Centers that recruited >400 sirolimus-eluting stent patients in this time period were considered high-volume centers; those with 150 to 400 patients were considered intermediate-volume centers; and those with <150 patients were designated as low-volume centers. The primary end point was all death, myocardial infarction, and target-vessel revascularization at 6 months. This end point occurred in 11.3%, 12.1%, and 9.0% of patients in the low-, intermediate-, and high-volume center groups, respectively (P=0.0001). There was no difference between groups in the rate of target-vessel revascularization (P=0.2) or cerebrovascular accidents (P=0.5). The difference in death/myocardial infarction remained significant after adjustment for baseline factors (odds ratio 1.85, 95% confidence interval 1.31 to 2.59, P<0.001 for low-volume centers; odds ratio 1.69, 95% confidence interval 1.29 to 2.21, P<0.001 for intermediate-volume centers). Patient and lesion selection, procedural features, and postprocedural medications differed significantly between groups. CONCLUSIONS: The volume of sirolimus-eluting stent procedures performed on an institutional level was inversely related to death and myocardial infarction but not to target-vessel revascularization at 6-month follow-up. Safety issues are better considered in high-volume centers. These findings have important public health policy implications.

Minimal local anesthetic volume for peripheral nerve block: a new ultrasound-guided, nerve dimension-based method

BACKGROUND AND OBJECTIVES: Nerve blocks using local anesthetics are widely used. High volumes are usually injected, which may predispose patients to associated adverse events. Introduction of ultrasound guidance facilitates the reduction of volume, but the minimal effective volume is unknown. In this study, we estimated the 50% effective dose (ED50) and 95% effective dose (ED95) volume of 1% mepivacaine relative to the cross-sectional area of the nerve for an adequate sensory block. METHODS: To reduce the number of healthy volunteers, we used a volume reduction protocol using the up-and-down procedure according to the Dixon average method. The ulnar nerve was scanned at the proximal forearm, and the cross-sectional area was measured by ultrasound. In the first volunteer, a volume of 0.4 mL/mm of nerve cross-sectional area was injected under ultrasound guidance in close proximity to and around the nerve using a multiple injection technique. The volume in the next volunteer was reduced by 0.04 mL/mm in case of complete blockade and augmented by the same amount in case of incomplete sensory blockade within 20 mins. After 3 up-and-down cycles, ED50 and ED95 were estimated. Volunteers and physicians performing the block were blinded to the volume used. RESULTS: A total 17 of volunteers were investigated. The ED50 volume was 0.08 mL/mm (SD, 0.01 mL/mm), and the ED95 volume was 0.11 mL/mm (SD, 0.03 mL/mm). The mean cross-sectional area of the nerves was 6.2 mm (1.0 mm). CONCLUSIONS: Based on the ultrasound measured cross-sectional area and using ultrasound guidance, a mean volume of 0.7 mL represents the ED95 dose of 1% mepivacaine to block the ulnar nerve at the proximal forearm.

The importance of growth hormone in the regulation of erythropoiesis, red cell mass, and plasma volume in adults with growth hormone deficiency

Total body water (TBW) is reduced in adult GH deficiency (GHD) largely due to a reduction of extracellular water. It is unknown whether total blood volume (TBV) contributes to the reduced extracellular water in GHD. GH and insulin-like growth factor I (IGF-I) have been demonstrated to stimulate erythropoiesis in vitro, in animal models, and in growing children. Whether GH has a regulatory effect on red cell mass (RCM) in adults is not known. We analyzed body composition by bioelectrical impedance and used standard radionuclide dilution methods to measure RCM and plasma volume (PV) along with measuring full blood count, ferritin, vitamin B12, red cell folate, IGF-I, IGF-binding protein-3, and erythropoietin in 13 adult patients with GHD as part of a 3-month, double blind, placebo-controlled trial of GH (0.036 U/kg.day). TBW and lean body mass significantly increased by 2.5 +/- 0.53 kg (mean +/- SEM; P < 0.004) and 3.4 +/- 0.73 kg (P < 0.004), respectively, and fat mass significantly decreased by 2.4 +/- 0.32 kg (P < 0.001) in the GH-treated group. The baseline RCM of all patients with GHD was lower than the predicted normal values (1635 +/- 108 vs. 1850 +/- 104 mL; P < 0.002). GH significantly increased RCM, PV, and TBV by 183 +/- 43 (P < 0.006), 350 +/- 117 (P < 0.03), and 515 +/- 109 (P < 0.004) mL, respectively. The red cell count increased by 0.36 +/- 0.116 x 10(12)/L (P < 0.03) with a decrease in ferritin levels by 39.1 +/- 4.84 micrograms/L (P < 0.001) after GH treatment. Serum IGF-I and IGF-binding protein-3 concentrations increased by 3.0 +/- 0.43 (P < 0.001) and 1.3 +/- 0.15 (P < 0.001) SD, respectively, but the erythropoietin concentration was unchanged after GH treatment. No significant changes in body composition or blood volume were recorded in the placebo group. Significant positive correlations could be established between changes in TBW and TBV, lean body mass and TBV (r = 0.78; P < 0.04 and r = 0.77; P < 0.04, respectively), and a significant negative correlation existed between changes in fat mass and changes in TBV in the GH-treated group (r = -0.95; P < 0.02). We conclude that 1) erythropoiesis is impaired in GHD; 2) GH stimulates erythropoiesis in adult GHD; and 3) GH increases PV and TBV, which may contribute to the increased exercise performance seen in these patients.

A Stereoscopic Volume Rendered Brain Atlas

The anatomy of the human brain is organized as a complex arrangement of interrelated structures in three dimensional space. To facilitate the understanding of both structure and function, we have created a volume rendered brain atlas (VRBA) with an intuitive interface that allows real-time stereoscopic rendering of brain anatomy. The VRBA incorporates 2-dimensional and 3-dimensional texture mapping to display segmented brain anatomy co-registered with a T1 MRI. The interface allows the user to remove and add any of the 62 brain structures, as well as control the display of the MRI dataset. The atlas also contains brief verbal and written descriptions of the different anatomical regions to correlate structure with function. A variety of stereoscopic projection methods are supported by the VRBA and provide an abstract, yet simple, way of visualizing brain anatomy and 3-dimensional relationships between different nuclei.