Applied Mathematics Commons^™

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 …

Application Of Symplectic Integration On A Dynamical System, William Frazier

Electronic Theses and Dissertations

Molecular Dynamics (MD) is the numerical simulation of a large system of interacting molecules, and one of the key components of a MD simulation is the numerical estimation of the solutions to a system of nonlinear differential equations. Such systems are very sensitive to discretization and round-off error, and correspondingly, standard techniques such as Runge-Kutta methods can lead to poor results. However, MD systems are conservative, which means that we can use Hamiltonian mechanics and symplectic transformations (also known as canonical transformations) in analyzing and approximating solutions. This is standard in MD applications, leading to numerical techniques known as symplectic …

Time Lags Associated With Effects Of Oceanic Conditions On Seabird Breeding In The Salish Sea Region Of The Northern California Current System, Rashida S. Smith, Lynelle M. Weldon, James L. Hayward, Shandelle M. Henson

Faculty Publications

No abstract provided.

Shear Driven Micro-Fluidic Pump For Cardiovascular Applications, Nihad E. Daidzic

Aviation Department Publications

A valveless shear-driven micro-fluidic pump design (SDMFP) for hemodynamic applications is presented in this work. One of the possible medical and biomedical applications is in-vivo hemodynamic (human blood circulation) support/assist. One or more SDMFPs can be inserted/implanted into vascular lumens in a form of a stent/duct in series and/or in parallel (bypass duct) to support blood circulation in-vivo. A comprehensive review of various micro-pump designs up to about mid 2000’s is given in [1,2]. Many of micropump designs considered are not suitable for in-vivo or even in-vitro medical/biomedical applications.

Operating principles, design, and SDMFP features are given in [3]. A …

An Rbf Interpolation Blending Scheme For Effective Shock-Capturing, M. Harris, Eduardo Divo, Alain J. Kassab

Publications

In recent years significant focus has been given to the study of Radial basis functions (RBF), especially in their use on solving partial differential equations (PDE). RBF have an impressive capability of inter- polating scattered data, even when this data presents localized discontinuities. However, for infinitely smooth RBF such as the Multiquadrics, inverse Multiquadrics, and Gaussian, the shape parameter must be chosen properly to obtain accurate approximations while avoiding ill-conditioning of the interpolating matrices. The optimum shape parameter can vary significantly depending on the field, particularly in locations of steep gradients, shocks, or discontinuities. Typically, the shape parameter is chosen …

A Coupled Localized Rbf Meshless/Drbem Formulation For Accurate Modeling Of Incompressible Fluid Flows, Leonardo Bueno, Eduardo Divo, Alain J. Kassab

Publications

Velocity-pressure coupling schemes for the solution of incompressible fluid flow problems in Computational Fluid Dynamics (CFD) rely on the formulation of Poisson-like equations through projection methods. The solution of these Poisson-like equations represent the pressure correction and the velocity correction to ensure proper satisfaction of the conservation of mass equation at each step of a time-marching scheme or at each level of an iteration process. Inaccurate solutions of these Poisson-like equations result in meaningless instantaneous or intermediate approximations that do not represent the proper time-accurate behavior of the flow. The fact that these equations must be solved to convergence at …

Haar Wavelet Solution Of Poisson’S Equation And Their Block Structures, The British University In Egypt, Ain Shams University

Basic Science Engineering

The structure of the algebraic system which results from the use of Haar wavelet when solving Poisson’s equation is studied. Haar wavelet technique is used to solve Poisson’s equation on a unit square domain. The form of collocation points are used at the mid points of the subintervals i.e at the odd multiple of the sub interval length labeling. It is proved that the coefficient matrix has symmetric block structure. Comparison with the tridagonal block structure obtained by the finite difference with the natural ordering is introduced. The numerical results have illustrated the superiority of the use of Haar wavelet …

Computational Algorithms For Improved Representation Of The Model Error Covariance In Weak-Constraint 4d-Var, Jeremy A. Shaw

Dissertations and Theses

Four-dimensional variational data assimilation (4D-Var) provides an estimate to the state of a dynamical system through the minimization of a cost functional that measures the distance to a prior state (background) estimate and observations over a time window. The analysis fit to each information input component is determined by the specification of the error covariance matrices in the data assimilation system (DAS). Weak-constraint 4D-Var (w4D-Var) provides a theoretical framework to account for modeling errors in the analysis scheme. In addition to the specification of the background error covariance matrix, the w4D-Var formulation requires information on the model error statistics and …

