Physical Insights, Steady Aerodynamic Effects, and a Design Tool for Low-Pressure Turbine Flutter

The successful, efficient, and safe turbine design requires a thorough understanding of the underlying physical phenomena. This research investigates the physical understanding and parameters highly correlated to flutter, an aeroelastic instability prevalent among low pressure turbine (LPT) blades in both aircraft engines and power turbines. The modern way of determining whether a certain cascade of LPT blades is susceptible to flutter is through time-expensive computational fluid dynamics (CFD) codes. These codes converge to solution satisfying the Eulerian conservation equations subject to the boundary conditions of a nodal domain consisting fluid and solid wall particles. Most detailed CFD codes are accompanied by cryptic turbulence models, meticulous grid constructions, and elegant boundary condition enforcements all with one goal in mind: determine the sign (and therefore stability) of the aerodynamic damping. The main question being asked by the aeroelastician, ``is it positive or negative?'' This type of thought-process eventually gives rise to a black-box effect, leaving physical understanding behind. Therefore, the first part of this research aims to understand and reveal the physics behind LPT flutter in addition to several related topics including acoustic resonance effects. A percentage of this initial numerical investigation is completed using an influence coefficient approach to study the variation the work-per-cycle contributions of neighboring cascade blades to a reference airfoil. The second part of this research introduces new discoveries regarding the relationship between steady aerodynamic loading and negative aerodynamic damping. Using validated CFD codes as computational wind tunnels, a multitude of low-pressure turbine flutter parameters, such as reduced frequency, mode shape, and interblade phase angle, will be scrutinized across various airfoil geometries and steady operating conditions to reach new design guidelines regarding the influence of steady aerodynamic loading and LPT flutter. Many pressing topics influencing LPT flutter including shocks, their nonlinearity, and three-dimensionality are also addressed along the way. The work is concluded by introducing a useful preliminary design tool that can estimate within seconds the entire aerodynamic damping versus nodal diameter curve for a given three-dimensional cascade.

Fluid resuscitation practice patterns in intensive care units of the USA: a cross-sectional survey of critical care physicians.

BACKGROUND: Fluid resuscitation is a cornerstone of intensive care treatment, yet there is a lack of agreement on how various types of fluids should be used in critically ill patients with different disease states. Therefore, our goal was to investigate the practice patterns of fluid utilization for resuscitation of adult patients in intensive care units (ICUs) within the USA. METHODS: We conducted a cross-sectional online survey of 502 physicians practicing in medical and surgical ICUs. Survey questions were designed to assess clinical decision-making processes for 3 types of patients who need volume expansion: (1) not bleeding and not septic, (2) bleeding but not septic, (3) requiring resuscitation for sepsis. First-choice fluid used in fluid boluses for these 3 patient types was requested from the respondents. Descriptive statistics were performed using a Kruskal-Wallis test to evaluate differences among the physician groups. Follow-up tests, including t tests, were conducted to evaluate differences between ICU types, hospital settings, and bolus volume. RESULTS: Fluid resuscitation varied with respect to preferences for the factors to determine volume status and preferences for fluid types. The 3 most frequently preferred volume indicators were blood pressure, urine output, and central venous pressure. Regardless of the patient type, the most preferred fluid type was crystalloid, followed by 5 % albumin and then 6 % hydroxyethyl starches (HES) 450/0.70 and 6 % HES 600/0.75. Surprisingly, up to 10 % of physicians still chose HES as the first choice of fluid for resuscitation in sepsis. The clinical specialty and the practice setting of the treating physicians also influenced fluid choices. CONCLUSIONS: Practice patterns of fluid resuscitation varied in the USA, depending on patient characteristics, clinical specialties, and practice settings of the treating physicians.

Perioperative fluid therapy: a statement from the international Fluid Optimization Group.

