Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fractional operators”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Characterization of W production during ICRF operations: experiments and modeling

For successfully heating plasma with waves in the ion cyclotron range of frequencies (ICRFs), mitigating impurity production is just as crucial as maximizing power coupling, especially in high-Z environments. ICRF can effectively deposit energy on ions, modify turbulence-driven transport, and enhance fusion reaction efficiency, but only when its power coupling has minimal impact on impurity production. To do so, one must rely on a toroidal array of at least three active elements excited with appropriate phasing and power ratio to reduce the currents induced on the antenna frame below levels critical for physical sputtering. In contrast to classic two-strap antennas, which are optimized for dipole phasing with equal power on both straps, three-strap antennas in ASDEX Upgrade (AUG)—but also four-strap antennas in JET, Alcator C-Mod, SPARC and ITER—offer the possibility to act also on the power ratio between the central and outer straps. With optimal settings, impurity production can be reduced substantially, making the ICRF compatible with the high-Z wall. This paper explores the characteristics of the AUG three-strap antennas in terms of impurity production, as well as the key role of plasma composition in this process. Numerical simulations were performed using SSWICH and Petra-M (finite element codes) to quantify impurity production and compare with experimental results. Energies of ions falling on antenna limiters (measured with probes) are well predicted by both codes. These tools are then used to further describe the source of the impurity, namely the gross erosion of tungsten from an ICRF antenna, for different plasma mixtures. Results are also compared to spectroscopy data. Ultimately, we show that deleterious effects of the ICRF on plasma surface interactions will be weaker in plasmas containing larger fractions of highly ionized heavier low-Z impurity, which is typically relevant for experiments relying on impurity seeding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Predicting Thermal Excursions During in situ Oxidative Regeneration of Packed Bed Catalytic Fast Pyrolysis Catalyst

Ex situ catalytic fast pyrolysis (CFP) uses a secondary reactor to upgrade biomass pyrolysis vapors to stabilized CFP oils with reduced oxygen content. In one configuration, the secondary reactor is operated as a packed-bed swing reactor system which allows coke-deactivated beds to be decarbonized in situ while other beds remain online for vapor upgrading. In situ decarbonization must be done carefully to avoid irreversible deactivation and/or physical degradation of catalyst pellets. Given that packed bed reactors are well known to have poor heat transfer characteristics, this is a critical issue impacting scaleability and commercial viability of the technology. To predict thermal excursions during regeneration, finite element computational models have been built to assist in scaling up oxidative decarbonization of a Pt/TiO2 CFP catalyst (0.5 mm spheres) from a bench scale packed bed with 100 g of catalyst to a pilot scale packed bed with 2 kg of catalyst and internal cooling tubes. Based on transient measurements of outlet temperature and effluent CO2 concentration, and using an assumed coke profile and activation energy, this paper demonstrates that specific combinations of effective thermal conductivity and wall heat transfer coefficient can fit bench scale oxidative regeneration data equally well. For the upscaled 2 kg bed, four bench-scale “best fit” parameter pairs give different predictions for location and magnitude of thermal excursions, with the maximum computed bed temperature gradients ranging from 30 °C cm-1 to as high as 3000 °C cm-1. The larger the fraction of heat removal by conduction through the cooling tubes, the greater the differences between the parameter pairs. The modelling results presented in this paper cast doubt on the industrial viability of the proposed combination of catalyst, bed and regeneration process, and point to the need for alternate reactor designs. However, there is considerable uncertainly in some of the key model parameters. The reliability of model predictions can be increased by adding more temperature measurements at key bed locations, testing additional variations in process conditions, performing careful bed dissections to determine the true coke profile, and perhaps most importantly, directly measuring the effective thermal conductivity of the catalyst pellets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

TRACER-VOC: Speciated Volatile Organic Compounds at the Tracking Aerosol Convection Interactions Experiment Field Campaign Report

Formation of secondary organic aerosol (SOA) is responsible for the growth of particles into cloud condensation nuclei (CCN)-relevant sizes and accounts for a substantial fraction of aerosol mass in the troposphere. Despite its wide abundance and importance, capturing SOA in regional and climate models is still challenging due to complex chemical transformations involved in production of SOA. Oxidation products of biogenic hydrocarbons, such as isoprene, monoterpenes, and sesquiterpenes account for the largest fraction of SOA mass over the continents. As such, biogenic SOA is an important source of continental CCN, and has been recently shown to nucleate ice in the deposition mode. To better account for the processes involved in SOA formation in global and regional climate models, we need better measurements of sources and transformations of volatile organic compounds (VOCs), the precursors to SOA formation in the atmosphere. Understanding aerosol-deep convection interactions is an important science driver of the Tracking Aerosol Convection Interactions Experiment (TRACER), requiring a characterization of the critical controls on aerosol properties and processes, and how those properties and processes vary spatially and temporally in the Houston area. In that environment, the atmospheric photo-oxidative chemistry that drives the evolution and transformation of aerosol properties and processes is strongly impacted by the emission of VOCs, which have both anthropogenic (e.g., industrial, urban, transportation, biomass burning) and biogenic (e.g., marine) sources in the region. Understanding (and disentangling) the complex role of VOCs in such a heterogenous region requires high-time-resolution, speciated VOC measurements that are spatially resolved. To meet these measurement requirements, the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s time-of-flight proton transfer reaction-mass spectrometer (ToF-PTR-MS) was deployed during the TRACER intensive operational period (IOP) from June to September, 2022 at the ancillary site (ANC), which was sited to provide a characterization of the rural atmospheric environment in contrast to the main site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Preheating and premixing effects on NO x emissions in a high-pressure axially staged combustor

