Air pollution regulations have driven modern power generation systems to move from diffusion to premixed combustion. However, these premixed combustion systems are prone to combustion instability, causing high fluctuations in pressure and temperature. This results in shortening of component life, system failure, or even catastrophic disasters. A large number of studies have been performed to understand and quantify the onset of combustion instability and the limit cycle amplitude. However, much work remains due to the complexity of the process associated with flow dynamics and chemistry. This thesis focuses on identifying, quantifying and predicting mechanisms of flame response subject to disturbances. A promising tool for predicting combustion instability is a flame transfer function. The flame transfer function is obtained by integrating unsteady heat release over the combustor domain. Thus, the better understanding of spatio-temporal characteristics of flame is required to better predict the flame transfer function. The spatio-temporal flame response is analyzed by the flame kinematic equation, so called G-equation. The flame is assumed to be a thin interface separating products and reactant, and the interface is governed by the local flow and the flame propagation. Much of the efforts were done to the flame response subject to the harmonic velocity disturbance. A key assumption allowing for analytic solutions is that the velocity is prescribed. For the mathematical tools, small perturbation theory, Hopf-Lax formula and numerical simulation were used. Solutions indicated that the flame response can be divided into three regions, referred to here as the near-field, mid-field, and farfield. In each regime, analytical expressions were derived, and those results were compared with numerical and experimental data. In the near field, it was shown that the flame response grows linearly with the normal component of the velocity disturbance. In the mid field, the flame response shows peaks in gain, and the axial location of these peaks can be predicted by the interference pattern by two characteristic waves. Lastly, in the far field where the flame response decreases, three mechanisms are studied; they are kinematic restoration, flame stretch, and turbulent flow effects. For each mechanism, key parameters are identified and their relative significances are compared.

The effect of freestream turbulence on generation of crossflow disturbances over swept wings is investigated through direct numerical simulations.  The set up follows  the  experiments  performed  by Downs  et  al.  in their  TAMU  experi- ment.  In this experiment the authors use ASU(67)-0315 wing geometry which promotes  growth  of crossflow  disturbances.   Distributed  roughness  elements are locally placed near the leading edge with a span-wise wavenumber, to ex- cite the corresponding crossflow vortices.  The response of boundary layer to external disturbances such as roughness heights, span-wise wavenumbers, Rey- nolds numbers and freestream turbulence characteristics are studied.  It must be noted that the experiments were conducted at a very low level of freestream turbulence  intensity  (T u).   In this  study,  we fully  reproduce the  freestream isotropic homogenous turbulence through a DNS code using detailed freestream spectrum data provided by the experiment. The generated freestream fields are then applied as the inflow boundary condition for direct numerical simulation of the wing. The geometrical set up is the same as the experiment along with application of distributed roughness elements near the leading edge to precipi- tate stationary crossflow disturbances.  The effects of the generated freestream turbulence are then studied on the initial amplitudes and growth of the bound- ary layer perturbations.  It appears that the freestream turbulence damps out the dominant stationary crossflow vortices. / <p>QC 20130604</p>

The current study addresses the stability of shock waves propagating through dissipative media, analogous to both granular media and molecular gases undergoing endothermic reactions. In order to investigate the stability, a simple molecular dynamics model was developed to observe shock waves and their structures with the inclusion of energy dissipation. For this, an Event Driven Molecular Dynamics model was implemented in a 2D environment, where a molecule is represented by a disk. The simulations addressed the formation of a shock wave in a gas by the sudden acceleration of a piston. Inelastic collisions were assumed to occur only if an impact velocity threshold is surpassed, representing the activation energy of the dissipative reactions. Parametric studies were conducted for this molecular model, by varying the strength of the shock wave, the activation threshold and the degree of inelasticity in the collisions. The resulting simulations showed that a shock structure does indeed become unstable with the presence of dissipative collisions. This instability manifests itself in the form of distinctive high density non-uniformities behind the shock wave, which take the form of convective rolls. The spacing and size of this ``finger-like" unstable pattern was shown to be dependent on the degree of inelasticity, the activation energy, and the strength of the driving piston. The mechanism responsible for the instability was addressed by studying the time evolution of the material undergoing the shock wave compression and further relaxation. It is found that the gas develops the instability on the same time scales as the clustering instability in homogeneous gases, first observed by Goldhirsch and Zanetti in granular gases. This confirmed that the clustering instability is the dominant mechanism.

