BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Bartosz Protas

Bartosz Protas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynmath.NAmath.OCphysics.comp-phmath.APmath.DSmath-phmath.CVmath.MPmath.PRNumerical Analysisphysics.chem-ph

Catalog footprint

What is connected

21works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

Data-Driven Optimal Closures for Mean-Cluster Models: Beyond the Classical Pair Approximation

This study concerns the mean-clustering approach to modelling the evolution of lattice dynamics. Instead of tracking the state of individual lattice sites, this approach describes the time evolution of the concentrations of different cluster types. It leads to an infinite hierarchy of ordinary differential equations which must be closed by truncation using a so-called closure condition. This condition approximates the concentrations of higher-order clusters in terms of the concentrations of lower-order ones. The pair approximation is the most common form of closure. Here, we consider its generalization, termed the "optimal approximation", which we calibrate using a robust data-driven strategy. To fix attention, we focus on a recently proposed structured lattice model for a nickel-based oxide, similar to that used as cathode material in modern commercial Li-ion batteries. The form of the obtained optimal approximation allows us to deduce a simple sparse closure model. In addition to being more accurate than the classical pair approximation, this ``sparse approximation'' is also physically interpretable which allows us to a posteriori refine the hypotheses underlying construction of this class of closure models. Moreover, the mean-cluster model closed with this sparse approximation is linear and hence analytically solvable such that its parametrization is straightforward. On the other hand, parametrization of the mean-cluster model closed with the pair approximation is shown to lead to an ill-posed inverse problem.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimal Eddy Viscosity in Closure Models for 2D Turbulent Flows

We consider the question of fundamental limitations on the performance of eddy-viscosity closure models for turbulent flows, focusing on the Leith model for 2D {Large-Eddy Simulation}. Optimal eddy viscosities depending on the magnitude of the vorticity gradient are determined subject to minimum assumptions by solving PDE-constrained optimization problems defined such that the corresponding optimal Large-Eddy Simulation best matches the filtered Direct Numerical Simulation. First, we consider pointwise match in the physical space and the main finding is that with a fixed cutoff wavenumber $k_c$, the performance of the Large-Eddy Simulation systematically improves as the regularization in the solution of the optimization problem is reduced and this is achieved with the optimal eddy viscosities exhibiting increasingly irregular behavior with rapid oscillations. Since the optimal eddy viscosities do not converge to a well-defined limit as the regularization vanishes, we conclude that in this case the problem of finding an optimal eddy viscosity does not in fact have a solution and is thus ill-posed. We argue that this observation is consistent with the physical intuition concerning closure problems. The second problem we consider involves matching time-averaged vorticity spectra over small wavenumbers. It is shown to be better behaved and to produce physically reasonable optimal eddy viscosities. We conclude that while better behaved and hence practically more useful eddy viscosities can be obtained with stronger regularization or by matching quantities defined in a statistical sense, the corresponding Large-Eddy Simulations will not achieve their theoretical performance limits.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Singularity Formation in the Deterministic and Stochastic Fractional Burgers Equation

This study is motivated by the question of how singularity formation and other forms of extreme behavior in nonlinear dissipative partial differential equations are affected by stochastic excitations. To address this question we consider the 1D fractional Burgers equation with additive colored noise as a model problem. This system is interesting, because in the deterministic setting it exhibits finite-time blow-up or a globally well-posed behavior depending on the value of the fractional dissipation exponent. The problem is studied by performing a series of accurate numerical computations combining spectrally-accurate spatial discretization with a Monte-Carlo approach. First, we carefully document the singularity formation in the deterministic system in the supercritical regime where the blow-up time is shown to be a decreasing function of the fractional dissipation exponent. Our main result for the stochastic problem is that there is no evidence for the noise to regularize the evolution by suppressing blow-up in the supercritical regime, or for the noise to trigger blow-up in the subcritical regime. However, as the noise amplitude becomes large, the blow-up times in the supercritical regime are shown to exhibit an increasingly non-Gaussian behavior. Analogous observations are also made for the maximum attained values of the enstrophy and the times when the maxima occur in the subcritical regime.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Systematic Search For Extreme and Singular Behavior in Some Fundamental Models of Fluid Mechanics

