Numerical Simulations of Shock and Rarefaction Waves Interacting With Interfaces in Compressible Multiphase Flows

Developing a highly accurate numerical framework to study multiphase mixing in high speed flows containing shear layers, shocks, and strong accelerations is critical to many scientific and engineering endeavors. These flows occur across a wide range of scales: from tiny bubbles in human tissue to massive stars collapsing. The lack of understanding of these flows has impeded the success of many engineering applications, our comprehension of astrophysical and planetary formation processes, and the development of biomedical technologies. Controlling mixing between different fluids is central to achieving fusion energy, where mixing is undesirable, and supersonic combustion, where enhanced mixing is important. Iron, found throughout the universe and a necessary component for life, is dispersed through the mixing processes of a dying star. Non-invasive treatments using ultrasound to induce bubble collapse in tissue are being developed to destroy tumors or deliver genes to specific cells. Laboratory experiments of these flows are challenging because the initial conditions and material properties are difficult to control, modern diagnostics are unable to resolve the flow dynamics and conditions, and experiments of these flows are expensive. Numerical simulations can circumvent these difficulties and, therefore, have become a necessary component of any scientific challenge. Advances in the three fields of numerical methods, high performance computing, and multiphase flow modeling are presented: (i) novel numerical methods to capture accurately the multiphase nature of the problem; (ii) modern high performance computing paradigms to resolve the disparate time and length scales of the physical processes; (iii) new insights and models of the dynamics of multiphase flows, including mixing through hydrodynamic instabilities. These studies have direct applications to engineering and biomedical fields such as fuel injection problems, plasma deposition, cancer treatments, and turbomachinery.PhDMechanical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/133458/1/marchdf_1.pd

Reactive Flows in Deformable, Complex Media

Many processes of highest actuality in the real life are described through systems of equations posed in complex domains. Of particular interest is the situation when the domain is changing in time, undergoing deformations that depend on the unknown quantities of the model. Such kind of problems are encountered as mathematical models in the subsurface, material science, or biological systems.The emerging mathematical models account for various processes at different scales, and the key issue is to integrate the domain deformation in the multi-scale context. The focus in this workshop was on novel techniques and ideas in the mathematical modelling, analysis, the numerical discretization and the upscaling of problems as described above

Moment Methods for Advection on Networks and an Application to Forest Pest Life Cycle Models

This paper develops low-dimensional moment methods for advective problems on networks of domains. The evolution of a density function is described by a linear advection-diffusion-reaction equation on each domain, combined via advective flux coupling across domains in the network graph. The PDEs' coefficients vary in time and across domains but they are fixed along each domain. As a result, the solution on each domain is frequently close to a Gaussian that moves, decays, and widens. For that reason, this work studies moment methods that track only three degrees of freedom per domain -- in contrast to traditional PDE discretization methods that tend to require many more variables per domain. A simple ODE-based moment method is developed, as well as an asymptotic-preserving scheme. We apply the methodology to an application that models the life cycle of forest pests that undergo different life stages and developmental pathways. The model is calibrated for the spotted lanternfly, an invasive species present in the Eastern USA. We showcase that the moment method, despite its significant low-dimensionality, can successfully reproduce the prediction of the pest's establishment potential, compared to much higher-dimensional computational approaches.Comment: 31 pages, 14 figure

Parallel and Adaptive Galerkin Methods for Radiative Transfer Problems

Finite element methods for radiative transfer problems are developed and their implementation on parallel computer architectures is discussed. A posteriori error estimates allow for efficient grid adaptation in two and three space dimensions