Application Of An Rbf Blending Interpolation Method To Problems With Shocks, Michael Harris, Eduardo Divo, Alain J. Kassab

Publications

Radial basis functions (RBF) have become an area of research in recent years, especially in the use of solving partial differential equations (PDE). Radial basis functions have an impressive capability in interpolating scattered data, even for data with discontinuities. Although, for infinitely smooth radial basis functions such as the multi-quadrics and inverse multi-quadrics, the shape parameter must be chosen properly to obtain accurate approximations while avoiding ill-conditioning of the interpolating matrices. The optimum shape parameter can vary depending on the field, such as in locations of sharp gradients or shocks. Typically, the shape parameter is chosen to maintain a high …

On The Three Dimensional Interaction Between Flexible Fibers And Fluid Flow, Bogdan Nita, Ryan Allaire

Department of Mathematics Faculty Scholarship and Creative Works

In this paper we discuss the deformation of a flexible fiber clamped to a spherical body and immersed in a flow of fluid moving with a speed ranging between 0 and 50 cm/s by means of three dimensional numerical simulation developed in COMSOL . The effects of flow speed and initial configuration angle of the fiber relative to the flow are analyzed. A rigorous analysis of the numerical procedure is performed and our code is benchmarked against well established cases. The flow velocity and pressure are used to compute drag forces upon the fiber. Of particular interest is the behavior …

Steady State Probabilities In Relation To Eigenvalues, Pellegrino Christopher

The Kabod

By using the methods of Hamdy Taha, eigenvectors can be used in solving problems to compute steady state probabilities, and they work every time.

Numerical Simulation Of Acoustic Emission During Crack Growth In 3-Point Bending Test, Mihhail Berezovski, Arkadi Berezovski

Publications

Numerical simulation of acoustic emission by crack propagation in 3-point bending tests is performed to investigate how the interaction of elastic waves generates a detectable signal. It is shown that the use of a kinetic relation for the crack tip velocity combined with a simple crack growth criterion provides the formation of waveforms similar to those observed in experiments.

Weihrauch Reducibility And Finite-Dimensional Subspaces, Sean Sovine

Theses, Dissertations and Capstones

In this thesis we study several principles involving subspaces and decompositions of vector spaces, matroids, and graphs from the perspective of Weihrauch reducibility. We study the problem of decomposing a countable vector space or countable matroid into 1-dimensional subspaces. We also study the problem of producing a finite-dimensional or 1-dimensional subspace of a countable vector space, and related problems for producing finite-dimensional subspaces of a countable matroid. This extends work in the reverse mathematics setting by Downey, Hirschfeldt, Kach, Lempp, Mileti, and Montalb´an (2007) and recent work of Hirst and Mummert (2017). Finally, we study the problem of producing a …

Kinetic Monte Carlo Methods For Computing First Capture Time Distributions In Models Of Diffusive Absorption, Daniel Schmidt

HMC Senior Theses

In this paper, we consider the capture dynamics of a particle undergoing a random walk above a sheet of absorbing traps. In particular, we seek to characterize the distribution in time from when the particle is released to when it is absorbed. This problem is motivated by the study of lymphocytes in the human blood stream; for a particle near the surface of a lymphocyte, how long will it take for the particle to be captured? We model this problem as a diffusive process with a mixture of reflecting and absorbing boundary conditions. The model is analyzed from two approaches. …

Dynamics Of Gene Networks In Cancer Research, Paul Scott

College of Graduate Studies: Theses & Dissertations

Cancer prevention treatments are being researched to see if an optimized treatment schedule would decrease the likelihood of a person being diagnosed with cancer. To do this we are looking at genes involved in the cell cycle and how they interact with one another. Through each gene expression during the life of a normal cell we get an understanding of the gene interactions and test these against those of a cancerous cell. First we construct a simplified network model of the normal gene network. Once we have this model we translate it into a transition matrix and force changes on …

Neural Network Predictions Of A Simulation-Based Statistical And Graph Theoretic Study Of The Board Game Risk, Jacob Munson

Murray State Theses and Dissertations

We translate the RISK board into a graph which undergoes updates as the game advances. The dissection of the game into a network model in discrete time is a novel approach to examining RISK. A review of the existing statistical findings of skirmishes in RISK is provided. The graphical changes are accompanied by an examination of the statistical properties of RISK. The game is modeled as a discrete time dynamic network graph, with the various features of the game modeled as properties of the network at a given time. As the network is computationally intensive to implement, results are produced …

