The effects of different shearing processes and sample sizes on the fluid flow anisotropy and its impact on particle transport process in rough rock fractures are significant factors that need to be considered in the performance and safety assessments of underground nuclear waste repositories. The subjects, however, have not been adequately investigated previously in either laboratory experiments or numerical modeling. This thesis addresses these problems using numerical modeling approaches. The modeling consists of two parts: 1) fluid flow simulations considering more complex but realistic flow boundary conditions during shear processes that cannot be realized readily in laboratory experiments, using digitalized fracture surfaces scanned in the laboratory, so that anisotropic fluid flow induced by shearing with channeling phenomenon can be directly simulated and quantified; 2) particle tracking simulations to demonstrate the impacts of such channeling effects on characteristic properties of particle transport. The numerical method chosen for the simulations is the Finite Element Method (FEM). Scale effects were considered in the simulations by using fracture surface samples of different sizes. The distributions of fracture aperture during shear were obtained by numerically generating relative translational and rotary movements between two digitalized surfaces of a rock fracture replica without considering normal loading. From the evolutions of the aperture distributions during the shearing processes, the evolutions of the transmissivity fields were determined by assuming the validity of the cubic law locally. A geostatistical approach was used to quantify the scale effects of the aperture and transmissivity fields. The fluid flow was simulated using different flow boundary conditions, corresponding to translational and rotary shear processes. Corresponding to translational shear (with a 1 mm shear displacement interval up to a maximum shear displacement of 20 mm), three different flow patterns, i.e., unidirectional (flow parallel with and perpendicular to the shear direction), bi-directional and radial, were taken into account. Corresponding to rotary shear (with a 0.5o shear angle interval up to 90o), only the radial flow pattern was considered. The particle transport was simulated using the Particle Tracking Method, with the particles motion following the fluid velocity fields during shear, as calculated by FEM. For the unidirectional particle transport, the breakthrough curves were analyzed by fitting to an analytical solution of 1-D advection-dispersion equation. The dispersivity, Péclet number and tracer velocity, as well as their evolutions during shear, were determined numerically. The results show that the fracture aperture increases anisotropically during translational shear, with a more pronounced increase in the direction perpendicular to the shear displacement, causing significant fluid flow channelling. A more significant increase of flow rate and decrease in travel time of the particles in the direction perpendicular to the shear direction is predicted. The particle travel time and characteristics are, correspondingly, much different when such effects caused by shear are included. This finding may have an important impact on the interpretation of the results of coupled hydro-mechanical and tracer experiments for measurements of hydraulic properties of rock fractures, because hydraulic properties are usually calculated from flow test results along the shear directions, with the effects of the significant anisotropic flow perpendicular to the shear direction ignored. The results also show that safety assessment of a nuclear repository, without considering the effects of stress/deformation of rocks on fluid flow and transport processes, may have significant risk potential. The results obtained from numerical simulations show that fluid flow through a single rough fracture changes with increasing sample size, indicating that representativehydro-mechanical properties of the fractures in the field can only be accurately determined using samples of representative sizes beyond their stationarity thresholds. / QC 20101207

This work is concerned with the development of computational methods for approximating sensitivities of solutions to boundary value problems. We focus on the continuous sensitivity equation method and investigate the application of adaptive meshing and smoothing projection techniques to enhance the basic scheme. The fundamental ideas are first developed for a one dimensional problem and then extended to 2-D flow problems governed by the incompressible Navier-Stokes equations. Numerical experiments are conducted to test the algorithms and to investigate the benefits of adaptivity and smoothing. / Ph. D.

