Applied Mathematics Commons^™

Wall Mechanical Properties And Hemodynamics Of Unruptured Intracranial Aneurysms, J. R. Cebral, X. Duan, Bong Jae Chung, C. Putman, Khaled Aziz, A. M. Robertson

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

BACKGROUND AND PURPOSE: Aneurysm progression and rupture is thought to be governed by progressive degradation and weakening of the wall in response to abnormal hemodynamics. Our goal was to investigate the relationship between the intra-aneurysmal hemodynamic conditions and wall mechanical properties in human aneurysms. MATERIALS AND METHODS: A total of 8 unruptured aneurysms were analyzed. Computational fluid dynamics models were constructed from preoperative 3D rotational angiography images. The aneurysms were clipped, and the domes were resected and mechanically tested to failure with a uniaxial testing system under multiphoton microscopy. Linear regression analysis was performed to explore possible correlations between hemodynamic …

A Framework For Inferring Unobserved Multistrain Epidemic Subpopulations Using Synchronization Dynamics, Eric Forgoston, Leah B. Shaw, Ira B. Schwartz

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

A new method is proposed to infer unobserved epidemic subpopulations by exploiting the synchronization properties of multistrain epidemic models. A model for dengue fever is driven by simulated data from secondary infective populations. Primary infective populations in the driven system synchronize to the correct values from the driver system. Most hospital cases of dengue are secondary infections, so this method provides a way to deduce unobserved primary infection levels. We derive center manifold equations that relate the driven system to the driver system and thus motivate the use of synchronization to predict unobserved primary infectives. Synchronization stability between primary and …

Penalized Function-On-Function Regression, Andrada Ivanescu, Ana Maria Staicu, Fabian Scheipl, Sonja Greven

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

A general framework for smooth regression of a functional response on one or multiple functional predictors is proposed. Using the mixed model representation of penalized regression expands the scope of function-on-function regression to many realistic scenarios. In particular, the approach can accommodate a densely or sparsely sampled functional response as well as multiple functional predictors that are observed on the same or different domains than the functional response, on a dense or sparse grid, and with or without noise. It also allows for seamless integration of continuous or categorical covariates and provides approximate confidence intervals as a by-product of the …

On The 3-Dimensional Fluid-Structure Interaction Of Flexible Fibers In A Flow, Ryan Howard Allaire

Theses, Dissertations and Culminating Projects

We discuss the equilibrium configurations 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. Experimental results are presented with both two-dimensional and three-dimensional numerical simulations used to model this problem. We present the effects of flow speed and initial configuration angle between the fiber and the direction of the flow. Investigations reveal that both the orientation of the fiber and the fiber length have a significant impact on the deformation of the fiber as well as on the forces it experiences. Specifically, we …

Predicting Successful Long-Term Weight Loss From Short-Term Weight-Loss Outcomes: New Insights From A Dynamic Energy Balance Model (The Pounds Lost Study), Diana Thomas, W Andrada Ivanescu, Corby K. Martin, Steven B. Heymsfield, Kaitlyn Marshall, Victoria E. Bodrato, Donald Williamson, Stephen Anton, Frank M. Sacks, Donna Ryan, George A. Bray

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Background: Currently, early weight-loss predictions of long-term weight-loss success rely on fixed percent-weight-loss thresholds.

Objective: The objective was to develop thresholds during the first 3 mo of intervention that include the influence of age, sex, baseline weight, percent weight loss, and deviations from expected weight to predict whether a participant is likely to lose 5% or more body weight by year 1.

Design: Data consisting of month 1, 2, 3, and 12 treatment weights were obtained from the 2-y Preventing Obesity Using Novel Dietary Strategies (POUNDS Lost) intervention. Logistic regression models that included covariates of age, height, sex, baseline weight, …

Modeling Contagion In The Eurozone Crisis Via Dynamical Systems, Giuseppe Castellacci, Youngna Choi

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

We recently (Castellacci and Choi, 2013) formulated a theoretical framework for the modeling of financial instability contagion using the theories of dynamical systems. Here, our main goal is to model the Eurozone financial crisis within that framework. The underlying system comprises many economic agents that belong to several subsystems. In each instantiation of this framework, the hierarchy and nesting of the subsystems is dictated by the nature of the problem at hand. We describe in great detail how a suitable model can be set up for the Eurozone crisis. The dynamical system is defined by the evolution of the wealths …

