Barriers to front propagation in ordered and disordered vortex flows

We present experiments on reactive front propagation in a two-dimensional (2D) vortex chain flow (both time-independent and time-periodic) and a 2D spatially disordered (time-independent) vortex-dominated flow. The flows are generated using magnetohydrodynamic forcing techniques, and the fronts are produced using the excitable, ferroin-catalyzed Belousov-Zhabotinsky chemical reaction. In both of these flows, front propagation is dominated by the presence of burning invariant manifolds (BIMs) that act as barriers, similar to invariant manifolds that dominate the transport of passive impurities. Convergence of the fronts onto these BIMs is shown experimentally for all of the flows studied. The BIMs are also shown to collapse onto the invariant manifolds for passive transport in the limit of large flow velocities. For the disordered flow, the measured BIMs are compared to those predicted using a measured velocity field and a three-dimensional set of ordinary differential equations that describe the dynamics of front propagation in advection-reaction-diffusion systems.

Assessing the Validity of Brain Glucose Concentration Measurement Using Microdialysis

Brain functions, such as learning, orchestrating locomotion, memory recall, and processing information, all require glucose as a source of energy. During these functions, the glucose concentration decreases as the glucose is being consumed by brain cells. By measuring this drop in concentration, it is possible to determine which parts of the brain are used during specific functions and consequently, how much energy the brain requires to complete the function. One way to measure in vivo brain glucose levels is with a microdialysis probe. The drawback of this analytical procedure, as with many steadystate fluid flow systems, is that the probe fluid will not reach equilibrium with the brain fluid. Therefore, brain concentration is inferred by taking samples at multiple inlet glucose concentrations and finding a point of convergence. The goal of this thesis is to create a three-dimensional, time-dependent, finite element representation of the brainprobe system in COMSOL 4.2 that describes the diffusion and convection of glucose. Once validated with experimental results, this model can then be used to test parameters that experiments cannot access. When simulations were run using published values for physical constants (i.e. diffusivities, density and viscosity), the resulting glucose model concentrations were within the error of the experimental data. This verifies that the model is an accurate representation of the physical system. In addition to accurately describing the experimental brain-probe system, the model I created is able to show the validity of zero-net-flux for a given experiment. A useful discovery is that the slope of the zero-net-flux line is dependent on perfusate flow rate and diffusion coefficients, but it is independent of brain glucose concentrations. The model was simplified with the realization that the perfusate is at thermal equilibrium with the brain throughout the active region of the probe. This allowed for the assumption that all model parameters are temperature independent. The time to steady-state for the probe is approximately one minute. However, the signal degrades in the exit tubing due to Taylor dispersion, on the order of two minutes for two meters of tubing. Given an analytical instrument requiring a five μL aliquot, the smallest brain process measurable for this system is 13 minutes.

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.

Do women have impaired regional systolic function in hypertensive heart disease? A 3-dimensional echocardiography study

AIMS: In pressure overload left ventricular (LV) hypertrophy, gender-related differences in global LV systolic function have been previously reported. The goal of this study was to determine regional systolic function of the left ventricle in male and female patients with hypertensive heart disease. METHODS AND RESULTS: Regional LV function was analyzed from multiplane transesophageal echocardiography with three-dimensional (3D) reconstruction of the left ventricle. In 24 patients (13 males and 11 females), four parallel (2 basal and 2 apical) equidistant short axis cross-sections from base to apex were obtained from the reconstructed LV. In each short axis 24 wall-thickness measurements were carried out at 15 degrees intervals at end-diastole and end-systole. Thus, a total of 192 measurements were obtained in each patient. Wall thickening was calculated as difference of end-diastolic and end-systolic wall thickness, and fractional thickening as thickening divided by end-diastolic thickness. Fractional thickening and wall stress were inversely related to end-diastolic wall thickness in both, males and females. Females showed less LV systolic function when compared to males (p<0.001). However, when corrected for wall stress, which was higher in females, there was no gender difference in systolic function. CONCLUSION: There are regional differences in LV systolic function in females and males which are directly related to differences in wall stress. Thus, gender-related differences in LV regional function are load-dependent and not due to structural differences.

Development of the premature infant nose throat-model (PrINT-Model): an upper airway replica of a premature neonate for the study of aerosol delivery