BACKGROUND: Perioperative fluid therapy remains a highly debated topic. Its purpose is to maintain or restore effective circulating blood volume during the immediate perioperative period. Maintaining effective circulating blood volume and pressure are key components of assuring adequate organ perfusion while avoiding the risks associated with either organ hypo- or hyperperfusion. Relative to perioperative fluid therapy, three inescapable conclusions exist: overhydration is bad, underhydration is bad, and what we assume about the fluid status of our patients may be incorrect. There is wide variability of practice, both between individuals and institutions. The aims of this paper are to clearly define the risks and benefits of fluid choices within the perioperative space, to describe current evidence-based methodologies for their administration, and ultimately to reduce the variability with which perioperative fluids are administered. METHODS: Based on the abovementioned acknowledgements, a group of 72 researchers, well known within the field of fluid resuscitation, were invited, via email, to attend a meeting that was held in Chicago in 2011 to discuss perioperative fluid therapy. From the 72 invitees, 14 researchers representing 7 countries attended, and thus, the international Fluid Optimization Group (FOG) came into existence. These researches, working collaboratively, have reviewed the data from 162 different fluid resuscitation papers including both operative and intensive care unit populations. This manuscript is the result of 3 years of evidence-based, discussions, analysis, and synthesis of the currently known risks and benefits of individual fluids and the best methods for administering them. RESULTS: The results of this review paper provide an overview of the components of an effective perioperative fluid administration plan and address both the physiologic principles and outcomes of fluid administration. CONCLUSIONS: We recommend that both perioperative fluid choice and therapy be individualized. Patients should receive fluid therapy guided by predefined physiologic targets. Specifically, fluids should be administered when patients require augmentation of their perfusion and are also volume responsive. This paper provides a general approach to fluid therapy and practical recommendations.

A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.

Seawater carbonate chemistry and hydrogen ions and carbonate chemistry in calcifying fluid of coral Astrangia poculata during experiments, 2011

A generalized physicochemical model of the response of marine organisms' calcifying fluids to CO2-induced ocean acidification is proposed. The model is based upon the hypothesis that some marine calcifiers induce calcification by elevating pH, and thus Omega aragonite, of their calcifying fluid by removing protons (H+). The model is explored through two end-member scenarios: one in which a fixed number of H+ is removed from their calcifying fluid, regardless of atmospheric pCO2, and another in which a fixed external-internal proton ratio ([H+]E/[H+]I) is maintained. The model is able to generate the full range of calcification response patterns observed in prior ocean acidification experiments and is consistent with the assertion that organisms' calcification response to ocean acidification is more negative for marine calcifiers that exert weaker control over their calcifying fluid pH. The model is empirically evaluated for the temperate scleractinian coral Astrangia poculata with in situ pH microelectrode measurements of the coral's calcifying fluid under control and acidified conditions. These measurements reveal that (1) the pH of the coral's calcifying fluid is substantially elevated relative to its external seawater under both control and acidified conditions, (2) the coral's [H+]E/[H+]I remains constant under control and acidified conditions, and (3) the coral removes fewer H+ from its calcifying fluid under acidified conditions than under control conditions. Thus, the carbonate system dynamics of A. poculata's calcifying fluid appear to be most consistent with the fixed [H+]E/[H+]I end-member scenario. Similar microelectrode experiments performed on additional taxa are required to assess the model's general applicability.

Stochastic Modeling of Reversible Biochemical Reaction-Diffusion Systems and High-Resolution Shock-Capturing Methods for Fluid Interfaces

Thesis (Ph.D.)--University of Washington, 2016-08

Analysis of heat transfer coefficient in silica nanofluids with water as base fluid under transitional flow

The effect of Reynolds number variation in a vertical double pipe counterflow heat exchanger due to the changes in viscosity can cause the change in flow regime, for instance, when heats up and cools down, it can convert from turbulent to laminar or inversely, that can have significant effect on heat transfer coefficient and pressure drop. Mainly, the range of transition phase has been studied in this study with the investigation of silica nanofluid dispersed in water in three different concentrations. The results have been compared with distilled water sample and showed a remarkable raise in heat transfer coefficient while pressure drop has been increased respectively, as well. Although pumping power has to go up at the same time and it is a drawback, heat transfer efficiency grows for diluted samples. On the other hand, for the most concentrated sample, effect of pressure drop dominates which leads to decline in the overall efficiency.

Tris(2-ethylhexyl) trimellitate (TOTM) as a potential industrial reference fluid for viscosity at high temperatures and high pressures: New viscosity, density and surface tension measurements