The Coupled Within- And Between-Host Dynamics In The Evolution Of Hiv/Aids In China, Jie Lou, Hongna Zhou, Dong Liang, Zhen Jin, Baojun Song

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

In this work, we develop and analyze mathematical models for the coupled within-host and between-host dynamics caricaturing the evolution of HIV/AIDS. The host population is divided into susceptible, the infected without receiving treatment and the infected receiving ART treatment in accordance with China’s Four-Free-One-Care Policy. The within-host model is a typical ODE model adopted from literatures. The between-host model incorporates age-since-infection described by a system of integrodifferential equations. The two models are coupled via the viral load and number of CD4+ T cells of within the hosts. For the between-host model with an arbitrarily selected HIV infected individual, we focus …

Magnetoviscous Effects Of Magnetized Particle Threads In Magnetized Ferrofluid, Alexander Francis Cali

Theses, Dissertations and Culminating Projects

The magnetoviscous effect of applied fields on ferrofluids has been utilized in many applications in which the ferrofluid must remain in a fixed position while this effect on ferrofluids in motion has yet to be rigorously explored. In light of potential biomedical applications such as drug targeting, experiments were conducted to probe the rheology of ferrofluids on the micrometer scale. A non-conducting glass sphere of diameter 550 μm is dropped into a cylindrical container of magnetized ferrofluid of inner diameter 5.2 mm. This was repeated for two applied field strengths (980 gauss and 480 gauss) and over multiple angles with …

Experimental Validation Of Robotic Manifold Tracking In Gyre-Like Flows, Matthew Michini, M. Ani Hsieh, Eric Forgoston, Ira B. Schwartz

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

In this paper, we present a first attempt toward experimental validation of a multi-robot strategy for tracking manifolds and Lagrangian coherent structures (LCS) in flows. LCS exist in natural fluid flows at various scales, and they are time-varying extensions of stable and unstable manifolds of time invariant dynamical systems. In this work, we present the first steps toward experimentally validating our previously proposed real-time manifold and LCS tracking strategy that relies solely on local measurements. Although we have validated the strategy in simulations using analytical flow models, experimental flow data, and actual ocean data, the strategy has never been implemented …

Transportation Of Perishable And Refrigerated Foods In Mylar Foil Bags And Insulated Containers: A Time-Temperature Study, Yanyan Li, John P. Schrade, Haiyan Su, John Specchio

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Data are lacking on the temperature changes of food during transport without the use of refrigerated trucks. The purpose of this study was to evaluate the ability of several insulated and noninsulated containers with or without frozen gel packs to keep perishable and refrigerated foods within the temperature safe zone in relationship to duration of transport. The study was designed to duplicate the practices exhibited by customers purchasing perishable food products from a cash-and-carry business. Approximately 40 perishable food items were evaluated. Four types of containers were tested: a mylar foil bag, a commercial insulated bag, a generic insulated bag, …

Different Types Of Backward Bifurcations Due To Density-Dependent Treatments, Baojun Song, Wen Du, Jie Lou

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

A set of deterministic SIS models with density-dependent treatments are studied to understand the disease dynamics when different treatment strategies are applied. Qualitative analyses are carried out in terms of general treatment functions. It has become customary that a backward bifurcation leads to bistable dynamics. However, this study finds that finds that bistability may not be an option at all; the disease-free equilibrium could be globally stable when there is a backward bifurcation. Furthermore, when a backward bifurcation occurs, the fashion of bistability could be the coexistence of either dual stable equilibria or the disease-free equilibrium and a stable limit …

Predicting Unobserved Exposures From Seasonal Epidemic Data, Eric Forgoston, Ira B. Schwartz

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

We consider a stochastic Susceptible-Exposed-Infected-Recovered (SEIR) epidemiological model with a contact rate that fluctuates seasonally. Through the use of a nonlinear, stochastic projection, we are able to analytically determine the lower dimensional manifold on which the deterministic and stochastic dynamics correctly interact. Our method produces a low dimensional stochastic model that captures the same timing of disease outbreak and the same amplitude and phase of recurrent behavior seen in the high dimensional model. Given seasonal epidemic data consisting of the number of infectious individuals, our method enables a data-based model prediction of the number of unobserved exposed individuals over very …

