Les Moteurs à Propergol Solide sont le siège d'instabilités de combustion qui se manifestent par des fluctuations de la pression interne et, en conséquence, de la poussée induite. Les phénomènes oscillatoires observés résultent d'un couplage entre l'acoustique de la chambre de combustion et une ou plusieurs instabilités hydrodynamiques de l'écoulement interne.Dans cette thèse, on s'intéresse à l'étude de l'écoulement induit par injection pariétale en géométrie axisymétrique. Plus particulièrement, on se concentre sur l'interaction qui peut avoir lieu entre cet écoulement et une structure : les Protections Thermiques de Face.Une étude expérimentale reposant sur le montage gaz froid VALDO de l'ONERA, modélisation représentative des MPS P230, a permis d'examiner l'influence d'obstacles, rigides et en élastomère, sur les instationnarités de l'écoulement. Le caractère instable de l'écoulement induit par injection pariétale est alors confirmé et le rôle primordial joué par les PTF sur les fréquences émergeant au sein de l'écoulement a été identifié.L'analyse de stabilité linéaire de l'écoulement, conduite en approche globale, permet l'étude de sa dynamique en tant qu'amplificateur de bruit. Les m¶mécanismes mis en jeu dans l'instabilité de l'écoulement induit par injection pariétale ainsi que son interaction complexe avec la couche de cisaillement issue de l'obstacle sont analysés. La réponse de l'écoulement à un forçage harmonique est alors examinée à la lumière des résultats expérimentaux.Cette étude a été complétée par une approche numérique de l'Interaction Fluide-Structure sur une configuration jugée critique. Une étape préliminaire à la simulation multi-physique est l'élaboration d'un modèle réaliste du comportement des PTF en élastomère. L'identification et la prise en considération des mécanismes propres aux élastomères dans la loi de comportement est ainsi détaillée. Les résultats issus des simulations sont confrontés aux résultats expérimentaux et rendent possible la validation d'un scénario de couplage dédié à l'Interaction Fluide-Structure au sein des MPS. / Fluid Structure Interaction of an inhibitor with the internal flow induced by wall injection was studied in an axisymmetric cold flow apparatus. Experiments were carried out over a wide range of injection velocities in order to underline how the obstacle not only modifies the mean flowfield but also its entire dynamic behaviour.The resulting instability (from the interaction of the unstable shear layer with the Taylor-Culick flow) exhibits a significant shift with respect to the Taylor-Culick instability and therefore emphasizes the strong impact of the inhibitor on hydrodynamics.The mecanisms responsible of such behaviour are studied in this thesis with the help of global linear stability analysis and multi-physics numerical computations.