Tris(2-ethylhexyl) trimellitate (TOTM) was recently suggested as a reference fluid for industrial use associated with high viscosity at elevated temperature and pressure. Viscosity and density data have already been published on one sample covering the temperature range (303-373) K and at pressures up to about 65 MPa. The viscosity covered a range from about (9 to 460) mPa s. In the present article we study several other characteristics of TOTM that must be available if it were to be adopted as a standard. First, we present values for the viscosity and density obtained with a different sample of TOTM to examine the important feature of consistency among different samples. Vibrating-wire viscosity measurements were performed at pressures from (5 to 100) MPa, along 6 isotherms between (303 and 373) K. Density measurements were carried out from (293 to 373) K up to 68 MPa, along 4 isotherms, using an Anton Paar DMA HP vibrating U-tube densimeter. Secondly, we report a study of the effect of water contamination on the viscosity of TOTM, performed using an Ubbelhode viscometer under atmospheric pressure. Finally, in order to support the use of TOTM as a reference liquid for the calibration of capillary viscometers, values of its surface tension, obtained by the pendant drop method, are provided. (C) 2016 Elsevier B.V. All rights reserved.

Computational modelling of non-equilibrium condensing steam flows in low-pressure steam turbines

The steam turbines play a significant role in global power generation. Especially, research on low pressure (LP) steam turbine stages is of special importance for steam turbine man- ufactures, vendors, power plant owners and the scientific community due to their lower efficiency than the high pressure steam turbine stages. Because of condensation, the last stages of LP turbine experience irreversible thermodynamic losses, aerodynamic losses and erosion in turbine blades. Additionally, an LP steam turbine requires maintenance due to moisture generation, and therefore, it is also affecting on the turbine reliability. Therefore, the design of energy efficient LP steam turbines requires a comprehensive analysis of condensation phenomena and corresponding losses occurring in the steam tur- bine either by experiments or with numerical simulations. The aim of the present work is to apply computational fluid dynamics (CFD) to enhance the existing knowledge and understanding of condensing steam flows and loss mechanisms that occur due to the irre- versible heat and mass transfer during the condensation process in an LP steam turbine. Throughout this work, two commercial CFD codes were used to model non-equilibrium condensing steam flows. The Eulerian-Eulerian approach was utilised in which the mix- ture of vapour and liquid phases was solved by Reynolds-averaged Navier-Stokes equa- tions. The nucleation process was modelled with the classical nucleation theory, and two different droplet growth models were used to predict the droplet growth rate. The flow turbulence was solved by employing the standard k-ε and the shear stress transport k-ω turbulence models. Further, both models were modified and implemented in the CFD codes. The thermodynamic properties of vapour and liquid phases were evaluated with real gas models. In this thesis, various topics, namely the influence of real gas properties, turbulence mod- elling, unsteadiness and the blade trailing edge shape on wet-steam flows, are studied with different convergent-divergent nozzles, turbine stator cascade and 3D turbine stator-rotor stage. The simulated results of this study were evaluated and discussed together with the available experimental data in the literature. The grid independence study revealed that an adequate grid size is required to capture correct trends of condensation phenomena in LP turbine flows. The study shows that accurate real gas properties are important for the precise modelling of non-equilibrium condensing steam flows. The turbulence modelling revealed that the flow expansion and subsequently the rate of formation of liquid droplet nuclei and its growth process were affected by the turbulence modelling. The losses were rather sensitive to turbulence modelling as well. Based on the presented results, it could be observed that the correct computational prediction of wet-steam flows in the LP turbine requires the turbulence to be modelled accurately. The trailing edge shape of the LP turbine blades influenced the liquid droplet formulation, distribution and sizes, and loss generation. The study shows that the semicircular trailing edge shape predicted the smallest droplet sizes. The square trailing edge shape estimated greater losses. The analysis of steady and unsteady calculations of wet-steam flow exhibited that in unsteady simulations, the interaction of wakes in the rotor blade row affected the flow field. The flow unsteadiness influenced the nucleation and droplet growth processes due to the fluctuation in the Wilson point.

A randomised controlled trial comparing weight adjusted dose versus fixed dose prophylactic phenylephrine infusion on maintaining systolic blood pressure during caesarean section under spinal anaesthesia.

