Full Content is available to subscribers

Subscribe/Learn More  >

A Numerical Study of Thermoacoustic Oscillations in a Rectangular Channel Using CMSIP Method

[+] Author Affiliations
K. M. Akyuzlu

University of New Orleans, New Orleans, LA

K. Albayrak, C. Karaeren

Middle East Technical University, Ankara, Turkey

Paper No. IMECE2009-13109, pp. 943-953; 11 pages
  • ASME 2009 International Mechanical Engineering Congress and Exposition
  • Volume 9: Heat Transfer, Fluid Flows, and Thermal Systems, Parts A, B and C
  • Lake Buena Vista, Florida, USA, November 13–19, 2009
  • Conference Sponsors: ASME
  • ISBN: 978-0-7918-4382-6 | eISBN: 978-0-7918-3863-1
  • Copyright © 2009 by ASME


This paper presents a mathematical model that was developed to study instabilities (primarily thermoacoustic oscillations) experienced inside a channel (with a rectangular cross section) heated symmetrically (from its top and bottom.) The heated channel is configured to simulate a combustion chamber of a rocket hybrid rocket motor and is connected to a converging–diverging nozzle in the downstream and to a plenum with a flow straightener in the upstream side. The working fluid is supplied from a pressurized storage tank to the upstream plenum through a throttle valve. A multi-component approach is used to model this test apparatus. In this integrated component model, the unsteady flow through the throttle valve and the nozzle is assumed to be one-dimensional and isentropic where as the flow in the forward plenum and the heated channel is assumed to be a two-dimensional, unsteady, compressible, turbulent, and subsonic. The physics based mathematical model of the flow in the channel consists of conservation of mass, momentum (two-dimensional Navier-Stokes) and energy equations subject to appropriate boundary conditions as defined by the physical problem stated above. The working fluid is assumed to be compressible where the density of the fluid is related to the pressure and temperature of the fluid through a simple ideal gas relation. The governing equations are discretized using second order accurate central differencing for spatial derivatives and second order accurate (based on Taylor expansion) finite difference approximations for temporal derivatives. The resulting nonlinear equations are then linearized using Newton’s linearization method. The set of algebraic equations that result from this process are then put into a matrix form and solved using a Coupled Modified Strongly Implicit Procedure (CMSIP) for the unknowns (primitive variables, i.e., pressure, temperature, and the velocity field) of the problem. The turbulence model equations and the unsteady flow equation for the throttle valve are solved using a second order accurate explicit finite difference technique. Convergence and grid independence studies were done to determine the optimum mesh size and computational time increment. Furthermore, two benchmark cases (unsteady driven cavity and laminar channel flows) were simulated using the developed numerical model to verify the accuracy of the proposed solution procedure. Numerical experiments were then carried out to simulate the thermoacoustic oscillations inside rectangular channels with various aspect ratios ranging from 5 to 20 for various operating conditions (i.e., for Re numbers between 102 and 106 ) and to determine the flow regions where these oscillations are sustained. The numerical simulation results indicate that the mathematical model for the gas flow in the heated channel predicts the expected unsteady temperature and pressure distributions, and the velocity field, successfully. Furthermore, it is concluded that the proposed integrated component model is successful in generating the characteristics of the instabilities associated with thermal, hydrodynamic, and thermoacoustic oscillations in heated channels.

Copyright © 2009 by ASME



Interactive Graphics


Country-Specific Mortality and Growth Failure in Infancy and Yound Children and Association With Material Stature

Use interactive graphics and maps to view and sort country-specific infant and early dhildhood mortality and growth failure data and their association with maternal

Citing articles are presented as examples only. In non-demo SCM6 implementation, integration with CrossRef’s "Cited By" API will populate this tab (http://www.crossref.org/citedby.html).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In