NO x emissions remain a primary concern for modern and future gas turbines, particularly as we progress towards a net-zero carbon future. It is critical to identify novel strategies to mitigate environmental pollutants for both modern turbines operated on natural gas and future turbines projected to use carbon-free fuels (such as hydrogen). In the current study, NO x emissions of a reacting methane-air jet in a vitiated crossflow are experimentally investigated in a model axially staged combustor at 5 atm. The combustion products from the main stage combustor flow into the axial stage at temperatures ranging from 1580 to 1650 °C. The axial stage contains a reacting jet in crossflow, which provides an overall temperature rise ranging from 50 to 220 °C to generate exit temperatures similar to gas turbine engines. The focus of this paper is to explore the effects of flame liftoff and ignition timescales on the NO x contribution of this secondary jet. This is done by varying the premixed level and preheat temperature of the axial jet. The effects are explored at different momentum flux ratios (J), jet equivalence ratios (φ jet ), exit temperatures (T exit ), and main stage (headend) temperatures (T HE ). High-speed CH* chemiluminescence imaging is employed for each test case to determine the flame stabilization point. For a given temperature rise and combustor exit temperature, an increase in ignition delay was observed with a decrease in jet temperature. This provided a NO x benefit, particularly for jets that are injected at or close to an ignitable mixture fraction. Similarly, the coaxial injector caused a significant ignition delay compared to the fully premixed injector; this also led to a NO x benefit. Overall, the results show that a thermally driven delay in the axial stage provides a greater NO x reduction than a mixing delay.

03 NATURAL GAS↗
DOE OSTI2022

Mind-the-gap - Part III: Doppler velocity measurements from space

Convective motions and hydrometeors microphysical properties are highly sought-after parameters for evaluating atmospheric numerical models. With most of the Earth's surface covered by water, space-borne Doppler radars are ideal for acquiring such measurements at a global scale. While these systems have proven to be useful tools for retrieving cloud microphysical and dynamical properties from the ground, their adequacy, and specific requirements for spaceborne operation still need to be evaluated. Comprehensive forward simulations enable us to assess the advantages and drawbacks of six different Doppler radar architectures currently planned or under consideration by space agencies for the study of cloud dynamics. Radar performance is examined against state-of-the-art numerical model simulations of well-characterized shallow and deep, continental, and oceanic convective cases. Mean Doppler velocity (MDV) measurements collected at multiple frequencies (13, 35 and 94 GHz) provide complementary information in deep convective cloud systems. The high penetration capability of the 13-GHz radar enables to obtain a complete, albeit horizontally under-sampled, view of deep convective storms. The smaller instantaneous field of view (IFOV) of the 35-GHz radar captures more precise information about the location and size of convective updrafts above 5-8 km height of most systems which was determined is the portion of storms where the mass flux peak is typically located. Finally, the lower mean Doppler velocity uncertainty of displaced phase center antenna (DPCA) radars makes them an ideal system for studying microphysics in shallow convection and frontal systems, as well as ice and mixed-phase clouds. It is demonstrated that a 94-GHz DCPA system can achieve retrieval errors as low as 0.05-0.15 mm for raindrop volume-weighted mean diameter and 25% for rime fraction (for a -10 dBZ echo).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Reducing the L-H transition power threshold in ITER-similar-shape DIII-D hydrogen plasmas

Recent dedicated DIII-D experiments in low-torque, ITER-similar-shape (ISS) hydrogen plasmas (at a plasma current I p ~ 1.5 MA and ITER-similar edge safety factor q 95 ~ 3.6) show that the L-H transition power threshold P LH can be reduced substantially (~30%) with L-mode helium admixtures n He /n e $\leqslant$ 25%. In the ensuing H-mode, helium ion fractions n He /n H remain below 25%. H-mode normalized pressure and confinement quality are only slightly affected by helium seeding, and Z eff $\leqslant$ 2.15 (including helium and carbon content). The plasmas investigated here are electron-heat dominated, with temperatures T e (0)/T i (0) $\geqslant$ 1 and edge heat flux ratio Q e /Q i (ρ = 0.95) ~ 1.2–1.5. Without mitigation, P LH is higher by a factor of 2–3 in comparison to similar ISS deuterium plasmas. ISS hydrogen plasmas with lower plasma current I p ~ 1 MA (increased edge safety factor q 95 ~ 5.1) exhibit a substantially lower power threshold. This plasma current dependence, also observed previously on ASDEX-U and in JET, is not accounted for by the commonly used 2008 ITPA multi-machine threshold scaling, but could potentially allow H-mode access at marginal heating power during the initial plasma current ramp-up. Attempts to reduce P LH with low-field- and high-field-side hydrogen pellet injection, using 1.7 mm diameter pellets, have not demonstrated a robust threshold reduction, in contrast to successful earlier experiments with larger 2.7 mm pellets. Finally, techniques for reducing P LH are very important for ITER, in particular for accessing H-mode in hydrogen plasmas during the Pre-Fusion Power Operation-1 (PFPO-1) campaign with marginal auxiliary heating power (20–30 MW of ECH).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Evaluation of New Additions to OLI Software in Predicting Mercuric and Mercurous Species in Liquid Waste Operations

