Source author record

Andrew Wynn

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

ResearcherUnclaimed source record

math.OC physics.flu-dyn math-ph math.AP math.MP math.FA

Catalog footprint

What is connected

7works

6topics

4close 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.

preprint2021arXiv

Analytical bounds on the heat transport in internally heated convection

We obtain an analytical bound on the mean vertical convective heat flux $\langle w T \rangle$ between two parallel boundaries driven by uniform internal heating. We consider two configurations, one with both boundaries held at the same constant temperature, and the other one with a top boundary held at constant temperature and a perfectly insulating bottom boundary. For the first configuration, Arslan et al. (J. Fluid Mech. 919:A15, 2021) recently provided numerical evidence that Rayleigh-number-dependent corrections to the only known rigorous bound $\langle w T \rangle \leq 1/2$ may be provable if the classical background method is augmented with a minimum principle stating that the fluid's temperature is no smaller than that of the top boundary. Here, we confirm this fact rigorously for both configurations by proving bounds on $\langle wT \rangle$ that approach $1/2$ exponentially from below as the Rayleigh number is increased. The key to obtaining these bounds are inner boundary layers in the background fields with a particular inverse-power scaling, which can be controlled in the spectral constraint using Hardy and Rellich inequalities. These allow for qualitative improvements in the analysis not available to standard constructions.

preprint2021arXiv

The background method: Theory and computations

The background method is a widely used technique to bound mean properties of turbulent flows rigorously. This work reviews recent advances in the theoretical formulation and numerical implementation of the method. First, we describe how the background method can be formulated systematically within a broader "auxiliary function" framework for bounding mean quantities, and explain how symmetries of the flow and constraints such as maximum principles can be exploited. All ideas are presented in a general setting and are illustrated on Rayleigh-Bénard convection between stress-free isothermal plates. Second, we review a semidefinite programming approach and a timestepping approach to optimizing bounds computationally, revealing that they are related to each other through convex duality and low-rank matrix factorization. Open questions and promising directions for further numerical analysis of the background method are also outlined.

preprint2019arXiv

Chordal decomposition in operator-splitting methods for sparse semidefinite programs

We employ chordal decomposition to reformulate a large and sparse semidefinite program (SDP), either in primal or dual standard form, into an equivalent SDP with smaller positive semidefinite (PSD) constraints. In contrast to previous approaches, the decomposed SDP is suitable for the application of first-order operator-splitting methods, enabling the development of efficient and scalable algorithms. In particular, we apply the alternating direction method of multipliers (ADMM) to solve decomposed primal- and dual-standard-form SDPs. Each iteration of such ADMM algorithms requires a projection onto an affine subspace, and a set of projections onto small PSD cones that can be computed in parallel. We also formulate the homogeneous self-dual embedding (HSDE) of a primal-dual pair of decomposed SDPs, and extend a recent ADMM-based algorithm to exploit the structure of our HSDE. The resulting HSDE algorithm has the same leading-order computational cost as those for the primal or dual problems only, with the advantage of being able to identify infeasible problems and produce an infeasibility certificate. All algorithms are implemented in the open-source MATLAB solver CDCS. Numerical experiments on a range of large-scale SDPs demonstrate the computational advantages of the proposed methods compared to common state-of-the-art solvers.

preprint2019arXiv

New bounds on the vertical heat transport for Bénard-Marangoni convection at infinite Prandtl number

We prove a new rigorous upper bound on the vertical heat transport for Bénard-Marangoni convection of a two- or three-dimensional fluid layer with infinite Prandtl number. Precisely, for Marangoni number $Ma \gg 1$ the Nusselt number $Nu$ is bounded asymptotically by $Nu \lesssim Ma^{2/7}(\ln Ma)^{-1/7}$. Key to our proof are a background temperature field with a hyperbolic profile near the fluid's surface, and new estimates for the coupling between temperature and vertical velocity.

preprint2017arXiv

Fast ADMM for Semidefinite Programs with Chordal Sparsity

Many problems in control theory can be formulated as semidefinite programs (SDPs). For large-scale SDPs, it is important to exploit the inherent sparsity to improve the scalability. This paper develops efficient first-order methods to solve SDPs with chordal sparsity based on the alternating direction method of multipliers (ADMM). We show that chordal decomposition can be applied to either the primal or the dual standard form of a sparse SDP, resulting in scaled versions of ADMM algorithms with the same computational cost. Each iteration of our algorithms consists of a projection on the product of small positive semidefinite cones, followed by a projection on an affine set, both of which can be carried out efficiently. Our techniques are implemented in CDCS, an open source add-on to MATLAB. Numerical experiments on large-scale sparse problems in SDPLIB and random SDPs with block-arrow sparse patterns show speedups compared to some common state-of-the-art software packages.

preprint2015arXiv

$β$-admissibility of observation and control operators for hypercontractive semigroups

We prove a Weiss conjecture on $β$-admissibility of control and observation operators for discrete and continuous $γ$-hypercontractive semigroups of operators, by representing them in terms of shifts on weighted Bergman spaces and using a reproducing kernel thesis for Hankel operators. Particular attention is paid to the case $γ=2$, which corresponds to the unweighted Bergman shift.

preprint2012arXiv

Observer design for systems with an energy-preserving non-linearity

Observer design is considered for a class of non-linear systems whose non-linear part is energy preserving. A strategy to construct convergent observers for this class of non-linear system is presented. The approach has the advantage that it is possible, via convex programming, to prove whether the constructed observer converges, in contrast to several existing approaches to observer design for non-linear systems. Finally, the developed methods are applied to the Lorenz attractor and to a low order model for shear fluid flow.