O fenômeno de vibração induzida por vórtices (VIV) é um problema fundamental dentro da Mecânica dos Fluidos e um exemplo importante de interação fluido-estrutura. Esta tese investiga fenômeno de VIV quando um cilindro rígido, submetido a escoamento uniforme, está livre para oscilar na direção transversal e alinhada com a corrente incidente. A tese foi estruturada ao redor de sete perguntas relacionadas ao fenômeno de VIV: 1) O fenômeno e os resultados experimentais são repetitivos? 2) Como ocorre a transição entre ramos de resposta? 3) Qual é o papel da inércia da estrutura oscilante? 4) Qual é o papel de sua rigidez? 5) Quais são as frequências naturais mais importantes da estrutura? 6) Quais padrões de esteira se desenvolvem para VIV com dois graus de liberdade? 7) Quais são os efeitos do movimento na direção alinhada com a corrente no processo de formação e desprendimento de vórtices? O fenômeno de VIV é estudado de maneira experimental em uma base elástica pendular capaz de oscilar com o mesmo momento de inércia e frequência natural nas duas direções. Os experimentos de VIV foram realizados em canal de água recirculante e com diferentes condições de inércia e rigidez. A técnica de velocimetria por imagem de partículas foi usada e permitiu identificar diferentes padrões de esteira de vórtices. Verificou-se que o VIV é repetitivo a nível de amplitudes médias e frequências dominantes. A transição dos ramos pode ocorrer de maneira intermitente ou com histerese. Os parâmetros de inércia e rigidez da estrutura são capazes de mudar o regime de oscilação e, para algumas condições, suprimir as vibrações alinhadas com a corrente. Dentre os padrões de esteira observados, um deles não havia sido relatado na literatura e é definido nesta tese. O novo modo de emissão apresenta dois vórtices com circulação oposta e elevada intensidade emitidos por ciclo. A influência da direção alinhada com o escoamento está relacionada a dois efeitos: a velocidade relativa entre o cilindro e o fluido, responsável pelo aumento da circulação dos vórtices na esteira, e o ângulo de fase do movimento nas direções alinhada e transversal, capaz de mudar o processo de formação dos vórtices. / Vortex-induced vibration (VIV) phenomenon is a fundamental problem of Fluid Mechanics and a typical example of fluid-structure interaction. This thesis explores the VIV phenomenon for a rigid circular cylinder immersed in a uniform fluid current. The cylinder is free to oscillate in two degrees of freedom (2dof): in-line and cross flow. The thesis was structured in order to answer seven questions regarding VIV: 1) Are the phenomenon and its experimental results repetitive? 2) How do the branches transition occur? 3) What is the role played by the inertia of the oscillating structure? 4) What is the role of its stiffness? 5) Which are the most relevant natural frequencies of the structure? 6) Which are the vortex wake patterns developed in VIV 2dof? 7) What are the influences of the in-line movement to the process of vortex formation and shedding? The phenomenon is experimentally investigated using an elastic base similar to a pendulum and able to oscillate with the same moment of inertia and natural frequencies in both directions. All the experiments were conducted in a recirculating water channel facility and with several combinations of moment of inertia and stiffness of the structure. Particle image velocimetry provided visualization of different vortex wake patterns. The phenomenon is repetitive in terms of its mean amplitudes and dominant frequencies. The transitions between dfferent branches can be hysteretic or intermittent. It is shown that both the moment of inertia and the stiffness of the structure are able to change the regime of oscillations and, for some cases, suppress the in-line movement. Among the different vortex wake patterns observed, one has not been reported previously in the literature. The new wake pattern shows two large vortices with high and opposite circulations shed per cycle. The influence of the displacement in the current direction is related to two different effects: the relative velocity between the incoming flow and the structure motion, responsible for the increase in the net circulation shed in the vortex wake, and the influence of the phase angle between the displacement in the in-line and cross-flow directions, capable of changing the vortex formation process.

Apresentam-se as bases teóricas do Método de Elementos Finitos (MEF) piezoelétrico, e a sua aplicação na modelagem de transdutores de ultra-som piezoelétricos, que consiste na determinação das características vibracionais (frequências de ressonância e anti-ressonância, modos de vibrar e coeficiente de acoplamento eletromecânico), obtenção da curva de admitancia, análise transiente da estrutura piezoelétrica sujeita a uma excitação pulsada e análise da influência da variação das constantes piezoelétricas do transdutor com o raio. Utilizando-se o MEF aplicado a acústica obteve-se o campo acústico gerado pelo transdutor operando em onda contínua, bem como iniciou-se o estudo da propagação de ondas num líquido, analisando-se as ondas geradas pela excitação pulsada de um pistão plano em contato com o fluido. Os modos de vibrar e os valores de frequências de ressonância obtidos para um transdutor, foram comparados com os resultados experimentais. / The theoretical basis of piezoelectric finite element method (FEM), and its application in piezoelectric ultrasonic transducer modelling is presented. Among these applications we have the calculation of resonance and antiresonance frequencies, vibration modes, piezoelectric coupling coefficient, admittance curve and transient analysis of piezoelectric structure excited by a short pulse. By means of piezoelectric FEM the influence of variation of piezoelectric constant with radius is analysed. It is discussed three kind of functions (linear, cosinoidal and Gaussian). This technique is called apodization. The acoustic filed generated by the transducer operating in continuous wave (CW) was calculated by using FEM applied to acoustic, considering the fluid-structure coupling. The study of wave propagation in liquids is started by using FEM, analyzing the waves generated by a plane piston in contact with the fluid, excited by a short pulse. For each case discussed above, all boundary conditions and hypothesis assumed in the construction of finite element models are discussed. Although the models considered are circular transducers, the concepts acquired can be expanded to other geometries. The vibrational modes were visualized by means of a laser interferometry technique (ESPI), and the admittance curves were measured by using an impedometer. These results were compared with the FEM results, and the models precision was discussed.

