Testing of subgrid scale (SGS) models for large-eddy simulation (LES) of turbulent channel flow

Sub-grid scale (SGS) models are required in order to model the influence of the unresolved small scales on the resolved scales in large-eddy simulations (LES), the flow at the smallest scales of turbulence. In the following work two SGS models are presented and deeply analyzed in terms of accuracy through several LESs with different spatial resolutions, i.e. grid spacings. The first part of this thesis focuses on the basic theory of turbulence, the governing equations of fluid dynamics and their adaptation to LES. Furthermore, two important SGS models are presented: one is the Dynamic eddy-viscosity model (DEVM), developed by \cite{germano1991dynamic}, while the other is the Explicit Algebraic SGS model (EASSM), by \cite{marstorp2009explicit}. In addition, some details about the implementation of the EASSM in a Pseudo-Spectral Navier-Stokes code \cite{chevalier2007simson} are presented. The performance of the two aforementioned models will be investigated in the following chapters, by means of LES of a channel flow, with friction Reynolds numbers $Re_\tau=590$ up to $Re_\tau=5200$, with relatively coarse resolutions. Data from each simulation will be compared to baseline DNS data. Results have shown that, in contrast to the DEVM, the EASSM has promising potentials for flow predictions at high friction Reynolds numbers: the higher the friction Reynolds number is the better the EASSM will behave and the worse the performances of the DEVM will be. The better performance of the EASSM is contributed to the ability to capture flow anisotropy at the small scales through a correct formulation for the SGS stresses. Moreover, a considerable reduction in the required computational resources can be achieved using the EASSM compared to DEVM. Therefore, the EASSM combines accuracy and computational efficiency, implying that it has a clear potential for industrial CFD usage.

Development and in-vitro characterization of an implantable flow sensing transducer for hydrocephalus

An implantable transducer for monitoring the flow of Cerebrospinal fluid (CSF) for the treatment of hydrocephalus has been developed which is based on measuring the heat dissipation of a local thermal source. The transducer uses passive telemetry at 13.56 MHz for power supply and read out of the measured flow rate. The in vitro performance of the transducer has been characterized using artificial Cerebrospinal Fluid (CSF) with increased protein concentration and artificial CSF with 10\% fresh blood. After fresh blood was added to the artificial CSF a reduction of flow rate has been observed in case that the sensitive surface of the flow sensor is close to the sedimented erythrocytes. An increase of flow rate has been observed in case that the sensitive surface is in contact with the remaining plasma/artificial CSF mix above the sediment which can be explained by an asymmetric flow profile caused by the sedimentation of erythrocythes having increased viscosity compared to artificial CSF. After removal of blood from artificial CSF, no drift could be observed in the transducer measurement which could be associated to a deposition of proteins at the sensitive surface walls of the packaged flow transducer. The flow sensor specification requirement of +-10\% for a flow range between 2 ml/h and 40 ml/h. could be confirmed at test conditions of 37 degrees C.

Performance Characteristics of Airlift Pumps with Vortex Induced by Tangential Fluid Injection

The effect of the swirl component of air injection on the performance of an airlift pump was examined experimentally. An airlift pump is a device that pumps a liquid or slurry using only gas injection. In this study, the liquid used was water and the injected gas was air. The effect of the air swirl was determined by measuring the water discharge from an airlift pump with an air injection nozzle in which the air flow had both axial and tangential components and then repeating the tests with a nozzle with only axial injection. The induced water flow was measured using an orifice meter in the supply pipeline. Tests were run for air pressures ranging from 10 to 30 pounds per square inch, gauge (psig), at flow rates from 5 standard cubic feet per minute (scfm) up the maximum values attainable at the given pressure (usually in the range from 20 to 35 scfm). The nozzle with only axial injection produced a water flow rate that wasequivalent to or better than that induced by the nozzle with swirl. The swirl component of air injection was found to be detrimental to pump performance for all but the smallest air injection flow rate. Optimum efficiency was found for air injection pressures of 10 psig to 15 psig. In addition, the effect of using auxiliary tangential injection of water to create a swirl component in the riser before air injection on the overall capacity (i.e., flow rate) and efficiencyof the pump was examined. Auxiliary tangential water injection was found to have no beneficial effect on the pump capacity or performance in the present system.

