Applied Mathematics Commons^™

Development And Internal Validation Of An Aneurysm Rupture Probability Model Based On Patient Characteristics And Aneurysm Location, Morphology, And Hemodynamics, Felicitas J. Detmer, Bong Jae Chung, Fernando Mut, Martin Slawski, Farid Hamzei-Sichani, Christopher Putman, Carlos Jiménez, Juan R. Cebral

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Purpose: Unruptured cerebral aneurysms pose a dilemma for physicians who need to weigh the risk of a devastating subarachnoid hemorrhage against the risk of surgery or endovascular treatment and their complications when deciding on a treatment strategy. A prediction model could potentially support such treatment decisions. The aim of this study was to develop and internally validate a model for aneurysm rupture based on hemodynamic and geometric parameters, aneurysm location, and patient gender and age. Methods: Cross-sectional data from 1061 patients were used for image-based computational fluid dynamics and shape characterization of 1631 aneurysms for training an aneurysm rupture probability …

Identification Of Hostile Hemodynamics And Geometries Of Cerebral Aneurysms: A Case-Control Study, Bong Jae Chung, F. Mut, C. M. Putman, F. Hamzei-Sichani, W. Brinjikji, D. Kallmes, C. M. Jimenez, Juan R. Cebral

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

BACKGROUND AND PURPOSE: Hostile hemodynamic conditions and geometries are thought to predispose aneurysms for instability and rupture. This study compares stable, unstable, and ruptured aneurysms while controlling for location and patient characteristics. MATERIALS AND METHODS: The hemodynamics and geometries of 165 stable, 65 unstable, and 554 ruptured aneurysms were compared. Hemodynamics was modeled using image-based computational fluid dynamics. Case-control pairs were selected matching aneurysm location, patient age, and sex. Paired Wilcoxon tests were used to compare hemodynamic and geometric variables among different aneurysm groups. The pairing was repeated 100 times, and the combined P values were calculated and adjusted for …

Development Of A Statistical Model For Discrimination Of Rupture Status In Posterior Communicating Artery Aneurysms, Felicitas J. Detmer, Bong Jae Chung, Fernando Mut, Michael Pritz, Martin Slawski, Farid Hamzei-Sichani, David Kallmes, Christopher Putman, Carlos Jimenez, Juan R. Cebral

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Background: Intracranial aneurysms at the posterior communicating artery (PCOM) are known to have high rupture rates compared to other locations. We developed and internally validated a statistical model discriminating between ruptured and unruptured PCOM aneurysms based on hemodynamic and geometric parameters, angio-architectures, and patient age with the objective of its future use for aneurysm risk assessment. Methods: A total of 289 PCOM aneurysms in 272 patients modeled with image-based computational fluid dynamics (CFD) were used to construct statistical models using logistic group lasso regression. These models were evaluated with respect to discrimination power and goodness of fit using tenfold nested …

Dynamics Of Coupled Particles In A Time-Dependent, Double-Gyre Flow, Manuel Albrizzio

Theses, Dissertations and Culminating Projects

We consider a time-dependent, wind-driven, stochastic double-gyre flow, and investigate the interaction between the flow and coupled particles operating within the flow. It is known that noise can cause individual particles to escape from one gyre to another gyre. By computing the Lagrangian coherent structures (LCS) of the system, one can determine low and high probability regions of particle escape. We adjust the coupling between two particles, and study the effect on particle escape for a variety of initial conditions and noise intensities.

Understanding The Flow Structure Of Low Reynolds Number Flows, Albert Jarvis

Theses, Dissertations and Culminating Projects

Ocean flows and the mechanisms by which their contents are organized has been a longstanding area of interest in applied mathematics. In recent years, a new theory has been developed to identify the structures responsible for the organization of fluid particles within complex geophysical flows. This theory is known as the theory of Lagrangian Coherent Structures (LCS) and details which structures are responsible for the organization of the flow and how to identify them. Being able to identify these LCS in real time has far reaching implications ranging from developing strategies for search and rescue missions to identifying the best …

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 …

A Practical Guide To Big Data, Ekaterina Smirnova, Andrada Ivanescu, Jiawei Bai, Ciprian M. Crainiceanu

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Big Data is increasingly prevalent in science and data analysis. We provide a short tutorial for adapting to these changes and making the necessary adjustments to the academic culture to keep Biostatistics truly impactful in scientific research.

Analysis Of Daily Precipitation Data From Selected Sites In The United States, Sahar Ahmed

Theses, Dissertations and Culminating Projects

Global warming is a contentious topic since modern climate records only exist for the last 100 years in contrast to ice-core analysis that establishes ice ages tens of thousands of years ago. Nevertheless, patterns associated with events such as El Niño Southern Oscillation (ENSO), precipitation, tornadoes, and snowfall amounts over the last century can provide a useful and objective indicator of climate “change”. This project focuses on daily precipitation totals for the state of New Jersey over the last 100 to 150 years from nineteen meteorological recording stations and involves large data sets with a million observations. This research utilizes …

