Physical Sciences and Mathematics Commons^™

An Application Of M-Matrices To Preserve Bounded Positive Solutions To The Evolution Equations Of Biofilm Models, Richard S. Landry Jr.

University of New Orleans Theses and Dissertations

In this work, we design a linear, two step implicit finite difference method to approximate the solutions of a biological system that describes the interaction between a microbial colony and a surrounding substrate. Three separate models are analyzed, all of which can be described as systems of partial differential equations (PDE)s with nonlinear diffusion and reaction, where the biological colony grows and decays based on the substrate bioavailability. The systems under investigation are all complex models describing the dynamics of biological films. In view of the difficulties to calculate analytical solutions of the models, we design here a numerical technique …

Feasible Computation In Symbolic And Numeric Integration, Robert H.C. Moir

Electronic Thesis and Dissertation Repository

Two central concerns in scientific computing are the reliability and efficiency of algorithms. We introduce the term feasible computation to describe algorithms that are reliable and efficient given the contextual constraints imposed in practice. The main focus of this dissertation then, is to bring greater clarity to the forms of error introduced in computation and modeling, and in the limited context of symbolic and numeric integration, to contribute to integration algorithms that better account for error while providing results efficiently.

Chapter 2 considers the problem of spurious discontinuities in the symbolic integration problem, proposing a new method to restore continuity …

Radial Basis Function Differential Quadrature Method For The Numerical Solution Of Partial Differential Equations, Daniel Watson

Dissertations

In the numerical solution of partial differential equations (PDEs), there is a need for solving large scale problems. The Radial Basis Function Differential Quadrature (RBFDQ) method and local RBF-DQ method are applied for the solutions of boundary value problems in annular domains governed by the Poisson equation, inhomogeneous biharmonic equation, and the inhomogeneous Cauchy-Navier equations of elasticity. By choosing the collocation points properly, linear systems can be obtained so that the coefficient matrices have block circulant structures. The resulting systems can be efficiently solved using matrix decomposition algorithms (MDAs) and fast Fourier transforms (FFTs). For the local RBFDQ method, the …

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 …

Non-Equispaced Fast Fourier Transforms In Turbulence Simulation, Aditya M. Kulkarni

Masters Theses

Fourier pseudo-spectral method on equispaced grid is one of the approaches in turbulence simulation, to compute derivative of discrete data, using fast Fourier Transform (FFT) and gives low dispersion and dissipation errors. In many turbulent flows the dynamically important scales of motion are concentrated in certain regions which requires a coarser grid for higher accuracy. A coarser grid in other regions minimizes the memory requirement. This requires the use of Non-equispaced Fast Fourier Transform (NFFT) to compute the Fourier transform, by solving a system of linear equations.

To achieve similar accuracy, the NFFT needs to return more Fourier coefficients than …

High-Order Relaxed Multirate Infinitesimal Step Methods For Multiphysics Applications, Jean Sexton

Mathematics Theses and Dissertations

In this work, we consider numerical methods for integrating multirate ordinary differential equations. We are interested in the development of new multirate methods with good stability properties and improved efficiency over existing methods. We discuss the development of multirate methods, particularly focusing on those that are based on Runge-Kutta theory. We introduce the theory of Generalized Additive Runge-Kutta methods proposed by Sandu and Günther. We also introduce the theory of Recursive Flux Splitting Multirate Methods with Sub-cycling described by Schlegel, as well as the Multirate Infinitesimal Step methods this work is based on. We propose a generic structure called Flexible …

Information Theoretic Study Of Gaussian Graphical Models And Their Applications, Ali Moharrer

LSU Doctoral Dissertations

In many problems we are dealing with characterizing a behavior of a complex stochastic system or its response to a set of particular inputs. Such problems span over several topics such as machine learning, complex networks, e.g., social or communication networks; biology, etc. Probabilistic graphical models (PGMs) are powerful tools that offer a compact modeling of complex systems. They are designed to capture the random behavior, i.e., the joint distribution of the system to the best possible accuracy. Our goal is to study certain algebraic and topological properties of a special class of graphical models, known as Gaussian graphs. First, …