Bone is the second most transplanted tissue in the body, with approximately 2.2 million bone graft procedures performed annually worldwide. Currently, autogenous bone is the gold standard for bone grafting due to its ability to achieve functional healing; however, it is limited in supply and results in secondary injury at the donor site. Tissue engineering has emerged as a promising means for the development of new bone graft substitutes in order to overcome the limitations of the current grafts. In this research project, the specific approach for bone tissue engineering involves seeding osteoprogenitor cells within a biomaterial scaffold then culturing this construct in a biodynamic bioreactor. The bioreactor imparts osteoinductive mechanical stimuli on the cells to stimulate the synthesis of an extracellular matrix rich in osteogenic and angiogenic factors that are envisioned to guide bone healing in vivo. Fluid flow, which exerts a hydrodynamic shear stress on adherent cells, has been identified as one of the strongest stimuli on bone cell behavior. It has been shown to enhance the deposition of osteoblastic matrix proteins in vitro, and is particularly important for the delivery of oxygen and nutrients to cells within large scaffolds suitable for bone tissue regeneration. In particular, dynamic flow profiles have been shown to be more efficient at initiating mechanotransductive signaling and enhancing gene expression of osteoblastic cells in vitro relative to steady flow. However, the molecular signaling mechanisms by which bone cells convert hydrodynamic shear stress into biochemical signals and express osteoblastic matrix proteins are not fully understood. Therefore, the overall goal of this research project was to determine the effect of dynamic fluid flow on mechanotransductive signaling and expression of bioactive factors and bone matrix proteins. In the first study, an intermittent flow regimen, in which 5 min rest periods were inserted during fluid flow, was examined. Results showed that signaling molecules, mitogen activated protein kinases (MAPKs) and prostaglandin E2, were modulated with the flow regimen, but that expression of bone matrix proteins, collagen 1α1, osteopontin, bone sialoprotein (BSP), and osteocalcin (OC), were similar under continuous and intermittent flow. Thus, this study suggested that variation of the flow regimen modulates mechanotranductive signaling. In the second study, four flow conditions were examined: continuous flow, 0.074 Hz, 0.044 Hz, and 0.015 Hz pulsatile flow. This study demonstrated that pulsatile flow enhances expression of BSP and OC over steady flow. Similarly, bone morphogenetic protein (BMP)-2 and -7 were enhanced with pulsatile flow, while BMP-4 was suppressed with all flow conditions, suggesting that the mechanism by which fluid flow enhances bone matrix proteins may involve the induction of BMP-2 and -7, but not BMP-4. In the third study, the molecular mechanism by which fluid flow simulates expression of BMPs was examined. Results from this study suggest that this mechanism may involve activation of MAPKs, but BMP-2, -4, and -7 are regulated through multiple different signaling pathways. Overall, the results from this research demonstrate that dynamic flow modulates mechanotransductive signaling and expression of osteoblastic matrix proteins by osteoblast cells. In particular, BMPs, important for formation in vivo, were shown to be induced by fluid flow. Therefore, this work may be beneficial in understanding and developing 3D perfusion culture systems for the creation of a clinically effective engineering bone tissue. / Ph. D.

