Physical Sciences and Mathematics Commons^™

Reducing Communication In The Solution Of Linear Systems, Neil S. Lindquist

Doctoral Dissertations

There is a growing performance gap between computation and communication on modern computers, making it crucial to develop algorithms with lower latency and bandwidth requirements. Because systems of linear equations are important for numerous scientific and engineering applications, I have studied several approaches for reducing communication in those problems. First, I developed optimizations to dense LU with partial pivoting, which downstream applications can adopt with little to no effort. Second, I consider two techniques to completely replace pivoting in dense LU, which can provide significantly higher speedups, albeit without the same numerical guarantees as partial pivoting. One technique uses randomized …

Space-Angle Discontinuous Galerkin Finite Element Method For Radiative Transfer Equation, Hang Wang

Doctoral Dissertations

Radiative transfer theory describes the interaction of radiation with scattering and absorbing media. It has applications in neutron transport, atmospheric physics, heat transfer, molecular imaging, and others. In steady state, the radiative transfer equation is an integro-differential equation of five independent variables, which are 3 dimensions in space and 2 dimensions in the angular direction. This high dimensionality and the presence of the integral term present serious challenges when solving the equation numerically. Over the past 50 years, several techniques for solving the radiative transfer equation (RTE) have been introduced. These include, but are certainly not limited to, Monte Carlo …

Towards Reduced-Order Model Accelerated Optimization For Aerodynamic Design, Andrew L. Kaminsky

Doctoral Dissertations

The adoption of mathematically formal simulation-based optimization approaches within aerodynamic design depends upon a delicate balance of affordability and accessibility. Techniques are needed to accelerate the simulation-based optimization process, but they must remain approachable enough for the implementation time to not eliminate the cost savings or act as a barrier to adoption.

This dissertation introduces a reduced-order model technique for accelerating fixed-point iterative solvers (e.g. such as those employed to solve primal equations, sensitivity equations, design equations, and their combination). The reduced-order model-based acceleration technique collects snapshots of early iteration (pre-convergent) solutions and residuals and then uses them to project …

Model Based Force Estimation And Stiffness Control For Continuum Robots, Vincent A. Aloi

Doctoral Dissertations

Continuum Robots are bio-inspired structures that mimic the motion of snakes, elephant trunks, octopus tentacles, etc. With good design, these robots can be naturally compliant and miniaturizable, which makes Continuum Robots ideal for traversing narrow complex environments. Their flexible design, however, prevents us from using traditional methods for controlling and estimating loading on rigid link robots.

In the first thrust of this research, we provided a novel stiffness control law that alters the behavior of an end effector during contact. This controller is applicable to any continuum robot where a method for sensing or estimating tip forces and pose exists. …

Numerical Methods For Stochastic Stokes And Navier-Stokes Equations, Liet Vo

Doctoral Dissertations

This dissertation consists of three main parts with each part focusing on numerical approximations of the stochastic Stokes and Navier-Stokes equations.

Part One concerns the mixed finite element methods and Chorin projection methods for solving the stochastic Stokes equations with general multiplicative noise. We propose a modified mixed finite element method for solving the Stokes equations and show that the numerical solutions converge optimally to the PDE solutions. The convergence is under energy norms (strong convergence) for the velocity and in a time-averaged norm (weak convergence) for the pressure. In addition, after establishing the error estimates in second moment, high …

A Weak Fractional Calculus Theory And Numerical Methods For Fractional Differential Equations, Mitchell D. Sutton

Doctoral Dissertations

This dissertation is comprised of four integral parts. The first part comprises a self-contained new theory of weak fractional differential calculus in one-dimension. The crux of this new theory is the introduction of a weak fractional derivative notion which is a natural generalization of integer order weak derivatives; it also helps to unify multiple existing fractional derivative definitions.

The second part of this work presents three new families of fractional Sobolev spaces and their accompanying theory in one-dimension. The new construction and theory are based on a newly developed notion of weak fractional derivatives, which are natural generalizations of the …

