Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow Conditions”

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.

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

At least 199 records · Page 11

Summary of DASSH-CFD Inter-Assembly Heat Transfer Comparison

Fast reactor core design requires integration of detailed power distributions obtained from neutronics models and thermal fluids calculations to determine temperature distributions and flow conditions. During the design phase, calculations may be repeated as specifications change and computational expense must be balanced with solution accuracy, creating an opportunity for the use of intermediate-fidelity codes. To that end, the Ducted Assembly Steady-State Heat Transfer Software (DASSH) has been developed at Argonne National Laboratory to perform full-core, steady state coolant and fuel pin temperature calculations in reactors consisting of hexagonal, ducted assemblies. DASSH is a first-order, finite volume subchannel code based on methodologies initially introduced in the legacy thermal hydraulics codes SUPERENERGY-2 and SE2-ANL. It can work coupled with the Argonne Reactor Computational (ARC) code suite to obtain neutron and gamma power distributions based on neutron transport solutions and features many improvements and new capabilities compared to its predecessors. DASSH calculates temperatures by balancing energy between subchannels and relies on correlations to account for flow distribution and coolant mixing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Experimental Analysis of Advanced Turbine for Supercritical CO2 Power Cycle

An advanced first stage high pressure turbine blade for a supercritical CO2 power cycle is tested in the Big Rig for Aerothermal Stationary Turbine Analysis (BRASTA) annular cascade at the Purdue Experimental Turbine Aerothermal Laboratory (PETAL) high pressure blow down facility alongside a baseline blade for comparison of aerodynamic performance. Both geometries are tested simultaneously using a novel off-axis rig design to allow for the blade geometries to be scaled up to achieve higher Reynolds numbers. The off-axis design necessitates the use of discrete sectors of airfoils, which are designed and additively manufactured in house using an mSLA printer. Printing the blades allows for unique routing of passage for blade surface static pressure taps that contour to the blade shape instead of requiring a straight line view from surface tap to egress, allowing all instrumentation for 15%, 50%, and 85% span to exist in the same passage. A flow conditioning gauze [1] is placed upstream of the blade passages to impart pressure, Mach, and swirl profiles to mimic the rotor relative frame inlet conditions to the blade. Downstream, a sonic valve is used to alter the backpressure to achieve different pressure ratios for different blowdown setpoints. The design of the rig, sonic wheel, and instrumentation are discussed along with the manufacturing and GD&T of the blades with respect to the build-up of this novel experimental apparatus. Initial commissioning data is obtained from inlet total pressure and temperature rakes, blade static pressure taps, and exit total pressure rakes. Oil visualization is performed using a silicone oil mixture containing titanium dioxide and pigmented for contrast of the blade suction side surface with the hub and shroud endwalls. The viscosity of the oil is tuned so that the oil does not fully thin and blow away during the time of the test, and endoscopic cameras are placed in the rig for live monitoring and recording. The preliminary oil visualization and rake profiles are presented, showcasing the operability of this additively manufactured off-axis flow path design.

Tuite, Logan [Purdue University]↗

Modeling Flow in a Microreactor Core: From High-Fidelity CFD to Subchannel Analysis

The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Surrogates for Valve-Controlled Pipe Flow: Accelerating Nuclear Reactor Design

Neural surrogate models are developed to replace expensive steady-state RANS CFD simulations for valve-controlled pipe flow in nuclear reactor design. Using parametric CFD data generated with MOOSE Pronghorn across a range of valve geometry and flow conditions, three approaches are compared: a POD-based reduced-order model, a structured UNet on a cylindrical grid, and unstructured models (DeepONet and BiStride MeshGraphNet) on nondimensionalized point clouds. POD achieves the highest accuracy (99%) with fast inference but requires storing all solution snapshots, while the DeepONet and BSMS-GNN both achieve ~89% accuracy at sub-second inference, with the BSMS-GNN offering superior geometric generalizability. These surrogates enable rapid ranking of candidate valve designs and can warm-start CFD solvers to accelerate convergence, supporting agentic design iteration on the Prometheus platform.

42 - ENGINEERING↗

Dynamic structural and microstructural responses of a metal–organic framework type material to carbon dioxide under dual gas flow and supercritical conditions

