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

Numerical Analysis and Computation

Series

Institution
Keyword
Publication Year
Publication

Articles 1 - 28 of 28

Full-Text Articles in Ordinary Differential Equations and Applied Dynamics

Game-Theoretic Approaches To Optimal Resource Allocation And Defense Strategies In Herbaceous Plants, Molly R. Creagar Dec 2023

Game-Theoretic Approaches To Optimal Resource Allocation And Defense Strategies In Herbaceous Plants, Molly R. Creagar

Department of Mathematics: Dissertations, Theses, and Student Research

Empirical evidence suggests that the attractiveness of a plant to herbivores can be affected by the investment in defense by neighboring plants, as well as investment in defense by the focal plant. Thus, allocation to defense may not only be influenced by the frequency and intensity of herbivory but also by defense strategies employed by other plants in the environment. We incorporate a neighborhood defense effect by applying spatial evolutionary game theory to optimal resource allocation in plants where cooperators are plants investing in defense and defectors are plants that do not. We use a stochastic dynamic programming model, along …


Lecture 14: Randomized Algorithms For Least Squares Problems, Ilse C.F. Ipsen Apr 2021

Lecture 14: Randomized Algorithms For Least Squares Problems, Ilse C.F. Ipsen

Mathematical Sciences Spring Lecture Series

The emergence of massive data sets, over the past twenty or so years, has lead to the development of Randomized Numerical Linear Algebra. Randomized matrix algorithms perform random sketching and sampling of rows or columns, in order to reduce the problem dimension or compute low-rank approximations. We review randomized algorithms for the solution of least squares/regression problems, based on row sketching from the left, or column sketching from the right. These algorithms tend to be efficient and accurate on matrices that have many more rows than columns. We present probabilistic bounds for the amount of sampling required to achieve a …


Lecture 13: A Low-Rank Factorization Framework For Building Scalable Algebraic Solvers And Preconditioners, X. Sherry Li Apr 2021

Lecture 13: A Low-Rank Factorization Framework For Building Scalable Algebraic Solvers And Preconditioners, X. Sherry Li

Mathematical Sciences Spring Lecture Series

Factorization based preconditioning algorithms, most notably incomplete LU (ILU) factorization, have been shown to be robust and applicable to wide ranges of problems. However, traditional ILU algorithms are not amenable to scalable implementation. In recent years, we have seen a lot of investigations using low-rank compression techniques to build approximate factorizations.
A key to achieving lower complexity is the use of hierarchical matrix algebra, stemming from the H-matrix research. In addition, the multilevel algorithm paradigm provides a good vehicle for a scalable implementation. The goal of this lecture is to give an overview of the various hierarchical matrix formats, such …


Lecture 07: Nonlinear Preconditioning Methods And Applications, Xiao-Chuan Cai Apr 2021

Lecture 07: Nonlinear Preconditioning Methods And Applications, Xiao-Chuan Cai

Mathematical Sciences Spring Lecture Series

We consider solving system of nonlinear algebraic equations arising from the discretization of partial differential equations. Inexact Newton is a popular technique for such problems. When the nonlinearities in the system are well-balanced, Newton's method works well, but when a small number of nonlinear functions in the system are much more nonlinear than the others, Newton may converge slowly or even stagnate. In such a situation, we introduce some nonlinear preconditioners to balance the nonlinearities in the system. The preconditioners are often constructed using a combination of some domain decomposition methods and nonlinear elimination methods. For the nonlinearly preconditioned problem, …


Lecture 02: Tile Low-Rank Methods And Applications (W/Review), David Keyes Apr 2021

Lecture 02: Tile Low-Rank Methods And Applications (W/Review), David Keyes

Mathematical Sciences Spring Lecture Series