The simultaneous role of an alveolus as flow mixer and flow feeder for the deposition of inhaled submicron particles

In an effort to understand the fate of inhaled submicron particles in the small sacs, or alveoli, comprising the gas-exchange region of the lung, we calculated the flow in three-dimensional (3D) rhythmically expanding models of alveolated ducts. Since convection toward the alveolar walls is a precursor to particle deposition, it was the goal of this paper to investigate the streamline maps' dependence upon alveoli location along the acinar tree. On the alveolar midplane, the recirculating flow pattern exhibited closed streamlines with a stagnation saddle point. Off the midplane we found no closed streamlines but nested, funnel-like, spiral, structures (reminiscent of Russian nesting dolls) that were directed towards the expanding walls in inspiration, and away from the contracting walls in expiration. These nested, funnel-like, structures were surrounded by air that flowed into the cavity from the central channel over inspiration and flowed from the cavity to the central channel over expiration. We also found that fluid particle tracks exhibited similar nested funnel-like spiral structures. We conclude that these unique alveolar flow structures may be of importance in enhancing deposition. In addition, due to inertia, the nested, funnel-like, structures change shape and position slightly during a breathing cycle, resulting in flow mixing. Also, each inspiration feeds a fresh supply of particle-laden air from the central channel to the region surrounding the mixing region. Thus, this combination of flow mixer and flow feeder makes each individual alveolus an effective mixing unit, which is likely to play an important role in determining the overall efficiency of convective mixing in the acinus.

High interstitial fluid pressure is associated with low tumour penetration of diagnostic monoclonal antibodies applied for molecular imaging purposes

The human epithelial cell adhesion molecule (EpCAM) is highly expressed in a variety of clinical tumour entities. Although an antibody against EpCAM has successfully been used as an adjuvant therapy in colon cancer, this therapy has never gained wide-spread use. We have therefore investigated the possibilities and limitations for EpCAM as possible molecular imaging target using a panel of preclinical cancer models. Twelve human cancer cell lines representing six tumour entities were tested for their EpCAM expression by qPCR, flow cytometry analysis and immunocytochemistry. In addition, EpCAM expression was analyzed in vivo in xenograft models for tumours derived from these cells. Except for melanoma, all cell lines expressed EpCAM mRNA and protein when grown in vitro. Although they exhibited different mRNA levels, all cell lines showed similar EpCAM protein levels upon detection with monoclonal antibodies. When grown in vivo, the EpCAM expression was unaffected compared to in vitro except for the pancreatic carcinoma cell line 5072 which lost its EpCAM expression in vivo. Intravenously applied radio-labelled anti EpCAM MOC31 antibody was enriched in HT29 primary tumour xenografts indicating that EpCAM binding sites are accessible in vivo. However, bound antibody could only be immunohistochemically detected in the vicinity of perfused blood vessels. Investigation of the fine structure of the HT29 tumour blood vessels showed that they were immature and prone for higher fluid flux into the interstitial space. Consistent with this hypothesis, a higher interstitial fluid pressure of about 12 mbar was measured in the HT29 primary tumour via "wick-in-needle" technique which could explain the limited diffusion of the antibody into the tumour observed by immunohistochemistry.

Biofilm-related infections of cerebrospinal fluid shunts

Cerebrospinal fluid (CSF) shunts carry a high risk of complications. Infections represent a major cause of shunt failure. Diagnosis and therapy of such infections are complicated by the formation of bacterial biofilms attached to shunt surfaces. This study correlated the pathophysiology and clinical course of biofilm infections with microscopical findings on the respective shunts. Surface irregularities, an important risk-factor for shunt colonisation with bacteria, were found to increase over time because of silicone degradation. Scanning electron-microscopy (SEM) documented residual biological material (dead biofilm), which can further promote extant bacterial adhesion, on newly manufactured shunts. Clinical course and SEM both documented bacterial dissemination against CSF flow and the monodirectional valve. In all cases, biofilms grew on both the inner and outer surfaces of the shunts. Microscopy and conventional culture detected all bacterial shunt infections. Analyses of 16S rDNA sequences using conserved primers identified bacteria in only one of three cases, probably because of previous formalin fixation of the samples.

Effect of a lung recruitment maneuver by high-frequency oscillatory ventilation in experimental acute lung injury on organ blood flow in pigs

