Applied Mathematics Commons^™

Hpcnmf: A High-Performance Toolbox For Non-Negative Matrix Factorization, Karthik Devarajan, Guoli Wang

COBRA Preprint Series

Non-negative matrix factorization (NMF) is a widely used machine learning algorithm for dimension reduction of large-scale data. It has found successful applications in a variety of fields such as computational biology, neuroscience, natural language processing, information retrieval, image processing and speech recognition. In bioinformatics, for example, it has been used to extract patterns and profiles from genomic and text-mining data as well as in protein sequence and structure analysis. While the scientific performance of NMF is very promising in dealing with high dimensional data sets and complex data structures, its computational cost is high and sometimes could be critical for …

Minimizing Differences Of Convex Functions With Applications To Facility Location And Clustering, Mau Nam Nguyen, R. Blake Rector, Daniel J. Giles

Mathematics and Statistics Faculty Publications and Presentations

In this paper we develop algorithms to solve generalized Fermat-Torricelli problems with both positive and negative weights and multifacility location problems involving distances generated by Minkowski gauges. We also introduce a new model of clustering based on squared distances to convex sets. Using the Nesterov smoothing technique and an algorithm for minimizing differences of convex functions called the DCA introduced by Tao and An, we develop effective algorithms for solving these problems. We demonstrate the algorithms with a variety of numerical examples.

Single Valued Neutrosophic Graphs, Florentin Smarandache, Said Broumi, Assia Bakali, Mohamed Talea

Branch Mathematics and Statistics Faculty and Staff Publications

The notion of single valued neutrosophic sets is a generalization of fuzzy sets, intuitionistic fuzzy sets. We apply the concept of single valued neutrosophic sets, an instance of neutrosophic sets, to graphs. We introduce certain types of single valued neutrosophic graphs (SVNG) and investigate some of their properties with proofs and examples.

Newton's Law Of Cooling, Caleb J. Emmons

Journal of Humanistic Mathematics

A poem reflecting three different viewpoints on Newton's Law of Cooling.

Signal Flow Graph Approach To Efficient Dst I-Iv Algorithms, Sirani M. Perera

Publications

In this paper, fast and efficient discrete sine transformation (DST) algorithms are presented based on the factorization of sparse, scaled orthogonal, rotation, rotation-reflection, and butterfly matrices. These algorithms are completely recursive and solely based on DST I-IV. The presented algorithms have low arithmetic cost compared to the known fast DST algorithms. Furthermore, the language of signal flow graph representation of digital structures is used to describe these efficient and recursive DST algorithms having (n�1) points signal flow graph for DST-I and n points signal flow graphs for DST II-IV.

Early Afternoon Concurrent Panel Sessions: Commercial Space Industry Snapshot: Presentation: Small Catapult-Assisted Horizontal-Launch Reusable Rbcc Ssto Spaceplane For Economical Short-Duration Leo Access, Nihad E. Daidzic,

Aviation / Aeronautics / Aerospace International Research Conference

This article discusses the conceptual design, flight trajectory calculations, and utilization of the possible future horizontally-launched reusable Single-Stage-to-Orbit (SSTO) spaceplane for small payload short-duration manned/unmanned access to Low-Earth-Orbit (LEO). The 10,000 lb spaceplane would use 5,000 ft catapult-assist horizontal-launch facility and conduct powered approach and landing on conventional horizontal paved runways following the gliding atmospheric re-entry. To increase the economy of operation, the launch facility located at high elevations (4,000+ ft) equatorial region is required, such as, the plateaus in Kenya and Tanzania in Africa and/or Ecuador in South America. A 500-lb payload, including pilot-commander, is envisioned. The propulsion cycle …

Thermally Induced Aggregation Of Rigid Spheres On A Liquid Surface, Eric Forgoston, Leo Hentschker, Siobhan Soltau, Patrick Truitt, Ashuwin Vaidya

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Fluids provide the optimal setting to explore natural patterns far from thermodynamic equilibrium. Experiments suggest that randomly dispersed particles on a liquid surface tend to aggregate on the surface of liquid over time, and the process is enhanced by an increase in the temperature of the liquid. We show that the agglomeration radii increases monotonically with temperature up until the point where all particles in the system form a single, large aggregate. The aggregation dynamics is related to changes in the material properties of the liquid including its viscosity and surface tension as well as the convection driven flow generated …

Parts Of The Whole: Teaching Quantitative Reasoning In The Predator-Prey Model, Dorothy Wallace

Numeracy