The tumour suppressor p53-binding protein 1 (53BP1) is phosphorylated following DNA double strand breaks (DSBs); however, little is understood about the upstream signaling pathways that control this phosphorylation. Additionally, it is not known how these processes combine with 53BP1 to control the survival of cells following DNA damage such as that imparted by ionizing radiation (IR), which is the basis of radiotherapy. In this thesis, I have shown that 53BP1 is phosphorylated specifically in S-phase cells, but not relocalized to intranuclear foci, in response to severe oxygen stress. This occurs with only partial dependence on the ATM kinase (Chapter 2). Following IR, I find that both ATM and DNA-PKcs contribute to intranuclear phosphorylated 53BP1 foci, but that this phosphorylation is independent of proximal signaling molecules that control the localization of 53BP1 to initial DSBs (Chapter 3). Furthermore, I show that 53BP1 loss confers sensitivity to IR and this can be further augmented by inhibition of ATM and DNA-PKcs kinases suggesting that there are both 53BP1-dependent and -independent pathways of survival from IR (Chapter 4). These findings may have important implications for molecular pathology and personalized medicine as 53BP1 has recently been found to be activated or lost in subsets of human tumours. I have collaborated to initiate the development of a novel system to interrogate the implications of 53BP1 loss as traditional siRNA approaches in human cancer cells were not feasible (Chapter 5 and Appendix 2). This system can be used in vivo as tumour xenografts to further understand how 53BP1 and the tumour microenvironment interact endogenously and in response to IR. I also present the possibility and proof of concept for the use of 53BP1 as a biomarker in primary human prostate cancer tissue where little is known about 53BP1 biology (Chapter 5).

The research described in this thesis began with a single long-term objective: modelling of the vertebral artery during chiropractic manipulation of the cervical spine. Although chiropractic treatment has become prevalent, the possible correlation between neck manipulation and subsequent stroke in patients has been the subject of debate without resolution. Past research has been qualitative or statistical, whereas resolution demands a fundamental understanding of the associated mechanics.<p> Analysis in the thesis begins with a study of the anatomy and properties pertinent to the chiropractic problem. This indicates that the complexity of the problem will necessitate a long-term multidisciplinary effort including a nonlinear finite element formulation effective in analysing image data for soft tissue modelled as nearly incompressible. This leads to an assessment of existing finite element methods and the conclusion that new equation solving techniques are needed to ensure numerical stability.<p> Three techniques for effectively eliminating the source of numerical instability are developed and demonstrated with the aid of original finite element codes. Two of the methods are derived as modifications of matrix decomposition algorithms, while the third method constitutes a new finite element formulation. In addition, the understanding gained in developing these methods is used to produce a theorem for assessing a different but related problem: deformation of a nearly incompressible material subjected to a single concentrated force. Throughout the thesis, an interdisciplinary path from chiropractic problem to numerical algorithms is outlined, and results are in the form of mathematical proofs and derivations of both existing and new methods.

Down-hole damages such as borehole collapse, circulation loss and rock tensile/compressive cracking in the open-hole system are well understood at drilling and well completion stages. However, less effort has been made to understand the instability of cemented sections in High Pressure High Temperature (HPHT) wells. The existing analysis shows that, in the perforation zones, casing/cement is subject to instability, particularly in the presence of cavities. This dissertation focuses on the instability mechanism of casing/cement in the non-perforated zones. We investigate the transient thermal behavior in the casing-cement-formation system resulting from the movement of wellbore fluid using finite element method. The critical value of down-hole stresses is identified in both wellbore heating and cooling effects. Differently with the heating effect, the strong cooling effect in a cased hole can produce significant tension inside casing/cement. The confining formation has an obvious influence on the stability of casing/cement. The proposed results reveal that the casing/cement system in the non-homogeneous formation behaves differently from that in homogeneous formation. With this in mind, a three-dimensional layered finite element model is developed to illustrate the casing/cement mechanical behavior in the non-homogeneous formation. The radial stress of cement sheath is found to be highly variable and affected by the contrast in Young’s moduli in the different formation layers. The maximum stress is predicted to concentrate in the casing-cement system confined by the sandstone. Casing wear in the cased-hole system causes significant casing strength reduction, possibly resulting in the casing-cement tangential collapse. In this study, an approach for calculating the stress concentration in the worn casing with considering temperature change is developed, based on boundary superposition. The numerical results indicate that the casing-cement system after casing wear will suffer from severe tangential instability due to the elevated compressive hoop stress. Gas migration during the cementing process results from the fluid cement’s inability to balance formation pore pressure. Past experience emphasized the application of chemical additives to reduce or control gas migration during the cementing process. This report presents the thermal and mechanical behaviors in a cased hole caused by created gas channels after gas migration. In conclusion, the size and the number of gas channels are two important factors in determining mechanical instability in a casing-cement system.

No description available.

No description available.

No description available.

