Radial Temperature Distribution Research Articles

Methane steam reforming with axial variable diameter particle structures in grille-sphere composite packed bed: A numerical study of hydrogen production performance

Methane steam reforming with packed bed is an important approach to get hydrogen in industrial production. Compared with the conventional packed bed, the grille-sphere composite packed bed has the advantages of lowering the pressure drop, increasing the uniformity of the radial temperature distribution and improving overall efficiency. In the present paper, the axial variable diameter particle structure of grille-sphere composite packed bed has been proposed, and four different structures are studied to compare their performances. Moreover, two simulation methods are used in study: (a) the solid particle method which can show details of the simulation, (b) the equivalent medium method which can show the macro situation of simulation. Through these two methods, performances of four structures in different parameters are compared and analyzed. It is found that the axial variable diameter particle structures have better performance than the typical structures in different parameters cases. Inlet temperature, inlet velocity and inlet steam-carbon ratio are set to parameters respectively in simulation. According to the simulations, higher inlet temperature and inlet steam-carbon ratio will facilitate the reactions, however, higher inlet velocity will hinder the reactions. Through the solid particle method, the reasons of better performance have been found, and two methods difference are also studied. Compared with typical structures, the specific axial variable diameter structure brings 4.5% higher outlet hydrogen mass flow and 5.2% higher outlet hydrogen selectivity. Analysis of this study can guide the design of packed beds in industrial production and the results of this study have significance in industrial hydrogen production.

Effect of Radiation on Heat Transfer Inside Aeroengine Compressor Rotors

AbstractThe blade clearance in aero-engine compressors is mainly controlled by the radial growth of the compressor discs, to which the blades are attached. This growth depends on the radial distribution of the disc temperature, which in turn is determined by the heat transfer inside the internal rotating cavity between adjacent discs. The buoyancy-induced convection inside the cavity is significantly weaker than that associated with the forced convection in the external mainstream flow, and consequently radiation between the cavity surfaces cannot be ignored in the calculation of the disc temperatures. In this paper, both the Monte Carlo Ray-Trace (MCRT) method and the view factor (VF) method are used to calculate the radiative flux when the temperatures of the discs, shroud, and inner shaft of the compressor vary radially and axially. The Monte Carlo Ray-Trace method is computationally expensive, but it is able to incorporate the effect of complex geometries on radiation. The view factor method is quick to compute and, although the derivation becomes complicated when geometrical details are considered, it can be used as a first check of the effect of radiation in compressor cavities. Given distributions of surface temperatures, the blackbody and gray body heat fluxes were calculated for the discs, shroud, and inner shaft in two experimental compressor rigs and in a simulated compressor stage. For the experimental rigs, although the effect of radiation was relatively small for the case of large Grashof numbers, the relative effect of radiation increases as Gr (and consequently the convective heat transfer) decreases. For the simulated compressor, with a pressure ratio of 50:1 for state-of-the-art aircraft engines, radiation could have a significant effect on the disc temperature and consequently on the blade clearance; the effect is predicted to be more prominent for the next generation of aircraft engines with pressure ratios up to 70:1.

Design and Testing of a Rig to Investigate Buoyancy-Induced Heat Transfer in Aero-Engine Compressor Rotors

Abstract The change in compressor blade-tip clearance across the flight cycle depends on the expansion of the rotor, which in turn depends on the temperature and stress in the disks. The radial distribution of temperature is directly coupled to the buoyancy-driven flow and heat transfer in the rotating disk cavities. This paper describes a new test rig specifically designed to investigate this conjugate phenomenon. The rig test section includes four rotating disks enclosing three cavities. Two disks in the central cavity are instrumented with thermocouples to provide the radial distribution of temperature; the two outer cavities are thermally insulated to create appropriate boundary conditions for the heat transfer analysis. An axial throughflow of air is supplied between a stationary shaft and the bore of the disks. The temperature of the throughflow air is measured by thermocouples in rakes upstream and downstream of the central cavity. For a cold throughflow, the outer shroud of the central cavity is heated. Two independently controlled radiant heaters allow differential shroud temperatures for the upstream and downstream disks, as found in aero-engine compressors. Alternatively, the throughflow can be heated above the shroud temperature to simulate the transient conditions during engine operation where stratified flow can occur inside the cavity. The rig is designed to operate in conditions where both convective and radiative heat transfer dominate; all internal surfaces of the cavity are painted matt black to allow the accurate calculation of the radiant heat transfer. Separate attachments can be fitted to the cobs of both central disks; the attachments reduce the axial gap between the cobs—reducing the gap to zero creates a closed cavity, which can occur in some compressor designs. Other instrumentation includes heat-flux gages on the shroud and high-frequency pressure transducers embedded into the disk diaphragm to capture unsteady flow structures. Attention has been given to experimental uncertainty, including the computation of the thermal-disturbance errors, caused by thermocouples embedded in the rotating disks; a Bayesian statistical model is used to reduce the effect of uncertainties in temperature measurements on the calculation of the Nusselt number. The effect of relevant nondimensional parameters on the radial distribution of the disk and throughflow temperatures has been shown for some typical cases.

