Physical Sciences and Mathematics Commons^™

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

University of 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 …

Feasible Computation In Symbolic And Numeric Integration, Robert H.C. Moir

Electronic Thesis and Dissertation Repository

Two central concerns in scientific computing are the reliability and efficiency of algorithms. We introduce the term feasible computation to describe algorithms that are reliable and efficient given the contextual constraints imposed in practice. The main focus of this dissertation then, is to bring greater clarity to the forms of error introduced in computation and modeling, and in the limited context of symbolic and numeric integration, to contribute to integration algorithms that better account for error while providing results efficiently.

Chapter 2 considers the problem of spurious discontinuities in the symbolic integration problem, proposing a new method to restore continuity …

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 …

Rate Of Convergence For Generalized Szász–Mirakyan Operators In Exponential Weighted Space, Sevilay K. Serenbay, Özge Dalmano˘Glu

Applications and Applied Mathematics: An International Journal (AAM)

In the present paper, generalized Szász–Mirakyan operators in exponential weighted space of functions of one variable are introduced. Using a method given by Rempulska and Walczak, some theorems on the degree of approximation are investigated. Furthermore, a numerical example with an illustrative graphic is given to show comparison for the error estimates of the operators.

Evaluation Of Some Reliability Characteristics Of A Single Unit System Requiring Two Types Of Supporting Device For Operations, Ibrahim Yusuf, Nura J. Fagge

Applications and Applied Mathematics: An International Journal (AAM)

This study presents the reliability assessment of a single unit connected to two types of external supporting devices for its operation. Each type of external supporting device has two copies I and II on standby. First order differential difference equations method is used to obtain the explicit expression for the steady state availability, busy period due to failure of type I and II supporting devices of repairmen, steady-state availability and profit function. Based on assumed numerical values given to system parameters, graphical illustrations are given to highlight important results. Comparisons are performed to highlight the impact of unit failure and …

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.

A New Hybrid Method For Solving Nonlinear Fractional Differential Equations, R. Delpasand, M. M. Hosseini, F. M. Maalek Ghaini

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, numerical solution of initial and boundary value problems for nonlinear fractional differential equations is considered by pseudospectral method. In order to avoid solving systems of nonlinear equations resulting from the method, the residual function of the problem is constructed, as well as a suggested unconstrained optimization model solved by PSOGSA algorithm. Furthermore, the research inspects and discusses the spectral accuracy of Chebyshev polynomials in the approximation theory. The following scheme is tested for a number of prominent examples, and the obtained results demonstrate the accuracy and efficiency of the proposed method.

Numerical Experiments For Finding Roots Of The Polynomials In Chebyshev Basis, M. S. Solary

Applications and Applied Mathematics: An International Journal (AAM)

Root finding for a function or a polynomial that is smooth on the interval [a; b], but otherwise arbitrary, is done by the following procedure. First, approximate it by a Chebyshev polynomial series. Second, find the zeros of the truncated Chebyshev series. Finding roots of the Chebyshev polynomial is done by eigenvalues of a nXn matrix such as companion or comrade matrices. There are some methods for finding eigenvalues of these matrices such as companion matrix and chasing procedures.We derive another algorithm by second kind of Chebyshev polynomials.We computed the numerical results of these methods for some special and ill-conditioned …

Effective Modified Hybrid Conjugate Gradient Method For Large-Scale Symmetric Nonlinear Equations, Jamilu Sabi'u, Mohammed Y. Waziri

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, we proposed hybrid conjugate gradient method using the convex combination of FR and PRP conjugate gradient methods for solving Large-scale symmetric nonlinear equations via Andrei approach with nonmonotone line search. Logical formula for obtaining the convex parameter using Newton and our proposed directions was also proposed. Under appropriate conditions global convergence was established. Reported numerical results show that the proposed method is very promising.

The Fmx/Fm/1 Queue With Multiple Working Vacation, G. Kannadasan, N. Sathiyamoorthi

Applications and Applied Mathematics: An International Journal (AAM)

This study investigates the batch arrival FM^X/FM/1 queue with multiple working vacation. For this fuzzy queuing model, this research obtains some performance measure of interest such as mean system length, mean system sojourn time, mean busy period for the server and working vacation period. Finally, numerical results are presented to show the effects of system parameters.

