Tsunamis are among the most terrifying and complex physical phenomena potentially affecting almost all coastal regions of the Earth. Tsunami waves propagate in the ocean over thousands of kilometres away from their generating source at considerable speeds. Among several other tsunamis that occurred during the past decade, the 2004 Indian Ocean Tsunami and the 2011 Tohoku Tsunami in Japan, considered to be the deadliest and costliest natural disasters in the history of mankind, respectively, have hit wide stretches of densely populated coastal areas. During these major events, severe destruction of inland structures resulted from the action of extreme hydrodynamic forces induced by tsunami flooding. Subsequent field surveys in which researchers from the University of Ottawa participated ultimately revealed that, in contrast to seismic forces, such hydrodynamic forces are not taken into proper consideration when designing buildings for tsunami prone areas. In view of these limitations, a novel interdisciplinary hydraulic-structural engineering research program was initiated at the University of Ottawa, in cooperation with the Canadian Hydraulic Centre of the National Research Council, to help develop guidelines for the sound design of nearshore structures located in such areas. The present study aims to simulate the physical laboratory experiments performed within the aforementioned research program using a single-phase three-dimensional weakly compressible Smoothed Particle Hydrodynamics (SPH) numerical model. These experiments consist in the violent impact of rapidly advancing tsunami-like hydraulic bores with individual slender structural elements. Such bores are emulated based on the classic dam-break problem. The quantitatively compared measurements include the time-history of the net base horizontal force and of the pressure distribution acting on columns of square and circular cross-sections, as well as flow characteristics such as bore-front velocity and water surface elevation. Good agreement was obtained. Results show that the magnitude and duration of the impulsive force at initial bore impact depend on the degree of entrapped air in the bore-front. The latter was found to increase considerably if the bed of the experimental flume is covered with a thin water layer of even just a few millimetres. In order to avoid large fluctuations in the pressure field and to obtain accurate simulations of the hydrodynamic forces, a Riemann solver-based formulation of the SPH method is utilized. However, this formulation induces excessive numerical diffusion, as sudden and large water surface deformations, such as splashing at initial bore impact, are less accurately reproduced. To investigate this particular issue, the small-scale physical experiment of Kleefsman et al. (2005) is also considered and modeled. Lastly, taking full advantage of the validated numerical model to better understand the underlying flow dynamics, the influence of the experimental test geometry and of the bed condition (i.e. dry vs. wet) is investigated. Numerical results show that when a bore propagates over a wet bed, its front is both deeper and steeper and it also has a lower velocity compared to when it propagates over a dry bed. These differences significantly affect the pressure distributions and resulting hydrodynamic forces acting on impacted structures.

Currently, the accurate prediction of the impact of an extreme wave on infrastructure located near shore is difficult to assess. There is a lack of established methods to accurately quantify these impacts. Extreme waves, such as tsunamis generate, through breaking, extremely powerful hydraulic bores that impact and significantly damage coastal structures and buildings located close to the shoreline. The damage induced by such hydraulic bores is often due to structural failure. Examples of devastating coastal disasters are the 2004 Indian Ocean Tsunami, 2005 Hurricane Katrina and most recently, the 2011 Tohoku Japan Tsunami. As a result, more advanced research is needed to estimate the magnitude of forces exerted on structures by such bores. This research presents results of a numerical model based on the Smoothed Particle Hydrodynamics (SPH) method which is used to simulate the impact of extreme hydrodynamic forces on shore protection walls. Typically, fluids are modeled numerically based on a Lagrangian approach, an Eulerian approach or a combination of the two. Many of the common problems that arise from using more traditional techniques can be avoided through the use of SPH-based models. Such challenges include the model computational efficiency in terms of complexity of implementation. The SPH method allows water particles to be individually modeled, each with their own characteristics, which then accurately depicts the behavior and properties of the flow field. An open source code, known as SPHysics, was used to run the simulations presented in this thesis. Several cases analysed consist of hydraulic bores impacting a flat vertical wall as well as a sloping seawall. The analysis includes comparisons of the numerical results with published experimental data. The model is shown to accurately reproduce the formation of solitary waves as well as their propagation and breaking. The impacting bore profiles as well as the resulting pressures are also efficiently simulated using the model.

