This thesis deals with the physical modeling of the initiation and propagation of upward positive leader discharges from grounded structures during lightning strikes. It includes the analysis of upward leaders initiated under the influence of the electric field produced by a dominant negative cloud charge and due to the combined action of a negative thundercloud and a descending downward stepped negative leader. Thus, a self-consistent model based on the physics of leader discharges is developed for the evaluation of the attachment of lightning flashes to any kind of grounded structure. The predictions of the model have been found to be in good agreement with the results of laboratory long air gap experiments and with classical and altitude rocket triggered lightning experiments. Due to the high application level and predictive power of the developed model, several contributions to the physical understanding of factors influencing the initiation and propagation of upward positive leaders during thunderstorms have been made. For instance, it has been found that the initiation of upward connecting leaders is strongly affected by the average velocity of the downward stepped leader. Similarly, it is shown that the switching voltage impulses used in the laboratory do not “fairly approximate” the electric fields produced by a descending downward leader, as claimed by supporters of Early Streamer Emission (ESE) devices. Furthermore, it is found that the space charge layer created by corona at ground level significantly increases the thundercloud electric fields required to initiate upward lightning leaders from tall objects. On the other hand, it is also shown that the upward leader velocity depends on the downward leader average velocity, the prospective return stroke current, the lateral distance of the downward leader channel and the ambient electric field. By implementing the model to the analysis of complex structures, it has been observed that the corners of actual buildings struck by lightning coincide rather well with the places characterized by low leader inception electric fields. Besides, it has been found that the leader inception zones of the corners of complex structures do not define symmetrical and circular regions as it is generally assumed.

The aim of this thesis work is to increase the knowledge of phenomena taking place during the initial stage in a top blown converter. The work has been done in a few steps resulting in four different supplements. Water model experiments have been carried out using particle image velocimetry (PIV) technology. The system investigated was a fundamental top blown converter where an air jet was set to impinge on a water surface. The flow field of the combined blown case, where an air jet was introduced through a bottom nozzle, was also captured by the PIV. The work clearly showed that the flow field caused by an impinging top blown jet alone could not match that of the bottom blown case. The main re-circulation loop (or vortex) was investigated with respect to position and it was found that an increased flow rate pushes the center of the re-circulation loop downwards into the bath. However, for the top-blown case there is a point when the flow rate is too large to cause a distinguishable re-circulation loop since the jet becomes more plunging (i.e. penetrates deep into the bath) than impinging, with large surface agitation and splashing as a result.A numerical model with the same dimensions as the experimental system was then created. Three different turbulence models from the same family were tested: standard-, realizable- and a modified-(slight modification of one of the coefficients in order to produce less spreading of the air jet) k-ε turbulence model. It could be shown that for the family of k-ε turbulence models the difference in penetration depth was small and that the values corresponded well to literature data. However, when it comes to the position of the re-circulation loop it was shown that the realizable k-ε model produced better results when comparing the results to the experimental data produced from the PIV measurements, mentioned earlier.It was then shown how the computational fluid dynamics (CFD) model could be coupled to thermodynamics databases in order to solve for both reactions and transport in the system. Instead of an air-water system, a gas-steel-slag system was created using the knowledge obtained in the previous simulation step described above. Reactions between gas-steel, gas-slag, steel-slag and gas-steel-slag were considered. Extrapolation of data from a few seconds of simulation was used for comparison to experimental data from the literature and showed reasonable agreement. The overall conclusion was that it is possible to make a coupling of the Thermo-Calc databases and a CFD software to make dynamic simulations of metallurgical processes such as a top-blown converter.A parametric study was then undertaken where two different steel grades were tested; one with high initial carbon content (3.85 mass-%) and one with lower carbon content (0.5 mass-%). The initial silicon content was held constant at 0.84 mass-%. Different initial temperatures were tested and also some variation in initial dissolved oxygen content was tried. It was found that the rate of decarburization/desiliconization was influenced by the temperature and carbon concentration in the melt, where a high temperature as well as a high carbon concentration favors decarburization over desiliconization. It was also seen that the region affected by a lower concentration of alloys (or impurities) was quite small close to the axis where the impinging jet hits the bath. Add the oscillating nature of the cavity and it was realized that sampling from this region during an experiment might be quite difficult. / QC 20100720