Speciation of mercury during the pretreatment steps of tank waste processing is critical to successful mercury removal prior to vitrification during Liquid Waste Operations (LWO) at SRS. OLI software has been used to predict mercury speciation and activity throughout LWO. The OLI software operates based on a thermodynamic framework called the Mixed Solvent Electrolyte (MSE) framework. The MSE framework allows prediction in theoretically infinitely dilute to concentrated mixtures (e.g., purely solute solutions). Before modification to the MSE framework databanks, certain critical mercury species were missing in the MSE databank, and some thermodynamic data needed to be updated for the OLI software to accurately predict mercury chemical species in SRS waste tanks. To better reflect streams across LWO, new mercury species were integrated into the MSE database. To evaluate the changes to the OLI MSE framework per the Technical Task Request (TTR) and the Task Technical and Quality Assurance Plan (TTQAP), waste stream compositions from Tanks 38, 43, and Tank 50 decontaminated salt solution (DSS) were used as model inputs. Models were developed and executed using both the old and new databases. Compositional analyses from caustic Tank 50 DSS and caustic Tanks 38 and 43 were used as the input streams. These streams represent the most comprehensive chemical data sets where both mercury and tank constituents were measured together. Results for Tank 50 DSS predict HgO as the predominant species in both databases. Both methyl and dimethyl Hg species are present when the new database is ‘on’ and are not predicted with the new database turned ‘off’. The new database predicts a greater amount of HgO and a greater fraction of it in the solid phase. Pourbaix diagrams (potential vs. pH) generated for each Tank 50 DSS were identical regardless of which database was used. Elemental Hg and HgO were predicted in the water stable region under basic conditions. Tanks 38 and 43 follow similar trends as the Tank 50 DSS models. Unlike Tanks 38 and 50 DSS, the Tank 43 Pourbaix plot shows a region of stability for an aqueous HgOHCO3 - species between approximately pH 7-11. In all streams, when MeHg+ is included in the inputs, the new database predicts aqueous MeHgOH as the dominant species. If elemental or dimethyl mercury is in the waste stream, the new database model predicts they are unchanged and remain in those states and quantities. Additionally, the total mercury values are reported for both the measured input data and the OLI output data for all considered tanks. The summary indicates that the percentage error between the measured and calculated values is less than 1% in all cases The reconciliations and generation of the Pourbaix diagrams for Tank 50 DSS took approximately ten times longer with the new database ‘on’. In addition, over the course of that time, models with the new database ‘on’ were more likely to crash or display an error. Some modest performance improvements were noted when modeling with an i7 processor versus an i5. An example error is found in Appendix A. Furthermore, Appendix B provides V&V for two chemical systems analyzed with the OLI software, results were satisfactory. It is recommended to utilize the new databases (i.e., HCO.ddb and SR-Hg.ddb) in future Savannah River Mission Completion applications of OLI to represent pseudo steady-state. Furthermore, the integration and utilization of the new databases (i.e., HCO.ddb and SR-Hg.ddb) in modeling applications (e.g., Aspen) is also recommended.

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

Minimal Energy Routing of a Leader and a Wingmate with Periodic Connectivity

We consider a route planning problem in which two unmanned vehicles are required to complete a set of tasks present at distinct locations, referred to as targets, with minimum energy consumption. The mission environment is hazardous, and to ensure a safe operation, the UVs are required to communicate with each other at every target they visit. The problem objective is to determine the allocation of the tasks to the UVs and plan tours for the UVs to visit the targets such that the weighted sum of the distances traveled by the UVs and the distances traveled by the communicating signals between them is minimized. We formulate this problem as an Integer program and show that naively solving the problem using commercially available off-the-shelf solvers is insufficient in determining scalable solutions efficiently. To address this computational challenge, we develop an approximation and a heuristic algorithm, and employ them to compute high-quality solutions to a special case of the problem where equal weights are assigned to the distances traveled by the vehicles and the communicating signals. For this special case, we show that the approximation algorithm has a fixed approximation ratio of 3.75. We also develop lower bounds to the optimal cost of the problem to evaluate the performance of these algorithms on large-scale instances. We demonstrate the performance of these algorithms on 500 randomly generated instances with the number of targets ranging from 6 to 100, and show that the algorithms provide high-quality solutions to the problem swiftly; the average computation time of the algorithmic solutions is within a fraction of a second for instances with at most 100 targets. Finally, we show that the approximation ratio has a variable ratio for the weighted case of the problem. Specifically, if ρ denotes the ratio of the weights assigned to the distances representing the communication and travel costs, the algorithm has an a posteriori ratio of $3 + \frac{3ρ}{4}$ when ρ ≥ 1, and $\frac{3}{ρ}$ + $\frac{3}{4}$ when ρ ≤ 1.

42 ENGINEERING↗
DOE OSTI2024

Modeling Flow and Particle Deposition in a Spent Nuclear Fuel Assembly