Submitted by Erika Demachki (erikademachki@gmail.com) on 2014-10-21T20:53:51Z No. of bitstreams: 5 Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 01.pdf: 20230731 bytes, checksum: 56aa4f90641e87e6081a5e3e9319bd24 (MD5) Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 02.pdf: 10802505 bytes, checksum: 4e34e2a114366cce38898dda48f31144 (MD5) Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 03.pdf: 14680983 bytes, checksum: e7a40859c186d5c2130a3debfec64edd (MD5) Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 04.pdf: 11913507 bytes, checksum: 73e7928267b0b902da703843e18b2919 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Approved for entry into archive by Jaqueline Silva (jtas29@gmail.com) on 2014-10-22T19:05:17Z (GMT) No. of bitstreams: 5 Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 01.pdf: 20230731 bytes, checksum: 56aa4f90641e87e6081a5e3e9319bd24 (MD5) Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 02.pdf: 10802505 bytes, checksum: 4e34e2a114366cce38898dda48f31144 (MD5) Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 03.pdf: 14680983 bytes, checksum: e7a40859c186d5c2130a3debfec64edd (MD5) Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 04.pdf: 11913507 bytes, checksum: 73e7928267b0b902da703843e18b2919 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Made available in DSpace on 2014-10-22T19:05:17Z (GMT). No. of bitstreams: 5 Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 01.pdf: 20230731 bytes, checksum: 56aa4f90641e87e6081a5e3e9319bd24 (MD5) Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 02.pdf: 10802505 bytes, checksum: 4e34e2a114366cce38898dda48f31144 (MD5) Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 03.pdf: 14680983 bytes, checksum: e7a40859c186d5c2130a3debfec64edd (MD5) Dissertação - Ana Larissa dal Piva Argenta - 2013 - Parte 04.pdf: 11913507 bytes, checksum: 73e7928267b0b902da703843e18b2919 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Previous issue date: 2013-06-13 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / Circular cylindrical shells are widely used structures in several engineering areas and have great capacity to withstand both axial and lateral loads. However, they may present a complex dynamic behavior. Thus, a detailed study of the behavior of cylindrical shells subjected to different loading and support conditions as well as the influence of material characteristics and geometric relations evaluation is justified. In this work the influence of geometry and orthotropy on the nonlinear dynamic behavior of orthotropic simply supported cylindrical shells subjected to both axial and lateral time depending loads is studied. To model the shell, the Donnell nonlinear shallow shell theory, neglecting the effects of shear deformations, is used. It is considered the shell in three different situations: empty, filled with static fluid and subjected to internal flow of incompressible and non-viscous fluid, whose motion is isentropic and irrotational. The radial displacements are described as an expansion with eight degrees of freedom which satisfies the boundary conditions. The Galerkin method is applied to obtain a set of nonlinear equations of motion, which are in turn solved by the Runge-Kutta method. A detailed analysis is performed to study the influence of material orthotropy and geometric relations such as length-radio (L/R) and radio-thickness (R/h) on the natural frequencies, critical loads, critical flow velocities, post-critical paths, frequency-amplitude relations, instability boundaries, bifurcation diagrams and resonance curves. Obtained results display the strong influence of both material orthotropy and geometric relations on the linear and nonlinear behavior of the shells and, depending on these characteristics, the shell can display softening or hardening behavior. / Cascas cilíndricas são estruturas com diversas aplicações em várias áreas da engenharia e têm grande capacidade para resistir a carregamentos axiais e a pressões laterais. Entretanto, são estruturas que podem apresentar um complexo comportamento dinâmico. Assim, um estudo detalhado do comportamento das cascas cilíndricas submetidas a diferentes condições de carregamento bem como a avaliação da influência do material e da geometria se justifica plenamente. Este trabalho tem como objetivo estudar a influência da ortotropia do material e da geometria no comportamento dinâmico não linear das cascas cilíndricas ortotrópicas simplesmente apoiadas e submetidas a carregamentos axiais e laterais variáveis com o tempo. Para modelar a casca é utilizada a teoria não linear de Donnell para cascas abatidas desprezando-se os efeitos das deformações cisalhantes. Considera-se a casca em três situações distintas: vazia, preenchida com fluido estático e submetida ao escoamento interno de fluido incompressível e não viscoso, cujo movimento é isentrópico e irrotacional. O deslocamento radial da casca é descrito, de maneira geral, por uma expansão com oito graus de liberdade que satisfaz as condições de contorno. O método de Galerkin é utilizado para obter o sistema discreto de equações diferencias não lineares de movimento, que são resolvidas através do método de Runge-Kutta de quarta ordem. Uma análise detalhada é realizada visando observar a influência das características do material e das relações geométricas comprimento-raio (L/R) e raio-espessura (R/h) nas frequências naturais, cargas críticas, velocidades críticas do fluido bem como nos caminhos pós-críticos, relações frequência-amplitude, fronteiras de instabilidade, diagramas de bifurcação e curvas de ressonância das cascas. Os resultados obtidos permitem observar a forte influência que a geometria e as propriedades dos materiais exercem no comportamento linear e não linear das cascas cilíndricas ortotrópicas e, dependendo de certas características, verifica-se que a casca pode ter comportamentos com ganho ou perda de rigidez.