The Northern Current flows cyclonically contouring the continental slope in the NW Mediterranean. At the entrance of the Gulf of Lions this current is about 20 -- 30 km wide and flows along the deepest half of the continental slope, i.e. over the 1000 to 2000 m isobaths approximately. Surface speeds are of 30 -- 50 cm s^{-1}. In the MATER HFF experiment (March -- May 1997) mesoscale variability of the Northern Current is observed from current meter records, SST images and hydrographic data. The HFF experimental box is 20 x 40 km, covering the upper half of the slope (i.e. covering from 250 m to 1250 m depth isobaths). Current meter and satellite data show that the site is embedded in a region of significant Northern Current meandering and eddy activity. From SST images, meander wavelengths are estimated larger than 60 km, embracing smaller structures. These flow patterns affect upper-layer waters down to at least 650 m depth. Current meter data distinguish two narrow energetic bands centred at 3.5 days and 7.5 days, respectively, in agreement with previous studies.Baroclinic instability is viewed as a possible mechanism to explain the generation of the Northern Current meanders. The analytical model of Tang (1975) predicts the development of unstable waves of wavelength (> 60 km) and periods compatible with the 7.5 day band recorded with current meter devices. The higher frequency band of 3.5 days is out of the frequency range predicted by the classical baroclinic instability theory and it is discussed as a restriction of quasi-geostrophic theory.Barotropic instability is studied using a laboratory model of a -westward' jet flowing over the lower half of the continental slope, which considers dynamic similarity with the Northern Current. The laboratory model is cross-validated with a corresponding numerical model. Jet instabilities of currents similar to the Northern Current (i.e. westward jets) occur at the edges of the jet, showing a clear meandering tendency over the mid-slope. Westward currents of Ro = 0.1 -- 0.2 develop instabilities of wavelengths (50 -- 75 km) similar to those observed from SST images, with periods (3.3 -- 3.8 days) compatible with the 3.5 days period band recorded with HFFE current meters.The laboratory and numerical experiments have reproduced westward jets (as the Northern Current), but also eastward jets, in order to have a full approach to better understand the role of the bottom topography on barotropic instabilities. The slope current instabilities are successfully explained by the Marcus and Lee theory (1998) of jets on a beta plane. This theory is valid for westward flows with Ro > 0.1 and for eastward flows with Ro > 0.2 (jets of the so-called Regime II flows in this thesis), and it states that the instabilities of each shear layer of the barotropic jet take the appearance of a Kelvin-Helmholtz-like pattern, associated with a Rossby wave (of topographic origin in our case). According to this theory, the differences between eastward and westward jets rely on the disposition of the Rossby waves --at the centre of the current in eastward flows and at the edges of the jet in westward currents. Jets over a sloping bottom with small Rossby numbers (Ro < 0.1 for westward jets; Ro < 0.2 for eastward jets) show a flow pattern (the so-called Regime I in this thesis) that has common characteristics for eastward and westward flows. In these -small'-Ro flows, Kelvin-Helmholtz-like instabilities dominate, whereas Rossby waves are too weak to produce any major difference between jets flowing in eastward or westward direction. This occurs when the topographic influence, assumed proportional to the Ro number of the jet, is small.The differences between eastward and westward slope currents observed in this work (and similar observations of jets on a beta-plane from previous works) are explained in this thesis by a simple scheme based on conservation of potential vorticity, considering there are two main components in balance: the shear-induced vorticity and the topographically induced vorticity. The signs of these two components are determined by the relative direction of the flow with respect to the inclination of the bottom topography. Once the critical Rossby number is overpassed so that the topographic effects are important (Ro > 0.1 for westward jets; Ro > 0.2 for eastward jets), conservation of potential vorticity tends to enhance vortices at the centre of eastward jets --eastward jets show meandering at the jet core. In westward jets, potential vorticity conservation is responsible of enhancing vortices at each edge of the jet. Thus, westward jets (as the Northern Current) are broad and meandering occurs at the jet edges.In Ro > 0.1 westward flows (i.e. Regime II westward jets) a topographic Rossby wave appears over the shelf break. This result is likely observed because of the specific topography used in this work --a continental slope and a continental shelf separated by a shelf break, producing a strong change in ambient potential vorticity. Numerical simulations reveal that this Rossby wave is triggered by the slope current. This topographic Rossby wave is a robust pattern, since it is independent of the position of the current over the slope, the shape of the velocity shear profile of the jet, and the jet width. Although this type of wave could not be inferred from the HFFE field data, it could be a focus of study in further field experiments. It also needs further analytical consideration. The general conclusion extracted from this thesis that tries to explain the mesoscale variability associated to the Northern Current is that both baroclinic and barotropic instability could explain part of the oceanic observations. As a consequence, mixed barotropic-baroclinic instability (which occurs at wavelengths which are between those corresponding to pure barotropic and pure baroclinic instability) is thought to play an important role on the observed mesoscale variability. The resulting wavelength would depend on the relative strength of both mechanisms.