When planning for a tunnel, the ground conditions in which the tunnel is going to be excavated through will be investigated to different extent. Lack of relevant pre-investigation data or misinterpretations of the available data can cause both economical and/or unexpected stability problems. Weakness zones that are expected to cross the tunnel could be investigated thoroughly with a variety of methods. Refraction seismicity survey and 2D resistivity survey are two geophysical methods that are common in Norway for obtaining information about the rock quality in weakness zones. In this work, a twin tunnel under construction in Bodø (northern Norway) called the Bodøtunnel is studied. The predictions based on the pre-investigation for crossing of some expected weakness zones are compared to the actual conditions encountered during tunneling. Tunneling observations (Geological mapping and photos), rock samples and measurement while drilling (MWD) were used to describe the weakness zones that were encountered during tunneling. Rock samples were collected from two weakness zones and the general rock mass. These samples were tested in a point bearing machine for determination of their uniaxial compressive strength (UCS). These results indicated that the rock samples gathered from the weakness zones had significantly lower UCS than the samples from the rock mass. This was exceedingly clear for the samples of fault rock gathered in connection with a shear zone. The results from this work demonstrate that refraction seismicity had a high success rate for locating weakness zones, with the exception for the crossed narrow zones that were interpreted lacking a shear component. Empirical formulas relating Q-value and UCS with the seismic wave speed were used for calculating these factors for some interesting locations. The empirically calculated UCS was similar to the obtained UCS from the point bearing tests, while the empirically calculated Q-value showed large deviations from the mapped Q-value. The resistivity measurements had a low success rate so far in this project; the reason for this could be disturbances in the ground and the location of the resistivity profiles, which had to adapted to the nearby railroad. It should be noted that only one full resistivity profile has been crossed and the rest of the profiles are expected to be more accurate. Based on the results from the crossed profile(s), the suitability of resistivity survey 2D in urban areas can be brought to question. This work also stumbled upon problems regarding the definition of weakness zones. Shear/fault zones are one of the more common type of weakness zones encountered in tunneling. These kind of zones often consists of different parts. Depending on which parts are regarded as a weakness zone by the responsible engineers, the Q-value might differ due to the SRF. Different scenarios were also evaluated with numerical modeling for the expected remaining major weakness zones. This analysis highlights the importance of differentiation between more fractured zones and zones containing fault rock, such as breccia. The width of the zone had a major impact on the stability while the dip for wide zones had a minor impact on the stability, as long the zones dip is not so small that both tunnels are intersected at the same time. The rock mechanical parameter of the weakness zones that had the most impact on the overall stability was the cohesion.

The widespread recognition of (ultra)high-pressure ((U)HP) metamorphic rocks in orogens worldwide suggests that subduction and exhumation of crustal rocks from mantle depths are normal processes at convergent plate margins. However, the dynamics of these processes, in particular the comparative roles of erosion and crustal extension, and the driving forces of extension during (U)HP rock exhumation, remain controversial. This thesis presents numerical modeling and field/analytical studies that address the geodynamics of crustal subduction and exhumation in two intensely studied orogens, the Alps and the Caledonides. The 2D numerical models show how different scales and durations of orogeny and plate motions can lead to marked contrasts in the style of orogenic growth, crustal subduction, and (U)HP exhumation. In the Western Alps, rapid exhumation (1-3 cm/a) can be explained by local, syn-orogenic extension driven by the buoyant ascent of deforming (U)HP crust from the subduction channel. Later trans-crustal exhumation probably resulted from the combined effects of syn-convergent thrusting, local extension, and erosion. The low temperatures (500-700°C) of Alpine (U)HP metamorphism are attributable to the small size of the orogen and short duration of subduction/exhumation. Contrary to recent suggestions, neither erosion nor absolute extension is required to explain (U)HP exhumation in the Alps. The Western Gneiss Region (WGR) (Norwegian Caledonides), in contrast, can be explained by subduction to (U)HP conditions followed by plate divergence. Gravitational spreading of a thick, hot orogenic wedge leads to a short period of coeval thrusting and extension. Exhumation of (U)HP crust from the subduction channel is achieved by normal-sense shearing along the top of the (U)HP terrane, with minor associated shortening. Trans-crustal exhumation by vertical thinning of the orogenic wedge results from continued absolute extension and erosion. The comparatively high temperatures (700-800°C) achieved by Caledonian (U)HP rocks reflect the orogen's greater size, slower exhumation rates, and possible stalling of the (U)HP terrane at depth. These contrasting models underscore the variety of possible mechanisms responsible for (U)HP exhumation, and represent new benchmarks in the understanding of Alpine and Caledonian tectonics and (U)HP rock exhumation in general.

No description available.