Problemas envolvendo interação entre fluido e estrutura são desafiadores para a engenharia e, ao mesmo tempo em que abrangem dois meios com características físicas distintas, demandam uma descrição matemática para cada um deles que seja compatível, de forma a permitir o acoplamento. Assim, este trabalho apresenta uma formulação em descrição Lagrangeana para análises dinâmicas de sólidos, fluidos incompressíveis e interação fluido-estrutura (IFE). Nos problemas de IFE é comum a estrutura apresentar grandes deslocamentos, o que torna imprescindível considerar o efeito da não-linearidade geométrica. Levando isso em consideração, é empregada uma formulação do método dos elementos finitos (MEF) baseada em posições, cuja aplicação em análises dinâmicas de estruturas em regime de grandes deslocamentos vem se mostrando bastante robusta. Já no âmbito da dinâmica dos fluidos, sabe-se que uma descrição Lagrangeana acaba por eliminar os termos convectivos das equações de Navier-Stokes, dispensando o uso de métodos estabilizantes nessas equações. Por outro lado, a dificuldade é então transferida para o uso de técnicas eficientes de remesh, preservação da qualidade da malha e de identificação do contorno, uma vez que os fluidos podem deformar-se indefinidamente quando submetidos a forças de cisalhamento. Assim, uma combinação do método dos elementos finitos e do método de partículas é utilizada, onde as forças de interação entre as partículas de fluido são calculadas por meio de uma malha de elementos finitos que é renovada para cada passo de tempo. Por meio de técnicas que reconstroem automaticamente o contorno, é possível simular problemas de superfície livre que sofram severas alterações e, até mesmo, uma eventual separação de partículas do domínio inicial, representando, por exemplo, a formação de gotas. Por fim, o sistema de acoplamento entre o fluido e o sólido é simplificado devido a ambos os domínios serem descritos através de um referencial Lagrangeano, não necessitando de métodos para a adaptação da malha do fluido de modo a acompanhar o movimento da estrutura. / Problems involving fluid-structure interaction are challenging for engineering and, while involving two different materials with distinct physical properties, they require a compatible mathematical description for both solid and fluid domain in order to allow the coupling. Thus, this work introduces a formulation, under Lagrangian description, for the solution of solid, incompressible fluid dynamics and fluid-structure interaction (FSI). In FSI problems, the structure usually presents large displacements thus making mandatory a geometric non-linear analysis. Considering it, we adopt a position based formulation of the finite element method (FEM) which has been shown to be very robust when applied to large displacement solid dynamics. For the fluid mechanics problem it is well known that a Lagrangian description eliminates the convective terms from the Navier-Stokes equations and thus, no stabilization technique is required. However, the difficulty is then transferred to the need of efficient re-meshing, mesh quality and external boundary identification techniques, since the fluid presents no resistance to shear stresses and may deform indefinitely. In this sense, we employ a combination of finite element and particle methods in which the particle interaction forces are computed by mean of a finite element mesh which is re-constructed at every time step. Free surface flows are simulated by a boundary recognition technique enabling large domain distortions or even the particles separation from the main domain, representing for instance a water drop. Finally, the fluid-structure coupling is simplified due to the Lagrangian description adopted for both materials, with no need for extra adaptive mesh-moving technique for the fluid computational domain to follow the structure motion.

