Multilevel Domain Decomposition Algorithms for Monolithic Fluid-Structure Interaction Problems with Application to Haemodynamics

Finite element techniques for solving the problem of fluid-structure interaction of an elastic solid material in a laminar incompressible viscous flow are described. The mathematical problem consists of the Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian formulation coupled with a non-linear structure model, considering the problem as one continuum. The coupling between the structure and the fluid is enforced inside a monolithic framework which computes simultaneously for the fluid and the structure unknowns within a unique solver. We used the well-known Crouzeix-Raviart finite element pair for discretization in space and the method of lines for discretization in time. A stability result using the Backward-Euler time-stepping scheme for both fluid and solid part and the finite element method for the space discretization has been proved. The resulting linear system has been solved by multilevel domain decomposition techniques. Our strategy is to solve several local subproblems over subdomain patches using the Schur-complement or GMRES smoother within a multigrid iterative solver. For validation and evaluation of the accuracy of the proposed methodology, we present corresponding results for a set of two FSI benchmark configurations which describe the self-induced elastic deformation of a beam attached to a cylinder in a laminar channel flow, allowing stationary as well as periodically oscillating deformations, and for a benchmark proposed by COMSOL multiphysics where a narrow vertical structure attached to the bottom wall of a channel bends under the force due to both viscous drag and pressure. Then, as an example of fluid-structure interaction in biomedical problems, we considered the academic numerical test which consists in simulating the pressure wave propagation through a straight compliant vessel. All the tests show the applicability and the numerical efficiency of our approach to both two-dimensional and three-dimensional problems.

Windkanalstudien zum Graupelwachstum und Bestimmung von Retentionskoeffizienten verschiedener Spurengase während der Bereifung

Am vertikalen Windkanal der Johannes Gutenberg-Universität Mainz wurden physikalische und chemische Bereifungsexperimente durchgeführt. Dabei lagen die Umgebungstemperaturen bei allen Experimenten zwischen etwa -15 und -5°C und der Flüssigwassergehalt erstreckte sich von 0,9 bis etwa 1,6g/m³, typische Bedingungen für Mischphasenwolken in denen Bereifung stattfindet. Oberflächentemperaturmessungen an wachsenden hängenden Graupelpartikeln zeigten, dass während der Experimente trockene Wachstumsbedingungen herrschten.rnZunächst wurde das Graupelwachstum an in einer laminaren Strömung frei schwebenden Eispartikeln mit Anfangsradien zwischen 290 und 380µm, die mit flüssigen unterkühlten Wolkentröpfchen bereift wurden, studiert. Ziel war es, den Kollektionskern aus der Massenzunahme des bereiften Eispartikels und dem mittleren Flüssigwassergehalt während des Wachstumsexperimentes zu bestimmen. Die ermittelten Werte für die Kollektionskerne der bereiften Eispartikel erstreckten sich von 0,9 bis 2,3cm³/s in Abhängigkeit ihres Kollektorimpulses (Masse * Fallgeschwindigkeit des bereifenden Graupels), der zwischen 0,04 und 0,10gcm/s lag. Bei den Experimenten zeigte sich, dass die hier gemessenen Kollektionskerne höher waren im Vergleich mit Kollektionskernen flüssiger Tropfen untereinander. Aus den aktuellen Ergebnissen dieser Arbeit und der vorhandenen Literaturwerte wurde ein empirischer Faktor entwickelt, der von dem Wolkentröpfchenradius abhängig ist und diesen Unterschied beschreibt. Für die untersuchten Größenbereiche von Kollektorpartikel und flüssigen Tröpfchen können die korrigierten Kollektionskernwerte in Wolkenmodelle für die entsprechenden Größen eingebunden werden.rnBei den chemischen Experimenten zu dieser Arbeit wurde die Spurenstoffaufnahme verschiedener atmosphärischer Spurengase (HNO3, HCl, H2O2, NH3 und SO2) während der Bereifung untersucht. Diese Experimente mussten aus technischen Gründen mit hängenden Eispartikeln, dendritischen Eiskristallen und Schneeflocken, bereift mit flüssigen Wolkenlösungströpfchen, durchgeführt werden.rnDie Konzentrationen der Lösung, aus der die Wolkentröpfchen mit Hilfe von Zweistoffdüsen erzeugt wurden, lagen zwischen 1 und 120mg/l. Für die Experimente mit Ammoniak und Schwefeldioxid wurden Konzentrationen zwischen 1 und 22mg/l verwendet. Das Schmelzwasser der bereiften hängenden Graupel und Schneeflocken wurden ionenchromatographisch analysiert und zusammen mit der bekannten Konzentration der bereifenden Wolkentröpfchen konnte der Retentionskoeffizient für jeden Spurenstoff bestimmt werden. Er gibt die Menge an Spurenstoff an, die bei der Phasenumwandlung von flüssig zu fest in die Eisphase übergeht. Salpetersäure und Salzsäure waren nahezu vollständig retiniert (Mittelwerte der gesamten Experimente entsprechend 99±8% und 100±9%). Für Wasserstoffperoxid wurde ein mittlerer Retentionskoeffizient von 65±17% bestimmt. rnDer mittlere Retentionskoeffizient von Ammoniak ergab sich unabhängig vom Flüssigwassergehalt zu 92±21%, während sich für Schwefeldioxid 53±10% für niedrige und 29±7% für hohe Flüssigphasenkonzentrationen ergaben. Bei einigen der untersuchten Spurenstoffe wurde eine Temperaturabhängigkeit beobachtet und wenn möglich durch Parametrisierungen beschrieben.rn