Background: Spinal anaesthesia is the standard of care for elective caesarean delivery. It has advantages over general anaesthesia. However the sympathetic blockade induced by spinal anaesthesia results in an 80 percent incidence of hypotension without prophylactic management. Current evidence supports co-loading with intravenous fluids in conjunction with the use of vasopressors as the most effective way to prevent and treat the hypotension. Phenylephrine is the accepted vasopressor of choice in the parturient. A prophylactic phenylephrine infusion combined with a fluid co-load is proven to be an effective and safe method of maintaining maternal hemodynamic stability. While most published studies have assessed the effectiveness of a prophylactic phenylephrine fixed dose infusion, few studies have assessed the effect of a prophylactic phenylephrine weight adjusted dose infusion on maintaining maternal hemodynamic stability following spinal anesthesia for a cesarean delivery. Objective: To compare the incidence of hypotension between women undergoing elective caesarean section under spinal anaesthesia, receiving prophylactic phenylephrine infusion at a fixed dose of 37.5 micrograms per minute versus a weight adjusted dose of 0.5 micrograms per kilogram per minute. Methods: One hundred and eight patients scheduled for non-urgent caesarean section under spinal anaesthesia were randomized into 2 groups; control group and intervention group using a computer generated table of numbers. Control group; Received prophylactic phenylephrine fixed dose infusion at 37.5 micrograms per minute. Intervention group; Received prophylactic phenylephrine weight adjusted dose infusion at 0.5 micrograms per kilogram per minute Results: The two groups had similar baseline characteristics in terms of ; Age, sex, weight and height. There was a 35.2% incidence of hypotension in the fixed dose group and an 18.6% incidence of hypotension in the weight adjusted dose group. This difference was found to be of borderline statistical significance p-value 0.05, and the difference in the incidence rates between the two groups was found to be statistically significant p= 0.03. The difference in the incidence of reactive hypertension and bradycardia between the two groups was not statistically significant: p-value of 0.19 for reactive hypertension and p-value of 0.42 for the incidence of bradycardia. There was also no statistically significant difference in the use of phenylephrine boluses, use of atropine, intravenous fluid used and the number of times the infusion was stopped. Conclusion: Among this population, the incidence of hypotension was significantly less in the weight adjusted dose group than in the fixed dose group. There was no difference in the number of physician interventions required to keep the blood pressure within 20% of baseline, and no difference in the proportion of reactive hypertension or bradycardia between the two groups. Administering prophylactic phenylephrine infusion at a weight adjusted dose of 0.5 micrograms per kilogram per minute results in a lower incidence of hypotension compared to its administration at a fixed dose of 37.5 micrograms per minute.

A Numerical study of fluid flow in peristaltic pumps and polymeric elastic pipes

This thesis work deals with a mathematical description of flow in polymeric pipe and in a specific peristaltic pump. This study involves fluid-structure interaction analysis in presence of complex-turbulent flows treated in an arbitrary Lagrangian-Eulerian (ALE) framework. The flow simulations are performed in COMSOL 4.4, as 2D axial symmetric model, and ABAQUS 6.14.1, as 3D model with symmetric boundary conditions. In COMSOL, the fluid and structure problems are coupled by monolithic algorithm, while ABAQUS code links ABAQUS CFD and ABAQUS Standard solvers with single block-iterative partitioned algorithm. For the turbulent features of the flow, the fluid model in both codes is described by RNG k-ϵ. The structural model is described, on the basis of the pipe material, by Elastic models or Hyperelastic Neo-Hookean models with Rayleigh damping properties. In order to describe the pulsatile fluid flow after the pumping process, the available data are often defective for the fluid problem. Engineering measurements are normally able to provide average pressure or velocity at a cross-section. This problem has been analyzed by McDonald's and Womersley's work for average pressure at fixed cross section by Fourier analysis since '50, while nowadays sophisticated techniques including Finite Elements and Finite Volumes exist to study the flow. Finally, we set up peristaltic pipe simulations in ABAQUS code, by using the same model previously tested for the fl uid and the structure.

A two-fluid model for tissue growth within a dynamic flow environment