Preconditioned Nesterov’S Accelerated Gradient Descent Method And Its Applications To Nonlinear Pde, Jea Hyun Park

Doctoral Dissertations

We develop a theoretical foundation for the application of Nesterov’s accelerated gradient descent method (AGD) to the approximation of solutions of a wide class of partial differential equations (PDEs). This is achieved by proving the existence of an invariant set and exponential convergence rates when its preconditioned version (PAGD) is applied to minimize locally Lipschitz smooth, strongly convex objective functionals. We introduce a second-order ordinary differential equation (ODE) with a preconditioner built-in and show that PAGD is an explicit time-discretization of this ODE, which requires a natural time step restriction for energy stability. At the continuous time level, we show …

Examining The Accumulation Statistics Of Index1 Saddle Points On The Potential Energy Surface And Imposing Early Termination On A Rejection Scheme For Off Lattice Kinetic Monte Carlo, Jonathan W. Hicks

Doctoral Dissertations

In the calculation of time evolution of an atomic system where a chemical reaction and/or diffusion occurs, off-lattice kinetic Monte Carlo methods can be used to overcome timescale and lattice based limitations from other methods such as Molecular Dynamics and kinetic Monte Carlo procedures. Off-lattice kinetic Monte Carlo methods rely on a harmonic approximation to Transition State Theory, in which the rate of the rare transitions from one energy minimum to a neighboring minimum require surmounting a minimum energy barrier on the Potential Energy Surface, which is found at an index-1 saddle point commonly referred to as a transition state. …

Mathematical Modeling Of Mixtures And Numerical Solution With Applications To Polymer Physics, John Timothy Cummings

Doctoral Dissertations

We consider in this dissertation the mathematical modeling and simulation of a general diffuse interface mixture model based on the principles of energy dissipation. The model developed allows for a thermodynamically consistent description of systems with an arbitrary number of different components, each of which having perhaps differing densities. We also provide a mathematical description of processes which may allow components to source or sink into other components in a mass conserving, energy dissipating way, with the motivation of applying this model to phase transformation. Also included in the modeling is a unique set of thermodynamically consistent boundary conditions which …

Efficient Methods For Multidimensional Global Polynomial Approximation With Applications To Random Pdes, Peter A. Jantsch

Doctoral Dissertations

In this work, we consider several ways to overcome the challenges associated with polynomial approximation and integration of smooth functions depending on a large number of inputs. We are motivated by the problem of forward uncertainty quantification (UQ), whereby inputs to mathematical models are considered as random variables. With limited resources, finding more efficient and accurate ways to approximate the multidimensional solution to the UQ problem is of crucial importance, due to the “curse of dimensionality” and the cost of solving the underlying deterministic problem.

The first way we overcome the complexity issue is by exploiting the structure of the …

Numerical Methods For Non-Divergence Form Second Order Linear Elliptic Partial Differential Equations And Discontinuous Ritz Methods For Problems From The Calculus Of Variations, Stefan Raymond Schnake

Doctoral Dissertations

This dissertation consists of three integral parts. Part one studies discontinuous Galerkin approximations of a class of non-divergence form second order linear elliptic PDEs whose coefficients are only continuous. An interior penalty discontinuous Galerkin (IP-DG) method is developed for this class of PDEs. A complete analysis of the proposed IP-DG method is carried out, which includes proving the stability and error estimate in a discrete W^2;p-norm [W^2,p-norm]. Part one also studies the convergence of the vanishing moment method for this class of PDEs. The vanishing moment method refers to a PDE technique for approximating these PDEs by a …

Jet-Hadron Correlations Relative To The Event Plane Pb--Pb Collisions At The Lhc In Alice, Joel Anthony Mazer

Doctoral Dissertations

