Dynamic rupture simulation integrated with earthquake observations

Dynamic rupture simulations are unique in their contributions to the study of earthquake physics. The current rapid development of dynamic rupture simulations poses several new questions: Do the simulations reflect the real world? Do the simulations have predictive power? Which one should we believe when the simulations disagree? This thesis illustrates how integration with observations can help address these questions and reduce the effects of non-uniqueness of both dynamic rupture simulations and kinematic inversion problems. Dynamic rupture simulations with observational constraints can effectively identify non-physical features inferred from observations. Moreover, the integrative technique can also provide more physical insights into the mechanisms of earthquakes. This thesis demonstrates two examples of such kinds of integration: dynamic rupture simulations of the M_w 9.0 2011 Tohoku-Oki earthquake and of earthquake ruptures in damaged fault zones:

(1) We develop simulations of the Tohoku-Oki earthquake based on a variety of observations and minimum assumptions of model parameters. The simulations provide realistic estimations of stress drop and fracture energy of the region and explain the physical mechanisms of high-frequency radiation in the deep region. We also find that the overridding subduction wedge contributes significantly to the up-dip rupture propagation and large final slip in the shallow region. Such findings are also applicable to other megathrust earthquakes.

(2) Damaged fault zones are usually found around natural faults, but their effects on earthquake ruptures have been largely unknown. We simulate earthquake ruptures in damaged fault zones with material properties constrained by seismic and geological observations. We show that reflected waves in fault zones are effective at generating pulse-like ruptures and head waves tend to accelerate and decelerate rupture speeds. These mechanisms are robust in natural fault zones with large attenuation and off-fault plasticity. Moreover, earthquakes in damaged fault zones can propagate at super-Rayleigh speeds that are unstable in homogeneous media. Supershear transitions in fault zones do not require large fault stresses. In the end, we present observations in the Big Bear region, where variability of rupture speeds of small earthquakes correlates with the laterally variable materials in a damaged fault zone.

Ferromagnetic Resonance Force Microscopy

A new approach to magnetic resonance was introduced in 1992 based upon detection of spin-induced forces by J. Sidles [1]. This technique, now called magnetic resonance force microscopy (MRFM), was first demonstrated that same year via electron paramagnetic resonance (EPR) by D. Rugar et al. [2]. This new method combines principles of magnetic resonance with those of scanned probe technology to detect spin resonance through mechanical, rather than inductive, means. In this thesis the development and use of ferromagnetic resonance force microscopy (FMRFM) is described. This variant of MRFM, which allows investigation of ferromagnetic samples, was first demonstrated in 1996 by Z. Zhang et al. [3]. FMRFM enables characterization of (a) the dynamic magnetic properties of microscale magnetic devices, and (b) the spatial dependence of ferromagnetic resonance within a sample. Both are impossible with conventional ferromagnetic resonance techniques.

Ferromagnetically coupled systems, however, pose unique challenges for force detection. In this thesis the attainable spatial resolution - and the underlying physical mechanisms that determine it - are established. We analyze the dependence of the magnetostatic modes upon sample dimensions using a series of microscale yttrium iron garnet (YIG) samples. Mapping of mode amplitudes within these sample is attained with an unprecedented spatial resolution of 15μm. The modes, never before analyzed on this scale, fit simple models developed in this thesis for samples of micron dimensions. The application of stronger gradient fields induces localized perturbation of the ferromagnetic resonance modes. The first demonstrations of this effect are presented in this study, and a simple theoretical model is developed to explain our observations. The results indicate that the characteristics of the locally-detected ferromagnetic modes are still largely determined by the external fields and dimensions of the entire sample, rather than by the localized interaction volume (i.e., the locale most strongly affected by the local gradient field). Establishing this is a crucial first step toward understanding FMRFM in the high gradient field limit where the dispersion relations become locally determined. In this high gradient field regime, FMRFM imaging becomes analogous with that of EPR MRFM.

FMRFM has also been employed to characterize magnetic multilayers, similar to those utilized in giant magnetoresistance (GMR) devices, on a lateral scale 40 x 40μm. This is orders of magnitude smaller than possible via conventional methods. Anisotropy energies, thickness, and interface qualities of individual layers have been resolved.