Analysis of evening and morning transitions in complex terrain: MATERHORN campaign 2013

A field study of thermal circulation over very gentle slope is described for a specific day characterised by weak synoptic conditions. The emphasis was on morning and evening transitions, but measurements cover the entire day; therefore a brief analysis is performed to represent the general thermal circulation pattern. Both transition periods are characterised by complex dynamic behaviours. During evening transition, the upslope flow has got through a stagnation condition characterised by wind velocity U<0:5 m=s. Only when the stagnating air has become negative buoyant, the flow is allowed to pour downslope like a slab. Some features of front formation has been found during the transition development, such as delay time of downslope flow start up along the slope, and the presence of positive turbulent kinetic energy at the onset of the motion. Eventually the observed evening transition has followed a mixed mechanisms, with features from different models. Therefore the Rayleigh number seems not to be a good criterion by which parametrise evening transition itself. Morning transition is characterised by destruction of nocturnal temperature inversion and the onset of upslope flow. Inversion destruction can be described in terms of CBL growth at surface and inversion decent from the top of the layer. KH has found to be a good indicator of inversion breakup, if used as parameter to study the inversion breakup in terms of temperature reversal. After the inversion breakup, buoyancy and mechanical productions supply the flow with the necessary energy to start the upslope wind. More quantitative analysis are provided by the study of stability parameters and turbulent kinetic energy budgets. Gradient Richardson number has been used in this terms, finding that a mixed SBL-CBL behaviour dominates the most of the observed layers. Tke budget has shown high turbulent behaviour during morning transition while the evening transition has developed entirely in laminar condition, apart from short intermittent turbulent events.

Lung on Chip: In vitro HGF effects on injured alveolar A549 epithelial cells in microfluidic system

Microfluidic systems have become competitive tools in the invitro modelling of diseases and promising alternatives to animal studies. They allow obtaining more invivo like conditions for cellular assays. Research in idiopathic pulmonary fibrosis could benefit from this novel methodological approach to understand the pathophysiology of the disease & develop efficient therapies. The use of hepatocyte growth factor (HGF) for alveolar reepithelisation is a promising approach. In this study, we show a new microfluidic system to analyse the effects of HGF on injured alveolar epithelial cells. Microfluidic systems in polydimethylsiloxane were fabricated by soft lithography. The alveolar A549 epithelial cells (10,000 cells) were seeded and studied in these microfluidic systems with media perfusion (1μl/30min). Injury tests were made on the cells by the perfusion with media containing H2O2 or bleomycin. The degree of injury was then assessed by a metabolic and an apoptotic assays. Wound assays were also performed with a central laminar flow of trypsin. Monitoring of wound closure with HGF vs control media was assessed. The alveolar A549 epithelial cells grew and proliferated in the microfluidic system. In the wound closure assay, the degree of wound closure after 5 hours was (53.3±1.3%) with HGF compared to (9.8±2.4%) without HGF (P <0.001). We present a novel microfluidic model that allows culture, injury and wounding of A549 epithelial cells and represents the first step towards the development of an invitro reconstitution of the alveolar-capillary interface. We were also able to confirm that HGF increased alveolar epithelial repair in this system.