Brief Report: Parental Child-Directed Speech As A Predictor Of Receptive Language In Children With Autism Symptomatology, Twyla Y. Perryman, Alice S. Carter, Daniel S. Messinger, Wendy L. Stone, Andrada Ivanescu, Paul J. Yoder

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Facilitative linguistic input directly connected to children's interest and focus of attention has become a recommended component of interventions for young children with autism spectrum disorder (ASD). This longitudinal correlational study used two assessment time points and examined the association between parental undemanding topic-continuing talk related to the child's attentional focus (i.e.; follow-in comments) and later receptive language for 37 parent-child dyads with their young (mean = 21 months, range 15-24 months) children with autism symptomology. The frequency of parental follow-in comments positively predicted later receptive language after considering children's joint attention skills and previous receptive language abilities.

Mathematical Modelling And Control Of Echinococcus In Qinghai Province, China, Liumei Wu, Baojun Song, Wen Du, Jie Lou

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

In this paper, two mathematical models, the baseline model and the intervention model, are proposed to study the transmission dynamics of echinococcus. A global forward bifurcation completely characterizes the dynamical behavior of the baseline model. That is, when the basic reproductive number is less than one, the disease-free equilibrium is asymptotically globally stable; when the number is greater than one, the endemic equilibrium is asymptotically globally stable. For the intervention model, however, the basic reproduction number alone is not enough to describe the dynamics, particularly for the case where the basic reproductive number is less then one. The emergence of …

An Experimental Testbed For Multi-Robot Tracking Of Manifolds And Coherent Structures In Flows, Matthew Michini, Kenneth Mallory, Dennis Larkin, M. Ani Hsieh, Eric Forgoston, Philip A. Yecko

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

In this paper, we describe the development of an experimental testbed capable of producing controllable ocean-like flows in a laboratory setting. The objective is to develop a testbed to evaluate multi-robot strategies for tracking manifolds and Lagrangian coherent structures (LCS) in the ocean. Recent theoretical results have shown that LCS coincide with minimum energy and minimum time optimal paths for autonomous vehicles in the ocean. Furthermore, knowledge of these structures enables the prediction and estimation of the underlying fluid dynamics. The testbed is a scaled flow tank capable of generating complex and controlled quasi-2D flow fields that exhibit wind-driven doublegyre …

Thermalization And Initial State-Recurrence In Discrete Kdv-Like Lattices, Garrett Taylor Nieddu

Theses, Dissertations and Culminating Projects

Three discretizations of the Korteweg de-Vries equation are studied; convergence rate, initial state-recurrence, and the energy distribution of the three schemes are all considered. For each discrete scheme over 300 lattices with varying grid sizes were investigated, and the solutions were compared with other lattices from the same scheme, as well as solutions from the other two. It is found that the two schemes that are least accurate display the best recurrence at intermediate grid sizes, away from convergence. This is a notable result because the best recurrence is expected to be found in the most accurate, and converged lattices. …

Modeling The Curvature Of A Ferrofluid Interface Using A Height Function Method, Holly Timme

Theses, Dissertations and Culminating Projects

The behavior of an interface embedded in a fluid is central to a wide range of biological, chemical, environmental and physical problems and engineering processes. Modeling the evolution of a fluid interface is thus a critical and important problem. In many instances, including two-phase (e.g. liquid-gas) flows, the interface is an internal boundary within a PDE model. A model of the interface properties and its evolution is then typically performed by numerical computation, within the framework of the PDE solution method, such as finite differences (FD). Volume of Fluid (VOF) is a simple FD based method which exhibits excellent volume …

Coherent Pattern Prediction In Swarms Of Delay-Coupled Agents, Luis Mier-Y-Teran-Romero, Eric Forgoston, Ira B. Schwartz

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

We consider a general swarm model of self-propelling agents interacting through a pairwise potential in the presence of noise and communication time delay. Previous work has shown that a communication time delay in the swarm induces a pattern bifurcation that depends on the size of the coupling amplitude. We extend these results by completely unfolding the bifurcation structure of the mean field approximation. Our analysis reveals a direct correspondence between the different dynamical behaviors found in different regions of the coupling-time delay plane with the different classes of simulated coherent swarm patterns. We derive the spatiotemporal scales of the swarm …