This initial work clearly demonstrates the immense and unique potential that FMRFM offers for characterizing advanced magnetic nanostructures and magnetic devices.

From molecules to organs : Microscopy and multi-scale nature of development

Morphogenesis is a phenomenon of intricate balance and dynamic interplay between processes occurring at a wide range of scales (spatial, temporal and energetic). During development, a variety of physical mechanisms are employed by tissues to simultaneously pattern, move, and differentiate based on information exchange between constituent cells, perhaps more than at any other time during an organism's life. To fully understand such events, a combined theoretical and experimental framework is required to assist in deciphering the correlations at both structural and functional levels at scales that include the intracellular and tissue levels as well as organs and organ systems. Microscopy, especially diffraction-limited light microscopy, has emerged as a central tool to capture the spatio-temporal context of life processes. Imaging has the unique advantage of watching biological events as they unfold over time at single-cell resolution in the intact animal. In this work I present a range of problems in morphogenesis, each unique in its requirements for novel quantitative imaging both in terms of the technique and analysis. Understanding the molecular basis for a developmental process involves investigating how genes and their products- mRNA and proteins-function in the context of a cell. Structural information holds the key to insights into mechanisms and imaging fixed specimens paves the first step towards deciphering gene function. The work presented in this thesis starts with the demonstration that the fluorescent signal from the challenging environment of whole-mount imaging, obtained by in situ hybridization chain reaction (HCR), scales linearly with the number of copies of target mRNA to provide quantitative sub-cellular mapping of mRNA expression within intact vertebrate embryos. The work then progresses to address aspects of imaging live embryonic development in a number of species. While processes such as avian cartilage growth require high spatial resolution and lower time resolution, dynamic events during zebrafish somitogenesis require higher time resolution to capture the protein localization as the somites mature. The requirements on imaging are even more stringent in case of the embryonic zebrafish heart that beats with a frequency of ~ 2-2.5 Hz, thereby requiring very fast imaging techniques based on two-photon light sheet microscope to capture its dynamics. In each of the hitherto-mentioned cases, ranging from the level of molecules to organs, an imaging framework is developed, both in terms of technique and analysis to allow quantitative assessment of the process in vivo. Overall the work presented in this thesis combines new quantitative tools with novel microscopy for the precise understanding of processes in embryonic development.

Wideband two-port beam splitter of a binary fused-silica phase grating

The usual beam splitter of multilayer-coated film with a wideband spectrum is not easy to achieve. We describe the realization of a wideband transmission two-port beam splitter based on a binary fused-silica phase grating. To achieve high efficiency and equality in the diffracted 0th and -1st orders, the grating profile parameters are optimized using rigorous coupled-wave analysis at a wavelength of 1550 nm. Holographic recording and the inductively coupled plasma dry etching technique are used to fabricate the fused-silica beam splitter grating. The measured efficiency of (45% x 2) = 90% diffracted into the both orders can be obtained with the fabricated grating under Littrow mounting. The physical mechanism of such a wideband two-port beam splitter grating can be well explained by the modal method based on two-beam interference of the modes excited by the incident wave. With the high damage threshold, low coefficient of thermal expansion, and wideband high efficiency, the presented beam splitter etched in fused silica should be a useful optical element for a variety of practical applications. (C) 2008 Optical Society of America.

Evidence from the Pacific Northwest for solar modulation of climate and of northern ecosystems

EXTRACT (SEE PDF FOR FULL ABSTRACT): The 250-year net annual snow accumulation, or mass balance, time series derived from the Mt. Logan (Yukon) ice core has been spectrally analyzed and is found to contain a nominal 11-year waveform. The stable isotope time series contains a significant amount of power between 9 and 13 years, although this record is evidently not a straightforward proxy for air temperatures. The signal in the mass balance time series exhibits a close relationship with the sunspot cycle waveform and is, therefore, assumed to be related to it. Waveforms showing a high correlation with the solar cycle are found in other climate data in the region. ... Taken collectively, the data point to a link between solar variability, atmospheric variability, climate, and selected ecological dynamics in the Pacific Northwest, but other data, not presented, indicate these relationships may hold elsewhere. So far, the evidence is empirical; complete details of the physical mechanisms involved have yet to be synthesized in a satisfactory way.

