The mountainous areas in the Intermountain West (IMW) of the North America are considered as the major water reservoir for the Western US. Summer evapotranspiration (ET) and soil moisture are key factors affecting the annual water yield in the montane region of the IMW. This research estimated ET of four common vegetation types (aspen, conifer, grass, and sage) and areal soil moisture in an advanced instrumentation site located at the T.W. Daniel Experimental Forest (TWDEF). Among instrumented forest research sites worldwide, TWDEF is one of a few with triplicate measures of meteorological parameters, radiation, and soil moisture within four common vegetation types in the IMW. This unique dataset enables study and understanding of the ecological and hydrological responses to climate change in Utah and the IMW region. In a second phase of this study, summer water uses from the four common vegetation types were simulated using a numerical simulation model, Hydrus-1D. The simulation was informed by soil moisture measurements at three depths (0.1 m, 0.25 m, and 0.5 m) and by ET measured from an eddy covariance tower. The results confirmed the value of numerical simulations as a viable alternate method to estimated ET where no direct ET measurements are available. It also provided comparison of water use by these vegetation species including both high and low water years. In the third phase of this study, a comparison was made between the intermediate-scale areal soil moisture measured by a Cosmic-ray neutron probe (CRNP) and the in situ TDT soil moisture network at the TWDEF site. Improved correlations were obtained, especially after shallow rainfall events, by including numerically simulated soil moisture above 0.1 m where no measurements were available. The original CRNP calibration exhibited a dry bias during spring/early summer, leading to the need for a site-specific enhanced calibration, which improved the accuracy of the CRNP soil moisture estimate at the TWDEF site.

Both primary and secondary liver cancers are common and the majority of patients are not eligible for surgical resection or a liver transplant, which are considered the only hope of cure. Mortality rates are high and there is a need for alternative treatment options. New forms of local treatment work best on small tumours; large ones, however, remain difficult to treat. Hyperthermia involves heating tumours to 40??-44?? C. The aim is to heat the entire tumour without damaging the surrounding normal tissue. Treating deep seated tumours is technically challenging. Ferromagnetic embolisation hyperthermia (FEH) is a novel method of treating liver tumours. Magnetic microspheres are infused into the hepatic artery and lodge primarily in the tumour periphery. An applied alternating-current magnetic field causes the microspheres to heat. Animal experiments have shown that this is a promising technique. There is a need for modelling of FEH prior to commencement of clinical trials. Analytical and numerical models of tumour heating during FEH treatment are presented here. The models help predict the temperature distributions that are likely to arise during treatment and give insight into the factors affecting tumour and liver heating. The models incorporate temperature-dependent thermal properties and blood perfusion rates of the tissues and a heterogeneous clustering of microspheres in the tumour periphery. Simulations show that the poorly perfused tumours heat preferentially while the liver is effectively cooled by blood flow from the portal vein. A peripheral distribution of heat sources produces a more even temperature field throughout the tumour, compared to a heat source that is centred within the tumour core. Large tumours reach higher temperatures and have higher heating rates, supporting experimental findings. Using temperature-dependent, rather than constant, values for thermal conductivities and blood perfusion rates results in higher temperatures within the tumour. The uneven clustering of microspheres in the tumour periphery leads to a more heterogeneous temperature distribution in the core, but it has less of an effect on the wellperfused liver. The results show that FEH has the potential to effectively treat liver tumours and the technique merits further investigation.