The structural and microstructural responses of a model metal–organic framework material, Ni(3-methyl-4,4'-bipyridine)[Ni(CN) 4 ] (Ni-BpyMe or PICNIC-21), to CO 2 adsorption and desorption are reported for in situ small-angle X-ray scattering and X-ray diffraction measurements under different gas pressure conditions for two technologically important cases. These conditions are single or dual gas flow (CO 2 with N 2 , CH 4 or H 2 at sub-critical CO 2 partial pressures and ambient temperatures) and supercritical CO 2 (with static pressures and temperatures adjusted to explore the gas, liquid and supercritical fluid regimes on the CO 2 phase diagram). The experimental results are compared with density functional theory calculations that seek to predict where CO 2 and other gas molecules are accommodated within the sorbent structure as a function of gas pressure conditions, and hence the degree of swelling and contraction in the associated structure spacings and void spaces. Furthermore, these predictions illustrate the insights that can be gained concerning how such sorbents can be designed or modified to optimize the desired gas sorption properties relevant to enhanced gas recovery or to addressing carbon dioxide reduction through carbon mitigation, or even direct air capture of CO 2 .

36 MATERIALS SCIENCE↗

Laser Measurements and Modeling of Shock Tunnel Freestream Velocity and Multispecies Thermal Nonequilibrium

Coherent anti-Stokes Raman scattering (CARS) and nitric oxide molecular tagging velocimetry (NO-MTV) were used to characterize the freestream in Sandia’s Hypersonic Shock Tunnel (HST) using a burst-mode laser operated at 100 kHz. Experiments were performed at nominal freestream velocities of 3 and 4 km/s using both air and N 2 test gas. The CARS diagnostic provides nonequilibrium characterization of the flow by measuring vibrational and rotational temperatures of N 2 and O 2 , which are compared to NO temperatures from separate laser absorption experiments. Simultaneous, collinear freestream velocities were measured using NO MTV along with pitot pressure measurements. Furthermore, this extensive freestream dataset is compared to nonequilibrium CFD that is capable of modeling species-specific, vibrational temperatures throughout the nozzle expansion. Significant nonequilibrium between vibrational and rotational temperatures was measured at each flow condition. N 2 exhibits the most nonequilibrium, followed by O 2 and NO. The CFD model captures this trend, although it consistently overpredicts N 2 and O 2 vibrational temperatures. At 3 km/s, the modeled NO nonequilibrium is underpredicted, whereas it is overpredicted at 4 km/s. Good agreement is seen between CFD and the velocity and rotational temperature measurements. Experiments with water added to the test gas yielded no discernable difference in vibrational relaxation.

Aerodynamics↗

Iron Single‐Atom Catalyzed N‐Alkylation of Amines with Alcohols via Solvent‐Free Borrowing Hydrogen Strategy

Abstract Industrial hydrogenation is a pivotal process in chemical synthesis. However, it has significant drawbacks, including high cost, safety risks associated with the use of molecular hydrogen gas, and substantial energy demands due to the need for elevated temperatures and pressures to achieve satisfactory yields. The borrowing hydrogen synthesis, which enables the transfer of hydrogen between molecules, offers a promising approach for green, one‐pot synthesis of industrially important chemicals and intermediates. Despite its potential, the broad application remains limited due to the reliance on toxic solvents, expensive noble metal catalysts, and the still restricted efficiency and substrate scope. In this study, the first solvent‐free strategy for theN‐alkylation of amines with alcohols is presented, employing an N‐doped graphene‐supported Fe single‐atom catalyst (Fe SA @N‐G; 1.06 wt.%). This approach achieves superior conversion and selectivity (up to 99%) along with record values for turnover number (TON, 1032.7) and turnover frequency (TOF, 413.1 h -1 ) for the coupling reaction of aniline with benzyl alcohol, surpassing all previously reported catalysts. DFT calculations, combined with experimental data, elucidated the reaction mechanism and identified the Fe 1 (III)‐N 4 active site participating in Fe‐H hydride transfer and containing two pyrrolic and two pyridinic nitrogens bound to the Fe center. The developed technology is further supported by the catalyst's excellent scalability, reusability, and performance under continuous‐flow conditions. Additionally, the exceptional efficiency of the single‐atom catalyst is demonstrated across more than 50 substrates, including reactions involving both aliphatic and aromatic amines with aliphatic and aromatic alcohols. The industrial applicability of this technology is validated through the synthesis of pharmaceutically relevant compounds, including stimulant drugs, antihistamines, and pharmaceutical intermediates.

Chemistry↗

Performance Comparisons for Artificially Propagated and Wild Pacific Lamprey Juveniles and Larvae

ABSTRACT Artificially propagated Pacific lamprey ( Entosphenus tridentatus ) are produced for restoration and for use in dam passage studies to reduce the demand for wild fish. Such uses require that animals are representative of their wild counterparts. Previous work indicated that this is true for Pacific lamprey larvae and juveniles reared in the hatchery with respect to the length of sustained swimming. However, more subtle differences in behaviour and performance that lamprey need to survive have not been assessed. In this study, artificially propagated and wild fish were compared in laboratory tests under no‐flow conditions to examine light avoidance, burrowing speed, burst swim speed, volitional routine swim speed and time to come to rest. Most larvae burrowed in less than a minute, and we found highly significant differences ( p < 0.001) between artificially propagated and wild larvae burrowing times, a critical escape behaviour. This could have implications for studies of larval entrainment at irrigation diversion canals or in turbine boils at dams. Interestingly, all of the wild juveniles tested came to rest quickly after introduction to the chamber (1.5 min), while artificially propagated lamprey swam robotically near the surface and 48% did not come to rest in the first 10 min (median time to rest = 9.5 min). In contrast, wild juveniles quickly (median = 1.47 min) sought areas near the tank bottom and attached strongly with their oral disc. Such behavioural differences could have important survival consequences for artificially propagated lamprey as they approach turbine intakes, bypass screens and irrigation diversion headgates. This study highlights the need to conduct behavioural assays that examine subtleties of fish behaviour that can be missed with traditional swim tunnel comparisons.

Frick, Kinsey [Fish Ecology Division, Northwest Fi↗

Oxygen consumption of sexually mature adult, first-feeding larval, and yearling Pacific Lampreys

Abstract Objective The Pacific Lamprey Entosphenus tridentatus is considered a first food of Native American tribes, such as the Confederated Tribes of the Umatilla Indian Reservation. Populations of this imperiled species have declined such that harvest for traditional use is limited. As a result, propagation and culture of Pacific Lampreys have been initiated to supply animals for restoration and recovery efforts. These efforts require transport and periodic holding of both adults and larvae under static (no-flow) conditions. Hence, guidelines for ensuring an adequate oxygen supply are needed for all life stages of this species. Methods We measured oxygen consumption rates of mature adults, first-feeding larvae, and yearlings under static conditions at 12–15°C. In addition, we recorded indicators of stress during hypoxia for adults and estimated routine respiration rates during and after larval feeding. Result Adults exhibited surprisingly high metabolic rates when at rest in static chambers. At 12°C, a single adult typically reduced dissolved oxygen levels in 10 L to <2 mg/L in just 90 min (oxygen consumption rates of 100–200 mg/kg/h). Adults often started to climb the walls of open chambers when dissolved oxygen levels approached 2 mg/L, allowing them to raise their branchiopores into air. Larvae remained quiescent, even when oxygen levels dropped below 1 mg/L, and costs of feeding increased routine respiration by 22–56%. Conclusion Based on these observations, we recommend that adult Pacific Lampreys always be transported in aerated containers and with access to air at the top of the tank. Although larvae exhibited hypoxia tolerance, care should be taken to ensure adequate oxygen availability, particularly during and immediately after feeding.

Moser, Mary↗

New vortex‐lift and tangential‐force models for HAWT aerodynamic load prediction

Abstract Horizontal axis wind turbines (HAWTs) experience three‐dimensional rotational and unsteady aerodynamic phenomena at the rotor blades sections. These highly unsteady three‐dimensional effects have a dramatic impact on the aerodynamic load distributions on the blades, in particular, when they occur at high angles of attack due to stall delay and dynamic stall. Unfortunately, there is no complete understanding of the flow physics yet at these unsteady 3D flow conditions, and hence, the existing published theoretical models are often incapable of modelling the impact on the turbine response realistically. The purpose of this paper is to provide an insight on the combined influence of the stall delay and dynamic stall on the blade load history of wind turbines in controlled and uncontrolled conditions. New dynamic stall vortex and nonlinear tangential force coefficient modules, which integrally take into account the three dimensional rotational effect, are also proposed in this paper. This module along with the unsteady influence of turbulent wind speed and tower shadow is implemented in a blade element momentum (BEM) model to estimate the aerodynamic loads on a rotating blade more accurately. This work presents an important step to help modelling the combined influence of the stall delay and dynamic stall on the load history of the rotating wind turbine blades which is vital to have lighter turbine blades and improved wind turbine design systems.

Elgammi, Moutaz↗

Control methods for mitigating flow oscillations in a supercritical CO 2 recompression closed Brayton cycle

A dynamic model of a 10 MWe supercritical CO 2 (sCO 2 ) recompression closed Brayton power generation cycle is used to investigate control methods for mitigating oscillations in flow conditions and power demand load, especially as the load ramp rate is increased and the load setpoint is more closely tracked. The focus is on control of the main compressor inlet temperature (MCIT) since oscillations are prompted by the strong nonlinearities in sCO 2 fluid properties near the critical point and by changes in system inventory for efficient load management. In this study, the oscillatory behavior during load ramping is first shown and then two potential control solutions for substantially reducing oscillations are given: 1) the use of the main compressor inlet guide vanes (IGV) and 2) bypassing a portion of the sCO 2 around the cooler to maintain the MCIT. These control methods reduce the effect of cycle feedback from interactions between the MCIT and load. When implementing either control method, a load setpoint ramp of 7.5%/min was maintainable while achieving closer load setpoint tracking with limited oscillations. Finally, the impact of sensor noise on the control was examined since derivative action was used in the IGV and cooler sCO 2 bypass control. Finally, while some effect was evident, signal noise was not problematic to the control.

