Physical Sciences and Mathematics Commons^™

An Optimal Decay Estimation Of The Solution To The Airy Equation, Ashley Scherf

Theses and Dissertations

In this thesis, we investigate the initial value problem to the Airy equation \begin{align} \partial_t u + \partial_{x}^3 u &= 0\\ u(0,x) &= f(x). \end{align}

Bloch Spectra For High Contrast Elastic Media, Jayasinghage Ruchira Nirmali Perera

LSU Doctoral Dissertations

The primary goal of this dissertation is to develop analytic representation formulas and power series to describe the band structure inside periodic elastic crystals made from high contrast inclusions. We use source free modes associated with structural spectra to represent the solution operator of the Lame' system inside phononic crystals. Then we obtain convergent power series for the Bloch wave spectrum using the representation formulas. An explicit bound on the convergence radius is given through the structural spectra of the inclusion array and the Dirichlet spectra of the inclusions. Sufficient conditions for the separation of spectral branches of the dispersion …

Sparse Spectral-Tau Method For The Two-Dimensional Helmholtz Problem Posed On A Rectangular Domain, Gabriella M. Dalton

Mathematics & Statistics ETDs

Within recent decades, spectral methods have become an important technique in numerical computing for solving partial differential equations. This is due to their superior accuracy when compared to finite difference and finite element methods. For such spectral approximations, the convergence rate is solely dependent on the smoothness of the solution yielding the potential to achieve spectral accuracy. We present an iterative approach for solving the two-dimensional Helmholtz problem posed on a rectangular domain subject to Dirichlet boundary conditions that is well-conditioned, low in memory, and of sub-quadratic complexity. The proposed approach spectrally approximates the partial differential equation by means of …

Diagonalization Of 1-D Schrodinger Operators With Piecewise Constant Potentials, Sarah Wright

Master's Theses

In today's world our lives are very layered. My research is meant to adapt current inefficient numerical methods to more accurately model the complex situations we encounter. This project focuses on a specific equation that is used to model sound speed in the ocean. As depth increases, the sound speed changes. This means the variable related to the sound speed is not constant. We will modify this variable so that it is piecewise constant. The specific operator in this equation also makes current time-stepping methods not practical. The method used here will apply an eigenfunction expansion technique used in previous …

Methods In Modeling Wildlife Disease From Model Selection To Parameterization With Multi-Scale Data, Ian Mcgahan

All Graduate Theses and Dissertations, Spring 1920 to Summer 2023

The effects of emerging wildlife diseases are global and profound, resulting in loss of human life, economic and agricultural impacts, declines in wildlife populations, and ecological disturbance. The spread of wildlife diseases can be viewed as the result of two simultaneous processes: spatial spread of wildlife populations and disease spread through a population. For many diseases these processes happen at different timescales, which is reflected in available data. These data come in two flavors: high-frequency, high-resolution telemetry data (e.g. GPS collar) and low-frequency, low-resolution presence-absence disease data. The multi-scale nature of these data makes analysis of such systems challenging. Mathematical …

Hadamard Well-Posedness For Two Nonlinear Structure Acoustic Models, Andrew Becklin

Department of Mathematics: Dissertations, Theses, and Student Research