A two-dimensional axisymmetric implementation of the smoothed particle hydrodynamics (SPH) technique, called for short AxisSPH, has been described in this thesis, along with a number of basic tests and realistic applications. The main goal of this work was to fill a gap on a topic which has been scarcely addressed in the published literature concerning SPH. Although the application of AxisSPH to the simulation of real problems is restricted to those systems which display the appropriate symmetry there are, however, many interesting examples of physical systems which evolve following the axisymmetric premise. These examples belong to a variety of scientific and technological areas such as, for example, astrophysics, laboratory astrophysics or inertial confinement fusion. Additionally AxisSPH can be also useful in convergence studies of the standard 3D-SPH technique because the higher resolution achieved in 2D can be used to benchmark the three-dimensional codes. The main improvements implemented in AxisSPH with respect existing axisymmetric SPH formulations are summarized as follows: 1) We have derived simple analytical expressions for correction factors which largely improves the calculation of density and velocity in the vicinity of the z-axis. These expressions and their derivatives were given as a function of an adimensional parameter and do not increase the computational load of the scheme. 2) We have obtained the appropriate expression of the fluid Euler equations containing the new correction functions and their derivatives. Far enough from the singular axis, the scheme reduces to the standard formulation discussed by Brookshaw (2003). 3) A novel expression for the heat conduction term, which has to be added to the energy equation was devised and checked. This new term improves the description of the heat flux for those particles located at the axis neighborhoods. 4) Until now axisymmetric SPH hydrocodes handle artificial viscosity using a crude approach because it was treated as a simple restriction of the standard 3D Cartesian viscosity to 2D. Here we propose to calculate the viscous pressure as a combination of two terms, the first one is the (standard) Cartesian part and the second is the axis-converging part of the viscosity respectively. As expected this last term is of special relevance to simulate implosions. 5) We have developed an original method to incorporate gravity into AxisSPH. First the direct ring to ring force was found as a function of the Euclidean distance between the 2D particles. In second place the gravitational force on a given particle was obtained by summing the contributions of all N particles. We have also developed a more efficient scheme to obtain the gravitational force calculating the potential of the ring, instead the force because it involves lesser algebraic operations. The scheme has been checked using a large number of tests cases. These tests range from very specific oriented to check a particular algorithm or a piece of physics, to rather complex ones intended to analyze the behavior of the scheme in potential real applications (ICF, jets, astrophysics). At least in one case, the head on collision of a pair of white dwarfs, the result of the simulations carried out using AxisSPH brings new, unpublished, scientific material. / En esta tesis se ha desarrollado un código, que hemos llamado AxisSPH, en dos dimensiones axisimétrico a partir de la técnica conocida como SPH (“smooothed particle hydrodynamics”). AxisSPH ha sido validado después de realizar una serie de tests básicos y algunas simulaciones de situaciones reales. El objetivo principal de este trabajo ha sido llenar, en parte, el vacío existente al respecto en la literatura sobre SPH. Aunque sólo se puede aplicar AxisSPH en problemas reales que presenten la apropiada simetría, existen muchos ejemplos interesantes de sistemas físicos que presentan la simetría axial demandada. Existen ejemplos en campos de aplicación tanto científica como tecnológica, por ejemplo en astrofísica, en el llamado laboratorio de astrofísica o en fusión por confinamiento inercial (ICF). Otra interesante aplicación de AxisSPH puede ser su utilización en estudios de convergencia con otros códigos 3D-SPH debido a su mayor resolución, al tratarse de un código 2D. Las mejoras implementadas en el código AxisSPH en comparación con otros códigos axisimétricos SPH existentes se pueden resumir en los siguientes puntos: 1) Hemos deducido expresiones analíticas simples para unos factores de corrección que mejoran el cálculo de la densidad y la velocidad en las proximidades del eje z. Dichas expresiones y sus derivadas dependen de un parámetro adimensional que no incrementa mucho el peso computacional del esquema propuesto. 2) Hemos obtenido las expresiones adecuadas de las ecuaciones de Euler que contienen estas nuevas funciones correctoras y sus derivadas. Lejos del eje de singularidad estas ecuaciones se transforman en las de la formulación estándar propuesta por Brookshaw (2003). 3) Una expresión novedosa del término de conducción, que debe de añadirse a la ecuación de la energía, se ha propuesto y validado. Este nuevo término mejora la evolución del flujo de calor de las partículas situadas en las proximidades del eje z. 4) Hasta el momento los códigos hidrodinámicos SPH axisimétricos existentes trabajaban con una aproximación poco elaborada de la viscosidad artificial ya que consistían en una restricción a dos dimensiones de la viscosidad estándar 3D. En este trabajo proponemos el cálculo de la presión debida a la viscosidad como combinación de dos términos, el primero reflejo de la parte cartesiana y la segunda da cuenta de la parte relacionada con la convergencia en el eje. Como era de esperar este último término es de relevante importancia en la simulación de implosiones. 5) Hemos desarrollado un método original para incorporar el cálculo de la gravedad en el código AxisSPH. En primer lugar la fuerza directa de anillo a anillo y en segundo lugar la fuerza de la gravedad que sufre una determinada partícula a partir de la contribución del resto de las N partículas existentes. También hemos desarrollado un esquema más eficiente para calcular la gravedad a partir del cálculo del potencial del anillo en lugar del cálculo directo de la fuerza ya que implica un menor número de operaciones algebraicas. El método ha sido verificado con un gran número de test numéricos. Desde los más específicos orientados a comprobar la validez de un algoritmo particular o la capacidad para simular un fenómeno físico en particular, hasta simulaciones bastante más complejas, con la intención de validar la capacidad de simular aplicaciones potencialmente más reales (ICF, jets, astrofísica). Así, en al menos un caso, en la colisión frontal de dos enanas blancas, los resultados de la simulación utilizando AxisSPH pueden aportar material científico publicable.