This review article offers a survey of the research program focused on a systematic computational search for extreme and potentially singular behavior in hydrodynamic models motivated by open questions concerning the possibility of a finite-time blow-up in the solutions of the Navier-Stokes system. Inspired by the seminal work of Lu & Doering (2008), we sought such extreme behavior by solving PDE optimization problems with objective functionals chosen based on certain conditional regularity results and a priori estimates available for different models. No evidence for singularity formation was found in extreme Navier-Stokes flows constructed in this manner in 3D. We also discuss the results obtained for 1D Burgers and 2D Navier-Stokes systems, and while singularities are ruled out in these flows, the results presented provide interesting insights about sharpness of different energy-type estimates known for these systems. Connections to other bounding techniques are also briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Maximum Amplification of Enstrophy in 3D Navier-Stokes Flows

This investigation concerns a systematic search for potentially singular behavior in 3D Navier-Stokes flows. Enstrophy serves as a convenient indicator of the regularity of solutions to the Navier Stokes system --- as long as this quantity remains finite, the solutions are guaranteed to be smooth and satisfy the equations in the classical (pointwise) sense. However, there are no estimates available with finite a priori bounds on the growth of enstrophy and hence the regularity problem for the 3D Navier-Stokes system remains open. In order to quantify the maximum possible growth of enstrophy, we consider a family of PDE optimization problems in which initial conditions with prescribed enstrophy $\mathcal{E}_0$ are sought such that the enstrophy in the resulting Navier-Stokes flow is maximized at some time $T$. Such problems are solved computationally using a large-scale adjoint-based gradient approach derived in the continuous setting. By solving these problems for a broad range of values of $\mathcal{E}_0$ and $T$, we demonstrate that the maximum growth of enstrophy is in fact finite and scales in proportion to $\mathcal{E}_0^{3/2}$ as $\mathcal{E}_0$ becomes large. Thus, in such worst-case scenario the enstrophy still remains bounded for all times and there is no evidence for formation of singularity in finite time. We also analyze properties of the Navier-Stokes flows leading to the extreme enstrophy values and show that this behavior is realized by a series of vortex reconnection events.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Optimal Closures in a Simple Model for Turbulent Flows

In this work we introduce a computational framework for determining optimal closures of the eddy-viscosity type for Large-Eddy Simulations (LES) of a broad class of PDE models, such as the Navier-Stokes equation. This problem is cast in terms of PDE-constrained optimization where an error functional representing the misfit between the target and predicted observations is minimized with respect to the functional form of the eddy viscosity in the closure relation. Since this leads to a PDE optimization problem with a nonstandard structure, the solution is obtained computationally with a flexible and efficient gradient approach relying on a combination of modified adjoint-based analysis and Sobolev gradients. By formulating this problem in the continuous setting we are able to determine the optimal closure relations in a very general form subject only to some minimal assumptions. The proposed framework is thoroughly tested on a model problem involving the LES of the 1D Kuramoto-Sivashinsky equation, where optimal forms of the eddy viscosity are obtained as generalizations of the standard Smagorinsky model. It is demonstrated that while the solution trajectories corresponding to the DNS and LES still diverge exponentially, with such optimal eddy viscosities the rate of divergence is significantly reduced as compared to the Smagorinsky model. By systematically finding {optimal forms of the eddy viscosity within a certain general class of closure} models, thisframework can thus provide insights about the fundamental performance limitations of these models.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Rotating Equilibria of Vortex Sheets