Clinical efficacy of aerosol therapy in premature newborns depends on the efficiency of delivery of aerosolized drug to the bronchial tree. To study the influence of various anatomical, physical, and physiological factors on aerosol delivery in preterm newborns, it is crucial to have appropriate in vitro models, which are currently not available. We therefore constructed the premature infant nose throat-model (PrINT-Model), an upper airway model corresponding to a premature infant of 32-wk gestational age by three-dimensional (3D) reconstruction of a three-planar magnetic resonance imaging scan and subsequent 3D-printing. Validation was realized by visual comparison and comparison of total airway volume. To study the feasibility of measuring aerosol deposition, budesonide was aerosolized through the cast and lung dose was expressed as percentage of nominal dose. The airway volumes of the initial magnetic resonance imaging and validation computed tomography scan showed a relative deviation of 0.94%. Lung dose at low flow (1 L/min) was 61.84% and 9.00% at high flow (10 L/min), p < 0.0001. 3D-reconstruction provided an anatomically accurate surrogate of the upper airways of a 32-wk-old premature infant, making the model suitable for future in vitro testing.

Reynolds-stress turbulence model in the KIVA code for engine simulation

A Reynolds-Stress Turbulence Model has been incorporated with success into the KIVA code, a computational fluid dynamics hydrocode for three-dimensional simulation of fluid flow in engines. The newly implemented Reynolds-stress turbulence model greatly improves the robustness of KIVA, which in its original version has only eddy-viscosity turbulence models. Validation of the Reynolds-stress turbulence model is accomplished by conducting pipe-flow and channel-flow simulations, and comparing the computed results with experimental and direct numerical simulation data. Flows in engines of various geometry and operating conditions are calculated using the model, to study the complex flow fields as well as confirm the model’s validity. Results show that the Reynolds-stress turbulence model is able to resolve flow details such as swirl and recirculation bubbles. The model is proven to be an appropriate choice for engine simulations, with consistency and robustness, while requiring relatively low computational effort.

Three members of the connective tissue growth factor family CCN are differentially regulated by mechanical stress

Expression of connective tissue growth factor (CTGF), a member of the CCN gene family, is known to be significantly induced by mechanical stress. We have therefore investigated whether other members of the CCN gene family, including Cyr61 and Nov, might reveal a similar stress-dependent regulation. Fibroblasts growing under stressed conditions within a three-dimensional collagen gel showed at least a 15 times higher level of Cyr61 mRNA than cells growing under relaxed conditions. Upon relaxation, the decline of the Cyr61 mRNA to a lower level occurred within 2 h, and was thus quicker than the response of CTGF. The regulation was fully reversible when stress was reapplied. Thus, Cyr61 represents another typical example of a stress-responsive gene. The level of the Nov mRNA was low in the stressed state, but increased in the relaxed state. This CCN gene therefore shows an inverted regulation relative to that of Cyr61 and CTGF. Inhibition of protein kinases by means of staurosporine suppressed the stress-induced expression of Cyr61 and CTGF. Elevated levels of cAMP induced by forskolin mimicked the effects of relaxation on the regulation of Cyr61, CTGF and Nov. Thus, adenylate cyclase as well as one or several protein kinases might be involved in the mechanoregulation of these CCN genes.

3D Ultrasound imaging--a useful non-invasive tool to detect AV fistulas in transplanted kidneys

BACKGROUND: A precise, non-invasive, non-toxic, repeatable, convenient and inexpensive follow-up of renal transplants, especially following biopsies, is in the interest of nephrologists. Formerly, the rate of biopsies leading to AV fistulas had been underestimated. Imaging procedures suited to a detailed judgement of these vascular malformations are to be assessed. METHODS: Three-dimensional (3D) reconstruction techniques of ultrasound flow-directed and non-flow-directed energy mode pictures were compared with a standard procedure, gadolinium-enhanced nuclear magnetic resonance imaging angiography (MRA) using the phase contrast technique. RESULTS: Using B-mode and conventional duplex information, AV fistulas were localized in the upper pole of the kidney transplant of the index patient. The 3D reconstruction provided information about the exact localization and orientation of the fistula in relation to other vascular structures, and the flow along the fistula. The MRA provided localization and orientation information, but less functional information. Flow-directed and non-flow-directed energy mode pictures could be reconstructed to provide 3D information about vascular malformations in transplanted kidneys. CONCLUSION: In transplanted kidneys, 3D-ultrasound angiography may be equally as effective as MRA in localizing and identifying AV malformations. Advantages of the ultrasound method are that it is cheaper, non-toxic, non-invasive, more widely availability and that it even provides more functional information. Future prospective studies will be necessary to evaluate the two techniques further.