Zugl.: Freiburg (Breisgau), Univ., Diss., 2009

Due to its adaptable nature in a broad range of problem domains, Smoothed Particle Hydrodynamics (SPH) is a popular numerical technique for computing solutions in astrophysics. This dissertation discusses the SPH technique and assesses its capabilities for reproducing steady-state spherically-symmetric accretion flow. The accretion scenario is of great interest for its applicability in a diverse array of astrophysical phenomena and, under certain assumptions, it also provides an accepted analytical solution against which the numerical method can be validated. After deriving the necessary equations from astrophysical fluid dynamics, giving a detailed review of solving the steady-state spherical accretion problem, and developing the SPH methodology, this work suggests solutions to the issues that must be overcome in order to successfully employ the SPH methodology to reproduce steady-state spherical accretion flow. Several techniques for setting initial data are addressed, resolution requirements are illustrated, inner and outer boundary conditions are discussed, and artificial dissipation parameters and methodologies are explored. / text

In this thesis, we apply smoothed empirical likelihood method to investigate confidence intervals for the receiver operating characteristic (ROC) curve with right censoring. As a particular application of comparison of distributions from two populations, the ROC curve is constructed by the combination of cumulative distribution function and quantile function. Under mild conditions, the smoothed empirical likelihood ratio converges to chi-square distribution, which is the well-known Wilks's theorem. Furthermore, the performances of the empirical likelihood method are also illustrated by simulation studies in terms of coverage probability and average length of confidence intervals. Finally, a primary biliary cirrhosis data is used to illustrate the proposed empirical likelihood procedure.

Currently, the accurate prediction of the impact of an extreme wave on infrastructure located near shore is difficult to assess. There is a lack of established methods to accurately quantify these impacts. Extreme waves, such as tsunamis generate, through breaking, extremely powerful hydraulic bores that impact and significantly damage coastal structures and buildings located close to the shoreline. The damage induced by such hydraulic bores is often due to structural failure. Examples of devastating coastal disasters are the 2004 Indian Ocean Tsunami, 2005 Hurricane Katrina and most recently, the 2011 Tohoku Japan Tsunami. As a result, more advanced research is needed to estimate the magnitude of forces exerted on structures by such bores. This research presents results of a numerical model based on the Smoothed Particle Hydrodynamics (SPH) method which is used to simulate the impact of extreme hydrodynamic forces on shore protection walls. Typically, fluids are modeled numerically based on a Lagrangian approach, an Eulerian approach or a combination of the two. Many of the common problems that arise from using more traditional techniques can be avoided through the use of SPH-based models. Such challenges include the model computational efficiency in terms of complexity of implementation. The SPH method allows water particles to be individually modeled, each with their own characteristics, which then accurately depicts the behavior and properties of the flow field. An open source code, known as SPHysics, was used to run the simulations presented in this thesis. Several cases analysed consist of hydraulic bores impacting a flat vertical wall as well as a sloping seawall. The analysis includes comparisons of the numerical results with published experimental data. The model is shown to accurately reproduce the formation of solitary waves as well as their propagation and breaking. The impacting bore profiles as well as the resulting pressures are also efficiently simulated using the model.

