Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

GeN-ROM—An OpenFOAM®-based multiphysics reduced-order modeling framework for the analysis of Molten Salt Reactors

This work presents a projection-based multiphysics Model Order Reduction (MOR) framework for the analysis of nuclear systems and its application to parametric simulations of Molten Salt Reactors (MSR). The framework, named GeN-ROM, is developed using OpenFOAM® and employs a Proper Orthogonal Decomposition aided Reduced-Basis technique (POD-RB). It can be used to reduce steady-state and transient multiphysics problems involving parametric fluid dynamics, heat exchange, and neutronics phenomena. For the treatment of structural elements in the hydraulic systems, a porous medium approach has been adopted. The reduction process is data-driven and snapshot information is extracted via POD to learn the solution manifold and to build global spatial basis functions. At the data collection phase, GeN-ROM makes use of the solvers available in GeN-Foam, a similarly OpenFOAM®-based multiphysics framework developed for the analysis of nuclear reactors. The global bases are used both to approximate the solution fields and to project the full-order equations onto lower-dimensional subspaces, thus considerably reducing the number of unknowns in a numerical system. This reduction leads to significant computational speedups, which is ideal for multi-query applications such as uncertainty quantification or design optimization. The developed tool has been tested using a 2D multiphysics model of the Molten Salt Fast Reactor (MSFR) with steady-state and transient scenarios, with speedups on the order of 10 – 10 5 .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Elucidating the Impact of Flow on Material-Sensitive Critical Heat Flux and Boiling Heat Transfer Coefficients: An Experimental Study with Various Materials

In this work, steady-state internal flow boiling critical heat flux (CHF) experiments were carried out at 10°C inlet subcooling and atmospheric pressure (84 kPa at an altitude of about 5300 feet) with mass flux ranging from 200 to 2,000 kg/m 2 -s on seven tube materials, including the accident-tolerant fuel cladding candidate FeCrAl (C26M) alloy. At a relatively low flow rate, the sensitivity of CHF to different materials was appreciable (the average difference among the tested materials was 19% at 200 kg/m 2 -s). As the absolute CHF differences stayed constant for higher flow rates, the increased flow rate led to reduced relative CHF differences among the tested materials (average relative difference of 15% at 1,000 kg/m 2 -s). FeCrAl alloy claddings that were oxidized under simulated Light Water Reactor (LWR) conditions yielded CHF comparable to as-received FeCrAl alloy cladding, with less than 1.8% difference. Experimental CHF results imply marginal material sensitivity impacts on the departure from nucleate boiling ratio in steady-state nuclear reactor cores. Yet, for accident progression analysis, material sensitivity may hold importance, as boiling behavior variation at low mass flux is appreciable. Macroscopic material parameters (surface wettability, surface roughness, and thermal properties) do not provide a suitable explanation for observed material-sensitive CHFs. Yet, surface roughness and thermal effusivity revealed power and linear relations, respectively, to the nucleate boiling heat transfer coefficients at low mass flux. At high mass flux, material sensitivity on nucleate boiling heat transfer significantly decreased due to the increased shear stress on boiling surface overwhelming material-related boiling implications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Perspective on the oil-dominated gas hydrate plugging conceptual picture as applied to transient Shut-In/Restart

Gas hydrate formation is rarely an issue during continuous production conditions. Unplanned and/or extended shut-downs comprise the majority of scenarios where overlapping concern and uncertainty for gas hydrate blockages exist. To date, most experimental studies and theoretical analysis of gas hydrate formation derive from the steady-state conceptual picture. Yet, significant distinctions induced by phase separation and prolonged cooling suggest an altered conceptual picture is required for transient shut-in/restart operations. To address this outstanding issue, this work presents a conceptual picture for the base case of a non-surface active hydrocarbon undergoing cold restart without enacted mitigation efforts. First, phase separation upon shut-in leads to the localization of gas hydrate ingredients due to fluids density differences and pipeline inclination. Cooling to the seafloor temperature, then, induces gas hydrate film formation at the hydrocarbon interface and maximizes subcooling, resulting in rapid and confined gas hydrate growth occurring upon restart. Furthermore, this ultimately leads to plugging dominated by the deposition of porous, low conversion gas hydrate masses. Direct visual observations and results from Raman spectroscopy, hydrate contact angle and film growth, micromechanical force apparatus, rocking cell, and flow loop studies support the mechanisms outlined within each step. Similar to the influence provided by the original Colorado School of Mines steady-state conceptual picture that was introduced about 15 years ago, this transient counterpart builds a foundation that could inspire new directions of hydrate research and analysis in areas where uncertainty abounds.

