Numerical Analysis and Computation Commons^™

Numerical Analysis Of A Model For The Growth Of Microorganisms, Alexander Craig Montgomery, Braden J. Carlson

Rose-Hulman Undergraduate Mathematics Journal

A system of first-order differential equations that arises in a model for the growth of microorganisms in a chemostat with Monod kinetics is studied. A new, semi-implicit numerical scheme is proposed to approximate solutions to the system. It is shown that the scheme is uniquely solvable and unconditionally stable, and further properties of the scheme are analyzed. The convergence rate of the numerical solution to the true solution of the system is given, and it is shown convergence of the numerical solutions to the true solutions is uniform over any interval [0, T ] for T > 0.

A Conservative Numerical Scheme For The Multilayer Shallow Water Equations, Evan Butterworth

All Theses

An energy-conserving numerical scheme is developed for the multilayer shallow water equations (SWE’s). The scheme is derived through the Hamiltonian formulation of the inviscid shallow water flows related to the vorticity-divergence variables. Through the employment of the skew-symmetric Poisson bracket, the continuous system for the multilayer SWE’s is shown to preserve an infinite number of quantities, most notably the energy and enstrophy. An energy-preserving numerical scheme is then developed through the careful discretization of the Hamiltonian and the Poisson bracket, ensuring the skew-symmetry of the latter. This serves as the groundwork for developing additional schemes that preserve other conservation properties …

On The Consistency Of Alternative Finite Difference Schemes For The Heat Equation, Tran April

Rose-Hulman Undergraduate Mathematics Journal

While the well-researched Finite Difference Method (FDM) discretizes every independent variable into algebraic equations, Method of Lines discretizes all but one dimension, leaving an Ordinary Differential Equation (ODE) in the remaining dimension. That way, ODE's numerical methods can be applied to solve Partial Differential Equations (PDEs). In this project, Linear Multistep Methods and Method of Lines are used to numerically solve the heat equation. Specifically, the explicit Adams-Bashforth method and the implicit Backward Differentiation Formulas are implemented as Alternative Finite Difference Schemes. We also examine the consistency of these schemes.

Multi-Valued Solutions For The Equation Of Motion, Darcy-Jordan Model, As A Cauchy Problem: A Shocking Event, Chandler Shimp

Master's Theses

Shocks are physical phenomenon that occur quite often around us. In this thesis we examine the occurrence of shocks in finite amplitude acoustic waves from a numerical perspective. These waves, or jump discontinuities, yield ill-behaved solutions when solved numerically. This study takes on the challenge of finding both single- and multi-valued solutions.

The previously unsolved problem in this study is the representation of the Equation of Motion (EoM) in the form of the Darcy-Jordan model (DJM) and expressed as a dimensionless IVP Cauchy problem. Prior attempts to solve have resulted only in implicit solutions or explicit solutions with certain initial …

Numerical Integration Through Concavity Analysis, Daniel J. Pietz

Rose-Hulman Undergraduate Mathematics Journal

We introduce a relationship between the concavity of a C2 func- tion and the area bounded by its graph and secant line. We utilize this relationship to develop a method of numerical integration. We then bound the error of the approximation, and compare to known methods, finding an improvement in error bound over methods of comparable computational complexity.

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 …

Analysis Of A Partial Differential Equation Model Of Surface Electromigration, Selahittin Cinar

Masters Theses & Specialist Projects

A Partial Differential Equation (PDE) based model combining surface electromigration and wetting is developed for the analysis of the morphological instability of mono-crystalline metal films in a high temperature environment typical to operational conditions of microelectronic interconnects. The atomic mobility and surface energy of such films are anisotropic, and the model accounts for these material properties. The goal of modeling is to describe and understand the time-evolution of the shape of film surface. I will present the formulation of a nonlinear parabolic PDE problem for the height function h(x,t) of the film in the horizontal …

A Posteriori Error Estimates For Surface Finite Element Methods, Fernando F. Camacho

Theses and Dissertations--Mathematics

Problems involving the solution of partial differential equations over surfaces appear in many engineering and scientific applications. Some of those applications include crystal growth, fluid mechanics and computer graphics. Many times analytic solutions to such problems are not available. Numerical algorithms, such as Finite Element Methods, are used in practice to find approximate solutions in those cases.

In this work we present L2 and pointwise a posteriori error estimates for Adaptive Surface Finite Elements solving the Laplace-Beltrami equation −△_Γ u = f . The two sources of errors for Surface Finite Elements are a Galerkin error, and a …

Analysis Of Roms Estimated Posterior Error Utilizing 4dvar Data Assimilation, Joseph Patrick Horton

Mathematics

The appropriateness of the approximate error calculated by the Regional Ocean Modeling System (ROMS) is analyzed using Four-Dimensional Data Assimilation (4DVAR) performed on a numerical model of the San Luis Obispo Bay. An effective method of sampling data to minimize the actual error associated with the assimilated numerical model is explored by using different data sampling methods. An idealized state of the SLO bay region ("Real Run") is created to be used as the real ocean, then a numerical model of this region is created approximating this Real Run; this is known as the "Simulated State". By taking samples from …