Numerical modeling of fluid flow and dissolved species transport in the subsurface is a challenging task, given variability and measurement uncertainty in the physical properties of the rock, the complexities of multi-fluid interaction, and limited computational resources. Nonetheless, this thesis seeks to expand our modeling capabilities in the context of contaminant hydrogeology. We describe the numerical simulator CompFlow Bio and use it to model invasion of a nonaqueous phase liquid (NAPL) contaminant through the vadose zone and below the water table in a fractured porous rock. CompFlow Bio is a three-phase, multicomponent, deterministic numerical model for fluid flow and dissolved species transport; it includes capillary pressure and equilibrium partitioning relationships. We have augmented the model to include randomly generated, axis-aligned, discrete fracture networks (DFNs). The DFN is coupled with the porous medium (PM) to form a single continuum. The domain is discretized using a finite-volume scheme in an unstructured mesh of rectilinear control volumes (CVs). Herein we present the governing equations, unstructured mesh creation scheme, algebraic development of fracture intersection CV elimination, and coupling of PM CVs over a fracture plane to permit asperity contact bridged flow. We include: small scale two-phase water-air and NAPL-water simulations to validate the practice of intersection CV elimination; small scale simulations with water-air, NAPL-water, and NAPL-water-air systems in a grid refinement exercise and to demonstrate the effect of asperity contact bridged flow; intermediate scale 3D simulations of NAPL invading the saturated zone, based on the Smithville, Ontario, site; intermediate scale 2D and 3D simulations of NAPL invading the vadose zone and saturated zone with transient recharge, based on the Santa Susana Field Laboratory site, California. Our findings indicate that: the formulation provides a practical and satisfactory way of modeling three-phase flow in discretely fractured porous rock; numerical error caused by spatial discretization manifests itself as several biases in physical flow processes; that asperity contact is important in establishing target water saturation conditions in the vadose zone; and simulation results are sensitive to relative permeability-saturation-capillary pressure relationships. We suggest a number of enhancements to CompFlow Bio to overcome certain computational limitations.

The Brunswick No. 6 area, which is located in the Bathurst Mining Camp, New Brunswick, Canada, is the focus of this thesis. Almost a decade ago, in order to improve the understanding of the crustal structures and explore for new mineral deposits at depth, three 2D seismic profiles totaling about 30 km and 3D seismic data covering an area of about 38 km2 were acquired from the study area. Petrophysical properties including compressional-wave velocity and density were also measured in two deep boreholes in the area. These data were recovered and reanalyzed, and the improved seismic images interpreted as the main part of this PhD thesis. A prestack DMO and poststack migration algorithm was considered for processing both 2D and 3D data. Processing of 2D data revealed shallow and deep reflections, which correlate well with surface geology. Steeply-dipping reflections, some of which could host mineral deposits, were imaged down to a depth of 6-7 km. Processing of 3D data showed similar results to the processed 2D profiles. Nevertheless, the non-orthogonal nature of the 3D survey, combined with irregular distribution of offsets, azimuths and trace midpoints, caused a severe acquisition footprint masking reflections in the DMO-corrected unmigrated stacked cube. An FK-dip filter in the wavenumber domain was designed to reduce the effects of the acquisition footprint. To better understand wave propagation and scattering effects, calculated acoustic impedance log from the available borehole data was used to estimate vertical scale length using a von Karman autocorrelation function. 2D synthetic models representative of heterogeneity in the area were generated accounting for the estimated scale length. Numerical modeling was used to study the scattering effects on the synthetic models, where some predefined targets were superimposed in the provided 2D heterogeneous medium. The effects of variable source frequency, longer horizontal scale length and petrophysical fluctuations of heterogeneous medium were also investigated. The modeling results indicate that, in the presence of large horizontal, but small vertical scale lengths (structural anisotropy), the identification of mineral deposits is possible in the unmigrated stacked sections, but can be challenging in the migrated sections.