The Persistence Of Infectious Diseases In Metapopulations, Jonathan Calvin Hayes

Theses, Dissertations and Culminating Projects

Mathematical models provide a great deal of information about the dynamics of disease spread. In this paper, we use stochastic simulation to investigate spontaneous disease extinction and réintroduction in a SIR model. We begin by investigating path to extinction and time to extinction in single population models, and then expand to a multipopulation model linked with linear migration. We have found that in a single population model, it is more effective to use random pulse vaccinations less per year at a higher removal rate. We have expanded this result by developing a vaccination strategy giving one large, well timed pulse …

Noise, Bifurcations, And Modeling Of Interacting Particle Systems, Luis Mier-Y-Teran-Romero, Eric Forgoston, Ira B. Schwartz

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

We consider the stochastic patterns of a system of communicating, or coupled, self-propelled particles in the presence of noise and communication time delay. For sufficiently large environmental noise, there exists a transition between a translating state and a rotating state with stationary center of mass. Time delayed communication creates a bifurcation pattern dependent on the coupling amplitude between particles. Using a mean field model in the large number limit, we show how the complete bifurcation unfolds in the presence of communication delay and coupling amplitude. Relative to the center of mass, the patterns can then be described as transitions between …

Complex Dynamics Of Discrete Seis Models With Simple Demography, Hui Cao, Yicang Zhou, Baojun Song

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

We investigate bifurcations and dynamical behaviors of discrete SEIS models with exogenous reinfections and a variety of treatment strategies. Bifurcations identified from the models include period doubling, backward, forward-backward, and multiple backward bifurcations. Multiple attractors, such as bistability and tristability, are observed. We also estimate the ultimate boundary of the infected regardless of initial status. Our rigorously mathematical analysis together with numerical simulations show that epidemiological factors alone can generate complex dynamics, though demographic factors only support simple equilibrium dynamics. Our model analysis supports and urges to treat a fixed percentage of exposed individuals.

Set-Based Corral Control In Stochastic Dynamical Systems: Making Almost Invariant Sets More Invariant, Eric Forgoston, Lora Billings, Philip Yecko, Ira Schwartz

Department of Mathematics Faculty Scholarship and Creative Works

We consider the problem of stochastic prediction and control in a time-dependent stochastic environment, such as the ocean, where escape from an almost invariant region occurs due to random fluctuations. We determine high-probability control-actuation sets by computing regions of uncertainty, almost invariant sets, and Lagrangian coherent structures. The combination of geometric and probabilistic methods allows us to design regions of control, which provide an increase in loitering time while minimizing the amount of control actuation. We show how the loitering time in almost invariant sets scales exponentially with respect to the control actuation, causing an exponential increase in loitering times …

Maximal Sensitive Dependence And The Optimal Path To Epidemic Extinction, Eric Forgoston, Simone Bianco, Leah B. Shaw, Ira B. Schwartz

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Extinction of an epidemic or a species is a rare event that occurs due to a large, rare stochastic fluctuation. Although the extinction process is dynamically unstable, it follows an optimal path that maximizes the probability of extinction. We show that the optimal path is also directly related to the finite-time Lyapunov exponents of the underlying dynamical system in that the optimal path displays maximum sensitivity to initial conditions. We consider several stochastic epidemic models, and examine the extinction process in a dynamical systems framework. Using the dynamics of the finite-time Lyapunov exponents as a constructive tool, we demonstrate that …

Epidemic Spread Of Influenza Viruses: The Impact Of Transient Populations On Disease Dynamics, Karen R. Ríos-Soto, Baojun Song, Carlos Castillo-Chavez

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

The recent H1N1 ("swine u") pandemic and recent H5N1 ("avian u") outbreaks have brought increased attention to the study of the role of animal populations as reservoirs for pathogens that could invade human populations. It is believed that pigs acquired u strains from birds and humans, acting as a mixing vessel in generating new influenza viruses. Assessing the role of animal reservoirs, particularly reservoirs involving highly mobile populations (like migratory birds), on disease dispersal and persistence is of interests to a wide range of researchers including public health experts and evolutionary biologists. This paper studies the interactions between transient and …