We consider relative equilibrium solutions of the two-dimensional Euler equations in which the vorticity is concentrated on a union of finite-length vortex sheets. Using methods of complex analysis, more specifically the theory of the Riemann-Hilbert problem, a general approach is proposed to find such equilibria which consists of two steps: first, one finds a geometric configuration of vortex sheets ensuring that the corresponding circulation density is real-valued and also vanishes at all sheet endpoints such that the induced velocity field is well-defined; then, the circulation density is determined by evaluating a certain integral formula. As an illustration of this approach, we construct a family of rotating equilibria involving different numbers of straight vortex sheets rotating about a common center of rotation and with endpoints at the vertices of a regular polygon. This equilibrium generalizes the well-known solution involving single rotating vortex sheet. With the geometry of the configuration specified analytically, the corresponding circulation densities are obtained in terms of a integral expression which in some cases lends itself to an explicit evaluation. It is argued that as the number of sheets in the equilibrium configuration increases to infinity, the equilibrium converges in a certain distributional sense to a hollow vortex bounded by a constant-intensity vortex sheet, which is also a known equilibrium solution of the two-dimensional Euler equations.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Bayesian Uncertainty Quantification in Inverse Modelling of Electrochemical Systems

This study proposes a novel approach to quantifying uncertainties of constitutive relations inferred from noisy experimental data using inverse modelling. We focus on electrochemical systems in which charged species (e.g., Lithium ions) are transported in electrolyte solutions under an applied current. Such systems are typically described by the Planck-Nernst equation in which the unknown material properties are the diffusion coefficient and the transference number assumed constant or concentration-dependent. These material properties can be optimally reconstructed from time- and space-resolved concentration profiles measured during experiments using the Magnetic Resonance Imaging (MRI) technique. However, since the measurement data is usually noisy, it is important to quantify how the presence of noise affects the uncertainty of the reconstructed material properties. We address this problem by developing a state-of-the-art Bayesian approach to uncertainty quantification in which the reconstructed material properties are recast in terms of probability distributions, allowing us to rigorously determine suitable confidence intervals. The proposed approach is first thoroughly validated using "manufactured" data exhibiting the expected behavior as the magnitude of noise is varied. Then, this approach is applied to quantify the uncertainty of the diffusion coefficient and the transference number reconstructed from experimental data revealing interesting insights.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Harnessing the Kelvin-Helmholtz Instability: Feedback Stabilization of an Inviscid Vortex Sheet

In this investigation we use a simple model of the dynamics of an inviscid vortex sheet given by the Birkhoff-Rott equation to obtain fundamental insights about the potential for stabilization of shear layers using feedback control. As actuation we consider two arrays of point sinks/sources located a certain distance above and below the vortex sheet and subject to the constraint that their mass fluxes separately add up to zero. First, we demonstrate using analytical computations that the Birkhoff-Rott equation linearized around the flat-sheet configuration is in fact controllable when the number of actuator pairs is sufficiently large relative to the number of discrete degrees of freedom present in the system, a result valid for generic actuator locations. Next we design a state-based LQR stabilization strategy where the key difficulty is the numerical solution of the Riccati equation in the presence of severe ill-conditioning resulting from the properties of the Birkhoff-Rott equation and the chosen form of actuation, an issue which is overcome by performing computations with a suitably increased arithmetic precision. Analysis of the linear closed-loop system reveals exponential decay of the perturbation energy and of the corresponding actuation energy in all cases. Computations performed for the nonlinear closed-loop system demonstrate that initial perturbations of nonnegligible amplitude can be effectively stabilized when a sufficient number of actuators is used. We also thoroughly analyze the sensitivity of the closed-loop stabilization strategies to the variation of a number of key parameters. Subject to the known limitations of inviscid vortex models, our findings indicate that, in principle, it may be possible to stabilize shear layers for relatively large initial perturbations, provided the actuation has sufficiently many degrees of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Drift Due to Two Obstacles in Different Arrangements

