BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Bartosz Protas

Bartosz Protas contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynmath.APmath.OCphysics.comp-phmath-phmath.CVmath.MPmath.NAmath.PRNumerical Analysisphysics.chem-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 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