A Numerical Study Of Construction Of Honey Bee Comb, Pamela Guerrero, Pamela C. Guerrero

Murray State Theses and Dissertations

We use finite difference methods in the treatment of an existing system of partial differential equations that captures the dynamics of parallel honeycomb construction in a bee hive. We conduct an uncertainty analysis by calculating the partial rank correlation coefficient for the parameters to find which are most important to the outcomes of the model. We then use an eFAST method to determine both the individual and total sensitivity index for the parameters. Afterwards we examine our numerical model under varying initial conditions and parameter values, and compare ratios found from local data with the golden mean by fitting images …

A Physics-Based Approach To Modeling Wildland Fire Spread Through Porous Fuel Beds, Tingting Tang

Theses and Dissertations--Mechanical and Aerospace Engineering

Wildfires are becoming increasingly erratic nowadays at least in part because of climate change. CFD (computational fluid dynamics)-based models with the potential of simulating extreme behaviors are gaining increasing attention as a means to predict such behavior in order to aid firefighting efforts. This dissertation describes a wildfire model based on the current understanding of wildfire physics. The model includes physics of turbulence, inhomogeneous porous fuel beds, heat release, ignition, and firebrands. A discrete dynamical system for flow in porous media is derived and incorporated into the subgrid-scale model for synthetic-velocity large-eddy simulation (LES), and a general porosity-permeability model is …

Long And Short-Range Air Navigation On Spherical Earth, Nihad E. Daidzic

International Journal of Aviation, Aeronautics, and Aerospace

Global range air navigation implies non-stop flight between any two airports on Earth. Such effort would require airplanes with the operational air range of at least 12,500 NM which is about 40-60% longer than anything existing in commercial air transport today. Air transportation economy requires flying shortest distance, which in the case of spherical Earth are Orthodrome arcs. Rhumb-line navigation has little practical use in long-range flights, but has been presented for historical reasons and for comparison. Database of about 50 major international airports from every corner of the world has been designed and used in testing and route validation. …

Paving The Randomized Gauss-Seidel, Wei Wu

Scripps Senior Theses

The Randomized Gauss-Seidel Method (RGS) is an iterative algorithm that solves overdetermined systems of linear equations Ax = b. This paper studies an update on the RGS method, the Randomized Block Gauss-Seidel Method. At each step, the algorithm greedily minimizes the objective function L(x) = kAx bk2 with respect to a subset of coordinates. This paper describes a Randomized Block Gauss-Seidel Method (RBGS) which uses a randomized control method to choose a subset at each step. This algorithm is the first block RGS method with an expected linear convergence rate which can be described by the properties of the matrix …

Modeling And Simulation Of The Peristaltic Flow Of Newtonian And Non-Newtonian Fluids With Application To The Human Body, Samer Alokaily

Dissertations, Master's Theses and Master's Reports

Computational models are developed to investigate peristaltic motion in the human gastro-intestinal tract. The peristaltic motion is simulated by means of traveling waves which deform the boundary of the tubes. An axisymmetric tube of uniform diameter is used to model the small intestines, and an axisymmetric conical geometry is developed to model the lower part of the human stomach. The conical geometry represents a simplification of the more complicated three-dimensional models of the human stomach that have been used in other studies. Also, they seeks to reduce computational costs and circumvent difficulties of mesh generation. The computations are performed within …

Numerical Simulation Of Viscoelastic Multiphase Flows Using An Improved Two-Phase Flow Solver, Olabanji Shonibare

Dissertations, Master's Theses and Master's Reports

The production of uniformly-sized droplets has numerous applications in various fields including the biotechnology and chemical industries. For example, in the separation of mixtures based on their relative absorbency, an optimal arrangement of monodispersed droplets in columns is desired for an effective separation. However, very few numerical studies on the formation of viscoelastic droplets via cross-flow shear are available, none of which have considered the case when the flow of the continuous phase is Couette. In this work, a new solver capable of automatic mesh refinement is developed for the OpenFOAM CFD toolbox to solve viscoelastic two-phase flow problems. The …

Daily Traffic Flow Pattern Recognition By Spectral Clustering, Matthew Aven

CMC Senior Theses

This paper explores the potential applications of existing spectral clustering algorithms to real life problems through experiments on existing road traffic data. The analysis begins with an overview of previous unsupervised machine learning techniques and constructs an effective spectral clustering algorithm that demonstrates the analytical power of the method. The paper focuses on the spectral embedding method’s ability to project non-linearly separable, high dimensional data into a more manageable space that allows for accurate clustering. The key step in this method involves solving a normalized eigenvector problem in order to construct an optimal representation of the original data.