02 PETROLEUM↗
DOE OSTI2023

Highly Coupled Heterobicycle‐Fused Porphyrin Dimers: Excitonic Coupling and Charge Separation with Coordinated Fullerene, C 60

Abstract Porphyrin dimers have been widely explored and studied owing to their importance in photosynthetic systems. A vast variety of dimers linked by different groups and at different angles have been synthesized and studied; however, the means by which to synthesize rigidly fused porphyrins with direct conjugation of the chromophores remains limited. Such a class of porphyrins may possess interesting properties that unconjugated or stacked dimers may not exhibit. In this study, bisbenzimidazole‐fused porphyrin dimers and their mono‐ and bis−zinc derivatives are synthesized and characterized. As a consequence of excitonic coupling, these dimers exhibit a split Soret band irrespective of the metal ion in the porphyrin cavity. Steady‐state fluorescence and excitation spectra followed by femtosecond transient absorption spectral studies of the heterometallated dimer, (free‐base and zinc porphyrin) reveals the occurrence of efficient singlet‐singlet energy transfer (>95 % efficiency and rate constant >10 12 s −1 ) within the dyad. Further, donor‐acceptor conjugates were formed by metal‐ligand axial coordination of phenyl imidazole functionalized C 60 and were characterized by a variety of physicochemical techniques. Excited state charge separation from both singlet and triplet excited states of ZnP in the conjugates has been established. The lifetime of the final charge‐separated state was in the 30–40 μs range revealing charge stabilization. Interestingly, no charge separation in the conjugate derived from the heterometallated dimer was observed wherein excitation transfer dominated the process. The present study brings out the importance of the rigid π‐spacer connecting porphyrin dimers in governing the energy and electron transfer events when coupled with an electron acceptor.

Chemistry↗
DOE OSTI2022

Mechanistic pathways and role of oxygen in oxidative coupling of methane derived from transient kinetic studies

Oxidative coupling of methane (OCM) is a promising industrial process to upgrade natural gas to high value chemicals. In this study, Temporal Analysis of Products (TAP) and steady-state experiments were conducted to distinguish how the composition of surface and gas phase oxygen influence mechanistic details of the selective conversion of CH 4 to C 2 H 4 over the Mn-Na 2 WO 4 /SiO 2 catalyst. The results from TAP studies indicate that methane activation on this catalyst proceeds predominantly via a short-lived, transient surface oxygen species and there is a competition for this species to form either CO or methyl radicals on the surface. This active species has a total lifetime of 3 s and is identified to have a dioxygen (e.g. O 2 2- or O 2 -) form. We show that the concentration of the transient surface oxygen species significantly impacts the OCM performance. Oxygen attributed to the catalyst lattice (in a singular form e.g., O - ), is found to activate methane to a lesser degree, but exclusively forms CO 2 . Evidence for these surface pathways for methyl radical, CO and CO 2 formation identified by TAP are also validated through steady-state experiments. Finally, by distinguishing different catalyst oxygen species and their role in selective/nonselective pathways, important screening criteria have been identified for the advancement of superior catalyst formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Learning the factors controlling mineral dissolution in three-dimensional fracture networks: applications in geologic carbon sequestration

