Open Access. Powered by Scholars. Published by Universities.®

Numerical Analysis and Computation Commons™

Open Access. Powered by Scholars. Published by Universities.®

Rose-Hulman Institute of Technology

Discipline
Keyword
Publication Year
Publication
Publication Type

Articles 1 - 26 of 26

Full-Text Articles in Numerical Analysis and Computation

Modeling An Infection Outbreak With Quarantine: The Sibkr Model, Mikenna Dew, Amanda Langosch, Theadora Baker-Wallerstein Apr 2024

Modeling An Infection Outbreak With Quarantine: The Sibkr Model, Mikenna Dew, Amanda Langosch, Theadora Baker-Wallerstein

Rose-Hulman Undergraduate Mathematics Journal

Influenza is a respiratory infection that places a substantial burden in the world population each year. In this project, we study and interpret a data set from a flu outbreak in a British boarding school in 1978 with mathematical modeling. First, we propose a generalization of the SIR model based on the quarantine measure in place and establish the long-time behavior of the model. By analyzing the model mathematically, we determine the analytic formulas of the basic reproduction number, the long-time limit of solutions, and the maximum number of infection population. Moreover, we estimate the parameters of the model based …

Go to article

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

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.

Go to article

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

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.

Go to article

Numerical Integration Through Concavity Analysis, Daniel J. Pietz Jan 2021

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.

Go to article

Algorithms To Approximate Solutions Of Poisson's Equation In Three Dimensions, Ray Dambrose Mar 2019

Algorithms To Approximate Solutions Of Poisson's Equation In Three Dimensions, Ray Dambrose

Rose-Hulman Undergraduate Mathematics Journal

The focus of this research was to develop numerical algorithms to approximate solutions of Poisson's equation in three dimensional rectangular prism domains. Numerical analysis of partial differential equations is vital to understanding and modeling these complex problems. Poisson's equation can be approximated with a finite difference approximation. A system of equations can be formed that gives solutions at internal points of the domain. A computer program was developed to solve this system with inputs such as boundary conditions and a nonhomogenous source function. Approximate solutions are compared with exact solutions to prove their accuracy. The program is tested with an …

Go to article

Analytical Solution Of The Symmetric Circulant Tridiagonal Linear System, Sean A. Broughton, Jeffery J. Leader Aug 2014

Analytical Solution Of The Symmetric Circulant Tridiagonal Linear System, Sean A. Broughton, Jeffery J. Leader

Mathematical Sciences Technical Reports (MSTR)

A circulant tridiagonal system is a special type of Toeplitz system that appears in a variety of problems in scientific computation. In this paper we give a formula for the inverse of a symmetric circulant tridiagonal matrix as a product of a circulant matrix and its transpose, and discuss the utility of this approach for solving the associated system.

Go to article

A Spectral Approach To Protein Structure Alignment, Yosi Shibberu, Allen Holder Aug 2010

A Spectral Approach To Protein Structure Alignment, Yosi Shibberu, Allen Holder

Mathematical Sciences Technical Reports (MSTR)

We present two algorithms that use spectral methods to align protein folds. One of the algorithms is suitable for database searches, the other for difficult alignments. We present computational results for 780 pairwise alignments used to classify 40 proteins as well as results for a separate set of 36 protein alignments used for comparison to four other alignment algorithms. We also provide a mathematically rigorous development of the intrinsic geometry underlying our spectral approach.

Go to article

Fast Protein Structure Alignment, Yosi Shibberu, Allen Holder, Kyla Lutz Feb 2010

Fast Protein Structure Alignment, Yosi Shibberu, Allen Holder, Kyla Lutz

Mathematical Sciences Technical Reports (MSTR)

We address the problem of aligning the 3D structures of two proteins. Our pairwise comparisons are based on a new optimization model that is succinctly expressed in terms of linear transformations and highlights the problem’s intrinsic geometry. The optimization problem is approximately solved with a new polynomial time algorithm. The worst case analysis of the algorithm shows that the solution is bounded by a constant depending only on the data of the problem.

Go to article

Electrothermal Imaging In One And Two Dimensions, Michael Janus, David Kibbling Dec 2008

Electrothermal Imaging In One And Two Dimensions, Michael Janus, David Kibbling

Mathematical Sciences Technical Reports (MSTR)

Developing methods for the nondestructive testing of materials is an important area of research for industry. Situations often arise in which the integrity of an object is questioned, but testing it is very difficult. For example, a support bar may be embedded in a larger structure so that testing the bar’s integrity directly would require the impractical task of breaking down the larger structure. Instead, the ends of the bar might be accessible without dismantling the enclosing structure. The goal of nondestructive testing is to use methods that require taking measurements at the ends of the bar alone to give …

Go to article

Determining The Shape Of A Resistor Grid, Esther Chiew, Vincent Selhorst-Jones Dec 2008

Determining The Shape Of A Resistor Grid, Esther Chiew, Vincent Selhorst-Jones

Mathematical Sciences Technical Reports (MSTR)

