An official website of the United States government
Although supercritical CO2 (sCO2) heat transfer has been employed in industrial process since the 1960s, the underlying transport phenomenon in high-flux microscale geometries, as could be employed in concentrating solar receivers, is poorly understood. To date, nearly all experimental studies and simulations of supercritical convective heat transfer have focused on large diameter vertical channel and tube bundle flows, which may differ dramatically from microscale supercritical convection. Computational studies have primarily employed Reynolds averaged (RANS) turbulence modeling approaches, which may not capture effects from the sharply varying property trends of supercritical fluids. In this study, large eddy simulation (LES) turbulence modeling techniques are employed to study heat transfer characteristics of sCO2 in microscale heat exchangers. The simulation geometry consists of a microchannel of 750×737 μm cross-section and 5 mm length, heated from all four sides. Simulation cases are evaluated at reduced pressure P_r = 1.1, mass flux G = 1000 kg/m^2-s, heat flux q'' = 1.7 − 8.9 W/cm^2 , and varying inlet temperature: 20 − 100℃. Computational results reveal thermal transport mechanisms specific to microscale sCO2 flows. Results have been compared with available supercritical convection correlations to identify the most applicable heat transfer models for engineering of microchannel sCO2 heat exchangers.
more »
« less
This content will become publicly available on February 1, 2026
Unified, Integral Approach to Modeling and Design of High-Pressure Pipelines: Assessment of Various Flow and Temperature Models for Supercritical Fluids
Abstract A unified one-dimensional (1D), steady-state flow and heat transfer model is presented for the pipeline transport of fluids at high pressures, including the supercritical (SC) conditions. The model includes a generalized temperature equation, presented here for the first time, and accounts for all of the important effects, including the property variation, viscous dissipation, Joule-Thomson (J-T) cooling, and heat exchange with the surrounding. With appropriate approximations, this model can yield all isothermal and nonisothermal pipe flow solutions reported thus far. A generalized multizone integral method is developed which solves the two resulting algebraic equations for pressure and temperature in conjunction with a property database, such as the National Institute of Standard and Technology (NIST) reference fluid thermodynamic and transport properties (REFPROP). With appropriately selected number and size of the zones and using property values at the mean temperature and pressure within each zone, this integral method can accurately predict the complex effects of the governing parameters, such as the pipe diameter and length, inlet and exit pressures, mass flowrate, J-T cooling, and inlet and surrounding temperatures. Its accuracy for small-to-large diameter pipes has been ascertained by a comparison with the numerical solutions of the differential form of governing equations that requires a large number of small grids along the pipe and the values of mean properties within each grid. Indeed, this integral model can be used for the pipeline transport at both subcritical and supercritical pressures as long as the fluid does not encounter its anomalous states and the phase-change.
more »
« less
- PAR ID:
- 10636311
- Publisher / Repository:
- ASME
- Date Published:
- Journal Name:
- Journal of Fluids Engineering
- Volume:
- 147
- Issue:
- 2
- ISSN:
- 0098-2202
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Miniaturization and high heat flux of power electronic devices have posed a colossal challenge for adequate thermal management. Conventional air-cooling solutions are inadequate for high-performance electronics. Liquid cooling is an alternative solution thanks to the higher specific heat and latent heat associated with the coolants. Liquid-cooled cold plates are typically manufactured by different approaches such as: skived, forged, extrusion, electrical discharge machining. When researchers are facing challenges at creating complex geometries in small spaces, 3D-printing can be a solution. In this paper, a 3D-printed cold plate was designed and characterized with water coolant. The printed metal fin structures were strong enough to undergo pressure from the fluid flow even at high flow rates and small fin structures. A copper block with top surface area of 1 inch by 1 inch was used to mimic a computer chip. Experimental data has good match with a simulation model which was built using commercial software 6SigmaET. Effects of geometry parameters and operating parameters were investigated. Fin diameter was varied from 0.3 mm to 0.5 mm and fin height was maintained at 2 mm. A special manifold was designed to maximize the surface contact area between coolant and metal surface and therefore minimize thermal resistance. The flow rate was varied from 0.75 L/min to 2 L/min and coolant inlet temperature was varied from 25 to 48 oC. It was observed that for the coolant inlet temperature 25 oC and aluminum cold plate, the junction temperature was kept below 63.2 oC at input power 350 W and pressure drop did not exceed 23 Kpa. Effects of metal materials used in 3D-printing on the thermal performance of the cold plate were also studied in detail.more » « less
-
All fluids exhibit large property-variations near the critical point in a region identified as the anomalous state. The anomaly starts in the liquid and extends well into the supercritical state, which can be identified thermodynamically using the Gibbs free energy (g). The specific heat, isobaric expansion, and isothermal compressibility parameters governing the transitions are: (cp/T), (vβ), and (vκ), rather cp, β, and κ. They are essentially the second-order derivatives of g and have two extrema (minimum, maximum); only maxima reported ever. When applied to the van der Waals fluid, these extrema exhibit closed loops on the phase-diagram to satisfy d3g = 0 and map the anomalous region. The predicted liquid-like to gas-like transitions are related to the ridges reported earlier, and the Widom delta falls between these loops. Evidently, in the anomalous region, both the liquid and the supercritical fluid need to be treated differently. Beyond the anomalous states, the supercritical fluids show monotonic, gradual changes in their properties. The analysis for argon, methane, nitrogen, carbon dioxide, and water validates the thermodynamic model, supports the stated observations, and identifies their delimiting pressures and temperatures for the anomalous states. It also demonstrates the applicability of the law of corresponding states. Notably, the critical point is a state where d3g = 0, the anomaly in the fluid's properties/behavior is maximal, and the governing parameters approach infinity. Also the following are presented: (a) the trajectory of the liquid–vapor line toward the melt-solid boundary and (b) a modified phase diagram (for water) exhibiting the anomalous region.more » « less
-
Abstract This paper proposes a computational fluid dynamics (CFD) simulation methodology for the multi-design variable optimization of heat sinks for natural convection single-phase immersion cooling of high power-density Data Center server electronics. Immersion cooling provides the capability to cool higher power-densities than air cooling. Due to this, retrofitting Data Center servers initially designed for air-cooling for immersion cooling is of interest. A common area of improvement is in optimizing the air-cooled component heat sinks for the fluid and thermal properties of liquid cooling dielectric fluids. Current heat sink optimization methodologies for immersion cooling demonstrated within the literature rely on a server-level optimization approach. This paper proposes a server-agnostic approach to immersion cooling heat sink optimization by developing a heat sink-level CFD to generate a dataset of optimized heat sinks for a range of variable input parameters: inlet fluid temperature, power dissipation, fin thickness, and number of fins. The objective function of optimization is minimizing heat sink thermal resistance. This research demonstrates an effective modeling and optimization approach for heat sinks. The optimized heat sink designs exhibit improved cooling performance and reduced pressure drop compared to traditional heat sink designs. This study also shows the importance of considering multiple design variables in the heat sink optimization process and extends immersion heat sink optimization beyond server-dependent solutions. The proposed approach can also be extended to other cooling techniques and applications, where optimizing the design variables of heat sinks can improve cooling performance and reduce energy consumption.more » « less
-
null (Ed.)An in-rack cooling system connected to an external vapor recompression loop can be an economical solution to harness waste heat recovery in data centers. Validated subsystem-level models of the thermosyphon cooling and recompression loops (evaporator, heat exchangers, compressor, etc.) are needed to predict overall system performance and to perform design optimization based on the operating conditions. This paper specifically focuses on the model of the evaporator, which is a finned-tube heat exchanger incorporated in a thermosyphon cooling loop. The fin-pack is divided into individual segments to analyze the refrigerant and air side heat transfer characteristics. Refrigerant flow in the tubes is modeled as 1-D flow scheme with transport equations solved on a staggered grid. The air side is modeled using differential equations to represent the air temperature and humidity ratio and to predict if moisture removal will occur, in which case the airside heat transfer coefficient is suitably reduced. The louver fins are modeled as individual hexagons and are treated in conjunction with the tube walls. A segment-by-segment approach is utilized for each tube and the heat exchanger geometry is subsequently evaluated from one end to the other, with air property changes considered for each subsequent row of tubes. Model predictions of stream outlet temperature and pressure, refrigerant outlet vapor quality and heat exchanger duty show good agreement when compared against a commercial software.more » « less
