Degree Of Subcooling Research Articles

Thermodynamic analysis and performance optimization on an ultra‐low‐temperature cascade refrigeration system using refrigerants R290 and R170

AbstractTo meet the strict requirement on whole chain of vaccine production, storage, transportation, and distribution, most researches have been done on cascade refrigeration systems to achieve high operation performance, while a research gap on the performance exploration of eco‐friendly refrigeration systems still exists. In this paper, a comprehensive thermodynamic model was established to analyze the operation performance of a pre‐cooled cascade refrigeration system with eco‐friendly refrigerants propane (R290) and ethane (R170). Based on the present thermodynamic model, the performance optimization on the R290‐R170 cascade refrigerator was made with considerations of degree of subcooling, degree of superheat, evaporation temperature, condensation temperature, cascade condensation temperature, and cascade temperature difference. Variations of the coefficient of performance, exergy destruction, and total exergy efficiency of the refrigeration cycle were analyzed. Two mathematical correlations yielding the optimal cascade condensation temperature and maximized coefficient of performance were developed by multilinear regression analysis. When the evaporation temperature is − 60°C, the maximized coefficient of performance and total exergy efficiency are 1.276 and 49.51%. This paper demonstrates the potential for improving the R290‐R170 cascade refrigeration system and furnishes the basis for further exploration on ultra‐low‐temperature refrigerators.

Energy, exergy, and economic analysis of a novel solar-assisted ejector subcooling CO2 transcritical refrigeration systemt

To further enhance the performance of transcritical CO2 refrigeration cycle (base cycle), a solar-assisted ejector subcooling transcritical CO2 refrigeration cycle (SESRS) is proposed. The base cycle is subcooled utilizing a solar-assisted ejector refrigeration cycle. This research performed energy, exergy, and economic evaluations of SESRS utilizing the eco-friendly refrigerant R152a in the ejector cycle based on the established model. The thermodynamic analysis results demonstrate that SESRS outperform the base cycle. At typical working conditions, the COPm of the SESRS demonstrates a 22.96 % increase compared to that of the base cycle. Although the total cost of the SESRS is 7.98 %–17.9 % higher than that of the base cycle when subcooling degree is 5 °C, the enhancement of COPm is larger and is increased by approximately 18.34 %–26.22 %. The impact of subcooling degree and other parameters on system performance are also investigated. Overall, this study confirms the potential of the SESRS system for utilization in the refrigeration and air-conditioning fields. Considering both the system’s performance and economy, further optimized analysis should be done on key operational parameters, such as subcooling degree and discharge pressure.

A Combined Experimental and Computational Study on the Effect of the Reactor Configuration and Operational Procedures on the Formation, Growth and Dissociation of Carbon Dioxide Hydrate

Clathrate hydrate-based technologies are considered promising and sustainable alternatives for the effective management of the climate change risks related to emissions of carbon dioxide produced by human activities. This work presents a combined experimental and computational investigation of the effects of the operational procedures and characteristics of the experimental configuration, on the phase diagrams of CO2-H2O systems and CO2 hydrates’ formation, growth and dissociation conditions. The operational modes involved (i) the incremental (step-wise) temperature cycling and (ii) the continuous temperature cycling processes, in the framework of an isochoric pressure search method. Also, two different high-pressure PVT configurations were used, of which one encompassed a stirred tank reactor and the other incorporated an autoclave of constant volume with magnetic agitation. The experimental results implied a dependence of the subcooling degree, (P, T) conditions for hydrate formation and dissociation, and thermal stability of the hydrate phase on the applied temperature cycling mode and the technical features of the utilized PVT configuration. The experimental findings were complemented by a thermodynamic simulation model and other calculation approaches, with the aim to resolve the phase diagrams including the CO2 dissolution over the entire range of the applied (P, T) conditions.

Mechanistic insights into pore water conversion to gas hydrates in clay minerals