The general objective of the research effort reported in this thesis was to develop the knowledge required to optimize the design and operation of solid and semi-solid manure handling and land application equipment. Selected physical and rheological properties of manure products deemed to have an influence on the performances of manure handling and land application equipment were measured and general trends were identified among the measured properties. Relationships were also established between the measured properties and the type of manure as well as its total solids concentration. Field experiments were carried out to evaluate the effects of selected mechanical configurations, operating parameters and product properties on the discharge of manure spreaders. The influence of the type of conveying system (scraper conveyor and system of four augers) and the velocity at which it is operated, the geometry of the holding system and the position of a flow-control gate were all included in the analysis. The discharge rates of the machines as well as the specific energy required by the unloading operations were measured. A numerical modeling method called discrete element method (DEM) was used to create virtual manure, a numerical model of the real product. The measured physical and flow properties were used to develop and validate the virtual manure models. It was found that manure products could successfully be represented in a DE framework and that several parameters defining the contact constitutive model in the DEM had an influence on the behaviour of the virtual products. The DEM was then used to study machine-product interactions taking place in handling and land application equipment. Results from field experiments carried out using various land application equipment were used in the development and validation of the interaction models. The predicted flow rates and power requirements were in good agreement with measured data. The results obtained allowed for a better understanding of the flow of manure products in manure handling and land application equipment. It was found that the constitutive model used for the product influenced the results of the machine-product interactions models. A precision banded applicator under development at the University of Saskatchewan was also modeled. The discharge rate of this equipment is influenced by a number of parameters. The predicted mass distribution across the width of the banded applicator was well correlated to the experimental results. The models developed in this thesis have the potential to become powerful engineering tools for the design of improved machines for the handling and land application of solid and semi-solid manure.

In this thesis, studies associated with the stability analyses of the box-shaped dump-site culvert constructed in Tinaz Surface Mine of Turkish Coal Enterprises (TKi) are presented. In addition, stability conditions of other culvert alternatives are evaluated. Existence of creeks in a surface mining area is a significant factor to be considered in selection of dump-site location. Since, the dumped overburden material on the valley acts as a barrier and behaves like a dam causing flood problem behind the dump-site. TKi engineers prevented the flood potential that might have occurred behind the dump-site by constructing a 480-meter long reinforced-concrete culvert on the downstream of Gevenez Creek Valley. However, considerable amount of deformations occurred in the first 100 meters of the culvert, as a result of overburden material being replaced on this structure. In order to determine the failure mechanism associated with the culvert, a series of numerical modeling analyses were carried out utilizing back analysis technique. The validity of the numerical model was justified by convergence measurements and observations carried out inside the culvert as overburden material being replaced on the stable part of this structure. Finally, based on the numerical model developed, the stability of other culvert alternatives that could be used in future projects were evaluated considering different embankment conditions (positive projecting and negative projecting), bedding conditions (impermissible, ordinary, first-class and concrete cradle), culvert shapes (box and circular) and dumping conditions.