On The Ramberg-Osgood Stress-Strain Model And Large Deformations Of Cantilever Beams, Ronald J. Giardina Jr

University of New Orleans Theses and Dissertations

In this thesis the Ramberg-Osgood nonlinear model for describing the behavior of many different materials is investigated. A brief overview of the model as it is currently used in the literature is undertaken and several misunderstandings and possible pitfalls in its application is pointed out, especially as it pertains to more recent approaches to finding solutions involving the model. There is an investigation of the displacement of a cantilever beam under a combined loading consisting of a distributed load across the entire length of the beam and a point load at its end and new solutions to this problem are …

Joint Inversion Of Compact Operators, James Ford

Boise State University Theses and Dissertations

The first mention of joint inversion came in [22], where the authors used the singular value decomposition to determine the degree of ill-conditioning in inverse problems. The authors demonstrated in several examples that combining two models in a joint inversion, and effectively stacking discrete linear models, improved the conditioning of the problem. This thesis extends the notion of using the singular value decomposition to determine the conditioning of discrete joint inversion to using the singular value expansion to determine the well-posedness of joint linear operators. We focus on compact linear operators related to geophysical, electromagnetic subsurface imaging.

The operators are …

Some Problems Arising From Mathematical Model Of Ductal Carcinoma In Situ., Heng Li

Electronic Theses and Dissertations

Ductal carcinoma in situ (DCIS) is the earliest form of breast cancer. Three mathematical models in the one dimensional case arising from DCIS are proposed. The first two models are in the form of parabolic equation with initial and known moving boundaries. Direct and inverse problems are considered in model 1, existence and uniqueness are proved by using tool from heat potential theory and Volterra integral equations. Also, we discuss the direct problem and nonlocal problem of model 2, existence and uniqueness are proved. And approximation solution of these problems are implemented by Ritz-Galerkin method, which is the first attempt …

Numerical Solution Of Partial Differential Equations Using Polynomial Particular Solutions, Thir R. Dangal

Dissertations

Polynomial particular solutions have been obtained for certain types of partial differential operators without convection terms. In this dissertation, a closed-form particular solution for more general partial differential operators with constant coefficients has been derived for polynomial basis functions. The newly derived particular solutions are further coupled with the method of particular solutions (MPS) for numerically solving a large class of elliptic partial differential equations. In contrast to the use of Chebyshev polynomial basis functions, the proposed approach is more flexible in selecting the collocation points inside the domain. Polynomial basis functions are well-known for yielding ill-conditioned systems when their …

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 …

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 …

Eignefunctions For Partial Differential Equations On Two-Dimensional Domains With Piecewise Constant Coefficients, Abdullah Muheel Momit Aurko

Master's Theses

In this thesis, we develop a highly accurate and efficient algorithm for computing the solution of a partial differential equation defined on a two-dimensional domain with discontinuous coefficients. An example of such a problem is for modeling the diffusion of heat energy in two space dimensions, in the case where the spatial domain represents a medium consisting of two different but homogeneous materials, with periodic boundary conditions.

Since diffusivity changes based on the material, it will be represented using a piecewise constant function, and this results in the formation of a complicated mathematical model. Such a model is impossible to …

Evaluation And Enhancement Of Clean Energy Systems: Analytical, Computational And Experimental Study Of Solar And Nuclear Cycles, Nima Fathi

Mechanical Engineering ETDs

Clean (and specifically renewable) energy is steadily improving its global share. However, finite availability of fossil fuels and the growing effects of climate change make it an urgent priority to convince the industry and governments to incentivize investment in the renewable energy field and to make it more attractive by decreasing the capital cost. Until recently, uncertainties in funding limited renewable energy development, especially in the US. That limitation has been one of the barriers to progress. Another limitation of many renewable energy systems is the variability in their output, which makes them unsuitable for baseline power production. Therefore, fossil …