CFD (Computational Fluid Dynamic) simulation of aerosol-laden natural convective flow and particle deposition in a spent fuel storage canister with 37 assemblies is currently computationally prohibitive. PWR (Pressurized Water Reactor) assemblies have up to 289 pins or tubes with several spacer grids to align the pins. Spacer grids with mixing vanes induce swirling during operation to increase heat transfer. Each spacer grid contains hundreds of small structures such as retaining clips, channel walls, and openings. The largest canisters store 37 PWR assemblies thus, there are numerous pins, tubes, and spacer grids for which the flow region between and around these structures needs to be determined along with the movement and deposition of aerosol particles. Because of the complicated geometry, modeling the intricate flow even for just one assembly is currently impractical. Nonetheless, we are developing techniques for a practical model to assess the natural aerosol particle deposition process in a canister in the event that a release occurs from one or more fuel pins. In the previous work it was demonstrated that CFD can model the flow through a PWR spacer grid with mixing vanes, including particle deposition, in a reasonable amount of time on a personal computer. In this work, the analysis is extended to include the bypass region between an assembly and the canister basket walls. It is shown that the flow velocity in the bypass region is about three times that of the interstitial region between the pins. The lengths before and after the spacer grid are also extended to determine when the flow becomes fully developed. In addition, the approach of computationally “stitching together” segments of an assembly is demonstrated with the plan to ultimately model a full assembly. The fraction of particles that are deposited in a segment with a spacer grid is determined as a function of particle size and flow velocity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Investigation of Thermal-Hydraulic Effects of Dry Storage Canister Helium Backfill Loss Using the Horizontal Dry Cask Simulator

A previous investigation produced data sets that can be used to benchmark the codes and best practices presently used to determine cladding temperatures and induced cooling air flows in modern horizontal dry storage systems. The horizontal dry cask simulator (HDCS) was designed to generate this benchmark data and add to the existing knowledge base. The objective of the previous HDCS investigation was to capture the dominant physics of a commercial dry storage system in a well-characterized test apparatus for a wide range of operational parameters. The close coupling between the thermal response of the canister system and the resulting induced cooling air flow rate was of particular importance. The previous investigation explored these parameters using helium backfill at 100 kPa and 800 kPa pressure as well as air backfill with a series of simulated decay heats. The helium tests simulated a horizontal dry cask storage system at normal storage conditions with either atmospheric or elevated backfill pressure, while the air tests simulated horizontal storage canisters following a complete loss of helium backfill, in which case the helium would be replaced by air. The present HDCS investigation adds to the previous investigation by exploring steady-state conditions at various stages of the loss of helium backfill from a horizontal dry cask storage system. This is achieved by using helium/air blends as a backfill in the HDCS and running a series of tests using various simulated decay heats to explore the effects of relative helium/air molar concentration on the thermal response of a simulated horizontal dry cask storage system. A total of twenty tests were conducted where the HDCS achieved steady state for various assembly powers, representative of decay heat. The power levels tested were 0.50, 1.00, 2.50, and 5.00 kW. All tests were run at 100 kPa vessel pressure. The backfill gases used in these tests are given in this report as a function of mole fraction of helium (He), balanced by air: 1.0, 0.9, 0.5, 0.1, and 0.0 He. Steady-state conditions (where the steady-state start condition is defined as where the change in temperature with respect to time for the majority of HDCS components is less than or equal to 0.3 K/h) were achieved for all test cases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Non-Fouling, Low Cost Electrolytic Coagulation & Disinfection for Treating Flowback and Produced Water for Reuse

Executive Summary This Final Report is composed of two major sections. The first section presents experimental results obtained in the laboratory during Budget Period 1. The second section presents results from the field test conducted during Budget Period 2. Laboratory Results This research investigated a novel electrochemical process for producing a ferric iron coagulant for use in treating flowback and produced water from hydraulic fracturing and oil production operations. The treatment system improves the effectiveness and lowers the cost of coagulation processes using Fe3+ as the coagulant. The electrolytic coagulant generation (ECG) system uses an electrochemical cell to produce acid and base from oilfield brine solutions. The acid is used to dissolve scrap iron to provide a Fe3+ coagulating agent. The base is used to neutralize the treated water. Compared to conventional electrocoagulation (EC), the main advantage of the ECG system is an order of magnitude lower cost for the source of iron. The second advantage over conventional EC is that it can deliver Fe3+ doses greater than 1 mM, since it is not limited by the amount of dissolved oxygen in the water required to oxidize ferrous to ferric iron. The capital costs for conventional EC and the ECG system are similar, but the operational costs for the ECG system are an order of magnitude lower than conventional EC. The combined costs for iron and electrical energy for treating 1 m3 of FPW with 1 mM Fe3+ is estimated to be $0.87 for conventional EC, and $0.087 for the ECG system. The estimated all-in cost for treating FPW with a 2 mM Fe3+ dose is $0.73/m3 ($0.12/bbl). Field Test Results The field test was conducted at the Paul Foster Central Tank Battery (CTB) in Lea County, New Mexico from November 10, 2022 through December 15, 2022. The feed water to the system was produced water from the Tatanka 1H formation. After approximately two weeks of testing, the initial batch of produced water had been treated and no untreated produced water was available. Thus, after this time, the feed water to the system consisted of previously treated water (i.e., recycled water). The recycled water had nearly all colloidal particles removed, and had a much lower alkalinity due to precipitation of carbonate minerals during the first pass through the system. Although the recycled water was not an ideal test solution due to its low particulate concentrations, its lower alkalinity did allow us to identify the main problem with the treatment system. The main problem with the treatment system was caused by the high alkalinity of the initial feed water (5.4 meq/L) that consumed a significant fraction of the electrochemically generated acid. This resulted in pH values exiting the iron contact tank that were too high to dissolve enough iron to effectively treat the produced water. Tests performed with recycled water with lower alkalinity did not have this issue, and dissolved iron concentrations greater than 20 mM could be achieved. One consistent observation was that effluent water from the iron contact tank was always free of particulates, even when fed with circumneutral solutions. This suggests that there is no need to dissolve high concentrations of iron if all the water to be treated is passed through the scrap iron canister. In this case, dissolved O2 and hypochlorous acid can promote sufficient iron corrosion to provide an effective coagulating agent – even in neutral pH water. This solves the problem resulting from highly alkaline produced water. Modifications to the design of the treatment system were made based on the field test results. These modifications will add minimal additional cost, and were tested in bench-scale laboratory experiments. In short-term testing, the modified treatment process was able to remove colloidal FeS particulates to levels below detection. Long-term, steady state testing will be required to determine whether the modified process is suitable for commercial treatment systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Overview of the ASDEX Upgrade results