Inertial Particle Transport By Lagrangian Coherent Structures In Geophysical Flows, Alexa Aucoin

Theses, Dissertations and Culminating Projects

Lagrangian Coherent Structures (LCS) provide a skeleton for the underlying structures in geophysical flows. It is known that LCS govern the movement of fluid particles within a flow, but it is not well understood how these same LCS influence the movement of inertial particles within a fluid flow. In this thesis, we consider two geophysical flows, the double-gyre model, and a single-layer quasi-geostrophic PDE model. In particular, we use finite-time Lyapunov exponents (FTLE) to characterize the attracting and repelling LCS for these models and show how inertial particles aggregate with respect to LCS. We numerically investigate the dynamics of inertial …

Seasonal Switching Affects Bacterial-Fungal Dominance In An Ecological System, Kristin Carfora

Theses, Dissertations and Culminating Projects

We consider a model inspired by producer-herbivore-decomposer soil food webs and determine the effect of ecological parameters on the decomposer pool. In particular, we observe how seasonal changes in the stoichiometric quality of the producer coupled with the efficiency of herbivory over the calendar year can induce a shift in the composition of the decomposer pool. Decomposers have a significant effect on the movement of essential nutrients throughout an ecosystem; we further determine how this shift between a bacterially dominated decomposer pool and a fungally dominated pool affects primary production and relative distribution of biomass of the other compartments.

Borrowing Capacity, Financial Instability, And Contagion: Case Study Of The U.S. Subprime Mortgage Crisis, Youngna Choi

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

We use the case of the 2007 United States subprime mortgage crisis to investigate the impact of borrowing capacity limitations on financial instability and contagion. We divide an economy into agents that interact via flow of funds and express the financial instability level of each agent as a function of time derivatives of its wealth, cash inflows, and borrowing capacity. We show that among these factors, the borrowing capacity, which is determined by other economic constraints, has the largest impact on financial instability. It is suggested that borrowing capacity limitations could even cause contagion through feedback loop formed by flow …

The Mechanics And Thermodynamics Of Amyloid Beta Protein Aggregation In Competing Pathways, Edward Steen

Theses, Dissertations and Culminating Projects

The primary purpose of this paper is to investigate the mechanics of Aβ protein aggregation within the brain through mathematical modeling and simulation. Aggregation of Aβ is the cause of plaques within the brain of Alzheimer’s Disease sufferers. Because the pathways of aggregation from monomer to oligomer to polymer are numerous and complex, we have had to simplify our model to a limited number of species. Of great concern, too, is the process by which Aβ can form as “off-pathway” species, which is when Aβ reacts with fatty acid micelles. It is this species of Aβ, which due to its …

A Primer On Noise-Induced Transitions In Applied Dynamical Systems, Eric Forgoston, Richard O. Moore

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Noise plays a fundamental role in a wide variety of physical and biological dynamical systems. It can arise from an external forcing or due to random dynamics internal to the system. It is well established that even weak noise can result in large behavioral changes such as transitions between or escapes from quasi-stable states. These transitions can correspond to critical events such as failures or extinctions that make them essential phenomena to understand and quantify, despite the fact that their occurrence is rare. This article will provide an overview of the theory underlying the dynamics of rare events for stochastic …

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 …

Hemodynamic Characteristics Of Ruptured And Unruptured Multiple Aneurysms At Mirror And Ipsilateral Locations, Ravi Doddasomayajula, Bong Jae Chung, Fernando Mut, Carlos M. Jimenez, Farid Hamzei-Sichani, Christopher M. Putman, Juan R. Cebral

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

BACKGROUND AND PURPOSE: Different hemodynamic patterns have been associated with aneurysm rupture. The objective was to test whether hemodynamic characteristics of the ruptured aneurysm in patients with multiple aneurysms were different from those in unruptured aneurysms in the same patient.

MATERIALS AND METHODS: Twenty-four mirror and 58 ipsilateral multiple aneurysms with 1 ruptured and the others unruptured were studied. Computational fluid dynamics models were created from 3D angiographies. Case-control studies of mirror and ipsilateral aneurysms were performed with paired Wilcoxon tests.

RESULTS: In mirror pairs, the ruptured aneurysm had more oscillatory wall shear stress (P = .007) than the …

Angioarchitectures And Hemodynamic Characteristics Of Posterior Communicating Artery Aneurysms And Their Association With Rupture Status, Bong Jae Chung, Ravi Doddasomayajula, Fernando Mut, Felicitas Detmer, Michael Pritz, Farid Hamzei-Sichani, Waleed Brinjikji, David F. Kallmes, Carlos M. Jimenez, Christopher Putman, Juan Cebral

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

