Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuation power 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.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The role of microstructural evolution during spark plasma sintering on the soft magnetic and electronic properties of a CoFe–Al2O3 soft magnetic composite

Abstract For transformers and inductors to meet the world’s growing demand for electrical power, more efficient soft magnetic materials with high saturation magnetic polarization and high electrical resistivity are needed. This work aimed at the development of a soft magnetic composite synthesized via spark plasma sintering with both high saturation magnetic polarization and high electrical resistivity for efficient soft magnetic cores. CoFe powder particles coated with an insulating layer of Al 2 O 3 were used as feedstock material to improve the electrical resistivity while retaining high saturation magnetic polarization. By maintaining a continuous non-magnetic Al 2 O 3 phase throughout the material, both a high saturation magnetic polarization, above 1.5 T, and high electrical resistivity, above 100 μΩ·m, were achieved. Through microstructural characterization of samples consolidated at various temperatures, the role of microstructural evolution on the magnetic and electronic properties of the composite was elucidated. Upon consolidation at relatively high temperature, the CoFe was to found plastically deform and flow into the Al 2 O 3 phase at the particle boundaries and this phenomenon was attributed to low resistivity in the composite. In contrast, at lower consolidation temperatures, perforation of the Al 2 O 3 phase was not observed and a high electrical resistivity was achieved, while maintaining a high magnetic polarization, ideal for more efficient soft magnetic materials for transformers and inductors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

SOFE2023 Presentation: "Tritium Safety and Sensitivity Analysis of Tritium Extraction eXperiment (TEX) PbLi Loop Using MELCOR-TMAP"

The efficient extraction of tritium from the breeder blanket (BB) in a fusion reactor is vital for maintaining the tritium economy for a sustainable and self-sufficient fusion power plant. The vacuum permeator (VP) design is one of several promising technology platforms for extracting tritium from the BB with the advantages of being capable of operating continuously at high temperatures with the aim of achieving high tritium extraction efficiencies. The Tritium Extraction eXperiment (TEX), under construction at the Safety and Tritium Applied Research (STAR) Facility at Idaho National Laboratory (INL), is designed to test the viability of VP technology by introducing atomic hydrogen (H, D, and T) from a source permeator (SP) into a forced convection lead lithium eutectic (PbLi) loop and determining the hydrogen extraction rates from the vanadium VP. Presented here is a computational analysis of the hydrogen inventory and permeation behavior in TEX using the hydrogen transport and accident analysis code MELCOR-TMAP developed at INL. The results of a sensitivity analysis based on MELCOR-TMAP simulations of hydrogen extraction efficiencies from the TEX test section using varying PbLi flow rates, hydrogen source rates from the SP, and material property uncertainties are also reported. Operation scenarios for TEX in consideration of the analyses based on results from MELCOR-TMAP are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Pulsed power to control growth of silicon nanoparticles in low temperature flowing plasmas

Low-temperature plasmas have seen increasing use for synthesizing high-quality, mono-disperse nanoparticles (NPs). Recent work has highlighted that an important process in NP growth in plasmas is particle trapping—small, negatively charged nanoparticles become trapped by the positive electrostatic potential in the plasma, even if only momentarily charged. In this article, results are discussed from a computational investigation into how pulsing the power applied to an inductively coupled plasma (ICP) reactor may be used for controlling the size of NPs synthesized in the plasma. The model system is an ICP at 1 Torr to grow silicon NPs from an Ar/SiH 4 gas mixture. This system was simulated using a two-dimensional plasma hydrodynamics model coupled to a three-dimensional kinetic NP growth and trajectory tracking model. The effects of pulse frequency and pulse duty cycle are discussed. We identified separate regimes of pulsing where particles become trapped for one pulsed cycle, a few cycles, and many cycles—each having noticeable effects on particle size distributions. For the same average power, pulsing can produce a stronger trapping potential for particles when compared to continuous wave power, potentially increasing particle mono-dispersity. Pulsing may also offer a larger degree of control over particle size for the same average power. Experimental confirmation of predicted trends is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

FEED Study of Carbon Capture Inc. DAC and CarbonCure Utilization Using United States Steel's Gary Works Waste Plant