The human brain is a complex organ that is responsible for regulating all the physiological processes in the body, ranging from memory to movement. As humans age, the brain goes through a variety of changes including a reduction in glymphatic waste clearance and increase in glial reactivity. Two neurological conditions that affect individuals over the age of 65 include glioblastoma (GBM) and Alzheimer's disease (AD). Interestingly, patients with GBM do not present with AD and vice versa. Both conditions are characterized by a disruption in interstitial fluid flow (IFF) and an increase in neuroinflammation. Throughout the following dissertation, we examined the role of IFF in AD and GBM progression using a three-sided approach (in vivo, in vitro, and in silico). Increased IFF underlies glioma invasion into the surrounding tumor microenvironment (TME) in GBM. We used a 3D hydrogel model of the GBM TME to examine potential pathways by which astrocytes and microglia contribute to glioma invasion. A reduction in IFF contributes to accumulation of the toxic protein amyloid beta (Aβ) in AD. We sought to create a novel, patient-inspired model of the AD hippocampus for examination of the relationship between IFF and Aβ clearance. Human AD and unaffected control hippocampal brain samples were stained for markers of neurons, astrocytes, microglia and Aβ. The percentage of each cell population in the CA1 region of the hippocampus was calculated. We also analyzed the amount and characteristics of the Aβ aggregates present in this hippocampal region. Pearson correlation analysis was completed to assess the relationships between the various cell populations, Aβ load, and patient descriptors. The cell ratios gleaned from the patient samples were incorporated into a novel, 3D hydrogel model of the AD hippocampus. This model features a hydrogel mixture like the native brain extracellular matrix (ECM) and allows for the application of IFF and Aβ. To our knowledge, we are the first group to create a patient-specific triculture model of the AD hippocampus, which is the main site of Aβ aggregation in the AD brain. We used this model to examine the relationship between IFF-mediated Aβ clearance and glial reactivity. The last aim of this dissertation was to create a computational model for examining Aβ binding within the ECM and the effects of IFF on Aβ clearance. In vitro experiments were conducted to generate 3D renderings of glial cells and to determine relevant parameters for our model. Throughout this work, we discuss the relationship between disruption in IFF and glial reactivity in the context of GBM and AD. / Doctor of Philosophy / The human brain is a complex organ that is responsible for regulating all the physiological processes in the body, ranging from memory to movement. As humans age, the brain goes through a variety of changes including a reduction in brain waste removal and an increase in inflammation. Two neurological conditions that affect individuals over the age of 65 include glioblastoma (GBM) and Alzheimer's disease (AD). Interestingly, patients with GBM do not present with AD and vice versa. Both conditions are characterized by a disruption in brain interstitial fluid flow (IFF) and an increase in neuroinflammation. Throughout the following dissertation, we examined the role of IFF in AD and GBM progression using a three-sided approach including analysis of mouse and human tissues, engineered cell models, and computational methods. Specific interactions between brain cell types and their relationships with glioma invasion were examined using a 3D cell model that mimics the brain. Through the work presented here, we also sought to create a novel cell model of the hippocampus region located in the AD brain. We quantified the various cell types in the hippocampus of AD patient samples and incorporated this information into our hydrogel model. The resulting model features three brain cell types (astrocytes, microglia, and neurons) that are added at patient relevant ratios, a matrix that mimics the native brain scaffold, and allows for the application of IFF. In the AD brain there is a reduction in brain waste removal that leads to accumulation of the toxic protein amyloid beta (Aβ). We were successfully able to incorporate this protein within our model so we could assess the relationship between IFF and Aβ removal from the brain. We further studied this relationship using a new computational model of Aβ accumulation in the brain. Throughout this work, we discuss the connection between disrupted IFF and neuroinflammation in the context of GBM and AD.