An Exploration Of Modeling Techniques For The Study Of The Dynamics Of E-Mail Viruses, Karin Weule

Theses, Dissertations and Culminating Projects

We analyze real data sets from two e-mail viruses, the Magistr.b and the Sircam.a to explore how we can use mathematical models to predict the behavior described by the data. Analysis of the data is conducted primarily with computer programming in MatLab. We focus mainly on the use of two continuous models commonly used in the study of biological diseases, the SIS and the SIR models. A discrete modeling approach using agent-based simulations is also explored and revealed to be potentially useful in developing a compartmentalized model that incorporates both SIS and SIR model behavior. The theory behind the continuous …

Switching Exponent Scaling Near Bifurcation Points For Non-Gaussian Noise, Lora Billings, Ira Schwartz, Mary Mccrary, A. Korotkrov, M. Dykman

Department of Mathematics Faculty Scholarship and Creative Works

We study noise-induced switching of a system close to bifurcation parameter values where the number of stable states changes. For non-Gaussian noise, the switching exponent, which gives the logarithm of the switching rate, displays a non-power-law dependence on the distance to the bifurcation point. This dependence is found for Poisson noise. Even weak additional Gaussian noise dominates switching sufficiently close to the bifurcation point, leading to a crossover in the behavior of the switching exponent.

An Empirical Likelihood-Based Method For Comparison Of Treatment Effects-Test Of Equality Of Coefficients In Linear Models, Haiyan Su, Hua Liang

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

To compare two treatment effects, which can be described as the difference of the parameters in two linear models, we propose an empirical likelihood-based method to make inference for the difference. Our method is free of the assumptions of normally distributed and homogeneous errors, and equal sample sizes. The empirical likelihood ratio for the difference of the parameters of interest is shown to be asymptotically chi-squared. Simulation experiments illustrate that our method outperforms the published ones. Our method is used to analyze a data set from a drug study.

A Time Series Analysis Of The New Jersey Meadowlands Weather And Air Quality Data, Steven Spero

Theses, Dissertations and Culminating Projects

This research applies time series methods to determine relationships among a set of weather variables which are continually monitored in the Hackensack Meadowlands region of northern New Jersey. Weather data includes chemical and atmospheric factors. Chemical factors are Nitrogen Oxide, atmospheric Ozone, Carbon Monoxide, and Carbon Dioxide. Weather factors are wind speed, barometric pressure, air temperature, humidity, and solar radiation. Additionally, traffic density and time of week are brought in as categorical factors. This research attempts to (a) introduce the reader to various time series methodologies, (b) find a significant and efficient model for forecasting Nitrogen Oxide levels, and (c) …

Association Of Workplace Chronic And Acute Stressors With Employee Weight Status: Data From Worksites In Turmoil, Isabel Diana Fernandez, Haiyan Su, Paul C. Winters, Hua Liang

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Objectives: To examine the independent and joint effects of psychosocial chronic and acute stressors with weight status and to report the intraclass correlation coefficient for body mass index (BMI). Methods: Baseline data on 2782 employees from a group-randomized weight gain prevention intervention were examined to investigate the effect of high job strain and job insecurity on BMI and on the odds of overweight/obesity including potential confounders and mediating variables. Data were analyzed using mixed models. Results: The mediating variables removed the effect of high job strain on weight (β = 0.68, P = 0.07; odds ratios = 1.34, confidence interval …

Computational Fluid Dynamics Of Aggregating Red Blood Cells In Postcapillary Venules, Bong Jae Chung, Sangho Kim, Paul C. Johnson, Aleksander S. Popel

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Aggregate formation of red blood cells (RBCs) in a postcapillary venular bifurcation is investigated with three-dimensional computer simulations using the Chimera grid method. Interaction energy between the RBCs is modelled by a depletion interaction theory; RBCs are modelled as rigid oblate ellipsoids. The cell-cell interactions of RBCs are strongly dependent on vessel geometry and shear rates. The experimental data on vessel geometry, pseudoshear rates, and Dextran concentration obtained in our previous in vivo RBC aggregation study in postcapillary venules of the rat spinotrapezius muscle were used to simulate RBC aggregation. The computational results were compared to the experimental results from …