As simulation and analytics enter the exascale era, numerical algorithms, particularly implicit solvers that couple vast numbers of degrees of freedom, must span a widening gap between ambitious applications and austere architectures to support them. We present fifteen universals for researchers in scalable solvers: imperatives from computer architecture that scalable solvers must respect, strategies towards achieving them that are currently well established, and additional strategies currently being developed for an effective and efficient exascale software ecosystem. We consider recent generalizations of what it means to “solve” a computational problem, which suggest that we have often been “oversolving” them at the …


Lecture 11: The Road To Exascale And Legacy Software For Dense Linear Algebra, Jack Dongarra Apr 2021

Lecture 11: The Road To Exascale And Legacy Software For Dense Linear Algebra, Jack Dongarra

Mathematical Sciences Spring Lecture Series

In this talk, we will look at the current state of high performance computing and look at the next stage of extreme computing. With extreme computing, there will be fundamental changes in the character of floating point arithmetic and data movement. In this talk, we will look at how extreme-scale computing has caused algorithm and software developers to change their way of thinking on implementing and program-specific applications.


Lecture 01: Scalable Solvers: Universals And Innovations, David Keyes Apr 2021

Lecture 01: Scalable Solvers: Universals And Innovations, David Keyes

Mathematical Sciences Spring Lecture Series

As simulation and analytics enter the exascale era, numerical algorithms, particularly implicit solvers that couple vast numbers of degrees of freedom, must span a widening gap between ambitious applications and austere architectures to support them. We present fifteen universals for researchers in scalable solvers: imperatives from computer architecture that scalable solvers must respect, strategies towards achieving them that are currently well established, and additional strategies currently being developed for an effective and efficient exascale software ecosystem. We consider recent generalizations of what it means to “solve” a computational problem, which suggest that we have often been “oversolving” them at the …


Optimal Allocation Of Two Resources In Annual Plants, David Mcmorris Jul 2020

Optimal Allocation Of Two Resources In Annual Plants, David Mcmorris

Department of Mathematics: Dissertations, Theses, and Student Research

The fitness of an annual plant can be thought of as how much fruit is produced by the end of its growing season. Under the assumption that annual plants grow to maximize fitness, we can use techniques from optimal control theory to understand this process. We introduce two models for resource allocation in annual plants which extend classical work by Iwasa and Roughgarden to a case where both carbohydrates and mineral nutrients are allocated to shoots, roots, and fruits in annual plants. In each case, we use optimal control theory to determine the optimal resource allocation strategy for the plant …


Teaching And Learning Of Fluid Mechanics, Ashwin Vaidya Apr 2020

Teaching And Learning Of Fluid Mechanics, Ashwin Vaidya

Department of Mathematics Facuty Scholarship and Creative Works

Fluid mechanics occupies a privileged position in the sciences; it is taught in various science departments including physics, mathematics, environmental sciences and mechanical, chemical and civil engineering, with each highlighting a different aspect or interpretation of the foundation and applications of fluids. Doll’s fluid analogy [5] for this idea is especially relevant to this issue: “Emergence of creativity from complex flow of knowledge—example of Benard convection pattern as an analogy—dissipation or dispersal of knowledge (complex knowledge) results in emergent structures, i.e., creativity which in the context of education should be thought of as a unique way to arrange information so …


Parameter Estimation And Optimal Design Techniques To Analyze A Mathematical Model In Wound Healing, Nigar Karimli Apr 2019

Parameter Estimation And Optimal Design Techniques To Analyze A Mathematical Model In Wound Healing, Nigar Karimli

Masters Theses & Specialist Projects

For this project, we use a modified version of a previously developed mathematical model, which describes the relationships among matrix metalloproteinases (MMPs), their tissue inhibitors (TIMPs), and extracellular matrix (ECM). Our ultimate goal is to quantify and understand differences in parameter estimates between patients in order to predict future responses and individualize treatment for each patient. By analyzing parameter confidence intervals and confidence and prediction intervals for the state variables, we develop a parameter space reduction algorithm that results in better future response predictions for each individual patient. Moreover, use of another subset selection method, namely Structured Covariance Analysis, that …