Interactions Of Thermoelastic Beam In Modified Couple Stress Theory, Rajneesh Kumar, Shaloo Devi

Applications and Applied Mathematics: An International Journal (AAM)

This paper is concerned with the study of thermoelastic beam in modified couple stress theory. The governing equations of motion for modified couple stress theory and heat conduction equation for non-Fourier (non-classical process) are investigated to model the vibrations in a homogeneous isotropic thin beam in a closed form by employing the Euler Bernoulli beam theory. The generalized theories of thermoelasticity with one and two relaxation times are used to model the problem. Both ends of the beam are simply supported. The Laplace transform technique applied to solve the system of equations which are written in dimensionless form. A general …

Analytical Solution For Determination The Control Parameter In The Inverse Parabolic Equation Using Ham, M. M. Khader

Applications and Applied Mathematics: An International Journal (AAM)

In this article, the homotopy analysis method (HAM) for obtaining the analytical solution of the inverse parabolic problem and computing the unknown time-dependent parameter is introduced. The series solution is developed and the recurrence relations are given explicitly. Special attention is given to satisfy the convergence of the proposed method. A comparison of HAM with the variational iteration method is made. In the HAM, we use the auxiliary parameter ~ to control with a simple way in the convergence region of the solution series. Applying this method with several

Mathematical Modeling Of Mixtures And Numerical Solution With Applications To Polymer Physics, John Timothy Cummings

Doctoral Dissertations

We consider in this dissertation the mathematical modeling and simulation of a general diffuse interface mixture model based on the principles of energy dissipation. The model developed allows for a thermodynamically consistent description of systems with an arbitrary number of different components, each of which having perhaps differing densities. We also provide a mathematical description of processes which may allow components to source or sink into other components in a mass conserving, energy dissipating way, with the motivation of applying this model to phase transformation. Also included in the modeling is a unique set of thermodynamically consistent boundary conditions which …

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 …

Stochastic Analysis Of A Mammalian Circadian Clock Model: Small Protein Number Effects, David W. Morgens, Blerta Shtylla

Spora: A Journal of Biomathematics

The circadian clock, responsible for coordinating organism function with daily and seasonal changes in the day-night cycle, is controlled by a complex protein network that constitutes a robust biochemical oscillator. Deterministic ordinary differential equation models have been used extensively to model the behavior of these central clocks. However, due to the small number of proteins involved in the circadian oscillations, mathematical models that track stochastic variations in the numbers of clock proteins may reveal more complex and biologically relevant behaviors. In this paper, we compare the response of a robust yet detailed deterministic model for the mammalian circadian clock with …

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 …

Non-Equispaced Fast Fourier Transforms In Turbulence Simulation, Aditya M. Kulkarni

Masters Theses

Fourier pseudo-spectral method on equispaced grid is one of the approaches in turbulence simulation, to compute derivative of discrete data, using fast Fourier Transform (FFT) and gives low dispersion and dissipation errors. In many turbulent flows the dynamically important scales of motion are concentrated in certain regions which requires a coarser grid for higher accuracy. A coarser grid in other regions minimizes the memory requirement. This requires the use of Non-equispaced Fast Fourier Transform (NFFT) to compute the Fourier transform, by solving a system of linear equations.

To achieve similar accuracy, the NFFT needs to return more Fourier coefficients than …

Heads And Tails, Julie Simons

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

A Cellular Automaton Modeling Approach To Chestnut Blight Canker Development, Samuel Iselin

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Distributed Evolution Of Spiking Neuron Models On Apache Mahout For Time Series Analysis, Andrew Palumbo

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.

High-Order Relaxed Multirate Infinitesimal Step Methods For Multiphysics Applications, Jean Sexton

Mathematics Theses and Dissertations

In this work, we consider numerical methods for integrating multirate ordinary differential equations. We are interested in the development of new multirate methods with good stability properties and improved efficiency over existing methods. We discuss the development of multirate methods, particularly focusing on those that are based on Runge-Kutta theory. We introduce the theory of Generalized Additive Runge-Kutta methods proposed by Sandu and Günther. We also introduce the theory of Recursive Flux Splitting Multirate Methods with Sub-cycling described by Schlegel, as well as the Multirate Infinitesimal Step methods this work is based on. We propose a generic structure called Flexible …

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 …