Continuous Formation and Separation of Calcium-Alginate Capsules Containing Viable Mammalian Cells in Microfluidic Devices

Microfluidic devices can be used for many applications, including the formation of well-controlled emulsions. In this study, the capability to continuously create monodisperse droplets in a microfluidic device was used to form calcium-alginate capsules.Calcium-alginate capsules have many potential uses, such as immunoisolation of cells and microencapsulation of active drug ingredients or bitter agents in food or beverage products. The gelation of calcium-alginate capsules is achieved by crosslinking sodiumalginate with calcium ions. Calcium ions dissociated from calcium carbonate due to diffusion of acetic acid from a sunflower oil phase into an aqueous droplet containing sodium-alginate and calcium carbonate. After gelation, the capsules were separated from the continuous oil phase into an aqueous solution for use in biological applications. Typically, capsules are separated bycentrifugation, which can damage both the capsules and the encapsulated material. A passive method achieves separation without exposing the encapsulated material or the capsules to large mechanical forces, thereby preventing damage. To achieve passiveseparation, the use of a microfluidic device with opposing channel wa hydrophobicity was used to stabilize co-laminar flow of im of hydrophobicity is accomplished by defining one length of the channel with a hydrogel. The chosen hydrogel was poly (ethylene glycol) diacrylate, which adheres to the glass surface through the use of self-assembled monolayer of 3-(trichlorosilyl)-propyl methacrylate. Due to the difference in surface energy within the channel, the aqueous stream is stabilized near a hydrogel and the oil stream is stabilized near the thiolene based optical adhesive defining the opposing length of the channel. Passive separation with co-laminar flow has shown success in continuously separating calcium-alginatecapsules from an oil phase into an aqueous phase. In addition to successful formation and separation of calcium alginate capsules,encapsulation of Latex micro-beads and viable mammalian cells has been achieved. The viability of encapsulated mammalian cells was determined using a live/dead stain. The co-laminar flow device has also been demonstrated as a means of separating liquid-liquidemulsions.

Efficient Liquid Liquid Extraction By Emulsion Formation and Separation in Robust Milli-Fluidic Devices

Conventional liquid liquid extraction (LLE) methods require large volumes of fluids to achieve the desired mass transfer of a solute, which is unsuitable for systems dealing with a low volume or high value product. An alternative to these methods is to scale down the process. Millifluidic devices share many of the benefits of microfluidic systems, including low fluid volumes, increased interfacial area-to-volume ratio, and predictability. A robust millifluidic device was created from acrylic, glass, and aluminum. The channel is lined with a hydrogel cured in the bottom half of the device channel. This hydrogel stabilizes co-current laminar flow of immiscible organic and aqueous phases. Mass transfer of the solute occurs across the interface of these contacting phases. Using a y-junction, an aqueous emulsion is created in an organic phase. The emulsion travels through a length of tubing and then enters the co-current laminar flow device, where the emulsion is broken and each phase can be collected separately. The inclusion of this emulsion formation and separation increases the contact area between the organic and aqueous phases, therefore increasing the area over which mass transfer can occur. Using this design, 95% extraction efficiency was obtained, where 100% is represented by equilibrium. By continuing to explore this LLE process, the process can be optimized and with better understanding may be more accurately modeled. This system has the potential to scale up to the industrial level and provide the efficient extraction required with low fluid volumes and a well-behaved system.

Hemodynamic Modeling of the Intrarenal Circulation