This thesis investigates tsunami impact on shore-based, low-rise structures in coastal areas. The aims are to investigate tsunami wave inundation in built-up coastal areas with reference to structural response to wave inundation, to assess the performance of current design codes in comparison with validated state-of-the-art numerical models and to improve structural design of residential buildings in tsunami risk areas. Tsunami events over the past few decades have shown that a significant proportion of fatalities can be attributed to the collapse of building infrastructure due to various actions of the incident waves. Although major tsunami events have demonstrated the potential catastrophic effects on built infrastructure, current building codes have no detailed or consistent guidance on designing structures in tsunami-prone regions. Furthermore, considerable differences in existing empirical formulae highlight that new research is necessary to appropriately address the particularities of the tsunami-induced forces and structure response into the design standards. In this thesis, numerical modelling methods are used to simulate hydrodynamic impact on shore-based coastal structures. The hydrodynamic simulations were conducted using a novel meshless numerical method, smoothed particle hydrodynamics (SPH), which is coupled with the finite element (FE) method to model structural behaviour. The SPH method was validated with experimental data for bore impact on an obstacle using a convergence study to identify the optimum particle size to capture the hydrodynamics. The FE model was validated against experimental data for plates under transient blast loads which have similar load characteristics with impulsive tsunami-induced bore impacts. One of the contributions of the thesis is the use of a new coupling method of the SPH-based software DualSPHysics and FE-based software ABAQUS. Using SPH particle spacing of the same size as the FE mesh size, enables the SPH output pressure to be directly applied as an input to the structural response model. Using this approach the effects of arrangement and orientation of single and multiple low rise structures are explored. Test cases were performed in 2-D and 3-D involving a discrete structure and multiple structures. The 3-D SPH simulations with single and multiple structures used an idealised coastal structure in the form of a cube with different on-plan orientations (0°, 30°, 45° and 60°) relative to the oncoming bore direction. The single structure cases were intended to study the improvement of the resilience of coastal structures by reducing the acting pressures on the vertical surfaces by changing the structure’s orientation. It was found the pressure exerted on the vertical surface of structure can be reduced by up to 50% with the 60° orientation case. The multiple structure models were conducted to examine shielding and flow focusing phenomena in tsunami events. The results reveal that the distance between two adjacent front structures can greatly influence the pressure exerted on the rear structure. This thesis also demonstrates the capability of SPH numerical method in simulating standard coastal engineering problems such as storm waves impact on a recurve wall in 2-D. The idealised structures were represented as standard timber construction and the finite element modelling was used to determine the corresponding stress distributions under tsunami impact. Following the comparison of the method used in this thesis with commonly used design equations based on the quasi-static approach, large differences in stress prediction were observed. In some cases the loads according to the design equations predicted maximum stresses almost one order of magnitude lower. This large discrepancy clearly shows the potential for non-conservative design by quasi-static approaches. The new model for the simulation of tsunami impact on discrete and multiple structures shows that the resilience of a coastal structure can be improved by changing the orientation and arrangement. The characteristics of tsunami waves during propagation and bore impact pressures on structures can be assessed in great detail with the combined SPH and FE modelling strategy. The techniques outlined in this thesis will enable engineers to gain a better insight into tsunami wave-structure interaction with a view towards resilience optimisation of structures vulnerable to tsunami impact events.