Impedance imaging has received a lot of attention in the past two decades, as a means for non-destructively imaging the interior of a conductive object. One injects a known electrical current pattern into an object at the exterior boundary, then measures the induced potential (voltage) on some portion of the boundary. The goal is to recover information about the interior conductivity of the object, which (we hope) influences the voltages we measure. Of course one can also use multiple input currents and measured voltages. A variation on this problem is that of "boundary identification," in which some portion of the …

Go to article

Utilizing Thermal Testing For Recovering, James Preciado, Thomas Werne May 2007

Utilizing Thermal Testing For Recovering, James Preciado, Thomas Werne

Mathematical Sciences Technical Reports (MSTR)

Given a two-dimensional region that contains one or more circular voids, we develop mathematical methods to locate the center and radius of the voids based on thermal boundary data. These methods can be readily applied in the field of non-destructive evaluation.

Go to article

Thermal Imaging To Recover A Defect In Three Dimensional Objects, Breanne Baker Jul 2006

Thermal Imaging To Recover A Defect In Three Dimensional Objects, Breanne Baker

Mathematical Sciences Technical Reports (MSTR)

This paper focuses on the inverse problem of identifying an internal void in a bounded two- or three-dimensional region. Information, in form of a heat flux and temperature, is assumed to be obtainable only on the external boundary of the region. The reciprocity gap approach with a suitable test functions is used in both the two- and three-dimensional cases.

Go to article

Non-Destructive Recovery Of Voids Within A Three Dimensional Domain Using Thermal Imaging, Victor B. Oyeyemi Jul 2006

Non-Destructive Recovery Of Voids Within A Three Dimensional Domain Using Thermal Imaging, Victor B. Oyeyemi

Mathematical Sciences Technical Reports (MSTR)

We develop an algorithm capable of detecting the presence of spherical voids in a thermally conducting object. In addition, the process recovers both the radii and locations of each void. Our method involves the application of a known steady state heat flux to the object's boundary and measurement of the induced steady state temperature on the boundary.

Go to article

Time-Dependent Thermal Imaging Of Circular Inclusions, Donald L. Brouwn, Mark Hubenthal Jul 2005

Time-Dependent Thermal Imaging Of Circular Inclusions, Donald L. Brouwn, Mark Hubenthal

Mathematical Sciences Technical Reports (MSTR)

This paper considers the inverse problem of locating one or more circular inclusions in a two-dimensional domain using thermal boundary data, specifically, the input heat flux and measured boundary temperature. The forward problem is governed by the heat equation. We show how the position and size of such defects can be recovered using the boundary data and various approximations of the solution to the forward problem. We also consider the stability of the algorithm involved to recover the defects.

Go to article

Reconstruction Of Partially Conductive Cracks Using Boundary Data, David Mccune, Janine Haugh Sep 2004

Reconstruction Of Partially Conductive Cracks Using Boundary Data, David Mccune, Janine Haugh

Mathematical Sciences Technical Reports (MSTR)

This paper develops an algorithm for finding one or more non-insulated, pair-wise disjoint, linear cracks in a two dimensional region using boundary measurements.

Go to article

Non-Destructive Testing Of Thermal Resistances For A Single Inclusion In A 2-Dimensional Domain, Nicholas Christian, Mathew A. Johnson Sep 2004

Non-Destructive Testing Of Thermal Resistances For A Single Inclusion In A 2-Dimensional Domain, Nicholas Christian, Mathew A. Johnson

Mathematical Sciences Technical Reports (MSTR)

In this paper we examine the inverse problem of determining the amount of corrosion/disbonding which has occurred on the boundary of a single circular (or nearly circular) inclusion D in a two dimensional domain W using Cauchy data for the steady-state heat equation. We develop an algorithm for reconsructing a function which qunatifies the level of corrosion/disbonding at each point in ¶W. We also address the issue of well-posedness and develop a simple regularization scheme. Then we provide several numerical examples. We shall show a simple procedure for recovering the center of D assuming that the boundary of W and …

Go to article

Reconstruction Of An Unknown Boundary Portion From Cauchy Data In N-Dimensions, Kurt M. Bryan, Lester Caudill Jul 2004

Reconstruction Of An Unknown Boundary Portion From Cauchy Data In N-Dimensions, Kurt M. Bryan, Lester Caudill

Mathematical Sciences Technical Reports (MSTR)

We consider the inverse problem of determining the shape of some inacces­ sible portion of the boundary of a region in n dimensions from Cauchy data for the heat equation on an accessible portion of the boundary. The inverse problem is quite ill-posed, and nonlinear. We develop a Newton-like algorithm for solving the problem, with a simple and efficient means for computing the required derivatives, develop methods for regularizing the process, and provide computational examples

Go to article

Determining The Length Of A One-Dimensional Bar, Natalya Yarlikina, Holly Walrath Jul 2004

Determining The Length Of A One-Dimensional Bar, Natalya Yarlikina, Holly Walrath

Mathematical Sciences Technical Reports (MSTR)

In this paper we examine the inverse problem of determining the length of a one-dimensional bar from thermal measurements (temperature and heat flux) at one end of the bar (the "accessible" end); the other inaccessible end of the bar is assumed to be moving. We develop two different approaches to estimating the length of the bar, and show how one approach can also be adapted to find unknown boundary conditions at the inaccessible end of the bar.

