Visualizing Lagrangian Coherent Structures from Realistic Flow Data

136845-Thumbnail Image.png
Description
The goal of this project was to examine the separatricies that define regions of distinct flow behaviors in realistic time-dependent dynamical systems. In particular, we adapted previously available methods for computing the Finite-Time Lyapunov Exponent (FTLE) to a set of…

The goal of this project was to examine the separatricies that define regions of distinct flow behaviors in realistic time-dependent dynamical systems. In particular, we adapted previously available methods for computing the Finite-Time Lyapunov Exponent (FTLE) to a set of measured wind velocity data in order to visualize the separatricies as ridges of the FTLE field in a section of the atmosphere. This visualization required a number of alterations to the original methods, including interpolation techniques and two different adaptive refinement schemes for producing more detailed results. Overall, there were two computations performed with the wind velocity data: once along a single spherical surface, on which the separatricies could be visualized as material lines, and then along a three-dimensional section of the atmosphere, for which the separatricies were material surfaces. The resulting figures provide an image of the Antarctic polar vortex from the wind velocity data, which is consistent with other data gathered on the same date.
Date Created
2014-05
Agent

Prediction of heat transport in multiple tokamak devices with neural networks

136271-Thumbnail Image.png
Description
The OMFIT (One Modeling Framework for Integrated Tasks) modeling environment and the BRAINFUSE module have been deployed on the PPPL (Princeton Plasma Physics Laboratory) computing cluster with modifications that have rendered the application of artificial neural networks (NNs) to the…

The OMFIT (One Modeling Framework for Integrated Tasks) modeling environment and the BRAINFUSE module have been deployed on the PPPL (Princeton Plasma Physics Laboratory) computing cluster with modifications that have rendered the application of artificial neural networks (NNs) to the TRANSP databases for the JET (Joint European Torus), TFTR (Tokamak Fusion Test Reactor), and NSTX (National Spherical Torus Experiment) devices possible through their use. This development has facilitated the investigation of NNs for predicting heat transport profiles in JET, TFTR, and NSTX, and has promoted additional investigations to discover how else NNs may be of use to scientists at PPPL. In applying NNs to the aforementioned devices for predicting heat transport, the primary goal of this endeavor is to reproduce the success shown in Meneghini et al. in using NNs for heat transport prediction in DIII-D. Being able to reproduce the results from is important because this in turn would provide scientists at PPPL with a quick and efficient toolset for reliably predicting heat transport profiles much faster than any existing computational methods allow; the progress towards this goal is outlined in this report, and potential additional applications of the NN framework are presented.
Date Created
2015-05
Agent

Prey-predator-parasite: an ecosystem model with fragile persistence

155408-Thumbnail Image.png
Description
Using a simple $SI$ infection model, I uncover the

overall dynamics of the system and how they depend on the incidence

function. I consider both an epidemic and endemic perspective of the

model, but in both cases, three classes of incidence

functions are identified.

In…

Using a simple $SI$ infection model, I uncover the

overall dynamics of the system and how they depend on the incidence

function. I consider both an epidemic and endemic perspective of the

model, but in both cases, three classes of incidence

functions are identified.

In the epidemic form,

power incidences, where the infective portion $I^p$ has $p\in(0,1)$,

cause unconditional host extinction,

homogeneous incidences have host extinction for certain parameter constellations and

host survival for others, and upper density-dependent incidences

never cause host extinction. The case of non-extinction in upper

density-dependent

incidences extends to the case where a latent period is included.

Using data from experiments with rhanavirus and salamanders,

maximum likelihood estimates are applied to the data.

With these estimates,

I generate the corrected Akaike information criteria, which

reward a low likelihood and punish the use of more parameters.

This generates the Akaike weight, which is used to fit

parameters to the data, and determine which incidence functions

fit the data the best.

From an endemic perspective, I observe

that power incidences cause initial condition dependent host extinction for

some parameter constellations and global stability for others,

