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

Articles 1 - 16 of 16

Full-Text Articles in Numerical Analysis and Computation

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.


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 …


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 …


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.


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.


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.


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 …


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


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.


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.


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.


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.


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.


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.


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.


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.