The conversion of pore water to gas hydrates is an important prerequisite for natural gas hydrates (NGHs) accumulation and hydrate-based CO2 sequestration. However, the mechanism of water-to-hydrate conversion in clays, which are prevalent in hydrate reservoirs or marine sediments serving as potential hosts for CO2 hydrates, still remains poorly understood, impeding further advancements in exploiting NGHs and developing CO2 sequestration strategies. In this paper, the characteristics of pore water conversion to hydrates in montmorillonite and kaolin were comprehensively examined through quantitative measurements and NMR analysis, and compared with those in silt and coarse sand. It revealed that the water-to-hydrate conversion was more restricted in clays than in sands. Especially in montmorillonite, since almost all pore water was detected as bound water, the final water conversion was extremely low and sensitive to the subcooling degree. Additionally, as the initial water content increased from 15 wt% to 35 wt%, the final water conversion in kaolin, silt, and coarse sand decreased exponentially due to intensified kinetic limitation, whereas that in montmorillonite initially rose and then fell, attaining the maximum at a threshold of 20 wt%. It is found that the competition between thermodynamic and kinetic limitations controls water-to-hydrate conversion in clays. Specifically, in kaolin, silt and coarse sand, the main controlling factor of water-to-hydrate conversion is the kinetic limitation, as their unconverted equilibrium water is negligible. However, in montmorillonite, the thermodynamic limitation dominates at low water contents, transitioning to kinetic limitation at high water contents. This study provides direct evidence and mechanistic insights into how clays affect water-to-hydrate conversion in porous media, which has significant implications for the geological accumulation of NGHs and the site selection for hydrate-based CO2 sequestration.

Ultra-stable counter flow diverging minichannel heat sink with integrated microstructures for superior cooling performance

Efficient thermal management of miniaturized electronic devices is quite challenging nowadays. It requires dissipation of very high heat flux without considerable elevation of surface temperature. To address this issue, this study proposes the counter flow diverging minichannel heat sink in a compact size of 3 cm2. Additionally, two innovative microstructures have been developed and integrated into the bottom surface of the proposed minichannel heat sink: one with gradient micro pin-fins and another combining gradient microcavities with micro pin-fins. Comprehensive flow boiling experiments are conducted under various operating conditions to identify the best configuration with the most appropriate operating condition to dissipate the highest heat flux with minimal pressure drop and stable two-phase flow patterns. The current study reveals that the synergistic effects of gradient micro pin-fins and optimized gradient microcavities achieve a very high heat flux of 746.14 W/cm2 and a two-phase heat transfer coefficient of 28.18 W/cm2K under the limitation of wall temperature lower than 132 °C. Additionally, it maintains an exceptionally low pressure drop of 2.35 kPa with stable two-phase flow patterns under a mass flux of 600 kg/m2s and a degree of subcooling of 40 °C. Compared to the bare surface, counter flow diverging minichannels with gradient micro pin-fins and optimized gradient microcavities show excellent enhancements of 201.12 % in the effective heat flux and 124 % in the two-phase heat transfer coefficient, suggesting it is the best design in the present study for advanced cooling applications.

Thermodynamic Comparative Analysis of Cascade Refrigeration System Pairing R744 with R404A, R448A and R449A with Internal Heat Exchanger: Part 2—Exergy Characteristics