homogeneous incidences have host extinction for certain parameter constellations and

host survival for others, and upper density-dependent incidences

never cause host extinction.

The dynamics when the incidence function is homogeneous are deeply explored.

I expand the endemic considerations in the homogeneous case

by adding a predator into the model.

Using persistence theory, I show the conditions for the persistence of each of the

predator, prey, and parasite species. Potential dynamics of the system include parasite mediated

persistence of the predator, survival of the ecosystem at high initial predator levels and

ecosystem collapse at low initial predator levels, persistence of all three species, and much more.
Date Created
2017
Agent

Rotating split-cylinder flows

155323-Thumbnail Image.png
Description
The three-dimensional flow contained in a rapidly rotating circular

split cylinder is studied numerically solving the Navier--Stokes

equations. The cylinder is completely filled with fluid

and is split at the midplane. Three different types of boundary

conditions were imposed, leading to a variety…

The three-dimensional flow contained in a rapidly rotating circular

split cylinder is studied numerically solving the Navier--Stokes

equations. The cylinder is completely filled with fluid

and is split at the midplane. Three different types of boundary

conditions were imposed, leading to a variety of instabilities and

complex flow dynamics.

The first configuration has a strong background rotation and a small

differential rotation between the two halves. The axisymmetric flow

was first studied identifying boundary layer instabilities which

produce inertial waves under some conditions. Limit cycle states and

quasiperiodic states were found, including some period doubling

bifurcations. Then, a three-dimensional study was conducted

identifying low and high azimuthal wavenumber rotating waves due to

G’ortler and Tollmien–-Schlichting type instabilities. Over most of

the parameter space considered, quasiperiodic states were found where

both types of instabilities were present.

In the second configuration, both cylinder halves are in exact

counter-rotation, producing an O(2) symmetry in the system. The basic state flow dynamic

is dominated by the shear layer created

in the midplane. By changing the speed rotation and the aspect ratio

of the cylinder, the flow loses symmetries in a variety of ways

creating static waves, rotating waves, direction reversing waves and

slow-fast pulsing waves. The bifurcations, including infinite-period

bifurcations, were characterized and the flow dynamics was elucidated.

Additionally, preliminary experimental results for this case are

presented.

In the third set up, with oscillatory boundary conditions, inertial

wave beams were forced imposing a range of frequencies. These beams

emanate from the corner of the cylinder and from the split at the

midplane, leading to destructive/constructive interactions which

produce peaks in vorticity for some specific frequencies. These

frequencies are shown to be associated with the resonant Kelvin

modes. Furthermore, a study of the influence of imposing a phase

difference between the oscillations of the two halves of the cylinder

led to the interesting result that different Kelvin

modes can be excited depending on the phase difference.
Date Created
2017
Agent

Extraction of coherent structures using direct numerical simulation in 3D turbulent flows and its effects on chemotaxis

155243-Thumbnail Image.png
Description
A numerical study of chemotaxis in 3D turbulence is presented here. Direct Numerical

Simulation were used to calculate the nutrient uptake for both motile and non-motile bacterial

species and by applying the dynamical systems theory the effect of flow topology…

A numerical study of chemotaxis in 3D turbulence is presented here. Direct Numerical

Simulation were used to calculate the nutrient uptake for both motile and non-motile bacterial

species and by applying the dynamical systems theory the effect of flow topology on the

variability of chemotaxis is analyzed. It is done by injecting a highly localized patch of nutrient

in the turbulent flow, and analyzing the evolution of reaction associated with the observed

high and low stretching regions. The Gaussian nutrient patch is released at different locations

and the corresponding nutrient uptake is obtained. The variable stretching characteristics of

the flow is depicted by Lagrangian Coherent Structures and the roles they play in affecting the

uptake are analyzed. The Lagrangian Coherent Structures are quantified by the Finite Time

Lyapunov Exponents which is a measure of the average stretching experienced by the flow in

finite time. It is found that in high stretching regions, the motile bacteria are attracted to the