This thesis uses results of systematic numerical models to argue that externally inactive salt structures, which are potential targets for radioactive waste disposal, might be internally active due to the presence of dense layers or blocks within a salt layer. The three papers that support this thesis use the Gorleben salt diapir (NW Germany), which was targeted as a future final repository for high-grade radioactive waste, as a general guideline. The first two papers present systematic studies of the parameters that control the development of a salt diapir and how it entrains a dense anhydrite layer. Results from these numerical models show that the entrainment of a dense anhydrite layer within a salt diapir depends on four parameters: sedimentation rate, viscosity of salt, perturbation width and the stratigraphic location of the dense layer. The combined effect of these four parameters, which has a direct impact on the rate of salt supply (volume/area of the salt that is supplied to the diapir with time), shape a diapir and the mode of entrainment. Salt diapirs down-built with sedimentary units of high viscosity can potentially grow with an embedded anhydrite layer and deplete their source layer (salt supply ceases). However, when salt supply decreases dramatically or ceases entirely, the entrained anhydrite layer/segments start to sink within the diapir. In inactive diapirs, sinking of the entrained anhydrite layer is inevitable and strongly depends on the rheology of the salt, which is in direct contact with the anhydrite layer. During the post-depositional stage, if the effective viscosity of salt falls below the threshold value of around 1018-1019 Pa s, the mobility of anhydrite blocks might influence any repository within the diapir. However, the internal deformation of the salt diapir by the descending blocks decreases with increase in effective viscosity of salt. The results presented in this thesis suggest that it is highly likely that salt structures where dense and viscous layer/blocks are present undergo an internal deformation processes when these dense blocks start sinking within the diapir. Depending on size and orientation of these blocks, deformation pattern is significantly different within the diapir. Furthermore, model results applied to the Gorleben diapir show that the rate of descent of the entrained anhydrite blocks differs on different sides of the diapir. This suggests that if the anhydrite blocks descent within the Gorleben diapir, they initiate an asymmetric internal flow within it.

Increasing attention has been directed towards the development of optically transparent and mechanically flexible thin film transistors (TFTs) and associated circuits based on the transition metal oxides. These flexible see-through structures offer reduced weight, potential low-cost fabrication, and high performance compared to commonly used hydrogenated amorphous silicon (a-Si:H) in applications for large-area electronics and displays. As these emerging technologies evolve towards commercialization, a thorough investigation of the impacts of the thermo-mechanical stress and strain and their effects on the electrical and mechanical stability of the flexible microelectronic devices have become increasingly necessary. However, not much progress has been reported in this area, and the numerical modeling of the flexible transistors with the Finite Element Method (FEM) would provide unique insight to the design and operation of the flexible TFTs. In this thesis, numerical models of flexible TFTs are built up by COMSOL Multiphysics and compared with analytical models to reach the best agreement between the experimental measurements and the numerical analyses. These simulations provide additional insight into the local stress induced strain within the device due to both intrinsic and applied stress. It was shown that the thermal and mechanical impacts on the TFT performance can be reduced by placing the vital active layer of the flexible device near the neutral mechanical plane or by proper designing the device structure and processing conditions based on the data derived from the numerical models. The mathematical analysis and numerical simulation will be used to improve the electrical and mechanical performance and the reliability of the transistors for flexible applications.