BACKGROUND AND PURPOSE: Intracranial aneurysms originating at the posterior communicating artery are known to have high rupture risk compared with other locations. We tested the hypothesis that different angioarchitectures (ie, branch point configuration) of posterior communicating artery aneurysms are associated with aneurysm hemodynamics, which in turn predisposes aneurysms to rupture.

MATERIALS AND METHODS: A total of 313 posterior communicating artery aneurysms (145 ruptured, 168 unruptured) were studied with image-based computational fluid dynamics. Aneurysms were classified into different angioarchitecture types depending on the location of the aneurysm with respect to parent artery bifurcation. Hemodynamic characteristics were compared between ruptured and unruptured …

Fluid Dynamics Of Watercolor Painting : Experiments And Modelling, David Edward Baron

Theses, Dissertations and Culminating Projects

In his classic study in 1908, A.M. Worthington gave a thorough account of splashes and their formation through visualization experiments. In more recent times, there has been renewed interest in this subject, and much of the underlying physics behind Worthington's experiments has now been clarified. One specific set of such recent studies, which motivates this thesis, concerns the fluid dynamics behind Jackson Pollock's drip paintings. The physical processes and the mathematical structures hidden in his works have received serious attention and have made the scientific pursuit of art a compelling area of exploration. Our current work explores the interaction of …

Understanding Angiography-Based Aneurysm Flow Fields Through Comparison With Computational Fluid Dynamics, Juan R. Cebral, F. Mut, Bong Jae Chung, L. Spelle, J. Moret, F. Van Nijnatten, D. Ruijters

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

BACKGROUND AND PURPOSE: Hemodynamics is thought to be an important factor for aneurysm progression and rupture. Our aim was to evaluate whether flow fields reconstructed from dynamic angiography data can be used to realistically represent the main flow structures in intracranial aneurysms. MATERIALS AND METHODS: DSA-based flow reconstructions, obtained during interventional treatment, were compared qualitatively with flow fields obtained from patient-specific computational fluid dynamics models and quantitatively with projections of the computational fluid dynamics fields (by computing a directional similarity of the vector fields) in 15 cerebral aneurysms. RESULTS: The average similarity between the DSA and the projected computational fluid …

Rainbow Perfect Matchings And Hamilton Cycles In The Random Geometric Graph, Deepak C. Bal, Patrick Bennett, Xavier Pérez‐Giménez, Paweł Prałat

Department of Mathematics Faculty Scholarship and Creative Works

