Tsunamis impose significant threat to human life and coastal infrastructure. The goal of my dissertation is to develop a robust, accurate, and computationally efficient numerical model for quantitative hazard assessment of tsunamis. The length scale of the physical domain of interest ranges from hundreds of kilometers, in the case of landslide-generated tsunamis, to thousands of kilometers, in the case of far-field tsunamis, while the water depth varies from couple of kilometers, in deep ocean, to few centimeters, in the vicinity of shoreline. The large multi-scale computational domain leads to challenging and expensive numerical simulations. I present and compare the numerical results for different important problems --- such as tsunami hazard mitigation due to presence of coastal vegetation, boulder dislodgement and displacement by long waves, and tsunamis generated by an asteroid impact --- in risk assessment of tsunamis. I employ depth-integrated shallow water equations and Serre-Green-Naghdi equations for solving the problems and compare them to available three-dimensional results obtained by mesh-free smoothed particle hydrodynamics and volume of fluid methods. My results suggest that depth-integrated equations, given the current hardware computational capacities and the large scales of the problems in hand, can produce results as accurate as three-dimensional schemes while being computationally more efficient by at least an order of a magnitude. / Ph. D.

Left ventricular diastolic dysfunction (LVDD) is any abnormality in the filling of the left ventricle (LV).  Despite the prevalence of this disease, it remains difficult to diagnose, mainly due to inherent compensatory mechanisms and a limited physical understanding of the filling process.  LV filling can be non-invasively imaged using color m-mode echocardiography which provides a spatio-temporal map of inflow velocity.  These filling patterns, or waves, are conventionally used to qualitatively assess the filling pattern, however, this work aims to physically quantify the filling waves to improve understanding of diastole and develop robust, reliable, and quantitative parameters. This work reveals that LV filling waves in a normal ventricle act as dispersive waves and not only propagate along the length of the LV but also spread and disperse in the direction of the apex.  In certain diseased ventricles, this dispersion is limited due to changes in LV geometry and wall motion.  This improved understanding could aid LVDD diagnostics not only for determining health and disease, but also for distinguishing between progressing disease states. This work also identifies a limitation in a current LVDD parameter, intra ventricular pressure difference (IVPD), and presents a new methodology to address this limitation.  This methodology is also capable of synthesizing velocity information from a series of heartbeats to generating one representative heartbeat, addressing inaccuracies due to beat-to-beat variations.  This single beat gives a comprehensive picture of that specific patient's filling pattern.  Together, these methods improve the clinical utility of IVPD, making it more robust and limiting the chance for a misdiagnosis. / Master of Science

We consider a gas of indistinguishable bosons, confined to a ring of radius R, and interacting via a pair-wise cosine potential. This may be thought of as the quantized Hamiltonian Mean Field (HMF) model for bosons originally introduced by Chavanis as a generalization of Antoni and Ruffo’s classical model. This thesis contains three parts: In part one, the dynamics of a Bose-condensate are considered by studying a generalized Gross-Pitaevskii equation (GGPE). Quantum effects due to the quantum pressure are found to substantially alter the system’s dynamics, and can serve to inhibit a pathological instability for repulsive interactions. The non-commutativity of the large-N , long-time, and classical limits is discussed. In part two, we consider the GGPE studied above and seek static solutions. Exact solutions are identified by solving a non-linear eigenvalue problem which is closely related to the Mathieu equation. Stationary solutions are identified as solitary waves (or solitons) due to their small spatial extent and the system’s underlying Galilean invariance. Asymptotic series are developed to give an analytic solution to the non- linear eigenvalue problem, and these are then used to study the stability of the solitary wave mentioned above. In part three, the exact solutions outlined above are used to study quantum fluctuations of gapless excitations in the HMF model’s symmetry broken phase. It is found that this phase is destroyed at zero temperature by large quantum fluctuations. This demonstrates that mean-field theory is not exact, and can in fact be qualitatively wrong, for long-range interacting quantum systems, in contrast to conventional wisdom. / Thesis / Doctor of Philosophy (PhD) / The Hamiltonian Mean Field (HMF) model was initially proposed as a simplified description of self-gravitating systems. Its simplicity shortens calculations and makes the underlying physics more transparent. This has made the HMF model a key tool in the study of systems with long-range interactions. In this thesis we study a quantum extension of the HMF model. The goal is to understand how quantum effects can modify the behaviour of a system with long-range interactions. We focus on how the model relaxes to equilibrium, the existence of special “solitary waves”, and whether quantum fluctuations can prevent a second order (quantum) phase transition from occurring at zero temperature.