We perform a set of high-fidelity simulations of geochemical reactions within three-dimensional discrete fracture networks (DFN) and use various machine learning techniques to determine the primary factors controlling mineral dissolution. The DFN are partially filled with quartz that gradually dissolves until quasi-steady state conditions are reached. At this point, we measure the quartz remaining in each fracture within the domain as our primary quantity of interest. We observe that a primary sub-network of fractures exists, where the quartz has been fully dissolved out. This reduction in resistance to flow leads to increased flow channelization and reduced solute travel times. However, depending on the DFN topology and the rate of dissolution, we observe substantial variability in the volume of quartz remaining within fractures outside of the primary subnetwork. This variability indicates an interplay between the fracture network structure and geochemical reactions. We characterize the features controlling these processes by developing a machine learning framework to extract their relevant impact. Specifically, we use a combination of high-fidelity simulations with a graph-based approach to study geochemical reactive transport in a complex fracture network to determine the key features that control dissolution. We consider topological, geometric and hydrological features of the fracture network to predict the remaining quartz in quasi-steady state. We found that the dissolution reaction rate constant of quartz and the distance to the primary sub-network in the fracture network are the two most important features controlling the amount of quartz remaining. This study is a first step towards characterizing the parameters that control carbon mineralization using an approach with integrates computational physics and machine learning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Barium stars as tracers of s -process nucleosynthesis in AGB stars

Barium (Ba) stars help to verify asymptotic giant branch (AGB) star nucleosynthesis models since they experienced pollution from an AGB binary companion and thus their spectra carry the signatures of the slow neutron capture process (s process). For a large number (180) of Ba stars, we searched for AGB stellar models that match the observed abundance patterns. We aim to uncover any systematic deviations of the sample abundances from the predictions of the nucleosynthesis models. We employed three machine learning algorithms as classifiers: a Random Forest method, developed for this work, and the two classifiers used in our previous study. Compared to that work, we also expanded our observational sample with 11 Ba stars available in the supersolar metallicity range. We studied the statistical behaviour of the different s-process elements in the observational sample to investigate if the AGB models systematically under- or overpredict the abundances observed in the Ba stars and show the results in the form of violin plots of the residuals between spectroscopic abundances and model predictions. We inspected the correlations between the observed [Fe/H], the s-process elemental abundances, and the residuals. We employed the [Zr/Fe] and [Nb/Fe] abundances as a thermometer to constrain the operational temperature that rules the production of these elements in the sample stars, assuming a steady-state s process. We also investigated the mass distribution of the identified polluter AGB stars and the behaviour of the δ parameter, which describes the fraction of accreted AGB material relative to the Ba star envelope. We find a significant trend in the residuals that implies an underproduction of the elements just after the first s-process peak (Nb, Mo, and Ru) in the models relative to the observations. This may originate from a neutron-capture process (e.g. the intermediate neutron-capture process, i process) not yet included in the AGB models of metallicity from solar to roughly 1/5 solar, corresponding to the range of the Ba stars. Correlations are found between the residuals of these peculiar elements, suggesting a common origin for the deviations from the models. In addition, there is a weak metallicity dependence of the residuals of these elements. The s-process temperatures derived with the [Zr/Fe] – [Nb/Fe] thermometer have an unrealistic value for the majority of our stars. The most likely explanation is that at least a fraction of these elements are not produced in a steady-state s process, and instead may be due to processes not included in the AGB models. The mass distribution of the identified models confirms that our sample of Ba stars was polluted by low-mass AGB stars (< 4 M ⊙ ). Most of the matching AGB models require low accreted mass, but a few systems with high accreted mass are needed to explain the observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Progress in the development and understanding of a high poloidal-beta tokamak operating scenario for an attractive fusion pilot plant

Abstract The high poloidal-beta ( $$\beta _{\textrm{P}}$$ β P ) regime was first proposed as a high bootstrap current scenario for a steady-state fusion pilot plant (FPP) in the 1990s (Kikuchi in Nucl Fusion 30:265, 1990). Since then, there have been many theoretical, modeling, and experimental research activities on this topic. A joint DIII-D/EAST research team began exploring the high- $$\beta _{\textrm{P}}$$ β P regime in 2013, focusing on addressing the needs of attractive FPP design by taking advantage of the extensive diagnostic set and sophisticated plasma control system on DIII-D and the well-developed integrated modeling capability at General Atomics. The ultimate goal is to demonstrate such a scenario on EAST with truly long pulse and metal wall compatibility. This paper summarizes the highlights of the research results on DIII-D by the joint team in the past decade. Experimental evidence and modeling analysis show the high- $$\beta _{\textrm{P}}$$ β P scenario has great advantages in addressing key needs for an attractive FPP design, such as high-energy confinement quality at low rotation, excellent core-edge integration, high line-averaged density above the Greenwald limit, low disruption risk, and high bootstrap current fraction for steady-state operation. This provides a relatively safe and economical option to base an FPP design on that will lead to commercial fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Techno-economic analysis of renewable energy generation at the South Pole