As explosive blasts continue to cause casualties in both civil and military environments, there is a need for increased understanding of the mechanisms of blast trauma at the organ level and a need for a more detailed predictive methodology. A fundamental understanding of blast injury will lead to the development of improved protective equipment and ultimately reduce the severity of injury. Models capable of predicting injury to varied blast loading will also reduce the emphasis on animal blast testing. To provide some historical context, this research was begun shortly after the U.S. led invasion of Iraq, and came to a close while there continues to be daily loss of life from blast injuries in the Middle East, as well as continued threats of terrorism throughout the world. In addition to industrial accidents, it is clear that blast injury is far more than just a military concern. Simplified finite element models of the human and sheep thoraces were created in order to provide practical and flexible models for the prediction of primary blast injury in simple and complex blast environments, and subsequently for the development of improved protective equipment. The models were created based on actual human and sheep geometries and published material properties. The fluid-structure interaction of the models compared well with experimental blast studies carried out during the course of the research, as shown by comparing actual and predicted overpressures in the free field and at the thorax. By comparing the models to published experimental data from simple blasts, trends in the results were verified and peak lung pressure was proposed as a trauma criterion. Local extent of injury in the lung is correlated to the peak pressure measured in each finite element, categorized as no injury (< 60 kPa), trace (60-100 kPa), slight (100-140 kPa), moderate (140-240 kPa) and severe (> 240 kPa). The calculation of the mean value of the peak lung pressures of all of the finite elements allows for an overall estimate of the injury level, with 35 kPa predicting threshold damage, 129 kPa for one percent lethality, and 186 kPa for fifty percent lethality. The simple blast results also compared well to the predictions of two previously validated mathematical models. Variation of predicted injury within a given loading severity was 15%, which is comparable to the model by Stuhmiller that had a variation of 20%. The model by Axelsson had very little variation (1.4%), but the differences between levels of severity were quite small, and often difficult to decipher. In addition to predicting consistent levels of injury, the finite element models were able to provide insight into the trauma mechanism, map the extent of injury through the lungs, and validate a local injury criterion. The models were then applied to predict injury under complex blast loading by subjecting the human finite element torso to a threshold level blast while located at varying distances from a wall or a corner. The results compared well to the validated mathematical models, showing a sharp increase in injury severity as the model approached the reflecting surface. When directly against the wall, the mean of the peak lung pressure values was 57 kPa, and in the corner, the mean value reached 69 kPa. Although these values did not reach the level representing one percent lethality, they do represent a significant increase in injury above threshold as a direct result of the surrounding geometry. Once again, the finite element models correctly showed injury trends and lung injury patterns reported in experiments. The models predicted the level of injury and were able to predict the time varying pattern of injury, which is something existing models cannot do. Having designed the models from physical principals, and having validated the models against published results, they can now be used in the continued development of protective equipment. Acknowledging that this model was the first iteration, the author believes that improvements in material properties, mesh refinement, and the investigation of other possible parameters for the prediction of injury will lead to substantial advances in the understanding of primary blast injury.

Basic oxygen furnace (BOF) slag media were studied through a series of laboratory, modeling and field studies as a potential treatment material for use in on site wastewater disposal systems. Microsphere enumeration methodology was examined in a factorial experiment to evaluate the minimum density and minimum number of microspheres that should be counted to ensure accurate and precise estimations of concentration. The results suggest that to minimize variability at least 350 microspheres should be counted and a microsphere density of 25-40 microspheres field-1 is necessary. A review of existing methodologies for high-titer bacteriophage production was conducted and an amalgamation of existing methodologies was chosen that reliably achieves elevated concentration and ensures a purified suspension. A combination of batch and column studies was conducted to evaluate the removal of the bacteriophage, PRD-1, and virus-sized fluorescent microspheres by BOF media, and to delineate the relative contributions of the two principle attenuation processes, inactivation and attachment. In the batch studies, substantial removal of PRD-1 does not occur in the pH 7.6 and 9.5 suspensions, but at pH 11.4, removal of the virus was 2.1 log C/C0 day-1 for the first two days, followed by 0.124 log C/C0 day-1 over the subsequent 10 days. Two column studies were conducted after 60 and 300 days of saturation with artificial groundwater at a flow rate of 1 pore volume day-1 using two BOF mixtures. After 300 days of column saturation, microsphere concentrations approached input levels, indicating a removal of 0.1-0.2 log C/C0 and suggesting attachment processes were negligible. PRD-1 removal was more pronounced (1.0-1.5 log C/C0). The reduction of PRD-1 is likely the result of a combination of virus inactivation at elevated pH (10.6-11.4), and attachment processes. Geochemical factors controlling microsphere attachment were compared between the two sets of experiments after 60 and 300 days of column saturation. Differences in attachment efficiency may be due to higher influent DOC concentration in the second experiment, conversion of amorphous iron phases to more crystalline forms over time, reductive dissolution of preferable attachment sites on iron phases, or precipitation of calcite. Hydrus-1D, a one-dimensional numerical model, was used to quantify transport processes, inactivation and attachment/detachment, occurring in the column experiments by model inversion. Fitted microsphere breakthrough closely matched observed data, whereas PRD-1 breakthrough with realistic parameter values does not closely match the peaked nature of the observed curves. The model achieved improved fits for microsphere and PRD-1 breakthrough when both strongly- and weakly-binding sites are represented. A unique set of parameter estimates could not be determined because of overparameterization of the inverse modeling for the experimental systems. An alternative latrine incorporating BOF slag media was constructed in a periurban community located near São Paulo, Brazil. Pathogen indicator removal is approximately 4-5 orders of magnitude in less than one meter of vertical transport through the BOF slag media. In a control latrine, constructed with similar hydraulic characteristics and inert materials, comparable reductions in pathogenic indicators were observed over three meters of vertical transport.

