Simulation of a GOx-GCH4 Rocket Combustor and the Effect of the GEKO Turbulence Model Coefficients

In this study, a single injector methane-oxygen rocket combustor is numerically studied. The simulations included in this study are based on the hardware and experimental data from the Technical University of Munich. The focus is on the recently developed generalized k–ω turbulence model (GEKO) and the effect of its adjustable coefficients on the pressure and on wall heat flux profiles, which are compared with the experimental data. It was found that the coefficients of ‘jet’, ‘near-wall’, and ‘mixing’ have a major impact, whereas the opposite can be deduced about the ‘separation’ parameter Csep, which highly influences the pressure and wall heat flux distributions due to the changes in the eddy-viscosity field. The simulation results are compared with the standard k–ε model, displaying a qualitatively and quantitatively similar behavior to the GEKO model at a Csep equal to unity. The default GEKO model shows a stable performance for three oxidizer-to-fuel ratios, enhancing the reliability of its use. The simulations are conducted using two chemical kinetic mechanisms: Zhukov and Kong and the more detailed RAMEC. The influence of the combustion model is of the same order as the influence of the turbulence model. In general, the numerical results present a good or satisfactory agreement with the experiment, and both GEKO at Csep = 1 or the standard k–ε model can be recommended for usage in the CFD simulations of rocket combustion chambers, as well as the Zhukov–Kong mechanism in conjunction with the flamelet approach.

Towards Mechanistic Wall Heat Flux Partitioning Model for Fully Developed Nucleate Boiling

Mechanistic models developed to predict partial nucleate boiling are not adequate for fully developed nucleate boiling due to differences in the prevailing heat transfer governing mechanisms. In place of the mechanistic model, several empirical correlations and semi-mechanistic models have been proposed over the years for the prediction of fully developed nucleate boiling as presented in this study but they are unsuitable for use in computational fluid dynamics (CFD) code. Recently, the simulation of fully developed nucleate boiling has become much more practical because of advancement in a computational method that involves the coupling of the interface capturing method (for slug bubbles) with the Eulerian multi-fluid model (for dispersed spherical bubbles). Nonetheless, there is a need for a mechanistic closure law for the fully developed nucleate boiling phenomenon that would complement this advancement in CFD. Towards this end, a mechanistic wall heat flux partitioning model for fully developed nucleate boiling is proposed in this study. This model is predicated on the hypothesis that a high heat flux nucleate boiling is distinguished by the existence of a liquid macro-layer between the heated wall and the slug or elongated bubbles and that the macro-layer is interspersed with numerous high frequency nucleate small bubbles. With this hypothesis, the heat flux generated on the heated wall is partitioned into two parts: conduction heat transfer across the macro-layer liquid film thickness and evaporation heat flux of the microlayer of the nucleating small bubbles. The proposed model is validated against experimental data.