We review and test twelve different approaches to the detection of finite-time coherent material structures in two-dimensional, temporally aperiodic flows. We consider both mathematical methods and diagnostic scalar fields, comparing their performance on three benchmark examples: the quasiperiodically forced Bickley jet, a two-dimensional turbulence simulation, and an observational wind velocity field from Jupiter's atmosphere. A close inspection of the results reveals that the various methods often produce very different predictions for coherent structures, once they are evaluated beyond heuristic visual assessment. As we find by passive advection of the coherent set candidates, false positives and negatives can be produced even by some of the mathematically justified methods due to the ineffectiveness of their underlying coherence principles in certain flow configurations. We summarize the inferred strengths and weaknesses of each method, and make general recommendations for minimal self-consistency requirements that any Lagrangian coherence detection technique should satisfy.
Extreme events, such as rogue waves, earthquakes, and stock market crashes, occur spontaneously in many dynamical systems. Because of their usually adverse consequences, quantification, prediction, and mitigation of extreme events are highly desirable. Here, we review several aspects of extreme events in phenomena described by high-dimensional, chaotic dynamical systems. We especially focus on two pressing aspects of the problem: (i) mechanisms underlying the formation of extreme events and (ii) real-time prediction of extreme events. For each aspect, we explore methods relying on models, data, or both. We discuss the strengths and limitations of each approach as well as possible future research directions.
We derive precursors of extreme dissipation events in a turbulent channel flow. Using a recently developed method that combines dynamics and statistics for the underlying attractor, we extract a characteristic state that precedes laminarization events that subsequently lead to extreme dissipation episodes. Our approach utilizes coarse statistical information for the turbulent attractor, in the form of second-order statistics, to identify high-likelihood regions in the state space. We then search within this high-probability manifold for the state that leads to the most finite-time growth of the flow kinetic energy. This state has both high probability of occurrence and leads to extreme values of dissipation. We use the alignment between a given turbulent state and this critical state as a precursor for extreme events and demonstrate its favorable properties for prediction of extreme dissipation events. Finally, we analyze the physical relevance of the derived precursor and show its robust character for different Reynolds numbers. Overall, we find that our choice of precursor works well at the Reynolds number it is computed at and at higher Reynolds number flows with similar extreme events.
11 January 2023 in Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
Compartmental models are an important quantitative tool in epidemiology, enabling us to forecast the course of a communicable disease. However, the model parameters, such as the infectivity rate of the disease, are riddled with uncertainties, which has motivated the development and use of stochastic compartmental models. Here, we first show that a common stochastic model, which treats the uncertainties as white noise, is fundamentally flawed since it erroneously implies that greater parameter uncertainties will lead to the eradication of the disease. Then, we present a principled modelling of the uncertainties based on reasonable assumptions on the contacts of each individual. Using the central limit theorem and Doob’s theorem on Gaussian Markov processes, we prove that the correlated Ornstein–Uhlenbeck (OU) process is the appropriate tool for modelling uncertainties in the infectivity rate. We demonstrate our results using a compartmental model of the COVID-19 pandemic and the available US data from the Johns Hopkins University COVID-19 database. In particular, we show that the white noise stochastic model systematically underestimates the severity of the Omicron variant of COVID-19, whereas the OU model correctly forecasts the course of this variant. Moreover, using an SIS model of sexually transmitted disease, we derive an exact closed-form solution for the final distribution of infected individuals. This analytical result shows that the white noise model underestimates the severity of the pandemic because of unrealistic noise-induced transitions. Our results strongly support the need for temporal correlations in modelling of uncertainties in compartmental models of infectious disease.
24 January 2022 in SIAM Journal on Scientific Computing
Reduced-order models of time-dependent partial differential equations (PDEs) where the solution is assumed as a linear combination of prescribed modes are rooted in a well-developed theory. However...