Skip to main content

A reduced-order partial differential equation model for the flow in a thermosyphon

  • E. A. BURROUGHS (a1), E. A. COUTSIAS (a2) and L. A. ROMERO (a3)

Flow in a closed-loop thermosyphon heated from below exhibits a sequence of bifurcations with increasing Grashof number. Using the Navier–Stokes equations in the Boussinesq approximation we have derived a model where, in the case of a slender circular loop, the first Fourier modes exactly decouple from all other Fourier modes, leaving a system of three coupled nonlinear partial differential equations that completely describes the flow in the thermosyphon. We have characterized the flow through two bifurcations, identifying stable periodic solutions for flows of Prandtl number greater than 18.5, a much lower value than predicted previously. Because of the quadratic nonlinearity in this system of equations, it is possible to find the global stability limit, and we have proved that it is identical to the first bifurcation point.

The numerical study of the model equations is based on a highly accurate Fourier–Chebyshev spectral method, combined with asymptotic analysis at the various bifurcation points. Three-dimensional computations with a finite element method computational fluid dynamics code (MPSalsa), are also pursued. All three approaches are in close agreement.

Recommend this journal

Email your librarian or administrator to recommend adding this journal to your organisation's collection.

Journal of Fluid Mechanics
  • ISSN: 0022-1120
  • EISSN: 1469-7645
  • URL: /core/journals/journal-of-fluid-mechanics
Please enter your name
Please enter a valid email address
Who would you like to send this to? *


Full text views

Total number of HTML views: 0
Total number of PDF views: 8 *
Loading metrics...

Abstract views

Total abstract views: 31 *
Loading metrics...

* Views captured on Cambridge Core between September 2016 - 18th November 2017. This data will be updated every 24 hours.