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

Partial Differential Equations

Institution
Keyword
Publication Year
Publication
Publication Type
File Type

Articles 241 - 265 of 265

Full-Text Articles in Numerical Analysis and Computation

Mathematical Methods In Composing Melodies, Thomas Brown Apr 2007

Mathematical Methods In Composing Melodies, Thomas Brown

Undergraduate Theses and Capstone Projects

This thesis, “Mathematical Methods in Composing Melodies,” explores the different ways in which mathematics can be used to create music. Some research has been done in this field already. Richard F. Voss and John Clarke used fractals and different frequencies of noise to create music. The Greek composer Iannis Xenakis used Markovian Stochastic trees to create some of his compositions. Explored in this thesis are seven different methods to compose melodies. After compiling the different melodies, they were categorized by different musical periods based on the musical characteristics found in the melody. This thesis differs from other research that deals …


Multiphoton Response Of Retinal Rod Photoreceptors, Vasilios Alexiades, Harihar Khanal Jan 2007

Multiphoton Response Of Retinal Rod Photoreceptors, Vasilios Alexiades, Harihar Khanal

Publications

Phototransduction is the process by which light is converted into an electrical response in retinal photoreceptors. Rod photoreceptors contain a stack of (about 1000) disc membranes packed with photopigment rhodopsin molecules, which absorb the photons. We present computational experiments which show the profound effect on the response of the distances (how many discs apart) photons happen to be absorbed at. This photon-distribution effect alone can account for much of the observed variability in response.


A Partially Discretized Age-Dependent Population Model With An Additional Stucture, Jean Tchuenche Dec 2006

A Partially Discretized Age-Dependent Population Model With An Additional Stucture, Jean Tchuenche

Applications and Applied Mathematics: An International Journal (AAM)

A semi-discretization method for solving an age-dependent population dynamics model with an additional structure is proposed. This method, unlike previous ones, considers the partial discretization which reduces the model equation into a first order ordinary differential equation. The latter is then solved explicitly and conditions under which second order accuracy arises are given. While the approach adopted is basically analytical, the main result shows that the sum of errors is bounded. An extension to the non-trivial case where growth depends on the additional parameter leads to a Riccati equation, and the existence and
convergence of solutions are proved.


Numerical Simulation Of Nonlinear Elastic Wave Propagation In Piecewise Homogeneous Media, Arkadi Berezovski, Mihhail Berezovski, Juri Engelbrecht Jan 2006

Numerical Simulation Of Nonlinear Elastic Wave Propagation In Piecewise Homogeneous Media, Arkadi Berezovski, Mihhail Berezovski, Juri Engelbrecht

Publications