The cascade refrigeration systems (CRS) used in hypermarkets and supermarkets, which are used by many people, have been employing R744 for the low-temperature cycle (LTC) and R404A for the high-temperature cycle (HTC) due to environmental and public safety issues. However, the use of R404A is limited due to its high GWP, and therefore research on alternative refrigerants is necessary. Nevertheless, there is no detailed study in the literature that compares and analyzes the three refrigerants for practical design by applying R744 for LTC and R404A, R448A, and R449A for HTC in CRS. Therefore, this study aims to provide data for the practical detailed design of an alternative system to R744/R404A CRS. Under standard conditions, we analyzed how the exergy destruction rate (EDR) and exergy efficiency (EE) of the system and the EDR of each component change when the important factors affecting CRS (degree of superheating (DSH), degree of subcooling (DSC), and internal heat exchanger (IHX) efficiency of HTC, DSH of LTC, condensation temperature (CT), evaporation temperature (ET), cascade evaporation temperature (CET), and temperature difference of CHX) are varied over a wide range. The main conclusions are as follows. (1) Under the given constant conditions, the smallest change in system EDR based on R448A is DSH of HTC (decreased by 0.07–0.1 kW), followed by IHX of HTC (decreased by 0.12–0.3 kW), DSH of LTC (increased by 0.19–0.25 kW), DSC of HTC (decreased by 0.59–0.69 kW), temperature difference of CHX (increased by 1.57–1.83 kW), CET (decreased and then increased by 0.67–4.43 kW), CT (increased by 1.49–3.9 kW), ET (decreased by 2.39–4.61 kW). (2) The highest change rate of system EE based on R448A is CET (increased and then decreased by 1.38–8.28%), followed by temperature difference of CHX (decreased by 2.96–3.16%), ET (increased and then decreased by 0.63–2.75%), DSC of HTC (increased by 1.26–1.34%), CT (increased and then decreased by 0.24–1.12%), IHX of HTC (increased by 0.11–1.02%), DSH of LTC (decreased by 0.35–0.49%), and DSH of HTC (increased by 0.14–0.19%).

Thermodynamic analysis of a modified two-stage transcritical CO2 refrigeration cycle with an ejector and a subcooler

In the application scope of large-scale supermarkets, the practicality of transcritical CO2 refrigeration device has been confirmed to be quite satisfying. A modified two-stage transcritical CO2 refrigeration cycle with an ejector and a subcooler (MTC) is proposed in this paper. In the modified cycle, the ejector recovers expansion works and reduces irreversible losses in the throttling process. The subcooler provides subcooling degree for the CO2 entering the low-temperature (LT) evaporator, thus increasing the refrigeration capacity and improving the COP of the modified cycle. Thermodynamic analysis has shown that the exergy efficiency (ηex) and COP of MTC under given operating condition has been enhanced by 13.7 % and 14.2 % compared to BTC. The discharge temperature at the outlet of high-pressure compressor in MTC is decreased by 8.3℃. Under typical operating condition, the optimal discharge pressure of MTC is 9.07 MPa, which is lower than that of BTC. The correlation to calculate optimal discharge pressure for single-stage CO2 refrigeration cycle is also suitable for MTC under given operating conditions. The MTC has also shown better performance under variable operating conditions. For MTC, when the temperature at the outlet of gas cooler increases from 35 – 45℃, the COP and ηex are enhanced by 14.1 % - 16.1 % and 12.8 % - 14.2 % compared to BTC, respectively. When gas cooler outlet pressure decreases from 11.0 to 7.5 MPa, the COP and ηex are enhanced by 14.0 % - 22.8 % and 13.4 % - 18.7 %. As the evaporating temperature at the cold side of subcooler increases from -25 to -15℃, the COP and ηex increase from 1.82 to 1.98 and 27.4 % to 29.7 %. With the ratio of refrigeration capacity (Rec) between medium-temperature evaporator and low-temperature evaporator varies from 0.7 to 1.2, the COP and ηex are improved by 10.7 % - 16.3 % and 10.7 % - 15.4 % compared to BTC. There is the maximum exergy loss at gas cooler in the MTC, whereas that of BTC is located in the expansion valve before LT evaporator. The economic analysis shows the cost per unit of exergy of MTC is decreased by 11.5 % under typical operation condition. According to simulation results, the modified cycle has better performance in severe working conditions such as high gas cooler outlet temperature and low gas cooler outlet pressure in the given range of working conditions compared to BTC.

4E Analysis of a new cogeneration system coupling heat pumps and PEMFC in summer and winter modes

