SEARCH · Search NASA

Results for “INLET FLOW”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Impact of particle size and particle-flow-wall coupling on pressurized oxy-combustion in the down-fired burner

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the burner of a pressurized oxy-combustor is close or even higher than the maximum limit that commercial CFD codes (e.g., ANSYS FLUENT) can handle. At this high particle volume fraction, the interactions among particles, fluid flow, and wall need to be re-evaluated. The present computational work is the first step of a systematic analysis of the particle influence, like releasing method, releasing locating, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing ANSYS FLUENT, using Reynolds-averaged Navier-Stokes (RANS) modeling. It is revealed that for this pilot-scale, pressurized burner, velocity profiles of the near-wall region in this POC facility exhibit some discrepancy against the near-wall turbulent flow velocity profile. In order to investigate the particle influence in the near-wall region, particle releasing location will be investigated. Numerical simulation results exhibit the coupling effect of turbulence flow and particles in this case. The particle size also demonstrates a great effect on particle trajectory, then further has an impact on flame stability and temperature.

Li, Lei↗

DOE OSTI2022

Thermo-economic analysis of a solid oxide fuel cell-gas turbine hybrid with commercial off-the-shelf gas turbine

The thermodynamic and economic performance of a natural gas-fueled solid oxide fuel cell (SOFC)-gas turbine (GT) hybrid system with a commercial off-the-shelf GT is investigated. Until today, commercial GTs are primarily engineered for the direct use of energy dense fuels, such as natural gas, and the integration of commercial GTs into a SOFC-GT hybrid remains challenging. In this work technically feasible and economically viable GT operating modes are identified to gauge the technology’s competitiveness on the free market. Steady state, full load operating conditions were established for various GT operating modes: I) constant spool speed operation, II) variable spool speed operation and III) partially closed compressor inlet guide vanes. For each GT operating mode, the performance was investigated over a range of SOFC fuel utilization factors, while considering physical constraints inside the SOFC, such as local temperature gradients in flow direction as well as overall cell temperature differences. The results of the thermodynamic evaluation served as inputs for the economic analysis of the SOFC-GT hybrid power plant. The results show that the integration of an off-the-shelf GT not only results in a significant derating of the GT, but also substantially impacts the SOFC operation, the main power producer in this SOFC-GT hybrid. To maximize the SOFC power output it is desirable operate the GT in a region of high air mass flow rates and low pressure ratios, which increases the number of stacks that can be accommodated and reduces the SOFC cooling requirement. The lowest costs of electricity are obtained at constant spool speed operation, while the highest efficiencies are reached at variable spool speed operation. Critical for the integration of off-the-shelf GTs, which historically have been designed for natural gas, is the surge margin. The largest surge margins are obtained by closing the compressor inlet guide vanes. Furthermore, higher SOFC fuel utilization factors are shown to increase the surge margin as the turbine firing temperature is decreased.

25 ENERGY STORAGE↗

DOE OSTI2020

Heat Transfer Coefficients of Additively Manufactured Tubes with Internal Pin Fins for Supercritical Carbon Dioxide Cycle Recuperators

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO 2 flowing in additively manufactured tubes with internal pin fins at the US DOE’s National Energy Technology Laboratory in Morgantown, WV. The measurement procedures were validated by conducting benchmark tests with smooth stainless-steel tube and comparing the results with published correlations for Nusselt number (Nu) and friction factor. Over Reynolds numbers (ReD) ranging from 5 × 10 4 to 2.5 × 10 5 , measured Nu was within 5% of the Dittus-Boelter correlation and measured friction factors were within 5% of the McAdams correlation for smooth tube flow. The candidate pin fin patterned pipes were additively manufactured (AM) at the Oak Ridge National Laboratory. The pins were circular or elliptical in cross-section and were printed at a 30⁰ angle relative to the inner wall (to meet AM constraints). The pin arrangement was helical to promote enhanced heat transfer due to swirl flow. Pin length to diameter aspect ratio was 1.33, 2, and 8, while the pin diameter to tube diameter ratio was 0.188, 0.125, and 0.063. Tests were performed for ReD varying from 6.9 × 10 4 to 2.2 × 10 5 and at conditions equivalent to the low pressure (LP) outlet (8.69 MPa, 361 K) and the high pressure (HP) inlet (20.7 MPa, 350 K) of the low temperature recuperator (LTR) in an indirect sCO 2 cycle. The Wilson plot technique was utilized to measure the bulk heat transfer coefficients. For the best performing design (tube A, pin length to tube diameter ratio: 1.33, pin diameter to tube diameter ratio: 0.19), the local heat transfer coefficient increased by 136% relative to the Dittus-Boelter correlation at the LTR low pressure outlet and 194% at the LTR high pressure inlet. These correspond to a 282% and a 271% increase in the product of the heat transfer coefficient and surface area (adjusted for fin efficiency) product, respectively. Large pressure drops across the test articles were observed. For Tube Design A, the average friction factor, across the range of ReD considered, was significantly larger than the McAdams correlation at both the LTR LP outlet and the LTR HP inlet. A thermal performance factor was utilized to express the ratio of material required to build a finned heat exchanger relative to a finless heat exchanger with the same heat duty and pumping power. Tube Design A was estimated to decrease the required heat exchanger material by 13%.