Systematic experimental work [S. Zhuang, G. Ravichandran, D. Grady, J. Mech. Phys. Solids 51 (2003) 245–265] on laminated composites subjected to high velocity impact loading exhibits the dispersed wave field and the oscillatory behavior of waves with respect to a mean value. Such a behavior is absent in homogeneous solids. An approximate solution to the plate impact in layered heterogeneous solids has been developed in [X. Chen, N. Chandra, A.M. Rajendran, Int. J. Solids Struct. 41 (2004) 4635–4659]. The influence of the particle velocity on many process characteristics was demonstrated. Based on earlier results [A. Berezovski, J. Engelbrecht, G.A. Maugin, …


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.


Transient Non-Linear Heat Conduction Solution By A Dual Reciprocity Boundary Element Method With An Effective Posteriori Error Estimator, Eduardo Divo, Alain J. Kassab Jan 2005

Transient Non-Linear Heat Conduction Solution By A Dual Reciprocity Boundary Element Method With An Effective Posteriori Error Estimator, Eduardo Divo, Alain J. Kassab

Publications

A Dual Reciprocity Boundary Element Method is formulated to solve non-linear heat conduction problems. The approach is based on using the Kirchhoff transform along with lagging of the effective non-linear thermal diffusivity. A posteriori error estimate is used to provide effective estimates of the temporal and spatial error. A numerical example is used to demonstrate the approach.


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.


Response Of Dark-Adapted Retinal Rod Photoreceptors, H. Khanal, V. Alexiades, E. Dibenedetto Jan 2004

Response Of Dark-Adapted Retinal Rod Photoreceptors, H. Khanal, V. Alexiades, E. Dibenedetto

Publications

The process of phototransduction, whereby light is converted into an electrical response, in rod and cone photoreceptors in the retina, involves as a key setp, the diffusion of the cytoplasmic, signaling molecules cGMP (cyclic guanosime monophosphate) and Ca2+ diffuse in the cytoplasm (the fluid surrounding the discs). the complex geometry of the rod creates computational difficulties. We present spatio-temporal compuational models for interacctions and diffusion of cGMP and Ca2+ in the cytoplasm of vertebrate rod photoreceptors, as well as numerical simulations fo the response to light of dark-adapted Salamander rods.


Computational Models For Diffusion Of Second Messengers In Visual Transduction, Harihar Khanal Aug 2003

Computational Models For Diffusion Of Second Messengers In Visual Transduction, Harihar Khanal

Publications

The process of phototransduction, whereby light is converted into an electrical response in retinal rod and cone photoreceptors, involves, as a crucial step, the diffusion of cytoplasmic signaling molecules, termed second messengers. A barrier to mathematical and computational modeling is the complex geometry of the rod outer segment which contains about 1000 thin discs. Most current investigations on the subject assume a well-stirred bulk aqueous environment thereby avoiding such geometrical complexity. We present theoretical and computational spatio-temporal models for phototransduction in vertebrate rod photoreceptors, which are pointwise in nature and thus take into account the complex geometry of the …


Simulation Of Dynamic Electrochemical Processes, John Cassidy Jan 2003

Simulation Of Dynamic Electrochemical Processes, John Cassidy

Articles

This work is designed to introduce electrochemists in a tutorial manner to the basics of modeling of electrochemical systems based primarily on diffusion equations. There is an introduction to analytical and numerical methods with examples taken from typical electrochemical experiments. The Laplace transform is used to derive the Cottrell equation and chronopotentiometry. The response of an electrode to a Gaussian concentration profile is detailed. Laplace’s equation is solved for a simple cell to determine the potential distribution. Discrete methods are employed to calculate the current time behavior following a potential step using the explicit finite difference method. Cyclic voltammetry is …


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.


Droplet Evaporation And Deformations In An Amplitude Modulated Ultrasonic Field, Nihad E. Daidzic, Rene Stadler, Adrian Melling Apr 1996

Droplet Evaporation And Deformations In An Amplitude Modulated Ultrasonic Field, Nihad E. Daidzic, Rene Stadler, Adrian Melling

Aviation Department Publications

The aim of the report presented is the measurements of droplet oscillations.


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.


Electrostatic Positioning Of Droplets In Turbulent Flows (Lstm 375/Te/93), Nihad E. Daidzic, Adrian Melling Apr 1993

Electrostatic Positioning Of Droplets In Turbulent Flows (Lstm 375/Te/93), Nihad E. Daidzic, Adrian Melling

Aviation Department Publications

Report LSTM 375/TE/93, Lehrstuhl fuer Stroemungsmechanik Universitaet Erlangen-Nuernberg Cauerstr. 4, 8520 Erlangen Germany.


Some Stability Problems In Droplet Formation And Breakup (Report Lstm 351/T/92)., Nihad E. Daidzic May 1992

Some Stability Problems In Droplet Formation And Breakup (Report Lstm 351/T/92)., Nihad E. Daidzic

Aviation Department Publications

In this study the instability of droplets and cylindrical jets is investigated. The understanding of these processes has both academic and practical value. Instability of cylindrical jets is theoretically investigated for infinitesimal and finite, but small initial amplitudes (linear and nonlinear stability). For droplets, only linear theory is presented. It is assumed that the capillary force play a dominant role. It is determined that the viscosity exerts a damping effect. In the first section we give an introduction, after which the linear stability theory of cylindrical liquid jets is presented. In the third section the nonlinear jet stability theory is …


A Mathematical Model And Simulation Of Natural Circulation Loop Thermalhydraulics, Nihad E. Daidzic Mar 1992

A Mathematical Model And Simulation Of Natural Circulation Loop Thermalhydraulics, Nihad E. Daidzic

Aviation Department Publications

In this work, a simple, nonlinear, seventh-order mathematical model with lumped parameter analysis of the thermalhydraulic processes in a natural circulation steam boiler loop is presented. Digital simulation has been performed in order to predict open loop transients in severe accidental situations. Parts of the circulation loop such as the steam drum, downcomer, and riser were modelled separately on the basis of the integral balance equations. The two-phase flow in the steam drum and in the riser was modelled with the aid of a homogeneous-equilibrium, two-phase flow model. Special care was taken to ensure the stability of the numerical scheme, …


A Numerical Method For The Solution Of The Schrödinger Equation By A Trial Wavefunction Improvement Formula, Chun-Sheng Ko Jan 1980

A Numerical Method For The Solution Of The Schrödinger Equation By A Trial Wavefunction Improvement Formula, Chun-Sheng Ko

Masters Theses

A numerical method and corresponding computer algorithm for solving the one-dimensional radial Schrödinger equation to any desired accuracy is developed. The method uses a finite difference scheme in which an initial trial wavefunction is digitalized over a lattice covering the region of integration. The values of a rough solution are then altered at each lattice point by a simple improvement formula decreasing the value of the variational energy until the desired minimum is reached. The accuracy of these solutions depends only on the grid size. This method is characterized and tested with a harmonic oscillator potential. Practical evaluations and applications …