24 POWER TRANSMISSION AND DISTRIBUTION↗

CFD modeling of turbulent air flow in self-heated gyroid TPMS structures: Thermal-hydraulic performance and validation

The application of mathematically derived geometries, such as triply periodic minimal surface (TPMS) lattices, has garnered significant interest across various fields, including the nuclear sector, due to their superior thermal-hydraulic characteristics for heat transfer compared to traditional plain or finned tubes. Here, this study validates a computational fluid dynamics (CFD) model, evaluates different turbulence models and CFD model settings, and performs uncertainty quantification to provide a comprehensive analysis. Despite extensive research on CFD modeling of TPMS lattices, such as gyroid and diamond geometries, there is a notable lack of publicly available literature providing comprehensive details on numerical analysis aspects, including convergence and methodological best practices. This study embarks on a benchmark analysis of a gyroid geometry to evaluate its thermal-hydraulic performance under turbulent flow conditions and scrutinize various CFD model configurations. The main contributions of this work include validating the CFD model, assessing and comparing different turbulence models, and enhancing pressure drop and temperature prediction capabilities. The results aim to support the development of methodologies needed to benchmark and enhance numerical analysis techniques for TPMS lattices. This work seeks to complement the existing body of knowledge, support the development of TPMS reactor concepts, and improve best practices for CFD modeling of TPMS lattices, ultimately advancing methodologies to support future applications in this domain.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of convection heat transfer on multiscale rough superhydrophobic and liquid infused surfaces