Numerical modelling of dynamic response of loss modulated monolithic multisection mode-locked semiconductor lasers

Monolithic multisection mode-locked semiconductor lasers with an integrated distributed Bragg reflector (DBR) have recently been demonstrated to generate stable picosecond pulses at high repetition rates suitable for optical communication systems. However, there has been very little theoretical work on understanding the physical mechanisms of the device and on optimisation of the absorber modulator design. This article presents numerical modeling of the loss modulated mode-locking process in these lasers. The model predicts most aspects experimentally observed within this type of device, and the results show the output waveform, optical spectrum, instantaneous frequency chirp, and stable operating range.

Applications of the dynamic mode decomposition

The decomposition of experimental data into dynamic modes using a data-based algorithm is applied to Schlieren snapshots of a helium jet and to time-resolved PIV-measurements of an unforced and harmonically forced jet. The algorithm relies on the reconstruction of a low-dimensional inter-snapshot map from the available flow field data. The spectral decomposition of this map results in an eigenvalue and eigenvector representation (referred to as dynamic modes) of the underlying fluid behavior contained in the processed flow fields. This dynamic mode decomposition allows the breakdown of a fluid process into dynamically revelant and coherent structures and thus aids in the characterization and quantification of physical mechanisms in fluid flow. © 2010 Springer-Verlag.

Real-time experimental demonstration of low-cost VCSEL intensity-modulated 11.25 Gb/s optical OFDM signal transmission over 25 km PON systems.

The feasibility of utilising low-cost, un-cooled vertical cavity surface-emitting lasers (VCSELs) as intensity modulators in real-time optical OFDM (OOFDM) transceivers is experimentally explored, for the first time, in terms of achievable signal bit rates, physical mechanisms limiting the transceiver performance and performance robustness. End-to-end real-time transmission of 11.25 Gb/s 64-QAM-encoded OOFDM signals over simple intensity modulation and direct detection, 25 km SSMF PON systems is experimentally demonstrated with a power penalty of 0.5 dB. The low extinction ratio of the VCSEL intensity-modulated OOFDM signal is identified to be the dominant factor determining the maximum obtainable transmission performance. Experimental investigations indicate that, in addition to the enhanced transceiver performance, adaptive power loading can also significantly improve the system performance robustness to variations in VCSEL operating conditions. As a direct result, the aforementioned capacity versus reach performance is still retained over a wide VCSEL bias (driving) current (voltage) range of 4.5 mA to 9 mA (275 mVpp to 320 mVpp). This work is of great value as it demonstrates the possibility of future mass production of cost-effective OOFDM transceivers for PON applications.

Cinematographic OH-PLIF measurements of two interacting turbulent premixed flames with and without acoustic forcing

This paper describes an experimental investigation into the interactions that occur between two lean turbulent premixed flames stabilised on conical bluff-bodies when they are moved closer together. Cinematographic OH-PLIF measurements were acquired to investigate adjacent flame front interactions as a function of flame separation distance (S). Flame surface density (FSD) and curvature were determined to characterise the unforced flames. Acoustic forcing was then applied to explore the amplitude dependent thermo-acoustic response. Phase-averaged FSD and global heat release measurements in the form of OH * chemiluminescence were obtained for a range of forcing frequencies (f) and amplitudes (A) as a function of S. As the flames were brought closer together the adjacent annular jets were found to merge into a single jet structure. This caused adjacent flame fronts to merge above the wake region between the two flames at a location determined by the jet efflux (flame angle) and S. This region of flame-flame interaction we refer to as 'interacting region'. In the unforced flames, a trend of increasingly negative curvature for decreasing S produced a small net increase in flame surface area via cusp formation. When subjected to acoustic forcing, S-dependent regimes were found in the global heat release response as a function A. The overall trend showed that the occurrence of jet/flame merging reduces the value of A at which non-linear response occurs. In support of previous findings for flames stabilised along shear layers, the phase-averaged FSD showed that the flame dynamics that drive the thermo-acoustic response result from the roll-up of vortices which generate large-scale vortex-flame interactions. Compared with axisymmetric flames, the occurrence of jet merging alters the vortex-flame interactions resulting in an asymmetric contribution to the heat release between the wall and interacting regions. The majority of the heat release was found to occur in the interacting region through the rapid production and destruction of flame surface area. The occurrence of jet merging and large-scale interactions between adjacent flames result in different physical mechanisms that drive the thermo-acoustic response compared with single axisymmetric flames. © 2011.

