Rigorous upper limits on the vertical heat transport in two-dimensional Rayleigh-Bénard convection between stress-free isothermal boundaries are derived from the Boussinesq approximation of the Navier-Stokes equations. The Nusselt number Nu is bounded in terms of the Rayleigh number Ra according to Nu≤0.2891Ra(5/12) uniformly in the Prandtl number Pr. This scaling challenges some theoretical arguments regarding asymptotic high Rayleigh number heat transport by turbulent convection.

Obtaining pressure field data from particle image velocimetry (PIV) is an attractive technique in fluid dynamics due to its noninvasive nature. The application of this technique generally involves integrating the pressure gradient or solving the pressure Poisson equation using a velocity field measured with PIV. However, very little research has been done to investigate the dynamics of error propagation from PIV-based velocity measurements to the pressure field calculation. Rather than measure the error through experiment, we investigate the dynamics of the error propagation by examining the Poisson equation directly. We analytically quantify the error bound in the pressure field, and are able to illustrate the mathematical roots of why and how the Poisson equation based pressure calculation propagates error from the PIV data. The results show that the error depends on the shape and type of boundary conditions, the dimensions of the flow domain, and the flow type.

Rigorous bounds on heat transport are derived for thermal convection between stress-free horizontal plates. For three-dimensional Rayleigh-Bénard convection at infinite Prandtl number (Pr), the Nusselt number (Nu) is bounded according to Nu 0.28764Ra 5/12 where Ra is the standard Rayleigh number. For convection driven by a uniform steady internal heat source between isothermal boundaries, the spatially and temporally averaged (non-dimensional) temperature is bounded from below by T 0.6910R −5/17 in three dimensions at infinite Pr and by T 0.8473R −5/17 in two dimensions at arbitrary Pr, where R is the heat Rayleigh number proportional to the injected flux.

We derive an improved rigorous lower bound on the space and time averaged temperature T of an infinite Prandtl number Boussinesq fluid contained between isothermal no-slip boundaries thermally driven by uniform internal heating. A novel approach is used wherein a singular stable stratification is introduced as a perturbation to a non-singular background profile, yielding the estimatewhere R is the heat Rayleigh number. The analysis relies on a generalized Hardy-Rellich inequality that is proved in the appendix.

We consider the Boussinesq approximation for Rayleigh-Bénard convection perturbed by an additive noise and with boundary conditions corresponding to heating from below. In two space dimensions, with sufficient stochastic forcing in the temperature component and large Prandtl number P r > 0, we establish the existence of a unique ergodic invariant measure. In three space dimensions, we prove the existence of a statistically invariant state, and establish unique ergodicity for the infinite Prandtl Boussinesq system. Throughout this work we provide streamlined proofs of unique ergodicity which invoke an asymptotic coupling argument, a delicate usage of the maximum principle, and exponential martingale inequalities. Lastly, we show that the background method of Constantin-Doering [CD96] can be applied in our stochastic setting, and prove bounds on the Nusselt number relative to the unique invariant measure.

The dynamical core of an atmospheric general circulation model is engineered to satisfy a delicate balance between numerical stability, computational cost, and an accurate representation of the equations of motion. It generally contains either explicitly added or inherent numerical diffusion mechanisms to control the buildup of energy or enstrophy at the smallest scales. The diffusion fosters computational stability and is sometimes also viewed as a substitute for unresolved subgrid-scale processes. A particular form of explicitly added diffusion is horizontal divergence damping.In this paper a von Neumann stability analysis of horizontal divergence damping on a latitude-longitude grid is performed. Stability restrictions are derived for the damping coefficients of both second-and fourthorder divergence damping. The accuracy of the theoretical analysis is verified through the use of idealized dynamical core test cases that include the simulation of gravity waves and a baroclinic wave. The tests are applied to the finite-volume dynamical core of NCAR's Community Atmosphere Model version 5 (CAM5). Investigation of the amplification factor for the divergence damping mechanisms explains how small-scale meridional waves found in a baroclinic wave test case are not eliminated by the damping.

This work applies a continuous data assimilation scheme-a particular framework for reconciling sparse and potentially noisy observations to a mathematical model-to Rayleigh-Bénard convection at infinite or large Prandtl numbers using only the temperature field as observables. These Prandtl numbers are applicable to the earth's mantle and to gases under high pressure. We rigorously identify conditions that guarantee synchronization between the observed system and the model, then confirm the applicability of these results via numerical simulations. Our numerical experiments show that the analytically derived conditions for synchronization are far from sharp; that is, synchronization often occurs even when the conditions of our theorems are not met. We also develop estimates on the convergence of an infinite Prandtl model to a large (but finite) Prandtl number generated set of observations. Numerical simulations in this hybrid setting indicate that the mathematically rigorous results are accurate, but of practical interest only for extremely large Prandtl numbers.

No abstract