Searle, Matthew↗

DOE OSTI2022

QUIC-URB and QUIC-fire extension to complex terrain: Development of a terrain-following coordinate system

Ensemble-based approaches to prescribed fire planning cannot be supported by CFD-based models like FIRETEC and WFDS because they are too computationally expensive and cannot leverage LES approaches like CAWFE and WRF-SFIRE because too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, typical for instance of the Western United States. In this paper, we describe the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. In particular, the paper presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captured the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. We report the model agreed well with measurements on the upwind side of the peak, but overestimated speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field was obtained in 6 s, making the solver appealing for coupled fire–atmosphere simulations. The Socorro Mountain was highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect was addressed in this work). Results highlight that future work should focus on the implementation of parameterizations of wake-eddies, similar to QUIC-URB’s building parameterizations, and on thermodynamic-driven flow.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗

DOE OSTI2021

Operation of Natural Gas Pipeline Networks With Storage Under Transient Flow Conditions

Here, we formulate a nonlinear optimal control problem for intraday operation of a natural gas pipeline network that includes storage reservoirs. The dynamics of compressible gas flow through pipes, compressors, reservoirs, and wells are considered. In particular, a reservoir is modeled as a rigid, hollow container that stores gas under isothermal conditions and uniform density, and a well is modeled as a vertical pipe. For each pipe, flow dynamics are described by a coupled partial differential equation (PDE) system in density and mass flux variables, with momentum dissipation modeled using the Darcy–Wiesbach friction approximation. Compressors are modeled as scaling up the pressure of gas between the inlet and outlet. The governing equations for all network components are spatially discretized and assembled into a nonlinear differential-algebraic equation (DAE) system, which synthesizes above-ground pipeline and subsurface reservoir dynamics into a single reduced-order model. We seek to maximize an objective function that quantifies economic profit and network efficiency subject to the flow equations and inequalities that represent operating limitations. The problem is solved using a primal–dual interior point solver, and the solutions are validated in computational experiments and simulations on several pipeline test networks to demonstrate the effectiveness of the proposed methodology.

03 NATURAL GAS↗

DOE OSTI2020

Radiant boiler for pressurized oxy-combustion and method of radiant trapping to control heat flux in high temperature particle-laden flows at elevated pressure

A boiler has a shell surrounding a vertical centerline. The shell defines an inner surface having an inner diameter and an inner length extending between an upper upstream end and a lower downstream end. The inner surface defines a hollow interior, the boiler having a pre-combustion zone, a combustion zone downstream from the pre-combustion zone, and a post-combustion zone downstream from the combustion zone. The shell is tapered outward along its length in at least a portion of the combustion zone. An oxidizer inlet is in fluid communication with the pre-combustion zone, and a fuel nozzle introduces fuel into the combustion zone. A tube assembly is mounted in the hollow interior of the shell for transferring heat to fluid flowing through the tube assembly. A flue duct is in fluid communication with the post-combustion zone for transporting flue gases from the hollow interior.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2024

Partially cracked ammonia mixture combustion in gas turbine combustors

