Comparison of turbulence statistics in isothermal and non-isothermal large eddy simulations of supercritical carbon dioxide jets

Supercritical carbon dioxide is of interest in a wide range of engineering problems, including carbon capture, utilization, and storage as well as advanced cycles for power generation. Non-ideal variations in physical properties of supercritical carbon dioxide impact the physics of these systems. It is important to understand how drastic changes in thermodynamic properties influence these flow physics to aid and optimize the design of future technologies related to carbon capture and sequestration. In this study, we simulate turbulent supercritical carbon dioxide jets to gain a better understanding of these physics. Of particular interest is the impact of pseudo-boiling on supercritical flow dynamics. We use a second-order finite volume discretization method with adaptive mesh refinement as implemented in the reacting flow solver, PeleC, to perform a large eddy simulation of three turbulent jets of supercritical carbon dioxide. We use the Soave-Redlich-Kwong equation of state to close the system and more accurately incorporate the departure from ideal gas behavior into the turbulent flow physics. We find that the isothermal supercritical jet exhibits many similar flow characteristics compared to ideal gas round turbulent jets, with minor differences seen in the decay and spreading rate of the jet and in a noticeable anisotropy between resolved turbulent kinetic energy components. The non-isothermal jet excluding the pseudo-boiling point exhibits only small difference compared to the isothermal case. Crossing the pseudo-boiling point results in markedly different behavior, with evidence indicative of increased Kelvin-Helmholtz-like instabilities and much faster jet decay and disintegration. These factors impact the degree of mixing in the transition region of the jet, indicating a potential for larger heat transfer and more rapid combustion dynamics.