In relativistic heavy ion collisions at the Large Hadron Collider (LHC), a hot, dense and strongly interacting medium known as the Quark Gluon Plasma (QGP) is produced. Quarks and gluons from incoming nuclei collide to produce partons at high momenta early in the collisions. By fragmenting into collimated sprays of hadrons, these partons form 'jets'. Within the framework of perturbative Quantum Chromodynamics (pQCD), jet production is well understood in pp collisions. We can use jets measured in pp interactions as a baseline reference for comparing to heavy ion collision systems to detect and study jet quenching. The jet quenching mechanism …

Accuracy And Stability Of Integration Methods For Neutrino Transport In Core Collapse Supernovae, Kyle A. Gregory

Chancellor’s Honors Program Projects

No abstract provided.

Surface Energy In Bond-Counting Models On Bravais And Non-Bravais Lattices, Tim Ryan Krumwiede

Doctoral Dissertations

Continuum models in computational material science require the choice of a surface energy function, based on properties of the material of interest. This work shows how to use atomistic bond-counting models and crystal geometry to inform this choice. We will examine some of the difficulties that arise in the comparison between these models due to differing types of truncation. New crystal geometry methods are required when considering materials with non-Bravais lattice structure, resulting in a multi-valued surface energy. These methods will then be presented in the context of the two-dimensional material graphene in a way that correctly predicts its equilibrium …

Using Poincaré And Coefficient Analyses To Assess Changes In Variability In Respiration As A Function Of Leptin Status, Sex, And Buprenorphine In Mice, Wateen Hussein Alami

Chancellor’s Honors Program Projects

No abstract provided.

Kinetic Monte Carlo Models For Crystal Defects, Kyle Louis Golenbiewski

Doctoral Dissertations

Kinetic Monte Carlo algorithms have become an increasingly popular means to simulate stochastic processes since their inception in the 1960's. One area of particular interest is their use in simulations of crystal growth and evolution in which atoms are deposited on, or hop between, predefined lattice locations with rates depending on a crystal's configuration. Two such applications are heteroepitaxial thin films and grain boundary migration. Heteroepitaxial growth involves depositing one material onto another with a different lattice spacing. This misfit leads to long-range elastic stresses that affect the behavior of the film. Grain boundary migration, on the other hand, describes …

Computational Simulation Of Mass Transport And Energy Transfer In The Microbial Fuel Cell System, Shiqi Ou

Doctoral Dissertations

This doctoral dissertation introduces the research in the computational modeling and simulation for the microbial fuel cell (MFC) system which is a bio-electrochemical system that drives a current by using bacteria and mimicking bacterial interactions found in nature. The numerical methods, research approaches and simulation comparison with the experiments in the microbial fuel cells are described; the analysis and evaluation for the model methods and results that I have achieved are presented in this dissertation.

The development of the renewable energy has been a hot topic, and scientists have been focusing on the microbial fuel cell, which is an environmentally-friendly …

Domain Decomposition Methods For Discontinuous Galerkin Approximations Of Elliptic Problems, Craig Dwain Collins

Doctoral Dissertations

The application of the techniques of domain decomposition to construct effective preconditioners for systems generated by standard methods such as finite difference or finite element methods has been well-researched in the past few decades. However, results concerning the application of these techniques to systems created by the discontinuous Galerkin method (DG) are much more rare.

This dissertation represents the effort to extend the study of two-level nonoverlapping and overlapping additive Schwarz methods for DG discretizations of second- and fourth-order elliptic partial differential equations. In particular, the general Schwarz framework is used to find theoretical bounds for the condition numbers of …

Numerical Methods For Deterministic And Stochastic Phase Field Models Of Phase Transition And Related Geometric Flows, Yukun Li

Doctoral Dissertations

This dissertation consists of three integral parts with each part focusing on numerical approximations of several partial differential equations (PDEs). The goals of each part are to design, to analyze and to implement continuous or discontinuous Galerkin finite element methods for the underlying PDE problem.