Measurement and Analysis of Buoyancy-Induced Heat Transfer in Aero-Engine Compressor Rotors

Abstract The flow inside cavities between corotating compressor disks of aero-engines is driven by buoyancy, with Grashof numbers exceeding 1013. This phenomenon creates a conjugate problem: the Nusselt numbers depend on the radial temperature distribution of the disks, and the disk temperatures depend on the Nusselt numbers. Furthermore, Coriolis forces in the rotating fluid generate cyclonic and anticyclonic circulations inside the cavity. Such flows are three-dimensional, unsteady, and unstable, and it is a challenge to compute and measure the heat transfer from the disks to the axial throughflow in the compressor. In this paper, Nusselt numbers are experimentally determined from measurements of steady-state temperatures on the surfaces of both disks in a rotating cavity of the Bath compressor-cavity rig. The data are collected over a range of engine-representative parameters and are the first results from a new experimental facility specifically designed to investigate buoyancy-induced flow. The radial distributions of disk temperature were collected under carefully controlled thermal boundary conditions appropriate for analysis using a Bayesian model combined with the equations for a circular fin. The Owen-Tang buoyancy model has been used to compare predicted radial distributions of disk temperatures and Nusselt numbers with some of the experimentally determined values, taking account of radiation between the interior surfaces of the cavity. The experiments show that the average Nusselt numbers on the disk increase as the buoyancy forces increase. At high rotational speeds, the temperature rise in the core, created by compressibility effects in the air, attenuates the heat transfer, and there is a critical rotational Reynolds number for which the Nusselt number is a maximum. In the cavity, there is an inner region dominated by forced convection and an outer region dominated by buoyancy-induced flow. The inner region is a mixing region, in which entrained cold throughflow encounters hot flow from the Ekman layers on the disks. Consequently, the Nusselt numbers on the downstream disk in the inner region tend to be higher than those on the upstream disk.

Particle Emissions and Disc Temperature Profiles from a Commercial Brake System Tested on a Dynamometer under Real-World Cycles

The particle emissions from a commercial brake system utilizing copper-free pads have been characterized on a brake dynamometer under two real-world driving cycles. These included a novel cycle developed from analysis of the database of the World Harmonized Test Procedure (WLTP-Brake) and a short version of the Los Angeles City Traffic cycle (3h-LACT) developed in the framework of the European LowBraSys project. Disc temperature measurements using an array of embedded thermocouples revealed a large temporal and spatial non-uniformity with the radial temperature distribution depending also on the test procedure. Averaging over the duration of the cycle, it effectively reduced the influence of thermocouple positioning, allowing for more reliable quantification of the effectiveness of convective cooling. Particulate Matter (PM) emissions were similar for both cycles with PM2.5 averaging at 2.2 (±0.2) mg/km over the WLTP-Brake and 2.2 (±0.2) mg/km over the 3h-LACT, respectively. The corresponding PM10 emissions were 5.6 (±0.2) mg/km and 8.6 (±0.7) mg/km, respectively. The measurements revealed the formation of nanosized particles peaking at 10 nm, which were thermally stable at 350 °C under both cycles. Volatile nanoparticles were observed over the more demanding 3h-LACT cycle, with their emission rates decreasing with increasing the tunnel flow, suggesting nucleation of organic vapors released during braking as a potential formation process.

Application of the triple Langmuir probe to study local plasma parameters in the radiofrequency discharge that is placed in an additional applied magnetostatic field

This paper considers the application of the triple Langmuir probe method to investigate the local plasma parameters of the radio frequency (RF) discharge that is placed in the applied magnetostatic field. The main criteria applicability for such a method for experimental diagnostics of RF plasma parameters is presented. As an example, this paper presents experimental results that were obtained using the proposed method on local plasma parameters in one of the cross sections of the discharge chamber of the RF ion source with an applied magnetostatic field. The influence of the magnetostatic field on the radial distributions of charged particle density and electron temperature is shown.