Given a graph on n vertices and an assignment of colours to the edges, a rainbow Hamilton cycle is a cycle of length n visiting each vertex once and with pairwise different colours on the edges. Similarly (for even n) a rainbow perfect matching is a collection of independent edges with pairwise different colours. In this note we show that if we randomly colour the edges of a random geometric graph with sufficiently many colours, then a.a.s. the graph contains a rainbow perfect matching (rainbow Hamilton cycle) if and only if the minimum degree is at least 1 (respectively, …

Differences In Hemodynamics And Rupture Rate Of Aneurysms At The Bifurcation Of The Basilar And Internal Carotid Arteries, Ravi Doddasomayajula, Bong Jae Chung, Farid Hamzei-Sichani, Christopher M. Putman, Juan Cebral

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

BACKGROUND AND PURPOSE: Cerebral aneurysms in the posterior circulation are known to have a higher rupture risk than those in the anterior circulation. We sought to test the hypothesis that differences in hemodynamics can explain the difference in rupture rates.

MATERIALS AND METHODS: A total of 117 aneurysms, 63 at the tip of the basilar artery (27 ruptured, 36 unruptured, rupture rate = 43%) and 54 at the bifurcation of the internal carotid artery (11 ruptured, 43 unruptured, rupture rate = 20%) were analyzed with image-based computational fluid dynamics. Several hemodynamic variables were compared among aneurysms at each location and …

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 …

Flow Conditions In The Intracranial Aneurysm Lumen Are Associated With Inflammation And Degenerative Changes Of The Aneurysm Wall, J. Cebral, E. Ollikainen, Bong Jae Chung, F. Mut, V. Sippola, B. R. Jahromi, R. Tulamo, J. Hernesniemi, M. Niemelä, A. Robertson, J. Frösen

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

BACKGROUND AND PURPOSE: Saccular intracranial aneurysm is a common disease that may cause devastating intracranial hemorrhage. Hemodynamics, wall remodeling, and wall inflammation have been associated with saccular intracranial aneurysm rupture.Weinvestigated how saccular intracranial aneurysm hemodynamics is associated with wall remodeling and inflammation of the saccular intracranial aneurysm wall. MATERIALS AND METHODS: Tissue samples resected during a saccular intracranial aneurysm operation (11 unruptured, 9 ruptured) were studied with histology and immunohistochemistry. Patient-specific computational models of hemodynamics were created from preoperative CT angiographies. RESULTS: More stable and less complex flows were associated with thick, hyperplastic saccular intracranial aneurysm walls, while slower flows …

Regional Mapping Of Flow And Wall Characteristics Of Intracranial Aneurysms, Juan R. Cebral, Xinjie Duan, Piyusha S. Gade, Bong Jae Chung, Fernando Mut, Khaled Aziz, Anne M. Robertson

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

The evolution of intracranial aneurysms (IAs) is thought to be driven by progressive wall degradation in response to abnormal hemodynamics. Previous studies focused on the relationship between global hemodynamics and wall properties. However, hemodynamics, wall structure and mechanical properties of cerebral aneurysms can be non-uniform across the aneurysm wall. Therefore, the aim of this work is to introduce a methodology for mapping local hemodynamics to local wall structure in resected aneurysm specimens. This methodology combines image-based computational fluid dynamics, tissue resection, micro-CT imaging of resected specimens mounted on 3D-printed aneurysm models, alignment to 3D vascular models, multi-photon microscopy of the …

Computing The Optimal Path In Stochastic Dynamical Systems, Martha Bauver, Eric Forgoston, Lora Billings

Department of Mathematics Faculty Scholarship and Creative Works

In stochastic systems, one is often interested in finding the optimal path that maximizes the probability of escape from a metastable state or of switching between metastable states. Even for simple systems, it may be impossible to find an analytic form of the optimal path, and in high- dimensional systems, this is almost always the case. In this article, we formulate a constructive methodology that is used to compute the optimal path numerically. The method utilizes finite-time Lyapunov exponents, statistical selection criteria, and a Newton-based iterative minimizing scheme. The method is applied to four examples. The first example is a …

A Two-Strain Tb Model With Multiple Latent Stages, Azizeh Jabbari, Carlos Castillo-Chavez, Fereshteh Nazari, Baojun Song, Hossein Kheiri

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

A two-strain tuberculosis (TB) transmission model incorporating antibiotic-generated TB resistant strains and long and variable waiting periods within the latently infected class is introduced. The mathematical analysis is carried out when the waiting periods are modeled via parametrically friendly gamma distributions, a reasonable alternative to the use of exponential distributed waiting periods or to integral equations involving "arbitrary" distributions. The model supports a globally-asymptotically stable disease-free equilibrium when the reproduction number is less than one and an endemic equilibriums, shown to be locally asymptotically stable, or l.a.s., whenever the basic reproduction number is greater than one. Conditions for the existence …

The Importance Of Prediction Model Validation And Assessment In Obesity And Nutrition Research, Andrada Ivanescu, P. Li, B. George, A. W. Brown, S. W. Keith, D. Raju, D. B. Allison

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Deriving statistical models to predict one variable from one or more other variables, or predictive modeling, is an important activity in obesity and nutrition research. To determine the quality of the model, it is necessary to quantify and report the predictive validity of the derived models. Conducting validation of the predictive measures provides essential information to the research community about the model. Unfortunately, many articles fail to account for the nearly inevitable reduction in predictive ability that occurs when a model derived on one data set is applied to a new data set. Under some circumstances, the predictive validity can …

Splash Dynamics Of Paint On Dry, Wet, And Cooled Surfaces, David Baron, Haiyan Su, Ashuwin Vaidya

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

In his classic study in 1908, A.M. Worthington gave a thorough account of splashes and their formation through visualization experiments. In more recent times, there has been renewed interest in this subject, and much of the underlying physics behind Worthington's experiments has now been clarified. One specific set of such recent studies, which motivates this paper, concerns the fluid dynamics behind Jackson Pollock's drip paintings. The physical processes and the mathematical structures hidden in his works have received serious attention and made the scientific pursuit of art a compelling area of exploration. Our current work explores the interaction of watercolors …

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 …

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 …

Hemodynamic Analysis Of Fast And Slow Aneurysm Occlusions By Flow Diversion In Rabbits, Bong Jae Chung, Fernando Mut, Ramanathan Kadirvel, Ravi Lingineni, David F. Kallmes, Juan R. Cebral

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Purpose: To assess hemodynamic differences between aneurysms that occlude rapidly and those occluding in delayed fashion after flow diversion in rabbits. Methods: Thirty-six elastase-induced aneurysms in rabbits were treated with flow diverting devices. Aneurysm occlusion was assessed angiographically immediately before they were sacrificed at 1 (n=6), 2 (n=4), 4 (n=8) or 8 weeks (n=18) after treatment. The aneurysms were classified into a fast occlusion group if they were completely or near completely occluded at 4 weeks or earlier and a slow occlusion group if they remained incompletely occluded at 8 weeks. The immediate post-treatment flow conditions in aneurysms of each …