Physics Commons^™

High-Order Positivity-Preserving L2-Stable Spectral Collocation Schemes For The 3-D Compressible Navier-Stokes Equations, Johnathon Keith Upperman

Mathematics & Statistics Theses & Dissertations

High-order entropy stable schemes are a popular method used in simulations with the compressible Euler and Navier-Stokes equations. The strength of these methods is that they formally satisfy a discrete entropy inequality which can be used to guarantee L₂ stability of the numerical solution. However, a fundamental assumption that is explicitly or implicitly used in all entropy stability proofs available in the literature for the compressible Euler and Navier-Stokes equations is that the thermodynamic variables (e.g., density and temperature) are strictly positive in the entire space{time domain considered. Without this assumption, any entropy stability proof for a numerical scheme …

Electrohydrodynamic Simulations Of Capsule Deformation Using A Dual Time-Stepping Lattice Boltzmann Scheme, Charles Leland Armstrong

Mathematics & Statistics Theses & Dissertations

Capsules are fluid-filled, elastic membranes that serve as a useful model for synthetic and biological membranes. One prominent application of capsules is their use in modeling the response of red blood cells to external forces. These models can be used to study the cell’s material properties and can also assist in the development of diagnostic equipment. In this work we develop a three dimensional model for numerical simulations of red blood cells under the combined influence of hydrodynamic and electrical forces. The red blood cell is modeled as a biconcave-shaped capsule suspended in an ambient fluid domain. Cell deformation occurs …

Topics In Electromagnetic, Acoustic, And Potential Scattering Theory, Umaporn Nuntaplook

Mathematics & Statistics Theses & Dissertations

With recent renewed interest in the classical topics of both acoustic and electromagnetic aspects for nano-technology, transformation optics, fiber optics, metamaterials with negative refractive indices, cloaking and invisibility, the topic of time-independent scattering theory in quantum mechanics is becoming a useful field to re-examine in the above contexts. One of the key areas of electromagnetic theory scattering of plane electromagnetic waves — is based on the properties of the refractive indices in the various media. It transpires that the refractive index of a medium and the potential in quantum scattering theory are intimately related. In many cases, understanding such scattering …

Perfectly Matched Layer Absorbing Boundary Conditions For The Discrete Velocity Boltzmann-Bgk Equation, Elena Craig

Mathematics & Statistics Theses & Dissertations

Perfectly Matched Layer (PML) absorbing boundary conditions were first proposed by Berenger in 1994 for the Maxwell's equations of electromagnetics. Since Hu first applied the method to Euler's equations in 1996, progress made in the application of PML to Computational Aeroacoustics (CAA) includes linearized Euler equations with non-uniform mean flow, non-linear Euler equations, flows with an arbitrary mean flow direction, and non-linear clavier-Stokes equations. Although Boltzmann-BGK methods have appeared in the literature and have been shown capable of simulating aeroacoustics phenomena, very little has been done to develop absorbing boundary conditions for these methods. The purpose of this work was …

The Straggling Green's Function Method For Ion Transport, Steven Andrew Walker

Mathematics & Statistics Theses & Dissertations

For many years work has been conducted on developing a concise theory and method for HZE ion transport capable of being validated in the laboratory. Previous attempts have ignored dispersion and energy downshift associated with nuclear fragmentation and energy and range straggling. Here we present a Green's function approach to ion transport that incorporates these missing elements. This work forms the basis for a new version of GRNTRN, a Green's function transport code. Comparisons of GRNTRN predictions and laboratory results for an ⁵⁶Fe ion beam with average energy at the target of one GeV/amu or more are presented for …

Mathematical Models Of Quiescent Solar Prominences, Iain Mckaig

Mathematics & Statistics Theses & Dissertations

Magnetic fields in the solar atmosphere suspend and insulate dense regions of cool plasma known as prominences. The convection zone may be the mechanism that both generates and expels this magnetic flux through the photosphere in order to make these formations possible. The connection is examined here by modeling the convection zone as both one-dimensional, then more realistically, two-dimensional.

First a Dirichlet problem on a semi-infinite strip is solved using conformal mapping and the method of images. The base of the strip represents the photosphere where a current distribution can be given as a boundary condition, and the strip extends …

Studies Of Mixing Processes In Gases And Effects On Combustion And Stability, Frank Paul Kozusko Jr.

Mathematics & Statistics Theses & Dissertations

Three physical models of laminar mixing of initially separated gases are studied. Two models study the effects of the mixing dynamics on the chemical reactions between the gases. The third model studies the structure and stability of a laminar mixing layer in a binary gas. The three models are:

1. Two ideal and incompressible gases representing fuel and oxidizer are initially at rest and separated across an infinite linear interface in a two dimensional system. Combustion, expected as the gases mix, will lead to a rapid rise in temperature in a localized area, i.e. ignition. The mixing of the gases …

On Shock Capturing For Liquid And Gas Media, Tze Jang Chen

Mathematics & Statistics Theses & Dissertations

The numerical investigation of shock phenomena in gas or liquid media where a specifying relation for internal energy is absent poses special problems. Classically, for gas dynamics the usual procedure is to employ a splitting scheme to remove the source terms from the Euler equations, then up-wind biased shock capturing algorithms are built around the Riemann problem for the system which remains. However, in the case where the Euler equations are formulated in the term of total enthalpy, a technical difficulty associated with equation splitting forces a pressure time derivative to be treated as a source term. This makes it …

A Mathematical Model Of The Dynamics Of An Optically Pumped Codoped Solid State Laser System, Thomas G. Wangler

Mathematics & Statistics Theses & Dissertations

This is a study of a mathematical model for the dynamics of an optically pumped codoped solid state laser system. The model comprises five first order, nonlinear, coupled, ordinary differential equations which describe the temporal evolution of the dopant electron populations in the laser crystal as well as the photon density in the laser cavity. The analysis of the model is conducted in three parts.

First, a detailed explanation of the modeling process is given and the full set of rate equations is obtained. The model is then simplified and certain qualitative properties of the solution are obtained.

In the …

On A Moving Boundary Problem Of Transitional Ballistics, Jen-Ing G. Hwang

Mathematics & Statistics Theses & Dissertations

A major problem which arises in computer simulation of the firing of a gun weapon is the development of numerical schemes which effectively account for the physics of projectile motion. The chief difficulty is that away from the projectile the calculation is ordinarily accomplished on a fixed numerical grid, whereas due to projectile movement some cells of the grid near the projectile undergo volume changes as the calculation proceeds. A local finite volume scheme is developed which accounts for the expansion or compression of cells fore-and-aft of the projectile. Through the process of numerical experiment, the effectiveness of the scheme …

Minimal Norm Constrained Interpolation, Larry Dean Irvine

Mathematics & Statistics Theses & Dissertations

In computational fluid dynamics and in CAD/CAM a physical boundary, usually known only discreetly (say, from a set of measurements), must often be approximated. An acceptable approximation must, of course, preserve the salient features of the data (convexity, concavity, etc.) In this dissertation we compute a smooth interpolant which is locally convex where the data are locally convex and is locally concave where the data are locally concave.

Such an interpolant is found by posing and solving a minimization problem. The solution is a piecewise cubic polynomial. We actually solve this problem indirectly by using the Peano kernel theorem to …