We study the drift induced by the passage of two cylinders through an unbounded extent of inviscid incompressible fluid under the assumption that the flow is two-dimensional and steady in the moving frame of reference. The goal is to assess how the resulting total particle drift depends on the parameters of the geometric configuration, namely, the distance between the cylinders and their angle with respect to the direction of translation. This problem is studied by numerically computing, for different cylinder configurations, the trajectories of particles starting at various initial locations. The velocity field used in these computations is expressed in closed form using methods of the complex function theory and the accuracy of calculations is carefully verified. We identify cylinder configurations which result in increased and decreased drift with respect to the reference case when the two cylinders are separated by an infinite distance. Particle trajectories shed additional light on the hydrodynamic interactions between the cylinders in configurations resulting in different drift values. This ensemble of results provides insights about the accuracy of models used to study biogenic transport.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Linear Stability of Hill's Vortex to Axisymmetric Perturbations

We consider the linear stability of Hill's vortex with respect to axisymmetric perturbations. Given that Hill's vortex is a solution of a free-boundary problem, this stability analysis is performed by applying methods of shape differentiation to the contour dynamics formulation of the problem in a 3D axisymmetric geometry. This approach allows us to systematically account for the effect of boundary deformations on the linearized evolution of the vortex under the constraint of constant circulation. The resulting singular integro-differential operator defined on the vortex boundary is discretized with a highly accurate spectral approach. This operator has two unstable and two stable eigenvalues complemented by a continuous spectrum of neutrally-stable eigenvalues. By considering a family of suitably regularized (smoothed) eigenvalue problems solved with a range of numerical resolutions we demonstrate that the corresponding eigenfunctions are in fact singular objects in the form of infinitely sharp peaks localized at the front and rear stagnation points. These findings thus refine the results of the classical analysis by Moffatt & Moore (1978).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computation of Steady Incompressible Flows in Unbounded Domains

In this study we revisit the problem of computing steady Navier-Stokes flows in two-dimensional unbounded domains. Precise quantitative characterization of such flows in the high-Reynolds number limit remains an open problem of theoretical fluid dynamics. Following a review of key mathematical properties of such solutions related to the slow decay of the velocity field at large distances from the obstacle, we develop and carefully validate a spectrally-accurate computational approach which ensures the correct behavior of the solution at infinity. In the proposed method the numerical solution is defined on the entire unbounded domain without the need to truncate this domain to a finite box with some artificial boundary conditions prescribed at its boundaries. Since our approach relies on the streamfunction-vorticity formulation, the main complication is the presence of a discontinuity in the streamfunction field at infinity which is related to the slow decay of this field. We demonstrate how this difficulty can be overcome by reformulating the problem using a suitable background "skeleton" field expressed in terms of the corresponding Oseen flow combined with spectral filtering. The method is thoroughly validated for Reynolds numbers spanning two orders of magnitude with the results comparing favourably against known theoretical predictions and the data available in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal Nonlinear Eddy Viscosity in Galerkin Models of Turbulent Flows

We propose a variational approach to identification of an optimal nonlinear eddy viscosity as a subscale turbulence representation for POD models. The ansatz for the eddy viscosity is given in terms of an arbitrary function of the resolved fluctuation energy. This function is found as a minimizer of a cost functional measuring the difference between the target data coming from a resolved direct or large-eddy simulation of the flow and its reconstruction based on the POD model. The optimization is performed with a data-assimilation approach generalizing the 4D-VAR method. POD models with optimal eddy viscosities are presented for a 2D incompressible mixing layer at $Re=500$ (based on the initial vorticity thickness and the velocity of the high-speed stream) and a 3D Ahmed body wake at $Re=300,000$ (based on the body height and the free-stream velocity). The variational optimization formulation elucidates a number of interesting physical insights concerning the eddy-viscosity ansatz used. The 20-dimensional model of the mixing-layer reveals a negative eddy-viscosity regime at low fluctuation levels which improves the transient times towards the attractor. The 100-dimensional wake model yields more accurate energy distributions as compared to the nonlinear modal eddy-viscosity benchmark {proposed recently} by Östh et al. (2014). Our methodology can be applied to construct quite arbitrary closure relations and, more generally, constitutive relations optimizing statistical properties of a broad class of reduced-order models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal Reconstruction of Inviscid Vortices