The classical predator-prey equations are in nearly every differential equations text and mathematical biology text. Usually they are presented fait accompli, leaving the student to analyze them or play with a computer program. Here we show that the process of fully understanding where these equations come from and how they are derived provides numerous opportunities to teach or reinforce quantitative reasoning skills necessary to future scientists. This example is used to invoke logic, systems thinking, causal reasoning, understanding functions of one or more variables, quantities versus rates of change, proportional reasoning, unit analysis, and comparison to data.

Modelling Subject-Specific Childhood Growth Using Linear Mixed-Effect Models With Cubic Regression Splines, Laura M. Grajeda, Andrada Ivanescu, Mayuko Saito, Ciprian Crainiceanu, Devan Jaganath, Robert H. Gilman, Jean E. Crabtree, Dermott Kelleher, Lilia Cabrera, Vitaliano Cama, William Checkley

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Background: Childhood growth is a cornerstone of pediatric research. Statistical models need to consider individual trajectories to adequately describe growth outcomes. Specifically, well-defined longitudinal models are essential to characterize both population and subject-specific growth. Linear mixed-effect models with cubic regression splines can account for the nonlinearity of growth curves and provide reasonable estimators of population and subject-specific growth, velocity and acceleration. Methods: We provide a stepwise approach that builds from simple to complex models, and account for the intrinsic complexity of the data. We start with standard cubic splines regression models and build up to a model that includes subject-specific …

An Elementary Proof Of Dodgson's Condensation Method For Calculating Determinants, R. Corban Harwood, Mitch Main, Micah Donor

Faculty Publications - Department of Mathematics

In 1866, Charles Ludwidge Dodgson published a paper concerning a method for evaluating determinants called the condensation method. His paper documented a new method to calculate determinants that was based on Jacobi's Theorem. The condensation method is presented and proven here, and is demonstrated by a series of examples. The condensation method can be applied to a number of situations, including calculating eigenvalues, solving a system of linear equations, and even determining the different energy levels of a molecular system. The method is much more efficient than cofactor expansions, particularly for large matrices; for a 5 x 5 matrix, the …

Effect Of Spalled Particles Thermal Degradation On A Hypersonic Flow Field Environment, Raghava S. C. Davuluri, Huaibao Zhang, Alexandre Martin

Mechanical Engineering Faculty Publications

Two-way coupling is performed between a spallation code and a hypersonic aerothermodynamics CFD solver to evaluate the effect of spalled particles on the flow field. Time accurate solutions are computed in argon and air flow fields. A single particle simulations and multiple particles simulations are performed and studied. The results show that the carbon vapor released by spalled particles tend to change the composition of the flow field, particularly the upstream region of the shock.

Models Of Internal Waves In The Presence Of Currents, Alan Compelli, Rossen Ivanov

Conference papers

A fluid system consisting of two domains is examined. The system is considered as being bounded at the bottom and top by a flatbed and wave-free surface respectively. An internal wave propagating in one direction, driven by gravity, acts as a free common interface between the fluids. Various current profiles are considered. The Hamiltonian of the system is determined and expressed in terms of canonical wave-related variables. Limiting behaviour is examined and compared to that of other known models. The linearised equations as well as long-wave approximations are formulated. The presented models provide potential applications to modelling of internal geophysical …

A Competitive Random Sequential Adsorption Model For Immunoassay Activity, Dana Mackey, Eilis Kelly, Robert Nooney

Conference papers

Immunoassays rely on highly specific reactions between antibodies and antigens and are used in biomedical diagnostics applications to detect biomarkers for a variety of diseases. Antibody immobilization to solid interfaces through random adsorption is a widely used technique but has the disadvantage of severely reducing the antigen binding activity and, consequently, the assay performance. This paper proposes a simple mathematical framework, based on the theory known as competitive random sequential adsorption (CRSA), for describing how the activity of immobilized antibodies depends on their orientation and packing density and generalizes a previous model by introducing the antibody aspect ratio as an …

Sphere Representations, Stacked Polytopes, And The Colin De Verdière Number Of A Graph, Lon Mitchell, Lynne Yengulalp

Mathematics Faculty Publications