Glioblastoma (GBM) is the most aggressive and malignant glioma. It accounts for 48.6% of all primary, malignant gliomas with a median survival of 15 months. Infiltration into the surrounding parenchyma is a hallmark of GBM. Radiotherapy is used to address the invasion; however, recent studies have implicated that radiation contributes to increased invasiveness of glioma. Although the effect of radiation on cells has been studied extensively, its effect on the transport of fluid is not well characterized. Transport in the brain which has significant roles in physiology, GBM pathophysiology, and GBM treatment. Thus, understanding the effect of radiation on transport within the lesion and surrounding interstitium will be beneficial in characterizing the effects of radiotherapy in GBM patients. This dissertation seeks to explore the relationship between radiation, transport, and movement of glioma cells and includes the following: 1) Characterizing in vitro motility metrics of glioma stem cell lines in and relating them to in vivo invasion. 2) Studying the effect of radiation on motility, flow-mediated invasion, extracellular matrix components, and transport within the lesion and interstitium. 3) Assessing transport in clinical images and relating transport parameters to progression of GBM. 4) Developing a novel pipeline for applying vector field topology to the study of interstitial fluid flow in glioma. Surprisingly, we found that motility metrics in vitro have a negative correlation trend with in vivo invasion. Next, we found that radiation causes a transient increase in advective flow, and a more sustained decrease in diffusivity in a murine glioma model. Tenascin C was found to correlate significantly with invasion and diffusivity, indicating that it might be a link between radiation, transport, and invasion. Furthermore, interstitial fluid flow was calculated and assessed in clinical images. This showed that interstitial fluid flow velocity magnitude in the tumor correlates with overall survival in GBM patients. Lastly, vector field topology was introduced as a novel method of studying transport that provides more detailed information to identify potential drivers of transport within a flow field. Altogether, this work presents novel insight into the effects of radiation on invasion and transport in GBM. Hopefully, this work can provide a foundation to build upon in efforts of improving treatment planning and clinical outcomes for GBM patients. / Doctor of Philosophy / Glioblastoma (GBM) is the most aggressive glioma. It accounts for 48.6% of all primary, malignant gliomas with a median survival of 15 months. The movement of cancer cells into the surrounding tissue is a defining factor of GBM. Radiotherapy is used after surgery to treat the remaining cancer cells in tissue surrounding the tumor; however, recent studies have implicated that radiation contributes to increased movement of glioma into surrounding tissue. Although the effect of radiation on cells has been studied extensively, its effect on transport of fluid is not well characterized. Interstitial fluid flow in the brain has significant roles in healthy bodily functions, GBM disease state, and GBM treatment. Thus, understanding the effect of radiation on transport within the tumor and surrounding tissue is beneficial in better characterizing the effects of radiotherapy. This dissertation seeks to explore the relationship between radiation, transport, and movement of glioma cells and includes the following: 1) Characterizing in vitro motility metrics of glioma cells in and relating them to in vivo movement into healthy tissue. 2) Studying the effect of radiation on motility, flow-mediated infiltration into healthy tissue, tissue matrix components, and fluid flow within the tumor and surrounding tissue. 3) Assessing transport in clinical images and relating transport parameters to progression of GBM. 4) Developing a novel pipeline for applying vector field topology to the study of interstitial fluid flow in glioma. Surprisingly, we found that motility metrics in vitro have a negative correlation trend with in vivo invasion. Next, we found that radiation causes a transient increase in flow velocity magnitude, and a more sustained decrease in diffusivity in a murine glioma model. Tenascin C, a component of the tissue matrix, was found to correlate significantly with invasion and diffusivity. This indicates that Tenascin C might be a link between radiation, transport, and invasion. Furthermore, interstitial fluid flow was calculated and assessed in clinical images which showed that interstitial fluid flow velocity magnitude in the tumor correlates with survival. Lastly, vector field topology was introduced as a novel method of studying fluid flow in glioma that provides more detailed information regarding the flow field. Altogether, this work presents novel insight into the effects of radiation on fluid flow and cellular movement in GBM. Hopefully, this work can provide a foundation to build upon in efforts of improving treatment planning and clinical outcomes for GBM patients.

A high resolution spectrally accurate three-dimensional flow solver is developed in order to simulate convection dominated fluid flows. The governing incompressible Navier Stokes equations along with the energy equation for temperature are discretized using a second-order accurate projection method which utilizes Adams Bashforth and Backward Differentiation formula for temporal discretization of the non-linear convective and linear viscous terms, respectively. Spatial discretization is performed using a Fourier/Chebyshev spectral method. Extensive tests on three-dimensional Taylor Couette flow are performed and it is shown that the method successfully captures the different states ranging from formation of Taylor vortices to wavy vortex regime. Next, the code is validated for convection dominated flows through a comprehensive comparison of the results for two dimensional Rayleigh Benard convection with the theoretical and experimental results from the literature. Finally, fully parallel simulations, with efficient utilization of computational resources and memory, are performed on a model three-dimensional axially homogeneous Rayleigh Benard convection problem in order to explore the high Rayleigh number flows and to test the scaling of global properties.