Part One studies discontinuous Galerkin (DG) approximations of two phase field models, namely, the Allen-Cahn and Cahn-Hilliard equations, and their related curvature-driven geometric problems, namely, the mean curvature flow and the Hele-Shaw flow. We derive two discrete spectrum estimates, which play an important role in proving the sharper error estimates which only depend on a …

Flexible Memory Allocation In Kinetic Monte Carlo Simulations, Aaron David Craig

Masters Theses

We introduce two new algorithms for Kinetic Monte Carlo simulations: the minimal and flexible allocation algorithms. The theory and computational challenges associated with K.M.C. simulations are briefly discussed. We outline the simple cubic, solid-on-solid model of epitaxial growth and analyze four methods for its simulation: the linear search, standard inverted list, minimal allocation, and flexible allocation algorithms. We then implement these algorithms, analyze their performances, and discuss implications of the results.

Numerical Methods For Solving Optimal Control Problems, Garrett Robert Rose

Masters Theses

There are many different numerical processes for approximating an optimal control problem. Three of those are explained here: The Forward Backward Sweep, the Shooter Method, and an Optimization Method using the MATLAB Optimization Tool Box. The methods will be explained, and then applied to three different test problems to see how they perform. The results show that the Forward Backward Sweep is the best of the three methods with the Shooter Method being a competitor.

Numerical Analysis Of Convex Splitting Schemes For Cahn-Hilliard And Coupled Cahn-Hilliard-Fluid-Flow Equations, Amanda Emily Diegel

Doctoral Dissertations

This dissertation investigates numerical schemes for the Cahn-Hilliard equation and the Cahn-Hilliard equation coupled with a Darcy-Stokes flow. Considered independently, the Cahn-Hilliard equation is a model for spinodal decomposition and domain coarsening. When coupled with a Darcy-Stokes flow, the resulting system describes the flow of a very viscous block copolymer fluid. Challenges in creating numerical schemes for these equations arise due to the nonlinear nature and high derivative order of the Cahn-Hilliard equation. Further challenges arise during the coupling process as the coupling terms tend to be nonlinear as well. The numerical schemes presented herein preserve the energy dissipative structure …

Impacts Of Climate Change On The Evolution Of The Electrical Grid, Melissa Ree Allen

Doctoral Dissertations

Maintaining interdependent infrastructures exposed to a changing climate requires understanding 1) the local impact on power assets; 2) how the infrastructure will evolve as the demand for infrastructure changes location and volume and; 3) what vulnerabilities are introduced by these changing infrastructure topologies. This dissertation attempts to develop a methodology that will a) downscale the climate direct effect on the infrastructure; b) allow population to redistribute in response to increasing extreme events that will increase under climate impacts; and c) project new distributions of electricity demand in the mid-21^st century.

The research was structured in three parts. The first …

Statistical Mechanics And Schramm-Loewner Evolution With Applications To Crack Propagation Processes, Christopher Borut Mesic

Masters Theses

Schramm-Loewner Evolution (SLE) has both mathematical and physical roots that extend as far back as the early 20th century. We present the progression of these humble roots from the Ideal Gas Law, all the way to the renormalization group and conformal field theory, to better understand the impact SLE has had on modern statistical mechanics. We then explore the potential application of the percolation exploration process to crack propagation processes, illustrating the interplay between mathematics and physics.

Numerical Methods And Algorithms For High Frequency Wave Scattering Problems In Homogeneous And Random Media, Cody Samuel Lorton

Doctoral Dissertations

This dissertation consists of four integral parts with a unified objective of developing efficient numerical methods for high frequency time-harmonic wave equations defined on both homogeneous and random media. The first part investigates the generalized weak coercivity of the acoustic Helmholtz, elastic Helmholtz, and time-harmonic Maxwell wave operators. We prove that such a weak coercivity holds for these wave operators on a class of more general domains called generalized star-shape domains. As a by-product, solution estimates for the corresponding Helmholtz-type problems are obtained.