Nonlinear optics in Silicon photonic crystal cavities

Silicon is known to be a very good material for the realization of high-Q, low-volume photonic cavities, but at the same it is usually considered as a poor material for nonlinear optical functionalities like second-harmonic generation, because its second-order nonlinear susceptibility vanishes in the dipole approximation. In this work we demonstrate that nonlinear optical effects in silicon nanocavities can be strongly enhanced and even become macroscopically observable. We employ photonic crystal nanocavities in silicon membranes that are optimized simultaneously for high quality factor and efficient coupling to an incoming beam in the far field. Using a low-power, continuous-wave laser at telecommunication wavelengths as a pump beam, we demonstrate simultaneous generation of second- and third harmonics in the visible region, which can be observed with a simple camera. The results are in good agreement with a theoretical model that treats third-harmonic generation as a bulk effect in the cavity region, and second-harmonic generation as a surface effect arising from the vertical hole sidewalls. Optical bistability is also observed in the silicon nanocavities and its physical mechanisms (optical, due to two-photon generation of free carriers, as well as thermal) are investigated. © 2011 IEEE.

Complex chemistry DNS of n-heptane spray autoignition at high pressure and intermediate temperature conditions

Direct Numerical Simulations (DNS) of turbulent n-heptane sprays autoigniting at high pressure (P=24bar) and intermediate air temperature (Tair=1000K) have been performed to investigate the physical mechanisms present under conditions where low-temperature chemistry is expected to be important. The initial turbulence in the carrier gas, the global equivalence ratio in the spray region, and the initial droplet size distribution of the spray were varied. Results show that spray ignition exhibits a spotty nature, with several kernels developing independently in those regions where the mixture fraction is close to its most reactive value ξMR (as determined from homogeneous reactor calculations) and the scalar dissipation rate is low. Turbulence reduces the ignition delay time as it promotes mixing between air and the fuel vapor, eventually resulting in lower values of scalar dissipation. High values of the global equivalence ratio are responsible for a larger number of ignition kernels, due to the higher probability of finding regions where ξ=ξMR. Spray polydispersity results in the occurrence of ignition over a wider range of mixture fraction values. This is a consequence of the inhomogeneities in the mixing field that characterize these sprays, where poorly mixed rich spots are seen to alternate with leaner ones which are well-mixed. The DNS simulations presented in this work have also been used to assess the applicability of the Conditional Moment Closure (CMC) method to the simulation of spray combustion. CMC is found to be a valid method for capturing spray autoignition, although care should be taken in the modelling of the unclosed terms appearing in the CMC equations. © 2013 The Combustion Institute.

Direct comparison of the gravimetric responsivities of ZnO-based FBARs and SMRs

Film bulk acoustic resonators (FBARs) and solidly mounted resonators (SMRs) have the potential to significantly improve upon the sensitivity and minimum detection limit of traditional gravimetric sensors based on quartz crystal microbalances (QCMs) and surface acoustic wave resonators (SAWs). To date, neither FBAR nor SMR devices have been demonstrated to be superior to the other; hence the choice between them depends primarily on the users' ability to design/fabricate membranes and/or Bragg reflectors. In this work, it is shown that identically designed FBAR and SMR devices resonating at the same frequency exhibit different responsivities to mass loadings, Rm, and that the SMRs are less responsive than the FBARs. For the specific device design and resonant frequency (~2 GHz) of the resonators presented here, the FBARs' mass responsivity is ~20% greater than that of the SMRs', and although this value is not universal for all possible device designs, it clearly shows that FBAR devices should be favoured over SMRs in gravimetric sensing applications where the FBARs' fragility is not an issue. Numerical calculations based on Mason's model offer an insight into the physical mechanisms behind the greater FBARs responsivity, and it was shown that the Bragg reflector has an effect on the acoustic load at one of the facets of the piezoelectric films which is in turn responsible for the SMRs' lower responsivity to mass loadings. © 2013 Elsevier B.V.