Cette thèse est consacrée au développement d'une méthode numérique visant à modéliser des écoulements dans des milieux poreux fissurés, ainsi qu'à déterminer leur perméabilité effective à partir des avancements théoriques récents. En parallèle, elle a été aussi l'occasion de continuer sur la voie théorique et d'obtenir de nouveaux résultats sur ce plan.Les équations générales gouvernant l'écoulement dans de tels matériaux sont rappelées, et plus particulièrement, la conservation de la masse à l'intersection entre fissures est établie explicitement [132, 139]. Des solutions générales du potentiel sont proposées sous la forme d'une équation intégrale singulière décrivant l'écoulement dans et autour des fissures entourées par une matrice infinie soumise à un champ lointain [136, 139]. Ces solutions représentent le champ de pression dans le milieu infini en fonction de l'infiltration dans les fissures qui prennent en compte complètement l'interaction et l'intersection entre fissures. En considérant le problème d'une fissure super-conductrice, des solutions analytiques de l'écoulement ont été développées et ces solutions sont utilisées comme cas de référence pour valider la solution numérique. De plus, la solution théorique de ce problème dans le cas 3D permet de comparer le modèle d'écoulement de Poiseuille dans une fissure elliptique d'épaisseur nulle et le modèle d'inclusion ellipsoïdale aplatie soumise à l'écoulement de Darcy [140]. Des outils numériques ont été développés en se basant sur la méthode des équations intégrales singulières afin de résoudre les équations générales du potentiel [132, 180]. Cela permet, d'une part, de modéliser l'écoulement stationnaire dans un domaine poreux contenant un grand nombre de fissures et, d'autre part, de proposer une solution semi-analytique de l'infiltration dans une fissure isolée dépendant de la perméabilité de la matrice, de la conductivité de la fissure et de la variable géométrique de la fissure. Cette dernière est l'élément essentiel pour déterminer de la perméabilité effective d'un milieu poreux fissuré en utilisant des schémas d'homogénéisation. Ce modèle auto-cohérent révèle un seuil de percolation qui est alors applicable pour l'estimation de la perméabilité effective d'un matériau contenant un grand nombre de fissures. L'approche par sous-structuration permet l'extension de la solution générale du potentiel, écrite pour un domaine infini, à celle d'un domaine fini [181]. Une solution analytique de l'écoulement dans et autour d'une fissure partiellement saturée est établie, fondée sur la solution analytique pour la fissure super-conductrice. Celle-ci est alors utilisée pour estimer la perméabilité effective d'un milieu poreux fissuré non-saturé [141]. Le modèle de la perméabilité effective est appliqué dans le contexte du stockage géologique du CO2 en vue d'étudier le comportement d'une zone de faille constituée par un noyau argileux et des zones fissurées. La pression d'injection provoque l'augmentation de la pression interstitielle dans le réservoir. Cette surpression engendrée dans le réservoir peut affecter la perméabilité de zones fissurées ce qui conduit à des phénomènes hydromécaniques couplés. Les résultats de simulations numériques permettent d'évaluer, d'une part, le risque de la remontée de la saumure à l'aquifère supérieur, et d'autre part, le risque de l'initiation d'une rupture sur le plan de la faille / This thesis aims to develop a method for numerical modelling of fluid flow through fractured porous media and for determination of their effective permeability by taking advantage of recent results based on formulation of the problem by Singular Integral Equations. In parallel, it was also an occasion to continue on the theoretical development and to obtain new results in this area. The governing equations for flow in such materials are reviewed first and mass conservation at the fracture intersections is expressed explicitly. Using the theory of potential, the general potential solutions are proposed in the form of a singular integral equation that describes the steady-state flow in and around several fractures embedded in an infinite porous matrix under a far-field pressure condition [136, 139]. These solutions represent the pressure field in the whole body as functions of the infiltration in the fractures, which fully take into account the fracture interaction and intersections. Closed-form solutions for the fundamental problem of fluid flow around a single fracture are derived, which are considered as the benchmark problems to validate the numerical solutions. In particular, the solution obtained for the case of an elliptical disc-shaped crack obeying to the Poiseuille's law has been compared to that obtained for ellipsoidal inclusions with Darcy's law [140].The numerical programs have been developed based on the singular integral equations method to resolve the general potential equations [132, 180]. These allow modeling the fluid flow through a porous medium containing a great number of fractures. Besides, this formulation of the problem also allows obtaining a semi-analytical infiltration solution over a single fracture depending on the matrice permeability, the fracture conductivity and the fracture geometry. This result is the important key to upscalling the effective permeability of a fractured porous medium by using different homogeneisation schemes. The results obtained by the self-consistent scheme have been in particular established. The multi-region approach can be used to extend the general potential solution written for the infinite domain to that for a finite domain [181]. A closed-form solution for flow in and around a single partially saturated fracture, surrounded by an infinite matrix subjected to a far-field condition, is also derived combining the solutions for a superconductive fracture and for an imprevious fracture. This solution is then employed to estimate the effective permeability of unsaturated fractured porous media [141].The effective permeability model is applied to study the hydromechanical behaviour of a fault zone constituted by a clay core surrounded by fractured zones in the context of CO2 geological storage. The pressure injection induces an overpressure in the reservoir that may affect the permeability of the fractured zones leading to complexe coupled hydromechanical phenomena. The simulation results allow evaluating the risk of leakage of the reservoir brine to higher aquifers as well as the risk of fault reactivation

