Memory-encoding vibrations in a disconnecting air bubble

Many nonlinear processes, such as the propagation of waves over an ocean or the transmission of light pulses down an optical fibre1, are integrable in the sense that the dynamics has as many conserved quantities as there are independent variables. The result is a time evolution that retains a complete memory of the initial state. In contrast, the nonlinear dynamics near a finite-time singularity, in which physical quantities such as pressure or velocity diverge at a point in time, is believed to evolve towards a universal form, one independent of the initial state2. The break-up of a water drop in air3 or a viscous liquid inside an immiscible oil4,5 are processes that conform to this second scenario. These opposing scenarios collide in the nonlinearity produced by the formation of a finite-time singularity that is also integrable. We demonstrate here that the result is a novel dynamics with a dual character

The g-mode Excitation in the Proto Neutron Star by the Standing Accretion Shock Instability

The so-called "acoustic revival mechanism" of core-collapse supernova proposed recently by the Arizona group is an interesting new possibility. Aiming to understand the elementary processes involved in the mechanism, we have calculated the eigen frequencies and eigen functions for the g-mode oscillations of a non-rotating proto neutron star. The possible excitation of these modes by the standing accretion shock instability, or SASI, is discussed based on these eigen functions. We have formulated the forced oscillations of $g$-modes by the external pressure perturbations exerted on the proto neutron star surface. The driving pressure fluctuations have been adopted from our previous computations of the axisymmetric SASI in the non-linear regime. We have paid particular attention to low l modes, since these are the modes that are dominant in SASI and that the Arizona group claimed played an important role in their acoustic revival scenario. Here l is the index of the spherical harmonic functions, $Y_l^m$. Although the frequency spectrum of the non-linear SASI is broadened substantially by non-linear couplings, the typical frequency is still much smaller than those of g-modes, the fact leading to a severe impedance mismatch. As a result, the excitations of various $g$-modes are rather inefficient and the energy of the saturated g-modes is $\sim 10^{50}$erg or smaller, with the g_2-mode being the largest in our model. Here the g_2-mode has two radial nodes and is confined to the interior of the convection region. The energy transfer rate from the g-modes to out-going sound waves is estimated from the growth of the g-modes and found to be $\sim 10^{51}$erg/s in the model studied in this paper.Comment: 24 pages, 6 figure

A compressible Lagrangian framework for the simulation of underwater implosion problems

The development of efficient algorithms to understand implosion dynamics presents a number of challenges. The foremost challenge is to efficiently represent the coupled compressible fluid dynamics of internal air and surrounding water. Secondly, the method must allow one to accurately detect or follow the interface between the phases. Finally, it must be capable of resolving any shock waves which may be created in air or water during the final stage of the collapse. We present a fully Lagrangian compressible numerical framework for the simulation of underwater implosion. Both air and water are considered compressible and the equations for the Lagrangian shock hydrodynamics are stabilized via a variationally consistent multiscale method. A nodally perfect matched definition of the interface is used and then the kinetic variables, pressure and density, are duplicated at the interface level. An adaptive mesh generation procedure, which respects the interface connectivities, is applied to provide enough refinement at the interface level. This framework is then used to simulate the underwater implosion of a large cylindrical bubble, with a size in the order of cm. Rapid collapse and growth of the bubble occurred on very small spatial (0.3mm), and time (0.1ms) scales followed by Rayleigh-Taylor instabilities at the interface, in addition to the shock waves traveling in the fluid domains are among the phenomena that are observed in the simulation. We then extend our framework to model the underwater implosion of a cylindrical aluminum container considering a monolithic fluid-structure interaction (FSI). The aluminum cylinder, which separates the internal atmospheric-pressure air from the external high-pressure water, is modeled by a three node rotation-free shell element. The cylinder undergoes fast transient deformations, large enough to produce self-contact along it. A novel elastic frictionless contact model is used to detect contact and compute the non-penetrating forces in the discretized domain between the mid-planes of the shell. Two schemes are tested, implicit using the predictor/multi-corrector Bossak scheme, and explicit, using the forward Euler scheme. The results of the two simulations are compared with experimental data

A compressible Lagrangian framework for the simulation of underwater implosion problems