Direct Air Capture (DAC) has been proposed as a means of reducing atmospheric concentrations of CO2. While DAC has been evaluated through lab-scale, bench-scale, and small-scale pilot units, large-scale deployment has not been achieved. Feasibility studies are one tool to understand the potential design, operation, performance, and impact of commercial-scale DAC. This paper presents the results of a feasibility study for a passive DAC system deployed at >100,000 tpy scale in three different regions across the U.S. and awarded to Carbon Collect by the U.S. Department of Energy National Energy Technology Laboratory (DOE NETL). The MechanicalTree™technology has been designed and engineered by Carbon Collect based on the concept initially developed at Arizona State University Center for Negative Carbon Emissions. It uses a tower of stacked, sorbentcontaining disks supported by a lifting mechanism, exposed to the air to capture CO2 during the adsorption phase. The disks are then lowered into a regeneration chamber for vacuum and steam regeneration to produce CO2 product during the desorption step. The modular tree design allows large installations with repeatable, mass-manufactured units connected to common utilities such as steam supply, vacuum, and CO2 processing for compression and geologic storage. The passive DAC system eliminates the equipment and energy of forced air fans by using natural air circulation to contact the sorbent with the CO2 in the air. Because of this, the performance is dependent on the wind speed in addition to the temperature and relative humidity. Performance and flow rate fluctuations were incorporated into equipment and facility design with considerations for turndown to 10% of maximum rated flow to allow operation in all seasons. The feasibility study was undertaken to evaluate the technology in different regions and climates. Three locations were selected for this study representing different climates: Alabama (hot and humid), California (hot and arid), and Wyoming (continental). For each region, the adsorption/desorption cycle was optimized including tuning the heat integration and cycle timings. The performance of individual trees was then scaled to the full facility with the same design of more than 20,000 trees common between all regions. The installations had expected average CO2 capture rates of between 330,000 and 485,000 tonnes of CO2 per year depending on the climate. The initial engineering design for the facility at each location was performed with cost and performance estimates for the trees, carbon purification and compression, and the balance of plant. To supply thermal and electrical energy for the facility, carbon-free or low-carbon power must be considered. Options for low-carbon, continuous thermal and electrical energy were considered. The best-performing option from those considered was identified to be an electrically-heated molten salt energy storage system, powered by an on-site photovoltaic field. During molten salt thermal discharge, steam is produced from heat exchange with the molten salt and used to generate power in a steam turbine as well as steam for regenerating the DAC carbon trees. The thermal and electrical supply and analysis is presented in the context of low-carbon power for carbon removal. Modelling and economics of transportation and geologic storage of the CO2 is considered and presented for each location. The results of the feasibility study are incorporated into the presented techno-economic and life-cycle assessments. This work is intended to provide an understanding of the performance, cost, and impact of capturing CO2 at the commercial scale and the impact of climate and regional siting on the considered passive DAC system.

42 ENGINEERING↗
DOE OSTI2022

Environmental Flow Requirements from FERC Licenses Across the US