A Mathematical Framework On Machine Learning: Theory And Application, Bin Shi Nov 2018

A Mathematical Framework On Machine Learning: Theory And Application, Bin Shi

FIU Electronic Theses and Dissertations

The dissertation addresses the research topics of machine learning outlined below. We developed the theory about traditional first-order algorithms from convex opti- mization and provide new insights in nonconvex objective functions from machine learning. Based on the theory analysis, we designed and developed new algorithms to overcome the difficulty of nonconvex objective and to accelerate the speed to obtain the desired result. In this thesis, we answer the two questions: (1) How to design a step size for gradient descent with random initialization? (2) Can we accelerate the current convex optimization algorithms and improve them into nonconvex objective? For application, …


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

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 …


Simplicity And Sustainability: Pointers From Ethics And Science, Mehrdad Massoudi, Ashwin Vaidya Apr 2018

Simplicity And Sustainability: Pointers From Ethics And Science, Mehrdad Massoudi, Ashwin Vaidya

Department of Mathematics Facuty Scholarship and Creative Works

In this paper, we explore the notion of simplicity. We use definitions of simplicity proposed by philosophers, scientists, and economists. In an age when the rapidly growing human population faces an equally rapidly declining energy/material resources, there is an urgent need to consider various notions of simplicity, collective and individual, which we believe to be a sensible path to restore our planet to a reasonable state of health. Following the logic of mathematicians and physicists, we suggest that simplicity can be related to sustainability. Our efforts must therefore not be spent so much in pursuit of growth but in achieving …


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

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 …


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

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 …


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

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

Department of Mathematics Facuty 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 …


Development Of A Two-Fluid Drag Law For Clustered Particles Using Direct Numerical Simulation And Validation Through Experiments, Ahmadreza Abbasi Baharanchi Nov 2015

Development Of A Two-Fluid Drag Law For Clustered Particles Using Direct Numerical Simulation And Validation Through Experiments, Ahmadreza Abbasi Baharanchi

FIU Electronic Theses and Dissertations

This dissertation focused on development and utilization of numerical and experimental approaches to improve the CFD modeling of fluidization flow of cohesive micron size particles. The specific objectives of this research were: (1) Developing a cluster prediction mechanism applicable to Two-Fluid Modeling (TFM) of gas-solid systems (2) Developing more accurate drag models for Two-Fluid Modeling (TFM) of gas-solid fluidization flow with the presence of cohesive interparticle forces (3) using the developed model to explore the improvement of accuracy of TFM in simulation of fluidization flow of cohesive powders (4) Understanding the causes and influential factor which led to improvements and …


A Class Of High-Order Runge-Kutta-Chebyshev Stability Polynomials, Stephen O'Sullivan Jan 2015

A Class Of High-Order Runge-Kutta-Chebyshev Stability Polynomials, Stephen O'Sullivan

Articles

The analytic form of a new class of factorized Runge-Kutta-Chebyshev (FRKC) stability polynomials of arbitrary order N is presented. Roots of FRKC stability polynomials of degree L = MN are used to construct explicit schemes comprising L forward Euler stages with internal stability ensured through a sequencing algorithm which limits the internal amplification factors to ~ L2. The associated stability domain scales as M2 along the real axis. Marginally stable real-valued points on the interior of the stability domain are removed via a prescribed damping procedure. By construction, FRKC schemes meet all linear order conditions; for nonlinear …


Evolution Of Perturbations In Flow Field Mechanics, Samantha R. Bell, David Forliti, Nils Sedano, Kriss Vanderhyde Jan 2014

Evolution Of Perturbations In Flow Field Mechanics, Samantha R. Bell, David Forliti, Nils Sedano, Kriss Vanderhyde

STAR Program Research Presentations

