Numerical Analysis and Computation Commons^™

Adaptive Meshfree Methods For Partial Differential Equations, Jaeyoun Oh

Dissertations

There are many types of adaptive methods that have been developed with different algorithm schemes and definitions for solving Partial Differential Equations (PDE). Adaptive methods have been developed in mesh-based methods, and in recent years, they have been extended by using meshfree methods, such as the Radial Basis Function (RBF) collocation method and the Method of Fundamental Solutions (MFS). The purpose of this dissertation is to introduce an adaptive algorithm with a residual type of error estimator which has not been found in the literature for the adaptive MFS. Some modifications have been made in developing the algorithm schemes depending …

Full Field Computing For Elastic Pulse Dispersion In Inhomogeneous Bars, A. Berezovski, R. Kolman, M. Berezovski, D. Gabriel, V. Adamek

Publications

In the paper, the finite element method and the finite volume method are used in parallel for the simulation of a pulse propagation in periodically layered composites beyond the validity of homogenization methods. The direct numerical integration of a pulse propagation demonstrates dispersion effects and dynamic stress redistribution in physical space on example of a one-dimensional layered bar. Results of numerical simulations are compared with analytical solution constructed specifically for the considered problem. Analytical solution as well as numerical computations show the strong influence of the composition of constituents on the dispersion of a pulse in a heterogeneous bar and …

Masked Instability: Within-Sector Financial Risk In The Presence Of Wealth Inequality, Youngna Choi

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

We investigate masked financial instability caused by wealth inequality. When an economic sector is decomposed into two subsectors that possess a severe wealth inequality, the sector in entirety can look financially stable while the two subsectors possess extreme financially instabilities of opposite nature, one from excessive equity, the other from lack thereof. The unstable subsector can result in further financial distress and even trigger a financial crisis. The market instability indicator, an early warning system derived from dynamical systems applied to agent-based models, is used to analyze the subsectoral financial instabilities. Detailed mathematical analysis is provided to explain what financial …

Transport Phenomena In Field Effect Transistors, Ryan M. Evans, Arvind Balijepalli, Anthony Kearsley

Biology and Medicine Through Mathematics Conference

No abstract provided.

Godunov-Type Upwind Flux Schemes Of The Two-Dimensional Finite Volume Discrete Boltzmann Method, Leitao Chen, Laura Schaefer

Publications

A simple unified Godunov-type upwind approach that does not need Riemann solvers for the flux calculation is developed for the finite volume discrete Boltzmann method (FVDBM) on an unstructured cell-centered triangular mesh. With piecewise-constant (PC), piecewise-linear (PL) and piecewise-parabolic (PP) reconstructions, three Godunov-type upwind flux schemes with different orders of accuracy are subsequently derived. After developing both a semi-implicit time marching scheme tailored for the developed flux schemes, and a versatile boundary treatment that is compatible with all of the flux schemes presented in this paper, numerical tests are conducted on spatial accuracy for several single-phase flow problems. Four major …

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.

Clark Measures And A Theorem Of Ritt, Isabelle Chalendar, Pamela Gorkin, Jonathan R. Partington, William T. Ross

Department of Math & Statistics Faculty Publications