<p>Co-axial dripping is one of the many ways to make drops with a core-shell structure for encapsulated materials. However, in systems where the capsule components are not density matched or surfactants are not used, the shell will eventually thin and break if not solidified in time. If the shell fails before solidifying, the core will leak out and result in a non-functional capsule. This study assumes that all capsules will fail once the core has reached 80% eccentricity, meaning a shell region has thinned to 20% of its original thickness (~70 µm). In reality, rupture of the shell depends more on stochastic defects and disturbances, but locally decreasing the shell thickness will increase the probability of capsule rupture. With this assumption, the survival time of a core-shell drop is inversely proportional to the relative velocity of the inner drop, where the greater this relative velocity, the faster the shell phase will thin. Stoke's law is generally used to approximate the speed of a sphere in a fluid. However, this study demonstrates that Stoke's law is insufficient for predicting the inner drop's motion for a compound drop. This is due to internal flows that develop within all fluid drops because of shear forces on the drop’s external face during freefall. For core-shell drops, prior studies report how the inner drop velocity can change in magnitude and direction as a function of its eccentricity, meaning its position within the outer drop. Since previous studies did not analyze this core-shell drop relationship with a 50 vol% core and a high viscosity shell, a model was built in COMSOL Multiphysics to understand how the claims from literature would apply to a previous encapsulation study (Betancourt, 2021). The model was also put through a series of validation tests that confirmed the model’s ability to accurately represent the speed and direction of inner drop motion. The final model configuration was then used to identify the transition point between buoyancy-driven and internal flow-driven failure modes observed during the production of core-shell drops in a previous encapsulation study for phase change materials (Betancourt, 2021). The model results showed how the estimated inner drop velocity was significantly reduced once accounting for the internal flows within the shell phase of a compound drop. While this study does help characterize the motion of an inner drop and could be used to find a material system with a favorable velocity profile, it is still recommended to use an in-air curing system to produce concentric capsules. Achieving a concentric capsule would still require this co-axial dripping setup to be modified significantly. </p> <p>Betancourt-Jimenez, D., Wells, B., Youngblood, J. P., & Martinez, C. J. (2021). Encapsulation of biobased fatty acid amides for phase change material applications. <em>Journal of Renewable and Sustainable Energy</em>, <em>13</em>(6), 064101. https://doi.org/10.1063/5.0072105</p>