Ce travail de thèse porte sur l'étude d'un système d'interaction fluide-structure. Nous en traitons de nombreux aspects allant de sa modélisation jusqu'à l'étude de sa stabilisation et de sa simulation numérique. Le premier chapitre du manuscrit aborde la modélisation du système ainsi que l'existence de solutions fortes en temps petits. Le fluide est représenté par les équations de Navier-Stokes incompressibles. La structure est déformable et dépend d'un nombre fini de paramètres. Nous obtenons ses équations en appliquant un principe des travaux virtuels. Le système d'équations final est non linéaire. Nous prouvons l'existence locale d'une solution à ce système, dans un premier temps sur le système linéarisé autour de l'état nul. Puis, nous prouvons l'existence de solutions en temps petits au système non linéaire grâce à un argument de point fixe. Le deuxième chapitre traite de la stabilisation par feedback autour d'un état stationnaire non nul du système présenté dans le Chapitre 1. L'opérateur de feedback est déterminé à partir de l'analyse du problème linéarisé autour de l'état stationnaire et de la résolution d'une équation de Riccati. Le résultat de stabilisation portant sur le système non linéaire requiert des données petites et est obtenu par un argument de point fixe. Le troisième chapitre se concentre sur les aspects numériques de ce problème. La construction de l'opérateur de feedback correspond à la version discrétisée de celle proposée dans le Chapitre 2. Le système fluide-structure est simulé en utilisant une méthode de domaines fictifs. / This PhD thesis deals with the study of a fluid-structure interaction system. We are interested in several aspects such as modelling, stabilization and numerical simulation. In the first chapter of the manuscript, we show the modelling of the system and prove the existence of strong solutions in small times. The fluid is modelled by the incompressible Navier- Stokes equations. The structure is deformable and depends on a finite number of parameters. The equations are obtained with a virtual work principle. The final system of equations is nonlinear. We prove local existence of a solution to this system, first on the linearized system. Then, existence of solutions in small times to the full nonlinear system is obtained with a fixed point argument. In the second chapter, we prove feedback stabilization of the problem around a non-null stationary state. The feedback operator is computed with the solution to a Riccati equation obtained by the analysis of the linearized problem around the stationary state. The stabilization result holds on the full nonlinear system and requires small data. It is proven by a fixed point argument. In the third chapter, we focus on the numerical aspects of the problem. The feedback operator used corresponds to a discretization of the feedback operator of Chapter 2. The solution to the full nonlinear system is computed by the use of a fictitious domain method.

