Source author record

Bruce Turkington

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

ResearcherUnclaimed source record

math-ph math.MP physics.flu-dyn cond-mat.stat-mech math.DS nlin.CD

Catalog footprint

What is connected

4works

6topics

3close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2016arXiv

Optimal thermalization in a shell model of homogeneous turbulence

We investigate the turbulence-induced dissipation of the large scales in a statistically homogeneous flow using an "optimal closure," which one of us (BT) has recently exposed in the context of Hamiltonian dynamics. This statistical closure employs a Gaussian model for the turbulent scales, with corresponding vanishing third cumulant, and yet it captures an intrinsic damping. The key to this apparent paradox lies in a clear distinction between true ensemble averages and their proxies, most easily grasped when one works directly with the Liouville equation rather than the cumulant hierarchy. We focus on a simple problem for which the optimal closure can be fully and exactly worked out: the relaxation arbitrarily far-from-equilibrium of a single energy shell towards Gibbs equilibrium in an inviscid shell model of 3D turbulence. The predictions of the optimal closure are validated against DNS and contrasted with those derived from EDQNM closure.

preprint2015arXiv

Coarse-graining two-dimensional turbulence via dynamical optimization

A model reduction technique based on an optimization principle is employed to coarse-grain inviscid, incompressible fluid dynamics in two dimensions. In this reduction the spectrally-truncated vorticity equation defines the microdynamics, while the macroscopic state space consists of quasi-equilibrium trial probability densities on the microscopic phase space, which are parameterized by the means and variances of the low modes of the vorticity. A macroscopic path therefore represents a coarse-grained approximation to the evolution of a nonequilibrium ensemble of microscopic solutions. Closure in terms of the vector of resolved variables, namely, the means and variances of the low modes, is achieved by minimizing over all feasible paths the time integral of their mean-squared residual with respect to the Liouville equation. The equations governing the optimal path are deduced from Hamilton-Jacobi theory. The coarse-grained dynamics derived by this optimization technique contains a scale-dependent eddy viscosity, modified nonlinear interactions between the low mode means, and a nonlinear coupling between the mean and variance of each low mode. The predictive skill of this optimal closure is validated quantitatively by comparing it against direct numerical simulations. These tests show that good agreement is achieved without adjusting any closure parameters.

preprint2012arXiv

An optimization principle for deriving nonequilibrium statistical models of Hamiltonian dynamics

A general method for deriving closed reduced models of Hamiltonian dynamical systems is developed using techniques from optimization and statistical estimation. As in standard projection operator methods, a set of resolved variables is selected to capture the slow, macroscopic behavior of the system, and the family of quasi-equilibrium probability densities on phase space corresponding to these resolved variables is employed as a statistical model. The macroscopic dynamics of the mean resolved variables is determined by optimizing over paths of these probability densities. Specifically, a cost function is introduced that quantifies the lack-of-fit of such paths to the underlying microscopic dynamics; it is an ensemble-averaged, squared-norm of the residual that results from submitting a path of trial densities to the Liouville equation. The evolution of the macrostate is estimated by minimizing the time integral of the cost function. The value function for this optimization satisfies the associated Hamilton-Jacobi equation, and it determines the optimal relation between the statistical parameters and the irreversible fluxes of the resolved variables, thereby closing the reduced dynamics. The resulting equations for the macroscopic variables have the generic form of governing equations for nonequilibrium thermodynamics, and they furnish a rational extension of the classical equations of linear irreversible thermodynamics beyond the near-equilibrium regime. In particular, the value function is a thermodynamic potential that extends the classical dissipation function and supplies the nonlinear relation between thermodynamics forces and fluxes.

preprint2010arXiv

Best-fit quasi-equilibrium ensembles: a general approach to statistical closure of underresolved Hamiltonian dynamics

A new method of deriving reduced models of Hamiltonian dynamical systems is developed using techniques from optimization and statistical estimation. Given a set of resolved variables that define a model reduction, the quasi-equilibrium ensembles associated with the resolved variables are employed as a family of trial probability densities on phase space. The residual that results from submitting these trial densities to the Liouville equation is quantified by an ensemble-averaged cost function related to the information loss rate of the reduction. From an initial nonequilibrium state, the statistical state of the system at any later time is estimated by minimizing the time integral of the cost function over paths of trial densities. Statistical closure of the underresolved dynamics is obtained at the level of the value function, which equals the optimal cost of reduction with respect to the resolved variables, and the evolution of the estimated statistical state is deduced from the Hamilton-Jacobi equation satisfied by the value function. In the near-equilibrium regime, or under a local quadratic approximation in the far-from-equilibrium regime, this best-fit closure is governed by a differential equation for the estimated state vector coupled to a Riccati differential equation for the Hessian matrix of the value function. Since memory effects are not explicitly included in the trial densities, a single adjustable parameter is introduced into the cost function to capture a time-scale ratio between resolved and unresolved motions. Apart from this parameter, the closed equations for the resolved variables are completely determined by the underlying deterministic dynamics.