Transitioning from fossil-fuel power generation to renewable energy generation and energy storage in remote locations has the potential to reduce both carbon emissions and cost. Here, this study presents a techno-economic analysis for implementation of a hybrid renewable energy system at the South Pole in Antarctica, which currently hosts several high-energy physics experiments with nontrivial power needs. A tailored model of resource availability and economics for solar photovoltaics, wind turbine generators, lithium-ion energy storage, and long-duration energy storage at this site is explored in different combinations with and without existing diesel energy generation. The Renewable Energy Integration and Optimization (REopt) platform is used to determine the optimal system component sizing and the associated system economics and environmental benefit. We find that the least-cost system includes all three energy generation sources and lithium-ion energy storage. For an example steady-state load of 170 kW, this hybrid system includes 180 kW-DC of photovoltaic panels, 570 kW of wind turbines, and a 3.4 MWh lithium-ion battery energy storage system. This system reduces diesel consumption by 95% compared to an all -diesel configuration, resulting in approximately 1200 metric tons of carbon footprint avoided annually. Over the course of a 15-year analysis period the reduced diesel usage leads to a net savings of 57 million United States dollars, with a time to payback of approximately two years. All the scenarios modeled show that the transition to renewables is highly cost effective under the unique economics and constraints of this extremely remote site.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Creating an Advanced Sensor Network to calculate real-time, mass-weighted flue gas composition and air heater leakage of a coal-fired utility boiler under dynamic operating conditions

Utilization of renewable energy sources to minimize the environmental impact of energy production has changed the way utility boilers operate, requiring frequent load cycling between full load and partial loads as low as 30%. Dynamic operation of coal-fired utility boilers significantly reduces boiler efficiency when compared to steady state at full load. Data-driven plant optimization has shown success with coal-fired utility boilers under dynamic operating conditions. The purpose of this work was to create an Advanced Sensor Network (ASN) to provide more extensive real-time data to inform dynamic plant optimization of Net Unit Heat Rate (NUHR). The ASN consists of gas sampling grids in the convective pass of the boiler and downstream of the air heater. These sampling grids allow for quantification of spatial variation of flue gas within the boiler and calculation of mass-weighted composition of flue gas through the combination of composition, velocity, and temperature measurements. The comparison of O 2 between the inlet and outlet of the air heater is used to calculate air leakage in real time. Flue gas composition and air heater leakage are both important factors in boiler efficiency and NUHR. Further, the results of this work support the value of mass-weighted averages for determining flue gas composition accurately. The measurements from the ASN show increased composition stratification during dynamic operation, with an average standard deviation 38% higher than observed during steady-state operation. Air heater leakage was also observed to increase from 2.8% to 5.1% following a load change. Prior to the installation of the ASN, these data would not have been available for dynamic control. These real-time data will be leveraged to calculate and optimize for NUHR during dynamic operation in future work.

42 ENGINEERING↗
DOE OSTI2020

Gas Production from Hot Water Circulation through Hydraulic Fractures in Methane Hydrate-Bearing Sediments: THC-Coupled Simulation of Production Mechanisms

Methane hydrates, widely found in permafrost and deep marine sediments, have great potential as a future energy 6 source. Conventional production schemes perform poorly for challenging hydrate reservoirs with low permeability. We propose an 7 efficient production scheme by combining hydraulic fracturing from horizontal wells and hot water circulation through fractures. A 8 fully coupled thermo-hydro-chemical (THC) model is developed to simulate the key physical processes during gas production from 9 a hydrate reservoir representative of typical geological settings in Shenhu, South China Sea. We found that the gas production 10 process has two distinct stages divided by thermal breakthrough: a relatively short prebreakthrough stage and a postbreakthrough 11 stage yielding stable gas production. Heat advection along and near the hydraulic fracture dominates the prebreakthrough stage, 12 whereas conduction-driven thermal recovery in the volume around fractures dominates the postbreakthrough stage. We identified 13 that the steady-state injection temperature has a strong effect on the performance of the proposed scheme while the fluid mass 14 circulation rate has a moderate impact beyond a threshold. The proposed scheme proves to be efficient and robust over a range of 15 reservoir conditions with respect to initial hydrate saturation and intrinsic permeability, including their spatial heterogeneities, 16 thereby offering a promising solution for challenging reservoir conditions.