In response to the ongoing optimization of energy systems and the promotion of clean energy, this paper introduces a new high-efficiency cogeneration system based on proton exchange membrane fuel cells (PEMFC). Named the Transcritical Combined Cooling Heating and Power with Subcooling Coupled ORC (CPTSO) system, it undergoes energy, exergy, economic, and environmental (4E) assessments to evaluate its performance in both summer and winter modes. This paper also analyzes the effects of key factors on the new system's performance. The results indicate that the new system performs more effectively overall and achieves faster cost recovery during winter operations. With varying degrees of subcooling (ΔTsubcooling), the system's average energy efficiency, fuel energy saving ratio (FESR), exergy efficiency, and pollutant reduction in winter are 60.97 %, 4.51 %, 17.16 %, and 33.35 % higher, respectively, than in summer. Changes in the main evaporation temperature (Te) result in winter improvements of 74.37 % in energy efficiency, 16.86 % in FESR, 17.6 % in exergy efficiency, and 45.5 % in pollutant reduction compared to summer. Similarly, adjustments to the outlet temperature of the gas cooler (Tg) lead to winter increases of 56.75 % in energy efficiency, 0.49 % in FESR, 17.1 % in exergy efficiency, and 29.3 % in pollutant reduction over summer values.

The effect of subcooling on hydrate generation and solid deposition with two-phase flow of pure water and natural gases

Deepwater oil and gas transportation systems play a crucial role in energy supply, but also face the challenges of hydrate generation and blockage. The study observed hydrate generation and deposition at different subcooling degrees and flow rates by using a high-pressure fully visualized recirculation pipeline. It was found that the subcooling degree has a nonlinear positive correlation with the total amount of hydrate generated. The hydrate formation at 11.7 °C subcooling was 4.5 times higher than that at 9 °C subcooling. According to the hydrate conversion rate analysis, the hydrate growth process under low subcooling conditions directly transitions from the bubble-promoting zone to the termination point of hydrate generation. Combined with the deposition state of the pipe wall and the differential pressure change curves, the risk of blockage at different subcooling levels and flow rates can be assessed.

Time-resolved modeling of dropwise condensation patterns formed on a nanopillared substrate

Instantaneous condensation patterns formed over a vertical nanopillared copper surface are modeled by including drop level details including nucleation, growth, coalescence, droplet jumping and gravitational instability. The vapor phase is taken to be saturated and condensation is driven by a colder wall that imposes a given degree of subcooling. The mathematical model is setup in the form of a condensation cycle, starting from drop nucleation all the way to instability followed by fresh nucleation. In view of the condensing surface being pillared, a variety of droplet morphologies are possible, including Cassie, Wenzel, and partially-wetting that depend on the prevailing condensation conditions. The degree of subcooling is varied over 1 – 28 K, covering a wide range of condensation characteristics from Cassie jumping to Wenzel flooding with increasing temperature difference. Coalescence-based droplet jumping for Cassie and partially-wetting droplet configurations are incorporated as well. Gravitational instability of the largest drop leads to sweeping along the substrate followed by surface renewal and renucleation. With this overall approach, it is possible to explicitly determine the drop-size distribution as a function of time. The entire model is multiscale in space and time, making it a computationally demanding simulation. A domain decomposition framework is used and computations are carried out in a high-performance computing system with a parallel architecture employing message passing interface (MPI). The time-resolved model has been extensively tested against experimental data reported in the literature. The partially-wetting configuration for a subcooling of 8 K has a larger average droplet size and fewer jumping events generating similar heat fluxes as for Cassie droplets at a lower degree of subcooling (ΔT = 5 K) where the average droplet size is smaller and jumping events are frequent. In addition, the Cassie state yields an enhancement factor of around 2.5 relative to a flat surface while Wenzel droplets cause flooding and need subcooling levels as high as 24 K for similar heat flux values.

Theoretical analysis of the initial nucleation characteristics of trace water vapor frosting on a cryogenic surface