This dissertation focuses on the Hadamard well-posedness of two nonlinear structure acoustic models, each consisting of a semilinear wave equation defined on a smooth bounded domain $\Omega\subset\mathbb{R}^3$ strongly coupled with a Berger plate equation acting only on a flat portion of the boundary of $\Omega$. In each case, the PDE is of the following form: \begin{align*} \begin{cases} u_{tt}-\Delta u +g_1(u_t)=f(u) &\text{ in } \Omega \times (0,T),\\[1mm] w_{tt}+\Delta^2w+g_2(w_t)+u_t|_{\Gamma}=h(w)&\text{ in }\Gamma\times(0,T),\\[1mm] u=0&\text{ on }\Gamma_0\times(0,T),\\[1mm] \partial_\nu u=w_t&\text{ on }\Gamma\times(0,T),\\[1mm] w=\partial_{\nu_\Gamma}w=0&\text{ on }\partial\Gamma\times(0,T),\\[1mm] (u(0),u_t(0))=(u_0,u_1),\hspace{5mm}(w(0),w_t(0))=(w_0,w_1), \end{cases} \end{align*} where the initial data reside in the finite energy space, i.e., $$(u_0, u_1)\in H^1_{\Gamma_0}(\Omega) \times L^2(\Omega) \, \text{ …

Mathematical Modeling Of Microemulsification Processes, Numerical Simulations And Applications To Drug Delivery, Ogochukwu Nneka Ifeacho

Open Access Theses & Dissertations

Microemulsion systems are a great pharmaceutical tool for the delivery of formulations containing multiple hydrophilic and hydrophobic ingredients of varying physicochemical properties. These systems are gaining popularity because of its long shelf life, improved drug solubilisation capacity, easy preparation and improvement of bioavailability. Despite the advantages associated with the use of microemulsion systems in pharmaceutical industries, the major challenge impeding their use has been and continues to be the lack of understanding of these systems.

Microemulsions can be mathematically modeled by an initial boundary value problem involving a sixth order nonlinear time dependent equation. In this Thesis, we present a …

The Direct Scattering Map For The Intermediate Long Wave Equation, Joel Klipfel

Theses and Dissertations--Mathematics

In the early 1980's, Kodama, Ablowitz and Satsuma, together with Santini, Ablowitz and Fokas, developed the formal inverse scattering theory of the Intermediate Long Wave (ILW) equation and explored its connections with the Benjamin-Ono (BO) and KdV equations. The ILW equation\begin{align*} u_t + \frac{1}{\delta} u_x + 2 u u_x + Tu_{xx} = 0, \end{align*} models the behavior of long internal gravitational waves in stratified fluids of depth $0< \delta < \infty$, where $T$ is a singular operator which depends on the depth $\delta$. In the limit $\delta \to 0$, the ILW reduces to the Korteweg de Vries (KdV) equation, and in the limit $\delta \to \infty$, the ILW (at least formally) reduces to the Benjamin-Ono (BO) equation.

While the KdV equation is very well understood, a rigorous analysis of inverse scattering for the ILW equation remains to be accomplished. There is currently no rigorous proof that the Inverse Scattering …

Enhancement Of Krylov Subspace Spectral Methods Through The Use Of The Residual, Haley Dozier

Dissertations

Depending on the type of equation, finding the solution of a time-dependent partial differential equation can be quite challenging. Although modern time-stepping methods for solving these equations have become more accurate for a small number of grid points, in a lot of cases the scalability of those methods leaves much to be desired. That is, unless the timestep is chosen to be sufficiently small, the computed solutions might exhibit unreasonable behavior with large input sizes. Therefore, to improve accuracy as the number of grid points increases, the time-steps must be chosen to be even smaller to reach a reasonable solution. …

Automatic Construction Of Scalable Time-Stepping Methods For Stiff Pdes, Vivian Ashley Montiforte

Master's Theses

Krylov Subspace Spectral (KSS) Methods have been demonstrated to be highly scalable time-stepping methods for stiff nonlinear PDEs. However, ensuring this scalability requires analytic computation of frequency-dependent quadrature nodes from the coefficients of the spatial differential operator. This thesis describes how this process can be automated for various classes of differential operators to facilitate public-domain software implementation.

Survey Of Results On The Schrodinger Operator With Inverse Square Potential, Richardson Saint Bonheur

Electronic Theses and Dissertations

In this paper we present a survey of results on the Schrodinger operator with Inverse ¨ Square potential, L_a= −∆ + a/|x|^2 , a ≥ −( d−2/2 )^2. We briefly discuss the long-time behavior of solutions to the inter-critical focusing NLS with an inverse square potential(proof not provided). Later we present spectral multiplier theorems for the operator. For the case when a ≥ 0, we present the multiplier theorem from Hebisch [12]. The case when 0 > a ≥ −( d−2/2 )^2 was explored in [1], and their proof will be presented for completeness. No improvements on the sharpness …

Nonlocal Electrostatics In Spherical Geometries, Andrew Bolanowski

Theses and Dissertations

Nonlocal continuum electrostatic models have been used numerically in protein simulations, but analytic solutions have been absent. In this paper, two modified nonlocal continuum electrostatic models, the Lorentzian Model and a Linear Poisson-Boltzmann Model, are presented for a monatomic ion treated as a dielectric continuum ball. These models are then solved analytically using a system of differential equations for the charge distributed within the ion ball. This is done in more detail for a point charge and a charge distributed within a smaller ball. As the solutions are a series, their convergence is verified and criteria for improved convergence is …

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 …

Steady State Solutions For A System Of Partial Differential Equations Arising From Crime Modeling, Bo Li

HMC Senior Theses

I consider a model for the control of criminality in cities. The model was developed during my REU at UCLA. The model is a system of partial differential equations that simulates the behavior of criminals and where they may accumulate, hot spots. I have proved a prior bounds for the partial differential equations in both one-dimensional and higher dimensional case, which proves the attractiveness and density of criminals in the given area will not be unlimitedly high. In addition, I have found some local bifurcation points in the model.

Spectral Properties Of Photonic Crystals: Bloch Waves And Band Gaps, Robert Paul Viator Jr

LSU Doctoral Dissertations

The author of this dissertation studies the spectral properties of high-contrast photonic crystals, i.e. periodic electromagnetic waveguides made of two materials (a connected phase and included phase) whose electromagnetic material properties are in large contrast. A spectral analysis of 2nd-order divergence-form partial differential operators (with a coupling constant k) is provided. A result of this analysis is a uniformly convergent power series representation of Bloch-wave eigenvalues in terms of the coupling constant k in the high-contrast limit k -> infinity. An explicit radius of convergence for this power series is obtained, and can be written explicitly in terms of the …

Application Of A Numerical Method And Optimal Control Theory To A Partial Differential Equation Model For A Bacterial Infection In A Chronic Wound, Stephen Guffey

Masters Theses & Specialist Projects

In this work, we study the application both of optimal control techniques and a numerical method to a system of partial differential equations arising from a problem in wound healing. Optimal control theory is a generalization of calculus of variations, as well as the method of Lagrange Multipliers. Both of these techniques have seen prevalent use in the modern theories of Physics, Economics, as well as in the study of Partial Differential Equations. The numerical method we consider is the method of lines, a prominent method for solving partial differential equations. This method uses finite difference schemes to discretize the …

Mapping Algorithms And Software Environment For Data Parallel, Nikos Chrisochoides, Elias Houstis, John Rice

John R Rice

We consider computations associated with data parallel iterative solvers used for the numerical solution of Partial Differential Equations (PDEs). The mapping of such computations into load balanced tasks requiring minimum synchronization and communication is a difficult combinatorial optimization problem. Its optimal solution is essential for the efficient parallel processing of PDE computations. Determining data mappings that optimize a number of criteria, like workload balance, synchronization and local communication, often involves the solution of an NP-Complete problem. Although data mapping algorithms have been known for a few years there is lack of qualitative and quantitative comparisons based on the actual performance …

Mathematical Modeling Of Competition For Light And Nutrients Between Phytoplankton Species In A Poorly Mixed Water Column, Thomas George Stojsavljevic

Theses and Dissertations

Phytoplankton live in a complex environment with two essential resources forming various gradients. Light supplied from above is never homogeneously distributed in a body of water due to refraction and absorption from biomass present in the ecosystem and from other sources. Nutrients in turn are typically supplied from below. In poorly mixed water columns phytoplankton can be heterogeneously distributed forming various layering patterns. The relationship between the location and the thickness of the layers is an open problem of interest. Here we present three models which study how competition for light and resources can form common layering patterns seen in …

A Maximum Principle In The Engel Group, James Klinedinst

USF Tampa Graduate Theses and Dissertations

In this thesis, we will examine the properties of subelliptic jets in the Engel group of step 3. Step-2 groups, such as the Heisenberg group, do not provide insight into the general abstract calculations. This thesis then, is the first explicit non-trivial computation of the abstract results.

Using Partial Differential Equations To Model And Analyze The Treatment Of A Chronic Wound With Oxygen Therapy Techniques, Brandon C. Russell

Mahurin Honors College Capstone Experience/Thesis Projects

Chronic wounds plague approximately 1.3-3 million Americans. The treatment of these wounds requires knowledge of the complex healing process of typical wounds. With a system of partial differential equations, this project attempts to model the intricate biological process and to describe oxygen levels, neutrophil and bacteria concentrations, and other biological parameters with respect to time and space. Analytical solutions for the model will be derived for various frames of time in the wound-healing process. The system of equations will be numerically solved using Matlab. Numerical simulations are performed to determine optimal treatment strategies for a chronic wound.

Regularity For Solutions To Parabolic Systems And Nonlocal Minimization Problems, Joe Geisbauer

Department of Mathematics: Dissertations, Theses, and Student Research

The goal of this dissertation is to contribute to both the nonlocal and local settings of regularity within the calculus of variations. We provide analogues of higher differentiability results in the context of Besov spaces for minimizers of nonlocal functionals. We also establish the Holder continuity of solutions to a system of parabolic partial differential equations.

Advisor: Mikil Foss

A Block Operator Splitting Method For Heterogeneous Multiscale Poroelasticity, Paul M. Delgado

Open Access Theses & Dissertations

Traditional models of poroelastic deformation in porous media assume relatively homogeneous material properties such that macroscopic constitutive relations lead to accurate results. Many realistic applications involve heterogeneous material properties whose oscillatory nature require multiscale methods to balance accuracy and efficiency in computation.

The current study develops a multiscale method for poroelastic deformation based on a fixed point iteration based operator splitting method and a heterogeneous multiscale method using finite volume and direct stiffness methods. To characterize the convergence

of the operator splitting method, we use a numerical root finding algorithm to determine a threshold surface in a non-dimensional parameter space …

Solving The Partial Differential Equation Of Vibrations With Interval Parameters Using The Interval Finite Difference Method, Brenda G. Medina

Open Access Theses & Dissertations

Accuracy and efficiency are among the main factors that drive today's innovative disciplines. As technology rapidly advances, efficiency takes on new meanings but what about accuracy? How accurate is accurate? Human error, uncertainties in measurement, and rounding errors are just some causes of inaccuracy. Interval Computations is an area that allows for such issues to be taken into account; for each measurement attained (for example), an interval can be built by considering the error associated with the measurement, and such an interval can be utilized in the mathematical computations of interest.

We consider the partial differential equation (PDE) of vibrations …

Symmetry Methods And Self-Similar Solutions To Curve Shortening, Peter Geertz-Larson

Summer Research

Curve shortening is a geometric process that continually evolves a curve based on its curvature.Self-similar solutions to the curve shortening equation maintain their form throughoutthe process, though they can be scaled, translated, or rotated. These self-similar solutionscorrespond to the invariant solutions of the symmetry method for solving differential equations.

Generalizations Of A Laplacian-Type Equation In The Heisenberg Group And A Class Of Grushin-Type Spaces, Kristen Snyder Childers

USF Tampa Graduate Theses and Dissertations

In [2], Beals, Gaveau and Greiner find the fundamental solution to a 2-Laplace-type equation in a class of sub-Riemannian spaces. This fundamental solution is based on the well-known fundamental solution to the p-Laplace equation in Grushin-type spaces [4] and the Heisenberg group [6]. In this thesis, we look to generalize the work in [2] for a p-Laplace-type equation. After discovering that the "natural" generalization fails, we find two generalizations whose solutions are based on the fundamental solution to the p-Laplace equation.

Anti-Cloaking: The Mathematics Of Disguise, Theresa C. Anderson, Brooke E. Phillips

Mathematical Sciences Technical Reports (MSTR)

Recent developments in cloaking, the ability to selectively bend electromagnetic waves so as to render an object invisible, have been abundant. Based on cloaking principles, we will describe several ways to mathematically disguise objects in the context of electrical impedance imaging. Through the use of a change-of-variables scheme we show how one can make an object appear enlarged, translated, or rotated by surrounding it with a suitable "metamaterial," a man-made material that selectively redirects current. Analysis of eigenvectors and eigenvalues, which describe how current flows, follow. We prove that in order to disguise an object, a metamaterial must encompass both …

Adomian Decomposition Method For Solving The Equation Governing The Unsteady Flow Of A Polytropic Gas, M. A. Mohamed

Applications and Applied Mathematics: An International Journal (AAM)

In this article, we have discussed a new application of Adomian decomposition method on nonlinear physical equations. The models of interest in physics are considered and solved by means of Adomian decomposition method. The behavior of Adomian solutions and the effects of different values of time are investigated. Numerical illustrations that include nonlinear physical models are investigated to show the pertinent features of the technique.

Asymptotic Stability Of A Fluid-Structure Semigroup, George Avalos

Department of Mathematics: Faculty Publications

The strong stability problem for a fluid-structure interactive partial differential equation (PDE) is considered. The PDE comprises a coupling of the linearized Stokes equations to the classical system of elasticity, with the coupling occurring on the boundary interface between the fluid and solid media. It is now known that this PDE may be modeled by a $C_{0}$-semigroup of contractions on an appropriate Hilbert space. However, because of the nature of the unbounded coupling between fluid and structure, the resolvent of the semigroup generator will \emph{not} be a compact operator. In consequence, the classical solution to the stability problem, by means …

A Comparison Of Optimization Heuristics For The Data Mapping Problem, Nikos Chrisochoides, Nashat Mansour, Geoffrey C. Fox

Northeast Parallel Architecture Center

In this paper we compare the performance of six heuristics with suboptimal solutions for the data distribution of two dimensional meshes that are used for the numerical solution of Partial Differential Equations (PDEs) on multicomputers. The data mapping heuristics are evaluated with respect to seven criteria covering load balancing, interprocessor communication, flexibility and ease of use for a class of single-phase iterative PDE solvers. Our evaluation suggests that the simple and fast block distribution heuristic can be as effective as the other five complex and computational expensive algorithms.

Mapping Algorithms And Software Environment For Data Parallel, Nikos Chrisochoides, Elias Houstis, John Rice

Northeast Parallel Architecture Center

We consider computations associated with data parallel iterative solvers used for the numerical solution of Partial Differential Equations (PDEs). The mapping of such computations into load balanced tasks requiring minimum synchronization and communication is a difficult combinatorial optimization problem. Its optimal solution is essential for the efficient parallel processing of PDE computations. Determining data mappings that optimize a number of criteria, like workload balance, synchronization and local communication, often involves the solution of an NP-Complete problem. Although data mapping algorithms have been known for a few years there is lack of qualitative and quantitative comparisons based on the actual performance …