gas hydrate, fracture flow, gas production, numeri↗
DOE OSTI2023

Modeling of standing wave RF cavities for tracking through multi-pass energy recovery linac

Short bunches, high current and multiple linac pass are all characteristics of Energy Recovery Linacs (ERLs), which may result in collective effects. They in turn, may affect the beam, degrading its quality, or even yield to instabilities causing a beam loss. To study and mitigate these effects one needs a numerical simulation code, that can take into account both the collective effects, as well as, particular ERL features, such as a multi-turn design that does not reach a steady state or the multiple passages of the beam through Radio-Frequency (RF) cavities at different energies. CODAL [1], a code developed by SOLEIL in collaboration with IJCLab, enables such studies. It is a 6 dimensional (6D) tracking code applying 'kicks' based on the integration of the local Hamiltonian for each element of the lattice. It is also capable of simulating space charge, wakefields and coherent synchrotron radiation. However, to correctly take into account the ERL dynamics, an upgrade had to be made to include the effect of a standing wave RF cavity in 6D. In this paper, we will concentrate on the implementation and benchmarking (with DESY’s tracking code ASTRA [2]) of both the longitudinal and the transverse models (by J.B. Rosenzweig and L. Serafini [3]), which we use to carry out tracking of fully analytical 6D RF cavity.

Bogacz, A.↗
DOE OSTI2022

Fast Quasi-Static Time-Series Simulation for Accurate PV Inverter Semiconductor Fatigue Analysis with a Long-Term Solar Profile

Power system simulations with long-term data typically have large time steps varying from one second to a few minutes. However, for PV inverter semiconductors, the minimum thermal stress cycle occurs over the fundamental grid frequency (50 or 60 Hz). This requires the time step of the fatigue simulation to be around 100 µs. This small time step requires long computation times to process yearly power production profiles. This paper proposes a fast fatigue simulation for inverter semiconductors using the quasi-static time series (QSTS) simulation concept. The proposed simulation calculates the steady state of the semiconductor junction temperature by using a Fast Fourier Transform (FFT). The small thermal cycling during a switching period and even over the fundamental waveform is disregarded to further accelerate the simulation speed. The resulting time step of the fatigue simulation is 15 minutes, which is consistent with the solar dataset. The error of the proposed simulation is 0.16% compared to the fatigue simulation results using the complete thermal stress profile. A PV inverter that responds to a Transactive Energy System (TES) is simulated to demonstrate the use of the proposed fatigue simulation. The proposed simulation has the potential to co-simulate with system level simulation tools that also adopt the QSTS concept.

Liu, Yunting↗
DOE OSTI2023

A Perturbative Solution for Nonlinear Stratified Upwelling over a Frictional Slope

Abstract A perturbative solution of simplified primitive equations for nonlinear weakly stratified upwelling over a frictional slope is found that resolves the vertical structure of velocity fields and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The solution uses assumptions consistent with the model proposed by Lentz and Chapman, including a steady-state, constant cross-shore density gradient, no alongshore gradients, laterally inviscid, and consideration of cross-shore advection of alongshore momentum. The solution resolves the vertical structure of velocity fields (including subsurface maxima of compensational flow, not resolved by Lentz and Chapman) and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The dynamics are similar to Lentz and Chapman; bottom stress balances alongshore wind stress in a homogeneous density ocean and is replaced by nonlinear cross-shore transport of alongshore momentum as the Burger number (S=αN/f, whereα,N, andfare the bottom slope, buoyancy frequency, Coriolis frequency, respectively) increases. When the solution uses the empirical relation between cross-shore and vertical density gradients proposed by Lentz and Chapman, vorticity conservation is not satisfied and the nonlinear momentum transport estimated by the solution linearly increases withS, asymptotically matching Lentz and Chapman forS< 1. When the solution conserves interior potential vorticity, the momentum transport is proportional toS 2 forS< 1 and is in better agreement with numerical simulations.