We determine when a finite Blaschke product B can be written, in a non-trivial way, as a composition of two finite Blaschke products (Ritt's problem) in terms of the Clark measure for B. Our tools involve the numerical range of compressed shift operators and the geometry of certain polygons circumscribing the numerical range of the relevant operator. As a consequence of our results, we can determine, in terms of Clark measures, when two finite Blaschke products commute.

Real Solution Of Dae And Pdae System, Zahra Mohammadi, Greg Reid

Western Research Forum

General systems of differential equations don't have restrictions on the number or type of equations. For example, they can be over or under-determined, and also contain algebraic constraints (e.g. algebraic equations such as in Differential-Algebraic equations (DAE) and Partial differential algebraic equations (PDAE). Increasingly such general systems arise from mathematical modeling of engineering and science problems such as in multibody mechanics, electrical circuit design, optimal control, chemical kinetics and chemical control systems. In most applications, only real solutions are of interest, rather than complex-valued solutions. Much progress has been made in exact differential elimination methods, which enable characterization of all …

Numerical Simulation Of Energy Localization In Dynamic Materials, Arkadi Berezovski, Mihhail Berezovski

Publications

Dynamic materials are artificially constructed in such a way that they may vary their characteristic properties in space or in time, or both, by an appropriate arrangement or control. These controlled changes in time can be provided by the application of an external (non-mechanical) field, or through a phase transition. In principle, all materials change their properties with time, but very slowly and smoothly. Changes in properties of dynamic materials should be realized in a short or quasi-nil time lapse and over a sufficiently large material region. Wave propagation is a characteristic feature for dynamic materials because it is also …

Numerical Approximation Of Diffusive Capture Rates By Planar And Spherical Surfaces With Absorbing Pores, Andrew J. Bernoff, Alan E. Lindsay

All HMC Faculty Publications and Research

In 1977 Berg and Purcell published a landmark paper entitled Physics of Chemore- ception, which examined how a bacterium can sense a chemical attractant in the fluid surrounding it [H. C. Berg and E. M. Purcell, Biophys J, 20 (1977), pp. 193–219]. At small scales the attrac- tant molecules move by Brownian motion and diffusive processes dominate. This example is the archetype of diffusive signaling problems where an agent moves via a random walk until it either strikes or eludes a target. Berg and Purcell modeled the target as a sphere with a set of small circular targets (pores) that …

Computing Eigenmodes Of Elliptic Operators On Manifolds Using Radial Basis Functions, Vladimir Delengov

CGU Theses & Dissertations

In this work, a numerical approach based on meshless methods is proposed to obtain eigenmodes of Laplace-Beltrami operator on manifolds, and its performance is compared against existing alternative methods. Radial Basis Function (RBF)-based methods allow one to obtain interpolation and differentiation matrices easily by using scattered data points. We derive expressions for such matrices for the Laplace-Beltrami operator via so-called Reilly’s formulas and use them to solve the respective eigenvalue problem. Numerical studies of proposed methods are performed in order to demonstrate convergence on simple examples of one-dimensional curves and two-dimensional surfaces.

Blow Up Solution Of Bose-Einstein Condensates With Anisotropic Trapping Potential And Rotation, Christopher Leonard

College of Graduate Studies: Theses & Dissertations

In this paper we will analyze a nonlinear Schrodinger equation with quadratic anisotropic trapping potential used to describe a Bose-Einstein condensate. We will find sufficient condition for the global existence of the solution and for the blow up result using physical properties associated to the equation such as the mass, energy, and angular momentum along with some other identities related to the equation. We will finish the thesis with showing some graphical representations describing the solution of the equation.

Evaporation Of A Sessile Droplet On A Slope, Mitch Timm

Dissertations, Master's Theses and Master's Reports

We theoretically examine the drying of a stationary liquid droplet on an inclined surface. Both analytical and numerical approaches are considered, while assuming that the evaporation results from a purely diffusive transport of the liquid vapor and that the contact line is a pinned circle. For the purposes of our analytical calculations, we suppose that the effect of gravity relative to the surface tension is weak, i.e. the Bond number (Bo) is small. Then, we express the shape of the drop and the vapor concentration field as perturbation expansions in terms of Bo. When the Bond number is zero, the …

High Accuracy Methods And Regularization Techniques For Fluid Flows And Fluid-Fluid Interaction, Mustafa Aggul

Dissertations, Master's Theses and Master's Reports

This dissertation contains several approaches to resolve irregularity issues of CFD problems, including a decoupling of non-linearly coupled fluid-fluid interaction, due to high Reynolds number. New models present not only regularize the linear systems but also produce high accurate solutions both in space and time. To achieve this goal, methods solve a computationally attractive artificial viscosity approximation of the target problem, and then utilize a correction approach to make it high order accurate. This way, they all allow the usage of legacy code | a frequent requirement in the simulation of fluid flows in complex geometries. In addition, they all …

An Application Of M-Matrices To Preserve Bounded Positive Solutions To The Evolution Equations Of Biofilm Models, Richard S. Landry Jr.

LSU New Orleans Theses and Dissertations

In this work, we design a linear, two step implicit finite difference method to approximate the solutions of a biological system that describes the interaction between a microbial colony and a surrounding substrate. Three separate models are analyzed, all of which can be described as systems of partial differential equations (PDE)s with nonlinear diffusion and reaction, where the biological colony grows and decays based on the substrate bioavailability. The systems under investigation are all complex models describing the dynamics of biological films. In view of the difficulties to calculate analytical solutions of the models, we design here a numerical technique …

Flow Anisotropy Due To Thread-Like Nanoparticle Agglomerations In Dilute Ferrofluids, Alexander Cali, Wah-Keat Lee, A. David Trubatch, Philip Yecko

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Improved knowledge of the magnetic field dependent flow properties of nanoparticle-based magnetic fluids is critical to the design of biomedical applications, including drug delivery and cell sorting. To probe the rheology of ferrofluid on a sub-millimeter scale, we examine the paths of 550 μm diameter glass spheres falling due to gravity in dilute ferrofluid, imposing a uniform magnetic field at an angle with respect to the vertical. Visualization of the spheres’ trajectories is achieved using high resolution X-ray phase-contrast imaging, allowing measurement of a terminal velocity while simultaneously revealing the formation of an array of long thread-like accumulations of magnetic …

Analytical Solutions For The Black-Scholes Equation, Jalil Manafian, Mahnaz Paknezhad

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, the Black-Sholes equation (BS) has been applied successfully with the Cauchy-Euler method and the method of separation of variables and new analytical solutions have been found. The linear partial differential equation (PDE) transformed to linear ordinary differential equation (ODE) as well. We acquired three types of solutions including hyperbolic, trigonometric and rational solutions. Descriptions of these methods are given and the obtained results reveal that three methods are tools for exploring partial differential models.

Radial Basis Function Differential Quadrature Method For The Numerical Solution Of Partial Differential Equations, Daniel Watson

Dissertations

In the numerical solution of partial differential equations (PDEs), there is a need for solving large scale problems. The Radial Basis Function Differential Quadrature (RBFDQ) method and local RBF-DQ method are applied for the solutions of boundary value problems in annular domains governed by the Poisson equation, inhomogeneous biharmonic equation, and the inhomogeneous Cauchy-Navier equations of elasticity. By choosing the collocation points properly, linear systems can be obtained so that the coefficient matrices have block circulant structures. The resulting systems can be efficiently solved using matrix decomposition algorithms (MDAs) and fast Fourier transforms (FFTs). For the local RBFDQ method, the …

Examining The Electrical Excitation, Calcium Signaling, And Mechanical Contraction Cycle In A Heart Cell, Kristen Deetz, Nygel Foster, Darius Leftwich, Chad Meyer, Shalin Patel, Carlos Barajas, Matthias K. Gobbert, Zana Coulibaly

Spora: A Journal of Biomathematics

As the leading cause of death in the United States, heart disease has become a principal concern in modern society. Cardiac arrhythmias can be caused by a dysregulation of calcium dynamics in cardiomyocytes. Calcium dysregulation, however, is not yet fully understood and is not easily predicted; this provides motivation for the subsequent research. Excitation-contraction coupling (ECC) is the process through which cardiomyocytes undergo contraction from an action potential. Calcium induced calcium release (CICR) is the mechanism through which electrical excitation is coupled with mechanical contraction through calcium signaling. The study of the interplay between electrical excitation, calcium signaling, and mechanical …

Heads And Tails, Julie Simons

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

A Method For Sensitivity Analysis And Parameter Estimation Applied To A Large Reaction-Diffusion Model Of Cell Polarization, Marissa Renardy, Tau-Mu Yi, Dongbin Xiu, Ching-Shan Chou

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Analysis And Implementation Of Numerical Methods For Solving Ordinary Differential Equations, Muhammad Sohel Rana

Masters Theses & Specialist Projects

Numerical methods to solve initial value problems of differential equations progressed quite a bit in the last century. We give a brief summary of how useful numerical methods are for ordinary differential equations of first and higher order. In this thesis both computational and theoretical discussion of the application of numerical methods on differential equations takes place. The thesis consists of an investigation of various categories of numerical methods for the solution of ordinary differential equations including the numerical solution of ordinary differential equations from a number of practical fields such as equations arising in population dynamics and astrophysics. It …

Eignefunctions For Partial Differential Equations On Two-Dimensional Domains With Piecewise Constant Coefficients, Abdullah Muheel Momit Aurko

Master's Theses

In this thesis, we develop a highly accurate and efficient algorithm for computing the solution of a partial differential equation defined on a two-dimensional domain with discontinuous coefficients. An example of such a problem is for modeling the diffusion of heat energy in two space dimensions, in the case where the spatial domain represents a medium consisting of two different but homogeneous materials, with periodic boundary conditions.

Since diffusivity changes based on the material, it will be represented using a piecewise constant function, and this results in the formation of a complicated mathematical model. Such a model is impossible to …

Some Problems Arising From Mathematical Model Of Ductal Carcinoma In Situ., Heng Li

Electronic Theses and Dissertations

Ductal carcinoma in situ (DCIS) is the earliest form of breast cancer. Three mathematical models in the one dimensional case arising from DCIS are proposed. The first two models are in the form of parabolic equation with initial and known moving boundaries. Direct and inverse problems are considered in model 1, existence and uniqueness are proved by using tool from heat potential theory and Volterra integral equations. Also, we discuss the direct problem and nonlocal problem of model 2, existence and uniqueness are proved. And approximation solution of these problems are implemented by Ritz-Galerkin method, which is the first attempt …

Numerical Solution Of Partial Differential Equations Using Polynomial Particular Solutions, Thir R. Dangal

Dissertations

Polynomial particular solutions have been obtained for certain types of partial differential operators without convection terms. In this dissertation, a closed-form particular solution for more general partial differential operators with constant coefficients has been derived for polynomial basis functions. The newly derived particular solutions are further coupled with the method of particular solutions (MPS) for numerically solving a large class of elliptic partial differential equations. In contrast to the use of Chebyshev polynomial basis functions, the proposed approach is more flexible in selecting the collocation points inside the domain. Polynomial basis functions are well-known for yielding ill-conditioned systems when their …

Inverse Problem For A Parabolic System, Reza Pourgholi, Amin Esfahani, Hassan D. Mazraeh

Applications and Applied Mathematics: An International Journal (AAM)

In this paper a numerical approach combining the least squares method and a genetic algorithm is proposed for the determination of the source term in an inverse parabolic system (IPS). A numerical experiment confirm the utility of this algorithm as the results are in good agreement with the exact data. Results show that a reasonable estimation can be obtained by the genetic algorithm within a CPU with clock speed 2.7 GHz.

Finite Difference Schemes For Variable Order Time-Fractional First Initial Boundary Value Problems, Gunvant A. Birajdar, M. M. Rashidi

Applications and Applied Mathematics: An International Journal (AAM)

The aim of the study is to obtain the numerical solution of first initial boundary value problem (IBVP) for semi-linear variable order fractional diffusion equation by using different finite difference schemes. We developed the three finite difference schemes namely explicit difference scheme, implicit difference scheme and Crank-Nicolson difference scheme, respectively for variable order type semi-linear diffusion equation. For this scheme the stability as well as convergence are studied via Fourier method. At the end, solution of some numerical examples are discussed and represented graphically using Matlab.

High Performance Computation Of Cardiac Models In Real-Time Using Webgl, Abouzar Kaboudian, Flavio H. Fenton

Biology and Medicine Through Mathematics Conference

No abstract provided.

A Large Reaction-Diffusion Model For Cell Polarization In Yeast, Marissa Renardy

Biology and Medicine Through Mathematics Conference

No abstract provided.

Efficient Denoising And Sharpening Of Color Images Through Numerical Solution Of Nonlinear Diffusion Equations, Linh T. Duong

Honors Theses

The purpose of this project is to enhance color images through denoising and sharpening, two important branches of image processing, by mathematically modeling the images. Modifications are made to two existing nonlinear diffusion image processing models to adapt them to color images. This is done by treating the red, green, and blue (RGB) channels of color images independently, contrary to the conventional idea that the channels should not be treated independently. A new numerical method is needed to solve our models for high resolution images since current methods are impractical. To produce an efficient method, the solution is represented as …