COMPOSITE Cu-Cr MATERIALS UNDER THERMAL ACTION OF ELECTRIC ARC DISCHARGE PLASMA

The results of optical emission spectroscopy (OES) investigation of plasma of electric arc discharges in steadystate mode between Cu-Cr composite electrodes, manufactured at different sintered temperatures: 750, 850, 950 or 1050 °C, is presented. In particular, the impact of sintering temperature on erosion resistanceof such composite materials, which was determined in indirect manner by estimation of metal vapours content in the midsection of discharge gaps, is studied by the analysis of plasma parameters. These contents were calculated in assumption of local thermodynamic equilibrium (LTE) on the base of experimentally obtained radial distributions of plasma temperature and electron density.

SPECTROSCOPY PECULIARITIES OF ELECTRIC ARC DISCHARGE PLASMA WITH IRON VAPOURS

This work is devoted to spectroscopy peculiarities of electric arc discharge plasma with iron vapours. The solution of the main issue of optical emission spectroscopy, namely, selection of iron spectral lines, to study the parameters of non-uniform and non-steady-state plasma source, was considered within this paper. Specifically, the Boltzmann plots technique was used for detailed analysing of application possibility of Fe I spectral lines as well as for determination of plasma temperature. The spatial profiles of selected spectral line emission intensities were used to measure the radial distributions of plasma temperature of free-burning arc discharge between consumable electrodes at 3.5 A.

Electrochemical-distributed thermal coupled model-based state of charge estimation for cylindrical lithium-ion batteries

State of Charge (SoC) is essential in a smart Battery Management System (BMS). Traditional SoC estimation methods consider the lithium battery cell as an isothermal body simplistically, which is not accurate. Due to the heat conduction and convection, the temperature distribution along the radius is nonuniform, which means the battery model parameters subject to temperature are not identical radially. This paper proposed a modified electrochemical-distributed thermal coupled model (MEDTM) for improving the accuracy of SoC estimation. The distributed thermal model is employed to design a robust H-infinity temperature observer for the estimation of the radial temperature distribution of battery cell, where a radial discretized method is used for the realization of the observer design. Based on the observed temperature, a SoC observer is developed with the backstepping method. Finally, experiments and simulation are implemented. The two constant discharging experiments indicate that MEDTM is more accurate than the single particle model with surface temperature (SPMST), and it can generate reliable temperature prediction. A comparative simulation and a constant discharging experiment verify the H-infinity temperature observer can obtain a more accurate estimation of the battery cell’s temperature than the luenberger observer, and the H-infinity temperature observer has an estimation error of 0.3 K below. Then, through a constant discharging experiment and a simulation with UDDS current, the proposed SoC observer is verified to accurately estimate the SoC by comparing with the true SoC, where the estimated error of SoC under the two cases is below 1.5% and below 0.4%, respectively. Therefore, the proposed SoC observer also has good performance.

Shock Propagation Following an Intense Explosion: Comparison Between Hydrodynamics and Simulations

The solution for the radial distribution of pressure, density, temperature and flow velocity fields in a blast wave propagating through a medium at rest, following an intense explosion, starting from hydrodynamic equations, is one of the classic problems in gas dynamics. However, there is very little direct verification of the theory and its assumptions from simulations of microscopic models. In this paper, we compare the results and assumptions of the hydrodynamic theory with results from large scale event driven molecular dynamics simulations of a hard sphere gas in three dimensions. We find that the predictions for the radial distribution of the thermodynamic quantities do not match well with the numerical data. We improve the theory by replacing the ideal gas law with a more realistic virial equation of state for the hard sphere gas. While this improves the theoretical predictions, we show that they still fail to describe the data well. To understand the reasons for this discrepancy, the different assumptions of the hydrodynamic theory are tested within the simulations. A key assumption of the theory is the existence of a local equation of state. We validate this assumption by showing that the local pressure, temperature and density obey the equation of state for a hard sphere gas. However, the probability distribution of the velocity fluctuations has non-gaussian tails, especially away from the shock front, showing that the assumption of local equilibrium is violated. This, along with neglect of heat conduction, could be the possible reasons for the mismatch between theory and simulations.

Spatiotemporal evolution of laser-induced plasmas in air: Influence of pressure