Go to article

A Review Of Selected Works On Crack Indentification, Kurt M. Bryan Oct 2002

A Review Of Selected Works On Crack Indentification, Kurt M. Bryan

Mathematical Sciences Technical Reports (MSTR)

We give a short survey of some of the results obtained within the last 10 years or so concerning crack identification using impedance imaging techniques. We touch upon uniqueness results, continuous dependence results, and computational algorithms.

Go to article

Characterizing A Defect In A One-Dimensional Bar, Cynthia Gangi, Sameer Shah Sep 2002

Characterizing A Defect In A One-Dimensional Bar, Cynthia Gangi, Sameer Shah

Mathematical Sciences Technical Reports (MSTR)

We examine the inverse problem of locating and describing an internal point defect in a one­ dimensional rod W by controlling the heat inputs and measuring the subsequent temperatures at the boundary of W. We use a variation of the forward heat equation to model heat flow through W, then propose algorithms for locating an internal defect and quantifying the effect the defect has on the heat flow. We implement these algorithms, analyze the stability of the procedures, and provide several computational examples.

Go to article

Reconstruction Of Cracks With Unknown Transmission Condition From Boundary Data, F Ronald Ogborne Iii, Melissa E. Vellela Sep 2002

Reconstruction Of Cracks With Unknown Transmission Condition From Boundary Data, F Ronald Ogborne Iii, Melissa E. Vellela

Mathematical Sciences Technical Reports (MSTR)

We examine the problem of Identifying both the location and constitutive law governing electrical current flow across a one-dimensional linear crack in a two dimensional region when the crack only partially blocks the flow of current. We develop a a constructive numerical procedure for solving the inverse problem and provide computational examples.

Go to article

Fast Reconstruction Of Cracks Using Boundary Measurements, Nicholas A. Trainor, Rachel M. Krieger Jul 2002

Fast Reconstruction Of Cracks Using Boundary Measurements, Nicholas A. Trainor, Rachel M. Krieger

Mathematical Sciences Technical Reports (MSTR)

This paper develops a fast algorithm for locating one or more perfectly insulating, pair-wise disjoint, linear cracks in a homogeneous two-dimensional electrical conductor, using boundary measurements.

Go to article

Stability And Reconstruction For An Inverse Problem For The Heat Equations, Kurt M. Bryan, Lester Caudill May 1998

Stability And Reconstruction For An Inverse Problem For The Heat Equations, Kurt M. Bryan, Lester Caudill

Mathematical Sciences Technical Reports (MSTR)

We examine the inverse problem of determining the shape of some unknown portion of the boundary of a region W from measurements of the Cauchy data for solutions to the heat equation on W. By suitably linearizing the inverse problem we obtain uniqueness and continuous dependence results. We propose an algorithm for recovering estimates of the unknown portion of the surface and use the insight gained from a detailed analysis of the inverse problem to regularize the inversion. Several computational examples are presented.

Go to article

Effective Behavior Of Clusters Of Microscopic Cracks Inside A Homogeneous Conductor, Kurt M. Bryan, Michael Vogelius Jan 1996

Effective Behavior Of Clusters Of Microscopic Cracks Inside A Homogeneous Conductor, Kurt M. Bryan, Michael Vogelius

Mathematical Sciences Technical Reports (MSTR)

We study the effective behaviour of a periodic array of microscopic cracks inside a homoge­neous conductor. Special emphasis is placed on a rigorous study of the case in which the corresponding effective conductivity becomes nearly singular, due to the fact that adjacent cracks nearly touch. It is heuristically shown how thin clusters of such extremely close cracks may macroscopically appear as a single crack. The results have implications for our earlier work on impedance imaging.

Go to article

A Variable Time-Step Midpoint Scheme For Hamiltonian Systems, Yosi Shibberu May 1995

A Variable Time-Step Midpoint Scheme For Hamiltonian Systems, Yosi Shibberu

Mathematical Sciences Technical Reports (MSTR)

A smooth time-step selection formula for the midpoint method is derived which minimize deviations in the Hamiltonian function along piecewise-linear phase space trajectories of autonomous Hamiltonian systems. The time-step formula is implemented in a second order pre­dictor/corrector scheme and applied to Kepler's problem. The formula significantly improves energy conservation as well as the accuracy of the configuration space trajectory. Peak errors in position and momentum coordinates are not significantly reduced, but the time behavior of the errors is markedly more regular.

Go to article

Time-Discretization Of Hamiltonian Dynamical Systems, Yosi Shibberu Jun 1993

Time-Discretization Of Hamiltonian Dynamical Systems, Yosi Shibberu

Mathematical Sciences Technical Reports (MSTR)

Difference equations for Hamiltonian systems are derived from a discrete variational principle. The difference equations completely determine piecewise-linear, continuous trajectories which exactly conserve the Hamiltonian function at the midpoints of each linear segment. A generating function exists for transformations between the vertices of the trajectories. Existence and uniqueness results are present as well as simulation results for a simple pendulum and an inverse square law system.

Go to article