We address the question of constructing simple inviscid vortex models which optimally approximate realistic flows as solutions of an inverse problem. Assuming the model to be incompressible, inviscid and stationary in the frame of reference moving with the vortex, the "structure" of the vortex is uniquely characterized by the functional relation between the streamfunction and vorticity. It is demonstrated how the inverse problem of reconstructing this functional relation from data can be framed as an optimization problem which can be efficiently solved using variational techniques. In contrast to earlier studies, the vorticity function defining the streamfunction-vorticity relation is reconstructed in the continuous setting subject to a minimum number of assumptions. To focus attention, we consider flows in 3D axisymmetric geometry with vortex rings. To validate our approach, a test case involving Hill's vortex is presented in which a very good reconstruction is obtained. In the second example we construct an optimal inviscid vortex model for a realistic flow in which a more accurate vorticity function is obtained than produced through an empirical fit. When compared to available theoretical vortex-ring models, our approach has the advantage of offering a good representation of both the vortex structure and its integral characteristics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Wake Effects on Drift in Two-Dimensional Inviscid Incompressible Flows

This investigation analyzes the effect of vortex wakes on the Lagrangian displacement of particles induced by the passage of an obstacle in a two-dimensional incompressible and inviscid fluid. In addition to the trajectories of individual particles, we also study their drift and the corresponding total drift areas in the Föppl and Kirchhoff potential flow models. Our findings, which are obtained numerically and in some regimes are also supported by asymptotic analysis, are compared to the wakeless potential flow which serves as a reference. We show that in the presence of the Föppl vortex wake some of the particles follow more complicated trajectories featuring a second loop. The appearance of an additional stagnation point in the Föppl flow is identified as a source of this effect. It is also demonstrated that, while the total drift area increases with the size of the wake for large vortex strengths, it is actually decreased for small circulation values. On the other hand, the Kirchhoff flow model is shown to have an unbounded total drift area. By providing a systematic account of the wake effects on the drift, the results of this study will allow for more accurate modeling of hydrodynamic stirring.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Framework for Linear Stability Analysis of Finite-Area Vortices

In this investigation we revisit the question of the linear stability analysis of 2D steady Euler flows characterized by the presence of compact regions with constant vorticity embedded in a potential flow. We give a complete derivation of the linearized perturbation equation which, recognizing that the underlying equilibrium problem is of the free-boundary type, is done systematically using methods of the shape-differential calculus. Particular attention is given to the proper linearization of the contour integrals describing vortex induction. The thus obtained perturbation equation is validated by analytically deducing from it the stability analyses of the circular vortex, originally due to Kelvin (1880), and of the elliptic vortex, originally due to Love (1893), as special cases. We also propose and validate a spectrally-accurate numerical approach to the solution of the stability problem for vortices of general shape in which all singular integrals are evaluated analytically.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Method for Geometry Optimization in a Simple Model of Two-Dimensional Heat Transfer

This investigation is motivated by the problem of optimal design of cooling elements in modern battery systems. We consider a simple model of two-dimensional steady-state heat conduction described by elliptic partial differential equations and involving a one-dimensional cooling element represented by a contour on which interface boundary conditions are specified. The problem consists in finding an optimal shape of the cooling element which will ensure that the solution in a given region is close (in the least squares sense) to some prescribed target distribution. We formulate this problem as PDE-constrained optimization and the locally optimal contour shapes are found using a gradient-based descent algorithm in which the Sobolev shape gradients are obtained using methods of the shape-differential calculus. The main novelty of this work is an accurate and efficient approach to the evaluation of the shape gradients based on a boundary-integral formulation which exploits certain analytical properties of the solution and does not require grids adapted to the contour. This approach is thoroughly validated and optimization results obtained in different test problems exhibit nontrivial shapes of the computed optimal contours.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Optimal Model Identification For Oscillatory Dynamics With a Stable Limit Cycle