nutrient patch very quickly, but also dispersed quickly; whereas in low stretching regions the

bacteria respond slower towards the nutrient patch. However the total uptake is intricately

determined by stretching history. These reaction characteristics are reflected in the several

realizations of simulations. This helps in understanding turbulence intensity and how it affects

the uptake of the nutrient.
Date Created
2017
Agent

The domain dependence of chemotaxis in a two-dimensional turbulent flow

153936-Thumbnail Image.png
Description
Presented is a study on the chemotaxis reaction process and its relation with flow topology. The effect of coherent structures in turbulent flows is characterized by studying nutrient uptake and the advantage that is received from motile bacteria over other…

Presented is a study on the chemotaxis reaction process and its relation with flow topology. The effect of coherent structures in turbulent flows is characterized by studying nutrient uptake and the advantage that is received from motile bacteria over other non-motile bacteria. Variability is found to be dependent on the initial location of scalar impurity and can be tied to Lagrangian coherent structures through recent advances in the identification of finite-time transport barriers. Advantage is relatively small for initial nutrient found within high stretching regions of the flow, and nutrient within elliptic structures provide the greatest advantage for motile species. How the flow field and the relevant flow topology lead to such a relation is analyzed.
Date Created
2015
Agent

Lagrangian Coherent Structure Analysis of Terminal Winds: Three-Dimensionality, Intramodel Variations, and Flight Analyses

129265-Thumbnail Image.png
Description

We present a study of three-dimensional Lagrangian coherent structures (LCS) near the Hong Kong International Airport and relate to previous developments of two-dimensional (2D) LCS analyses. The LCS are contrasted among three independent models and against 2D coherent Doppler light…

We present a study of three-dimensional Lagrangian coherent structures (LCS) near the Hong Kong International Airport and relate to previous developments of two-dimensional (2D) LCS analyses. The LCS are contrasted among three independent models and against 2D coherent Doppler light detection and ranging (LIDAR) data. Addition of the velocity information perpendicular to the LIDAR scanning cone helps solidify flow structures inferred from previous studies; contrast among models reveals the intramodel variability; and comparison with flight data evaluates the performance among models in terms of Lagrangian analyses. We find that, while the three models and the LIDAR do recover similar features of the windshear experienced by a landing aircraft (along the landing trajectory), their Lagrangian signatures over the entire domain are quite different—a portion of each numerical model captures certain features resembling those LCS extracted from independent 2D LIDAR analyses based on observations.

Date Created
2014-11-30
Agent

Application of Short-Range LIDAR in Early Alerting for Low-Level Windshear and Turbulence at Hong Kong International Airport

129480-Thumbnail Image.png
Description

Hong Kong Observatory currently uses a series of meteorological instruments, including long-range LIDAR (light detection and ranging) systems, to provide alerting services of low-level windshear and turbulence for Hong Kong International Airport. For some events that are smaller in spatial…

Hong Kong Observatory currently uses a series of meteorological instruments, including long-range LIDAR (light detection and ranging) systems, to provide alerting services of low-level windshear and turbulence for Hong Kong International Airport. For some events that are smaller in spatial dimensions and are rapidly changing, such as low altitude windshear and turbulence associated with buildings or man-made structures, it would be necessary to involve meteorological instruments that offer greater spatial resolution. Therefore, the Observatory has set up a short-range LIDAR on the roof of the AsiaWorld-Expo during the summers over the past several years, conducting field research on the feasibility of strengthening early alerting for windshear and turbulence over the north runway’s eastern arrival runway (Runway 25RA) and developing an automated early alerting algorithm. This paper takes the pilot reports for Runway 25RA during the 2013 field research as verification samples, using different thresholds for radial wind velocity spatial and temporal changes detected by the short-range LIDAR to calculate the relative operating characteristic (ROC) curve, and analyzes its early alerting performance.

Date Created
2013-11-30
Agent

Mathematical and statistical insights in evaluating state dependent effectiveness of HIV prevention interventions