Oceanography↗
DOE OSTI2026

An Updated Picture of Presolar History from Short-lived Radioactive Isotopes and Inferences on the Birth of the Sun

We examine the origin of the short-lived radionuclides (SLRs; defined as having half-lives between 0.1 and 100 Ma) present in the early solar system (ESS) by investigating how predictions of their abundances in the interstellar medium (ISM) from steady-state equilibrium relate to their ESS values. For this, we take into account the non-negligible time t iso elapsed between the isolation of the presolar molecular cloud and the formation of the ESS, during which the SLRs decayed freely. We also consider the alternative scenario in which the presolar molecular cloud remained partially mixed with the ISM, with a mixing timescale t mix . We find that the ESS abundances of 107 Pd and 182 Hf produced by slow neutron captures (s-process) and of 53 Mn and 60 Fe produced by explosive nucleosynthesis can be consistently explained within these scenarios. Their required tiso is 9–12 Ma and their required tmix is 11–14 Ma (with one potential exception of t mix = 38 Ma), depending on galactic uncertainties, such as the galactic star formation history and efficiency and the star-to-gas mass ratio. Another s -process SLR, 205 Pb, has a more uncertain ESS value and falls within only some of these time values. The same applies to the SLRs produced by the p -process ( 92 Nb and 146 Sm), depending on the latter’s half-life. In agreement with previous studies, we find that the ESS abundances of the rapid neutron-capture isotopes ( 129 I, 244 Pu, and 247 Cm) and of the most short-lived radionuclides ( 26 Al, 36 Cl, and 41 Ca) cannot be explained by assuming steady-state equilibrium in the ISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Analysis of Design Alternatives of Actively Cooled RF Window for MPEX

The Materials Plasma Exposure eXperiment (MPEX) has been designed as a linear plasma divertor simulator in order to address plasma material interaction (PMI) science for next-generation fusion devices. It will have the capability to test neutron irradiated samples with plasma fluxes of greater than 10 24 m -2 s -1 . It is expected to operate steady state for up to 10 6 s to consider PMI affects through reactor end of life. The conceptual design of MPEX was completed in 2019, with preliminary design having begun in 2020. The plasma source for MPEX is a helicon antenna, where the energized helical antenna sits outside of the vacuum in order to minimize impurities in the plasma. It is expected to receive up to 200 kW of continuous power, and so the antenna and the window must be actively cooled. The water-cooled copper antenna has been operated at full power on the Proto-MPEX device (which is a test facility to demonstrate the plasma source and heating systems). The water-cooled window, however, is a novel component that must meet numerous competing requirements. It requires a low dielectric loss to allow the Radio Frequency (RF) power to create the plasma within the vacuum boundary. It must be structurally robust to handle the significant heat flux from the plasma and any heat from dielectric coupling. It must be compatible with the coolant (preferably water). It requires a vacuum seal that minimizes impurities into the plasma and does not compromise the structural integrity of the window. Two window designs have been tested. Results from these tests, where temperatures are measured and heat fluxes inferred from infrared camera data, have been correlated with thermal-structural simulations. When these simulations are extrapolated to the full power steady-state heat fluxes that are expected in MPEX, the designs do not appear to have the necessary structural robustness. This study explores design alternatives for the MPEX helicon antenna window, presents analysis results for several of the alternatives, and shows a viable solution that satisfies the requirements for MPEX operation.

Tipton Jr, Joseph↗
DOE OSTI2023

Heat Transfer Characteristics of Particle and Air Flow Through Additively Manufactured Lattice Frame Material Based on Octet-Shape Topology