Filtered Subspace Iteration For Selfadjoint Operators, Jay Gopalakrishnan, Luka Grubišić, Jeffrey S. Ovall

Portland Institute for Computational Science Publications

We consider the problem of computing a cluster of eigenvalues (and its associated eigenspace) of a (possibly unbounded) selfadjoint operator in a Hilbert space. A rational function of the operator is constructed such that the eigenspace of interest is its dominant eigenspace, and a subspace iteration procedure is used to approximate this eigenspace. The computed space is then used to obtain approximations of the eigenvalues of interest. An eigenvalue and eigenspace convergence analysis that considers both iteration error and dis- cretization error is provided. A realization of the proposed approach for a model second-order elliptic operator is based on a …

Convergence Of Media Attention Across 129 Countries, Jisun An, Hassan Aldarbesti, Haewoon Kwak

Research Collection School Of Computing and Information Systems

The objective of this study is to assess the longitudinal trends of media similarity and dissimilarity on the international scale. As news value has well-established political, cultural, and economic consequences, the degree to which media coverage and content is converging across countries has implications for international relations. To study this convergence, we use the daily data of the 100 topics that were over-reported in each country, compared to other countries, from March 7 to October 9, 2016. The results of this analysis indicate that two complementary patterns–globalization and domestication–explain the media attention across the countries. We conclude that this attention …

Information Theoretic Study Of Gaussian Graphical Models And Their Applications, Ali Moharrer

LSU Doctoral Dissertations

In many problems we are dealing with characterizing a behavior of a complex stochastic system or its response to a set of particular inputs. Such problems span over several topics such as machine learning, complex networks, e.g., social or communication networks; biology, etc. Probabilistic graphical models (PGMs) are powerful tools that offer a compact modeling of complex systems. They are designed to capture the random behavior, i.e., the joint distribution of the system to the best possible accuracy. Our goal is to study certain algebraic and topological properties of a special class of graphical models, known as Gaussian graphs. First, …

On The Ramberg-Osgood Stress-Strain Model And Large Deformations Of Cantilever Beams, Ronald J. Giardina Jr

University of New Orleans Theses and Dissertations

In this thesis the Ramberg-Osgood nonlinear model for describing the behavior of many different materials is investigated. A brief overview of the model as it is currently used in the literature is undertaken and several misunderstandings and possible pitfalls in its application is pointed out, especially as it pertains to more recent approaches to finding solutions involving the model. There is an investigation of the displacement of a cantilever beam under a combined loading consisting of a distributed load across the entire length of the beam and a point load at its end and new solutions to this problem are …

Large-Scale Online Feature Selection For Ultra-High Dimensional Sparse Data, Yue Wu, Steven C. H. Hoi, Tao Mei, Nenghai Yu

Research Collection School Of Computing and Information Systems

Feature selection (FS) is an important technique in machine learning and data mining, especially for large scale high-dimensional data. Most existing studies have been restricted to batch learning, which is often inefficient and poorly scalable when handling big data in real world. As real data may arrive sequentially and continuously, batch learning has to retrain the model for the new coming data, which is very computationally intensive. Online feature selection (OFS) is a promising new paradigm that is more efficient and scalable than batch learning algorithms. However, existing online algorithms usually fall short in their inferior efficacy. In this article, …

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 …

Joint Inversion Of Compact Operators, James Ford

Boise State University Theses and Dissertations

The first mention of joint inversion came in [22], where the authors used the singular value decomposition to determine the degree of ill-conditioning in inverse problems. The authors demonstrated in several examples that combining two models in a joint inversion, and effectively stacking discrete linear models, improved the conditioning of the problem. This thesis extends the notion of using the singular value decomposition to determine the conditioning of discrete joint inversion to using the singular value expansion to determine the well-posedness of joint linear operators. We focus on compact linear operators related to geophysical, electromagnetic subsurface imaging.

The operators are …