Automatic quantitation of multiple sclerosis lesions on MR images

A two-pronged approach for the automatic quantitation of multiple sclerosis (MS) lesions on magnetic resonance (MR) images has been developed. This method includes the design and use of a pulse sequence for improved lesion-to-tissue contrast (LTC) and seeks to identify and minimize the sources of false lesion classifications in segmented images. The new pulse sequence, referred to as AFFIRMATIVE (Attenuation of Fluid by Fast Inversion Recovery with MAgnetization Transfer Imaging with Variable Echoes), improves the LTC, relative to spin-echo images, by combining Fluid-Attenuated Inversion Recovery (FLAIR) and Magnetization Transfer Contrast (MTC). In addition to acquiring fast FLAIR/MTC images, the AFFIRMATIVE sequence simultaneously acquires fast spin-echo (FSE) images for spatial registration of images, which is necessary for accurate lesion quantitation. Flow has been found to be a primary source of false lesion classifications. Therefore, an imaging protocol and reconstruction methods are developed to generate "flow images" which depict both coherent (vascular) and incoherent (CSF) flow. An automatic technique is designed for the removal of extra-meningeal tissues, since these are known to be sources of false lesion classifications. A retrospective, three-dimensional (3D) registration algorithm is implemented to correct for patient movement which may have occurred between AFFIRMATIVE and flow imaging scans. Following application of these pre-processing steps, images are segmented into white matter, gray matter, cerebrospinal fluid, and MS lesions based on AFFIRMATIVE and flow images using an automatic algorithm. All algorithms are seamlessly integrated into a single MR image analysis software package. Lesion quantitation has been performed on images from 15 patient volunteers. The total processing time is less than two hours per patient on a SPARCstation 20. The automated nature of this approach should provide an objective means of monitoring the progression, stabilization, and/or regression of MS lesions in large-scale, multi-center clinical trials. ^

A longitudinal study of Type-B aortic dissection and endovascular repair scenarios: computational analyses

Conservative medical treatment is commonly first recommended for patients with uncomplicated Type-B aortic dissection (AD). However, if dissection-related complications occur, endovascular repair or open surgery is performed. Here we establish computational models of AD based on radiological three-dimensional images of a patient at initial presentation and after 4-years of best medical treatment (BMT). Computational fluid dynamics analyses are performed to quantitatively investigate the hemodynamic features of AD. Entry and re-entries (functioning as entries and outlets) are identified in the initial and follow-up models, and obvious variations of the inter-luminal flow exchange are revealed. Computational studies indicate that the reduction of blood pressure in BMT patients lowers pressure and wall shear stress in the thoracic aorta in general, and flattens the pressure distribution on the outer wall of the dissection, potentially reducing the progressive enlargement of the false lumen. Finally, scenario studies of endovascular aortic repair are conducted. The results indicate that, for patients with multiple tears, stent-grafts occluding all re-entries would be required to effectively reduce inter-luminal blood communication and thus induce thrombosis in the false lumen. This implicates that computational flow analyses may identify entries and relevant re-entries between true and false lumen and potentially assist in stent-graft planning.

Estimation of aortic pressure waveforms from 4D phase-contrast MRI

Several approaches for the non-invasive MRI-based measurement of the aortic pressure waveform over the heart cycle have been proposed in the last years. These methods are normally based on time-resolved, two-dimensional phase-contrast sequences with uni-directionally encoded velocities (2D PC-MRI). In contrast, three-dimensional acquisitions with tridirectional velocity encoding (4D PC-MRI) have been shown to be a suitable data source for detailed investigations of blood flow and spatial blood pressure maps. In order to avoid additional MR acquisitions, it would be advantageous if the aortic pressure waveform could also be computed from this particular form of MRI. Therefore, we propose an approach for the computation of the aortic pressure waveform which can be completely performed using 4D PC-MRI. After the application of a segmentation algorithm, the approach automatically computes the aortic pressure waveform without any manual steps. We show that our method agrees well with catheter measurements in an experimental phantom setup and produces physiologically realistic results in three healthy volunteers.

A Numerical Investigation of the Gulf Stream and Its Meanders in Response to Cold Air Outbreaks