Multiscale rough superhydrophobic or slippery liquid infused porous surfaces have gained much interest in recent years for their improved transport phenomena properties. While there have been several studies on drag reduction and condensation on non-wetting surfaces, convection heat transfer that is important in many thermal and thermochemical applications has not been addressed systematically. Here, this article utilizes a fractal description of rough surface topographies to develop analytical models for the Nusselt number and the thermal hydraulic factor for fluid flow and heat transfer inside a cylinder with non-wetting surfaces. For air-infused superhydrophobic surfaces, the model considers the dynamic stability of the air/fluid interface in the asperities. Using the analytical formulations and the stability criteria, systematic studies are presented on the effects of the fractal surface parameters, cylinder radius and Reynolds number on the convective heat transfer characteristics, from which surface texture design maps are developed for maximizing the convection heat transfer. It is shown that multiscale non-wetting surfaces are most effective in the range of lower Reynolds number and small cylinder radius for achieving the best convective heat transfer and thermal hydraulic performance. Applying the models to actual non-wetting surface topographies fabricated using electrodeposition and chemical etching, it is shown that contrary to prevailing notion, superhydrophobicity, characterized by the highest contact angles, does not always lead to the maximum convective heat transfer performance, and that under certain fluid flow conditions, hydrophobic surfaces may offer a greater thermal performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling

Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .

Close-clearance impellers↗

Rheological and pumping behaviour of 3D printable cementitious materials with varying aggregate content

3D printing with cementitious materials is an emerging technology and provides an active area of research. Despite this, there are only limited studies that examine the pumping behaviour of high yield stress fluids such as printable concrete. We present a study on the rheological behaviour of 3D printable cementitious materials with different aggregate-to-binder ratios (a/b) and analyse the pumping and printing characteristics with the measured rheological parameters. We observed that an increase in a/b resulted in a significant increase of plastic viscosity, and a nominal increase in yield stress and storage modulus. The variation in plastic viscosity and yield stress with increasing aggregate content was analysed with the Krieger-Dougherty and Chateau-Ovarlez-Trung model, respectively. Also, the rheological properties of the paste rich lubrication layer during pumping were determined with the help of a tribometer. The yield stress and viscous constant of the lubrication layer increase with the increase in a/b. The pumping pressure predictions were carried out by a suitable theoretical model by considering appropriate flow conditions for different mixtures. The variation in critical height of printed elements with an increase in a/b was correlated with the observed rheological parameters.

36 MATERIALS SCIENCE↗

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗

Modeling and mitigating airborne pathogen risk factors in school buses

Computational fluid dynamics (CFD) models were developed to simulate the impact of different ventilation scenarios on airborne exposure risks in a 72-passenger school bus. Here, scenarios and factors that were investigated included a moving vs. stationary bus, impacts of a heating unit within the bus, and impacts of alternative ventilation scenarios with different combinations of openings (e.g., windows, door, emergency hatch). Results of the simulations showed that when the bus was stationary, use of the heater increased receptor concentrations unless there was another opening. When the bus was moving, simulations with at least two sets of openings separated from each other in the forward and aft directions produced a through-flow condition that reduced concentrations via dilution from outside air by a factor of ten or more. A single opening in a moving bus generally increased concentrations throughout the cabin due to increased mixing with minimal ventilation. The cumulative exposure risk (time-averaged concentrations) was found to be inversely correlated to the air exchange rate. Stationary and moving-bus scenarios that yielded above ~20 air changes per hour resulted in the lowest cumulative exposures. Recommendations from this study were implemented in new safety and operating procedures by the Albuquerque Public Schools Transportation Center.

60 APPLIED LIFE SCIENCES↗

A fast reduced model for a shell-and-tube based latent heat thermal energy storage heat exchanger and its application for cost optimal design by nonlinear programming

Numerical simulation of latent heat thermal energy storage (LHTES) systems plays a fundamental role in studying the physical process and guiding the engineering design. Discretization of the PDEs describing the nonlinear solidification/melting process of phase change materials (PCMs) leads to a large-scale complex dynamics system, where the system behavior depends on a set of parameters. In a design setting, repeated model evaluations are required over the set of parameters results in significant computational burden. In this paper, an explicit analytic solution was built for the propagation of the solidification front in a cylindrical coordinate. The analytic solution approach is further employed to develop a low computational reduced model (RM) as a module for a shell-and-tube based LHTES heat exchanger. The levelized Cost of Energy (LCOE) is used as a design metric and the RM model is used to apply system-level constraints in the nonlinear programming formulation that facilitates efficient global optimal design of the PCM properties, flow conditions and tube geometries. The use of LCOE as the design metric prevents over design of the heat transfer rate and also establishes a fair ground for evaluation of different thermal storage technologies and their integrated applications with other systems. Optimal results showed that a higher effectiveness results in a higher LCOE; the velocity of the HTF and the length of the channel are highly correlated with each other; both larger PCM latent energy and conductivity result in lower LCOE.

25 ENERGY STORAGE↗