In certain extreme conditions characterized by ultra-low water vapor content and ultra-low temperatures (e.g. cryogenic wind tunnel), trace water vapor frosting can also pose significant hazards. Considering the crucial role of initial nucleation on subsequent frosting processes, this study firstly investigated and compared the nucleation characteristics of frosting at different water vapor contents (0.1−1000 ppmv) from humid air to trace water vapor based on classical nucleation theory. The preferred nucleation phase diagrams and critical nucleation conditions were analyzed. Sensitivity of the nucleation mass transfer rate, as well as the frosting pathways were further identified. Results show that the nucleation characteristics of trace water vapor frosting (<100 ppmv) is significantly different from that of humid air frosting (>1000 ppmv). Trace water vapor frosting is more inclined towards desublimation nucleation due to its lower nucleation temperature and larger critical contact angle. The critical contact angles for 1000 ppmv, 100 ppmv and 0.1 ppmv are 49°, 75° and 120°, respectively. Furthermore, trace water vapor nucleation requires a greater subcooling degree, has a smaller critical nucleation radius, and is highly sensitive to surface contact angle and further-subcooling degree. At a surface contact angle of 120°, the nucleation subcooling degree of 0.1 ppmv is 2.9 K greater than that of 1000 ppmv. This study helps understanding the nucleation characteristics under different water vapor content conditions, which indicates that reducing subcooling degree and increasing contact angle are more effective anti-frosting methods than reducing water vapor content for trace water vapor frosting.

Experimental Investigation of R404A Indirect Refrigeration System Applied Internal Heat Exchanger: Part 2—Exergy Characteristics

Although the R404A indirect refrigeration system (IRS) with an internal heat exchanger (IHX) and R744 as the secondary fluid has potential applications in supermarkets and hypermarkets, the exergy characteristics of this IRS have not been extensively investigated. In this study, the factors affecting the R744 exergy characteristics (degree of subcooling (DSB) and degree of superheating (DSP) of the R404A cycle, DSP of the R744 cycle, condensation temperature (CT) and cascade evaporation temperature (CET), and IHX efficiency) were experimentally evaluated to obtain basic data for the design of R404A IRS with R744 as the optimal secondary fluid. The main results can be summarized as follows: (1) Under given conditions, the smallest change in the system exergy destruction rate (EDR) according to the change in each parameter is the DSP of the R744 cycle (0.3–1%), followed by the DSB of the R404A cycle (6.1–8.8%), the IHX efficiency of the R404A cycle (3.8–14.3%), the DSP of the R404A cycle (11.7–15.9%), the CET (29.4–41.9%), and the CT (35–47%). (2) Also, in terms of the exergy efficiency of system (EES), the largest value was obtained for the DSP of the R404A cycle (2.4–12.7%), followed by the IHX efficiency of the R404A cycle (3–10.2%), the CET (2.2–8.7%), the CT (4–6.9%), the DSB of the R404A cycle (2.7–6.2%), and the DSP of the R744 cycle (0.04–1.2%). (3) In order to lower the system EDR, DSP, DSB, and IHX efficiency of the R404A cycle, the CET must be increased to the maximum, and to lower the DSP of the R744 cycle, the CT must be reduced to the minimum.

Determining the critical time points for hydrate formation in the presence of kinetic hydrate inhibitors: A rheological and kinetic study