The three-dimensional Princeton Ocean Model is used to examine the modification of the Gulf Stream and its meanders by cold air outbreaks. Two types of Gulf Stream meanders are found in the model. Meanders on the shoreward side of the Gulf Stream are baroclinically unstable. They are affected little by the atmospheric forcing because their energy source is stored at the permanent thermocline, well below the influence of the surface forcing. Meanders on the seaward side of the stream are both barotropically and baroclinically unstable. The energy feeding these meanders is stored at the surface front separating the Gulf Stream and the Sargasso Seal which is greatly reduced in case of cold air outbreaks. Thus, meanders there reduce strength and also seem to slow their downstream propagation due to the southward Ekman flow. Heat budget calculations suggest two almost separable processes. The oceanic heal released to the atmosphere during these severe cooling episodes comes almost exclusively from the upper water column. Transport of heat by meanders from the Gulf Stream to the shelf, though it is large, does not disrupt the principal balance. It is balanced nicely with the net heat transport in the downstream direction.

ANALYTICAL MODELS OF EXOPLANETARY ATMOSPHERES. I. ATMOSPHERIC DYNAMICS VIA THE SHALLOW WATER SYSTEM

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical, and spherical), rotation, magnetic tension, and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag, and magnetic drag), and magnetic tension are included. The global atmospheric structure is largely controlled by a single key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations, and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to develop decisively physical intuition and as a reference point for three-dimensional magnetohydrodynamic simulations of atmospheric circulation.

ON THE EXISTENCE OF SHOCKS IN IRRADIATED EXOPLANETARY ATMOSPHERES

Supersonic flows are expected to exist in the atmospheres of irradiated exoplanets, but the question of whether shocks develop lingers. Specifically, it reduces to whether continuous flow in a closed loop may become supersonic and if some portions of the supersonic flow steepen into shocks. We first demonstrate that continuous, supersonic flow may exist in two flavors: isentropic and non-isentropic, with shocks being included in the latter class of solutions. Supersonic flow is a necessary but insufficient condition for shocks to develop. The development of a shock requires the characteristics of neighboring points in a flow to intersect. We demonstrate that the intersection of characteristics may be quantified via the knowledge of the Mach number. Finally, we examine three-dimensional simulations of hot Jovian atmospheres and demonstrate that shock formation is expected to occur mostly on the dayside hemisphere, upstream of the substellar point, because the enhanced temperatures near the substellar point provide a natural pressure barrier for the returning flow. Understanding the role of shocks in irradiated exoplanetary atmospheres is relevant to correctly modeling observables such as the peak offsets of infrared phase curves.

THE EFFECTS OF IRRADIATION ON HOT JOVIAN ATMOSPHERES: HEAT REDISTRIBUTION AND ENERGY DISSIPATION

Hot Jupiters, due to the proximity to their parent stars, are subjected to a strong irradiating flux that governs their radiative and dynamical properties. We compute a suite of three-dimensional circulation models with dual-band radiative transfer, exploring a relevant range of irradiation temperatures, both with and without temperature inversions. We find that, for irradiation temperatures T irr lsim 2000 K, heat redistribution is very efficient, producing comparable dayside and nightside fluxes. For T irr ≈ 2200-2400 K, the redistribution starts to break down, resulting in a high day-night flux contrast. Our simulations indicate that the efficiency of redistribution is primarily governed by the ratio of advective to radiative timescales. Models with temperature inversions display a higher day-night contrast due to the deposition of starlight at higher altitudes, but we find this opacity-driven effect to be secondary compared to the effects of irradiation. The hotspot offset from the substellar point is large when insolation is weak and redistribution is efficient, and decreases as redistribution breaks down. The atmospheric flow can be potentially subjected to the Kelvin-Helmholtz instability (as indicated by the Richardson number) only in the uppermost layers, with a depth that penetrates down to pressures of a few millibars at most. Shocks penetrate deeper, down to several bars in the hottest model. Ohmic dissipation generally occurs down to deeper levels than shock dissipation (to tens of bars), but the penetration depth varies with the atmospheric opacity. The total dissipated Ohmic power increases steeply with the strength of the irradiating flux and the dissipation depth recedes into the atmosphere, favoring radius inflation in the most irradiated objects. A survey of the existing data, as well as the inferences made from them, reveals that our results are broadly consistent with the observational trends.