Modeling interstitial fluid flow involves processes such as fluid diffusion, convective transport in the extracellular matrix, and extravasation from blood vessels. In all of these processes, computational fluid dynamics can play a crucial role in elucidating the mechanisms of fluid flow in solid tumors and surrounding tissues. To date, microvasculature flow modeling has been most extensively studied with simple tumor shapes and their capillaries at different levels and scales. With our proposed numerical model, however, more complex and realistic tumor shapes and capillary networks can be studied. First, a mathematical model of interstitial fluid flow is developed, based on the application of the governing equations for fluid flow, i.e., the conservation laws for mass and momentum, to physiological systems containing solid tumors. Simulations of interstitial fluid transport in a homogeneous solid tumor demonstrate that, in a uniformly perfused tumor, i.e., one with no necrotic region, the interstitial pressure distribution results in a non-uniform distribution of drug particles. Pressure distribution for different values of necrotic radii is examined, and two new parameters, the critical tumor radius and critical necrotic radius, are defined. In specific ranges of these critical dimensions the interstitial fluid pressure is relatively lower, which in turn leads to a diminished opposing force against drug movement and a subsequently higher drug concentration and potentially enhanced therapeutic effects. In this work, the numerical model of fluid flow in solid tumors is further developed to incorporate and investigate non-spherical tumor shapes such as prolate and oblate ones. Using this enhanced model, tumor shape and size effects on drug delivery to solid tumors are then studied. Based on the assumption that drug particles flow with the interstitial fluid, the pressure and velocity maps of the latter are used to illustrate the drug delivery pattern in a solid tumor. Additionally, the effects of the surface area per unit volume of the tissue, as well as vascular and interstitial hydraulic conductivity on drug delivery efficiency, are investigated. Using a tumor-induced microvasculature architecture instead of a uniform distribution of vessels provides a more realistic model of solid tumors. To this end, continuous and discrete mathematical models of angiogenesis were utilized to observe the effect of matrix density and matrix degrading enzymes on capillary network formation in solid tumors. Additionally, the interactions between matrix-degrading enzymes, the extracellular matrix and endothelial cells are mathematically modeled. Existing continuous and discrete models of angiogenesis were modified to impose the effect of matrix density on the solution. The imposition has been performed by a specific function in movement potential. Implementing realistic boundary and initial conditions showed that, unlike in previous models, the endothelial cells accelerate as they migrate toward the tumor. Now, the tumor-induced microvasculature network can be applied to the model developed in Chapters 2 and 3. Once the capillary network was set up, fluid flow in normal and cancerous tissues was numerically simulated under three conditions: constant and uniform distribution of intravascular pressure in the whole domain, a rigid vascular network, and an adaptable vascular network. First, governing equations of sprouting angiogenesis were implemented to specify the different domains for the network and interstitium. Governing equations for flow modeling were introduced for different domains. The conservation laws for mass and momentum, Darcy’s equation for tissue, and a simplified Navier Stokes equation for blood flow through capillaries were then used for simulating interstitial and intravascular flows. Finally, Starling’s law was used to close this system of equations and to couple the intravascular and extravascular flows. The non-continuous behavior of blood and the adaptability of capillary diameter to hemodynamics and metabolic stimuli were considered in blood flow simulations through a capillary network. This approach provided a more realistic capillary distribution network, very similar to that of the human body. This work describes the first study of flow modeling in solid tumors to realistically couple intravascular and extravascular flow through a network generated by sprouting angiogenesis, consisting of one parent vessel connected to the network. Other key factors incorporated in the model for the first time include capillary adaptation, non-continuous viscosity blood, and phase separation of blood flow in capillary bifurcation. Contrary to earlier studies which arbitrarily assumed veins and arteries to operate on opposite sides of a tumor network, the present approach requires the same vessel to run and from the network. Expanding the earlier models by introducing the outlined components was performed in order to achieve a more-realistic picture of blood flow through solid tumors. Results predict an almost doubled interstitial pressure and are in better agreement with human biology compared to the more simplified models generally in use today.

Of all the trace elements emitted during coal combustion, mercury is most problematic. Mercury from the atmosphere enters into oceanic and terrestrial waters. Part of the inorganic Hg in water is converted into organic Hg (CH3Hg), which is toxic and bioaccumulates in human and animal tissue. The largest source of human-caused mercury air emissions in the U.S is from combustion coal, a dominant fuel used for power generation. The Hg emitted from plants primarily occurs in two forms: elemental Hg and oxidized Hg (Hg2+). The coal chlorine content and ash composition, gas temperature, residence time and presence of different gases will decide the speciation of Hg into Hg0 and Hg2+. For Wyoming coal the concentrations of mercury and chlorine in coal are 120ppb and 140ppb. In order to understand the basic process of formulation of HgCl2 and Hg0 a numerical model is developed in the current work to simulate in the detail i) heating ii) transient pyrolysis of coal and evolution of mercury and chlorine, iii) gas phase oxidation iv) reaction chemistry of Hg and v) heterogeneous oxidation of carbon during isolated coal particle combustion. The model assumes that mercury and chlorine are released as a part of volatiles in the form of elemental mercury and HCl. Homogenous reaction are implemented for the oxidation of mercury. Heterogeneous Hg reactions are ignored. The model investigates the effect of different parameters on the extent of mercury oxidation; particle size, ambient temperature, volatile matter, blending coal with high chlorine coal and feedlot biomass etc,. Mercury oxidation is increased when the coal is blended with feedlot biomass and high chlorine coal and Hg % conversion to HgCl2 increased from 10% to 90% when 20% FB is blended with coal. The ambient temperature has a negative effect on mercury oxidation, an increase in ambient temperature resulted in a decrease in the mercury oxidation. The percentage of oxidized mercury increases from 9% to 50% when the chlorine concentration is increased from 100ppm to 1000ppm. When the temperature is decreased from 1950 K to 950 K, the percentage of mercury oxidized increased from 3% to 27%.