After a 26-month vent ASDEX Upgrade (AUG) went back in operation with a newly designed upper W-divertor suitable for alternative divertor configurations (featuring in-vessel coils and cryo-pump). Parameter scans and an extensive set of measurements were obtained and their interpretation is ongoing. Prompted by the ITER wall change, dedicated experiments on non-boronized plasma startup were contrasted to that employing asymmetric and more symmetric boronizations. The asymmetric boronization proved to be as beneficial as the more symmetric one, which is in contrast to previous model calculations assuming perfect sticking of boron (measurements suggest sticking ≈ 0.3). In the startup phase also the impurity influxes at the outboard limiters were investigated contrasting the unboronized case featuring cold edges (low- Z radiation) to the boronized case, in which the lifetime of the boron layers could be estimated. Pedestal stability investigations revealed that the quasi-continuous exhaust (QCE) regime is obtained when ballooning modes are active in the vicinity of the separatrix and the global peeling-ballooning stability is high enough. Thus, at high enough shaping and high gas flux both can be achieved and QCE is a consequence. The closely related enhanced D-alpha (EDA) mode is not clearly distinguishable from QCE, e.g. the quasi-coherent mode characteristic for EDA also shows up in QCE. In QCE the impurity transport is behaving benign as could be measured for Ne with a novel analysis method making use of a comprehensive set of CXRS measurements. For high radiative fractions the regime of the X-point radiator (XPR) is accessible at AUG and the understanding of its access conditions and behaviour is further developed. Due to the localized radiative cooling at the X-point the XPR can be well diagnosed and thus controlled. For negative triangularity shapes, further experiments at increased shaping resulted in strongly heated L-mode plasmas avoiding ELMs. Two integrated modelling approaches towards ITER suggested that core W-accumulation will be no issue for ITER and that the fusion yield in ITER may be Q = 12 (i.e. ITPA20-IL scaling is too pessimistic). Further, investigations of the ITER ramp-down in AUG provide insights into maintaining position control. Various aspects of shattered pellet injection were investigated in AUG and one of the results show that with increasing Ne fraction the radiation during the current quench increases and the current decay becomes faster.

ASDEX Upgrade↗
DOE OSTI2023

Core-Pedestal Plasma Configurations in Advanced Tokamaks

Here, several configurations for the core and pedestal plasma are examined for a predefined tokamak design by implementing multiple heating/current drive (H/CD) sources to achieve an optimum configuration of high fusion power in a noninductive operation while maintaining an ideally magnetohydrodynamic (MHD) stable core plasma using the IPS-FASTRAN framework. IPS-FASTRAN is a component-based lightweight coupled simulation framework that is used to simulate magnetically confined plasma by integrating a set of high-fidelity codes to construct the plasma equilibrium (EFIT, TOQ, and CHEASE), calculate the turbulent heat and particle transport fluxes (TGLF), model various H/CD systems (TORIC, TORAY, GENRAY, and NUBEAM), model the pedestal pressure and width (EPED), and estimate the ideal MHD stability (DCON). The TGLF core transport model and EPED pedestal model are used to self-consistently predict plasma profiles consistent with ideal MHD stability and H/CD (and bootstrap) current sources. In order to evaluate the achievable and sustainable plasma beta, varying configurations are produced ranging from the no-wall stability to with-wall stability regimes, simultaneously subject to the self-consistent TGLF, EPED, and H/CD source profile predictions that optimize configuration performance. The pedestal density, plasma current, and total injected power are scanned to explore their impact on the target plasma configuration, fusion power, and confinement quality. A set of fully noninductive scenarios are achieved by employing ion-cyclotron, neutral beam injection, helicon, and lower-hybrid H/CDs to provide a broad profile for the total current drive in the core region for a predefined tokamak design. These noninductive scenarios are characterized by high fusion gain (Q ~ 4) and power (P fus ~ 600 MW), optimum confinement quality (H 98 ~ 1.1), and high bootstrap current fraction (f BS ~ 0.7) for Greenwald fraction below unity. The broad current profile configurations identified are stable to low-n kink modes either because the normalized pressure β N is below the no-wall limit or a wall is present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Methodology to evaluate design modifications intended to eliminate frosting and high discharge temperatures in air-source heat pumps (ASHPs) in cold climates