No presente trabalho, foi implementada uma formulação de otimização topológica com o objetivo de encontrar o mínimo volume de estruturas contínuas bidimensionais, em estado plano de tensão, sujeitas à restrição de tensão de von Mises. Foi utilizado o Método dos Elementos Finitos para discretizar o domínio, com o elemento não conforme de Taylor. A tensão foi suavizada, calculando-se um valor de tensão para cada nó do elemento. O fenômeno da singularidade foi contornado através do método de relaxação da tensão, penalizando-se o tensor constitutivo. Foi usada uma única medida de tensão global, a normap, resultando na redução do custo computacional do cálculo das sensibilidades. As sensibilidades da função objetivo e da restrição de tensão foram calculadas analiticamente. O problema de otimização topológica foi resolvido por um algoritmo de Programação Linear Sequencial. Os fenômenos da instabilidade de tabuleiro e da dependência da malha foram contornados pela utilização de um filtro de densidade linear. A formulação desenvolvida foi testada em 3 casos clássicos. No primeiro deles, foi testada uma viga curta em balanço, submetida a 3 diferentes tipos de penalização da função objetivo, obtendo-se uma estrutura com 27% do volume inicial, com reduzido número de elementos com densidades intermediárias. No segundo caso, foi testada a mesma estrutura submetida à flexão, chegandose a uma topologia bem definida no formato de duas barras, com 16,25% do volume inicial. No terceiro caso, em que foi utilizado um componente estrutural em formato de “L”, justamente por favorecer o surgimento de concentração de tensão em sua quina interna, o otimizador gerou uma estrutura bem definida, permanecendo, contudo, uma pequena região de concentração de tensão na topologia final. / A topology optimization formulation to search for the minimum volume of twodimensional linear elastic continuous structures in plane stress, subject to a von Mises stress constraint, was implemented in this study. The extended domain was discretized using Taylor nonconforming finite element. Nodal values of the stress tensor field were computed by global smoothing. A penalized constitutive tensor stress relaxation method bypassed the stress singularity problem. A single p-norm global stress measure was used to speed up the sensitivity analysis. The sensitivities of the objective function and stress constraints were derived analytically. The topology optimization problem was solved by a Sequential Linear Programming algorithm. A linear density filter avoided the checkerboard and the mesh dependence phenomena. The formulation was tested with three benchmark cases. In the first case, a tip loaded short cantilever beam was optimized using a sequence of three different objective function penalizations. The converged design had approximately 27% of the initial volume, with a small proportion of intermediate densities areas. In the second case, the same domain was subjected to shear, resulting a well defined two-bar design, with 16.25% of the initial volume. In the third case, an L-shape structure was studied, because it has a stress concentration at the reentrant corner. In this last case, the final topology was well-defined, but the stress concentration was not completely removed.

No presente trabalho, foi implementada uma formulação de otimização topológica com o objetivo de encontrar o mínimo volume de estruturas contínuas bidimensionais, em estado plano de tensão, sujeitas à restrição de tensão de von Mises. Foi utilizado o Método dos Elementos Finitos para discretizar o domínio, com o elemento não conforme de Taylor. A tensão foi suavizada, calculando-se um valor de tensão para cada nó do elemento. O fenômeno da singularidade foi contornado através do método de relaxação da tensão, penalizando-se o tensor constitutivo. Foi usada uma única medida de tensão global, a normap, resultando na redução do custo computacional do cálculo das sensibilidades. As sensibilidades da função objetivo e da restrição de tensão foram calculadas analiticamente. O problema de otimização topológica foi resolvido por um algoritmo de Programação Linear Sequencial. Os fenômenos da instabilidade de tabuleiro e da dependência da malha foram contornados pela utilização de um filtro de densidade linear. A formulação desenvolvida foi testada em 3 casos clássicos. No primeiro deles, foi testada uma viga curta em balanço, submetida a 3 diferentes tipos de penalização da função objetivo, obtendo-se uma estrutura com 27% do volume inicial, com reduzido número de elementos com densidades intermediárias. No segundo caso, foi testada a mesma estrutura submetida à flexão, chegandose a uma topologia bem definida no formato de duas barras, com 16,25% do volume inicial. No terceiro caso, em que foi utilizado um componente estrutural em formato de “L”, justamente por favorecer o surgimento de concentração de tensão em sua quina interna, o otimizador gerou uma estrutura bem definida, permanecendo, contudo, uma pequena região de concentração de tensão na topologia final. / A topology optimization formulation to search for the minimum volume of twodimensional linear elastic continuous structures in plane stress, subject to a von Mises stress constraint, was implemented in this study. The extended domain was discretized using Taylor nonconforming finite element. Nodal values of the stress tensor field were computed by global smoothing. A penalized constitutive tensor stress relaxation method bypassed the stress singularity problem. A single p-norm global stress measure was used to speed up the sensitivity analysis. The sensitivities of the objective function and stress constraints were derived analytically. The topology optimization problem was solved by a Sequential Linear Programming algorithm. A linear density filter avoided the checkerboard and the mesh dependence phenomena. The formulation was tested with three benchmark cases. In the first case, a tip loaded short cantilever beam was optimized using a sequence of three different objective function penalizations. The converged design had approximately 27% of the initial volume, with a small proportion of intermediate densities areas. In the second case, the same domain was subjected to shear, resulting a well defined two-bar design, with 16.25% of the initial volume. In the third case, an L-shape structure was studied, because it has a stress concentration at the reentrant corner. In this last case, the final topology was well-defined, but the stress concentration was not completely removed.