We study the growth of a tissue construct in a perfusion bioreactor, focussing on its response to the mechanical environment. The bioreactor system is modelled as a two-dimensional channel containing a tissue construct through which a flow of culture medium is driven. We employ a multiphase formulation of the type presented by G. Lemon, J. King, H. Byrne, O. Jensen and K. Shakesheff in their study (Multiphase modelling of tissue growth using the theory of mixtures. J. Math. Biol. 52(2), 2006, 571–594) restricted to two interacting fluid phases, representing a cell population (and attendant extracellular matrix) and a culture medium, and employ the simplifying limit of large interphase viscous drag after S. Franks in her study (Mathematical Modelling of Tumour Growth and Stability. Ph.D. Thesis, University of Nottingham, UK, 2002) and S. Franks and J. King in their study Interactions between a uniformly proliferating tumour and its surrounding: Uniform material properties. Math. Med. Biol. 20, 2003, 47–89). The novel aspects of this study are: (i) the investigation of the effect of an imposed flow on the growth of the tissue construct, and (ii) the inclusion of a chanotransduction mechanism regulating the response of the cells to the local mechanical environment. Specifically, we consider the response of the cells to their local density and the culture medium pressure. As such, this study forms the first step towards a general multiphase formulation that incorporates the effect of mechanotransduction on the growth and morphology of a tissue construct. The model is analysed using analytic and numerical techniques, the results of which illustrate the potential use of the model to predict the dominant regulatory stimuli in a cell population.

Specific shear-dependent viscoelastic third-grade fluid model

A specific modified constitutive equation for a third-grade fluid is proposed so that the model be suitable for applications where shear-thinning or shear-thickening may occur. For that, we use the Cosserat theory approach reducing the exact three-dimensional equations to a system depending only on time and on a single spatial variable. This one-dimensional system is obtained by integrating the linear momentum equation over the cross-section of the tube, taking a velocity field approximation provided by the Cosserat theory. From this reduced system, we obtain the unsteady equations for the wall shear stress and mean pressure gradient depending on the volume flow rate, Womersley number, viscoelastic coefficient and flow index over a finite section of the tube geometry with constant circular cross-section.

Vapor pressure and liquid density of fluorinated alcohols:Experimental, simulation and GC-SAFT-VR predictions

The vapor pressure of four liquid 1H,1H-perfluoroalcohols (CF3(CF2)n(CH2)OH, n ¼ 1, 2, 3, 4), often called odd-fluorotelomer alcohols, was measured as a function of temperature between 278 K and 328 K. Liquid densities were also measured for a temperature range between 278 K and 353 K. Molar enthalpies of vaporization were calculated from the experimental data. The results are compared with data from the literature for other perfluoroalcohols as well as with the equivalent hydrogenated alcohols. The results were modeled and interpreted using molecular dynamics simulations and the GC-SAFT-VR equation of state.

Wind loading predictions using computational fluid dynamics

Using Computational Wind Engineering, CWE, for solving wind-related problems is still a challenging task today, mainly due to the high computational cost required to obtain trustworthy simulations. In particular, the Large Eddy Simulation, LES, has been widely used for evaluating wind loads on buildings. The present thesis assesses the capability of LES as a design tool for wind loading predictions through three cases. The first case is using LES for simulating the wind field around a ground-mounted rectangular prism in Atmospheric Boundary Layer (ABL) flow. The numerical results are validated with experimental results for seven wind attack angles, giving a global understanding of the model performance. The case with the worst model behaviour is investigated, including the spatial distribution of the pressure coefficients and their discrepancies with respect to experimental results. The effects of some numerical parameters are investigated for this case to understand their effectiveness in modifying the obtained numerical results. The second case is using LES for investigating the wind effects on a real high-rise building, aiming at validating the performance of LES as a design tool in practical applications. The numerical results are validated with the experimental results in terms of the distribution of the pressure statistics and the global forces. The mesh sensitivity and the computational cost are discussed. The third case is using LES for studying the wind effects on the new large-span roof over the Bologna stadium. The dynamic responses are analyzed and design envelopes for the structure are obtained. Although it is a numerical simulation before the traditional wind tunnel tests, i.e. the validation of the numerical results are not performed, the preliminary evaluations can effectively inform later investigations and provide the final design processes with deeper confidence regarding the absence of potentially unexpected behaviours.