Air-source heat pumps (ASHPs) operating in cold climates experience problems with frosting and high refrigerant temperatures. These problems increase energy consumption, and their severity depends on the climatic conditions. In the present paper, a methodology for identifying the prevailing problem between frosting and high discharge temperatures is presented. Three performance indices, the frosting index (FI), the discharge index (DI), and the total loss index (TLI), are proposed to quantify the impacts of frosting and high discharge temperatures on the annual performance of ASHPs in different climatic conditions. The FI and DI show which problem (frosting or high discharge temperature) dominates, and the TLI indicates the combined effect of frosting and high discharge temperatures on the performance of an ASHP. A thermodynamic model of an ASHP coupled with the TRNSYS building simulation tool is used to estimate the performance of an ASHP and the proposed loss indices to estimate the impact of both frosting and high discharge temperatures for 45 cities in Canada. The results can be extended to other parts of the world that experience similar climatic conditions The results reveal that in cities in ASHRAE climatic zones 5 and 6 (classified as cold regions) where the ambient air temperatures are predominantly between -15 °C to 6 °C, ASHPs are heavily impacted by frosting. The problem of high discharge temperatures in ASHPs is predominant in cities in climate zones 7 and 8 (classified as very cold and subarctic regions) where the temperatures are frequently below -20 °C in winter. Among the cities considered, St. John, NL has the highest fraction of heating hours experiencing frosting (90 %), where the annual increase in energy consumption due to frosting is 13.5 % of the annual heating energy consumption. The highest annual increase in energy consumption due to high discharge temperatures is in Isachsen, NU (zone 8), where the increase is 30 % of the annual heating energy consumption. Based on the proposed indices, another index called the performance gain index (PGI) is created, which can be used as a first step to assess the energy-saving potential of design modifications applied to ASHPs to solve the problems of frosting and high discharge temperatures. The PGI will aid in developing climate specific ASHPs. One possible design modification is the use of a two-stage ASHP with an economizer. It is observed that the two-stage ASHP with economizer can mitigate high discharge temperatures and improve performance in very cold and subarctic regions (zones 7 and 8). However, it is not as beneficial in zones 5 and 6, where the impact of high discharge temperatures on performance is minimal and frosting dominates. Finally, a case study, using the PGI to evaluate the economic and environmental effectiveness of a two-stage ASHP with economizer is presented for the city of Saskatoon.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Production of Fischer-Tropsch Synfuels at Nuclear Plants

A case study analysis was performed to evaluate nuclear-powered synthetic fuel production in the midwestern United States (U.S.). A Fischer-Tropsch (FT) fuel synthesis plant design was used as the basis for the analysis. The FT plant design was configured to produce a product slate consisting of diesel fuel, jet fuel, and motor gasoline blend stocks from carbon dioxide (CO 2 ) and hydrogen (H 2 ) feedstocks. The CO 2 feedstock for the FT plant was assumed to be sourced from biorefineries in the region around a Midwest light water reactor (LWR) nuclear power plant (NPP). The analysis specifies that power from the LWR is used to produce H 2 via high-temperature steam electrolysis and to operate the FT synfuel production plant. Capital costs were estimated for the FT plant while capital costs for the electrolysis plant were based on previous Idaho National Laboratory (INL) studies. In addition to labor and maintenance costs for the FT and electrolysis plants, operating costs also include the costs for CO 2 feedstock transport. An analysis was performed to determine the cost of transporting CO 2 from the distributed biorefinery sources to the centralized fuel synthesis plant as a function of the synfuel plant capacity and corresponding CO 2 demand. The primary revenue streams are associated with sales of the synthetic fuel products. The synthetic fuel products will likely follow the same market trends as the conventional fuel products. The synfuel price data was thus based on projections made by the U.S. Energy Information Administration (EIA) 2021 Annual Energy Outlook (AEO) for conventional fuel products minus federal and state taxes, as well as marketing and distribution costs. The economic analysis also considered cases that included and excluded revenues from the 2022 Inflation Reduction Act (IRA) clean hydrogen production tax credit (PTC) of $\$ $3.00/kg for the first ten years of operation. The economic analysis calculated the net present value (NPV) for cases involving steady-state synfuel production for comparison with the NPV for a business-as-usual case in which NPP continues to sell only electric power to the grid. A synfuel production “Reference Case” was considered in addition to sensitivity cases in which the plant capacity, electricity price, and synthetic fuel product prices were perturbed. The synfuel production Reference Case considered a scenario in which the electrolysis and synfuel plants utilized a combined electrical load of 1000 megawatt electrical (MWe) from the LWR with the balance of the LWR power output being sold to the electric grid. The economic analysis suggests that the synfuel production Reference Case evaluated in this analysis would lead to considerable economic potential for near-term deployment of a nuclear-based synfuel production plant. Specifically, the economic analysis suggests that the deployment of a 1000 megawatt (MW) nuclear-powered synfuel plant could result in a NPV increase of approximately $\$ $1.7 billion for a case with no clean synfuel price premium relative to conventional petroleum fuels when accounting for the additional revenues from the 2022 IRA clean hydrogen PTCs of $\$ $3/kg. Sensitivity analysis was performed to evaluate the effect of perturbation of selected model input parameters on the NPV for the synfuel production Reference Case. The sensitivity analysis indicates that the plant capacity has the largest impact on the differential NPV, with a smaller synfuel production capacity resulting in a decrease in revenue when a larger fraction of the power from the NPP is sold to the grid and a smaller fraction of the power is used to produce synthetic fuel products. The synfuel product pricing has the next largest impact on the differential NPV, with lower synfuel prices resulting in decreased NPV from decreased synfuel sales revenue while higher synfuel prices result in increased NPV from increased synfuel sales revenue. Electricity pricing has a smaller effect on the NPV than the fuel sales price since, in the Reference Case, most of the energy from the NPP is used for synfuel production and a smaller amount of the system revenues are associated with electrical power sales. However, the electricity price sensitivity does indicate that the Synfuel Integrated Energy System (IES) would have a greater NPV than the business-as-usual case (e.g., grid power sales only) when electricity market prices are low, suggesting that synfuel production could provide a strategy for decreasing the economic risks to NPPs posed by a loss of revenues attributed to falling electricity market prices.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Low-Temperature Plasma Treatment for Enhanced Recovery of Highly Valued Critical REEs from Coal (Final Report)

