Numerical Analysis and Computation Commons^™

Sloane’S Gap: Do Mathematical And Social Factors Explain The Distribution Of Numbers In The Oeis?, Nicolas J.-P. Gauvrit, Jean-Paul Delahaye, Hector Zenil

Journal of Humanistic Mathematics

The Online Encyclopedia of Integer Sequences (OEIS) is a catalog of integer sequences. We are particularly interested in the number of occurrences of N(n) of an integer n in the database. This number N(n) marks the importance of n and it varies noticeably from one number to another, and from one number to the next in a series. “Importance” can be mathematically objective (2^10 is an example of an “important” number in this sense) or as the result of a shared mathematical culture (10^9 is more important than 9^10 because we use a decimal notation). The concept of algorithmic complexity …

Generalized Least-Squares Regressions Ii: Theory And Classification, Nataniel Greene

Publications and Research

In the first paper of this series, a variety of known and new symmetric and weighted least-squares regression methods were presented with efficient derivations. This paper continues and generalizes the previous work with a theory for deriving, analyzing, and classifying all symmetric and weighted least-squares regression methods.

Generalized Least-Squares Regressions I: Efficient Derivations, Nataniel Greene

Publications and Research

Ordinary least-squares regression suffers from a fundamental lack of symmetry: the regression line of y given x and the regression line of x given y are not inverses of each other. Alternative symmetric regression methods have been developed to address this concern, notably: orthogonal regression and geometric mean regression. This paper presents in detail a variety of least squares regression methods which may not have been known or fully explicated. The derivation of each method is made efficient through the use of Ehrenberg's formula for the ordinary least-squares error and through the extraction of a weight function g(b) which characterizes …

Near-Optimal Compressed Sensing Guarantees For Anisotropic And Isotropic Total Variation Minimization, Deanna Needell, Rachel Ward

CMC Faculty Publications and Research

Consider the problem of reconstructing a multidimensional signal from partial information, as in the setting of compressed sensing. Without any additional assumptions, this problem is ill-posed. However, for signals such as natural images or movies, the minimal total variation estimate consistent with the measurements often produces a good approximation to the underlying signal, even if the number of measurements is far smaller than the ambient dimensionality. Recently, guarantees for two-dimensional images were established. This paper extends these theoretical results to signals of arbitrary dimension and to both the anisotropic and isotropic total variation problems. To be precise, we show that …

Simulations Of Surfactant Driven Thin Film Flow, Shreyas Kumar

HMC Senior Theses

This thesis is intended to fulfill the requirements of the Math and Physics departments at Harvey Mudd College. We begin with a brief introduction to the study of surfactant dynamics followed by some background on the experimental framework our work is related to. We then go through a derivation of the model we use, and explore in depth the nature of the Equation of State (EoS), the relationship between the surface tension on a fluid and the surfactant concentration. We consider the effect of using an empirical equation of state on the results of the simulations and compare the new …

Two-Dimensional Hydrodynamic Modeling Of Two-Phase Flow For Understanding Geyser Phenomena In Urban Stormwater System, Zhiyu S. Shao

Theses and Dissertations--Civil Engineering

During intense rain events a stormwater system can fill rapidly and undergo a transition from open channel flow to pressurized flow. This transition can create large discrete pockets of trapped air in the system. These pockets are pressurized in the horizontal reaches of the system and then are released through vertical vents. In extreme cases, the transition and release of air pockets can create a geyser feature.

The current models are inadequate for simulating mixed flows with complicated air-water interactions, such as geysers. Additionally, the simulation of air escaping in the vertical dropshaft is greatly simplified, or completely ignored, in …

Time-Stepping For Laser Ablation, Harihar Khanal, David Autrique, Vasilios Alexiades

Publications

Nanosecond laser ablation is a popular technique, applied in many areas of science and technology such as medicine, archaeology, chemistry, environmental and materials sciences. We outline a computational model for radiative and collisional processes occurring during ns-laser ablation, and compare the performance of various low and high order time-stepping algorithms.

Hydrodynamic Modeling Of Ns-Laser Ablation, David Autrique, Vasilios Alexiades, Harihar Khanal

Publications

Laser ablation is a versatile and widespread technique, applied in an increasing number of medical, industrial and analytical applications. A hydrodynamic multiphase model describing nanosecond-laser ablation (ns- LA) is outlined. The model accounts for target heating and mass removal mechanisms as well as plume expansion and plasma formation. A copper target is placed in an ambient environment consisting of helium and irradiated by a nanosecond-laser pulse. The effect of variable laser settings on the ablation process is explored in 1-D numerical simulations.

Sampling From The Hardcore Process, William C. Dodds

CMC Senior Theses

Partially Recursive Acceptance Rejection (PRAR) and bounding chains used in conjunction with coupling from the past (CFTP) are two perfect simulation protocols which can be used to sample from a variety of unnormalized target distributions. This paper first examines and then implements these two protocols to sample from the hardcore gas process. We empirically determine the subset of the hardcore process's parameters for which these two algorithms run in polynomial time. Comparing the efficiency of these two algorithms, we find that PRAR runs much faster for small values of the hardcore process's parameter whereas the bounding chain approach is vastly …

Pressure Poisson Method For The Incompressible Navier-Stokes Equations Using Galerkin Finite Elements, John Cornthwaite

College of Graduate Studies: Theses & Dissertations

In this thesis we examine the Navier-Stokes equations (NSE) with the continuity equation replaced by a pressure Poisson equation (PPE). Appropriate boundary conditions are developed for the PPE, which allow for a fully decoupled numerical scheme to recover the pressure. The variational form of the NSE with PPE is derived and used in the Galerkin Finite Element discretization. The Galerkin finite element method is then used to solve the NSE with PPE. Moderate accuracy is shown.

Full Newton Step Interior Point Method For Linear Complementarity Problem Over Symmetric Cones, Andrii Berdnikov

College of Graduate Studies: Theses & Dissertations

In this thesis, we present a new Feasible Interior-Point Method (IPM) for Linear Complementarity Problem (LPC) over Symmetric Cones. The advantage of this method lies in that it uses full Newton-steps, thus, avoiding the calculation of the step size at each iteration. By suitable choice of parameters we prove the global convergence of iterates which always stay in the the central path neighborhood. A global convergence of the method is proved and an upper bound for the number of iterations necessary to find ε-approximate solution of the problem is presented.

Computational Fluid Dynamics (Cfd) Modeling Of A Laboratory Scale Coal Gasifier, Kiel S. Schultheiss

College of Graduate Studies: Theses & Dissertations

Furthering gasification technology is an essential part of advancing clean coal technologies. In order to seek insight into the appropriate operations for the formation of synthetic gas (syngas) a numerical simulation was performed to predict the phenomena of coal gasification in a laboratory scale entrained-flow coal gasifier. The mesh for the model was developed with ICEM CFD software and the chemical and physical phenomena were modeled using the fluid flow solver ANSYS FLUENT. Mesh independence was verified. The model was validated with experimental data from several studies performed on a laboratory scale gasifier.

Systematic examination of the model was performed …