<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1129469</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2023.1129469</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Heat transfer of supercritical pressure water in helical tubes</article-title>
<alt-title alt-title-type="left-running-head">Liu et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenrg.2023.1129469">10.3389/fenrg.2023.1129469</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Mengmeng</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/2179607/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhang</surname>
<given-names>Zhen</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1980696/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yang</surname>
<given-names>Xingtuan</given-names>
</name>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jiang</surname>
<given-names>Shengyao</given-names>
</name>
</contrib>
</contrib-group>
<aff>
<institution>Institute of Nuclear and New Energy Technology</institution>, <institution>Collaborative Innovation Center of Advanced Nuclear Energy Technology</institution>, <institution>Key Laboratory of Advanced Nuclear Reactor Engineering and Safety of Ministry of Education</institution>, <institution>Tsinghua University</institution>, <addr-line>Beijing</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/874360/overview">Deqi Chen</ext-link>, Chongqing University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1967053/overview">Hongna Zhang</ext-link>, Tianjin University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/78833/overview">Jiang Qin</ext-link>, Harbin Institute of Technology, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Zhen Zhang, <email>zhangzhen13@tsinghua.edu.cn</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Nuclear Energy, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>15</day>
<month>02</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1129469</elocation-id>
<history>
<date date-type="received">
<day>22</day>
<month>12</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>27</day>
<month>01</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Liu, Zhang, Yang and Jiang.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Liu, Zhang, Yang and Jiang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Subcritical steam generators in the current high-temperature gas-cooled reactor pebble-bed module in China can be replaced by supercritical steam generators to better coordinate the reactors and supercritical steam turbine unit to improve the thermal efficiency with no fluid phase change at supercritical pressures. The heat transfer of supercritical-pressure water flow in helical tubes with various helix diameters was simulated using the re-normalization group <italic>k</italic>&#x2013;<italic>&#x3b5;</italic> model to investigate the effects of finite-thickness walls, inlet mass flux, and helix diameter. The finite-thickness wall barely affected local heat transfer in the straight tube with inlet mass flux of 1,260&#xa0;kg/(m<sup>2</sup> s) whereas the local heat transfer in the helical tube was very different with relative difference in <italic>q</italic>/<italic>q</italic>
<sub>av</sub> between the inside and bottom, caused by the centrifugal force, is 36%. In helical tubes, the heat transfer coefficients in the inner bar were the lowest, with a low fluid velocity under the effect of the centrifugal force. The buoyancy and centrifugal force in the helical tubes caused non-uniformity in the wall temperature; however, when the fluid temperature increased below the pseudo-critical temperature, the temperature distributions showed an increasing trend. The non-uniformities in the wall temperature distributions decreased as the centrifugal force increased, and they were suppressed by the finite-thickness wall.</p>
</abstract>
<kwd-group>
<kwd>supercritical-pressure water</kwd>
<kwd>helical tube</kwd>
<kwd>finite-thickness wall</kwd>
<kwd>helix diameter</kwd>
<kwd>non-uniform heat transfer</kwd>
</kwd-group>
<contract-sponsor id="cn001">State Administration for Science, Technology and Industry for National Defense<named-content content-type="fundref-id">10.13039/501100010878</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">China National Nuclear Corporation<named-content content-type="fundref-id">10.13039/501100014883</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>The high-temperature gas-cooled reactor pebble-bed module (HTR-PM) being built at Tsinghua University&#x2019;s Institute of Nuclear and New Energy Technology in China is a typical Generation IV nuclear power system. The subcritical steam generators in the HTR-PM can be replaced by supercritical steam generators to work with the current reactors and a supercritical steam turbine to improve the thermal efficiency, as envisaged for the third technical phase of the HTR-PM program in China (<xref ref-type="bibr" rid="B22">Zhang et al., 2019a</xref>). By raising the pressure of the helical coil steam generator at the outlet to supercritical levels, the thermal efficiency of the HTR-PM can be significantly improved. Small variations in fluid temperature and pressure, especially when water is at supercritical pressures in a steam generator, can cause drastic changes in the thermophysical properties, of which heat transfer characteristics (HTCs) are critical.</p>
<p>
<xref ref-type="fig" rid="F1">Figure 1</xref> shows the variations in the thermophysical properties of supercritical pressure water (SCW) at 24.5&#xa0;MPa, where the specific heat, <italic>C</italic>
<sub>
<italic>p</italic>
</sub>, reaches a sharp peak at a specific temperature (656.2&#xa0;K), defined as the pseudo-critical temperature. The pseudo-critical enthalpy was 2,145&#xa0;kJ/kg at the pseudo-critical temperature. At 24.5 MPa, the specific heat performance, <italic>C</italic>
<sub>
<italic>p</italic>
</sub>, and volume expansivity, <italic>&#x3b1;</italic>
<sub>
<italic>p</italic>
</sub>, were non-monotonic with change in temperature, varying dramatically near the pseudo-critical temperature and reaching the maximum at that point. The dynamic viscosity, <italic>&#x3bc;</italic>, and density, <italic>&#x3c1;</italic>, gradually decreased with increasing temperature, and rapidly decreased near the pseudo-critical temperature, after which the rate of decrease tended to be gentle. The variation in the thermal conductivity, <italic>&#x3bb;</italic>, was similar to those of <italic>&#x3bc;</italic> and <italic>&#x3c1;</italic>, except for a slight fluctuation near the pseudo-critical temperature.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Variations in thermophysical properties of water (<italic>p</italic> &#x3d; 24.5&#xa0;MPa, <italic>T</italic>
<sub>pc</sub> &#x3d; 656.2&#xa0;K).</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g001.tif"/>
</fig>
<p>A drastic change in the properties can significantly affect the HTCs, complicating the related processes. Many researchers have conducted experiments and numerical simulations to investigate the HTCs of supercritical fluids in tubes.</p>
<p>Experimental studies can provide reliable data on supercritical fluid heat transfer. <xref ref-type="bibr" rid="B18">Yamagata et al. (1972)</xref> investigated the HTCs of SCW in horizontal and vertical tubes at pressures ranging from 22.6 to 29.4 MPa, with a heat fluxes ranging from 116 to 930&#xa0;kW/m<sup>2</sup> and inlet mass fluxes ranging from 310 to 1,830&#xa0;kg/(m<sup>2</sup>&#xb7;s). The extensive experimental data reported by them provide the foundation for further research. <xref ref-type="bibr" rid="B3">Fewster and Jackson. (2004)</xref> conducted experiments in a vertical tube with an inner diameter of 5&#xa0;mm and found that when the heat flux was high, there was local heat transfer deterioration at each position in the experimental section. They identified three heat transfer states of supercritical fluid: normal, heat transfer deterioration, and enhancement. The current experimental data were summarized to obtain an empirical correlation under specific working conditions. <xref ref-type="bibr" rid="B13">Mokry et al. (2010)</xref> measured the temperature field of SCW under 24&#xa0;MPa, with a mass flux range of 200&#x2013;1,500&#xa0;kg/(m<sup>2</sup>&#xb7;s), heat flux range of 160&#x2013;900&#xa0;kW/m<sup>2</sup> and inlet temperature range of 593.15&#x2013;613.15&#xa0;K. Their experiments also revealed the three heat-transfer states reported by <xref ref-type="bibr" rid="B3">Fewster and Jackson. (2004)</xref>; <xref ref-type="bibr" rid="B20">Zhang et al. (2015a)</xref> conducted experiments to study the HTCs of supercritical CO<sub>2</sub> in vertical helical tubes and analyzed the enhancement and inhibition mechanisms of the non-dimensional buoyancy number, Bo&#x2a;. <xref ref-type="bibr" rid="B15">Sheeba et al. (2019)</xref> experimentally investigated the HTCs of SCW in a helical tube.</p>
<p>Numerical simulations provide more detailed temperature field than experiments. Some researchers (<xref ref-type="bibr" rid="B21">Zhang et al., 2015b</xref>; <xref ref-type="bibr" rid="B16">Wang et al., 2015</xref>; <xref ref-type="bibr" rid="B6">Huang et al., 2018</xref>; <xref ref-type="bibr" rid="B19">Yang et al., 2021</xref>) studied the HTCs of supercritical-pressure fluid in helical tubes through simulations, in which the wall thickness of the tube was ignored and a uniform heat flux was loaded on the wall. <xref ref-type="bibr" rid="B21">Zhang et al. (2015b)</xref> calculated the HTCs of SCW using a numerical model with the re-normalization group (RNG) <italic>k</italic>&#x2013;<italic>&#x3b5;</italic> turbulent model coupled to an enhanced wall function. They compared the results with the experimental data reported by <xref ref-type="bibr" rid="B18">Yamagata et al. (1972)</xref>, and investigated the effects of the wall (including those of the wall properties and thicknesses) and fluid properties. <xref ref-type="bibr" rid="B19">Yang et al. (2021)</xref> used the RNG k&#x2013;&#x3b5; turbulent model with an enhanced wall function to investigate the mechanisms of deterioration and oscillation of the inner wall temperature during SCW cooling in horizontal tubes with sections of various non-circular shapes. <xref ref-type="bibr" rid="B6">Huang et al. (2018)</xref> studied the HTCs of SCW and CO<sub>2</sub> in a helical tube using a numerical calculation model, and discovered that the SST k&#x2013;&#x3c9; turbulent model agreed well with the experimental results of <xref ref-type="bibr" rid="B17">Wang et al. (2014)</xref>, who demonstrated that the maximum wall temperature was at a certain angle in the circumferential direction under the effects of centrifugal force and buoyancy in the helical tube. <xref ref-type="bibr" rid="B16">Wang et al. (2015)</xref> studied the applicability of various turbulence models for heat transfer calculations in helical tubes using the experimental data from convective heat transfer of supercritical CO<sub>2</sub> in helical tubes.</p>
<p>Ignoring the wall thickness and loading a uniform heat flux makes the analytical conditions different from those of the actual operation. Moreover, the circumferential non-uniformities at the cross section may deviate from the actual value when the wall thickness is ignored. To study the effect of heat transfer on SCW under non-uniform heating, researchers (<xref ref-type="bibr" rid="B11">Li et al., 2014</xref>; <xref ref-type="bibr" rid="B14">Rowinski et al., 2017</xref>; <xref ref-type="bibr" rid="B10">Li and Bai, 2018</xref>) loaded axial and circumferential non-uniform heat fluxes on the wall. It was found that a non-uniform heat flux on the wall affected the distribution of buoyancy, leading to the earlier appearance of wall temperature peaks and affecting the local HTCs.</p>
<p>In summary, the effects of the wall on the HTCs of helical tubes need further investigation. Here, the local HTCs of SCW in helical tubes with various helix diameters under 24.5&#xa0;MPa were investigated and the effects of finite-thickness wall, inlet mass flux, and helix diameter were analyzed using the computational fluid dynamics software package ANSYS Fluent. The results could provide technical guidance for the design and safe operation of supercritical steam generators in HTR-PM.</p>
</sec>
<sec id="s2">
<title>2 Numerical simulation</title>
<sec id="s2-1">
<title>2.1 Physical model</title>
<p>
<xref ref-type="fig" rid="F2">Figure 2</xref> shows the tube structures and circumferential angle <italic>&#x3b8;</italic> of the cross section in the helical tubes, where the values of <italic>&#x3b8;</italic> at the top, inside, bottom, and outside are 0&#xb0;, 90&#xb0;, 180&#xb0;, and 270&#xb0;, respectively. <xref ref-type="fig" rid="F2">Figure 2A</xref> shows a straight tube inclined at 7.16&#xb0;, with inner diameter of 7.5&#xa0;mm, heated length of 6,000&#xa0;mm, and adiabatic sections (500&#xa0;mm long) at the inlet and outlet. The helical tubes shown in <xref ref-type="fig" rid="F2">Figure 2B</xref> have a heated length of 6,000&#xa0;mm and adiabatic sections (500&#xa0;mm long) at the inlet and outlet. The helical tubes with inner diameter of 7.5&#xa0;mm are labeled A, B, and C, and the corresponding helix diameters <italic>D</italic> are 70.5&#xa0;mm, 105.5&#xa0;mm, and 140.5&#xa0;mm; all tubes have the same inclination angle of 7.16&#xb0;. The finite-thickness walls of the tubes were of 2&#xa0;mm thickness and made of Incoloy 800H, which is the key material of the HTR-PM subcritical steam generator tubes; the material density <italic>&#x3c1;</italic>
<sub>w</sub> was 8,030&#xa0;kg/m<sup>3</sup>, specific heat <italic>C</italic>
<sub>p,w</sub> was 502&#xa0;J/(kg&#xb7;K), and thermal conductivity <italic>&#x3bb;</italic>
<sub>w</sub> ranged between 11.6&#x2013;25.1&#xa0;W/(m&#xb7;K) depending on the temperature. The fluid in all tubes flowed from bottom to top.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Structures and angular coordinates of the investigated tubes. <bold>(A)</bold> 7.16&#xb0; inclined straight tube; <bold>(B)</bold> Helical tubes.</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g002.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F2">Figure 2</xref> shows the circumferential angle <italic>&#x3b8;</italic> of the cross section, where the value of <italic>&#x3b8;</italic> at the top was 0&#xb0; and those at the sides were 90&#xb0; and 270&#xb0; (owing to the symmetry of straight tubes); in the helical tubes, the values of <italic>&#x3b8;</italic> on the inside, outside, and at the bottom were 90&#xb0;, 270&#xb0;, and 180&#xb0;, respectively.</p>
<p>The parameters applied in the simulation cases are with various inlet mass flux in straight and helical tubes with zero-thickness (no wall)/2&#xa0;mm-thickness wall. SCW flowed upward in the tubes at a pressure of 24.5&#xa0;MPa. The inlet temperature was 586.02&#xa0;K, and the heat fluxes on the inner wall were all 100&#xa0;kW/m<sup>2</sup>. In all the simulations of the present work, the heat flux (heat flux in tubes with no wall or volumetric heat source in tubes with 2-mm-thickness wall) applied on the walls are assumed constant and uniform. The inlet mass flux under Con. I, Con. II, and Con. III were 225&#xa0;kg/(m<sup>2</sup>&#xb7;s), 338&#xa0;kg/(m<sup>2</sup>&#xb7;s), and 1,260&#xa0;kg/(m<sup>2</sup>&#xb7;s), respectively.</p>
</sec>
<sec id="s2-2">
<title>2.2 Numerical method</title>
<p>This study utilizes the finite volume method with the RNG <italic>k-&#x3b5;</italic> model with enhanced wall treatment (see our previous research (<xref ref-type="bibr" rid="B9">Lemmon et al., 2007</xref>; <xref ref-type="bibr" rid="B21">Zhang et al., 2015b</xref>)), taking various thermophysical properties of SCW into account. The NIST standard database (<xref ref-type="bibr" rid="B23">Zhang et al., 2019b</xref>) is used to determine the physical properties of the materials during the calculation process.</p>
<p>The mass, momentum, and energy conservation equations are as follows.</p>
<p>The mass-conservation is calculated as<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3c1;</mml:mi>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>The momentum -conservation is calculated as<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3c1;</mml:mi>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3c1;</mml:mi>
<mml:msub>
<mml:mi>g</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>where <inline-formula id="inf1">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>T</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf2">
<mml:math id="m4">
<mml:mrow>
<mml:mi>&#x3bc;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the molecular viscosity, <inline-formula id="inf3">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>T</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the turbulent viscosity, <inline-formula id="inf4">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>T</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3c1;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>&#x3bc;</mml:mi>
</mml:msub>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>/</mml:mo>
<mml:mi>&#x3b5;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, and <italic>C</italic>
<sub>
<italic>&#x3bc;</italic>
</sub> &#x3d; 0.0845.</p>
<p>The energy conservation equation is<disp-formula id="e3">
<mml:math id="m7">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>E</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>T</mml:mi>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>Pr</mml:mi>
<mml:mi>T</mml:mi>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>where <italic>E</italic> is the total energy, and <inline-formula id="inf5">
<mml:math id="m8">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mi>T</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>p</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi>&#x3c1;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>The RNG model is used with the enhanced wall treatment to study 3-D HTCs in straight and helical tubes. The turbulent transport equations are as follows.</p>
<p>Turbulent kinetic energy, <italic>k</italic>, transport is calculated as (<xref ref-type="bibr" rid="B2">Durbin, 1991</xref>; <xref ref-type="bibr" rid="B1">Behnia et al., 1998</xref>).<disp-formula id="e4">
<mml:math id="m9">
<mml:mrow>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>k</mml:mi>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>T</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>T</mml:mi>
</mml:msub>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>&#x3b5;</mml:mi>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>The transport equation of turbulent kinetic energy dissipation rate, <italic>&#x3b5;</italic> (<xref ref-type="bibr" rid="B2">Durbin, 1991</xref>; <xref ref-type="bibr" rid="B1">Behnia et al., 1998</xref>), is<disp-formula id="e5">
<mml:math id="m10">
<mml:mrow>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>&#x3b5;</mml:mi>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>T</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3b5;</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>T</mml:mi>
</mml:msub>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3b5;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3c1;</mml:mi>
<mml:msup>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mi>k</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>where<disp-formula id="e6">
<mml:math id="m11">
<mml:mrow>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1.42</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b7;</mml:mi>
<mml:mo>/</mml:mo>
<mml:mn>4.38</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.015</mml:mn>
<mml:msup>
<mml:mi>&#x3b7;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
<disp-formula id="e8">
<mml:math id="m13">
<mml:mrow>
<mml:mi>&#x3b7;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mi>&#x3b5;</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
<disp-formula id="e9">
<mml:math id="m14">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
<disp-formula id="e10">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
<disp-formula id="equ1">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3b5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1.68</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi mathvariant="normal">&#x3b5;</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.7179</mml:mn>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:math>
</disp-formula>
</p>
<p>The SIMPLEC algorithm was used in the simulation to determine the pressure-velocity coupling in the domain. For spatial discretization, the discretization settings were set to second-order upwind. The convergence criterion was set to 1 &#xd7; 10<sup>&#x2212;6</sup> for continuity, velocity, energy, turbulent kinetic energy, and turbulent dissipation rate.</p>
</sec>
<sec id="s2-3">
<title>2.3 Grid independence verification</title>
<p>The structured 3-D meshes are refined for the simulation. Away from the near-wall region, the grid size is increased by 1.1 toward the tube center. The distance between the first node and the wall is determined using the RNG <italic>k-&#x3b5;</italic> model with enhanced wall treatment, and the non-dimensional wall distance, y<sup>
<italic>&#x2b;</italic>
</sup>, is set to less than 0.2 to ensure y<sup>
<italic>&#x2b;</italic>
</sup> independence.</p>
<p>Because the helical tube structures are more complex than straight tube structures, the 7.16&#xb0; inclined helical tube is selected for grid independence analysis. To obtain grid-independent results, the number of circumferential and axial grids is refined, with the number of elements ranging from 1,026,300 to 7,694,400, numbered as Grid-1 to Grid-4. The outlet heat transfer coefficients are compared using various grids. The enthalpy <italic>i</italic>
<sub>b</sub>(<italic>l</italic>), bulk temperature <italic>T</italic>
<sub>b</sub>(<italic>l</italic>), and local heat transfer coefficients <italic>h</italic>(<italic>l</italic>) can be calculated as follows.</p>
<p>The enthalpy equation is expressed as<disp-formula id="e11">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>l</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>where <italic>i</italic>
<sub>b,in</sub> is the enthalpy at the heating-section inlet, J/kg, <italic>q</italic>
<sub>w</sub> is the inner-wall heat flux, W/m<sup>2</sup>, <italic>l</italic> is the length of the heating section, m, and <italic>G</italic> is the mass flux, kg/(m<sup>2</sup>&#xb7;s).</p>
<p>
<italic>T</italic>
<sub>b</sub> is calculated based on the local fluid enthalpy and pressure using the NIST standard database (<xref ref-type="bibr" rid="B23">Zhang et al., 2019b</xref>). The local heat transfer coefficients can be expressed as follows.<disp-formula id="e12">
<mml:math id="m18">
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi mathvariant="normal">b</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>where <italic>T</italic>
<sub>
<italic>w</italic>
</sub> is the local wall temperature, K.</p>
<p>The average heat transfer coefficients at the 7.16&#xb0; inclined helical tube outlet of Grid-1 to Grid-3 are compared with that of Grid-4 with the maximum number of grids. The relative error, <italic>&#x3b6;</italic>, is defined as follows.<disp-formula id="e13">
<mml:math id="m19">
<mml:mrow>
<mml:mi>&#x3b6;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mo>,</mml:mo>
<mml:mn>4</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mo>,</mml:mo>
<mml:mn>4</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<p>The relative errors <italic>&#x3b6;</italic> between Grid-1, Grid-2, Grid-3 and Grid-4 are 26%, 7% and 0.8%, respectively. With the encryption of circumferential and axial grids, the value of <italic>&#x3b6;</italic> decreases. Considering accuracy and calculation time, Grid-3 is selected, with 5,665,950 elements.</p>
</sec>
<sec id="s2-4">
<title>2.4 Model verification</title>
<p>
<xref ref-type="bibr" rid="B18">Yamagata et al. (1972)</xref> investigated the changes at the inner wall temperature with water enthalpy in a vertical tube (diameter: 7.5&#xa0;mm) with uniform heat flux of 233 and 465&#xa0;kW/m<sup>2</sup> at a pressure of 24.5&#xa0;MPa and mass flux of 1,260&#xa0;kg/(m<sup>2</sup>&#xb7;s). Here, the numerical model is tested on a 1,500&#xa0;mm heated vertical 3-D tube with adiabatic sections (500&#xa0;mm long) at the inlet and outlet to verify that it accurately predicts HTCS of SCW. SCW flows from the bottom to the top with a given inlet mass flux and temperature.</p>
<p>
<xref ref-type="fig" rid="F3">Figure 3</xref> depicts the wall temperature, <italic>T</italic>
<sub>w</sub>, in the experiment by <xref ref-type="bibr" rid="B18">Yamagata et al. (1972)</xref> and the numerical simulations under various uniform wall heat fluxes <italic>versus</italic> the fluid bulk enthalpy, <italic>i</italic>
<sub>b</sub>. The local wall temperature results exhibited a maximum relative deviation of &#x3c;2%. It is proved that the model is suitable for heat transfer calculation with SCW.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Simulated results compared to experimental results.</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g003.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="results|discussion" id="s3">
<title>3 Results and discussion</title>
<sec id="s3-1">
<title>3.1 Effects of finite-thickness wall on HTCs</title>
<p>This section presents the axial and circumferential variations in the HTCs in straight and helical tubes under the effects of the finite-thickness wall.</p>
<sec id="s3-1-1">
<title>3.1.1 Effects of finite-thickness wall in straight tubes on HTCs</title>
<p>This section presents the variation in the HTCs in straight tubes with inclination of 7.16&#xb0; due to the effects of a finite-thickness wall under Con. III.</p>
<p>
<xref ref-type="fig" rid="F4">Figure 4</xref> shows the axial dimentionless heat flux (<italic>q</italic>/<italic>q</italic>
<sub>av</sub>, ratio of local heat flux to the average heat flux on the inner wall) of three regions along the straight tube with wall thickness of 2&#xa0;mm under Con. III. Unlike a uniformly heated wall, the heat capacity along three regions (top, side, and bottom) in the tube with a finite-thickness wall changes with non-uniform turbulent kinetic energy, which results in differences in the circumferential dimentionless heat flux. At the lowest density, the dimentionless heat flux along the top region are the minimum, followed by those on the side.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Axial dimentionless heat flux along the straight tube with 2 mm wall thickness (Con. III: <italic>G</italic>
<sub>in</sub> &#x3d; 1, 260 kg).</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g004.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F5">Figure 5</xref> shows the axial heat transfer coefficients <italic>h</italic> of SCW along straight tubes under Con. III. The values of <italic>h</italic> in tubes with and without the finite-thickness wall are compared. The heat transfer coefficients at the bottom are the highest, followed by those at the sides and top. Unlike the differences in the heat flux factor along the three regions in the straight tube with wall thickness of 2&#xa0;mm, the heat transfer coefficients in straight tubes appear less sensitive to the wall.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Axial heat transfer coefficients along straight tubes with various wall thicknesses (Con. III: <italic>G</italic>
<sub>in</sub> &#x3d; 1, 260 kg).</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g005.tif"/>
</fig>
</sec>
<sec id="s3-1-2">
<title>3.1.2 Effects of finite-thickness wall in helical tubes on HTCs</title>
<p>This section presents the variation in the HTCs in helical tubes with helix diameter of 105.5 mm, caused by the effects of finite-thickness wall.</p>
<p>
<xref ref-type="fig" rid="F6">Figure 6</xref> shows the axial dimentionless heat flux, <italic>q</italic>/<italic>q</italic>
<sub>av</sub>, of four regions along the helical tube with wall thickness of 2&#xa0;mm under Con. III. In the helical tube with zero-thickness wall, a uniform heat flux was imposed on the tube wall, so the effects of the wall have been ignored. While in the helical tube with wall thickness of 2&#xa0;mm given volumetric heat source, the heat flux along four regions is significantly different. With the outlet enthalpy much lower than <italic>i</italic>
<sub>pc</sub> (thermophysical properties change dramatically near <italic>i</italic>
<sub>pc</sub> in <xref ref-type="fig" rid="F1">Figure 1</xref>), the dimentionless heat flux of the four regions change gently along the helical tube. The values of <italic>q</italic>/<italic>q</italic>
<sub>av</sub> represent the heat transfer capacity. The average values of <italic>q</italic>/<italic>q</italic>
<sub>av</sub> at the top and bottom are the same (1.06), whereas those on the inside and outside are 0.80 and 1.09, respectively. It can be concluded that the heat transfer capacity on the outside is the best, followed by those at the top and bottom, whereas that on the inside is the worst. The relative difference in <italic>q</italic>/<italic>q</italic>
<sub>av</sub> between the top and bottom, caused by buoyancy, is lower than 1%, whereas that between the inside and bottom, caused by the centrifugal force, is 36% in <xref ref-type="fig" rid="F6">Figure 6</xref>. Under a large inlet mass flux, the relative value of the centrifugal force is much higher than that of buoyancy in the helical tube.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Axial dimentionless heat flux along the helical tube with 2 mm wall thickness (<italic>D</italic> &#x003D; 105.5 mm, Con. III: <italic>G</italic>
<sub>in</sub> &#x3d; 1, 260 kg).</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g006.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F7">Figure 7</xref> shows the axial wall temperatures, <italic>T</italic>
<sub>w</sub>, along the four regions in helical tubes with helix diameter of 105.5&#xa0;mm under Con. III; the values of <italic>T</italic>
<sub>w</sub> in tubes with and without finite-thickness walls are compared. In <xref ref-type="fig" rid="F7">Figure 7A</xref>, the difference in the wall temperature between the inside and outside is 6.2&#xa0;K in the tube without the finite-thickness wall. In <xref ref-type="fig" rid="F7">Figure 7B</xref>, the wall temperature difference between the inside and outside is 4.4&#xa0;K in the tube with wall thickness of 2&#xa0;mm. The finite-thickness wall contributes to the circumferential heat exchange, significantly reducing the wall temperature difference between the inner and outer regions by 29%.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Axial wall temperature along helical tubes with various wall thicknesses (<italic>D</italic> &#x3d; 105.5 mm, Con. III: <italic>G</italic>
<sub>in</sub> &#x003D; 1, 260 kg).</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g007.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F8">Figure 8</xref> shows the axial heat transfer coefficients, <italic>h</italic>, of SCW along four regions in the helical tube with 105.5&#xa0;mm helix diameter under Con. I and Con. II. When the fluid enthalpy is close to <italic>i</italic>
<sub>pc</sub>, the heat transfer capacity increases with increase in the fluid velocity. However, as shown in <xref ref-type="fig" rid="F8">Figure 8</xref>, (<italic>&#x3c1;</italic>&#xb7;<italic>C</italic>
<sub>p</sub>), which affects the ability to accommodate and carry heat, reaches the maximum before <italic>i</italic>
<sub>pc</sub>. Correspondingly, under Con. I and Con. II, the maximum heat transfer coefficients at the four regions appear before <italic>i</italic>
<sub>pc</sub>. In <xref ref-type="fig" rid="F8">Figure 8A</xref> without the finite-thickness wall, the difference in the average heat transfer coefficient between the outer and inner regions is 6.5&#xa0;kW/(m<sup>2</sup>&#xb7;K) (relative difference 131%) under Con. I and 11.7&#xa0;kW/(m<sup>2</sup>&#xb7;K) (relative difference 142%) under Con. II. In <xref ref-type="fig" rid="F8">Figure 8B</xref> with wall thickness of 2&#xa0;mm, the difference in the average heat transfer coefficient between the outer and inner regions is 6.2&#xa0;kW/(m<sup>2</sup>&#xb7;K) (relative difference 128%) under Con. I and 11.4&#xa0;kW/(m<sup>2</sup>&#xb7;K) (relative difference 137%) under Con. II. The finite-thickness wall contributes to reducing the difference in <italic>h</italic> between the outer and inner regions.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Axial heat transfer coefficients along helical tubes with various wall thicknesses (<italic>D</italic> &#x3d; 105.5&#xa0;mm). <bold>(A)</bold> No wall-Con. I; <bold>(B)</bold> 2 mm wall thickness-Con. I; <bold>(C)</bold> No wall-Con. II; <bold>(D)</bold> 2 mm wall thickness-Con. II.</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g008.tif"/>
</fig>
<p>According to the development of the circumferential maximum wall temperature along the helical tubes, six typical temperature regions were selected: (i) the maximum wall temperature is lower than <italic>T</italic>
<sub>pc</sub> (<italic>T</italic>
<sub>w, max</sub> &#x3c; <italic>T</italic>
<sub>pc</sub>), (ii) the maximum wall temperature is higher but the average wall temperature is lower than <italic>T</italic>
<sub>pc</sub> (<italic>T</italic>
<sub>w, av</sub> &#x3c; <italic>T</italic>
<sub>pc</sub> &#x3c; <italic>T</italic>
<sub>w, max</sub>), (iii) the average wall temperature is higher but the minimum wall temperature is lower than <italic>T</italic>
<sub>pc</sub> (<italic>T</italic>
<sub>w, min</sub> &#x3c; <italic>T</italic>
<sub>pc</sub> &#x3c; <italic>T</italic>
<sub>w, av</sub>), (iv) the minimum wall temperature is higher but the fluid temperature is lower than <italic>T</italic>
<sub>pc</sub> (<italic>T</italic>
<sub>b</sub> &#x3c; <italic>T</italic>
<sub>pc</sub> &#x3c; <italic>T</italic>
<sub>w, min</sub>), (v) the fluid temperature reaches <italic>T</italic>
<sub>pc</sub> (<italic>T</italic>
<sub>b</sub> &#x3d; <italic>T</italic>
<sub>pc</sub>), (vi) and the fluid temperature exceeds <italic>T</italic>
<sub>pc</sub> (<italic>T</italic>
<sub>pc</sub> &#x3c; <italic>T</italic>
<sub>b</sub>).</p>
<p>
<xref ref-type="fig" rid="F9">Figure 9</xref> shows the development of the dimensionless circumferential wall temperature (<italic>T</italic>
<sub>w</sub>/<italic>T</italic>
<sub>w,max</sub>) in the helical tube with 105.5&#xa0;mm helix diameter under Con. I. <xref ref-type="fig" rid="F9">Figure 9A</xref> depicts the case of the tube without the finite-thickness wall; <xref ref-type="fig" rid="F9">Figure 9B</xref> depicts the case of the tube with wall thickness of 2&#xa0;mm. The maximum wall temperature appears when <italic>&#x3b8;</italic> increases from 0&#xb0; to 360&#xb0;, with the preceding wall temperature changing drastically before reaching <italic>T</italic>
<sub>w,max</sub>. The non-uniformities in the wall temperature shown in <xref ref-type="fig" rid="F9">Figure 9</xref> result from the combined effects of buoyancy and centrifugal forces. <italic>&#x3b8;</italic> of <italic>T</italic>
<sub>w,max</sub> aids the comparison of the relative scale between the centrifugal force and buoyancy. In <xref ref-type="fig" rid="F9">Figure 9A</xref>, for the temperature regions from <italic>T</italic>
<sub>w, max</sub> &#x3c; <italic>T</italic>
<sub>pc</sub> to <italic>T</italic>
<sub>b</sub> &#x3c; <italic>T</italic>
<sub>pc</sub> &#x3c; <italic>T</italic>
<sub>w,min</sub>, the maximum wall temperature moves from the inside to the top (from 64&#xb0; to 54&#xb0;), which proves that the buoyancy effect is enhanced. For the temperature region of <italic>T</italic>
<sub>w,max</sub> &#x3d; <italic>T</italic>
<sub>pc</sub>, the maximum wall temperature moves from the top to the inside (from 54&#xb0; to 72&#xb0;), which proves that with increase in the fluid temperature, the density decreases whereas the fluid velocity increases, resulting in increased centrifugal force. For the temperature region of <italic>T</italic>
<sub>pc</sub> &#x3c; <italic>T</italic>
<sub>b</sub>, the maximum wall temperature appears at <italic>&#x3b8;</italic>
<sub>Tw, max</sub> &#x3d; 90&#xb0; (inner region). The circumferential (<italic>T</italic>
<sub>w</sub>/<italic>T</italic>
<sub>w,max</sub>) differences in helical tubes with and without a finite-thickness wall are depicted in <xref ref-type="fig" rid="F10">Figure 10</xref>, where the value is 37&#xa0;K in <xref ref-type="fig" rid="F9">Figure 9A</xref> and 13&#xa0;K in <xref ref-type="fig" rid="F9">Figure 9B</xref>. The finite-thickness wall contributes to circumferential heat diffusion, which homogenizes the circumferential wall temperature distribution. The intervals of <italic>&#x3b8;</italic>
<sub>Tw, max</sub> in helical tubes with and without the finite-thickness wall are 36&#xb0; in <xref ref-type="fig" rid="F9">Figure 9A</xref> and 27&#xb0; in <xref ref-type="fig" rid="F9">Figure 9B</xref>. The finite-thickness wall reduces the interval of <italic>&#x3b8;</italic>
<sub>Tw, max</sub>.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Development of circumferential <italic>T</italic>
<sub>w</sub>/<italic>T</italic>
<sub>w,max</sub> in the helical tube versus <italic>&#x3b8;</italic> (<italic>D</italic> &#x3d; 105.5 mm, Con. I: <italic>G</italic>
<sub>in</sub> &#x3d; 225 kg).</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g009.tif"/>
</fig>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Development of circumferential heat transfer coefficients in the helical tube (Con. I: <italic>G</italic>
<sub>in</sub> &#x3d; 225 kg).</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g010.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F10">Figure 10</xref> shows the development of the circumferential heat transfer coefficient <italic>h</italic> of SCW in the helical tube with 105.5&#xa0;mm helix diameter under Con. I. <xref ref-type="fig" rid="F10">Figure 10A</xref> depicts the results for the tube without a finite-thickness wall, and <xref ref-type="fig" rid="F10">Figure 10B</xref> depicts the results for the tube with wall thickness of 2&#xa0;mm. In <xref ref-type="fig" rid="F10">Figure 10A</xref>, for the temperature regions from <italic>T</italic>
<sub>w, max</sub> &#x3c; <italic>T</italic>
<sub>pc</sub> to <italic>T</italic>
<sub>b</sub> &#x3d; <italic>T</italic>
<sub>pc</sub>, circumferential <italic>h</italic>
<sub>max</sub> increases from 7.7&#xa0;kW/(m<sup>2</sup>&#xb7;K) to 21.7&#xa0;kW/(m<sup>2</sup>&#xb7;K). For the temperature region of <italic>T</italic>
<sub>pc</sub> &#x3d; <italic>T</italic>
<sub>b</sub>, the average circumferential <italic>h</italic> is the maximum (17.4&#xa0;kW/(m<sup>2</sup>&#xb7;K)), with the maximum relative difference of 216%. For the temperature region of <italic>T</italic>
<sub>pc</sub> &#x3c; <italic>T</italic>
<sub>b</sub>, circumferential <italic>h</italic> drops sharply in the range of 1.9&#x2013;5.9&#xa0;kW/(m<sup>2</sup>&#xb7;K), which deteriorates the heat transfer. The maximum circumferential <italic>h</italic> in helical tubes without and with finite-thickness walls are 14.8&#xa0;kW/(m<sup>2</sup>&#xb7;K) (<xref ref-type="fig" rid="F10">Figure 10A</xref>) and 11.1&#xa0;kW/(m<sup>2</sup>&#xb7;K) (<xref ref-type="fig" rid="F10">Figure 10B</xref>), respectively. The finite-thickness wall contributes to the circumferential heat diffusion, which homogenizes the circumferential <italic>h</italic> distribution.</p>
</sec>
</sec>
<sec id="s3-2">
<title>3.2 HTCs in helical tubes with finite-thickness wall</title>
<p>This section presents the axial and circumferential variations in HTCs in helical tubes with wall thickness of 2&#xa0;mm under the effects of buoyancy, mass flux, and helix diameter.</p>
<sec id="s3-2-1">
<title>3.2.1 Effect of buoyancy</title>
<p>As shown in <xref ref-type="fig" rid="F1">Figure 1</xref>, the thermophysical properties of SCW change dramatically with increase in temperature. These changes induce buoyancy in the cross section and axial thermal acceleration. To quantitatively analyze the effect of buoyancy and thermal acceleration on the HTCs in the tubes, Hall and Jackson (<xref ref-type="bibr" rid="B4">Hall and Jackson, 1969</xref>; <xref ref-type="bibr" rid="B5">Hall, 1971</xref>) proposed the buoyancy factor, Bo&#x2a;. <xref ref-type="bibr" rid="B12">McEligot et al. (1970)</xref>; (<xref ref-type="bibr" rid="B7">Jackson et al., 1989</xref>) proposed the acceleration factor, Kv.<disp-formula id="e14">
<mml:math id="m26">
<mml:mrow>
<mml:mtext>Bo</mml:mtext>
<mml:mo>&#x2a;</mml:mo>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mtext>Gr</mml:mtext>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi mathvariant="normal">R</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
</mml:mrow>
<mml:mn>3.425</mml:mn>
</mml:msup>
<mml:msup>
<mml:mi>Pr</mml:mi>
<mml:mn>0.8</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>where,<disp-formula id="e15">
<mml:math id="m27">
<mml:mrow>
<mml:mtext>Gr</mml:mtext>
<mml:mo>&#x2a;</mml:mo>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>&#x3b2;</mml:mi>
<mml:msup>
<mml:mi>d</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
<mml:mi>q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:msup>
<mml:mi>&#x3bd;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
<disp-formula id="e16">
<mml:math id="m28">
<mml:mrow>
<mml:mi mathvariant="normal">R</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mi>&#x3bc;</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>
<disp-formula id="e17">
<mml:math id="m29">
<mml:mrow>
<mml:mi>Pr</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3bc;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>where, Gr<sup>&#x2a;</sup> is the Grashof number, Re is the Reynolds number, Pr is the Prandtl number, <italic>g</italic> denotes the acceleration due to gravity, <italic>d</italic> denotes the inner diameter, <italic>q</italic> denotes the heat flux, <italic>&#x3b2;</italic> denotes the isothermal compression coefficient, <italic>&#x3bb;</italic> denotes the thermal conductivity, <italic>&#x3c5;</italic> denotes the kinematic viscosity, <italic>G</italic> denotes the mass flux, <italic>&#x3bc;</italic> denotes the dynamic viscosity, and <italic>C</italic>
<sub>p</sub> denotes the specific heat.</p>
<p>
<xref ref-type="bibr" rid="B7">Jackson et al. (1989)</xref>; (<xref ref-type="bibr" rid="B8">Koshizuka et al., 1995</xref>) proposed a theoretical model for the relationship between Bo&#x2a; and heat transfer coefficients under heating. When Bo&#x2a; &#x3c;6 &#xd7; 10<sup>&#x2212;7</sup>, the buoyancy effect on the HTCs is relatively weak; when 6 &#xd7; 10<sup>&#x2212;7</sup> &#x3c; Bo&#x2a; &#x3c;1.2 &#xd7; 10<sup>&#x2212;6</sup>, the buoyancy degrades the heat transfer; when 1.2 &#xd7; 10<sup>&#x2212;6</sup> &#x3c; Bo&#x2a; &#x3c;8 &#xd7; 10<sup>&#x2212;6</sup>, the buoyancy helps to recover the heat transfer; when Bo&#x2a; &#x3e;6 &#xd7; 10<sup>&#x2212;7</sup>, the buoyancy enhances the heat transfer.<disp-formula id="e18">
<mml:math id="m30">
<mml:mrow>
<mml:mtext>Kv</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mi mathvariant="normal">p</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi mathvariant="normal">R</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3bc;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(18)</label>
</disp-formula>where, <italic>&#x3b1;</italic>
<sub>p</sub> refers to volume expansivity.</p>
<p>
<xref ref-type="bibr" rid="B12">McEligot et al. (1970)</xref> proposed that when Kv &#x3c; 3 &#xd7; 10<sup>&#x2212;6</sup>, the thermal acceleration effect on the HTCs is relatively weak, when Kv &#x2265; 3 &#xd7; 10<sup>&#x2212;6</sup>, the thermal acceleration deteriorates the heat transfer.</p>
<p>
<xref ref-type="fig" rid="F11">Figure 11</xref> shows the maximum Bo&#x2a; and Kv calculated using Equation <xref ref-type="disp-formula" rid="e14">14</xref>&#x2212;(<xref ref-type="disp-formula" rid="e18">18</xref>) for the nine cases with various inlet mass flux and helix diameter. The effects of the helix diameter on the maximum Bo&#x2a; and Kv are relatively weak, whereas obvious differences in the maximum Bo&#x2a; and Kv are observed under various inlet mass fluxes. In <xref ref-type="fig" rid="F11">Figure 11A</xref>, for the range of 1.2 &#xd7; 10<sup>&#x2212;6</sup> &#x3c; Bo&#x2a;<sub>max</sub> &#x3c; 8 &#xd7; 10<sup>&#x2212;6</sup> under Con. I and Con. II, the buoyancy affects the heat transfer and helps in recovering it. Under Con. III, for Bo&#x2a; &#x3c;6 &#xd7; 10<sup>&#x2212;7</sup>, the buoyancy effect on the HTCs is relatively weak. In <xref ref-type="fig" rid="F11">Figure 11B</xref>, Kv &#x3c; 3 &#xd7; 10<sup>&#x2212;6</sup>, which indicates that the thermal acceleration effect on the HTCs is relatively weak.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Maximum Bo&#x2a;and Kv in helical tubes with various helix diameter under various conditions (<italic>&#x3b4;</italic> &#x3d; 2&#xa0;mm). <bold>(A)</bold> No wall; <bold>(B)</bold> 2 mm wall thickness.</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g011.tif"/>
</fig>
</sec>
<sec id="s3-2-2">
<title>3.2.2 Effect of mass flux</title>
<p>In <xref ref-type="sec" rid="s3-1-2">Section 3.1.2</xref>, it can be concluded that the finite-thickness wall contributes to the circumferential heat exchange, significantly reducing the wall temperature and <italic>h</italic> difference between the inner and outer regions. This section mainly focuses on analyzing the effects of mass flux on heat transfer in the helical tubes with wall thickness of 2&#xa0;mm.</p>
<p>
<xref ref-type="fig" rid="F12">Figure 12</xref> shows the axial heat transfer coefficients, <italic>h</italic>, of SCW along four regions in the helical tube with 105.5&#xa0;mm helix diameter and wall thickness of 2&#xa0;mm under Con. III. Comparing <xref ref-type="fig" rid="F12">Figure 12</xref> with <xref ref-type="fig" rid="F8">Figures 8B, D</xref>, the relative differences in the average heat transfer coefficient in the top-to-bottom direction are 5%, 3%, and much less than 1% under Con. I, Con. II, and Con. III, respectively, which indicates that the buoyancy effect is weakened with increasing mass flux. In the inner-to-outer direction, the relative differences in the average heat transfer coefficient were 128%, 137%, and 211% under Con. I, Con. II, and Con. III, respectively, which indicates that with increasing mass flux, the heat transfer capacity in the outer region is enhanced and the difference in the heat transfer capacity between the inner and outer regions widens.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Axial heat transfer coefficients along the helical tube (<italic>D</italic> &#x3d; 105.5 mm, <italic>&#x3b4;</italic> &#x3d; 2 mm, Con. III: <italic>G</italic>
<sub>in</sub> &#x003D; 1, 260 kg).</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g012.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F13">Figure 13</xref> shows the development of the dimensionless circumferential wall temperature (<italic>T</italic>
<sub>w</sub>/<italic>T</italic>
<sub>w,max</sub>) in the helical tube with 105.5&#xa0;mm helix diameter and wall thickness of 2&#xa0;mm under various inlet mass fluxes. The maximum differences in the wall temperature under Con. I, Con. II, and Con. III are 11&#xa0;K, 9&#xa0;K, and 4&#xa0;K, respectively. With increase in the mass flux, Re increases and the heat transfer capacity is enhanced, which reduces the difference in the wall temperature. In addition, <italic>&#x3b8;</italic>
<sub>Tw, max</sub> shifts from 70&#xb0; to 90&#xb0;, and the effect of buoyancy is weakened whereas that of the centrifugal force is strengthened.</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>Development of circumferential <italic>T</italic>
<sub>w</sub>/<italic>T</italic>
<sub>w,max</sub> in helical tubes <italic>versus &#x3b8;</italic> under various conditions (<italic>D</italic> &#x3d; 105.5&#xa0;mm, <italic>&#x3b4;</italic> &#x3d; 2&#xa0;mm).</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g013.tif"/>
</fig>
</sec>
<sec id="s3-2-3">
<title>3.2.3 Effect of helix diameter</title>
<p>
<xref ref-type="fig" rid="F14">Figure 14</xref> shows the development of the circumferential dimensionless wall temperature (<italic>T</italic>
<sub>w</sub>/<italic>T</italic>
<sub>w,max</sub>) in helical tubes with various helix diameters and wall thickness of 2&#xa0;mm under Con. I. <xref ref-type="fig" rid="F14">Figures 14A, B</xref> depict the results for the cases of helix diameters of 70.5&#xa0;mm and 140.5&#xa0;mm, respectively. With <italic>&#x3b8;</italic> increasing from 0&#xb0; to 360&#xb0;, <italic>T</italic>
<sub>w, max</sub> occurs between the top region (0&#xb0;) and inner region (90&#xb0;); the wall temperature before the maximum value varies more dramatically than that after the maximum value. Comparing <xref ref-type="fig" rid="F14">Figure 14</xref> with <xref ref-type="fig" rid="F9">Figure 9B</xref>, with the helix diameter increasing from 70.5&#xa0;mm to 140.5&#xa0;mm, <italic>&#x3b8;</italic>
<sub>Tw, max</sub> decreases from 79&#xb0; to 70&#xb0;&#x2013;66&#xb0; in the <italic>T</italic>
<sub>w, max</sub> &#x3c; T<sub>pc</sub> section, shifting from the inner region to the top region; in the section <italic>T</italic>
<sub>pc</sub> &#x3c; <italic>T</italic>
<sub>b</sub>, the maximum difference in <italic>T</italic>
<sub>w</sub>/<italic>T</italic>
<sub>w, max</sub> decreases from 0&#xb0; to 360&#xb0; (0.019, 0.018, 0.017), and <italic>T</italic>
<sub>w</sub> gradually moves to the high-temperature zone.</p>
<fig id="F14" position="float">
<label>FIGURE 14</label>
<caption>
<p>Development of circumferential <italic>T</italic>
<sub>w</sub>/<italic>T</italic>
<sub>w,max</sub> in helical tubes versus <italic>&#x3b8;</italic> (<italic>&#x3b4;</italic> &#x3d; 2 mm, Con. I: <italic>G</italic>
<sub>in</sub> &#x3d; 225 kg). <bold>(A)</bold> <italic>D</italic> &#x3d; 70.5 mm; <bold>(B)</bold> <italic>D</italic> &#x3d; 140.5 mm.</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g014.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F15">Figure 15</xref> shows the development of the circumferential heat flux factor <italic>q</italic>/<italic>q</italic>
<sub>av</sub> along the helical tube with wall thickness of 2&#xa0;mm under Con. I. <xref ref-type="fig" rid="F15">Figures 15A, B, C</xref> depict the results for helix diameters of 70.5 mm, 105.5&#xa0;mm, and 140.5&#xa0;mm, respectively. Under the effect of centrifugal force, the mass flux in the inner region is lower than that in the outer region, resulting in poor heat transfer capacity. In <xref ref-type="fig" rid="F15">Figure 15B</xref>, for the case of <italic>T</italic>
<sub>w, max</sub> &#x3c; <italic>T</italic>
<sub>pc</sub> to <italic>T</italic>
<sub>b</sub> &#x3c; <italic>T</italic>
<sub>pc</sub> &#x3c; <italic>T</italic>
<sub>w, min</sub>, the buoyancy drains the fluid flow close to the wall from bottom to top, which shifts <italic>&#x3b8;</italic>
<sub>Tw, max</sub> from 70&#xb0; to 54&#xb0;. With increasing fluid temperature, the fluid velocity increases and density decreases, which enhances the centrifugal force, resulting in lower heat transfer capacity in the inner region. In the case of <italic>T</italic>
<sub>pc</sub> &#x3c; <italic>T</italic>
<sub>b</sub>, (<italic>q</italic>/<italic>q</italic>
<sub>av</sub>)<sub>min</sub> appears at <italic>&#x3b8;</italic>
<sub>Tw, max</sub> &#x3d; 90&#xb0; (inner region). <italic>q</italic>/<italic>q</italic>
<sub>av</sub> is greatly affected by the helix diameter before the fluid temperature reaches <italic>T</italic>
<sub>pc</sub>. In the region of <italic>T</italic>
<sub>w, max</sub> &#x3c; <italic>T</italic>
<sub>pc</sub>, the maximum difference in <italic>q</italic>/<italic>q</italic>
<sub>av</sub> for helix diameters of 70.5&#xa0;mm, 105.5 mm, and 140.5&#xa0;mm are 0.40 (<xref ref-type="fig" rid="F15">Figure 15A</xref>), 0.34 (<xref ref-type="fig" rid="F15">Figure 15B</xref>), and 0.33 (<xref ref-type="fig" rid="F15">Figure 15C</xref>), respectively. With increase in the helix diameter, the difference in <italic>q</italic>/<italic>q</italic>
<sub>av</sub> decreases.</p>
<fig id="F15" position="float">
<label>FIGURE 15</label>
<caption>
<p>Development of circumferential dimentionless heat flux in helical tubes versus <italic>&#x3b8;</italic> (<italic>&#x3b4;</italic> &#x3d; 2 mm, Con. I: <italic>G</italic>
<sub>in</sub> &#x3d; 225 kg). <bold>(A)</bold> <italic>D</italic> &#x3d; 70.5 mm; <bold>(B)</bold> <italic>D</italic> &#x3d; 105.5 mm; <bold>(C)</bold> <italic>D</italic> &#x3d; 140.5 mm.</p>
</caption>
<graphic xlink:href="fenrg-11-1129469-g015.tif"/>
</fig>
</sec>
</sec>
</sec>
<sec sec-type="conclusion" id="s4">
<title>4 Conclusion</title>
<p>The axial and circumferential non-uniform HTCs of SCW in helical tubes with various helix diameters under 24.5&#xa0;MPa were investigated and the effects of finite-thickness wall, inlet mass flux, and helix diameter were analyzed. The numerical simulation results were first compared with the experimental data from Yamagata et al. (<xref ref-type="bibr" rid="B18">Yamagata et al., 1972</xref>) to validate our numerical model, which uses the RNG <italic>k&#x2013;&#x3b5;</italic> turbulent model coupled with an enhanced wall function. The axial distribution of fluid temperature, wall temperature, and heat transfer coefficients, maximum Bo&#x2a; and Kv and development of circumferential <italic>T</italic>
<sub>w</sub>/<italic>T</italic>
<sub>w,max</sub>, heat transfer coefficients, <italic>q</italic>/<italic>q</italic>
<sub>av</sub> were studied. Based on these results, the following conclusions can be drawn:<list list-type="simple">
<list-item>
<p>(1) The finite-thickness wall barely affects local heat transfer in the straight tube with inlet mass flux of 1,260&#xa0;kg/(m<sup>2</sup> s) whereas the finite-thickness wall contributes to the circumferential heat exchange, significantly reducing the wall temperature difference between the inner and outer regions by 29% and with relative difference in <italic>q</italic>/<italic>q</italic>
<sub>av</sub> between the inside and bottom by 36%.</p>
</list-item>
<list-item>
<p>(2) The effects of the helix diameter on the maximum Bo&#x2a; and Kv are relatively weak, whereas obvious differences in the maximum Bo&#x2a; and Kv are observed under various inlet mass fluxes.</p>
</list-item>
<list-item>
<p>(3) With increase in the mass flux, the difference in the wall temperature gets reduced; the heat transfer capacity in the outer region is enhanced and the difference in the heat transfer capacity between the inner and outer regions widens.</p>
</list-item>
<list-item>
<p>(4) With increase in helix diameter, from 70.5&#xa0;mm to 140.5 mm, <italic>&#x3b8;</italic>
<sub>Tw, max</sub> decreases from 79&#xb0; to 70&#xb0;&#x2013;66&#xb0; in the <italic>T</italic>
<sub>w, max</sub> &#x3c; T<sub>pc</sub> section, the reason of which is that the buoyancy force strengthens first and then with increasing flow velocity, the centrifugal force increases; in the section <italic>T</italic>
<sub>pc</sub> &#x3c; <italic>T</italic>
<sub>b</sub>, the maximum difference in <italic>T</italic>
<sub>w</sub>/<italic>T</italic>
<sub>w, max</sub> decreases from 0&#xb0; to 360&#xb0; (0.019, 0.018, 0.017), and <italic>T</italic>
<sub>w</sub> gradually moves to the high-temperature zone. With increase in the helix diameter, the difference in <italic>q</italic>/<italic>q</italic>
<sub>av</sub> decreases.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s5">
<title>Data availability statement</title>
<p>The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s6">
<title>Author contributions</title>
<p>ML: numerical simulation; application of statistical, mathematical, and computational to analyze the data. ZZ: review and editing; acquisition of the financial support. XY: ideas; evolution of overarching research goals and aims. SJ: review.</p>
</sec>
<sec id="s7">
<title>Funding</title>
<p>The Analyses and Tests for the Key Safety Technology of Nuclear Heating Reactor NHR200-II; Young Talent Project of China National Nuclear Corporation.</p>
</sec>
<ack>
<p>The authors are grateful for the support of the Analyses and Tests for the Key Safety Technology of Nuclear Heating Reactor NHR200-II and Young Talent Project of China National Nuclear Corporation.</p>
</ack>
<sec sec-type="COI-statement" id="s8">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Behnia</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Parneix</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Durbin</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>Prediction of heat transfer in an axisymmetric turbulent jet impinging on a flat plate</article-title>. <source>Int. J. Heat. Mass Transf.</source> <volume>41</volume> (<issue>12</issue>), <fpage>1845</fpage>&#x2013;<lpage>1855</lpage>. <pub-id pub-id-type="doi">10.1016/s0017-9310(97)00254-8</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Durbin</surname>
<given-names>P. A.</given-names>
</name>
</person-group> (<year>1991</year>). <article-title>Near-wall turbulence closure modeling without &#x201c;damping functions&#x201d;</article-title>. <source>Theor. Comp. Fluid Dyn.</source> <volume>3</volume>, <fpage>1</fpage>&#x2013;<lpage>13</lpage>. <pub-id pub-id-type="doi">10.1007/bf00271513</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Fewster</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Jackson</surname>
<given-names>J. D.</given-names>
</name>
</person-group> (<year>2004</year>). &#x201c;<article-title>Experiments on supercritical pressure convective heat transfer having relevance to SCWR</article-title>,&#x201d; in <conf-name>Proceedings of the International Congress on Advances in Nuclear Power Plants</conf-name>, <conf-loc>Pittsburgh, USA</conf-loc>, <conf-date>24 August 2022</conf-date> (<publisher-name>ICAPP &#x2019;4-4342</publisher-name>).</citation>
</ref>
<ref id="B4">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Hall</surname>
<given-names>W. B.</given-names>
</name>
<name>
<surname>Jackson</surname>
<given-names>J. D.</given-names>
</name>
</person-group> (<year>1969</year>). <source>Laminarisation of a turbulent pipe flow by buoyancy forces, 11th national heat transfer conf., special session on laminarization of turbulent flows</source>. <publisher-loc>Mimmeapolis</publisher-loc>: <publisher-name>ASME</publisher-name>. <comment>Paper No. 69-HT-55</comment>.</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hall</surname>
<given-names>W. B.</given-names>
</name>
</person-group> (<year>1971</year>). <article-title>Heat transfer near the critical point</article-title>. <source>Adv. Heat Transf.</source> <volume>7</volume>, <fpage>1</fpage>&#x2013;<lpage>86</lpage>.</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>X. R.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X. T.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>S. Y.</given-names>
</name>
<name>
<surname>Tu</surname>
<given-names>J. Y.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Numerical study on heat transfer characteristic of CO<sub>2</sub> in helical tube at supercritical pressure</article-title>. <source>A. Energy Sci. Technol.</source> <volume>52</volume> (<issue>05</issue>), <fpage>769</fpage>&#x2013;<lpage>775</lpage>. <pub-id pub-id-type="doi">10.7538/yzk.2017.youxian.0605</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jackson</surname>
<given-names>J. D.</given-names>
</name>
<name>
<surname>Cotton</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Axcell</surname>
<given-names>B. P.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>Studies of mixed convection in vertical tubes</article-title>. <source>Int. J. Heat Fluid Flow</source> <volume>10</volume> (<issue>01</issue>), <fpage>2</fpage>&#x2013;<lpage>15</lpage>. <pub-id pub-id-type="doi">10.1016/0142-727x(89)90049-0</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Koshizuka</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Takano</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Oka</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title>Numerical analysis of deterioration phenomena in heat transfer to supercritical water</article-title>. <source>Int. J. Heat. Mass Transf.</source> <volume>38</volume> (<issue>16</issue>), <fpage>3077</fpage>&#x2013;<lpage>3084</lpage>. <pub-id pub-id-type="doi">10.1016/0017-9310(95)00008-w</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lemmon</surname>
<given-names>E. W.</given-names>
</name>
<name>
<surname>Huber</surname>
<given-names>M. L.</given-names>
</name>
<name>
<surname>McLinden</surname>
<given-names>M. O.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Reference fluid thermodynamic and transport properties, NIST standard reference database 23</article-title>.</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>F. B.</given-names>
</name>
<name>
<surname>Bai</surname>
<given-names>B. F.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Flow and heat transfer of supercritical water in the vertical helically-coiled tube under half-side heating condition</article-title>. <source>Appl. Therm. Eng.</source> <volume>133</volume>, <fpage>512</fpage>&#x2013;<lpage>519</lpage>. <pub-id pub-id-type="doi">10.1016/j.applthermaleng.2018.01.047</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Z. H.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Y. X.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>J. F.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>D. L. H.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Heat transfer to supercritical water in circular tubes with circumferentially non-uniform heating</article-title>. <source>Appl. Therm. Eng.</source> <volume>70</volume> (<issue>01</issue>), <fpage>190</fpage>&#x2013;<lpage>200</lpage>. <pub-id pub-id-type="doi">10.1016/j.applthermaleng.2014.05.013</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McEligot</surname>
<given-names>D. M.</given-names>
</name>
<name>
<surname>Coon</surname>
<given-names>C. W.</given-names>
</name>
<name>
<surname>Perkins</surname>
<given-names>H. C.</given-names>
</name>
</person-group> (<year>1970</year>). <article-title>Relaminarization in tubes</article-title>. <source>Int. J. Heat. Mass Transf.</source> <volume>13</volume> (<issue>2</issue>), <fpage>431</fpage>&#x2013;<lpage>433</lpage>. <pub-id pub-id-type="doi">10.1016/0017-9310(70)90118-3</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mokry</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Pioro</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Kirillov</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Gospodinov</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Supercritical-water heat transfer in a vertical bare tube</article-title>. <source>Nucl. Eng. Des.</source> <volume>240</volume> (<issue>3</issue>), <fpage>568</fpage>&#x2013;<lpage>576</lpage>. <pub-id pub-id-type="doi">10.1016/j.nucengdes.2009.09.003</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rowinski</surname>
<given-names>M. K.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>J. Y.</given-names>
</name>
<name>
<surname>White</surname>
<given-names>T. J.</given-names>
</name>
<name>
<surname>Soh</surname>
<given-names>Y. C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Numerical investigation of supercritical water flow in a vertical pipe under axially non-uniform heat flux</article-title>. <source>Prog. Nucl. Energy</source> <volume>97</volume>, <fpage>11</fpage>&#x2013;<lpage>25</lpage>. <pub-id pub-id-type="doi">10.1016/j.pnucene.2016.12.009</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sheeba</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Abhijith</surname>
<given-names>C. M.</given-names>
</name>
<name>
<surname>Jose Prakash</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Experimental and numerical investigations on the heat transfer and flow characteristics of a helical coil heat exchanger</article-title>. <source>Int. J. Refrig.</source> <volume>99</volume>, <fpage>490</fpage>&#x2013;<lpage>497</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijrefrig.2018.12.002</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>K. Z.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>X. X.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Y. Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Dang</surname>
<given-names>C. B.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Numerical investigation on heat transfer of supercritical CO<sub>2</sub> in heated helically coiled tubes</article-title>. <source>J. Supercrit. Fluid</source> <volume>99</volume>, <fpage>112</fpage>&#x2013;<lpage>120</lpage>. <pub-id pub-id-type="doi">10.1016/j.supflu.2015.02.001</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>S. X.</given-names>
</name>
<name>
<surname>Niu</surname>
<given-names>Z. Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>J. L.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Experimental investigations on heat transfer of supercritical CO<sub>2</sub> flowing through a helically coiled tube</article-title>. <source>Nucl. Power Eng.</source> <volume>35</volume> (<issue>01</issue>), <fpage>28</fpage>&#x2013;<lpage>31</lpage>.</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yamagata</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Nishikawa</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Hasegawa</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Fujii</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Yoshida</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>1972</year>). <article-title>Forced convective heat transfer to supercritical water flowing in tubes</article-title>. <source>Int. J. Heat. Mass Transf.</source> <volume>15</volume> (<issue>12</issue>), <fpage>2575</fpage>&#x2013;<lpage>2593</lpage>. <pub-id pub-id-type="doi">10.1016/0017-9310(72)90148-2</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G. R.</given-names>
</name>
<name>
<surname>Liao</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J. D.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Numerical study of characteristic influence on heat transfer of supercritical CO<sub>2</sub> in helically coiled tube with non-circular cross section</article-title>. <source>Int. J. Heat. Mass Transf.</source> <volume>176</volume>, <fpage>121511</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijheatmasstransfer.2021.121511</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S. X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C. D.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>J. L.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Mixed convective heat transfer of CO<sub>2</sub> at supercritical pressures flowing upward through a vertical helically coiled tube</article-title>. <source>Appl. Therm. Eng.</source> <volume>88</volume>, <fpage>61</fpage>&#x2013;<lpage>70</lpage>. <pub-id pub-id-type="doi">10.1016/j.applthermaleng.2014.10.031</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X. T.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>P. X.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Study on heat transfer and flow instability of supercritical water in vertical tube</article-title>. <source>A. Energy Sci. Technol.</source> <volume>49</volume> (<issue>11</issue>), <fpage>2011</fpage>&#x2013;<lpage>2016</lpage>. <pub-id pub-id-type="doi">10.7538/yzk.2015.49.11.2011</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Ye</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X. T.</given-names>
</name>
<name>
<surname>Ju</surname>
<given-names>H. M.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>S. Y.</given-names>
</name>
<name>
<surname>Tu</surname>
<given-names>J. Y.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Supercritical steam generator design and thermal analysis based on HTR-PM</article-title>. <source>Ann. Nucl. Energy</source> <volume>132</volume>, <fpage>311</fpage>&#x2013;<lpage>321</lpage>. <pub-id pub-id-type="doi">10.1016/j.anucene.2019.04.049</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>C. R.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X. T.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>P. X.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>S. Y.</given-names>
</name>
<name>
<surname>Tu</surname>
<given-names>J. Y.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Influences of tube wall on the heat transfer and flow instability of various supercritical pressure fluids in a vertical tube</article-title>. <source>Appl. Therm. Eng.</source> <volume>147</volume>, <fpage>242</fpage>&#x2013;<lpage>250</lpage>. <pub-id pub-id-type="doi">10.1016/j.applthermaleng.2018.10.024</pub-id>
</citation>
</ref>
</ref-list>
<sec id="s10">
<title>Nomenclature</title>
<def-list>
<def-item>
<term id="G1-fenrg.2023.1129469">
<bold>3-D</bold>
</term>
<def>
<p>Three-dimensional</p>
</def>
</def-item>
<def-item>
<term id="G2-fenrg.2023.1129469">
<bold>
<italic>A</italic>
</bold>
</term>
<def>
<p>Area, m<sup>2</sup>
</p>
</def>
</def-item>
<def-item>
<term id="G3-fenrg.2023.1129469">
<bold>bot</bold>
</term>
<def>
<p>Bottom</p>
</def>
</def-item>
<def-item>
<term id="G4-fenrg.2023.1129469">
<bold>Bo&#x2a;</bold>
</term>
<def>
<p>Buoyancy factor</p>
</def>
</def-item>
<def-item>
<term id="G5-fenrg.2023.1129469">
<bold>Con.</bold>
</term>
<def>
<p>Condition</p>
</def>
</def-item>
<def-item>
<term id="G6-fenrg.2023.1129469">
<bold>
<italic>C</italic>
</bold>
<sub>
<bold>p</bold>
</sub>
</term>
<def>
<p>Specific heat, J/(kg&#x2022;K)</p>
</def>
</def-item>
<def-item>
<term id="G7-fenrg.2023.1129469">
<bold>
<italic>D</italic>
</bold>
</term>
<def>
<p>Helical diameter, m</p>
</def>
</def-item>
<def-item>
<term id="G8-fenrg.2023.1129469">
<bold>
<italic>d</italic>
</bold>
</term>
<def>
<p>Diameter, m</p>
</def>
</def-item>
<def-item>
<term id="G9-fenrg.2023.1129469">
<bold>
<italic>E</italic>
</bold>
</term>
<def>
<p>Total energy, J/kg</p>
</def>
</def-item>
<def-item>
<term id="G10-fenrg.2023.1129469">
<bold>
<italic>G</italic>
</bold>
</term>
<def>
<p>Mass flux, kg/(m<sup>2</sup>&#x2022;s)</p>
</def>
</def-item>
<def-item>
<term id="G11-fenrg.2023.1129469">
<bold>Gr</bold>
</term>
<def>
<p>Grashof number</p>
</def>
</def-item>
<def-item>
<term id="G12-fenrg.2023.1129469">
<bold>
<italic>g</italic>
</bold>
</term>
<def>
<p>Acceleration of gravity, 9.81&#xa0;m/s<sup>2</sup>
</p>
</def>
</def-item>
<def-item>
<term id="G13-fenrg.2023.1129469">
<bold>HTC</bold>
</term>
<def>
<p>Heat transfer characteristic</p>
</def>
</def-item>
<def-item>
<term id="G14-fenrg.2023.1129469">
<bold>HTR-PM</bold>
</term>
<def>
<p>High-temperature gas-cooled reactor pebble-bed module</p>
</def>
</def-item>
<def-item>
<term id="G15-fenrg.2023.1129469">
<bold>
<italic>h</italic>
</bold>
</term>
<def>
<p>Heat transfer coefficient, W/(m<sup>2</sup>&#x2022;K)</p>
</def>
</def-item>
<def-item>
<term id="G16-fenrg.2023.1129469">
<bold>
<italic>i</italic>
</bold>
</term>
<def>
<p>Enthalpy, J/kg</p>
</def>
</def-item>
<def-item>
<term id="G17-fenrg.2023.1129469">
<bold>Kv</bold>
</term>
<def>
<p>Acceleration factor</p>
</def>
</def-item>
<def-item>
<term id="G18-fenrg.2023.1129469">
<bold>
<italic>k-&#x3b5;</italic>
</bold>
</term>
<def>
<p>Kappa-epsilon</p>
</def>
</def-item>
<def-item>
<term id="G19-fenrg.2023.1129469">
<bold>
<italic>L</italic>
</bold>
</term>
<def>
<p>Tube full length, m</p>
</def>
</def-item>
<def-item>
<term id="G20-fenrg.2023.1129469">
<bold>
<italic>l</italic>
</bold>
</term>
<def>
<p>Length, m</p>
</def>
</def-item>
<def-item>
<term id="G21-fenrg.2023.1129469">
<bold>NIST</bold>
</term>
<def>
<p>National Institute of Standards and Technology</p>
</def>
</def-item>
<def-item>
<term id="G22-fenrg.2023.1129469">
<bold>Pr</bold>
</term>
<def>
<p>Prandtl number</p>
</def>
</def-item>
<def-item>
<term id="G23-fenrg.2023.1129469">
<bold>
<italic>p</italic>
</bold>
</term>
<def>
<p>Pressure, Pa</p>
</def>
</def-item>
<def-item>
<term id="G24-fenrg.2023.1129469">
<bold>
<italic>q</italic>
</bold>
</term>
<def>
<p>Heat flux, W/m<sup>2</sup>
</p>
</def>
</def-item>
<def-item>
<term id="G25-fenrg.2023.1129469">
<bold>
<italic>q</italic>/<italic>q</italic>
</bold>
<sub>
<bold>av</bold>
</sub>
</term>
<def>
<p>Dimentionless heat flux</p>
</def>
</def-item>
<def-item>
<term id="G26-fenrg.2023.1129469">
<bold>Re</bold>
</term>
<def>
<p>Reynolds number</p>
</def>
</def-item>
<def-item>
<term id="G27-fenrg.2023.1129469">
<bold>RNG</bold>
</term>
<def>
<p>Renormalization group</p>
</def>
</def-item>
<def-item>
<term id="G28-fenrg.2023.1129469">
<bold>SCW</bold>
</term>
<def>
<p>Supercritical-pressure water</p>
</def>
</def-item>
<def-item>
<term id="G29-fenrg.2023.1129469">
<bold>
<italic>T</italic>
</bold>
</term>
<def>
<p>Temperature, K</p>
</def>
</def-item>
<def-item>
<term id="G30-fenrg.2023.1129469">
<bold>
<italic>T</italic>
</bold>
<sub>
<bold>w</bold>
</sub>
<bold>/<italic>T</italic>
</bold>
<sub>
<bold>w, max</bold>
</sub>
</term>
<def>
<p>Dimensionless circumferential wall temperatures</p>
</def>
</def-item>
<def-item>
<term id="G31-fenrg.2023.1129469">
<bold>
<italic>u</italic>
</bold>
</term>
<def>
<p>Velocity, m/s</p>
</def>
</def-item>
<def-item>
<term id="G32-fenrg.2023.1129469">
<bold>
<italic>v</italic>
</bold>
</term>
<def>
<p>Kinematic viscosity, m<sup>2</sup>/s</p>
</def>
</def-item>
<def-item>
<term id="G33-fenrg.2023.1129469">
<bold>
<italic>x</italic>
</bold>
</term>
<def>
<p>Axial direction</p>
</def>
</def-item>
<def-item>
<term id="G34-fenrg.2023.1129469">
<bold>y</bold>
<sup>
<bold>&#x2b;</bold>
</sup>
</term>
<def>
<p>Non-dimensional wall distance</p>
</def>
</def-item>
</def-list>
<def-list>
<title>Greek symbols</title>
<def-item>
<term id="G35-fenrg.2023.1129469">
<bold>
<italic>&#x3b1;</italic>
</bold>
</term>
<def>
<p>Volume expansivity, K<sup>&#x2212;1</sup>
</p>
</def>
</def-item>
<def-item>
<term id="G36-fenrg.2023.1129469">
<bold>
<italic>&#x3b2;</italic>
</bold>
</term>
<def>
<p>Isothermal compression coefficient, MPa<sup>-1</sup>
</p>
</def>
</def-item>
<def-item>
<term id="G37-fenrg.2023.1129469">
<bold>
<italic>&#x3b4;</italic>
</bold>
</term>
<def>
<p>Wall thickness, mm</p>
</def>
</def-item>
<def-item>
<term id="G38-fenrg.2023.1129469">
<bold>
<italic>&#x3b5;</italic>
</bold>
</term>
<def>
<p>Turbulent dissipation rate, m<sup>2</sup>/s<sup>3</sup>
</p>
</def>
</def-item>
<def-item>
<term id="G39-fenrg.2023.1129469">
<bold>
<italic>&#x3b6;</italic>
</bold>
</term>
<def>
<p>Solution relative error between different grid</p>
</def>
</def-item>
<def-item>
<term id="G40-fenrg.2023.1129469">
<bold>
<italic>&#x3b8;</italic>
</bold>
</term>
<def>
<p>Circumferential angle, &#xb0;</p>
</def>
</def-item>
<def-item>
<term id="G41-fenrg.2023.1129469">
<bold>
<italic>&#x3bb;</italic>
</bold>
</term>
<def>
<p>Thermal conductivity, W/(m&#x2022;K)</p>
</def>
</def-item>
<def-item>
<term id="G42-fenrg.2023.1129469">
<bold>
<italic>&#x3bc;</italic>
</bold>
</term>
<def>
<p>Dynamic viscosity, Pa&#x2022;s</p>
</def>
</def-item>
<def-item>
<term id="G43-fenrg.2023.1129469">
<bold>
<italic>&#x3c1;</italic>
</bold>
</term>
<def>
<p>Density, kg/m<sup>3</sup>
</p>
</def>
</def-item>
<def-item>
<term id="G44-fenrg.2023.1129469">
<bold>
<italic>&#x3c3;</italic>
</bold>
<sub>
<bold>k</bold>
</sub>
</term>
<def>
<p>Turbulent Prandtl number of turbulent kinetic energy</p>
</def>
</def-item>
<def-item>
<term id="G45-fenrg.2023.1129469">
<bold>
<italic>&#x3c3;</italic>
</bold>
<sub>
<bold>&#x3b5;</bold>
</sub>
</term>
<def>
<p>Turbulent Prandtl number of turbulent dissipation rate</p>
</def>
</def-item>
<def-item>
<term id="G46-fenrg.2023.1129469">
<bold>&#x3c4;</bold>
</term>
<def>
<p>Stress tensor, N/m<sup>2</sup>
</p>
</def>
</def-item>
</def-list>
<def-list>
<title>Subscripts</title>
<def-item>
<term id="G47-fenrg.2023.1129469">
<bold>4</bold>
</term>
<def>
<p>Grid 4</p>
</def>
</def-item>
<def-item>
<term id="G48-fenrg.2023.1129469">
<bold>av</bold>
</term>
<def>
<p>Area average</p>
</def>
</def-item>
<def-item>
<term id="G49-fenrg.2023.1129469">
<bold>b</bold>
</term>
<def>
<p>Bulk fluid</p>
</def>
</def-item>
<def-item>
<term id="G50-fenrg.2023.1129469">
<bold>i, j</bold>
</term>
<def>
<p>Coordinate direction</p>
</def>
</def-item>
<def-item>
<term id="G51-fenrg.2023.1129469">
<bold>in</bold>
</term>
<def>
<p>Inlet</p>
</def>
</def-item>
<def-item>
<term id="G52-fenrg.2023.1129469">
<bold>max</bold>
</term>
<def>
<p>Maximum</p>
</def>
</def-item>
<def-item>
<term id="G53-fenrg.2023.1129469">
<bold>min</bold>
</term>
<def>
<p>Minimum</p>
</def>
</def-item>
<def-item>
<term id="G54-fenrg.2023.1129469">
<bold>o</bold>
</term>
<def>
<p>Outlet</p>
</def>
</def-item>
<def-item>
<term id="G55-fenrg.2023.1129469">
<bold>pc</bold>
</term>
<def>
<p>Pseudo-critical</p>
</def>
</def-item>
<def-item>
<term id="G56-fenrg.2023.1129469">
<bold>T</bold>
</term>
<def>
<p>Turbulent</p>
</def>
</def-item>
<def-item>
<term id="G57-fenrg.2023.1129469">
<bold>w</bold>
</term>
<def>
<p>Wall</p>
</def>
</def-item>
</def-list>
</sec>
</back>
</article>