Environmental flow requirements included in Federal Energy Regulatory Commission (FERC) hydropower licenses are important for balancing natural properties and benefits of river ecosystems (e.g., healthy species, recreation, water supply, flood control) supporting hydropower production. In some cases, environmental flow requirements may limit operational flexibility given current operational schemes and make a hydropower plant less able to provide power to the electric grid on demand. Hydropower plants may gain some flexibility as hydropower scheduling time periods are made to be more responsive to the short-term needs of an energy grid increasingly reliant on intermittent renewables. However, many flow requirements focus on the daily, monthly, or seasonal flow fluctuations which matches the time scale of most paradigms linking flow alterations to the health of river ecosystems. This dataset seeks to provide a greater understanding of how flexibility in environmental requirements can be leveraged to create positive outcomes for both the power system and the environment. It contains information on environmental flow requirements from the Protection, Mitigation, and Enhancement section of 50 randomly selected FERC licenses: 25 issued from 1998-2013 that were also included in the ORNL Mitigation Database (Schramm et al. 2015) and 25 issued from 2014-present. The information on environmental flow requirements was extracted from the PM&E section of 50 randomly selected FERC licenses: 25 issued from 1998-2013 that were also included in the ORNL Mitigation Database (Schramm et al. 2015) and 25 issued from 2014-present. The flow requirements were then categorized into flow augmentation categories based on whether the license stated a specific water management purpose for the given requirement called augmentation categories (i.e., fisheries or habitat, recreation or boating, industry, and general or unspecified; Table B). Requirements were also grouped into flow type categories (e.g., minimum flow rate, maximum flow rate, ramping rate). Additional information related to flow requirements such as the augmentation time-period and whether the flow rate was continuous (i.e., condition must be present at all-times) or instantaneous (i.e., condition present at a point in time) was also extracted from the licenses. Some licenses had specific flow requirements based on whether the project was in a wet, dry, or normal water year. If that information was presented in the license, it was also included in the data set. The location within the project was noted, hereafter, zone, in the dataset for flow requirements relating to specific areas of hydropower projects (Dam, Powerhouse, Bypass Reach). Maximum discharge capacities of hydropower facilities were also extracted from both the Existing Hydropower Assets (EHA) data set and the National Inventory of Dams (NID) databases. Each facility was coded with project identification codes from the EHA dataset to facilitate cross-referencing between datasets.

13 HYDRO ENERGY↗
DOE OSTI2024

High Performance Heat Pipe Power Transient Testing at SPHERE Facility

Microreactors are being researched, designed, and built at Idaho National Laboratory (INL). Microreactors are small reactors defined at less than 20MW of power. These reactor concepts are also being looked at throughout industry for various applications. An important aspect of these reactor designs is economic feasibility i.e. lower overnight capital cost. The driving factors for implementing microreactors are quick setup and takedown, minimal operators, and the ability to manufacture them readily and to fit in mid-sized containers for transport. A specific area of research to aid in successful integration of these factors within the designs is passive heat removal of the core’s thermal power. Interest in heat pipes to achieve this passive heat removal has been shown across multiple industry partners. Because of this interest, INL has developed a test facility to facilitate experimental tests for sodium filled heat pipes. INL has developed the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility to run experiments on high performance, sodium filled heat pipes. As mentioned above, heat pipes are passive heat transfer devices. Radially, heat pipes are broken up into an outer wall, a small annular gap, a wick structure, and a centerline gap. They function by utilizing latent heat transfer. Heat pipes are traditionally separated into three regions, an evaporator (heat input), an adiabatic region, and finally a condenser region (heat removal). As heat is being applied to the evaporator, the working fluid undergoes a phase change to a vapor. This phase change causes a differential pressure across the axial length of the pipe driving flow down the center gap of the heat pipe. The vapor flows down past the adiabatic region to the condenser where the heat is removed. This heat removal forces the working fluid to phase change back to a liquid. The wick structure is then utilized to drive the flow back towards the evaporator by capillary forces. This backflow is aided by the annular gap. Because this heat transfer mechanism functions with latent heat transfer, the heat pipe is close to isothermal down the axial length. Heat pipes can operate under a wide range of working fluids. Considerations for these working fluids are primarily driven by operating temperatures amongst other important factors based around overall performance. Sodium filled heat pipes operate from 450°C up to 900°C. This temperature range works well for the current microreactor designs. In conjunction with this experimental capability, INL has developed a modeling software to simulate heat pipe physics within reactor cores. This modeling software is called Sockeye and functions under the established INL Multiphysics Object Oriented Simulation Environment (MOOSE). SPHERE also supports Sockeye development by providing the modeling team with experimental data on an array of setups and operating parameters to support validation efforts. A power transient experiment was performed utilizing the SPHERE facility to continue to aid with Sockeye development. The testing followed a proposed test plan to ramp up and down the temperature of the heat pipe. Sockeye models steady state heat pipe operation with high accuracy, the data provided by the power transient testing aims to assist with the validation efforts and further enhance transient modeling capability of the tool [2].

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Leveraging Optimal Sparse Sensor Placement to Aggregate a Network of Digital Twins for Nuclear Subsystems