The interaction of riblets with wall-bounded turbulence

The purpose of this thesis is to give answer to the question: why do riblets stop working for a certain size? Riblets are small surface grooves aligned in the mean direction of an overlying turbulent flow, designed specifically to reduce the friction between the flow and the surface. They were inspired by biological surfaces, like the oriented denticles in the skin of fastswimming sharks, and were the focus of a significant amount of research in the late eighties and nineties. Although it was found that the drag reduction depends on the riblet size scaled in wall units, the physical mechanisms implicated have not been completely understood up to now. It has been explained how riblets of vanishing size interact with the turbulent flow, producing a change in the drag proportional to their size, but that is not the regime of practical interest. The optimum performance is achieved for larger sizes, once that linear behavior has broken down, but before riblets begin adopting the character of regular roughness and increasing drag. This regime, which is the most relevant from a technological perspective, was precisely the less understood, so we have focused on it. Our efforts have followed three basic directions. First, we have re-assessed the available experimental data, seeking to identify common characteristics in the optimum regime across the different existing riblet geometries. This study has led to the proposal of a new length scale, the square root of the groove crosssection, to substitute the traditional peak-to-peak spacing. Scaling the riblet dimension with this length, the size of breakdown of the linear behavior becomes roughly universal. This suggests that the onset of the breakdown is related to a certain, fixed value of the cross-section of the groove. Second, we have conducted a set of direct numerical simulations of the turbulent flow over riblets, for sizes spanning the full drag reduction range. We have thus been able to reproduce the gradual transition between the different regimes. The spectral analysis of the flows has proven particularly fruitful, since it has made possible to identify spanwise rollers immediately above the riblets, which begin to appear when the riblet size is close to the optimum. This is a quite surprising feature of the flow, not because of the uniqueness of the phenomenon, which had been reported before for other types of complex and porous surfaces, but because most previous studies had focused on the detail of the flow above each riblet as a unit. Our novel approach has provided the adequate tools to capture coherent structures with an extended spanwise support, which interact with the riblets not individually, but collectively. We have also proven that those spanwise structures are responsible for the increase in drag past the viscous breakdown. Finally, we have analyzed the stability of the flow with a simplified model that connects the appearance of rollers to a Kelvin–Helmholtz-like instability, as is the case also for the flow over plant canopies and porous surfaces. In spite of the model emulating the presence of riblets only in an averaged, general fashion, it succeeds to capture the essential attributes of the breakdown, and provides a theoretical justification for the scaling with the groove cross-section.

Sensors based on SAW and FBAR technologies

Over the last few years a number of sensing platforms are being investigated for their use in drug development, microanalysis or medical diagnosis. Lab-on-a-chip (LOC) are devices integrating more than one laboratory functions on a single device chip of a very small size, and typically consist of two main components: microfluidic handling systems and sensors. The physical mechanisms that are generally used for microfluidics and sensors are different, hence making the integration of these components difficult and costly. In this work we present a lab-on-a-chip system based on surface acoustic waves (for fluid manipulation) and film bulk acoustic resonators (for sensing). Coupling surface acoustic waves into liquids induces acoustic streaming and motion of micro-droplets, whilst it is well-known that bulk acoustic waves can be used to fabricate microgravimetric sensors. Both technologies offer exceptional sensitivity and can be fabricated from piezoelectric thin films deposited on Si substrates, reducing the fabrication time/cost of the LOC devices. © 2013 SPIE.

Effects of indium ion implantation on the total dose hardness of fully depleted SOI NMOSFET

In this work, we investigate the effects of the indium ion implantation towards the back-channel interface on the total dose hardness of the n-channel SOI MOSFET. The results show that the indium implant has slight impact on the normal threshold voltage while preserving low leakage current after irradiation. The advantage is attributed to the narrow as-implanted and postanneal profile of the indium implantation. Two-dimensional simulations have been used to understand the physical mechanisms of the effects.