Ammonia is gaining popularity for single fuel use in existing gas turbines, as it is carbonless, has higher hydrogen content per unit volume than liquid hydrogen, and facilitates easy storage and handling. However, the low flame speed of ammonia makes it prone to flame blow off and remains the chief concern for its use in power generation. Following previous studies that have shown improved combustion performance by partially cracking the fuel, in this computational study, efficacy of partially cracked ammonia as a replacement fuel for natural gas is evaluated. A high-pressure optically accessible combustor operated at ~10 bar pressure with non-vitiated heated air was used as the test platform. This combustor was equipped with a novel multifuel multi-tube micromixing injector (M3 injector). High resolution simulations of the flow fields both in the injector as well as the combustor were performed to understand the fuel mixing in the injector and in the inlet region of the combustor. These modeling results showed that the injector was highly effective for achieving mixture homogeneity at the entrance of the combustor for ammonia and natural gas injection, whereas, for hydrogen an earlier injection facilitating longer residence time would be beneficial. Also, these simulations show that while harnessing waste heat the ideal fuel dissociation fraction is between 0.4 and 0.55 for optimal operation.

Ammonia↗

DOE OSTI2020

Three-Dimensional CFD Analysis of Construction Design Alternatives for an ERDC Coastal and Hydraulics Laboratory Flow Accelerator

The U.S. Army Engineer Research and Development Center, ERDC, Coastal and Hydraulics Laboratory in Vicksburg Mississippi is planning to build a large flume with a 10-foot-wide flow channel. Three-dimensional Computational Fluid Dynamics (CFD) analysis was used in the design of the flume entry and flow accelerator sections that feed into the 10-foot-wide channel. The design is based on the plans of the existing fiberglass flume inlet at the Turner-Fairbank Highway Research Center (TFHRC) J. Sterling Jones Hydraulics Research Laboratory in McLean Virginia. The TFHRC flume inlet feeds water to a 6-foot-wide flume channel, and therefore a scale up factor of 10/6 in the width of the inlet section of the TFHRC flume would meet the width requirement of the ERDC flume channel. The CFD analysis was performed using STAR-CCM+ commercial CFD software. Additional analysis was done with the open-source CFD software, OpenFOAM, to cross-check the analysis of one of the cases with a code-to-code comparison.

42 ENGINEERING↗

DOE OSTI2024

Comparison of Sub-Ppm Instrument Response Suggests Higher Detection Limits Could Be Used to Quantify Methane Emissions from Oil and Gas Infrastructure

Quantifying and controlling fugitive methane emissions from oil and gas facilities remains essential for addressing climate goals, but the costs associated with monitoring millions of production sites remain prohibitively expensive. Current thinking, supported by measurement and simple dispersion modelling, assumes single-digit parts-per-million instrumentation is required. To investigate instrument response, the inlets of three trace-methane (sub-ppm) analyzers were collocated on a facility designed to release gas of known composition at known flow rates between 0.4 and 5.2 kg CH 4 h –1 from simulated oil and gas infrastructure. Methane mixing ratios were measured by each instrument at 1 Hertz resolution over nine hours. While mixing ratios reported by a cavity ring-down spectrometer (CRDS)-based instrument were on average 10.0 ppm (range 1.8 to 83 ppm), a mid-infrared laser absorption spectroscopy (MIRA)-based instrument reported short-lived mixing ratios far larger than expected (range 1.8 to 779 ppm) with a similar nine-hour average to the CRDS (10.1 ppm). We suggest the peaks detected by the MIRA are likely caused by a micrometeorological phenomenon, where vortex shedding has resulted in heterogeneous methane plumes which only the MIRA can observe. Further analysis suggests an instrument like the MIRA (an optical-cavity-based instrument with cavity size ≤10 cm3 measuring at ≥2 Hz with air flow rates in the order of ≤0.3 slpm at distances of ≤20 m from the source) but with a higher detection limit (25 ppm) could detect enough of the high-concentration events to generate representative 20 min-average methane mixing ratios. Even though development of a lower-cost, high-precision, high-accuracy instrument with a 25 ppm detection threshold remains a significant problem, this has implications for the use of instrumentation with higher detection thresholds, resulting in the reduction in cost to measure methane emissions and providing a mechanism for the widespread deployment of effective leak detection and repair programs for all oil and gas infrastructure.