Nuclear power plants (NPPs) require continuous monitoring of various systems, structures, and components to ensure safe and efficient operations. The critical safety testing of new fuel compositions and the analysis of the effects of power transients on core temperatures can be achieved through modeling and simulations. They capture the dynamics of the physical phenomenon associated with failure modes and facilitate the creation of digital twins (DTs). Accurate reconstruction of fields of interest (e.g., temperature, pressure, velocity) from sensor measurements is crucial to establish a two-way communication between physical experiments and models. Sensor placement is highly constrained in most nuclear subsystems due to challenging operating conditions and inherent spatial limitations. This study develops optimized data-driven sensor placements for full-field reconstruction within reactor and steam generator subsystems of NPPs. Optimized constrained sensors reconstruct field of interest within a tri-structural isotropic (TRISO) fuel irradiation experiment, a lumped parameter model of a nuclear fuel test rod and a steam generator. The optimization procedure leverages reduced-order models of flow physics to provide a highly accurate full-field reconstruction of responses of interest, noise-induced uncertainty quantification and physically feasible sensor locations. Accurate sensor-based reconstructions establish a foundation for the digital twinning of subsystems, culminating in a comprehensive DT aggregate of an NPP.

42 ENGINEERING↗
DOE OSTI2021

General relativistic MHD simulations of non-thermal flaring in Sagittarius A*

Sgr A* exhibits regular variability in its multiwavelength emission, including daily X-ray flares and roughly continuous near-infrared (NIR) flickering. The origin of this variability is still ambiguous since both inverse Compton and synchrotron emission are possible radiative mechanisms. The underlying particle distributions are also not well constrained, particularly the non-thermal contribution. In this work, we employ the GPU-accelerated general relativistic magnetohydrodynamics code H-AMR to perform a study of flare flux distributions, including the effect of particle acceleration for the first time in high-resolution 3D simulations of Sgr A*. For the particle acceleration, we use the general relativistic ray-tracing code bhoss to perform the radiative transfer, assuming a hybrid thermal+non-thermal electron energy distribution. We extract ~60 h light curves in the sub-millimetre, NIR and X-ray wavebands, and compare the power spectra and the cumulative flux distributions of the light curves to statistical descriptions for Sgr A* flares. Our results indicate that non-thermal populations of electrons arising from turbulence-driven reconnection in weakly magnetized accretion flows lead to moderate NIR and X-ray flares and reasonably describe the X-ray flux distribution while fulfilling multiwavelength flux constraints. These models exhibit high rms per cent amplitudes, $\gtrsim 150{{\ \rm per\ cent}}$ both in the NIR and the X-rays, with changes in the accretion rate driving the 230 GHz flux variability, in agreement with Sgr A* observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

In Situ Wind and Turbulence Measurements in a Field of Full-Size Parabolic Trough Solar Collectors

Concentrated Solar Power (CSP) is a promising method for using Solar power for electricity generation with thermal energy storage. One of the primary drivers of structural design costs of CSP collector structures is wind loading. To date, the design of these structures has relied on data from wind tunnels that do not adequately capture the dynamic effects observed at scale. NREL initiated a field measurement campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs, providing insights on structural dynamic response, and generating a first-of-a-kind, comprehensive, high-resolution wind-loading dataset available for validating simulations of wind loading on collector structures. The measurements at NSO consist of Sonic anemometers on masts at different heights to characterize the incoming flow and conditions within and above the trough field, in combination with a Doppler Lidar scanning the horizontal plane above the troughs. The wind measurements at NSO have been continuously operating since October 2021 and provide a year-long dataset characterizing wind and turbulence conditions. The structural load measurements start in November 2022 and will complement the wind measurements. In this poster, we present first results of the wind measurement campaign by highlighting days with different atmospheric flow conditions. We identified three main factors altering the flow over the parabolic troughs: Wind speed, wind direction, and the angle of the parabolic troughs. The highest loads are expected when the wind blows perpendicular to the trough rows. In this case, the first rows experience the highest loads and block the subsequent rows, creating conditions with decreased wind speed and enhanced turbulence within the trough field. Also, turbulent length scales change. These affected conditions produce unique load cases on the structures that will be captured by the load measurements. If the wind blows along the trough rows, wind and turbulence conditions are less influenced by the troughs.