Particle-to-supercritical carbon dioxide (sCO 2 ) heat exchanger is a critical component in next-generation concentrating solar power (CSP) plants. The inherently low heat transfer between falling particles and sCO 2 imposes a challenge toward economic justification of levelized cost of electricity produced through solar energy. Introduction of integrated porous media with the walls bounding particle flow has the potential to enhance the overall particle-to-sCO 2 heat exchanger performance. Here, this paper presents an experimental study on heat transfer characterization of additively manufactured lattice frame material based on Octet-shaped unit cell with particles and air as working fluids. The lattice structures were additively manufactured in stainless steel (SS) 316L and SS420 (with 40% bronze infiltration) via Binder jetting process, where the lattice porosities were varied between 0.75 and 0.9. The mean particle diameters were varied from 266 μm to 966 μm. The effective thermal conductivity and averaged heat transfer coefficient were determined through steady-state experiments. It was found that the presence of lattice enhances the effective thermal conductivity by 2–4 times when compared to packed bed of particles alone. Furthermore, for gravity-assisted particle flow through lattice panel, significantly high convective heat transfer coefficients ranging from 200 W/m2K to 400 W/m2K were obtained for the range of particle diameters tested. The superior thermal transport properties of Octet-shape-based lattice frame for particle flow makes it a very promising candidate for particle-to-sCO 2 heat exchanger for CSP application.

14 SOLAR ENERGY↗
DOE OSTI2022

Assessment of Core Physics Characteristics of Extended Enrichment and Higher Burnup LWR Fuels using the Polaris/PARCS Two-Step Approach. Vol. I: PWR Fuel

Nuclear fuel with extended enrichment (235U enrichment within 5-8 wt%) is one of the evolutionary changes that have been pursued in recent years by commercial light water reactor operators and fuel vendors to improve the fuel cycle economy and operation performance of a nuclear plant. This work assesses the performance of the Polaris/PARCS two-step approach in core physics modeling of the pressurized water reactor cores with extended enrichment fuel, referred to as “LEU+” in this report. A representative LEU+ core with a 24-month fuel cycle developed by Southern Nuclear Company (SNC) was modeled using this two-step approach. A representative LEU core with an 18-month fuel cycle was also modeled to provide a reference for the LEU+ core. As expected, significantly more burnable poison absorbers were used in the LEU+ core to accomodate its higher fuel enrichment. Nine different fuel assembly types were modeled using Polaris for each core to generate the assembly cross sections, which were then processed by GenPMAX to prepare the cross-section data for PARCS. The average specific powers of each fuel batch in each core were derived from VERA results and higher specific powers in fresh assemblies were found in the LEU core due to its less total uranium loading included in the VERA LEU model, given that the total core power was assumed to be the same for both cores. PARCS models were developed to simulate the steady-state operations of both cores. PARCS results on the LEU+ core were first compared with the VERA results for verification purpose; good agreements were seen in soluble boron and burnup distribution results, indicating that the Polaris/PARCS modeling and simulation were correctly implemented. Core physics parameters calculated by PARCS, at zero power physics tests, beginning of cycle (BOC), and end of cycle conditions (EOC), were compared between the LEU+ core and the LEU core, including soluble boron concentration, burnup distributions, assembly and pin power peaking factors, fuel temperature reactivity coefficients, moderator temperature and density reactivity coefficients, control rod worth, and shut down margin. The main differences in PARCS results between the LEU+ and the LEU cores are summarized below: 1)The critical boron concentrations were found to be much higher in the LEU+ core than in the LEU core (1582 vs. 1335 ppm for peak values). 2)Higher assembly radial power peaking factors (1.4 vs. 1.3 for peak values), 2D pin peaking factors (1.53 vs. 1.42 for peak values), and 3D pin peaking factors (1.89 vs. 1.81 for peak values) were found in the LEU+ core than in the LEU core. 3)Significantly higher reactivity coefficients of moderator temperature (and density) were found in LEU+ than in LEU.4)Significantly lower control rod worth at EOC were found in LEU+ than LEU for all but one control banks.5)Significantly lower shut down margins were found in the LEU+ core than in the LEU core. The Polaris/GenPMAX/PARCS code suite was found to be capable of modeling the LEU+ PWR core for steady-state operations and no unexpected results in core physics parameters were observed, in spite of that a) several bugs in PARCS were identified and workarounds were used; b) several features were found lacking in the current version of PARCS that would be useful for core modeling. A list of requests for bug fixes and feature upgrades for PARCS originated from this work were transmitted to the code developers. The assessments on the performance of the Polaris/PARCS two-step approach in core modeling for boiling water reactor with LEU+ fuel is ongoing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