Kinetic hydrate inhibitors have the characteristics of delaying hydrate nucleation and reducing hydrate growth rate. The utilized of kinetic hydrate inhibitors for hydrate blockage risk management is a promising technology for hydrate prevention in oil and gas industry. In this study, the rheological and kinetic properties of gas–water systems contained pure water, polyvinyl pyrrolidone and newly-synthesized inhibitor during hydrate formation at different subcooling degrees were in situ characterized using a high-pressure rheometer. By analyzing the experimental results, the concept of ‘critical time’ have been proposed to divide the time-varying process of hydrate formation into multiple stages. The three critical times indicated distinct trends in viscosity and pressure changes, as well as varying risks of hydrate blockage in hydrate prevention. The results further compared the differences in determining critical time through pressure and viscosity. The critical time of hydrate nucleation determined by viscosity was earlier than pressure due to the hypothesis of lag in gas consumption. The critical time of rapid hydrate growth determined by viscosity was later than gas consumption, which may be due to that the sharp increase in viscosity requires the aggregation of hydrates. According to the inhibition effect evaluation results determined by three critical times, it was found that the newly-synthesized inhibitor was more effective than polyvinyl pyrrolidone in inhibiting nucleation, delaying rapid hydrate growth time and reaching high viscosity time. This study can provide a comprehensive evaluation method for assessing the performance of kinetic hydrate inhibitors from both kinetic and rheological perspectives.

Experimental investigation on hydrated salt phase change material for lithium-ion battery thermal management and thermal runaway mitigation

To efficiently ensure the suitable temperature for lithium-ion batteries, battery thermal management system based on phase change material (PCM) have gained considerable attention. A shape-stabilized composite PCM was prepared using the mixture of Na2SO4·10H2O and KAl(SO4)2·12H2O hydrated salts as phase change substrate, expanded graphite as thermal conductive filler and open-cell polyurethane foam as support material. The composite PCM containing 4 wt.% expanded graphite displayed a satisfied phase change temperature of 30.8 °C and 69.8 °C, high latent heat of 226.9 J/g, and good thermal conductivity of 1.39 W/(m·K). Meanwhile, it obtained a low subcooling degree, good leakage resistance, and excellent flame retardancy properties. Under PCM cooling in sandwich structure, during 0.5C charging and 1.5C discharging cycles, the maximum temperature and maximum temperature difference of battery were 52.47 °C and 2.20 °C, decreased by 20.4 % and 59.7 % respectively compared with nature air cooling. Furthermore, the composite PCM could significantly relieve the thermal runaway propagation between lithium-ion batteries. This work provides a feasible design and application strategy of hydrated salt based PCMs for battery thermal management and thermal runaway mitigation.

Benchmark DEBORA: Assessment of MCFD compared to high-pressure boiling pipe flow measurements

A benchmark activity on two-fluid simulations of high-pressure boiling upward flows in a pipe is performed by 12 participants using different MCFD (Multiphase Computational Fluid Dynamics) codes and closure relationships. More than 30 conditions from DEBORA experiment conducted by CEA are considered. Each case is characterised by the flow rate, inlet temperature, wall heat flux and outlet pressure. High-pressure Freon (R12) at 14bar and 26bar is boiled in a 19.2mm pipe heated over 3.5m. Flow rates range from 2000kgm−2s−1 to 5000kgm−2s−1 and exit quality x ranges from single-phase conditions to x=0.1 which leads to a peak void fraction of α=70%. In these high pressure conditions, bubbles remain small and there is no departure from the bubbly flow regime (François et al., 2011; Hösler, 1968). However, different kind of bubbly flows are observed: wall-peak, intermediate peak or core-peak, depending on the case considered. Measurements along the pipe radius near the end of the heated section are compared to code predictions. They include void fraction, bubble mean diameter, vapour velocity and liquid temperature. The benchmark covered two phases. In the first phase of the benchmark activities, experimental data were given to the participants, allowing to compare the simulation results and to develop, to select or to adjust the models in the CMFD codes. The second phase included blind cases where the participants could not compare to the measurements. In between the two phases, possible additional model adjustments or calibrations were performed.Overall, the benchmark involved very different closures and a wide range of models’ complexity was covered. Yet, it is extremely difficult to have a robust closure for all conditions considered, even knowing experimental measurements. The wall-to-core peak transition is not captured consistently by the models. The degree of subcooling and the void fraction level are also difficult to assess. We were not capable of showing superiority of some physical closures, even for part of the model. The interaction between mechanisms and their hierarchy are extremely difficult to understand.Although departure from nucleate boiling (DNB) was not considered in this benchmarking exercise, it is expected that DNB predictions at high-pressure conditions depend strongly on the near-wall flow, temperature, and void fraction distributions. Therefore, the suitability of the closures also limits the accuracy of DNB predictions. The benchmark also demonstrated that in order to progress further in models development and validation, it is compulsory to have new measurements that include simultaneously as many variables as possible (including liquid temperature, velocity, cross-correlations and wall temperature); also, a better knowledge of the local bubble sizes distributions is the key to discriminate performances of interfacial area modelling (IATE, MUSIG or iMUSIG models, considering for instance the possibility of two classes of bubbles having totally different behaviour regarding the lift force).Following this benchmark impulse, we hope that future activities will be engaged on high-pressure boiling water experiments with a continuation of models’ comparisons and development.