The formation of aeolian ripples has been modeled, quite successfully, using discrete approaches like cellular automaton models. Numerical analysis of continuum models to obtain similar success in modeling ripple evolution, however, has not been studied extensively. A numerical model based on continuum theories expedites calculations, as opposed to discrete approaches which model trajectory of each and every sand grain, and are hence relatively more economical. The numerical analysis strives to contribute to the field of study of aeolian ripple migration by an extensive comparison and discussion of modeled ripple evolution results with those of a particular laboratory based wind-tunnel experiment. This research also endeavors to under- stand the physics behind ripple generation and what parameters to be modified to account for multiple grain sizes. Incorporation of multiple grain sizes would enable us to study the stratigraphy of the generated bed-forms. To obtain smoother and realistic ripple surfaces, a sixth-order compact finite difference numerical scheme is used for spatial derivates and fourth-order Runge-Kutta scheme for time derivates. The boundary conditions incorporated are periodic and the initial condition employed to generate ripple is a rough sand surface. The numerical model is applied to study the effect of varying the angle, at which the sand bed gets impacted by sand grains, on the evolution of ripples. Ripples are analyzed qualitatively and quantitatively by considering the contribution of processes involved in the evolution process. The ripple profiles and the time taken to reach equilibrium state, obtained by numerical experiments, are in close agreement with the ones obtained by the wind-tunnel experiment.

In the present study, methodological approaches to assess the impact of climate change on the thermal and thermo-hydro-mechanical (THM) regimes of the ground in some selected Canadian no-permafrost areas (Ottawa, Sudbury, Toronto) is proposed. A modeling study to evaluate ground temperature variations due to global warming is conducted using TEMP/W software from Geoslope ltd. The effect of future climate change projections, up to 2100, on the ground freeze-thaw cycle frequencies, frost penetration depth, and frost duration is assessed in some selected sites located in the Canadian no-permafrost region. Moreover, three softwares (TEMP/W, SEEP/W and SIGMA/W from Geoslope Ltd) have been used to establish a numerical tool that enable to assess the effect of global warming on THM response of the grounds in the selected Canadian no-permafrost areas. TEMP/W and SEEP/W were coupled in a thermo-hydraulic analysis to assess the impact of global warming on the hydraulic regime of the ground. Afterwards, SEEP/W and SIGMA/W were coupled in a hydraulic-mechanical analysis to study the impact of climate change induced porewater pressures change on the mechanical regime of the ground in some no-permafrost regions. Simulation study to assess the effect ground temperature changes on key geotechnical properties of the soils in the selected sites is conducted by using the aforementioned numerical tool. The change of the porewater pressure changes and distributions in the soil induced by global warming is studied. The effect of climate change on the ground consolidation or settlement in the selected no-permafrost sites is also investigated. Finally, this study provides a simulation of a bridge pile foundation ground to detect the THM changes around the pile structure due to climate warming. The results indicate that climate change will affect the thermal regime of the ground in the selected Canadian no permafrost areas. Ground temperature in the studied no-permafrost regions will likely increase by 2 to 4 C by 2100 due to global warming. Furthermore, the frost penetration depth will be significantly reduced in all study areas. It is also found that the frost duration will experience a gradual reduction with time up to 2100. In addition, the simulation results showed minimal influence of global warming on the porewater pressure distribution and magnitude in the studied grounds. Aligned to this, climate change did not seem to have a significant effect on the consolidation behavior or settlement of the ground in the studied no-permafrost areas. The simulation of the foundation ground confirms the results mentioned above, as temperature changes around the pile structure falls within the same range found in the thermal analysis. Porewater pressure distributions and ground settlement are not significantly affected along the pile perimeter. Overall, the design of pile foundation in the Canadian no-permafrost region will not be significantly affected by climate change up 2100. The tools developed and results obtained will be useful for the geotechnical design of climate-adaptive civil engineering or transportation structures in Canadian no permafrost areas.

The main objective of this research was to create a two dimensional and three dimensional Environmental Fluid Dynamics Code (EFDC) model using Aberdeen Pool of the Tennessee-Tombigbee Waterway for the purpose of determining the differences in their ability to address sediment transport. These objectives were reached in the results with comparisons of water levels, sediment concentrations, shear stress, and bed change. The models produced very similar results for the majority of the sediment transport throughout both models with the overall trend being deposition except in the upper limits of the Tombigbee River. The main differences between the two models are produced from the 2D model being depth averaged and the 3D being able to transport sediment vertically. The results show the 2D model tends to erode less and deposit more whereas the 3D model tends to follow the same pattern except for less deposition with more erosion.

Underground mining has one of the highest fatal injury rates among any of the industries in the United States, which is more than five times the national average of the other industries (MSHA). Many of these incidents take place due to stress redistribution resulting from mine workings. Thus it is very important to develop some tools to predict this failure in advance and prevent any fatalities arising from the failure. The current study uses two tools — numerical modeling and laboratory tomography - to image the stress distribution in laboratory rock samples as they are uniaxially loaded. The discrete element code, PFC3D, is used. The laboratory properties of the rock sample need to be converted to the micro-properties of the particles in the model. Currently no theory exists for this conversion. In the current study an equation has been developed for this process. Based on the users' input, the equation determines the micro-properties for the model. Further, various techniques to study the stress redistribution from these models at the particle level are discussed. Tomography is a non-destructive technique through which the interior of a body can be imaged without penetrating the surface by any physical means. In the current study sensors were attached around the rock sample and tomograms were obtained at certain intervals of the load. Initially, an indentation load was applied on a rectangular block to study the comparison between the stress and the velocity in two dimensions. In the last part of the study three-dimensional tomograms were obtained from the rock samples as they were loaded to failure. / Ph. D.