The principle objective of this project was to develop a novel process flowsheet using low-temperature plasma treatment (LTP) integrated with conventional hydrometallurgical processes to recover rare earth elements (REEs), especially highly valued REEs (i.e. scandium and critical REEs), from coal and coal byproducts. LTP treatment was used to pretreat the solid feed material prior to leaching in an effort to improve the leaching characteristics of the REEs and potentially reduce acid consumption. Initially, a detailed experimental program was performed on a laboratory LTP furnace to identify a set of process parameter values to be used for assessing the enhanced leaching performance achieved on four specific gravity fractions of West Kentucky No. 13 and Fire Clay coarse coal refuse materials. High-temperature oxidation (HTO)) was also evaluated as a means for improving REE leachability on the same fractions of both coal sources. The lixiviants used to assess leachability improvements included: i) untreated water, ii) 0.1 M ammonium sulfate solution and iii) 1.2 M sulfuric acid solution at a temperature of 75°C. The laboratory data was used to develop multiple process flow sheets using HTO only to produce a rare earth oxide mix concentrate or utilizing both LTP and HTO treatment processes to generate a high purity heavy REO mix concentrate and a light REO mix concentrate. A techno-economic analysis was performed on each circuit. Low-temperature plasma oxidation was found to enhance REE leaching characteristics for all specific gravity (SG) fractions of the West Kentucky No. 13 coarse refuse material. The leaching recovery improvement was preferentially achieved for the heavy REEs with higher values obtained when treating the lower density fraction material. For example, the leaching recovery for the heavy REEs in the untreated 1.6 SG Float fraction material was 8% in a 0.1 M ammonium sulfate solution. After LTP treatment for 32 hours, heavy REE recovery increased to 33%. The corresponding recovery value for the light REEs was 10%. These findings indicate a strong affinity of the heavy REEs with the organic material, which is the component being reduced in content during the LTP treatment. Scanning electron micrographs showed the development of a honeycomb structure with increased LTP treatment with higher surface area as confirmed from BET measurements. Higher surface area promotes improved leaching efficiency due to enhance chemical interactions. When using 1.2 M sulfuric acid, the LTP treatment was found to provide a heavy REE recovery increase from 23% to 53% for the 1.6 SG Float fraction which corresponded to an ash content increase from around 20% to 35%. However, LTP treatment provided no improvement in leaching characteristic for the 2.2 SG sink material in the acid solution, which was likely due to the low organic content. The LTP treatment of the Fire Clay SG fractions provided no improvement in REE leaching characteristic, which may be due to the association with REE minerals such as monazite and REE-associated mineral matter such as clays. High temperature oxidation using temperatures in the range of 600-750OC was found to significantly enhance the leaching characteristics of both the light and heavy REEs for all density fractions in both coarse coal refuse materials. Generally, the increase in light REE recovery for the West Kentucky material was from 15% for the untreated material to 90% for the HTO treated material at 600°C. Recovery values dropped significantly when using temperatures below 600°C and above 750°C. The impact of the excessive temperatures was most likely a result of sintering. For the heavy REEs, HTO significantly improved recovery in the lighter SG fractions, which may be due to an affinity with the organic content. Only limited heavy REE recovery improvement was obtained for the denser SG fractions. Another benefit of HTO is the conversion of calcium-based minerals to a very soluble calcium oxide, which is removable using a mild acid solution, and the conversion of Fe minerals to relatively insoluble minerals such as hematite and magnetite. The pregnant leaching solution produced from the leaching of all SG fractions obtained from both coal sources was treated using two stages of solvent extraction procedurally modeling a rougher-cleaner circuit followed by oxalic acid precipitation. Rare earth oxide products having a purity greater than 97% was produced from nearly all of the SG fractions. A novel process flowsheet was developed based on the selective leaching of heavy REEs that was made possible through LTP treatment of West Kentucky No. 13 coarse refuse followed by salt leaching. In the flowsheet, coal refuse material is ground to a top size of 1mm and then treated by LTP. The treated material is mixed into a 0.1 M ammonium sulfate solution and then filtered to produce a heavy REE enriched pregnant leach solution. The heavy REEs are concentrated by two stages of selective precipitation to produce a heavy REE mix oxalate and then roasted to obtain a heavy REO mix concentrate. The solid residues remaining after salt leaching is dried and then subjected to HTO treatment followed by concentration using selective precipitation and roasting to produce a light REO mix concentrate. An alternative circuit was developed using HTO treatment only to treat 100% of the coarse refuse material followed by selective precipitation and roasting to generate a high purity REO mix product. A techno-economic analysis performed for the treatment of West Kentucky No. 13 and Fire Clay coarse refuse based on a hypothetical 500-tph plant revealed that the net present value for all scenarios was significantly negative. Primary reasons include: i) a substantial decrease in the basket value of the final product relative to the feed and ii) high chemical costs relative to the amount of product generated. The basket value reduction was primarily due to little to no scandium recovery due to poor leaching characteristics when using low acid concentrations as necessitated by the low feed grades. Operating costs were primarily influenced by the chemical costs relative to the kilograms of concentrate produced and feed grade. For the Fire Clay coarse refuse material, the total operating cost was estimated to be $354 per kg of REO produced which was higher than the basket value of the REO product ($276/kg). Approximately 37% of the total cost was the sulfuric acid used during the leaching process. A comparative analysis performed to assess the economic viability of using low-temperature plasma furnaces found that market value enhancement for the generation of both heavy REO and light REO concentrates will be needed in addition to significant technology advancements to reduce the cost of plasma furnaces.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Can a Nuclear-Assisted Biofuels System Enable Liquid Biofuels as the Economic Low-carbon Replacement for All Liquid Fossil Fuels and Hydrocarbon Feedstocks and Enable Negative Carbon Emissions?