03 NATURAL GAS↗

DOE OSTI2022

Scalable Simulation of Pressure Gradient-Driven Transport of Rarefied Gases in Complex Permeable Media Using Lattice Boltzmann Method

Accurate representations of slip and transitional flow regimes present a challenge in the simulation of rarefied gas flow in confined systems with complex geometries. In these regimes, continuum-based formulations may not capture the physics correctly. This work considers a regularized multi-relaxation time lattice Boltzmann (LB) method with mixed Maxwellian diffusive and halfway bounce-back wall boundary treatments to capture flow at high Kn. The simulation results are validated against atomistic simulation results from the literature. We examine the convergence behavior of LB for confined systems as a function of inlet and outlet treatments, complexity of the geometry, and magnitude of pressure gradient and show that convergence is sensitive to all three. The inlet and outlet boundary treatments considered in this work include periodic, pressure, and a generalized periodic boundary condition. Compared to periodic and pressure treatments, simulations of complex domains using a generalized boundary treatment conserve mass but require more iterations to converge. Convergence behavior in complex domains improves at higher magnitudes of pressure gradient across the computational domain, and lowering the porosity deteriorates the convergence behavior for complex domains.

42 ENGINEERING↗

DOE OSTI2021

Additively Manufactured Compliant Hybrid Gas Thrust Bearing for Supercritical Carbon Dioxide Turbomachinery: Experimental Evaluation and Fluid–Structure Model Predictions

This paper presents rotating test results and advances an analytical predictive fluid–structure model for a new type of gas-lubricated thrust bearing fabricated using direct metal laser melting. The bearing concept in this study is a compliant hybrid gas thrust bearing that uses external pressurization to increase load carrying capacity, where the testing campaign in this study was only focused on steady-state static performance. The need for the bearing concept comes from enabling highly efficient supercritical carbon dioxide (sCO2) turbomachinery by replacing oil-lubricated bearings with process gas lubrication. Leveraging the process gas of the turbomachine for bearing lubrication results in lowered bearing power loss, simplified mechanical design, and allows for novel oil-free hermetic drivetrains resulting in an efficient emission-free system. The new concept utilizes hydrostatic pressurization on individual tilting pads flexibly mounted with hermetic squeeze film dampers (HSFDs). This paper focuses on rotating tests of a 173 mm outer diameter gas thrust bearing in air up to 10 krpm and hydrostatic inlet pressures to 365 psi (2.52 MPa). The influence of thrust runner speed and bearing inlet pressure on force deflection characteristics and load carrying capability of the gas film were experimentally evaluated. This work also advances a predictive fluid–structure thrust bearing model using an isothermal ideal-gas-based compressible Reynolds flow equation directly coupled to a lumped stiffness element possessing axial and rotational degrees-of-freedom. The rotating testing demonstrated load capability of 1816 lbs (8.1 kN), which equates to a thrust bearing unit load of 67 psi (0.46 MPa). Gas film force–deflection curves reveal a nonlinear relationship between thrust load and film clearance. Comparison of film thickness values with the predictive model show good agreement under high load and inlet pressure, however deviate as load and pressure decrease. Load capability was shown to increase with increasing hydrostatic inlet pressure, while the increase in thrust runner speed revealed a small decrease in load capacity.

Engineering↗

DOE OSTI2023

Mapping the performance of a versatile water-based condensation particle counter (vWCPC) with numerical simulation and experimental study