Carbon-capture and sequestration (CCS) in geologic reservoirs is one strategy for reducing anthropogenic CO2 emissions from large-scale point source emitters. Recent developments have shown that basalt reservoirs are highly effective for permanent mineral trapping on the basis of CO2-water-rock interactions, which result in the formation of carbonate minerals. However, the injection of super-critical CO2 into the subsurface causes a disturbance in the pressure, temperature, and chemical systems within the target reservoir. How the ambient conditions change in response to a CO2 injection ultimately affects the transport and fate of the injected CO2. Understanding the behavior and transport of CO2 within a geologic reservoir is a difficult problem that is only exacerbated by heterogeneities within the reservoir; for example, permeability can be highly heterogeneous and exhibits significant control on the movement of CO2. This work is focused on constraining the permeability uncertainty within a flood basalt reservoir, specifically the Columbia River Basalt Group (CRBG). In order to do so, this dissertation is a culmination of four projects: (1) a geostatistical analysis resulting in a spatial correlation model of regional scale permeability within the CRBG, (2) a Monte Carlo-type modeling studying investigating the effects that permeability uncertainty has on the injectivity and storativity of the CRBG as a storage reservoir, (3) a modeling study utilizing 1-, 2-, and 3-D numerical models to investigate how the thermal signature of the CO2-water system evolves during a CO2 injection, and (4) a Monte Carlo-type modeling study focused on the integrity of the CRBG as a CO2 storage reservoir through a probabilistic assessment of static threshold criteria. / Doctor of Philosophy / The process of capturing CO2 from point-source emitters, such as power plants and injecting that CO2 into a geologic formation is one way to reduce anthropogenic CO2 emissions. Recent field studies have shown that basalt reservoirs may be very effective at permanently storing the injected CO2 making them a secure geologic formation to store the CO2. However, basalt reservoirs can be highly fractured, which causes the properties of the reservoir (e.g. permeability, porosity, etc.) to be nonuniform. Having nonuniform reservoir properties creates uncertainty when planning a large-scale CO2 injection. This research is focused on understanding and constraining the uncertainty of nonuniform reservoir properties associated with a large-scale CO2 injection. The work presented utilizes a geostatistical analysis of permeability to inform a variety of numerical models to study how nonuniform reservoir properties affect CO2 injection rate, how much CO2 can be stored, how the pressure and temperature of the reservoir changes, and how secure the storage reservoir is during a CO2 injection.

To date, research has not focused on the behavior of concrete block stoppings subjected to excessive vertical loading due to roof to floor convergence. For this reason, the failure mechanism of stoppings under vertical loading has not been fully understood. Numerical models were used in combination with physical testing to study the failure mechanisms of concrete block stoppings. Initially, the behavior of a single standard CMU block was observed and simulated using FLAC. Full-scale stoppings were then tested in the Mine Roof Simulator and modeled using UDEC. Through a combination of physical testing and numerical modeling a failure mechanism for concrete block stoppings was established. This failure mechanism consists of development of stress concentrations where a height difference as small as 1/32â exists between adjacent blocks. These stress concentrations lead to tensile cracking and, ultimately, premature failure of the wall. / Master of Science

Modeling of the electromagnetic, structural, thermal, or acoustic response of the human body to various external and internal stimuli is limited by the availability of anatomically accurate and numerically efficient computational models. The models currently approved for use are generally of proprietary or fixed format, preventing new model construction or customization. 1. This dissertation develops a new Visible Human Project - Female (VHP-F) computational phantom, constructed via segmentation of anatomical cryosection images taken in the axial plane of the human body. Its unique property is superior resolution on human head. In its current form, the VHP-F model contains 33 separate objects describing a variety of human tissues within the head and torso. Each obejct is a non-intersecting 2-manifold model composed of contiguous surface triangular elements making the VHP-F model compatible with major commercial and academic numerical simulators employing the Finite Element Method (FEM), Boundary Element Method (BEM), Finite Volume Method (FVM), and Finite-Difference Time-Domain (FDTD) Method. 2. This dissertation develops a new workflow used to construct the VHP-F model that may be utilized to build accessible custom models from any medical image data source. The workflow is customizable and flexible, enabling the creation of standard and parametrically varying models facilitating research on impacts associated with fluctuation of body characteristics (for example, skin thickness) and dynamic processes such as fluid pulsation. 3. This dissertation identifies, enables, and quantifies three new specific computational bioelectromagnetic problems, each of which is solved with the help of the developed VHP-F model: I. Transcranial Direct Current Stimulation (tDCS) of human brain motor cortex with extracephalic versus cephalic electrodes; II. RF channel characterization within cerebral cortex with novel small on-body directional antennas; III. Body Area Network (BAN) characterization and RF localization within the human body using the FDTD method and small antenna models with coincident phase centers. Each of those problems has been (or will be) the subject of a separate dedicated MS thesis.