While this …

Direct Sampling Methods For Inverse Scattering Problems, Ala Mahmood Nahar Al Zaalig

Dissertations, Master's Theses and Master's Reports

Recently, direct sampling methods became popular for solving inverse scattering problems to estimate the shape of the scattering object. They provide a simple tool to directly reconstruct the shape of the unknown scatterer. These methods are based on choosing an appropriate indicator function f on R^d, d=2 or 3, such that f(z) decides whether z lies inside or outside the scatterer. Consequently, we can determine the location and the shape of the unknown scatterer.

In this thesis, we first present some sampling methods for shape reconstruction in inverse scattering problems. These methods, which are described in Chapter 1, …

An Introduction To The Lebesgue Integral, Ikhlas Adi

All Student Theses and Dissertations

The Riemann integral is the simplest integral to define, and it allows one to integrate every continuous function. It is really important to have a definition of the integral that allows a wider class of functions to be integrated. However, there are many other types of integrals, the most important of which is the Lebesgue integral. The Lebesgue integral allows one to integrate unbounded or discontinuous functions whose Riemann integral does not exist, and it has mathematical properties that the Riemann integral does not. The definition of the Lebesgue integral requires the use of measure theory since picking out a …

Fitting A Linear Regression Model And Forecasting In R In The Presence Of Heteroskedascity With Particular Reference To Advanced Regression Technique Dataset On Kaggle.Com., Samuel Mbah Nde

All Student Theses and Dissertations

Since ancient times, men have built and sold houses. But just how much is a house worth? The challenge is to be able to use information about a house such as its location, and the area on which it is built to predict its price. Such predicted prices can be of great importance to any participant in the real estate business be it an agent, a buyer, seller or a bank to make intelligent decisions and the profit that come with such decisions. Since every company’s success depends on its ability to accurately predict financial outcomes, its profitability will depend …

Computational Fluid Dynamics In A Terminal Alveolated Bronchiole Duct With Expanding Walls: Proof-Of-Concept In Openfoam, Jeremy Myers

Theses and Dissertations

Mathematical Biology has found recent success applying Computational Fluid Dynamics (CFD) to model airflow in the human lung. Detailed modeling of flow patterns in the alveoli, where the oxygen-carbon dioxide gas exchange occurs, has provided data that is useful in treating illnesses and designing drug-delivery systems. Unfortunately, many CFD software packages have high licensing fees that are out of reach for independent researchers. This thesis uses three open-source software packages, Gmsh, OpenFOAM, and ParaView, to design a mesh, create a simulation, and visualize the results of an idealized terminal alveolar sac model. This model successfully demonstrates that OpenFOAM can be …

Multi-Scale Cardiovascular Flow Analysis By An Integrated Meshless-Lumped Parameter Model, Leonardo A. Bueno, Eduardo A. Divo, Alain J. Kassab

Publications

A computational tool that integrates a Radial basis function (RBF)-based Meshless solver with a Lumped Parameter model (LPM) is developed to analyze the multi-scale and multi-physics interaction between the cardiovascular flow hemodynamics, the cardiac function, and the peripheral circulation. The Meshless solver is based on localized RBF collocations at scattered data points which allows for automation of the model generation via CAD integration. The time-accurate incompressible flow hemodynamics are addressed via a pressure-velocity correction scheme where the ensuing Poisson equations are accurately and efficiently solved at each time step by a Dual-Reciprocity Boundary Element method (DRBEM) formulation that takes advantage …

Estimates Of Life Span Of Solutions Of A Cauchy Problem, Joon Hyuk Kang

Faculty Publications

In this paper we get estimates of life span of a Cauchy problem ut(x, t) = ∆ u(x, t) +u(x, t)p, x∈Rn, t >0,u(x,0) =λφ(x), x∈Rn in terms of the positive constant parameterλ whenφ(x)∈Lq is a nonnegative bounded continuous function in Rn but not identically zero, where q is large enough. The technique we used in this paper is the Comparison Principle.

Region Of Smooth Functions For Positive Solutions To An Elliptic Biological Model, Joon Hyuk Kang, Timothy Robertson

Faculty Publications

The non-existence and existence of the positive solution to the generalized elliptic model ∆u+g(u v) = 0 in Ω, ∆v+h(u, v) = 0 in Ω, u=v= 0 on∂Ω, were investigated.