Numerical investigation on dynamic flow characteristics of methane condensation in microchannels

Microchannel condensers play an essential role in cryogenic two-phase heat management systems due to their efficient heat transfer characteristics. Thus, it is worth conducting an in-depth study on microscale condensation characteristics of cryogenic fluids. This paper delves into the flow condensation process of methane in microchannels. A two-dimensional transient model with high accuracy for cryogenic fluids has been developed by combining a self-defined program for the source term of the phase transition model. The model fully considers the boundary layer thickness and accurately explores the mesh accuracy. The complete condensation flow patterns are captured for various vapor quality, mass flux, and wall subcooling degrees. The injection flow is a unique flow regime for condensation in microchannels. The decrease in wall subcooling degree and increase in mass flux leads to the separation point at the neck of the injected flow moving towards the exit, while the annular flow region is expanding and the flow pattern transition is lagging. The mass flux improves the heat transfer coefficient more significantly at high vapor quality. During injection and bubble flow, the wall shear stress and local heat transfer coefficient are subject to bouncing and oscillations, which may induce fluctuations in the upstream annular flow. The prediction performance of six classical heat transfer correlations is evaluated. The results indicate that the Nie et al. correlation has the highest comprehensive prediction accuracy with MRD and MARD of -5.00 % and 15.83 %, respectively.

CFD elucidation of high-pressure subcooled boiling flow towards effects of variable refrigerant properties using OpenFOAM empirical closures

Boiling flow presents a significant concern, especially when a liquid surpasses its boiling point, potentially leading to catastrophic consequences. This research utilizes a two-phase code in the OpenFOAM software to investigate bubble formation during flow boiling. The well-established empirical models for calculating wall heat components were selected based on the operating conditions. The study incorporates experimental data from high-pressure boiling flow (10–30 bars) with variable properties of refrigerant R-12. The predictions reveal underpredictions in void fraction and liquid temperature compared to experimental observations. Significantly, the impact of the subcooling degree on void fraction behaviour is emphasized, and a potential underprediction of the evaporation portion is highlighted, particularly near the wall. Challenges in modelling bubble size distribution are evident through discrepancies in bubble diameter and velocity data, indicating the necessity for further advancements in the code. In summary, this numerical study provides valuable insights into the intricate dynamics of high-pressure subcooled boiling flow, especially when considering variable working fluid properties. Future efforts will focus on refining models for nucleation site density, bubble departure size, and lift-off frequency to enhance prediction accuracy.

Effect of Plasma Electric Field Assistance on the Freezing Process and Ice Crystal Formation of Bactrian Camel Meat