atmospheric turbulence↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2022

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Heat transfer from glass melt to cold cap: Computational fluid dynamics study of cavities beneath cold cap

Efficient glass production depends on the continuous supply of heat from the glass melt to the floating layer of batch, or cold cap. Computational fluid dynamics (CFD) are employed to investigate the formation and behavior of gas cavities that form beneath the batch by gases released from the collapsing primary foam bubbles, ascending secondary bubbles, and in the case of forced bubbling, from the rising bubbling gas. The gas phase fraction, temperature, and velocity distributions below the cold cap are used to calculate local and average heat transfer rates as a function of the bubbling rate. It is shown that the thickness of the cavities is nearly independent of the cold cap shape and the amount of foam evolved during batch conversion. Furthermore, it is ~7 mm and up to ~15 mm for the cases without and with forced bubbling used to promote circulation within the melt, respectively. Using computed velocity and temperature profiles, the melting rate of the simulated high-level nuclear waste glass batch was estimated to increase with the bubbling rate to the power of ~0.3 to 0.9, depending on the flow pattern. The simulation results are in good agreement with experimental data from laboratory- and pilot-scale melter tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Microfluidic Liquid-Liquid Extraction Chip with Integrated Raman Sensors (Phase I and Phase II Final Report)

As demand for electricity continues to increase worldwide, the world’s nuclear power generating capacity will continue to grow. Nuclear power is the most environmentally benign way of producing electricity on a large scale. The long-term successful use of nuclear power, however, is critically dependent upon adequate and safe processing and disposal of spent nuclear fuels. A very important feature of nuclear energy is that spent fuels can be reprocessed to recover fissile and fertile materials that can then be used as fresh fuel for nuclear power plants. The DOE-NE Fuel Cycle Research and Development (FCR&D) is currently developing nuclear material reprocessing technologies. In a typical nuclear fuel reprocessing system, centrifugal contactors are used as liquid-liquid extraction devices where two immiscible liquids are mixed at high speeds using a rotor, which creates a fine dispersion of droplets of organic phase in an aqueous phase that contains the analyte. Understanding the extraction efficiency at these contactors through modeling and simulation is important in the development of reprocessing technologies and the optimization of current technologies such as the PUREX process. The outcome of this program will be a liquid-liquid microfluidic flow cell chip with embedded spectroscopic sensors for the analysis of extracted analytes. This device will be useful in the optimization of new fuel reprocessing schemes as well as existing reprocessing processes by providing a microfluidic modeling platform to optimize extraction parameters. The Phase I and Phase II work developed a microfluidic chip design that allows integration of Raman and absorption fiber optic probes and successfully demonstrated the feasibility of using Raman and absorption probes to detect key analytes that partition into the aqueous and organic phases.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Produced Water and Waste Heat-aided Blowdown Water Treatment: Using Chemical and Energy Synergisms for Value Creation