In this paper, we focused on the influence of gas pressure on the dynamic characteristics of laser-induced plasma (LIP) in air. The energy-dependent transmittances of a 1064 nm laser to LIPs generated at the pressure of 0.2–5 atm were measured. Laser Rayleigh scattering and Thomson scattering methods were applied to investigate the evolution of laser-produced shock wave and radial distribution of electron number density (ne) and electron temperature (Te) of LIP in air. The fraction of the shock wave energy to the total laser energy absorbed is about 40%–52% and increases upon increasing pressure. The higher the gas pressure, the greater ne and Te within the central of LIP and a faster decaying rate of Te and ne. The electron number density decays exponentially with time, and the decay index is between −0.891 and −1.177. A torus structure in LIPs is generated during the later stage of plasma decay and significantly affects the distribution of Te, ne and the intensity of plasma emission.

Magnetic properties of an ACSR conductor steel core at temperatures up to 230 ∘C and their impact on the transformer effect

The current density distribution between layers of steel-reinforced conductors having an odd number of aluminum layers is modified by the transformer effect, as a result of the steel core magnetisation. This redistribution of current affects the conductor alternating current resistance, the transitory temperature radial distribution, and depends on the total current intensity, and the core permeability, which varies with temperature. Although some steel-cored conductors are designed to operate at high temperatures to maximise their current capacity, the core magnetic properties in this scenario need to be further investigated. This paper presents the results of experimental work on an ACSR conductor, and on its steel core to investigate the steel core magnetic properties over an extended temperature range, and how they affect the transformer effect. Results show an unexpected hump-shaped variation of the permeability with temperature in a specific range of magnetic field strength, with maximum permeability values occurring at a temperature ranging from 150 to 170 ∘ C , which were not covered in previous studies. A similar behaviour is seen in the variation of the axial magnetic flux with increasing conductor temperature, causing the flux to decrease at high temperatures, and thus weakening the transformer effect during the conductor operation.

Ethanol supplement increases soot yields in nitrogen-diluted laminar ethylene diffusion flames at pressures from 3 to 5 bar

In spite of widely use of ethanol, mostly as a fuel extender in ground transportation engines, its sooting propensity with pressure and with ethanol content in the base fuel has not been clarified. Although, information on ethanol’s sooting and other combustion characteristics at atmospheric conditions is extensive, scaling this information to elevated pressures is problematic. Information that could be gathered from tractable laminar diffusion flames on the soot formation of ethanol’s response to pressure would be beneficial in assessing the soot emissions from engines fuelled with hydrocarbons supplemented with ethanol. To address this lack of information at pressures above atmospheric, high pressure soot formation in laminar co-flow diffusion flames fuelled by ethanol was investigated using nitrogen-diluted ethylene as the base fuel on a burner of 3 mm diameter fuel nozzle. Base fuel nitrogen to ethylene mass ratio was kept fixed at 6 for all experimental cases. In terms of total carbon mass flow, which was kept constant at 0.41 mg/s, ethanol’s contribution was varied as 0%, 10%, 30% and 40%, and the experiments were conducted at 3, 4, and 5 bar pressures. The main aim of the investigation was to evaluate the soot production at elevated pressures with increasing ethanol fraction in the fuel stream treating nitrogen-diluted ethylene as the other fuel component. Flame heights defined by the luminous visible tips did not display any significant changes as the ethanol fraction or pressure was changed. Spectrally-resolved line-of-sight soot radiation collected along the radial distance from the flame centreline at various heights above the burner exit with 1 mm increments was converted to radial soot temperature and volume fraction distributions through an Abel-type inversion algorithm. Within the bounds of pressure considered, soot production increased with the amount of ethanol added to the fuel stream. Largest incremental increase in soot production over ethylene diluted with nitrogen was observed at the condition with 10% carbon from ethanol in the fuel blend; further increases in ethanol fractions kept increasing the soot but at relatively lower rates. This non-linear effect with ethanol addition was discussed and it was argued that one of the reasons for this behavior is the influence of methyl radical produced by radical-induced decomposition of ethanol enhancing soot production for 10% ethanol case. With further increases in ethanol amount, this effect seems to slow down because of reduction in ethylene fraction and the increase in hydrogen to carbon ratio in the fuel stream at fixed carbon flow rate.

Evolution characteristic of axial magnetic field and Nernst effect in magnetized liner inertial fusion