The physical processes in lakes remain only partially understood despite successful data collection from a variety of sources spanning several decades. Although numerical models are already frequently employed to simulate the physics of lakes, especially in the context of water quality management, improved methods are necessary to better capture the wide array of dynamically important physical processes, spanning length scales from ~ 10 km (basin-scale oscillations) - 1 m (short internal waves). In this thesis, high-order numerical methods are explored for specialized model equations of lakes, so that their use can be taken into consideration in the next generation of more sophisticated models that will better capture important small scale features than their present day counterparts. The full three-dimensional incompressible density-stratified Navier-Stokes equations remain too computationally expensive to be solved for situations that involve both complicated geometries and require resolution of features at length-scales spanning four orders of magnitude. The main source of computational expense lay with the requirement of having to solve a three-dimensional Poisson equation for pressure at every time-step. Simplified model equations are thus the only way that numerical lake modelling can be carried out at present time, and progress can be made by seeking intelligent parameterizations as a means of capturing more physics within the framework of such simplified equation sets. In this thesis, we employ the long-accepted practice of sub-dividing the lake into vertical layers of different constant densities as an approximation to continuous vertical stratification. We build on this approach by including weakly non-hydrostatic dispersive correction terms in the model equations in order to parameterize the effects of small vertical accelerations that are often disregarded by operational models. Favouring the inclusion of weakly non-hydrostatic effects over the more popular hydrostatic approximation allows these models to capture the emergence of small-scale internal wave phenomena, such as internal solitary waves and undular bores, that are missed by purely hydrostatic models. The Fourier and Chebyshev pseudospectral methods are employed for these weakly non-hydrostatic layered models in simple idealized lake geometries, e.g., doubly periodic domains, periodic channels, and annular domains, for a set of test problems relevant to lake dynamics since they offer excellent resolution characteristics at minimal memory costs. This feature makes them an excellent benchmark to compare other methods against. The Discontinuous Galerkin Finite Element Method (DG-FEM) is then explored as a mid- to high-order method that can be used in arbitrary lake geometries. The DG-FEM can be interpreted as a domain-decomposition extension of a polynomial pseudospectral method and shares many of the same attractive features, such as fast convergence rates and the ability to resolve small-scale features with a relatively low number of grid points when compared to a low-order method. The DG-FEM is further complemented by certain desirable attributes it shares with the finite volume method, such as the freedom to specify upwind-biased numerical flux functions for advection-dominated flows, the flexibility to deal with complicated geometries, and the notion that each element (or cell) can be regarded as a control volume for conserved fluid quantities. Practical implementation details of the numerical methods used in this thesis are discussed, and the various modelling and methodology choices that have been made in the course of this work are justified as the difficulties that these choices address are revealed to the reader. Theoretical calculations are intermittently carried out throughout the thesis to help improve intuition in situations where numerical methods alone fall short of giving complete explanations of the physical processes under consideration. The utility of the DG-FEM method beyond purely hyperbolic systems is also a recurring theme in this thesis. The DG-FEM method is applied to dispersive shallow water type systems as well as incompressible flow situations. Furthermore, it is employed for eigenvalue problems where orthogonal bases must be constructed from the eigenspaces of elliptic operators. The technique is applied to the problem calculating the free modes of oscillation in rotating basins with irregular geometries where the corresponding linear operator is not self-adjoint.

Pressure swirl atomizers are commonly used in IC, aero-engines, and liquid propellant rocket combustion. Understanding the atomization process is important in order to enhance vaporization, mitigate soot formation, design of combustion chambers, and improve overall combustion efficiency. This work utilizes non-invasive techniques such as ultra -speed imaging, and Phase Doppler Particle Anemometry (PDPA) in order to investigate the cascade atomization process of pressure-swirl atomizers by examining swirling liquid film dynamics and the localized droplet characteristics of the resulting hollow cone spray. Specifically, experiments were conducted to examine these effects for three different nozzles with orifice diameters .3mm, .5mm, and .97mm. The ultra-speed imaging allowed for both visualization and interface tracking of the swirling conical film which emanated from each nozzle. Moreover, this allowed for the measurement of the radial fluctuations, film length, cone angle and maximum wavelength. Radial fluctuations are found to be maximum near the breakup or rupture of a swirling film. Film length decreases as Reynolds number increases. Cone angle increases until a critical Reynolds number is reached, beyond which it remains constant. A new approach to analyze the temporally unstable waves was developed and compared with the measured maximum wavelengths. The new approach incorporates the attenuation of a film thickness, as the radius of a conical film expands, with the classical dispersion relationship for an inviscid moving liquid film. This approach produces a new long wave solution which accurately matches the measured maximum wavelength swirling conical films generated from nozzles with the smallest orifice diameter. For the nozzle with the largest orifice diameter, the new long wave solution provides the upper bound limit, while the long wave solution for a constant film thickness provides the lower bound limit. These results indicate that temporal instability is the dominating mechanism which generates long Kelvin Helmholtz waves on the surface of a swirling liquid film. The PDPA was used to measure droplet size and velocity in both the near field and far field of the spray. For a constant Reynolds number, an increase in orifice diameter is shown to increase the overall diameter distribution of the spray. In addition, it was found that the probability of breakup, near the axis, decreases for the largest orifice diameter. This is in agreement with the cascading nature of atomization.