This project explores the stability analysis of a given flow field. Specifically, where the peak disturbance occurs in a flow as this is the disturbance that is most likely to occur. In rocket combustion, it is important to understand where the maximum disturbance occurs so that the mixing of fuel can be stabilized. The instabilities are the results of frequencies in the area surrounding the flow field. The linear stability governing equations are employed to better understand the disturbance. The governing equations for continuity and momentum in the x and y directions are used to form an equation for the …


Integrability, Recursion Operators And Soliton Interactions, Boyka Aneva, Georgi Grahovski, Rossen Ivanov, Dimitar Mladenov Jan 2014

Integrability, Recursion Operators And Soliton Interactions, Boyka Aneva, Georgi Grahovski, Rossen Ivanov, Dimitar Mladenov

Book chapter/book

This volume contains selected papers based on the talks,presentedat the Conference Integrability, Recursion Operators and Soliton Interactions, held in Sofia, Bulgaria (29-31 August 2012) at the Institute for Nuclear Research and Nuclear Energy of the Bulgarian Academy of Sciences. Included are also invited papers presenting new research developments in the thematic area. The Conference was dedicated to the 65-th birthday of our esteemed colleague and friend Vladimir Gerdjikov. The event brought together more than 30 scientists, from 6 European countries to celebrate Vladimir's scientific achievements. All participants enjoyed a variety of excellent talks in a friendly and stimulating atmosphere. …


Singular Solutions Of Coss-Coupled Epdiff Equations: Waltzing Peakons And Compacton Pairs, Colin Cotter, Darryl Holm, Rossen Ivanov, James Percival Jan 2012

Singular Solutions Of Coss-Coupled Epdiff Equations: Waltzing Peakons And Compacton Pairs, Colin Cotter, Darryl Holm, Rossen Ivanov, James Percival

Conference papers

We introduce EPDiff equations as Euler-Poincare´ equations related to Lagrangian provided by a metric, invariant under the Lie Group Diff(Rn). Then we proceed with a particular form of EPDiff equations, a cross coupled two-component system of Camassa-Holm type. The system has a new type of peakon solutions, 'waltzing' peakons and compacton pairs.


A Generalized Nonlinear Model For The Evolution Of Low Frequency Freak Waves, Jonathan Blackledge Jan 2010

A Generalized Nonlinear Model For The Evolution Of Low Frequency Freak Waves, Jonathan Blackledge

Articles

This paper presents a generalized model for simulating wavefields associated with the sea surface. This includes the case when `freak waves' may occur through an effect compounded in the nonlinear (cubic) Schrodinger equation. After providing brief introductions to linear sea wave models, `freak waves' and the linear and nonlinear Schrodinger equations, we present a unified model that provides for a piecewise continuous transition from a linear to a nonlinear state. This is based on introducing a fractional time derivative to develop a fractional nonlinear partial differential equation with a stochastic source function. In order to explore the characteristics of this …


Approximations Of Continuous Newton's Method: An Extension Of Cayley's Problem, Jon T. Jacobsen, Owen Lewis '05, Bradley Tennis '06 Feb 2007

Approximations Of Continuous Newton's Method: An Extension Of Cayley's Problem, Jon T. Jacobsen, Owen Lewis '05, Bradley Tennis '06

All HMC Faculty Publications and Research

Continuous Newton's Method refers to a certain dynamical system whose associated flow generically tends to the roots of a given polynomial. An Euler approximation of this system, with step size h=1, yields the discrete Newton's method algorithm for finding roots. In this note we contrast Euler approximations with several different approximations of the continuous ODE system and, using computer experiments, consider their impact on the associated fractal basin boundaries of the roots


Reconstruction Of An Unknown Boundary Portion From Cauchy Data In N- Dimensions, Kurt Bryan, Lester Caudill Feb 2005

Reconstruction Of An Unknown Boundary Portion From Cauchy Data In N- Dimensions, Kurt Bryan, Lester Caudill

Department of Math & Statistics Faculty Publications

We consider the inverse problem of determining the shape of some inaccessible 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.


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.


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, …