Engineering 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 …

Investigating The Feasibility And Stability For Modeling Acoustic Wave Scattering Using A Time-Domain Boundary Integral Equation With Impedance Boundary Condition, Michelle E. Rodio

Mathematics & Statistics Theses & Dissertations

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet and environmentally friendly aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction and are characterized by a frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners

This work considers using either an impedance or an admittance (inverse …

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 …

A Least Squares Closure Approximation For Liquid Crystalline Polymers, Traci Ann Sievenpiper

Mathematics & Statistics Theses & Dissertations

An introduction to existing closure schemes for the Doi-Hess kinetic theory of liquid crystalline polymers is provided. A new closure scheme is devised based on a least squares fit of a linear combination of the Doi, Tsuji-Rey, Hinch-Leal I, and Hinch-Leal II closure schemes. The orientation tensor and rate-of-strain tensor are fit separately using data generated from the kinetic solution of the Smoluchowski equation. The known behavior of the kinetic solution and existing closure schemes at equilibrium is compared with that of the new closure scheme. The performance of the proposed closure scheme in simple shear flow for a variety …

Diffusion Problems In Wound Healing And A Scattering Approach To Immune System Interactions, Julia Suzanne Arnold

Mathematics & Statistics Theses & Dissertations

A theoretical model for the existence of a Critical Size Defect (CSD) in certain animals is the focus of the majority of this dissertation. Adam [1] recently developed a one-dimensional model of this phenomenon, and chapters I–V address the exist the CSD in a two-dimensional model and a three-dimensional model. The two dimensional (or 1-d circular) model is the more appropriate for a study of CSD's. In that model we assume a circular wound of uniform depth and develop a time-independent form of the diffusion equation relevant to the study of the CSD phenomenon. It transpires that the range of …

The Solution Of Hypersingular Integral Equations With Applications In Acoustics And Fracture Mechanics, Richard S. St. John

Mathematics & Statistics Theses & Dissertations

The numerical solution of two classes of hypersingular integral equations is addressed. Both classes are integral equations of the first kind, and are hypersingular due to a kernel containing a Hadamard singularity. The convergence of a Galerkin method and a collocation method is discussed and computationally efficient algorithms are developed for each class of hypersingular integral equation.

Interest in these classes of hypersingular integral equations is due to their occurrence in many physical applications. In particular, investigations into the scattering of acoustic waves by moving objects and the study of dynamic Griffith crack problems has necessitated a computationally efficient technique …

High-Order Finite-Difference Schemes And Their Application To Computational Acoustics, Joe Leo Manthey

Mathematics & Statistics Theses & Dissertations

The primary focus of this study is upon the numerical stability of high-order finite-difference schemes and their application to duct acoustics. Since acoustic waves are known to be non-dissipative and non-dispersive, high-order schemes are favored for their low dissipation and low dispersion relative to the low-order schemes. The primary obstacle to the the development of explicit high-order finite-difference schemes is the construction of boundary closures which simultaneously maintain the formal order of accuracy and the numerical stability of the overall scheme. In this thesis a hybrid seven-point, fourth-order stencil for computing spatial derivatives is presented and the time stability is …

Reverse Engineering Of Aircraft Wing Data Using A Partial Differential Equation Surface Model, Jacalyn M. Huband

Mathematics & Statistics Theses & Dissertations

Reverse engineering is a multi-step process used in industry to determine a production representation of an existing physical object. This representation is in the form of mathematical equations that are compatible with computer-aided design and computer-aided manufacturing (CAD/CAM) equipment. The four basic steps to the reverse engineering process are data acquisition, data separation, surface or curve fitting, and CAD/CAM production. The surface fitting step determines the design representation of the object, and thus is critical to the success or failure of the reverse engineering process. Although surface fitting methods described in the literature are used to model a variety of …

Exact Solutions For Orthogonal And Non-Orthogonal Magnetohydrodynamic Stagnation-Point Flow, Shahrooz Moosavizadeh

Mathematics & Statistics Theses & Dissertations

The viscous plane flow of an electrically conducting fluid towards an infinite wall is solved in the presence of a magnetic field which is aligned with the flow far from the wall. The problem has two dimensionless parameters-- ε, the magnetic Prandtl number, and β, the square of the ratio of Alfven velocity to fluid velocity far from the wall. The problem has a similarity solution which reduces the governing equations to a system of coupled ordinary differential equations which can be solved numerically. For extreme values of ε, both large and small, singular perturbation techniques are used to derive …

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 …

Nozzle Flow With Vibrational Nonequilibrium, John Gary Landry

Mathematics & Statistics Theses & Dissertations

Flow of nitrogen gas through a converging-diverging nozzle is simulated. The flow is modeled using the Navier-Stokes equations that have been modified for vibrational nonequilibrium. The energy equation is replaced by two equations. One equation accounts for energy effects due to the translational and rotational degrees of freedom, and the other accounts for the affects due to the vibrational degree of freedom. The energy equations are coupled by a relaxation time which measures the time required for the vibrational energy component to equilibrate with the translational and rotational energy components. An improved relaxation time is used in this thesis. The …

Boundary Value Problems In Rectilinearly Anisotropic Thermoelastic Solids, Gilbert Kerr

Mathematics & Statistics Theses & Dissertations

The boundary value problems which are considered are the type that arise due to the presence of a Griffith crack (or cracks) in an anisotropic thermoelastic solid. The thermoelastic field, in such materials, when the infinitesimal theory is employed, is governed by a set of elliptic partial differential equations. The general solution of these equations is expressed in terms of arbitrary analytic functions whose real parts, in turn, are expressed in terms of Fourier type integrals or Fourier series. Integral transform techniques are then used to determine the stress intensity factors (and other pertinent information) for various crack geometries. In …

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 …

The Solution Of A Singular Integral Equation Arising From A Lifting Surface Theory For Rotating Blades, Mark H. Dunn

Mathematics & Statistics Theses & Dissertations

A technique is presented for the solution of a linear, two dimensional, singular, Volterra integral equation of the first kind. The integral equation, originally developed by Farassat and Myers, is derived from the basic equations of linearized acoustics and models the lifting force experienced by an infinitesimally thin surface moving tangent to itself. As a particular application, the motion of modern high speed aircraft propellers (Advanced Technology Propellers) is considered. The unknown propeller blade surface pressure distribution is approximated by a piecewise constant function and the integral equation is solved numerically by the method of collocation. Certain simplifying assumptions applied …

An Extension Of Essentially Non-Oscillatory Shock-Capturing Schemes To Multi-Dimensional Systems Of Conservation Laws, Jay Casper

Mathematics & Statistics Theses & Dissertations

In recent years, a class of numerical schemes for solving hyperbolic partial differential equations has been developed which generalizes the first-order method of Godunov to arbitrary order of accuracy. High-order accuracy is obtained, wherever the solution is smooth, by an essentially non-oscillatory (ENO) piecewise polynomial reconstruction procedure, which yields high-order pointwise information from the cell averages of the solution at a given point in time. When applied to piecewise smooth initial data, this reconstruction enables a flux computation that provides a time update of the solution which is of high-order accuracy, wherever the function is smooth, and avoids a Gibbs …