[논문 리뷰] Exhausting the background approach for bounding the heat transport in Rayleigh-Bénard convection

Zijing Ding, Rich R. Kerswell|arXiv (Cornell University)|2019. 06. 08.

Fluid Dynamics and Turbulent Flows참고 문헌 27인용 수 2

한 줄 요약

이 논문은 2D 레일리-베나르 대류에서 열전달 수(Nu)의 엄밀한 상한을 구하기 위해 배경 방법을 개선하여 열방정식과 시간 평균 브류시에 방정식을 완전히 충족시키는 2D 온도 및 속도 배경장에 기반한 방법을 사용한다. 비록 비대칭성과 코리올리력 제약 조건을 포함한 확장된 방법을 적용했음에도 불구하고, 상한은 여전히 Nu ≤ 0.055Ra^{1/2}에 머물러 있으며, 관측된 Ra^{1/3} 스케일링에 도달하지 못함으로써 이 문제에 대한 배경 방법의 한계에 도달했다.

ABSTRACT

We revisit the optimal heat transport problem for Rayleigh-B\'enard convection in which a rigorous upper bound on the Nusselt number, $Nu$, is sought as a function of the Rayleigh number $Ra$. Concentrating on the 2-dimensional problem with stress-free boundary conditions, we impose the full heat equation as a constraint for the bound using a novel 2-dimensional background approach thereby complementing the `wall-to-wall' approach of Hassanzadeh \etal \,(\emph{J. Fluid Mech.} extbf{751}, 627-662, 2014). Imposing the same symmetry on the problem, we find correspondence with their result for $Ra \leq Ra_c:=4468.8$ but, beyond that, the optimal fields complexify to produce a higher bound. This bound approaches that by a 1-dimensional background field as the length of computational domain $L ightarrow\infty$. On lifting the imposed symmetry, the optimal 2-dimensional temperature background field reverts back to being 1-dimensional giving the best bound $Nu\le 0.055Ra^{1/2}$ compared to $Nu \le 0.026Ra^{1/2}$ in the non-slip case. % We then show via an inductive bifurcation analysis that imposing the full time-averaged Boussinesq equations as constraints (by introducing 2-dimensional temperature {\em and} velocity background fields) is also unable to lower this bound. This then exhausts the background approach for the 2-dimensional (and by extension 3-dimensional) Rayleigh-Benard problem with the bound remaining stubbornly $Ra^{1/2}$ while data seems more to scale like $Ra^{1/3}$ for large $Ra$. % Finally, we show that adding a velocity background field to the formulation of Wen \etal\, (\emph{Phys. Rev. E.} extbf{92}, 043012, 2015), which is able to use an extra vorticity constraint due to the stress-free condition to lower the bound to $ Nu \le O(Ra^{5/12})$, also fails to improve the bound.

연구 동기 및 목표

• 2D 레일리-베나르 대류에서 개선된 배경 방법을 사용하여 열전달 수(Nu)의 엄밀한 상한을 구하는 것.
• 배경 방법에 전체 열방정식과 시간 평균 브류시 방정식을 제약 조건으로 도입함으로써 Nu의 상한을 향상시킬 수 있는지 테스트하는 것.
• 대칭성 가정이나 코리올리력 제약 조건을 도입함으로써 Nu ≤ 0.055Ra^{1/2} 이하의 상한을 낮출 수 있는지 조사하는 것.
• 배경 방법이 난류 대류에서 관측된 Nu ∼ Ra^{1/3} 스케일링을 예측하는 데 있어 본질적인 한계에 도달했는지 판단하는 것.

제안 방법

• 이 연구는 표준 투영을 초월하여 전체 열방정식을 제약 조건으로 삼는 2D 배경 방법을 사용한다.
• 최적의 유동 구조에서의 공간적 복잡성을 포착하기 위해 2D 온도 및 속도 배경장을 도입한다.
• 변분 최적화 프레임워크 내에서 운동량 및 열방정식을 제약 조건으로 삼기 위해 라그랑주 승수를 사용한다.
• 복잡성 증가에 따른 해의 안정성과 최적성 평가를 위해 인덕티브 분기 분석을 적용한다.
• 웬 등(2015)의 연구를 영감으로 삼아 수정된 라그랑주 승수 구조를 통해 코리올리력 제약 조건을 포함한 확장된 방법을 도입한다.
• 수치적 시간 적분 방법을 사용하여 변화를 시뮬레이션하지만, 2D에서는 다수의 안정 상태가 존재함에 따라 국소 최적해로의 수렴이 보장되지 않음을 입증한다.

실험 결과

연구 질문

• RQ1 배경 방법에 전체 열방정식을 제약 조건으로 도입함으로써 2D 레일리-베나르 대류에서 Nu의 상한을 향상시킬 수 있는가?
• RQ2 배경장에서 대칭성 가정을 제거하면 대칭 케이스보다 더 낮은 상한을 얻을 수 있는가?
• RQ3 확장된 배경 형식을 통해 코리올리력 제약 조건을 포함함으로써 Nu ≤ 0.055Ra^{1/2} 이하의 상한을 낮출 수 있는가?
• RQ4 배경 방법의 Ra^{1/2} 스케일링 한계는 본질적인가, 아니면 추가 제약 조건이 이를 깨뜨릴 수 있는가?
• RQ5 고도로 발전된 수식 구조를 사용했음에도 불구하고 배경 방법이 실험에서 관측된 Nu ∼ Ra^{1/3} 스케일링에 도달하지 못하는 이유는 무엇인가?

주요 결과

• Ra ≤ Rac ≈ 4468.8 범위에서는 대칭 2D 배경 방법이 하사잔자 등(2014)의 결과를 재현하지만, 이를 초월하면 최적 해의 복잡성 증가로 인해 상한이 증가한다.
• 도메인 길이 L → ∞일 때, 2D 상한은 1D 배경장의 상한에 점점 수렴하며, 이는 한계에서 1D 구조의 지배적 영향을 확인한다.
• 대칭성을 제거하면 최적의 2D 온도 배경장은 1D 프로파일로 복귀하여 Nu ≤ 0.055Ra^{1/2}의 상한을 얻는다.
• 시간 평균 브류시 방정식을 2D 속도 및 온도 배경장으로 완전히 충족시키는 조건을 도입한 후에도 상한은 여전히 Nu ≤ 0.055Ra^{1/2}에 머물러 있으며, 이는 방법의 고갈을 시사한다.
• 코리올리력 제약 조건을 포함한 속도 배경장을 도입함으로써 이전에 Nu ≤ O(Ra^{5/12})로 상한을 낮춘 바 있었지만, 추가로 개선되지 않아 Ra^{1/2} 스케일링이 방법의 경계임을 확인한다.
• 다수의 안정 상태가 초깃값에 따라 존재함에 따라 2D에서는 수치적 시간 적분이 최적 해로의 수렴을 보장하지 못한다.