A combination of in situ sampling and numerical modeling was used to investigate the effects of mesoscale (<50 km) circulation patterns and stratification on the evolution of hypoxia on the Louisiana Shelf. Temperature, salinity, and dissolved oxygen concentrations records reveal the presence of an alongshelf meander, which is manifested vertically and horizontally as a wave-like distribution of the properties in the water column. The observations suggest the meander is a ubiquitous characteristic of the shelf with alongshore spatial scale approximately 50 km and less, which is consistent with the locations of sandy shoals along the coast and the local deformation radius. Twelve numerical experiments using an idealized three-dimensional shelf circulation model were performed to evaluate the relative importance of the variable bottom topography and freshwater forcing on the development, evolution, and scales of the dynamic instabilities. The inclusion of the shoals into the bottom topography showed the development of the dynamic instabilities as the flow passed over the shoals and downstream. Introduction of fresh water onto the shelf resulted in greater salinity differences, and, as a consequence in the formation of the dynamically unstable salinity fronts along the plume edge. The combination of the freshwater forcing and shoaling topography produced competing and complex interactions. Six numerical experiments were analyzed in order to investigate the effect of dynamic instabilities on spatial and temporal patterns of dissolved oxygen concentrations along the shelf. Although a linear relationship between Brunt-Väisälä frequency and dissolved oxygen deficit was expected, a nonlinear loop-like relationship was discovered that reflects the response of biochemical properties to the alongshelf variability of the density field. Comparison of the numerical modeling runs to observations of density and dissolved oxygen concentrations on the Louisiana Shelf reinforces the importance of physical processes such as topographic steering and/or freshwater forcing on the alongshore distribution of physical and biochemical properties. It suggests that the time scales of respiration (~3 days) and buoyancy transfer processes (~5-7 days), associated with the physical processes that are responsible for water column stability and ventilation, are similar to the time scales associated with the benthic respiration rates.

Carbon Dioxide (CO₂) sequestration requires that the injected CO₂ be permanently trapped in the subsurface and not leak from the target location. To accomplish this, it is important to understand the main mechanisms associated with CO₂ flow and transport in the subsurface once CO₂ is injected. In this work CO₂ seeps at the Crystal Geyser site were studied using modeling and simulation to determine how CO₂ geochemically reacts with formation brines and how these interactions impact the migration of CO₂. Furthermore different scenarios for CO₂ migration and seepage along the Grand Wash fault are studied and the possible outcomes for these different scenarios are documented. The GEM (Generalized Equation-of-State Model) from CMG Ltd. was used to perform the simulation studies. A 2-D model was built without geochemical reactions to mainly study the mechanism associated with dissolution of CO₂ gas. The process of CO₂ release from the brine as the fluid mixture flows up along the fault was modeled. Then, 3-D models with geochemical reactions were built for CO₂ migration corresponding to two different sources of CO₂ - deep crustal ₂ and CO₂-dissolved in groundwater. In both these cases, CO₂ reacted with the aqueous components and minerals of the formation and caused carbonate mineralization. In the case of deep crustal CO₂ source, there were vertical patterns of calcite mineralization simulated along the fault that indicated that calcite mineralization might be localized to isolated vertical flow paths due to vertical channeling of CO₂ from the crust. In the case of CO₂-dissolved groundwater flowing along the sandstone layers, calcite mineralization is spread over the entire fault surface. In this case, the groundwater flow is interrupted by the fault and there is vertical flow along the fault until a permeable sandstone layer is encountered on the other side of the fault. This vertical migration of CO₂-saturated brine causes a release in pressure and subsequent ex-solution of CO₂. As a result, modeling allowed us to establish difference in surface expression of CO₂ leakage due to two different CO₂ migrations scenarios along the fault and helped develop a scheme for selecting appropriate model for CO₂ leakage based on surface observation of travertine mounds. A key observation at the Crystal Geyser site is the lateral migration of CO₂ seep sites over time. These migrations have been confirmed by isotope studies. In this modeling study, the mechanism for migration of seep sites was studied. A model for permeability reduction due to precipitation of calcite was developed. It is shown using percolation calculations that flow re-routing due to permeability alterations can result in lateral migration of CO₂ seeps at rates comparable to those established by isotope dating. / text