We propose a general framework for parameter-free identification of a class of dynamical systems. Here, the propagator is approximated in terms of an arbitrary function of the state, in contrast to a polynomial or Galerkin expansion used in traditional approaches. The proposed formulation relies on variational data assimilation using measurement data combined with assumptions on the smoothness of the propagator. This approach is illustrated using a generalized dynamic model describing oscillatory transients from an unstable fixed point to a stable limit cycle and arising in nonlinear stability analysis as an example. This 3-state model comprises an evolution equation for the dominant oscillation and an algebraic manifold for the low- and high-frequency components in an autonomous descriptor system. The proposed optimal model identification technique employs mode amplitudes of the transient vortex shedding in a cylinder wake flow as example measurements. The reconstruction obtained with our technique features distinct and systematic improvements over the well-known mean-field (Landau) model of the Hopf bifurcation. The computational aspect of the identification method is thoroughly validated showing that good reconstructions can also be obtained in the absence of of accurate initial approximations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Maximum Palinstrophy Growth in 2D Incompressible Flows

In this study we investigate vortex structures which lead to the maximum possible growth of palinstrophy in two-dimensional incompressible flows on a periodic domain. The issue of palinstrophy growth is related to a broader research program focusing on extreme amplification of vorticity-related quantities which may signal singularity formation in different flow models. Such extreme vortex flows are found systematically via numerical solution of suitable variational optimization problems. We identify several families of maximizing solutions parameterized by their palinstrophy, palinstrophy and energy and palinstrophy and enstrophy. Evidence is shown that some of these families saturate estimates for the instantaneous rate of growth of palinstrophy obtained using rigorous methods of mathematical analysis, thereby demonstrating that this analysis is in fact sharp. In the limit of small palinstrophies the optimal vortex structures are found analytically, whereas for large palinstrophies they exhibit a self-similar multipolar structure. It is also shown that the time evolution obtained using the instantaneously optimal states with fixed energy and palinstrophy as the initial data saturates the upper bound for the maximum growth of palinstrophy in finite time. Possible implications of this finding for the questions concerning extreme behavior of flows are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Optimal Reconstruction of Material Properties in Complex Multiphysics Phenomena

We develop an optimization-based approach to the problem of reconstructing temperature-dependent material properties in complex thermo-fluid systems described by the equations for the conservation of mass, momentum and energy. Our goal is to estimate the temperature dependence of the viscosity coefficient in the momentum equation based on some noisy temperature measurements, where the temperature is governed by a separate energy equation. We show that an elegant and computationally efficient solution of this inverse problem is obtained by formulating it as a PDE-constrained optimization problem which can be solved with a gradient-based descent method. A key element of the proposed approach, the cost functional gradients are characterized by mathematical structure quite different than in typical problems of PDE-constrained optimization and are expressed in terms of integrals defined over the level sets of the temperature field. Advanced techniques of integration on manifolds are required to evaluate numerically such gradients, and we systematically compare three different methods. As a model system we consider a two-dimensional unsteady flow in a lid-driven cavity with heat transfer, and present a number of computational tests to validate our approach and illustrate its performance.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Vortices, Maximum Growth and the Problem of Finite-Time Singularity Formation

In this work we are interested in extreme vortex states leading to the maximum possible growth of palinstrophy in 2D viscous incompressible flows on periodic domains. This study is a part of a broader research effort motivated by the question about the finite-time singularity formation in the 3D Navier-Stokes system and aims at a systematic identification of the most singular flow behaviors. We extend the results reported in Ayala & Protas (2013) where extreme vortex states were found leading to the growth of palinstrophy, both instantaneously and in finite-time, which saturates the estimates obtained with rigorous methods of mathematical analysis. Here we uncover the vortex dynamics mechanisms responsible for such extreme behavior in time-dependent 2D flows. While the maximum palinstrophy growth is achieved at short times, the corresponding long-time evolution is characterized by some nontrivial features, such as vortex scattering events.

0 citations0 reviewsNo code linkedOpen access