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

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 posteriori error estimation and modeling of unsaturated flow in fractured porous media

This doctoral thesis focuses on three topics: (1) modeling of unsaturated flow in fractured porous media, (2) a posteriori error estimation for mixed-dimensional elliptic equations, and (3) contributions to open-source software for complex multiphysics processes in porous media. In our first contribution, following a Discrete-Fracture Matrix (DFM) approach, we propose a model where Richards' equation governs the water flow in the matrix, whereas fractures are represented as lower-dimensional open channels, naturally providing a capillary barrier to the water flow. Therefore, water in the matrix is only allowed to imbibe the fracture if the capillary barrier is overcome. When this occurs, we assume that the water inside the fracture flows downwards without resistance and, therefore, is instantaneously at hydrostatic equilibrium. This assumption can be justifiable for fractures with sufficiently large apertures, where capillary forces play no role. Mathematically, our model can be classified as a coupled PDE-ODE system of equations with variational inequalities, in which each fracture is considered a potential seepage face. Our second contribution deals with error estimation for mixed-dimensional (mD) elliptic equations, which, in particular, model single-phase flow in fractured porous media. Here, based on the theory of functional a posteriori error estimates, we derive guaranteed upper bounds for the mD primal and mD dual variables, and two-sided bounds for the mD primal-dual pair. Moreover, we improve the standard results of the functional approach by proposing four ways of estimating the residual errors based on the conservation properties of the approximations, that is, (1) no conservation, (2) subdomain conservation, (3) local conservation, and (4) pointwise conservation. This results in sharper and fully-computable bounds when mass is conserved either locally or exactly. To our knowledge, to date, no error estimates have been available for fracture networks, including fracture intersections and floating subdomains. Our last contribution is related to the development of open-source software. First, we present the implementation of a new multipoint finite-volume-based module for unsaturated poroelasticity, compatible with the Matlab Reservoir Simulation Toolbox (MRST). Second, we present a new Python-based simulation framework for multiphysics processes in fractured porous media, named PorePy. PorePy, by design, is particularly well-suited for handling mixed-dimensional geometries, and thus optimal for DFM models. The first two contributions discussed above were implemented in PorePy.Denne avhandlingen tar for seg tre emner: (1) modellering av flyt i umettet porøst medium med sprekker, (2) a posteriori feilestimater for blandet-dimensjonale elliptiske ligninger, og (3) bidrag til åpen kildekode for komplekse multifysikk-prosesser i porøse medier. I det første bidraget anvender vi en Discrete-Fracture Matrix (DFM) (Diskret-Sprekk Matrise) metode til å sette opp en modell hvor Richard's ligning modellerer vann-flyt i matrisen, og sprekkene representeres som lavere-dimensjonale åpne kanaler, som naturlig virker som kapillærbarrierer til vann-flyten. Derfor vil vann i matrisen kun få tilgang til sprekken når kapillærbarrieren blir brutt. Når det inntreffer, antar vi at vannet i sprekken flyter nedover uten motstand, og at hydrostatisk ekvilibrium derfor inntreffer øyeblikkelig. Slike antakelser kan rettferdiggjøres for sprekker med tilstrekkelig stor apertur (åpning), hvor kapillærkrefter ikke har noen innvirkning. Fra et matematisk standpunkt kan modellen klassifiseres som en sammenkoblet PDE-ODE med variasjonelle ulikheter hvor hver sprekk behandles som en filtreringsfase. Det andre bidraget tar for seg feilestimater for blandet-dimensjonale elliptiske ligninger, som modellerer en-fase flyt i porøse medier med sprekker. Her anvender vi teorien for "funksjonal a posteriori feilestimater" til å finne øvre skranker for primær og dual variablene, samt øvre og nedre skranker for primær-dual paret. Dessuten viser vi at vi kan forbedre standardresultatene fra "funksjonal a posteriori feilestimater" ved å foreslå fire måte å estimere residualfeilen basert på bevaringsegenskapene til diskretiseringen. De fire forskjellige bevaringsegenskapene er; ingen bevaringsegenskap, under- domene bevaring, lokal bevaring og punktvis bevaring. Dette fører til skarpere skranker som er mulige å beregne når masse er bevart enten lokalt, eller eksakt. Vi kjenner ikke til andre tilgjengelige feilestimater for sprekknettverk som inkluderer snitt av sprekker og sprekkrender som ligger innenfor domenets rand. Det siste bidraget omhandler utvikling av åpen kildekode. Først presenterer vi imple- menteringen av en multipunktfluks-basert modul for flyt i umettet deformerbart porøst medium som er kompatibelt med "Matlab Reservoir Simulation Toolbox"(MRST). I tillegg presenterer vi et nytt Python-basert rammeverk for simulering av multifysikkprosesser i porøse medier med sprekker, som heter PorePy. Dette rammeverket er designet for å håndtere geometrier med blandede dimensjoner og er derfor optimalt for DFM modeller. De to første bidragene i avhandlingen (nevnt over) er implementert i PorePy.Doktorgradsavhandlin

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described

Software for Exascale Computing - SPPEXA 2016-2019

This open access book summarizes the research done and results obtained in the second funding phase of the Priority Program 1648 "Software for Exascale Computing" (SPPEXA) of the German Research Foundation (DFG) presented at the SPPEXA Symposium in Dresden during October 21-23, 2019. In that respect, it both represents a continuation of Vol. 113 in Springer’s series Lecture Notes in Computational Science and Engineering, the corresponding report of SPPEXA’s first funding phase, and provides an overview of SPPEXA’s contributions towards exascale computing in today's sumpercomputer technology. The individual chapters address one or more of the research directions (1) computational algorithms, (2) system software, (3) application software, (4) data management and exploration, (5) programming, and (6) software tools. The book has an interdisciplinary appeal: scholars from computational sub-fields in computer science, mathematics, physics, or engineering will find it of particular interest