INTRODUCTION: The objective was to study the effects of a lung recruitment procedure by stepwise increases of mean airway pressure upon organ blood flow and hemodynamics during high-frequency oscillatory ventilation (HFOV) versus pressure-controlled ventilation (PCV) in experimental lung injury. METHODS: Lung damage was induced by repeated lung lavages in seven anesthetized pigs (23-26 kg). In randomized order, HFOV and PCV were performed with a fixed sequence of mean airway pressure increases (20, 25, and 30 mbar every 30 minutes). The transpulmonary pressure, systemic hemodynamics, intracranial pressure, cerebral perfusion pressure, organ blood flow (fluorescent microspheres), arterial and mixed venous blood gases, and calculated pulmonary shunt were determined at each mean airway pressure setting. RESULTS: The transpulmonary pressure increased during lung recruitment (HFOV, from 15 +/- 3 mbar to 22 +/- 2 mbar, P < 0.05; PCV, from 15 +/- 3 mbar to 23 +/- 2 mbar, P < 0.05), and high airway pressures resulted in elevated left ventricular end-diastolic pressure (HFOV, from 3 +/- 1 mmHg to 6 +/- 3 mmHg, P < 0.05; PCV, from 2 +/- 1 mmHg to 7 +/- 3 mmHg, P < 0.05), pulmonary artery occlusion pressure (HFOV, from 12 +/- 2 mmHg to 16 +/- 2 mmHg, P < 0.05; PCV, from 13 +/- 2 mmHg to 15 +/- 2 mmHg, P < 0.05), and intracranial pressure (HFOV, from 14 +/- 2 mmHg to 16 +/- 2 mmHg, P < 0.05; PCV, from 15 +/- 3 mmHg to 17 +/- 2 mmHg, P < 0.05). Simultaneously, the mean arterial pressure (HFOV, from 89 +/- 7 mmHg to 79 +/- 9 mmHg, P < 0.05; PCV, from 91 +/- 8 mmHg to 81 +/- 8 mmHg, P < 0.05), cardiac output (HFOV, from 3.9 +/- 0.4 l/minute to 3.5 +/- 0.3 l/minute, P < 0.05; PCV, from 3.8 +/- 0.6 l/minute to 3.4 +/- 0.3 l/minute, P < 0.05), and stroke volume (HFOV, from 32 +/- 7 ml to 28 +/- 5 ml, P < 0.05; PCV, from 31 +/- 2 ml to 26 +/- 4 ml, P < 0.05) decreased. Blood flows to the heart, brain, kidneys and jejunum were maintained. Oxygenation improved and the pulmonary shunt fraction decreased below 10% (HFOV, P < 0.05; PCV, P < 0.05). We detected no differences between HFOV and PCV at comparable transpulmonary pressures. CONCLUSION: A typical recruitment procedure at the initiation of HFOV improved oxygenation but also decreased systemic hemodynamics at high transpulmonary pressures when no changes of vasoactive drugs and fluid management were performed. Blood flow to the organs was not affected during lung recruitment. These effects were independent of the ventilator mode applied.

Rapid restoration of microcirculatory blood flow with hyperviscous and hyperoncotic solutions lowers the transfusion trigger in resuscitation from hemorrhagic shock

Resuscitation from hemorrhagic shock relies on fluid retransfusion. However, the optimal properties of the fluid have not been established. The aim of the present study was to test the influence of the concentration of hydroxyethyl starch (HES) solution on plasma viscosity and colloid osmotic pressure (COP), systemic and microcirculatory recovery, and oxygen delivery and consumption after resuscitation, which were assessed in the hamster chamber window preparation by intravital microscopy. Awake hamsters were subjected to 50% hemorrhage and were resuscitated with 25% of the estimated blood volume with 5%, 10%, or 20% HES solution. The increase in concentration led to an increase in COP (from 20 to 70 and 194 mmHg) and viscosity (from 1.7 to 3.8 and 14.4 cP). Cardiac index and microcirculatory and metabolic recovery were improved with HES 10% and 20% when compared with 5% HES. Oxygen delivery and consumption in the dorsal skinfold chamber was more than doubled with HES 10% and 20% when compared with HES 5%. This was attributed to the beneficial effect of restored or increased plasma COP and plasma viscosity as obtained with HES 10% and 20%, leading to improved microcirculatory blood flow values early in the resuscitation period. The increase in COP led to an increase in blood volume as shown by a reduction in hematocrit. Mean arterial pressure was significantly improved in animals receiving 10% and 20% solutions. In conclusion, the present results show that the increase in the concentration of HES, leading to hyperoncotic and hyperviscous solutions, is beneficial for resuscitation from hemorrhagic shock because normalization of COP and viscosity led to a rapid recovery of microcirculatory parameters.