<p>  </p> <p>Heat pipes and vapor chambers are passive two-phase heat transport devices that are used for thermal management in electronics. The passive operation of a heat pipe is facilitated by capillary wicking of the working fluid through a porous wick, and thus is subject to an operational limit in terms of the maximum pressure head that the wick can provide. This operational limit, often termed as the capillary limit, is the maximum heat input at which the pressure drop in the wick is balanced by the maximum capillary pressure head; operating a heat pipe or a vapor chamber above the capillary limit at steady-state leads to dryout. It thus becomes important to predict the performance of heat pipes and vapor chambers and explore the parametric design space to provide guidelines for minimized thermal resistance while satisfying this capillary limit. An increasingly critical aspect is to predict the transient thermal response of vapor chambers. Moreover, heat pipes and vapor chambers are extensively being used in electronic systems where the power input is dictated by the end-user activity and is expected to even exceed the capillary limit for brief time intervals. Thus, it is imperative to understand the behavior of heat pipes and vapor chambers when operated at steady and transient heat loads above the capillary limit as dryout occurs. However, review of the literature on heat pipe performance characterization reveals that the regime of dryout operation has been virtually unexplored, and thus this thesis aims to fill this critical gap in understanding.</p> <p>The design for minimized thermal resistance of a vapor chamber or a heat pipe is guided by the relative contribution of thermal resistance due to conduction across the evaporator wick and the saturation temperature gradient in the vapor core. In the limit of very thin form factors, the contribution from the vapor core thermal resistance dominates the overall thermal resistance of the vapor chamber; recent work has focused on working fluid selection to minimize overall thermal resistance in this limit. However, the wick thermal resistance becomes increasingly significant as its thickness increases to support higher heat inputs while avoiding the capillary limit. A thermal resistance network model is thus utilized to investigate the importance of simultaneously considering the contributions of the wick and vapor core thermal resistances. A generalized approach is proposed for vapor chamber design which allows <em>simultaneous</em> selection of the working fluid and wick that provides minimum overall thermal resistance for a given geometry and operating condition. While the thermal resistance network model provides a convenient method for exploring the design space, it cannot be used to predict 3-D temperature fields in the vapor chamber. Moreover, such thermal resistance network models cannot predict transient performance and temperature evolution for a vapor chamber. Therefore, an easy-to-use approach is proposed for mapping of vapor chamber transport to the heat diffusion equation using a set of appropriately defined effective anisotropic thermophysical properties, thus allowing simulation of vapor chamber as a sold conduction block. This effective anisotropic properties approach is validated against a time-stepping analytical model and is shown to have good match for both spatial and transient temperature predictions.</p> <p>Moving the focus from steady-state and transient operation of vapor chambers, a comprehensive characterization of heat pipe operation above capillary limit is performed. Different user needs and device workloads can lead to highly transient heat loads which could exceed the notional capillary limit for brief time intervals. Experiments are performed to characterize the transient thermal response of a heat pipe subjected to heat input pulses of varying duration that exceed the capillary limit. Transient dryout events due to a wick pressure drop exceeding the maximum available capillary pressure can be detected from an analysis of the measured temperature signatures. It is discovered that under such transient heating conditions, a heat pipe can sustain heat loads higher than the steady-state capillary limit for brief periods of time without experiencing dryout. If the heating pulse is sufficiently long as to induce transient dryout, the heat pipe may experience an elevated steady-state temperature even after the heat load is reduced back to a level lower than the capillary limit. The steady-state heat load must then be reduced to a level much below the capillary limit to fully recover the original thermal resistance of the heat pipe. The recovery process of heat pipes is further investigated, and a mechanism is proposed for the thermal hysteresis observed in heat pipe performance after dryout. A model for <em>steady-state</em> heat pipe transport is developed based on the proposed mechanism to predict the parametric trends of thermal resistance following recovery from dryout-induced thermal hysteresis, and the model is mechanistically validated against experiments. The experimental characterization of the recovery process demonstrates the existence of a maximum hysteresis curve, which serves as the worst-case scenario for thermal hysteresis in heat pipe after dryout. Based on the learnings from the experimental characterization, a new procedure is introduced to experimentally characterize the steady-state dryout performance of a heat pipe.</p> <p>To recover the heat pipe performance under steady-state, it has been shown that the heat input needs to be lowered down or <em>throttled</em> significantly below the capillary limit. However, due to the highly transient nature of power dissipation from electronic devices, it becomes imperative to characterize heat pipe recovery from dryout under transient operations. Hence, power-throttling assisted recovery of heat pipe from dryout has been characterized under transient conditions. A minimum throttling time interval, defined as time-to-rewet, is identified to eliminate dryout induced thermal hysteresis using power throttling. Dependence of time-to-rewet on throttling power is explored, and guidelines are presented to advise the throttling need and choice of throttling power under transient conditions. </p> <p>The experimental characterization of heat pipe operation at pulse loads above the capillary limit and power throttling following the pulse load helped define the dryout and recovery performance of a heat pipe. Next, a physics-based model is developed to predict the heat pipe <em>transient</em> thermal response under dryout-inducing pulse load and power throttling assisted recovery. This novel model considers wick as a partially saturated media with spatially and temporally varying liquid saturation, and accounts for the effect of wick partial saturation in heat pipe transport. The model prediction are validated against experiments with commercial heat pipe samples, and it is shown that the model can accurately predict dryout and recovery characteristics, namely time-to-dryout, time-to-rewet, and dryout-induced thermal hysteresis, for heat pipes with a range of wick types, heat pipe lengths and pulse loads above the capillary limit. </p> <p>The work discussed in this thesis opens certain questions that are expected to guide further research in this area. First, the thermal hysteresis mechanism proposed could be further validated with direct visualization of the liquid in a vapor chamber. To achieve this, X-ray microscopy is proposed as a viable option for the imaging <em>in situ</em> wetting dynamics in a vapor chamber. Second, the model developed to predict the dryout and recovery characteristics of the heat pipe can be used to design heat pipe with improved performance under pulse loads and power throttling. Third, novel wick designs can be explored that utilize the understanding developed of governing mechanisms for recovery from dryout, and can eliminate thermal hysteresis at powers closer to capillary limit. Fourth, the modeling approach can be extended to predict dryout and recovery trends in vapor chamber since the heat transfer pathways in a vapor chamber are different than those of a heat pipe. Fifth, and lastly it was observed several times during experiments that some of the heat pipe samples would exhibit complete dryout (sudden catastrophic rise in temperature and thermal resistance at the point of dryout) whereas other samples would exhibit partial dryout (noticeable but small increase in thermal resistance at dryout) at operating powers just above the capillary limit. Exploring and explaining the cause of complete dryout, in particular, would be an extremely valuable contribution to the heat pipe research. </p> <p>The work discussed in this thesis has led to the comprehensive development of a functional and mechanistic understanding of heat pipe operation above the notional capillary limit. The experimental procedures developed in this work are utilized to characterize a heat pipe performance under dryout and recovery. The models based on the mechanistic understanding developed from experimental characterization of dryout and recovery operation of a heat pipe have been experimentally validated and are useful for predicting heat pipe performance under dryout-inducing pulse loads and power-throttling.   </p>