153290-Thumbnail Image.png
Description
Pre-Exposure Prophylaxis (PrEP) is any medical or public health procedure used before exposure to the disease causing agent, its purpose is to prevent, rather than treat or cure a disease. Most commonly, PrEP refers to an experimental HIV-prevention strategy that…

Pre-Exposure Prophylaxis (PrEP) is any medical or public health procedure used before exposure to the disease causing agent, its purpose is to prevent, rather than treat or cure a disease. Most commonly, PrEP refers to an experimental HIV-prevention strategy that would use antiretrovirals to protect HIV-negative people from HIV infection. A deterministic mathematical model of HIV transmission is developed to evaluate the public-health impact of oral PrEP interventions, and to compare PrEP effectiveness with respect to different evaluation methods. The effects of demographic, behavioral, and epidemic parameters on the PrEP impact are studied in a multivariate sensitivity analysis. Most of the published models on HIV intervention impact assume that the number of individuals joining the sexually active population per year is constant or proportional to the total population. In the second part of this study, three models are presented and analyzed to study the PrEP intervention, with constant, linear, and logistic recruitment rates. How different demographic assumptions can affect the evaluation of PrEP is studied. When provided with data, often least square fitting or similar approaches can be used to determine a single set of approximated parameter values that make the model fit the data best. However, least square fitting only provides point estimates and does not provide information on how strongly the data supports these particular estimates. Therefore, in the third part of this study, Bayesian parameter estimation is applied on fitting ODE model to the related HIV data. Starting with a set of prior distributions for the parameters as initial guess, Bayes' formula can be applied to obtain a set of posterior distributions for the parameters which makes the model fit the observed data best. Evaluating the posterior distribution often requires the integration of high-dimensional functions, which is usually difficult to calculate numerically. Therefore, the Markov chain Monte Carlo (MCMC) method is used to approximate the posterior distribution.
Date Created
2014
Agent

Stoichiometric producer-grazer models, incorporating the effects of excess food-nutrient content on grazer dynamics

152845-Thumbnail Image.png
Description
There has been important progress in understanding ecological dynamics through the development of the theory of ecological stoichiometry. This fast growing theory provides new constraints and mechanisms that can be formulated into mathematical models. Stoichiometric models incorporate the effects of…

There has been important progress in understanding ecological dynamics through the development of the theory of ecological stoichiometry. This fast growing theory provides new constraints and mechanisms that can be formulated into mathematical models. Stoichiometric models incorporate the effects of both food quantity and food quality into a single framework that produce rich dynamics. While the effects of nutrient deficiency on consumer growth are well understood, recent discoveries in ecological stoichiometry suggest that consumer dynamics are not only affected by insufficient food nutrient content (low phosphorus (P): carbon (C) ratio) but also by excess food nutrient content (high P:C). This phenomenon, known as the stoichiometric knife edge, in which animal growth is reduced not only by food with low P content but also by food with high P content, needs to be incorporated into mathematical models. Here we present Lotka-Volterra type models to investigate the growth response of Daphnia to algae of varying P:C ratios. Using a nonsmooth system of two ordinary differential equations (ODEs), we formulate the first model to incorporate the phenomenon of the stoichiometric knife edge. We then extend this stoichiometric model by mechanistically deriving and tracking free P in the environment. This resulting full knife edge model is a nonsmooth system of three ODEs. Bifurcation analysis and numerical simulations of the full model, that explicitly tracks phosphorus, leads to quantitatively different predictions than previous models that neglect to track free nutrients. The full model shows that the grazer population is sensitive to excess nutrient concentrations as a dynamical free nutrient pool induces extreme grazer population density changes. These modeling efforts provide insight on the effects of excess nutrient content on grazer dynamics and deepen our understanding of the effects of stoichiometry on the mechanisms governing population dynamics and the interactions between trophic levels.
Date Created
2014
Agent
Subscribe to Tang, Wenbo