It has long been known that any activity that results in changes in subsurface pressure, such as hydrocarbon production or waste or water reinjection, also causes underground deformations and movement, which can be described in terms of volumetric changes. Such deformations induce surface movement, which has a significant environmental impact. Induced surface deformations are measurable as vertical displacements; horizontal displacements; and tilts, which are the gradient of the surface deformation. The initial component of this study is a numerical model developed in C++ to predict and calculate surface deformations based on assumed subsurface volumetric changes occurring in a reservoir. The model is based on the unidirectional expansion technique using equations from Okada’s theory of dislocations (Okada, 1985). A second numerical model calculates subsurface volumetric changes based on surface deformation measurements, commonly referred to as solving for the inverse case. The inverse case is an ill-posed problem because the input is comprised of measured values that contain error. A regularization technique was therefore developed to help solve the ill-posed problem. A variety of surface deformation data sets were analyzed in order to determine the surface deformation input data that would produce the best solution and the optimum reconstruction of the initial subsurface volumetric changes. Tilt measurements, although very small, were found to be much better input than vertical displacement data for finding the inverse solution. Even in an ideal case with 0 % error, tilts result in a smaller RMSE (about 12 % smaller in the case studied) and thus a better resolution. In realistic cases with error, adding only 0.55 % of the maximum random error in the surface displacement data affects the back-calculated results to a significant extent: the RMSE increased by more than 13 times in the case studied. However, in an identical case using tilt measurements as input, adding 20 % of the maximum surface tilt value as random error increased the RMSE by 7 times, and remodelling the initial distribution of the volumetric changes in the subsurface was still possible. The required area of observation can also be reduced if tilt measurements are used. The optimal input includes tilt measurements in both directions: dz/dx and dz/dy. iv With respect to the number of observation points chosen, when tilts are used with an error of 0 %, very good resolution is obtainable using only 0.4 % of the unknowns as the number of benchmarks. For example, using only 10 observation points for a reservoir with 2500 elements, or unknowns resulted in an acceptable reconstruction. With respect to the sensitivity of the inverse solution to the depth of the reservoir and to the geometry of the observation grid, the deeper the reservoir, the more ill-posed the problem. The geometry of the benchmarks also has a significant effect on the solution of the inverse problem.

The initiation of subduction requires the formation of lithospheric plates which mostly deform at their edges. Shear heating is a possible candidate for producing such localized deformation. In this thesis we employ a 2D model of the mantle with a visco-elasto-plastic rheology and enabled shear heating. We are able to create a shear heating instability both in a constant strain rate and a constant stress boundary condition setup. For the rst case, localized deformation in our specic setup is found for strain rates of 10-15 1/s and mantle temperatures of 1300°C. For constant stress boundaries, the conditions for a setup to localize are more restrictive. Mantle motion is induced by large cold and hot temperature perturbations. Lithospheric stresses scale with the size of these perturbations; maximum stresses are on the order of the yield stress (1 GPa). Adding topography or large inhomogeneities does not result in lithospheric-scale fracture in our model. However, localized deformation does occur for a restricted parameter choice presented in this thesis. The perturbation size has little effect on the occurrence of localization, but large perturbations shorten its onset time.