Background: The longissimus dorsi muscle of the camel was chosen as the experimental material in order to investigate the impact of plasma electric field assisted freezing on the quality and freezing process of Bactrian camel meat. Methods: The effect of the auxiliary plasma electric field of different field strengths (1 kV/m, 3 kV/m and 5 kV/m) on the freezing velocity of the camel meat and the distribution of ice crystals as well as the diameter and roundness of ice crystal in the frozen meat were studied and analyzing them with the help of the thermodynamic theory. Result: The findings demonstrated that the meat samples in the auxiliary 5 kV/m plasma field environment took 60 min and 30 min, respectively, at -20°C and -35°C, to pass the maximal ice crystal production zone. The tissue sections revealed that at -20°C and -35°C assisted 5 kV/m plasma field freezing, the smallest ice crystal area and grain size were formed, the roundness of the ice crystals was highest. Using Boltzmann’s entropic equation to calculate the system’s entropy, the results revealed that the meat samples frozen by the aided 5 kV/m plasma electric field had the lowest entropy at -20°C and -35°C. The electric field assisted freezing requires less energy to improve the subcooling degree and accelerate the freezing rate of meat samples, proving that the auxiliary 5 kV/m plasma electric field can significantly improve the system frozen meat moisture vitrification degree and can make the crystallization process more controllable.

Composite microstructured surface with micro-cavities and micro-ditches on micro-pin-fins for enhancing pool boiling heat transfer

In this study, a novel composite microstructured surface with micro-cavities and micro-ditches on micro-pin-fins was designed to enhance pool boiling heat transfer. To reveal the enhancement mechanism of the newly-designed composite microstructured surface, the pool boiling heat transfer performance on four types of microstructured surfaces was investigated in HFE-7100 under different liquid subcooling degrees. The experimental results show that the micro-pin-finned surfaces with micro-ditches at the side walls of each micro-pin–fin can enhance the rewetting ability of the surfaces, thereby enhancing the boiling heat transfer at high heat fluxes and improving the critical heat flux (CHF) in comparison with the base micro-pin-finned surface. Meanwhile, the wickability of the micro-pin-finned surfaces with micro-ditches was found to increase with increasing the number of micro-ditches. Furthermore, by fabricating both micro-cavities and micro-ditches on each micro-pin–fin, the composite microstructured surface with micro-cavities and micro-ditches is capable of not only utilizing the advantage of increasing the nucleation site density through the micro-cavity structures, but also enhancing the rewetting ability of the surface through the micro-ditch structures. As a result, the newly-designed multiscale composite surface can significantly improve the heat transfer coefficient (HTC) and the CHF simultaneously. Quantitatively, the achieved largest HTC and highest CHF are 5.12 W/(cm2·K) and 75.2 W/cm2, respectively, which are 6.24 times and 2.82 times, respectively, as high as those on a smooth surface at the liquid subcooling degree of 35 K.

Experiment on the Excitation of Vapor Behavior and the Vibration by Thermal Flow Under C-Tube Pool Boiling Model

Abstract As an important safety equipment in nuclear system, the passive residual heat removal heat exchanger (PRHR HX) with C-tube is responsible for timely removing heat from the reactor core in accident. A C-tube heating device was setup to investigate the flow characteristics of the in-containment refueling water storage tank (IRWST), which is of the PRHR HX shell-side flow field and the vibration characteristics of the heat exchanger tube. The temperature of the flow field and the acceleration signal of the C-tube were measured. The C-tube frequency in the pool boiling condition is 0–0.660 Hz. There is no resonance with the first natural frequency of 8.98 Hz. The IRWST temperature field can be divided into slow-heating period, fast-heating period, and boiling period. Steam vapor exists in three periods, condensation and collapse of vapor concentrate on horizontal tube. The fast-heating period exists the subcooling degree ΔTsub for the most intense vapor excitation, ΔTsub is in the range of 13–16 °C. In term of thermal flow excitation vibration, the existence of the cross flow with C-plane direction is in the transition period from the slow-heating period to the fast-heating period, which excites the maximum root mean squares (RMS) amplitude 0.039 mm and 0.034 mm at the midpoint of upper horizontal tube and vertical tube. When the fluid enters fast-heating period, the direction of flow changes into the normal C-plane. Normal RMS amplitude reached maximum 0.02 mm at upper 1/4 of the vertical tube in boiling period. The above research provides important data support for the safety and reliability of PRHR HX.