We prove that a k-tree can be viewed as a subgraph of a special type of (k + 1)- tree that corresponds to a stacked polytope and that these “stacked” (k + 1)-trees admit representations by orthogonal spheres in R k+1. As a result, we derive lower bounds for Colin de Verdi`ere’s µ of complements of partial k-trees and prove that µ(G) + µ(G) > |G| − 2 for all chordal G.

Factorized Runge-Kutta-Chebyshev Methods, Stephen O'Sullivan

Conference papers

The second-order extended stability Factorized Runge-Kutta-Chebyshev (FRKC2) class of explicit schemes for the integration of large systems of PDEs with diffusive terms is presented. FRKC2 schemes are straightforward to implement through ordered sequences of forward Euler steps with complex stepsizes, and easily parallelised for large scale problems on distributed architectures.

Preserving 7 digits for accuracy at 16 digit precision, the schemes are theoretically capable of maintaining internal stability at acceleration factors in excess of 6000 with respect to standard explicit Runge-Kutta methods. The stability domains have approximately the same extents as those of RKC schemes, and are a third longer …

General Solution Of The Wind Triangle Problem And The Critical Tailwind Angle, Nihad E. Daidzic

Aviation Department Publications

A general analytical solution of the navigational wind-triangle problem and the calculation of the critical tailwind angle are presented in this study among other findings. Any crosswind component will effectively create a headwind component on fixed course tracks. The meaning of a route track is lost with excessive crosswinds representing the bifurcation point between the possible and the impossible navigational solutions. Any wind of constant direction and speed will effectively reduce groundspeed and increase time-of-flight on closed-loop multi-segment flights. Effective wind track component consists, in general, of true and induced components. The average groundspeed of multiple-leg flights is a harmonic …

Pattern Recognition In High-Dimensional Data, Matthew Dannenberg

HMC Senior Theses

Vast amounts of data are produced all the time. Yet this data does not easily equate to useful information: extracting information from large amounts of high dimensional data is nontrivial. People are simply drowning in data. A recent and growing source of high-dimensional data is hyperspectral imaging. Hyperspectral images allow for massive amounts of spectral information to be contained in a single image. In this thesis, a robust supervised machine learning algorithm is developed to efficiently perform binary object classification on hyperspectral image data by making use of the geometry of Grassmann manifolds. This algorithm can consistently distinguish between a …

The Global Stability Of The Solution To The Morse Potential In A Catastrophic Regime, Weerapat Pittayakanchit

HMC Senior Theses

Swarms of animals exhibit aggregations whose behavior is a challenge for mathematicians to understand. We analyze this behavior numerically and analytically by using the pairwise interaction model known as the Morse potential. Our goal is to prove the global stability of the candidate local minimizer in 1D found in A Primer of Swarm Equilibria. Using the calculus of variations and eigenvalues analysis, we conclude that the candidate local minimizer is a global minimum with respect to all solution smaller than its support. In addition, we manage to extend the global stability condition to any solutions whose support has a single …

Mathematical Modeling Of Blood Coagulation, Joana L. Perdomo

HMC Senior Theses

Blood coagulation is a series of biochemical reactions that take place to form a blood clot. Abnormalities in coagulation, such as under-clotting or over- clotting, can lead to significant blood loss, cardiac arrest, damage to vital organs, or even death. Thus, understanding quantitatively how blood coagulation works is important in informing clinical decisions about treating deficiencies and disorders. Quantifying blood coagulation is possible through mathematical modeling. This review presents different mathematical models that have been developed in the past 30 years to describe the biochemistry, biophysics, and clinical applications of blood coagulation research. This review includes the strengths and limitations …

4. Dragging Along, Ruth Dover

Differential Equations

More information on air drag.

The Kretschmann Scalar, Charles G. Torre

How to... in 10 minutes or less

On a pseudo-Riemannian manifold with metric g, the "Kretschmann scalar" is a quadratic scalar invariant of the Riemann R tensor of g, defined by contracting all indices with g. In this worksheet we show how to calculate the Kretschmann scalar from a metric.

Steady State Solutions For A System Of Partial Differential Equations Arising From Crime Modeling, Bo Li

HMC Senior Theses

I consider a model for the control of criminality in cities. The model was developed during my REU at UCLA. The model is a system of partial differential equations that simulates the behavior of criminals and where they may accumulate, hot spots. I have proved a prior bounds for the partial differential equations in both one-dimensional and higher dimensional case, which proves the attractiveness and density of criminals in the given area will not be unlimitedly high. In addition, I have found some local bifurcation points in the model.

Topological Data Analysis For Systems Of Coupled Oscillators, Alec Dunton

HMC Senior Theses

Coupled oscillators, such as groups of fireflies or clusters of neurons, are found throughout nature and are frequently modeled in the applied mathematics literature. Earlier work by Kuramoto, Strogatz, and others has led to a deep understanding of the emergent behavior of systems of such oscillators using traditional dynamical systems methods. In this project we outline the application of techniques from topological data analysis to understanding the dynamics of systems of coupled oscillators. This includes the examination of partitions, partial synchronization, and attractors. By looking for clustering in a data space consisting of the phase change of oscillators over a …

1-65-S-Algal Blooms: Algal Blooms Threatening Lake Chapala, R. Corban Harwood

Faculty Publications - Department of Mathematics

This modeling scenario investigates the massive algal blooms that struck Lake Chapala, Mexico, starting in 1994. After reading a summary of articles written on the incidents, students are guided through the process of creating a first order differential equation from a verbal model of the factors and analyze the nonautonomous ODE using direction field, parameter evaluation, and exact solution computation to fully describe the population behavior. Students are expected to be familiar with the separable method and direction fields. Students will learn building and improving a model from qualitative descriptions, nondimensionalization, evaluating parameters, and how to use DFIELD software to …

Oscillation-Free Method For Semilinear Diffusion Equations Under Noisy Initial Conditions, R. Corban Harwood, Likun Zhang, V. S. Manoranjan

Faculty Publications - Department of Mathematics

Noise in initial conditions from measurement errors can create unwanted oscillations which propagate in numerical solutions. We present a technique of prohibiting such oscillation errors when solving initial-boundary-value problems of semilinear diffusion equations. Symmetric Strang splitting is applied to the equation for solving the linear diffusion and nonlinear remainder separately. An oscillation-free scheme is developed for overcoming any oscillatory behavior when numerically solving the linear diffusion portion. To demonstrate the ills of stable oscillations, we compare our method using a weighted implicit Euler scheme to the Crank-Nicolson method. The oscillation-free feature and stability of our method are analyzed through a …

Logistics Of Mathematical Modeling-Focused Projects, R. Corban Harwood

Faculty Publications - Department of Mathematics

Projects provide tangible connections to course content and can motivate students to learn at a deeper level. This article focuses on the implementation of projects in both lower and upper division math courses which develop and analyze mathematical models of a problem based upon known data and real-life situations. Logistical pitfalls and insights are highlighted as well as several key implementation resources. Student feedback demonstrate a positive correlation between the use of projects and an enhanced understanding of the course topics when the impact of logistics is reduced. Best practices learned over the years are given along with example project …

Image Harvest: An Open-Source Platform For High-Throughput Plant Image Processing And Analysis, Avi C. Knecht, Malachy T. Campbell, Adam Caprez, David R. Swanson, Harkamal Walia

Holland Computing Center: Faculty Publications

High-throughput plant phenotyping is an effective approach to bridge the genotype-to-phenotype gap in crops. Phenomics experiments typically result in large-scale image datasets, which are not amenable for processing on desktop computers, thus creating a bottleneck in the image-analysis pipeline. Here, we present an open-source, flexible image-analysis framework, called Image Harvest (IH), for processing images originating from high-throughput plant phenotyping platforms. Image Harvest is developed to perform parallel processing on computing grids and provides an integrated feature for metadata extraction from large-scale file organization. Moreover, the integration of IH with the Open Science Grid provides academic researchers with the computational resources …

A Model For Spheroid Versus Monolayer Response Of Sk-N-Sh Neuroblastoma Cells To Treatment With 15-Deoxy-Pgj2, Dorothy I. Wallace, Ann Dunham, Paula X. Chen, Michelle Chen, Milan Huynh, Evan Rheingold, Olivia F. Prosper

Mathematics Faculty Publications

Researchers have observed that response of tumor cells to treatment varies depending on whether the cells are grown in monolayer, as in vitro spheroids or in vivo. This study uses data from the literature on monolayer treatment of SK-N-SH neuroblastoma cells with 15-deoxy-PGJ₂ and couples it with data on growth rates for untreated SK-N-SH neuroblastoma cells grown as multicellular spheroids. A linear model is constructed for untreated and treated monolayer data sets, which is tuned to growth, death, and cell cycle data for the monolayer case for both control and treatment with 15-deoxy-PGJ₂. The monolayer …

Multiple Problem-Solving Strategies Provide Insight Into Students’ Understanding Of Open-Ended Linear Programming Problems, Marla A. Sole

Publications and Research

Open-ended questions that can be solved using different strategies help students learn and integrate content, and provide teachers with greater insights into students’ unique capabilities and levels of understanding. This article provides a problem that was modified to allow for multiple approaches. Students tended to employ high-powered, complex, familiar solution strategies rather than simpler, more intuitive strategies, which suggests that students might need more experience working with informal solution methods. During the semester, by incorporating open-ended questions, I gained valuable feedback, was able to better model real-world problems, challenge students with different abilities, and strengthen students’ problem solving skills.

An Elliptic Exploration, David Curtis Swart

Online Theses and Dissertations

In this paper I will be giving an introduction to an interesting kind of equation called elliptic curves, and how they can be used to protect our national security through Cryptology. We will explore the unique operation for adding points on elliptic curves and the group structure that it creates, as well as the ECC method, which stands among the RSA and AES methods as one of the modern day's most secure systems of cryptography. In addition, I will also introduce several algorithms and methods that are useful for working with ECC such as Schoof's Algorithm, and I will also …