Three dimensional, time dependent numerical simulations of healthy and pathological conditions in a model kidney were performed. Blood flow in a kidney is not commonly investigated by computational approach, in contrast for example, to the flow in a heart. The flow in a kidney is characterized by relatively small Reynolds number (100 < Re < 0.01-laminar regime). The presented results give insight into the structure of such flow, which is hard to measure in vivo. The simulations have suggested that venous thrombosis is more likely than arterial thrombosis-higher shear rate observed. The obtained maximum velocity, as a result of the simulations, agrees with the observed in vivo measurements. The time dependent simulations show separation regimes present in the vicinity of the maximum pressure value. The pathological constriction introduced to the arterial geometry leads to the changes in separation structures. The constriction of a single vessel affects flow in the whole kidney. Pathology results in different flow rate values in healthy and affected branches, as well as, different pulsate cycle characteristic for the whole system.

Alterations of the thoracic spine in Marfan's syndrome

OBJECTIVE: The purpose of this study was to determine if the thoracic vertebral elements are altered in patients with Marfan's syndrome. MATERIALS AND METHODS: Thirty patients underwent helical CT of the thorax because of suspected thoracic aortic dilatation and acute dissection. Thirteen had Marfan's syndrome and 17 did not. Two reviewers, unaware of the final diagnosis, evaluated the images by consensus for laminar thickness, foraminal width, dural sac ratios, and vertebral scalloping for T2-T12. RESULTS: At T9-T12, dural sac ratios at the midcorpus level (p = 0.031) and foraminal width (p = 0.0124) were significantly greater in the patients with Marfan's syndrome than in the patients without. Dural sac ratios at lower endplate levels (p = 0.0685), laminar thickness (p = 0.951), and vertebral scalloping (p = 0.24) were not significantly greater in the patients with Marfan's syndrome than in the patients without. CONCLUSION: Because the phenotypic expression of Marfan's syndrome is variable, information on the spine from thoracic studies in combination with major criteria may be helpful clinically.

Multiple thromboembolic events in fetofetal transfusion syndrome in triplets contributing to the understanding of pathogenesis of hydranencephaly in combination with polymicrogyria

Over the last 180 years, several theories concerning the origin of hydranencephaly have been proposed with an emphasis on infectious, aplastic, and vascular etiologies. In this report, we present a case of triplets with fetofetal transfusion syndrome of which 2 fetuses (1 and 2) developed almost similar hydranencephaly, whereas the third exhibited the features of a fetus papyraceus (3). In the monochorial triamniotic placenta, multiple arteriovenous anastomoses were detected, representing a probable route for the transmission of thrombi originating from fetus 3 causing visceral lesions in fetus 2. Hydranencephaly was histologically characterized by necrosis, macrophage invasion, and endothelial proliferation. In addition, polymicrogyria was seen in fetuses 1 and 2. The combination of multiple visceral thromboembolic events and the death of fetus 3 approximately in the 11th week of gestation suggested a vascular thrombotic pathogenesis of hydranencephaly. Polymicrogyria can be considered as postmigratory laminar necrosis. Our findings contribute to the pathogenetic understanding of combined hydranencephaly and polymicrogyria.

Development of one-dimensional and two-dimensional computational tools that simulate steady internal condensing flows in terrestrial and zero-gravity environments

This dissertation presents an effective quasi one-dimensional (1-D) computational simulation tool and a full two-dimensional (2-D) computational simulation methodology for steady annular/stratified internal condensing flows of pure vapor. These simulation tools are used to investigate internal condensing flows in both gravity as well as shear driven environments. Through accurate numerical simulations of the full two dimensional governing equations, results for laminar/laminar condensing flows inside mm-scale ducts are presented. The methodology has been developed using MATLAB/COMSOL platform and is currently capable of simulating film-wise condensation for steady (and unsteady flows). Moreover, a novel 1-D solution technique, capable of simulating condensing flows inside rectangular and circular ducts with different thermal boundary conditions is also presented. The results obtained from the 2-D scientific tool and 1-D engineering tool, are validated and synthesized with experimental results for gravity dominated flows inside vertical tube and inclined channel; and, also, for shear/pressure driven flows inside horizontal channels. Furthermore, these simulation tools are employed to demonstrate key differences of physics between gravity dominated and shear/pressure driven flows. A transition map that distinguishes shear driven, gravity driven, and “mixed” driven flow zones within the non-dimensional parameter space that govern these duct flows is presented along with the film thickness and heat transfer correlations that are valid in these zones. It has also been shown that internal condensing flows in a micro-meter scale duct experiences shear driven flow, even in different gravitational environments. The full 2-D steady computational tool has been employed to investigate the length of annularity. The result for a shear driven flow in a horizontal channel shows that in absence of any noise or pressure fluctuation at the inlet, the onset of non-annularity is partly due to insufficient shear at the liquid-vapor interface. This result is being further corroborated/investigated by R. R. Naik with the help of the unsteady simulation tool. The condensing flow results and flow physics understanding developed through these simulation tools will be instrumental in reliable design of modern micro-scale and spacebased thermal systems.