Axial magnetic field is one of the main parameters of magnetized liner inertial fusion (MagLIF), which is greatly different from other traditional inertial confinement fusion configurations. The introduce of axial magnetic field dramatically increases energy deposit efficiency of alpha particles, when initial <i>B</i><i><sub>z</sub></i> increases from 0 to 30 T, the ratio of deposited alpha energy rises from 7% to 53%. In the MagLIF process, the evolvement of magnetic flux in fuel can be roughly divided into three main stages: undisturbed, oscillation, and equilibrium. The distributions and evolution characteristic of axial magnetic field are both determined by the liner conductivity, fuel conductivity, and the fluid dynamics. The pressure imbalance between fuel and liner, caused by laser injection, is the source of fluid oscillation, which is an intrinsic disadvantage of laser preheating method. This fluid oscillation does not lead the magnetic flux to decrease monotonically in the fuel during implosion process, but oscillate repeatedly, even increase in a short time. Nernst effect plays a negative role in MagLIF process. As initial axial magnetic field decreases from 30 to 20 to 10 T, the Nernst effect causes magnetic flux loss to increase from 28% to 44% to 73% correspondingly, and the deposited alpha energy ratio drops from 44% to 27% to 4% respectively. So the initial magnetic field is supposed to be moderately high. The radial distribution of temperature in fuel should be as uniform as possible after preheating, which is helpful in reducing the influence of Nernst effect. Compared with Nernst effect, the end loss effect is much responsible for rapid drawdown of fusion yield. A large number of physical images are acquired and summarized through this work, which are helpful in understanding the process of magnetic flux compression and diffusion in MagLIF process. The simulation can act as a powerful tool and the simulation results can serve as a useful guidance for the future experimental designs.

Control-oriented modeling of diesel oxidation catalyst

Diesel oxidation catalyst outlet temperature control is crucial for heat management to realize diesel particulate filter active regenerative control. In order to control the temperature of the active regeneration process in the filter, the temperature response process of the semi-physical oxidation catalyst model structure is proposed as a multi-stage inertia plus delay, and the equivalent inlet temperature step of the fuel oxidation reaction of the exhaust pipe. Combined with the test test, the control oriented oxidation catalyst model is established.A control-oriented oxidation catalyst model was constructed. By analysed the oxidation catalyst working process, the main chemical reactions, heat and mass transfer processes occurring inside the carrier were analyzed. Three-dimensional CFD model and one-dimensional chemical reaction kinetics model were established respectively. The radial and axial temperature distribution of the carrier was analyzed by model simulation. Based on the analysis of the system characteristics, the multi-step inertia plus delay semi-physical model structure was proposed. Combined with the test, the control oriented oxidation catalyst model is established. Select the appropriate working conditions to identify and verify the model parameters. The results show that the third order model can well indicate the temperature response characteristics of the oxidation catalyst outlet temperature. Considering the complexity of the system, the first-order and third-order model are selected as the basis of the control system design.

Combustion performance with ap and boron addition in staged hybrid rocket

This study investigated the effect of ammonium perchlorate (AP) and boron addition on controllability of the combustion temperature in fuel-rich effluent and system specific impulse. In addition, a series of combustion tests was carried out to analyse and evaluate the effect of boron addition on fuel regression rate and radial temperature profiles. To keep the hybrid rocket engine advantages, upper limit of AP and boron contents in solid fuel was set to be 15 and 10wt% respectively. The results suggested that AP addition proved to be very effective in controlling the effluent temperature. Not only that, it also offers less injected oxidizer mass flow, and resulted in less O/F variation, and less combustion pressure while producing the same combustion in temperature fuel-rich effluent. With boron addition, it also provided more uniform radial temperature distribution. Also, boron addition increased specific impulse by 6.8-13.7%. Secondary combustion using different solid fuels showed a decrease in c* efficiency with increasing boron content in solid fuel. Combustor length and combustion pressure were introduced to improve the combustion of boron particles, however, combustion efficiency was decreased with increasing combustion pressure and combustor length for given test conditions.

Solution of the gradient thermoelasticity problem for a cylinder with a heat-protected coating