The project objective was to develop a cooling blowdown water (BDW) treatment process utilizing produced water (PW) and low-grade heat to maximize water reuse and saleable by-product generation while reducing chemical and energy footprints of the treatment. The proposed treatment process consists of mixing, softening, organics and suspended solids removal, reverse osmosis (RO), thermal desalination, and brine electrolysis. BDW samples collected from a local coal-fired power plant and PW samples from two shale gas production wells were used in this study. Each treatment unit was first designed and tested to quantify its treatment efficiency, and its chemical and energy requirements. In addition, a process model was developed and model simulations were conducted based on the experimental results and literature data to optimize the treatment process. A techno-economic analysis was conducted to quantify chemical and energy savings as well as production of 10-lb brine as a saleable product. With the field-collected BDW and PW samples, mixing experiments determined a volumetric mixing ratio 10:1 (BDW:PW) resulted in the best performance of multivalent ions removal and largest chemical savings for softening. Softening of the BDW/PW mixtures using alkaline chemicals (Na 2 CO 3 and NaOH) achieved 95%-100% removal of scaling-forming cations (Ca, Mg, Fe, Ba, Sr) and 60% of silicon, and 10% of total organic carbon (TOC). The mixing and softening treatments yielded an effluent with total dissolved solids (TDS) concentration of 23 g/L. Activated carbon (AC) filtration removed TOC to a low level (< 3 mg/L) and further removed remaining scale-forming divalent metals and silica from the softened water. The AC filtration resulted in a slight reduction of TDS from 23 g/L to 20 g/L, leaving behind only mostly monovalent ions (i.e., sodium and chloride) in the filtered water. These pretreatments yielded a feed water that met the criteria of the downstream reverse osmosis (RO) to prevent membrane fouling. A cross-flow RO system was used to further concentrate the TDS of the AC effluent. Various factors including TDS, pH, and applied pressure were examined and optimal conditions were determined for the co-treatment process. An integrated process consisting of mixing, softening, AC filtration and RO was used to treat a continuous flow (0.25 – 1.2 L/min, or 0.07 – 0.32 gpm) and successfully generated RO permeate as product water (TDS < 0.5 g/L) for reuse in cooling operation, and a concentrate (TDS ~ 45 g/L) to be further treated in a thermal desalination unit. These flow rates meet the FOA’s criterion of 0.01 – 1 gpm. Overall, the co-treatment of BDW/PW allowed shorter ramp-up time compared to treatment of BDW alone. It resulted in 40% and 55% savings of Na 2 CO 3(s) and NaOH, respectively, compared to treating the BDW and PW individually for the same level of softening. The co-treatment also resulted in a 29% energy saving compared to treatment of BDW only for the level of TDS concentration. A thermal desalination system was designed using CFD simulations and manufactured in the WVU Innovation Hub for further treatment of the RO concentrate to generate 10-lb brine. The system has a design flow rate of 2 gpm and has been successfully tested. A bench-scale brine electrolysis system was developed for on-site generation of chlorine/hypochlorite (Cl 2 /OCl - ) and caustic soda (NaOH) as useful chemicals for the co-treatment process. Using salt solutions (0.5 M and 1 M), the system achieved faradaic efficiencies of 93%-97% and 70%-77% for caustic soda and chlorine/hypochlorite generation, respectively. An economic analysis showed that the electricity costs for on-site generation of these chemicals were significantly lower than the chemical prices offered by suppliers. An industrial-scale process model consisting of mixing, softening, AC filtration, RO, thermal desalination, and brine electrolysis was developed using the Aspen Plus V9 in conjunction with Aspen Custom Modeler V9. The model serves as a solvable Aspen Plus model and as basis to form the costing infrastructure. In addition, techno-economic analysis considering capital, operating, and transportation costs was conducted. An optimization solution showed that produced water for mixing is still advantageous in low quantities. The optimum solution approaches a leveled cost of water (LCW) of 2 $/m 3 which becomes cost competitive with nominal water treatment prices.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Relativistic gas accretion onto supermassive black hole binaries from inspiral through merger

Accreting supermassive black hole binaries are powerful multimessenger sources emitting both gravitational and electromagnetic (EM) radiation. Understanding the accretion dynamics of these systems and predicting their distinctive EM signals is crucial to informing and guiding upcoming efforts aimed at detecting gravitational waves produced by these binaries. To this end, accurate numerical modeling is required to describe both the spacetime and the magnetized gas around the black holes. In this paper, we present two key advances in this field of research. First, we have developed a novel 3D general relativistic magnetohydrodynamics (GRMHD) framework that combines multiple numerical codes to simulate the inspiral and merger of supermassive black hole binaries starting from realistic initial data and running all the way through merger. Throughout the evolution, we adopt a simple but functional prescription to account for gas cooling through photon emission. Next, we have applied our new computational method to follow the time evolution of a circular, equal-mass, nonspinning black hole binary for ∼200 orbits, starting from a separation of 20⁢𝑟 𝑔 and reaching the postmerger evolutionary stage of the system. We have shown how mass continues to flow toward the binary even after the binary “decouples” from its surrounding disk, but the accretion rate onto the black holes diminishes. We have identified how the minidisks orbiting each black hole are slowly drained and eventually dissolve as the binary compresses. We confirm previous findings that the system’s luminosity decreases by a factor of a few during inspiral; however, we observe an abrupt increase by ∼50% in this quantity at the time of merger, likely accompanied by an equally abrupt change in spectrum. Lastly, we have demonstrated that during the inspiral, fluid ram pressure regulates the fraction of the magnetic flux transported to the binary that attaches to the black holes’ horizons.