Computational study of internal and external condensing flows and experimental synthesis to investigate their attainability and stability in ground-based and space-based environments

This doctoral thesis presents the computational work and synthesis with experiments for internal (tube and channel geometries) as well as external (flow of a pure vapor over a horizontal plate) condensing flows. The computational work obtains accurate numerical simulations of the full two dimensional governing equations for steady and unsteady condensing flows in gravity/0g environments. This doctoral work investigates flow features, flow regimes, attainability issues, stability issues, and responses to boundary fluctuations for condensing flows in different flow situations. This research finds new features of unsteady solutions of condensing flows; reveals interesting differences in gravity and shear driven situations; and discovers novel boundary condition sensitivities of shear driven internal condensing flows. Synthesis of computational and experimental results presented here for gravity driven in-tube flows lays framework for the future two-phase component analysis in any thermal system. It is shown for both gravity and shear driven internal condensing flows that steady governing equations have unique solutions for given inlet pressure, given inlet vapor mass flow rate, and fixed cooling method for condensing surface. But unsteady equations of shear driven internal condensing flows can yield different “quasi-steady” solutions based on different specifications of exit pressure (equivalently exit mass flow rate) concurrent to the inlet pressure specification. This thesis presents a novel categorization of internal condensing flows based on their sensitivity to concurrently applied boundary (inlet and exit) conditions. The computational investigations of an external shear driven flow of vapor condensing over a horizontal plate show limits of applicability of the analytical solution. Simulations for this external condensing flow discuss its stability issues and throw light on flow regime transitions because of ever-present bottom wall vibrations. It is identified that laminar to turbulent transition for these flows can get affected by ever present bottom wall vibrations. Detailed investigations of dynamic stability analysis of this shear driven external condensing flow result in the introduction of a new variable, which characterizes the ratio of strength of the underlying stabilizing attractor to that of destabilizing vibrations. Besides development of CFD tools and computational algorithms, direct application of research done for this thesis is in effective prediction and design of two-phase components in thermal systems used in different applications. Some of the important internal condensing flow results about sensitivities to boundary fluctuations are also expected to be applicable to flow boiling phenomenon. Novel flow sensitivities discovered through this research, if employed effectively after system level analysis, will result in the development of better control strategies in ground and space based two-phase thermal systems.

Molecular Modeling of PMR-15 Polyimide

PMR-15 polyimide is a polymer that is used as a matrix in composites. These composites with PMR-15 matrices are called advanced polymer matrix composite that is abundantly used in the aerospace and electronics industries because of its high temperature resistivity. Apart from having high temperature sustainability, PMR-15 composites also display good thermal-oxidative stability, mechanical properties, processability and low costs, which makes it a suitable material for manufacturing aircraft structures. PMR-15 uses the reverse Diels-Alder (RDA) method for crosslinking which provides it with the groundwork for its distinctive thermal stability and a range of 280-300 degree Centigrade use temperature. Regardless of such desirable properties, this material has a number of limitations that compromises its application on a large scale basis. PMR-15 composites has been known to be very vulnerable to micro-cracking at inter and intra-laminar cracking. But the major factor that hinders its demand is PMR-15's carcinogenic constituent, methylene dianilineme (MDA), also a liver toxin. The necessity of providing a safe working environment during its production adds up to the cost of this material. In this study, Molecular Dynamics and Energy Minimization techniques are utilized to simulate a structure of PMR-15 at a given density of 1.324 g/cc and an attempt to recreate the polyimide to reduce the number of experimental testing and hence subdue the health hazards as well as the cost involved in its production. Even though this study does not involve in validating any mechanical properties of the model, it could be used in future for the validation of its properties and further testing for different properties like aging, microcracking, creep etc.

