English

Thermal Convection over Fractal Surfaces

Fluid Dynamics 2020-11-25 v3

Abstract

We use well resolved numerical simulations with the Lattice Boltzmann Method to study Rayleigh-B\'enard convection in cells with a fractal boundary in two dimensions for Pr=1Pr = 1 and Ra[107,1010]Ra \in \left[10^7, 10^{10}\right]. The fractal boundaries are functions characterized by power spectral densities S(k)S(k) that decay with wavenumber, kk, as S(k)kpS(k) \sim k^{p} (p<0p < 0). The degree of roughness is quantified by the exponent pp with p<3p < -3 for smooth (differentiable) surfaces and 3p<1-3 \le p < -1 for rough surfaces with Hausdorff dimension Df=12(p+5)D_f=\frac{1}{2}(p+5). By computing the exponent β\beta in power law fits NuRaβNu \sim Ra^{\beta}, where NuNu and RaRa are the Nusselt and the Rayleigh numbers for Ra[108,1010]Ra \in \left[10^8, 10^{10}\right], we observe that heat transport scaling increases with roughness over the top two decades of Ra[108,1010]Ra \in \left[10^8, 10^{10}\right]. For pp =3.0= -3.0, 2.0-2.0 and 1.5-1.5 we find β=0.288±0.005,0.329±0.006\beta = 0.288 \pm 0.005, 0.329 \pm 0.006 and 0.352±0.0110.352 \pm 0.011, respectively. We also observe that the Reynolds number, ReRe, scales as ReRaξRe \sim Ra^{\xi}, where ξ0.57\xi \approx 0.57 over Ra[107,1010]Ra \in \left[10^7, 10^{10}\right], for all pp used in the study. For a given value of pp, the averaged NuNu and ReRe are insensitive to the specific realization of the roughness.

Keywords

Cite

@article{arxiv.1908.10194,
  title  = {Thermal Convection over Fractal Surfaces},
  author = {Srikanth Toppaladoddi and Andrew J. Wells and Charles R. Doering and John S. Wettlaufer},
  journal= {arXiv preprint arXiv:1908.10194},
  year   = {2020}
}

Comments

15 pages, 13 figures

R2 v1 2026-06-23T10:57:57.203Z