Accurate airborne aerosol instrumentation is required to determine the spatial distribution of ambient aerosol particles, particularly when dealing with the complex vertical profiles and horizontal variations of atmospheric aerosols. A versatile water-based condensation particle counter (vWCPC) has been developed to provide aerosol concentration measurements under various environments with the advantage of reducing the health and safety concerns associated with using butanol or other chemicals as the working fluid. However, the airborne deployment of vWCPCs is relatively limited due to the lack of characterization of vWCPC performance at reduced pressures. Given the complex combinations of operating parameters in vWCPCs, modeling studies have advantages in mapping vWCPC performance. In this work, we thoroughly investigated the performance of a laminar-flow vWCPC using COMSOL Multiphysics® simulation coupled with MATLAB™. We compared it against a modified vWCPC (vWCPC model 3789, TSI, Shoreview, MN, USA). Our simulation determined the performance of particle activation and droplet growth in the vWCPC growth tube, including the supersaturation, $D$ p,kel,0 (smallest size of particle that can be activated), $D$ p,kel,50 (particle size activated with 50 % efficiency) profile, and final growth particle size D d under wide operating temperatures, inlet pressures $P$ (30–101 kPa), and growth tube geometry (diameter $D$ and initiator length $L$ ini ). The effect of inlet pressure and conditioner temperature on vWCPC 3789 performance was also examined and compared with laboratory experiments. The COMSOL simulation result showed that increasing the temperature difference (Δ$T$) between conditioner temperature $T$ con and initiator $T$ ini will reduce $D$ p,kel,0 and the cut-off size $D$ p,kel,50 of the vWCPC. In addition, lowering the temperature midpoint ($T$ mid = $\frac{T_{con}+T_{ini}}{2}$ increases the supersaturation and slightly decreases the $D$ p,kel . The droplet size at the end of the growth tube is not significantly dependent on raising or lowering the temperature midpoint but significantly decreases at reduced inlet pressure, which indirectly alters the vWCPC empirical cut-off size. Our study shows that the current simulated growth tube geometry ($D$=6.3 mm and $L$ ini =30 mm) is an optimized choice for current vWCPC flow and temperature settings. The current simulation can more realistically represent the $D$ p,kel for 7 nm vWCPC and also achieved good agreement with the 2 nm setting. Using the new simulation approach, we provide an optimized operation setting for the 7 nm setting. This study will guide further vWCPC performance optimization for applications requiring precise particle detection and atmospheric aerosol monitoring.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Thermal-Hydraulic Performance of a Zeotropic Refrigerant in a Brazed Plate Evaporator

In this work we report the results of a preliminary investigation focused on evaluating the thermal-hydraulic performance of R455A (CO2, R32, R1234yf – 3/21.5/75.5) in a commercial brazed plate heat exchanger. Three key parameters were measured which were refrigerant flow maldistribution, overall refrigerant side heat transfer coefficient, and total refrigerant side pressure drop. These same parameters were also measured for R134a which served as a base-line refrigerant for comparison. The experimental test matrix was set up to investigate the influence of refrigerant mass flux (3 ≤ Gref ≤ 5 kg m-2 s-1), refrigerant inlet vapor quality (0.1 ≤ xin ≤ 0.4), and refrigerant exit superheat (5 ≤ Tsh ≤ 10 K) on the above mentioned thermal-hydraulic parameters. For these experimental conditions, on average, the overall refrigerant side heat transfer coefficients for R455A were 45.2% lower than those of R134a. Refrigerant flow maldistribution for R455A, on average, was observed to be 62.6 % higher than that measured for R134a, and the total pressure drop across the heat exchanger (minor and frictional losses) for R455A was on average 26.5% higher than that measured for R134a. These initial results set the stage for a comprehensive investigation focused on further investigating the thermal-hydraulic performance of zeotropic refrigerants in brazed plate heat exchangers.

Jajja, Saad [ORNL]↗

DOE OSTI2022

Combining complex flow manifold with three dimensional woven lattices as a thermal management unit

The present invention is directed to a manifold for directing cooling fluid and/or gas to a heat exchanger in a flow configuration designed to optimize heat transfer from the heat exchanger. The manifold can take many different forms such as a layered construction with distributed inlet paths, local outlet paths, a central collection changer and a path for fluid removal. The manifold can be formed from a metal, plastic, rubber, ceramic, or other heat resistant material known to or conceivable by one of skill in the art. The manifold can also be combined with any type of heat exchanger known to or conceivable by one of skill in the art to form a thermal management unit. To optimize overall properties such as low pressure drop, high heat transfer, and excellent temperature uniformity of the thermal management unit, the manifold can be graded, expanded and scaled as needed.

Hemker, Kevin J.↗

DOE OSTI2024

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗

DOE OSTI2020

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