Non-destructive testing (NDT) techniques are critical to ensure integrity and safety of engineered structures. Structural health monitoring (SHM) is considered as the next step in the field enabling continuous monitoring of structures. The first part of the thesis concerns NDT and SHM using guided waves in plates, or Lamb waves, to perform imaging of plate structures. The imaging is performed using a fixed active array setup covering a larger area of a plate. Current methods are based on conventional beamforming techniques that do not efficiently exploit the available data from the small arrays used for the purpose. In this thesis an adaptive signal processing approach based on the minimum variance distortionless response (MVDR) method is proposed to mitigate issues related to guided waves, such as dispersion and the presence of multiple propagating modes. Other benefits of the method include a significant increase in resolution. Simulation and experimental results show that the method outperforms current standard processing techniques. The second part of the thesis addresses transducer design issues for resonant ultrasound inspections. Resonant ultrasound methods utilize the shape and frequency of the object's natural modes of vibration to detect anomalies. The method considered in the thesis uses transducers that are acoustically coupled to the inspected structures. Changes in the transducer's electrical impedance are used to detect defects. The sensitivity that can be expected from such a setup is shown to highly depend on the transducer resonance frequency, as well as the working frequency of the instrument. Through simulations and a theoretical argumentation, optimal conditions to achieve high sensitivity are given.

The central objective of this thesis is to suggest and develop one simple, unified method for communication over optical fiber networks, valid for all values of dispersion and nonlinearity parameters, and for a single-user channel or a multiple-user network. The method is based on the nonlinear Fourier transform (NFT), a powerful tool in soliton theory and exactly solvable models for solving integrable partial differential equations governing wave propagation in certain nonlinear media. The NFT decorrelates signal degrees of freedom in such models, in much the same way that the Fourier transform does for linear systems. In this thesis, this observation is exploited for data transmission over integrable channels such as optical fibers, where pulse propagation is governed by the nonlinear Schr\"odinger (NLS) equation. In this transmission scheme, which can be viewed as a nonlinear analogue of orthogonal frequency-division multiplexing commonly used in linear channels, information is encoded in the nonlinear spectrum of the signal. Just as the (ordinary) Fourier transform converts a linear convolutional channel into a number of parallel scalar channels, the nonlinear Fourier transform converts a nonlinear dispersive channel described by a \emph{Lax convolution} into a number of parallel scalar channels. Since, in the spectral coordinates the NLS equation is multiplicative, users of a network can operate in independent nonlinear frequency bands with no deterministic inter-channel interference. Unlike most other fiber-optic transmission schemes, this technique deals with both dispersion and nonlinearity directly and unconditionally without the need for dispersion or nonlinearity compensation methods. This thesis lays the foundations of such a nonlinear frequency-division multiplexing system.

The central objective of this thesis is to suggest and develop one simple, unified method for communication over optical fiber networks, valid for all values of dispersion and nonlinearity parameters, and for a single-user channel or a multiple-user network. The method is based on the nonlinear Fourier transform (NFT), a powerful tool in soliton theory and exactly solvable models for solving integrable partial differential equations governing wave propagation in certain nonlinear media. The NFT decorrelates signal degrees of freedom in such models, in much the same way that the Fourier transform does for linear systems. In this thesis, this observation is exploited for data transmission over integrable channels such as optical fibers, where pulse propagation is governed by the nonlinear Schr\"odinger (NLS) equation. In this transmission scheme, which can be viewed as a nonlinear analogue of orthogonal frequency-division multiplexing commonly used in linear channels, information is encoded in the nonlinear spectrum of the signal. Just as the (ordinary) Fourier transform converts a linear convolutional channel into a number of parallel scalar channels, the nonlinear Fourier transform converts a nonlinear dispersive channel described by a \emph{Lax convolution} into a number of parallel scalar channels. Since, in the spectral coordinates the NLS equation is multiplicative, users of a network can operate in independent nonlinear frequency bands with no deterministic inter-channel interference. Unlike most other fiber-optic transmission schemes, this technique deals with both dispersion and nonlinearity directly and unconditionally without the need for dispersion or nonlinearity compensation methods. This thesis lays the foundations of such a nonlinear frequency-division multiplexing system.