At present, most of the energy power consumed in the world is produced by fossil fuel combustion, which has raised increasing interest in renewable energy technologies, non-conventional oil and gas reservoirs, and nuclear power. Innovative nuclear fuels and reactors depend on the economical and environmental impacts of waste management. Disposals in mined geological formations are viewed as potential consolidated storage facilities before final disposition. Different stress paths during construction result in different kinds of failure mechanisms, which alter rock strength and induce anisotropy of rock elastic properties. Crack propagation in rock can be originated by these engineering activities (excavation, drilling, mining, building overburden), or by changes of the natural environment (tectonic processes, erosion or weathering). Damage is a mathematical variable that can represent a variety of microstructure changes, such as crack density, length, aspect ratio and orientation. The framework of Continuum Damage Mechanics allows modeling the resulting reduction in strength and stiffness, as well as the associated stress-induced anisotropy and irreversible deformation. This work presents a modeling framework for anisotropic crack propagation in rock, in conditions of stress typical of geological storage and oil and gas extraction. Emphasis is put on the prediction of the damage zone around cavities and ahead of pressurized fracture tips. An original model of anisotropic damage, the Differential Stress Induced Damage (DSID) model, is explained. The Drucker-Prager yield function is adapted to make the damage threshold depend on damage energy release rate and to distinguish between tension and compression strength. Flow rules are derived with the energy release rate conjugate to damage, which is thermodynamically consistent. The positivity of dissipation is ensured by using a non-associate flow rule for damage, while nonelastic deformation due to damage is computed by an associate flow rule. Stress paths simulated at the material point illustrate damaged stiffness and deformation variations in classical rock mechanics tests. The maximum likelihood method was employed to calibrate and verify the DSID model against stress-strain curves obtained during triaxial compression tests and uniaxial compression tests performed on clay rock and shale. Logarithmic transformation, normalization and forward deletion allowed optimizing the formulation of the DSID model, and reduce the number of damage constitutive parameters from seven to two for clay rock. The DSID model was implemented in ABAQUS Finite Element (FE) software. The iterative scheme was adapted in order to account for the non-linearities induce both by damage and damage-induced deformation. FE simulations of laboratory tests capture size an intrinsic anisotropy effects on the propagation of damage in rock. Smeared DSID zones representing shale delamination planes avoid some convergence problems encountered when modeling discontinuities with debonded contact surface elements. FE simulations of tunnel excavation, fracture propagation and borehole pressurization were performed to illustrate the evolution of the damage zone and the impact on energy dissipation, anisotropy of deformation, and loss of stiffness. Future work will focus on coupling the propagation of fractures with the evolution of the damage process zone, and on the transition from continuum damage to discrete fracture upon crack coalescence.

While several lines of evidence strongly hint at the biological importance of impact-induced hydrothermal systems during the impact cataclysm at ~3.9 Ga, these systems are not well understood. There is unambiguous evidence of hydrothermal activity at many terrestrial craters, but the available samples represent a very limited number of crater diameters and locations within the crater. Therefore, computer models are crucial for learning how impact-induced hydrothermal systems work, how long they last, and whether they provide suitable environments for thermophilic microorganisms. This dissertation presents detailed simulations of hydrothermal activity at the terrestrial craters Chicxulub and Sudbury, as well as at range of crater sizes on early Mars. A well-established computer code HYDROTHERM was used. The models for terrestrial craters were constrained by seismic, magnetic, and gravity surveys, as well as petrological, mineralogical, and chemical analyses of samples (by others).Sudbury crater is ~180 km in diameter, and 1.85 Ga. Simulation results indicate that a hydrothermal system at Sudbury crater remained active for several hundred thousand to several million years, depending on assumed permeability, and produced habitable volumes of up to ~20,000 km^3.Chicxulub crater is also ~180-km in diameter, but only 65 Ma. The lifetime of the hydrothermal system ranges from 1.5 Ma to 2.3 Ma depending on assumed permeability. The temperatures and fluxes observed in the model are consistent with alteration patterns observed by others in borehole samples.Another set of simulations modeled post-impact cooling of hypothetical craters with diameters of 30, 100, and 180 km in an early Martian environment. System lifetimes, averaged for all permeability cases examined, were 67,000 years for the 30-km crater, 290,000 years for the 100-km crater, and 380,000 for the 180-km crater. Also, an ap-proximation of the thermal evolution of a Hellas-sized basin (~2000 km) suggests poten-tial for hydrothermal activity for ~10 Myr after the impact. The habitable volume reached a maximum of ~6,000 km^3 in the 180-km crater model.Possible morphological and mineralogical signs of hydrothermal activity in Martian craters were observed, both in this work and by others. These observations, while by no means definitive, are generally consistent with model predictions.