skip to main content

Title: Algebraic bounds on the Rayleigh–Bénard attractor
Abstract The Rayleigh–Bénard system with stress-free boundary conditions is shown to have a global attractor in each affine space where velocity has fixed spatial average. The physical problem is shown to be equivalent to one with periodic boundary conditions and certain symmetries. This enables a Gronwall estimate on enstrophy. That estimate is then used to bound the L 2 norm of the temperature gradient on the global attractor, which, in turn, is used to find a bounding region for the attractor in the enstrophy–palinstrophy plane. All final bounds are algebraic in the viscosity and thermal diffusivity, a significant improvement over previously established estimates. The sharpness of the bounds are tested with numerical simulations.
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Page Range or eLocation-ID:
509 to 531
Sponsoring Org:
National Science Foundation
More Like this
  1. Motivated by the search for sharp bounds on turbulent heat transfer as well as the design of optimal heat exchangers, we consider incompressible flows that most efficiently cool an internally heated disc. Heat enters via a distributed source, is passively advected and diffused, and exits through the boundary at a fixed temperature. We seek an advecting flow to optimize this exchange. Previous work on energy-constrained cooling with a constant source has conjectured that global optimizers should resemble convection rolls; we prove one-sided bounds on energy-constrained cooling corresponding to, but not resolving, this conjecture. In the case of an enstrophy constraint,more »our results are more complete: we construct a family of self-similar, tree-like ‘branching flows’ whose cooling we prove is within a logarithm of globally optimal. These results hold for general space- and time-dependent source–sink distributions that add more heat than they remove. Our main technical tool is a non-local Dirichlet-like variational principle for bounding solutions of the inhomogeneous advection–diffusion equation with a divergence-free velocity. This article is part of the theme issue ‘Mathematical problems in physical fluid dynamics (part 1)’.« less
  2. In inertial-range turbulence, structure functions can diagnose transfer or dissipation rates of energy and enstrophy, which are difficult to calculate directly in flows with complex geometry or sparse sampling. However, existing relations between third-order structure functions and these rates only apply under isotropic conditions. We propose new relations to diagnose energy and enstrophy dissipation rates in anisotropic two-dimensional (2-D) turbulence. These relations use second-order advective structure functions that depend on spatial increments of vorticity, velocity, and their advection. Numerical simulations of forced-dissipative anisotropic 2-D turbulence are used to compare new and existing relations against model-diagnosed dissipation rates of energy andmore »enstrophy. These simulations permit a dual cascade where forcing is applied at an intermediate scale, energy is dissipated at large scales, and enstrophy is dissipated at small scales. New relations to estimate energy and enstrophy dissipation rates show improvement over existing methods through increased accuracy, insensitivity to sampling direction, and lower temporal and spatial variability. These benefits of advective structure functions are present under weakly anisotropic conditions, and increase with the flow anisotropy as third-order structure functions become increasingly inappropriate. Several of the structure functions also show promise for diagnosing the forcing scale of 2-D turbulence. Velocity-based advective structure functions show particular promise as they can diagnose both enstrophy and energy cascade rates, and are robust to changes in the effective resolution of local derivatives. Some existing and future datasets that are amenable to advective structure function analysis are discussed.« less
  3. The large deflections of panels in subsonic flow are considered, specifically a fully clamped von Karman plate accounting for both rotational inertia in plate filaments and (mild) structural damping. The panel is taken to be embedded in the boundary of the positive half-space in ℝ3 containing a linear, subsonic potential flow. Solutions are constructed via a semigroup approach despite the lack of natural dissipativity associated with the generator of the linear dynamics. The flow–plate dynamics are then reduced—via an explicit Neumann-to-Dirichlet (downwash-to-pressure) solver for the flow—to a memory-type dynamical system for the plate. For the non-conservative plate dynamics, a globalmore »attractor is explicitly constructed via Lyapunov and recent quasi-stability methods. Finally, it is shown that, via the compactness of the attractor and finiteness of the dissipation integral, all trajectories converge strongly to the set of stationary states.« less
  4. Abstract. The high-latitude atmosphere is a dynamic region with processes that respond to forcing from the Sun, magnetosphere, neutral atmosphere, andionosphere. Historically, the dominance of magnetosphere–ionosphere interactions has motivated upper atmospheric studies to use magneticcoordinates when examining magnetosphere–ionosphere–thermosphere coupling processes. However, there are significant differences between thedominant interactions within the polar cap, auroral oval, and equatorward of the auroral oval. Organising data relative to these boundaries hasbeen shown to improve climatological and statistical studies, but the process of doing so is complicated by the shifting nature of the auroral ovaland the difficulty in measuring its poleward and equatorward boundaries. This studymore »presents a new set of open–closed magnetic field line boundaries (OCBs) obtained from Active Magnetosphere and Planetary ElectrodynamicsResponse Experiment (AMPERE) magnetic perturbation data. AMPERE observations of field-aligned currents (FACs) are used to determine the location ofthe boundary between the Region 1 (R1) and Region 2 (R2) FAC systems. This current boundary is thought to typically lie a few degrees equatorwardof the OCB, making it a good candidate for obtaining OCB locations. The AMPERE R1–R2 boundaries are compared to the Defense MeteorologicalSatellite Program Special Sensor J (DMSP SSJ) electron energy flux boundaries to test this hypothesis and determine the best estimate of thesystematic offset between the R1–R2 boundary and the OCB as a function of magnetic local time. These calibrated boundaries, as well as OCBsobtained from the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) observations, are validated using simultaneous observations of theconvection reversal boundary measured by DMSP. The validation shows that the OCBs from IMAGE and AMPERE may be used together in statisticalstudies, providing the basis of a long-term data set that can be used to separate observations originating inside and outside of the polar cap.« less
  5. We construct a determining form for the 2D Rayleigh-Benard (RB) system in a strip with solid horizontal boundaries, in the cases of no-slip and stress-free boundary conditions. The determining form is an ODE in a Banach space of trajectories whose steady states comprise the long-time dynamics of the RB system. In fact, solutions on the global attractor of the RB system can be further identified through the zeros of a scalar equation to which the ODE reduces for each initial trajectory. The twist in this work is that the trajectories are for the velocity field only, which in turn determinesmore »the corresponding trajectories of the temperature.« less