Funktionalisierte Bauteile durch Selektives Maskensintern

Das Selektive Maskensintern ist ein neues, pulverbasiertes Additives Fertigungsverfahren. Das schichtweise aufgebrachte Kunststoffpulver wird hier flächig über einen Infrarotstrahler belichtet und aufgeschmolzen. Das Verfahren bietet ein großes Potential dreidimensionale, wärmeleitfähige Bauteile mit beliebiger Geometrie herzustellen. In diesem Beitrag wird darauf eingegangen, wie kommerziell erhältliches Polyamid 12-Pulver mit thermisch leitfähigen Füllstoffen, wie Aluminiumgrieß und Kupferkugeln, modifiziert und funktionalisierte Bauteile hergestellt werden können. Prozessrelevante Materialeigenschaften werden mittels Differential Scanning Kalorimetrie, Rotationsviskosimetrie und der Wärmeleitfähigkeit der modifizierten Pulver bestimmt. An den gefertigten Bauteilen wird die ausgebildete Morphologie, die mechanischen Eigenschaften als auch die Bauteilwärmeleitfähigkeit untersucht.

Mechanische Eigenschaften und Bruchverhalten maskengesinterter Kunststoffbauteile

Das selektive Maskensintern von Kunststoffen ermöglicht die flächige Belichtung des Bauraums, wodurch sich konstante, von der zu belichtenden Geometrie/Fläche unabhängige Zykluszeiten pro Schicht ergeben. Durch den Einsatz eines, über dem Bauraum platzierten, Infrarotstrahlerfeldes wird eine Modifikation des verarbeiteten Polyamid 12-Pulvers mit einem Absorber, hier Flammruß, notwendig. Bisher konnte gezeigt werden, dass Prototypen sowie wärmeleitfähige Kunststoffbauteile hergestellt werden können. Im Rahmen dieses Beitrags sollen die mechanischen Eigenschaften von SMS-Bauteilen betrachtet werden. Die Beeinflussung der mechanischen Kennwerte, durch variierende Materialeigenschaften sowie unterschiedliche Prozessparameter, werden ebenso wie die Richtungs-, Temperatur- und Belastungsartabhängigkeit, bei konstanten Material- und Prozessparametern, untersucht. Zur Charakterisierung des Bauteilversagens wurden Methoden wie die Lichtmikroskopie und Rasterelektronenmikroskopie eingesetzt und Bruchmechanismen abgeleitet.

Second Annual Kansas State–Nebraska Biochemical Engineering Symposium (April 29, 1972)

This booklet contains abstracts of papers presented at a biochemical engineering symposium conducted at the University of Nebraska-Lincoln on April 29, 1972. This was the second annual symposium on this subject, the first having been held at Kansas State University on June 4, 1971. It is expected that future symposia will alternate between the two campuses. ContentsS.H. Lin, Kansas State University, "Enzyme Reaction in a Tubular Reactor with Laminar Flow" Gregory C. Martin, University of Nebraska, "Estimation of Parameters in Population Models for Schizosaccharomyces pombe from Chemostat Data" Jaiprakash S. Shastry and Prakash N. Mishra, Kansas State University, "Immobilized Enzymes: Analysis of Ultrafiltration Reactors" Mark D. Young, University of Nebraska, "Modelling Unsteady-State Two-Species Data Using Ramkrishna's Staling Model" G.C.Y. Chu, Kansas State University, "Optimization of Step Aeration Waste Treatment Systems Using EVOP" Shinji Goto, University of Nebraska, "Growth of the Blue-Green Alga Microcytis aeruginosa under Defined Conditions" Prakash N. Mishra and Thomas M.C. Kuo, Kansas State University, "Digital Computer Simulation of the Activated Sludge System: Effect of Primary Clarifier on System Performance" Mark D. Young, University of Nebraska, "Aerobic Fermentation of Paunch Liquor"