The second part of the dissertation develops an absolutely stable (i.e. stable in all mesh regimes) interior …

Finite Difference And Discontinuous Galerkin Finite Element Methods For Fully Nonlinear Second Order Partial Differential Equations, Thomas Lee Lewis

Doctoral Dissertations

The dissertation focuses on numerically approximating viscosity solutions to second order fully nonlinear partial differential equations (PDEs). The primary goals of the dissertation are to develop, analyze, and implement a finite difference (FD) framework, a local discontinuous Galerkin (LDG) framework, and an interior penalty discontinuous Galerkin (IPDG) framework for directly approximating viscosity solutions of fully nonlinear second order elliptic PDE problems with Dirichlet boundary conditions. The developed frameworks are also extended to fully nonlinear second order parabolic PDEs. All of the proposed direct methods are tested using Monge-Ampere problems and Hamilton-Jacobi-Bellman (HJB) problems. Due to the significance of HJB problems …

Validation Of Weak Form Thermal Analysis Algorithms Supporting Thermal Signature Generation, Elton Lewis Freeman

Masters Theses

Extremization of a weak form for the continuum energy conservation principle differential equation naturally implements fluid convection and radiation as flux Robin boundary conditions associated with unsteady heat transfer. Combining a spatial semi-discretization via finite element trial space basis functions with time-accurate integration generates a totally node-based algebraic statement for computing. Closure for gray body radiation is a newly derived node-based radiosity formulation generating piecewise discontinuous solutions, while that for natural-forced-mixed convection heat transfer is extracted from the literature. Algorithm performance, mathematically predicted by asymptotic convergence theory, is subsequently validated with data obtained in 24 hour diurnal field experiments for …

Hard And Soft Error Resilience For One-Sided Dense Linear Algebra Algorithms, Peng Du

Doctoral Dissertations

Dense matrix factorizations, such as LU, Cholesky and QR, are widely used by scientific applications that require solving systems of linear equations, eigenvalues and linear least squares problems. Such computations are normally carried out on supercomputers, whose ever-growing scale induces a fast decline of the Mean Time To Failure (MTTF). This dissertation develops fault tolerance algorithms for one-sided dense matrix factorizations, which handles Both hard and soft errors.

For hard errors, we propose methods based on diskless checkpointing and Algorithm Based Fault Tolerance (ABFT) to provide full matrix protection, including the left and right factor that are normally seen in …

Numerical Analysis Of First And Second Order Unconditional Energy Stable Schemes For Nonlocal Cahn-Hilliard And Allen-Cahn Equations, Zhen Guan

Doctoral Dissertations

This PhD dissertation concentrates on the numerical analysis of a family of fully discrete, energy stable schemes for nonlocal Cahn-Hilliard and Allen-Cahn type equations, which are integro-partial differential equations (IPDEs). These two IPDEs -- along with the evolution equation from dynamical density functional theory (DDFT), which is a generalization of the nonlocal Cahn-Hilliard equation -- are used to model a variety of physical and biological processes such as crystallization, phase transformations, and tumor growth. This dissertation advances the computational state-of-the-art related to this field in the following main contributions: (I) We propose and analyze a family of two-dimensional unconditionally energy …

Analysis Of Solvability And Applications Of Stochastic Optimal Control Problems Through Systems Of Forward-Backward Stochastic Differential Equations., Kirill Yevgenyevich Yakovlev

Doctoral Dissertations

A stochastic metapopulation model is investigated. The model is motivated by a deterministic model previously presented to model the black bear population of the Great Smoky Mountains in east Tennessee. The new model involves randomness and the associated methods and results differ greatly from the deterministic analogue. A stochastic differential equation is studied and the associated results are stated and proved. Connections between a parabolic partial differential equation and a system of forward-backward stochastic differential equations is analyzed.

A "four-step" numerical scheme and a Markovian type iterative numerical scheme are implemented. Algorithms and programs in the programming languages C and …