Accretion disk & black-hole plasma↗
DOE OSTI2025

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗
DOE OSTI2020

Deuterium retention characteristics in Li film by coating and during flowing liquid Li limiter operation in experimental advanced superconducting tokamak

Lithium (Li) is a promising low-Z material for particle recycling and impurity control to improve plasma performance in fusion devices. In the experimental advanced superconducting tokamak (EAST), Li coating has become a routine method for wall conditioning, and a flowing liquid Li (FLiLi) limiter has been successfully tested several times. Deuterium retention characteristics in the Li film coated on the international thermonuclear experimental reactor-like tungsten divertor and FLiLi during plasma discharges, which is important for the utilization of Li in future fusion devices, were investigated in EAST. It is found that the absorption of the fuel particles by Li coatings decreases gradually, and recycling gradually increases over a series of discharges. The maximum net amount of deuterium retained reached ~0.8 g, corresponding to 12% deuterium in the Li. Furthermore, this corresponds to a whole day's worth of shots with a total of 87 plasma discharges and a total of ~640s plasma time after 11.75 g Li was deposited. Compared to the Li coating, it is shown that FLiLi continuously traps fuel particles and achieves a higher deuterium retention ratio over both the short- and long-term, leading to lower recycling. Meanwhile, it is also observed that the fuel particle retention ratio increases when FLiLi is closer to the plasma. This result is likely due to higher plasma heating power and limiter temperature, which cause an increased Li efflux from FLiLi due to a higher limiter temperature, which then redeposits on the other plasma-facing surfaces and increases fuel particle absorption. It is estimated that ≥ 80% of the retained D particles are captured by the continual renewal of the Li redeposition film during the FLiLi operation. This investigation would also be useful for D/T retention in future fusion devices if Li is used as a plasma-facing component.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

DS-GL: Advancing Graph Learning via Harnessing the Power of Nature within Dynamic Systems

With the rapid digitization of the world, an increasing number of real-world applications are turning to nonEuclidean data, modeled as graphs. Due to their intrinsic high complexity and irregularity, learning from graph data demands tremendous computational power. Recently, CMOS-compatible Ising machines, i.e., dynamic systems composed of CMOS components, have emerged as a new approach that harnesses the inherent power of natural annealing within dynamic systems to efficiently resolve binary optimization problems and have been adopted for traditional graph computation, such as max-cut. However, when performing complex Graph Learning (GL) tasks, Ising machines face significant hurdles: (i) they are inherently binary and thus ill-suited for real-valued problems; (ii) their expensive all-to-all coupling network that guarantees effective natural annealing poses daunting scalability concerns. To address these challenges, this paper proposes a nature-powered graph learning framework dubbed DS-GL, which is the first effort to transform the process of solving graph learning problems into the natural annealing process within a parameterized dynamic system embodied as a CMOS chip. To tackle the two major hurdles, DS-GL first augments the Ising machine architecture to modify the self-reaction term of its Hamiltonian function from linear to quadratic, effectively serving as an energy regulator. This adjustment maintains the system’s original physical interpretation while enabling it to process continuous, real-valued data. Second, to address the scaling issue, DS-GL further upgrades the real-valued dense Ising machine by decomposing it into a mesh-based multi-PE dynamic system that supports efficient distributed spatial-temporal co-annealing across different PEs through sparse interconnects. By exploiting the inherent sparsity and component structures in real-world graphs, DS-GL is able to map complex graph learning tasks onto the scalable dynamic system while maintaining high accuracy. Evaluations with three diverse GL applications across six real-world datasets, including traffic flow and COVID-19 prediction, show that DS-GL can deliver from 102× to 106× speedups and 500× energy reduction over Graph Neural Networks on GPUs, with 5% - 20% accuracy enhancement.

Song, Ruibing↗