This report integrates the results of a series of studies and workshops to address the question: Can a nuclear-assisted biofuels system enable liquid biofuels as the economic low-carbon replacement for all liquid fossil fuels and hydrocarbon feedstocks and simultaneously enable negative carbon emissions? “Economic” is defined as economically competitive relative to other low-carbon replacement options for crude oil. “All” refers to the capability to replace the 18 million barrels of oil per day used by the United States. “Nuclear-assisted” refers to the provision of massive quantities of low-carbon heat and hydrogen at the refinery to fully utilize the carbon content of the cellulosic biomass feedstocks for liquid hydrocarbon production. A system is proposed that decarbonizes about half of the U.S. economy while improving long-term soil productivity and sequestering carbon from the atmosphere. In the U.S. almost half the energy consumed by the final customer is in the form of liquid hydrocarbons produced from crude oil. Liquid hydrocarbons are used as (1) an energy source, (2) a method for daily-to-seasonal energy storage, (3) a chemical feedstock, (4) a chemical reducing agent, (5) a method to enhance high-temperature heat transfer in many furnaces and industrial processes and (6) other purposes. As a consequence, the U.S. consumes ~18 million barrels of crude oil per day to produce liquid hydrocarbons. While there are substitutes for liquid hydrocarbons for some applications, our assessment is that the costs and difficulty will dramatically increase if liquid hydrocarbon use goes much below the equivalent of 10 million barrels per day of crude oil. New uses of liquid hydrocarbons to partly replace coal and natural gas could increase demand beyond the equivalent of 10 million barrels of oil per day. One can produce and burn liquid hydrocarbons from biomass without any net addition of carbon dioxide to the atmosphere. Plants grow by removing carbon dioxide from the atmosphere; thus, burning biomass does not result in any net addition of carbon dioxide to the atmosphere. Biomass is typically 40% oxygen. To remove this oxygen to create hydrocarbon liquids, there are two options. The first option is to use biomass as (1) a feedstock, (2) an energy source to operate the process and (3) a supply of carbon to remove the biomass oxygen as carbon dioxide. The second option is to use external heat and hydrogen to remove the oxygen as water and produce liquid hydrocarbons. The use of massive quantities of external heat and hydrogen for hydrocarbon liquid fuels production reduces the biomass feedstock per unit of liquid hydrocarbon product by more than a factor of two reducing land use by more than a factor of two. Many cellulosic feedstocks unsuitable for liquid hydrocarbon production are viable feedstocks with external heat and hydrogen inputs. As a result, there is sufficient cellulosic feedstocks to meet U.S. and global liquid fuels hydrocarbon demand without significant impacts on food and fiber prices. The heat and hydrogen (heat plus electricity) are produced using base-load nuclear power plants—the most economic form of nuclear energy. If low-price natural gas, there is also the option of producing hydrogen from natural gas with sequestration of the carbon dioxide. The biomass is locally processed in depots to produce commodity feedstocks that can be shipped long distances to large biorefineries (250,000 barrels per day, oil equivalent) similar to existing refineries except modifications of front-end processes. The depot system enables local recycle of nutrients back to the soil. The biorefinery can produce variable quantities of liquid hydrocarbon fuels and carbon dioxide for sequestration enabling removal of carbon dioxide from the atmosphere. Preliminary estimates are that the liquid hydrocarbon costs are equivalent to crude oil at between $60 and $70 per barrel. The largest cost is for hydrogen with biomass feedstocks and refinery costs a smaller fraction of the total cost.

09 BIOMASS FUELS↗