The investigation of the stress-strain state of an infinitely long thermoelastic cylinder is carried out taking into account the scale effects. A thermal protective coating is applied to the outer side surface of the cylinder, the thermomechanical characteristics of which are functions of the radial coordinate. Thermal boundary conditions of the first kind are set on the stress-free side surfaces of the cylinder. To take into account the scale effects, the one-parameter gradient theory of thermoelasticity proposed by Aifantis is used. Additional boundary conditions and conjugation conditions for couple stresses are specified. The displacements and stresses are represented as the sum of the solutions of the thermoelasticity problem in the classical formulation and the gradient parts. After finding the radial temperature distribution, the thermoelasticity problem in the classical formulation with respect to radial displacements and stresses is solved numerically by the shooting method. Additional boundary layer terms for radial displacements at small values of the gradient parameter are found using the asymptotic method for solving linear differential equations with spatially varying coefficients (the Wentzel-Kramers-Brillouin method - WKB method). Calculations of radial displacements, Cauchy stresses, couple stresses and total stresses in the case of both homogeneous and inhomogeneous coatings are carried out using specific examples. It has been found that the Cauchy stresses and total stresses experience a jump at the boundary between the cylinder and the coating. The couple stresses at small gradient parameters are much less than the total stresses. An increase of the scale parameter reduces the values of radial displacements and total stresses. Deformations are continuous in both cases of coating (homogeneous and inhomogeneous). A comparative study of the influence of the value of the inhomogeneity parameter on the distribution of displacements and total stresses is carried out.

Evolution of the Quiescent Disk Surrounding a Superoutburst of the Dwarf Nova TW Virginis

Abstract In this paper, we investigate portions of the Kepler K2 Short Cadence light curve of the dwarf nova (DN) TW Vir at quiescence, using light-curve modeling. The light curve was separated into 24 sections, each with a data length of ∼0.93 days, comprising 4 sections before, and 20 after a superoutburst (SO). Due to morphological differences, the quiescent orbital modulation is classified into three types. Using a fixed disk radius and the two component stellar parameters, all 24 synthetic disk models from these sections show a consistent configuration, consisting of a disk and two hotspots: one at the vertical side of the edge of the disk and the other on the surface of the disk. Before the SO, the disk and a ringlike surface-hotspot are suddenly enhanced, triggering a precursor, and then the SO. At the end of the quiescent period following the SO and before the first normal outburst, the edge-hotspot becomes hotter, while the surface-hotspot switches into a “coolspot” with a coverage of nearly half of the disk’s surface. During quiescence, the surface-hotspot is always located at the outer part of the disk, with a constant radial width. A flat radial temperature distribution of the disk is found, and appears flatter when approaching the outburst. Like many U Gem-type DN with orbital periods of 3–5 hr, the mass transfer rate is significantly lower than the predictions of the standard/revised models of CV evolution.

Condensation heat transfer for downward flows of superheated steam-air mixture in a circular pipe

We measured the radial and axial temperature distributions of a saturated steam-air mixture in a vertical circular pipe (diameter, 49.5 mm; cooling height, 610 mm) in our previous study. Those of a superheated steam-air mixture in the same experimental setup were measured in this study. The condensation heat flux qc and heat transfer coefficient hc for the superheated mixture were evaluated, and they were compared with those for the saturated mixture. There was no large difference in qc between the saturated and superheated mixtures. The heat transfer coefficient hc (which was evaluated for the saturation temperature) for the superheated mixture agreed with that for the saturated mixture. hc computed with the diffusion layer model for saturated temperature agreed well with the hc measured for the superheated mixture.

Uncertainty Analysis on k-ε Turbulence Model in the Prediction of Subcooled Boiling in Vertical Pipes

Computational fluid dynamics (CFD) has become an effective method for researching two-phase flow in reactor systems. However, the uncertainty analysis of Computational fluid dynamics simulation is still immature. The effects of uncertainties from two-phase models and boundary conditions have been analyzed in our previous work. In this work, the uncertainties from a turbulence model on the prediction of subcooled boiling flow were analyzed with the DEBORA benchmark experiments by a deterministic sampling method. Seven parameters in the standard k-ε model, which interrelated momentum, energy, turbulent kinetic energy, and dissipation rate, were studied as uncertainty sources, including Cμ, Cμ,g, C1ε, C2ε, σk, σε, and Prt. Radial parameters were calculated to study the effects of uncertainties from the turbulence model. The contributions of each uncertainty source on void fraction and liquid temperature were also analyzed. It was found that the models can simulate subcooled boiling flow accurately and uncertainty analysis by deterministic sampling can give a reference interval to increase the reliability of results. The C2ε and C1ε, parameters in the production term and dissipation term of transport equations, dominate the radial distributions of void fraction and liquid temperature.