Time-stepping approach for solving upper-bound problems: Application to two-dimensional Rayleigh-Benard convection
MetadataShow full item record
An alternative computational procedure for numerically solving a class of variational problems arising from rigorous upper-bound analysis of forced-dissipative infinite-dimensional nonlinear dynamical systems, including the Navier-Stokes and Oberbeck-Boussinesq equations, is analyzed and applied to Rayleigh-Benard convection. A proof that the only steady state to which this numerical algorithm can converge is the required global optimal of the relevant variational problem is given for three canonical flow configurations. In contrast with most other numerical schemes for computing the optimal bounds on transported quantities (e.g., heat or momentum) within the "background field" variational framework, which employ variants of Newton's method and hence require very accurate initial iterates, the new computational method is easy to implement and, crucially, does not require numerical continuation. The algorithm is used to determine the optimal background-method bound on the heat transport enhancement factor, i.e., the Nusselt number (Nu), as a function of the Rayleigh number (Ra), Prandtl number (Pr), and domain aspect ratio L in two-dimensional Rayleigh-Benard convection between stress-free isothermal boundaries (Rayleigh's original 1916 model of convection). The result of the computation is significant because analyses, laboratory experiments, and numerical simulations have suggested a range of exponents alpha and beta in the presumed Nu similar to (PrRa beta)-Ra-alpha scaling relation. The computations clearly show that for Ra <= 10(10) at fixed L = 2 root 2, Nu <= 0.106Pr(0)Ra(5/12), which indicates that molecular transport cannot generally be neglected in the "ultimate" high-Ra regime.
CitationWen, Baole, Gregory P. Chini, Rich R. Kerswell, and Charles R. Doering. "Time-stepping approach for solving upper-bound problems: Application to two-dimensional Rayleigh-Bénard convection." Physical Review E 92, no. 4 (2015): 043012.
Showing items related by title, author, creator and subject.
Montgomery, M. H. (2010-08)We provide multiple lines of evidence that GD 358's convection zone was much thinner in 1996 than in 2006 and we interpret this as a result of a short-lived increase in its surface temperature. In addition, we find strong ...
Montgomery, M. H.; Hermes, J. J.; Winget, D. E. (2010-08)This poster examines the nonlinear combination frequencies found in the DAV G117-B15A. Using multiple epochs of observations we show that the ratio of the amplitude of the combination frequencies to those of the parent ...
Evidence For Temperature Change And Oblique Pulsation From Light Curve Fits Of The Pulsating White Dwarf GD 358 Montgomery, Michael H.; Provencal, J. L.; Kanaan, A.; Mukadam, Anjum S.; Thompson, Susan E.; Dalessio, J.; Shipman, H. L.; Winget, D. E.; Kepler, S. O.; Koester, Detlev (2010-06)Convective driving, the mechanism originally proposed by Brickhill for pulsating white dwarf stars, has gained general acceptance as the generic linear instability mechanism in DAV and dbV white dwarfs. This physical ...