Information Metrics For Predictive Modeling And Machine Learning, Kostantinos Gourgoulias

Doctoral Dissertations

The ever-increasing complexity of the models used in predictive modeling and data science and their use for prediction and inference has made the development of tools for uncertainty quantification and model selection especially important. In this work, we seek to understand the various trade-offs associated with the simulation of stochastic systems. Some trade-offs are computational, e.g., execution time of an algorithm versus accuracy of simulation. Others are analytical: whether or not we are able to find tractable substitutes for quantities of interest, e.g., distributions, ergodic averages, etc. The first two chapters of this thesis deal with the study of the …

Fast Algorithms On Random Matrices And Structured Matrices, Liang Zhao

Dissertations, Theses, and Capstone Projects

Randomization of matrix computations has become a hot research area in the big data era. Sampling with randomly generated matrices has enabled fast algorithms to perform well for some most fundamental problems of numerical algebra with probability close to 1. The dissertation develops a set of algorithms with random and structured matrices for the following applications: 1) We prove that using random sparse and structured sampling enables rank-r approximation of the average input matrix having numerical rank r. 2) We prove that Gaussian elimination with no pivoting (GENP) is numerically safe for the average nonsingular and well-conditioned matrix preprocessed with …

Building And Validating A Model For Investigating The Dynamics Of Isolated Water Molecules, Grant Cates

Senior Theses

Understanding how water molecules behave in isolation is vital to understand many fundamental processes in nature. To that end, scientists have begun studying crystals in which single water molecules become trapped in regularly occurring cavities in the crystal structure. As part of that investigation, numerical models used to investigate the dynamics of isolated water molecules are sought to help bolster our fundamental understanding of how these systems behave. To that end, the efficacy of three computational methods—the Euler Method, the Euler-Aspel Method and the Beeman Method—is compared using a newly defined parameter, called the predictive stability coefficient ρ. This …

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

Chancellor’s Honors Program Projects

No abstract provided.

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 …

Efficient Denoising And Sharpening Of Color Images Through Numerical Solution Of Nonlinear Diffusion Equations, Linh T. Duong

Honors Theses

The purpose of this project is to enhance color images through denoising and sharpening, two important branches of image processing, by mathematically modeling the images. Modifications are made to two existing nonlinear diffusion image processing models to adapt them to color images. This is done by treating the red, green, and blue (RGB) channels of color images independently, contrary to the conventional idea that the channels should not be treated independently. A new numerical method is needed to solve our models for high resolution images since current methods are impractical. To produce an efficient method, the solution is represented as …

Electrodynamical Modeling For Light Transport Simulation, Michael G. Saunders

Undergraduate Honors Theses

Modernity in the computer graphics community is characterized by a burgeoning interest in physically based rendering techniques. That is to say that mathematical reasoning from first principles is widely preferred to ad hoc, approximate reasoning in blind pursuit of photorealism. Thereby, the purpose of our research is to investigate the efficacy of explicit electrodynamical modeling by means of the generalized Jones vector given by Azzam [1] and the generalized Jones matrix given by Ortega-Quijano & Arce-Diego [2] in the context of stochastic light transport simulation for computer graphics. To augment the status quo path tracing framework with such a modeling …

Application Of Symplectic Integration On A Dynamical System, William Frazier

Electronic Theses and Dissertations

Molecular Dynamics (MD) is the numerical simulation of a large system of interacting molecules, and one of the key components of a MD simulation is the numerical estimation of the solutions to a system of nonlinear differential equations. Such systems are very sensitive to discretization and round-off error, and correspondingly, standard techniques such as Runge-Kutta methods can lead to poor results. However, MD systems are conservative, which means that we can use Hamiltonian mechanics and symplectic transformations (also known as canonical transformations) in analyzing and approximating solutions. This is standard in MD applications, leading to numerical techniques known as symplectic …

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 …