The development of efficient algorithms to understand implosion dynamics presents a number of challenges. The foremost challenge is to efficiently represent the coupled compressible fluid dynamics of internal air and surrounding water. Secondly, the method must allow one to accurately detect or follow the interface between the phases. Finally, it must be capable of resolving any shock waves which may be created in air or water during the final stage of the collapse. We present a fully Lagrangian compressible numerical framework for the simulation of underwater implosion. Both air and water are considered compressible and the equations for the Lagrangian shock hydrodynamics are stabilized via a variationally consistent multiscale method [109]. A nodally perfect matched definition of the interface is used [57, 25] and then the kinetic variables, pressure and density, are duplicated at the interface level. An adaptive mesh generation procedure, which respects the interface connectivities, is applied to provide enough refinement at the interface level. This framework is then used to simulate the underwater implosion of a large cylindrical bubble, with a size in the order of cm. Rapid collapse and growth of the bubble occurred on very small spatial (0.3mm), and time (0.1ms) scales followed by Rayleigh-Taylor instabilities at the interface, in addition to the shock waves traveling in the fluid domains are among the phenomena that are observed in the simulation. We then extend our framework to model the underwater implosion of a cylindrical aluminum container considering a monolithic fluid-structure interaction (FSI). The aluminum cylinder, which separates the internal atmospheric-pressure air from the external high-pressure water, is modeled by a three node rotation-free shell element. The cylinder undergoes fast transient deformations, large enough to produce self-contact along it. A novel elastic frictionless contact model is used to detect contact and compute the non-penetrating forces in the discretized domain between the mid-planes of the shell. Two schemes are tested, implicit using the predictor/multi-corrector Bossak scheme, and explicit, using the forward Euler scheme. The results of the two simulations are compared with experimental data.El desarrollo de métodos eficientes para modelar la dinámica de implosión presenta varios desafíos. El primero es una representación eficaz de la dinámica del sistema acoplado de aire-agua. El segundo es que el método tiene que permitir una detección exacta o un seguimiento adecuado de la interfase entre ambas fases. Por último el método tiene que ser capaz de resolver cualquier choque que podría generar en el aire o en el agua, sobre todo en la última fase del colapso. Nosotros presentamos un método numérico compresible y totalmente Lagrangiano para simular la implosión bajo el agua. Tanto el aire como el agua se consideran compresibles y las ecuaciones Lagrangianos para la hidrodinámica del choque se estabilizan mediante un método multiescala que es variacionalmente consistente [109]. Se utiliza una definición de interfase que coincide perfectamente con los nodos [57, 25]. Ésta, nos facilita duplicar eficazmente las variables cinéticas como la presión y la densidad en los nodos de la interfase. Con el fin de obtener suficiente resolución alrededor de la interfase, la malla se genera de forma adaptativa y respetando la posición de la interfase. A continuación el método desarrollado se utiliza para simular la implosión bajo el agua de una burbuja cilíndrica del tamaño de un centímetro. Varios fenómenos se han capturado durante el colapso: un ciclo inmediato de colapso-crecimiento de la burbuja que ocurre en un espacio (0.3mm) y tiempo (0.1ms) bastante limitado, aparición de inestabilidades de tipo Rayleigh-Taylor en la interfase y formaron de varias ondas de choque que viajan tanto en el agua como en el aire. Después, seguimos el desarrollo del método para modelar la implosión bajo el agua de un contenedor metálico considerando una interacción monolítica de fluido y estructura. El cilindro de aluminio, que a su vez contiene aire a presión atmosférica y está rodeada de agua en alta presión, se modelando con elementos de lámina de tres nodos y sin grados de libertad de rotación. El cilindro se somete a deformaciones transitorias suficientemente rápidos y enormes hasta llegar a colapsar. Un nuevo modelo elástico de contacto sin considerar la fricción se ha desarrollado para detectar el contacto y calcular las fuerzas en el dominio discretizado entre las superficies medianas de las laminas. Dos esquemas temporales están considerados, uno es implícito utilizando el método de Bossak y otro es explícito utilizando Forward Euler. Al final los resultados de ambos casos se comparan con los resultados experimentales

Magneto-Rayleigh-Taylor Instability: Theory and Simulation in Planar and Cylindrical Pulsed Power Targets.

Cylindrical liner implosions in the Magnetized Liner Inertial Fusion (MagLIF) concept are susceptible to the magneto-Rayleigh-Taylor instability (MRT). The danger of MRT enters in two phases, (1) during the main implosion, the outer surface of the liner is MRT unstable, and (2) during the short time period when the liner decelerates onto the hot fuel, the inner surface becomes unstable. Growth of MRT on the outer surface may also feedthrough, which may seed the inner surface leading to high MRT growth in the second phase. If MRT growth becomes large enough, confinement of the fuel is lost. To characterize MRT, we solve the linearized, ideal MHD equations in both planar and cylindrical geometries, including an axial magnetic field and the effects of sausage and kink modes. To evaluate our analytic growth rates, 1D HYDRA MHD simulations are used to generate realistic, evolving profiles (in density, pressure, and magnetic field) during the implosion. In general, the total instability growth rates in cylindrical geometry are larger than those in planar geometry. MRT and feedthrough are suppressed by strong magnetic field line bending (tension). We apply our analytic MRT growth rates to experiments on the Z-machine at Sandia National Laboratories. Analytic MRT growth rates for a typical magnetized MagLIF-like implosion show the kink mode to be the fastest growing early and very late in the liner implosion (during deceleration). Sophisticated 2D HYDRA simulations show that highly compressed axial magnetic fields can reduce the growth of perturbations at the fuel/liner interface during the implosion phase, enhancing the stability of the implosion. HYDRA 2D simulations also show that a non-uniform shock, driven from the liner exterior, can seed the liner interior, leading to substantial growth of instability far in excess of feedthrough. Large-scale perturbations on the liner interior may also feedout to the liner exterior when a shock wave interacts with the surface, which further destabilizes the liner. These effects are reduced when shock compression is minimized or significant perturbations are not present during shock compression. The feedthrough effects then dominate.PhDNuclear ScienceUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/113414/1/weisy_1.pd