Endothelin inhibition improves cerebral blood flow and is neuroprotective in pneumococcal meningitis

By using an infant rat model of pneumococcal meningitis, we determined whether endothelins contribute to neuronal damage in this disease. Cerebrospinal fluid analysis demonstrated a significant increase of endothelin-1 in infected animals compared with uninfected controls. Histopathological examination 24 hours after infection showed brain damage in animals treated with ceftriaxone alone (median, 9.2% of cortex; range, 0-49.1%). In infected animals treated intraperitoneally with the endothelin antagonist bosentan (30 mg/kg, every 12 hours) also, injury was reduced to 0.5% (range, 0-8.6%) of cortex. Cerebral blood flow was reduced in infected animals (6.5 +/- 4.0 ml/min/100 g of brain vs 14.9 +/- 9.1 ml/min/100 g in controls. Treatment with bosentan restored cerebral blood flow to levels similar to controls (12.8 +/- 5.3 ml/min/100 g). Improved blood flow was not mediated by nitric oxide production, because bosentan had no effect on cerebrospinal fluid or plasma nitrite/nitrate concentrations at 6, 12, or 18 hours. These data indicate that endothelins contribute to neuronal injury in this model of pneumococcal meningitis by causing cerebral ischemia.

Rapid diagnosis of experimental meningitis by bacterial heat production in cerebrospinal fluid

BACKGROUND: Calorimetry is a nonspecific technique which allows direct measurement of heat generated by biological processes in the living cell. We evaluated the potential of calorimetry for rapid detection of bacterial growth in cerebrospinal fluid (CSF) in a rat model of bacterial meningitis. METHODS: Infant rats were infected on postnatal day 11 by direct intracisternal injection with either Streptococcus pneumoniae, Neisseria meningitidis or Listeria monocytogenes. Control animals were injected with sterile saline or heat-inactivated S. pneumoniae. CSF was obtained at 18 hours after infection for quantitative cultures and heat flow measurement. For calorimetry, 10 microl and 1 microl CSF were inoculated in calorimetry ampoules containing 3 ml trypticase soy broth (TSB). RESULTS: The mean bacterial titer (+/- SD) in CSF was 1.5 +/- 0.6 x 108 for S. pneumoniae, 1.3 +/- 0.3 x 106 for N. meningitidis and 3.5 +/- 2.2 x 104 for L. monocytogenes. Calorimetric detection time was defined as the time until heat flow signal exceeded 10 microW. Heat signal was detected in 10-microl CSF samples from all infected animals with a mean (+/- SD) detection time of 1.5 +/- 0.2 hours for S. pneumoniae, 3.9 +/- 0.7 hours for N. meningitidis and 9.1 +/- 0.5 hours for L. monocytogenes. CSF samples from non-infected animals generated no increasing heat flow (<10 microW). The total heat was the highest in S. pneumoniae ranging from 6.7 to 7.5 Joules, followed by L. monocytogenes (5.6 to 6.1 Joules) and N. meningitidis (3.5 to 4.4 Joules). The lowest detectable bacterial titer by calorimetry was 2 cfu for S. pneumoniae, 4 cfu for N. meningitidis and 7 cfu for L. monocytogenes. CONCLUSION: By means of calorimetry, detection times of <4 hours for S. pneumoniae and N. meningitidis and <10 hours for Listeria monocytogenes using as little as 10 microl CSF were achieved. Calorimetry is a new diagnostic method allowing rapid and accurate diagnosis of bacterial meningitis from a small volume of CSF.

Perioperative fluid management: comparison of high, medium and low fluid volume on tissue oxygen pressure in the small bowel and colon

BACKGROUND AND OBJECTIVE: Insufficient blood flow and oxygenation in the intestinal tract is associated with increased incidence of postoperative complications after bowel surgery. High fluid volume administration may prevent occult regional hypoperfusion and intestinal tissue hypoxia. We tested the hypothesis that high intraoperative fluid volume administration increases intestinal wall tissue oxygen pressure during laparotomy. METHODS: In all, 27 pigs were anaesthetized, ventilated and randomly assigned to one of the three treatment groups (n = 9 in each) receiving low (3 mL kg-1 h-1), medium (7 mL kg-1 h-1) or high (20 mL kg-1 h-1) fluid volume treatment with lactated Ringer's solution. All animals received 30% and 100% inspired oxygen in random order. Cardiac index was measured with thermodilution and tissue oxygen pressure with a micro-oximetry system in the jejunum and colon wall and subcutaneous tissue. RESULTS: Groups receiving low and medium fluid volume treatment had similar systemic haemodynamics. The high fluid volume group had significantly higher mean arterial pressure, cardiac index and subcutaneous tissue oxygenation. Tissue oxygen pressures in the jejunum and colon were comparable in all three groups. CONCLUSIONS: The three different fluid volume regimens tested did not affect tissue oxygen pressure in the jejunum and colon, suggesting efficient autoregulation of intestinal blood flow in healthy subjects undergoing uncomplicated abdominal surgery.

Modeling of transfer kinetics at the serum-cerebrospinal fluid barrier in rabbits with experimental meningitis: application to grepafloxacin

The goals of the present study were to model the population kinetics of in vivo influx and efflux processes of grepafloxacin at the serum-cerebrospinal fluid (CSF) barrier and to propose a simulation-based approach to optimize the design of dose-finding trials in the meningitis rabbit model. Twenty-nine rabbits with pneumococcal meningitis receiving grepafloxacin at 15 mg/kg of body weight (intravenous administration at 0 h), 30 mg/kg (at 0 h), or 50 mg/kg twice (at 0 and 4 h) were studied. A three-compartment population pharmacokinetic model was fit to the data with the program NONMEM (Nonlinear Mixed Effects Modeling). Passive diffusion clearance (CL(diff)) and active efflux clearance (CL(active)) are transfer kinetic modeling parameters. Influx clearance is assumed to be equal to CL(diff), and efflux clearance is the sum of CL(diff), CL(active), and bulk flow clearance (CL(bulk)). The average influx clearance for the population was 0.0055 ml/min (interindividual variability, 17%). Passive diffusion clearance was greater in rabbits receiving grepafloxacin at 15 mg/kg than in those treated with higher doses (0.0088 versus 0.0034 ml/min). Assuming a CL(bulk) of 0.01 ml/min, CL(active) was estimated to be 0.017 ml/min (11%), and clearance by total efflux was estimated to be 0.032 ml/min. The population kinetic model allows not only to quantify in vivo efflux and influx mechanisms at the serum-CSF barrier but also to analyze the effects of different dose regimens on transfer kinetic parameters in the rabbit meningitis model. The modeling-based approach also provides a tool for the simulation and prediction of various outcomes in which researchers might be interested, which is of great potential in designing dose-finding trials.

Crystalloids versus colloids for goal-directed fluid therapy in major surgery

INTRODUCTION: Perioperative hypovolemia arises frequently and contributes to intestinal hypoperfusion and subsequent postoperative complications. Goal-directed fluid therapy might reduce these complications. The aim of this study was to compare the effects of goal-directed administration of crystalloids and colloids on the distribution of systemic, hepatosplanchnic, and microcirculatory (small intestine) blood flow after major abdominal surgery in a clinically relevant pig model. METHODS: Twenty-seven pigs were anesthetized and mechanically ventilated and underwent open laparotomy. They were randomly assigned to one of three treatment groups: the restricted Ringer lactate (R-RL) group (n = 9) received 3 mL/kg per hour of RL, the goal-directed RL (GD-RL) group (n = 9) received 3 mL/kg per hour of RL and intermittent boluses of 250 mL of RL, and the goal-directed colloid (GD-C) group (n = 9) received 3 mL/kg per hour of RL and boluses of 250 mL of 6% hydroxyethyl starch (130/0.4). The latter two groups received a bolus infusion when mixed venous oxygen saturation was below 60% ('lockout' time of 30 minutes). Regional blood flow was measured in the superior mesenteric artery and the celiac trunk. In the small bowel, microcirculatory blood flow was measured using laser Doppler flowmetry. Intestinal tissue oxygen tension was measured with intramural Clark-type electrodes. RESULTS: After 4 hours of treatment, arterial blood pressure, cardiac output, mesenteric artery flow, and mixed oxygen saturation were significantly higher in the GD-C and GD-RL groups than in the R-RL group. Microcirculatory flow in the intestinal mucosa increased by 50% in the GD-C group but remained unchanged in the other two groups. Likewise, tissue oxygen tension in the intestine increased by 30% in the GD-C group but remained unchanged in the GD-RL group and decreased by 18% in the R-RL group. Mesenteric venous glucose concentrations were higher and lactate levels were lower in the GD-C group compared with the two crystalloid groups. CONCLUSIONS: Goal-directed colloid administration markedly increased microcirculatory blood flow in the small intestine and intestinal tissue oxygen tension after abdominal surgery. In contrast, goal-directed crystalloid and restricted crystalloid administrations had no such effects. Additionally, mesenteric venous glucose and lactate concentrations suggest that intestinal cellular substrate levels were higher in the colloid-treated than in the crystalloid-treated animals. These results support the notion that perioperative goal-directed therapy with colloids might be beneficial during major abdominal surgery.

Characterizing the effects of turbulence on sediment oxygen and sediment nitrate demand using flow-through reactors

Onondaga Lake has received the municipal effluent and industrial waste from the city of Syracuse for more than a century. Historically, 75 metric tons of mercury were discharged to the lake by chlor-alkali facilities. These legacy deposits of mercury now exist primarily in the lake sediments. Under anoxic conditions, methylmercury is produced in the sediments and can be released to the overlying water. Natural sedimentation processes are continuously burying the mercury deeper into the sediments. Eventually, the mercury will be buried to a depth where it no longer has an impact on the overlying water. In the interim, electron acceptor amendment systems can be installed to retard these chemical releases while the lake naturally recovers. Electron acceptor amendment systems are designed to meet the sediment oxygen demand in the sediment and maintain manageable hypolimnion oxygen concentrations. Historically, designs of these systems have been under designed resulting in failure. This stems from a mischaracterization of the sediment oxygen demand. Turbulence at the sediment water interface has been shown to impact sediment oxygen demand. The turbulence introduced by the electron amendment system can thus increase the sediment oxygen demand, resulting in system failure if turbulence is not factored into the design. Sediment cores were gathered and operated to steady state under several well characterized turbulence conditions. The relationship between sediment oxygen/nitrate demand and turbulence was then quantified and plotted. A maximum demand was exhibited at or above a fluid velocity of 2.0 mm•s-1. Below this velocity, demand decreased rapidly with fluid velocity as zero velocity was approached. Similar relationships were displayed by both oxygen and nitrate cores.

Analyzing the effects of vertical acceleration and deceleration in multi-phase, debris-flow sediment-column experiments / by Christine Marie Williams.

Water-saturated debris flows are among some of the most destructive mass movements. Their complex nature presents a challenge for quantitative description and modeling. In order to improve understanding of the dynamics of these flows, it is important to seek a simplified dynamic system underlying their behavior. Models currently in use to describe the motion of debris flows employ depth-averaged equations of motion, typically assuming negligible effects from vertical acceleration. However, in many cases debris flows experience significant vertical acceleration as they move across irregular surfaces, and it has been proposed that friction associated with vertical forces and liquefaction merit inclusion in any comprehensive mechanical model. The intent of this work is to determine the effect of vertical acceleration through a series of laboratory experiments designed to simulate debris flows, testing a recent model for debris flows experimentally. In the experiments, a mass of water-saturated sediment is released suddenly from a holding container, and parameters including rate of collapse, pore-fluid pressure, and bed load are monitored. Experiments are simplified to axial geometry so that variables act solely in the vertical dimension. Steady state equations to infer motion of the moving sediment mass are not sufficient to model accurately the independent solid and fluid constituents in these experiments. The model developed in this work more accurately predicts the bed-normal stress of a saturated sediment mass in motion and illustrates the importance of acceleration and deceleration.