<p>  </p> <p>Cavitation is one of the most common causes of failures in axial piston machines. Due to the detrimental effects that cavitation has on unit performance, it is of important consideration both in the design of new units and in defining the operational limits of existing market products. The work in this thesis aimed to contribute to the current knowledge in both areas, with a focus on design considerations with respect to cavitation scalability, and on operating conditions by measuring cavitation severity under separate and combined inciting parameters. Though the application of unit scaling is common in industry for the design of pump families, there have been no comprehensive attempts to quantify whether cavitation in fluid power units may be adequately accounted for in published scaling laws. In this thesis, the scalability of cavitation phenomena was examined through a CFD scaling study performed using a modified version of the Full Cavitation Model.  Results indicate that linear scaling is consistent in maintaining volumetric efficiency performance within 1% across scaled units up to eight times larger or smaller than the baseline. However, the gas and vapor volume distributions vary significantly between scaled units, due largely to the linear non-scalability of fluid inertia and turbulent factors. Physical exchange between phases within a working fluid was shown to be time-dependent, such that the scaled-down unit exhibits bubble collapse rates up to 30% and 150% greater than the baseline and scaled-up units, respectfully. Considering these effects, the presented work demonstrates a potential for increased cavitation damage area when downscaling a unit and reduced risk in upscaling, despite the scaling law being a reliable indicator for volumetric efficiency. </p> <p>To define a more complete study of cavitation under a variety of operating conditions and inciting parameters, this a new experimental procedure and testing circuit was proposed with focus on repeatability by controlled pressure drops and preliminary quantification of inlet fluid quality. By measuring cavitation conditions under pressure starvation, incomplete filling, and combinations thereof, the direct effect of different inception methods on unit performance was shown to be readily identifiable. Through visualization of the inlet flow, reduction in inlet pressure levels was correlated to fluid cloudiness levels and bubble size, with transparency loss at 0.0 bar_g and transition from bubbly to plug flow at -0.4 bar_g. Incomplete filling-induced cavitation was also shown to be detectable by inlet flow conditions, with a distinct change in bubble coalescence rate when operating under shaft speeds greater than or equal to fill speed for a given inlet pressure. </p>