Dans cette thèse, nous étudions les propriétés des écoulements de fluides qui interagissent avec un corps rigide ou avec une source et un puits. Dans le cas d'un fluide visqueux incompressible qui satisfait les équations de Navier Stokes dans un domaine borné 2D, les solutions faibles de Leray-Hopf sont bien comprises. L'existence et l'unicité sont prouvées. De plus, les solutions sont continues en temps `a valeurs dans L 2 (Omega) et satisfont l’égalité d'énergie classique. Plus récemment, le problème d'un corps rigide en mouvement dans un fluide visqueux incompressible modélisé par les équations de Navier-Stokes couplées aux lois de Newton qui décrivent le mouvement du solide a également été abordé dans le cas où des conditions aux limites sans glissement ont été prescrites. Des résultats analogues concernant les solutions de Leray-Hopf ont également été démontrés dans ce contexte. Dans ce manuscrit, nous étudions le cas de conditions aux limites de Navier-Slip. Dans ce cadre, le résultat d'existence pour le système couplé a été prouvée par G'erard-Varet et Hillairet en 2014. Ici, nous montrons que les solutions sont continues en temps, qu'elles satisfont l’égalité d'énergie et qu’elles sont uniques. De plus, nous montrons un résultat d'existence des solutions faibles dans le cas d'un fluide incompressible visqueux auquel s'ajoute un corps rigide dans le cas où la vitesse du fluide a une partie orthonormale d'énergie infinie.Pour un fluide incompressible non visqueux modélisé par les équations d'Euler dans un domaine borné 2D, le cas où le fluide est autorisé à entrer et à sortir de la frontière a été traité par Judovic qui a introduit certaines conditions limites consistant à prescrire la composante normale de la vitesse et de la vorticité entrante. Dans ce manuscrit, nous considérons un domaine borné qui possède deux trous. L'un d'eux est une source, ce qui signifie que le fluide est autorisé à entrer dans le domaine et l'autre est un puits où le fluide peut sortir. En particulier, nous établissons les équations limites vérifiées par le fluide lorsque la source et le puits se contractent en deux points différents. Le système limite est caractérisé par un point source/puits et un point vortex en chacun des deux points où les trous se sont contractés. / In this thesis, we investigate properties of incompressible flows that interact with a rigid body or a source and a sink. In the case of an incompressible viscous fluid that satisfies the Navier Stokes equations in a 2D bounded domain well-posedness of Leray-Hopf weak solutions is well-understood. Existence and uniqueness are proved. Moreover solutions are continuous in time with values in L 2 (Omega) and they satisfy the energy equality. Recently the problem of a rigid body moving in a viscous incompressible fluid modeled by the Navier-Stokes equations coupled with the Newton laws that prescribe the motion of the solid, was also tackled in the case where the no-slip boundary conditions were imposed. And the correspondent well-posedness result for Leray-Hopf type weak solutions was proved. In this manuscript we consider the case of the Navier-slip boundary conditions. In this setting, the existence result for the coupled system was proved by G'erard-Varet and Hillairet in 2014. Here, we prove that solutions are continuous in time, that they satisfy the energy equality and that they are unique. Moreover we show an existence result for weak solutions of a viscous incompressible fluid plus rigid body system in the case where the fluid velocity has an orthoradial part of infinite energy.For an inviscid incompressible fluid modelled by the Euler equations in a 2D bounded domain, the case where the fluid is allowed to enter and to exit from the boundary was tackled by Judovic who introduced some conditions which consist in prescribing the normal component of the velocity and the entering vorticity. In this manuscript we consider a bounded domain with two holes, one of them is a source which means that the fluid is allowed to enter in the domain and the other is a sink from where the fluid can exit. In particular we find the limiting equations satisfied by the fluid when the source and the sink shrink to two different points. The limiting system is characterized by a point source/sink and a point vortex in each of the two points where the holes shrunk.

Since 2004 there have been multiple devastating tsunamis around the globe triggered by large magnitude earthquakes; with the most recent being the Tohoku, Japan tsunami in March 2011. These tsunamis have caused significant loss of life and damage to the coastal communities impacted by these powerful waves. The resulting devastation has raised awareness of the dangers of tsunamis and the Network for Earthquake Engineering Simulation (NEES) housesmash project (NEEShousesmash), was started to investigate several different areas of tsunami inundation. The work presented in the following two manuscripts was performed at the O.H. Hinsdale Wave Lab and Gene D. Knudson Wood Engineering Lab, which are located at Oregon State University. This work represents a small portion of the total NEEShousesmash project, and is focused on improving the knowledge and predictability of tsunami loading and structural performance. The first manuscript investigates tsunami wave impact on full scale light-frame wood walls, and compares the measured forces to calculated values using the linear momentum equation, previously evaluated by Cross (1967). The results show for each wave height tested a peak transient force followed by a sustained quasi-static force, with a ratio of transient force to quasi-static force of 2.2. The results also show that the linear momentum equation did an acceptable job of predicting the measured transient forces on the walls to within ±10%, and that increased wall flexibility, 2x4 vs. 2x6 dimensional lumber, resulted in lower measured transient forces when subjected to similar tsunami wave heights. These results are important for practical use because the linear momentum equation is a simple equation to use, that only requires a couple of site specific input variables. The second manuscript is a continuation of the work done in the wave lab for the first manuscript. These experiments provide a starting point for expanding the testing of the structural response and performance of larger scale structures subjected to tsunami wave loads. By simulating tsunami loading in a traditional structures laboratory, the inherent limits of testing structural performance in small scale tsunami laboratory facilities is removed. The results show that a light-frame wood shear wall, built to current standards, is susceptible to premature failures from concentrated impact loads at intermediate heights compared to the design strength at full height. It is also shown that the out-of-plane walls subjected to both elastic and inelastic loads behave like a one way slab with minimal load sharing between adjacent studs. The failures observed during the hydrodynamic wave testing of the nailed connection between the bottom plate and studs was successfully reproduced, and shows that current construction standards are not fully utilizing the available capacity of each stud when subjected to tsunami waves. The reinforcement of this connection with traditional metal brackets would help increase the capacity of the out-of-plane wall to resist tsunami wave loads. / Graduation date: 2012