Computational Algorithms For Improved Representation Of The Model Error Covariance In Weak-Constraint 4d-Var, Jeremy A. Shaw

Dissertations and Theses

Four-dimensional variational data assimilation (4D-Var) provides an estimate to the state of a dynamical system through the minimization of a cost functional that measures the distance to a prior state (background) estimate and observations over a time window. The analysis fit to each information input component is determined by the specification of the error covariance matrices in the data assimilation system (DAS). Weak-constraint 4D-Var (w4D-Var) provides a theoretical framework to account for modeling errors in the analysis scheme. In addition to the specification of the background error covariance matrix, the w4D-Var formulation requires information on the model error statistics and …

Weihrauch Reducibility And Finite-Dimensional Subspaces, Sean Sovine

Theses, Dissertations and Capstones

In this thesis we study several principles involving subspaces and decompositions of vector spaces, matroids, and graphs from the perspective of Weihrauch reducibility. We study the problem of decomposing a countable vector space or countable matroid into 1-dimensional subspaces. We also study the problem of producing a finite-dimensional or 1-dimensional subspace of a countable vector space, and related problems for producing finite-dimensional subspaces of a countable matroid. This extends work in the reverse mathematics setting by Downey, Hirschfeldt, Kach, Lempp, Mileti, and Montalb´an (2007) and recent work of Hirst and Mummert (2017). Finally, we study the problem of producing a …

Computational Fluid Dynamics In A Terminal Alveolated Bronchiole Duct With Expanding Walls: Proof-Of-Concept In Openfoam, Jeremy Myers

Theses and Dissertations

Mathematical Biology has found recent success applying Computational Fluid Dynamics (CFD) to model airflow in the human lung. Detailed modeling of flow patterns in the alveoli, where the oxygen-carbon dioxide gas exchange occurs, has provided data that is useful in treating illnesses and designing drug-delivery systems. Unfortunately, many CFD software packages have high licensing fees that are out of reach for independent researchers. This thesis uses three open-source software packages, Gmsh, OpenFOAM, and ParaView, to design a mesh, create a simulation, and visualize the results of an idealized terminal alveolar sac model. This model successfully demonstrates that OpenFOAM can be …

A Numerical Study Of Construction Of Honey Bee Comb, Pamela Guerrero, Pamela C. Guerrero

Murray State Theses and Dissertations

We use finite difference methods in the treatment of an existing system of partial differential equations that captures the dynamics of parallel honeycomb construction in a bee hive. We conduct an uncertainty analysis by calculating the partial rank correlation coefficient for the parameters to find which are most important to the outcomes of the model. We then use an eFAST method to determine both the individual and total sensitivity index for the parameters. Afterwards we examine our numerical model under varying initial conditions and parameter values, and compare ratios found from local data with the golden mean by fitting images …

Neural Network Predictions Of A Simulation-Based Statistical And Graph Theoretic Study Of The Board Game Risk, Jacob Munson

Murray State Theses and Dissertations

We translate the RISK board into a graph which undergoes updates as the game advances. The dissection of the game into a network model in discrete time is a novel approach to examining RISK. A review of the existing statistical findings of skirmishes in RISK is provided. The graphical changes are accompanied by an examination of the statistical properties of RISK. The game is modeled as a discrete time dynamic network graph, with the various features of the game modeled as properties of the network at a given time. As the network is computationally intensive to implement, results are produced …

Daily Traffic Flow Pattern Recognition By Spectral Clustering, Matthew Aven

CMC Senior Theses

This paper explores the potential applications of existing spectral clustering algorithms to real life problems through experiments on existing road traffic data. The analysis begins with an overview of previous unsupervised machine learning techniques and constructs an effective spectral clustering algorithm that demonstrates the analytical power of the method. The paper focuses on the spectral embedding method’s ability to project non-linearly separable, high dimensional data into a more manageable space that allows for accurate clustering. The key step in this method involves solving a normalized eigenvector problem in order to construct an optimal representation of the original data.

While this …