For the numerical simulation of many problems of engineering interest, it is desirable to have an automated mesh adaption tool. This is important especially for problems characterized by anisotropic features and require mesh clustering in the direction of high gradients. Another significant issue in meshing emerges in unsteady simulations with moving boundaries, where the boundary motion has to be accommodated by deforming the computational grid. Similarly, there exist problems where current mesh needs to be adapted to get more accurate solutions. To solve these problems, we propose three novel procedures. In the first part of this work, we present an optimization procedure for three-dimensional anisotropic tetrahedral grids based on metric-driven h-adaptation. Through the use of topological and geometrical operators, the mesh is iteratively adapted until the final mesh minimizes a given objective function. We propose an optimization process based on an ad-hoc application of the simulated annealing technique, which improves the likelihood of removing poor elements from the grid. Moreover, a local implementation of the simulated annealing is proposed to reduce the computational cost. Many challenging unsteady multi-physics problems are characterized by moving boundaries and/or interfaces. When the boundary displacements are large, degenerate elements are easily formed in the grid such that frequent remeshing is required. We propose a new r-adaptation technique that is valid for all types of elements (e.g., triangle, tet, quad, hex, hybrid) and deforms grids that undergo large imposed displacements at their boundaries. A grid is deformed using a network of linear springs composed of edge springs and a set of virtual springs. The virtual springs are constructed in such a way as to oppose element collapsing. Both frequent remeshing, and exact-pinpointing of clustering locations are great challenges of numerical simulations, which can be overcome by adaptive meshing algorithms. Therefore, we conclude this work by defining a novel mesh adaptation technique where the entire mesh is adapted upon application of a force field in order to comply with the target mesh or to get more accurate solutions. The method has been tested for two-dimensional problems of a-priori metric definitions as well as for oblique shock clusterings.

Liquid fragmentation phenomenon is explored from both a fundamental (fully resolved) and an engineering (modeled) perspective. The dual objectives compliment each other by providing an avenue to gain further understanding into fundamental processes of atomization as well as to use the newly acquired knowledge to address practical concerns. A compressible five-equation interface model based on a Roe-type scheme for the simulation of material boundaries between immiscible fluids with arbitrary equation of state is developed and validated. The detailed simulation model accounts for surface-tension, viscous, and body-force effects, in addition to acoustic and convective transport. The material interfaces are considered as diffused zones and a mixture model is given for this transition region. The simulation methodology combines a high-resolution discontinuity capturing method with a low-dissipation central scheme resulting in a hybrid approach for the solution of time- and space-accurate interface problems. Several multi-dimensional test cases are considered over a wide range of physical situations involving capillary, viscosity, and gravity effects with simultaneous presence of large viscosity and density ratios. The model is shown to accurately capture interface dynamics as well as to deal with dynamic appearance and disappearance of material boundaries. Simulation of atomization processes and its interaction with the flow field in practical devices is the secondary objective of this study. Three modeling requirements are identified to perform Large-Eddy Simulation (LES) of spray combustion in engineering devices. In concurrence with these requirements, LES of an experimental liquid-fueled Lean Direct Injection (LDI) combustor is performed using a subgrid mixing and combustion model. This approach has no adjustable parameters and the entire flow-path through the inlet swirl vanes is resolved. The inclusion of the atomization aspects within LES eliminates the need to specify dispersed-phase size-velocity correlations at the inflow boundary. Kelvin-Helmholtz (or aerodynamic) breakup model by Reitz is adopted for the combustor simulation. Two simulations (with and without breakup) are performed and compared with measurements of Cai et al. Time-averaged velocity prediction comparison for both gas- and liquid-phase with available data show reasonable agreement. The major impact of breakup is on the fuel evaporation in the vicinity of the injector. Further downstream, a wide range of drop sizes are recovered by the breakup simulation and produces similar spray quality as in the no-breakup case.