Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source Term Analysis”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

CO 2 to Bioplastics: Beneficial Re-use of Carbon Emissions from Coal-fired Power Plants using Microalgae

This project sought to address the technical and economic barriers to carbon dioxide (CO 2 ) capture and utilization using microalgae. Specifically, a dual photobioreactor (PBR)/open raceway pond (ORP) cultivation system was evaluated with respect to capital and operational costs, productivity, and culture health, and compared to an ORP-only cultivation system. In the dual (or “hybrid”) system, a cyclic flow PBR was used to provide inoculum for two open raceway ponds after each harvest, with the hypothesis that this PBR + pond system should result in improved performance compared to traditional raceway ponds. Two other raceway ponds were operated conventionally as a control, the experiments being performed at Duke Energy’s East Bend Station in northern KY using coal-derived flue gas as the CO 2 source. Although algae productivity achieved was mediocre due to the fact that the experiments had to be performed late in the growing season when climatic conditions were not optimal for high growth, the hybrid cultivation system showed a higher level of algae productivity (statistically significant) compared to the traditional ponds. In order to realize the maximum value of the algal biomass produced, fractionation of the biomass was examined. Bioplastic compounding and material characterization was subsequently performed using three feed stocks: whole biomass, lipid-extracted biomass and the proteinaceous residue from full fractionation. In general, the whole and lipid-extracted biomass gave similar results in terms of the properties of the resulting bioplastics, while the proteinaceous residue showed promise for the production of a polybutylene adipate terephthalate (PBAT) blend. Finally, sustainability assessment of a putative biorefinery system was conducted, including techno-economic and life cycle impact assessment. Nine different production scenarios were considered, comprising combinations of the three different biomass growth architectures (ORP, PBR and dual PBR-ORP systems) coupled with three different algae biomass processing pathways (drying only, lipid extraction and fractionation). Results show that the minimum selling price of the bioplastic feed stock (BPFS) is within the realm of economic competition with prices as low as $970 USD tonne -1 . Additionally, life cycle impact assessment results indicate drastic improvements in performance of the produced BPFS, with reductions in greenhouse gas emissions ranging between 67 and 116% compared to a petroleum based plastic feedstock.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS2015

The Compositional Evolution of C/2012 S1 (ISON) from Ground-Based High-Resolution Infrared Spectroscopy as Part of a Worldwide Observing Campaign

Volatile production rates, relative abundances, rotational temperatures, and spatial distributions in the coma were measured in C/2012 S1 (ISON) using long-slit high-dispersion (lambda/delta lambda approximately 2.5 times 10 (sup 4)) infrared spectroscopy as part of a worldwide observing campaign. Spectra were obtained on Universal Time 2013 October 26 and 28 with NIRSPEC (Near Infrared Spectrometer) at the W.M. Keck Observatory, and Universal Time 2013 November 19 and 20 with CSHELL (Cryogenic Echelle Spectrograph) at the NASA IRTF (Infrared Telescope Facility). H2O was detected on all dates, with production rates increasing markedly from (8.7 plus or minus 1.5) times 10 (sup 27) molecules per second on October 26 (Heliocentric Distance = 1.12 Astronomical Units) to (3.7 plus or minus 0.4) times 10 (sup 29) molecules per second on November 20 (Heliocentric Distance = 0.43 Astronomical Units). Short-term variability of H2O production is also seen as observations on November 19 show an increase in H2O production rate of nearly a factor of two over a period of about 6 hours. C2H6, CH3OH and CH4 abundances in ISON (International Scientific Optical Network) are slightly depleted relative to H2O when compared to mean values for comets measured at infrared wavelengths. On the November dates, C2H2, HCN and OCS abundances relative to H2O appear to be within the range of mean values, whereas H2CO and NH3 were significantly enhanced. There is evidence that the abundances with respect to H2O increased for some species but not others between October 28 (Heliocentric Distance = 1.07 Astronomical Units) and November 19 (Heliocentric Distance = 0.46 Astronomical Units). The high mixing ratios of H2CO to CH3OH and C2H2 to C2H6 on November 19, and changes in the mixing ratios of some species with respect to H2O between October 28 to November 19, indicates compositional changes that may be the result of a transition from sampling radiation-processed outer layers in this dynamically new comet to sampling more pristine natal material as the outer processed layer was increasingly eroded and the thermal wave propagated into the nucleus as the comet approached perihelion for the first time. On November 19 and 20, the spatial distribution for dust appears asymmetric and enhanced in the antisolar direction, whereas spatial distributions for volatiles (excepting CN) appear symmetric with their peaks slightly offset in the sunward direction compared to the dust. Spatial distributions for H2O, HCN, C2H6, C2H2, and H2CO on November 19 show no definitive evidence for significant contributions from extended sources; however, broader spatial distributions for NH3 and OCS may be consistent with extended sources for these species. Abundances of HCN and C2H2 on November 19 and 20 are insufficient to account for reported abundances of CN and C2 in ISON near this time. Differences in HCN and CN spatial distributions are also consistent with HCN as only a minor source of CN in ISON on November 19 as the spatial distribution of CN in the coma suggests a dominant distributed source that is correlated with dust and not volatile release. The spatial distributions for NH3 and NH2 are similar, suggesting that NH3 is the primary source of NH2 with no evidence of a significant dust source of NH2; however, the higher production rates derived for NH3 compared to NH2 on November 19 and 20 remain unexplained. This suggests a more complete analysis that treats NH2 as a distributed source and accounts for its emission mechanism is needed for future work.

Russo, N. Dello↗
DOE OSTI2024

Minimizing grid energy consumption in wastewater treatment plants: Towards green energy solutions, water sustainability, and cleaner environment

Wastewater treatment plants (WWTPs) consume significant amount of energy to sustain their operation. From this point, the current study aims to enhance the capacity of these facilities to meet their energy needs by integrating renewable energy sources. The study focused on the investigation of two primary solar energy systems in As Samra WWTP in Jordan. The first system combines parabolic trough collectors (PTCs) with thermal energy storage (TES). This system primarily serves to fulfill the thermal energy demands of the plant by reducing the demands from boiler units, which allows more biogas for electricity generation. The second system is a photovoltaic (PV) system with Lithium-Ion batteries, which directly produces electricity that will be used to cover part of the electrical energy demands of plant. To assess the optimal configuration, two distinct scenarios have been formulated and compared to the current case scenario (SC#1). The first scenario focuses on maximizing the net present value (NPV) and minimizing the levelized cost of electricity (LCOE). The second scenario is centred on minimizing the levelized cost of heat (LCOH). The findings indicate that both scenarios succeeded in reducing the reliance on the grid to a value that reach 1 %. Moreover, they both reduced biogas percentage in energy production from 88 % to approximately 65 % through the integration of the PV system. In terms of thermal demand, SC#2 reduced the reliance on biogas boiler units from 100 % to 25 %, while SC#3 achieved an even more impressive reduction to just 8 %. The best LCOE value was attained in SC#2, at 0.0895 USD/kWh, with an NPV of 10.54 million USD. Conversely, SC# 3 yielded an LCOH value of 0.0432 USD/kWh th compared to 0.0534 USD/kWh th USD for SC#2. In conclusion, despite their relatively high capital and operating costs, SC#2 and SC#3 managed to substantially decrease the annual electricity expenditure from approximately 2 million USD to 86,000 USD and 0 USD, respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Collaborative Research Proposal: Time Resolved Optical Emissions Spectroscopy and Laser Induced Fluorescence Spectroscopy of Nanosecond Pulsed Discharges in a Gas-Liquid Water Film Reactor

Electrical discharge plasma formed in contact with liquid water is of interest for a wide range of applications in chemical, biomedical, agricultural, electrical, and materials science and engineering. Such plasma reactors are of very timely importance since they also have significant disinfection capability by inactivating bacteria, viruses and other pathogens. Many types of plasma sources including those driven by AC, DC, RF, microwave, and pulsed electrical power supplies coupled to a wide range of different electrode configurations and reactor designs have been developed and explored which contact the plasma with liquid water. Recent roadmaps have recommended that further work is needed to develop our understanding of the fundamental chemical and physical process which occur at the interface of non-thermal plasma with liquid water solutions in order to advance this large diversity of applications which ultimately depend upon efficient production of key reactive chemical species. There is a wide range of interacting factors that affect the chemical reactions that occur in the plasma, in the liquid phase, and at the interface. These factors ultimately govern the key reactive chemical species formed and used in the various applications and they include a) the reactor design and input parameters, b) the discharge and transport processes, c) the plasma properties, and d) the resulting chemical reactions. For pulsed discharges, the power supply design and output parameters control the applied voltage, frequency, rise time, and width (duration) of the applied pulses. The reactor design involves specification of the gas-liquid contacting methods, electrode gap distance, reactor volume and shape, and gas and liquid flow rates achievable. The gas and liquid compositions as well as the liquid properties such as pH and conductivity are also of key importance in determination of the resulting chemical reactions. In addition, the important plasma properties include plasma gas temperature, electron density, electron energy (distribution), and size of the plasma channels which are all affected by the discharge and transport properties. Many studies have focused on specific aspects of these various processes. Of particular important and relevance to the proposed work is the utilization of nanosecond pulses to generate plasma in gas-liquid systems, liquid bubbles, underwater, and in gases. Recent advances in nanosecond pulses provide significant advantages in utilization with liquid water. For example, fast rise time and short pulses are less sensitive to water conductivity, such short pulses may provide advantages in fast temporal quenching of the plasma, and fundamental analysis of pulse properties, including pulse shape and width, may be facilitated by investigation of single filamentary fast pulses. Many studies have also dealt with the role of the gas composition on formation of reactive oxygen species (ROS) (i.e., hydroxyl radicals – ·OH, hydrogen peroxide – H 2 O 2 , various atomic oxygen species, ozone-O 3 , hydroperoxyl radicals HO 2 ·) and reactive nitrogen species (RNS) (i.e., nitrogen oxides – NO, NO 2 , N 2 O - collectively termed NO X , nitrite-NO 2 - , nitrate-NO 3 - , peroxynitrite-ONOO - ). The hydroxyl radical is the critical species in many chemical oxidation reactions for chemical degradation of toxic compounds in water and gases and for synthesis of some compounds. The mixture of various nitrogen oxide species is important for many biochemical and biological processes involved in biomedical and agricultural applications including disinfection and fertilizer production. The present proposal focuses on determination of the effects of time resolved electron density and hydroxyl radicals on plasma chemical reactions through collaboration with the Princeton Collaborative Low Temperature Plasma Research Facility (PCRF) at the Princeton Plasma Physics Laboratory (PPPL). In order to further investigate the role of the plasma generated electrons and hydroxyl radicals on the overall formation of the key species including hydrogen peroxide, hydroxyl radicals, and nitrogen oxides, the proposed work, thus seeks to determine high resolution time resolved electron density by optical emissions spectroscopy and time resolved hydroxyl radicals using laser induced fluorescence in the nanosecond discharge reactor. The combination of data on electron density and hydroxyl radical concentration will be utilized in the present work to more fully characterize the chemical reaction processes in this system and to advance the design, development, and operation of such chemical reactors for a wide range of applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

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

Nuclear Thermal Energy Storage Configurations for Industrial Combined Heat and Power Supply: Conceptual Study and Engineering Designs

The industries examined in this report primarily rely on moderate-temperature heat provided by gas- or coal-fired boilers and combined heat and power (CHP) plants, delivered through standard process steam systems. High-temperature energy demands are often industry-specific and typically exceed the capabilities of high-temperature gas-cooled reactors (HTGRs). While it is technically feasible to replace process steam from fossil-based heat sources with nuclear energy, certain industries, such as methanol production and pulp and paper, face technoeconomic challenges in integrating nuclear energy without major changes or a technological shift. This is mainly due to the limited external energy demand remaining after the use of internal byproducts, waste heat recovery, and simple efficiency improvements. Achieving full decarbonization of these processes with nuclear energy would require significant technological advancements, involving experimental technology and substantial investments, making widespread adoption in existing industrial plants unlikely in the near term. This study reviews TES options in the context of enabling a flexible CHP supply while maintaining a steady nuclear heat input. Heat storage systems that interface between the reactor primary fluid and the CHP system offer superior performance and flexibility. Specifically, steam extraction downstream of the reheater with a two-tank molten-salt TES appears as the best solution regarding thermodynamic system benefits and system drawbacks. Using selected system configurations, a conceptual design of an industrial energy park was developed for industries with varying energy demands, such as steel production plants utilizing electric arc furnaces (EAFs) and chemical plants, as well as for those with constant energy demands, like petroleum refineries. This design highlights the capabilities of TES and explores its potential business cases. The study also conceptually develops the potential for integrating additional energy sources with nuclear systems through the implementation of TES. The potential of the HTGR-TES-CHP system was also evaluated considering key uncertainties such as industrial demand profiles, external grid access availability, and eligible tax credit levels, using the Holistic Energy Resource Optimization Network. Sensitivity of net present value to these uncertainties was analyzed to determine the optimal number of nuclear reactors (and CHP systems) and the suitable TES capacity. The results were interpreted from a decision-maker’s perspective, focusing on three key areas: deployment strategy (oversized units vs. undersized units with TES support), industrial process characteristics (thermal-intensive single profiles vs. electricity-intensive combined profiles), and operational goals (maximizing profits vs. minimizing natural gas (NG) consumption or external grid dependence). The optimization results indicate that the HTGR-TES-CHP system significantly reduces reliance on NG boilers for individual industrial processes by 9-60% (in NG capacity factor), with an average reduction of 38%, compared to standalone NG boiler operation case (Business As Usual [BAU]). For combined industrial processes, the reduction ranges from 37-77%, with an average of 60%. Additionally, the system greatly reduces dependence on external grids. In meeting industrial electrical demands, a 33-100% self-sufficient internal electricity supply is achieved for single industrial process, with an average of 74%, compared to the BAU scenario, where 100% of electricity is imported. For combined processes, 35-100% of internal electricity demands are met by the reactor, with an average of 73%. At last, the relative NG price levels at which the proposed HTGR-TES-CHP system can cost-effectively enter the market currently dominated by existing NG boilers were estimated. For a moderate HTGR CAPEX level ($\$$2500/kWth, $\$$6329/kWe), the analysis suggests that NG prices must be 2.5 to 7 times higher than HTGR variable operating and maintenance costs for single industrial process, and 5.5 to 9.5 times higher for a combined process scenario. Tax credit modeling shows that the Investment Tax Credit significantly reduces the price threshold needed to break even, making the system competitive with NG boilers in certain cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Exploring filamentous fungi depolymerization of corn stover in the context bioenergy queuing operations

Recalcitrance of lignocellulosic feedstocks to depolymerization is a significant barrier for bioenergy production approaches that require conversion of monomeric carbohydrates to renewable energy sources. This study assesses how low-cost modifications in the feedstock supply chain can be transformed into targeted pretreatments in the context of the entire bioenergy supply chain. The aim of this research is to overcome the physiochemical barriers in corn stover that necessitate increased severity in conversion in terms of chemical loading, temperature, and residence time. Corn stover samples were inoculated with a selective (Ceriporiopsis subvermispora) and non-selective (Phaenarochaete chrysosporium) lignin degrading filamentous fungal strains, then stored aerobically to determine the working envelope for fungal pretreatment to achieve lignin degradation. Dry matter loss and gross chemical makeup of corn stover varied by the length of treatment (2 and 4 weeks) and by the moisture content of the treated corn stover samples (40 and 60%, wet basis). Dry matter loss in P. chrysosporium inoculated biomass was elevated compared to the C. subvermispora inoculated biomass; however, treatment also induced additional chemical composition changes suggestive of depolymerization. Scanning electron microscope images reveal hyphae attached within cell lumen and suggest structural changes within P. chrysosporium treated corn stover after 60% moisture storage. These results highlight that fungal treatment approaches must balance loss of convertible material with the potential for reduction in recalcitrance. Techno-economic assessment (TEA) of fungal pretreatment in a short-term queuing system indicated the viability of this approach compared to conventional queuing operations. The total queuing system cost was estimated at $\$$1.65/tonne of biomass stored. After applying the credit of $\$$1.48/tonne from energy savings in the conversion phase using fungal pretreated biomass, the total system cost was $0.80 lower than traditional biomass queueing approach. While the TEA results suggested that treating biomass with C. subvermispora is the most economically viable storage method in the designed fungal-assisted queuing system, future research should focus on additional fungal depolymerization such as those observed in the P. chrysosporium inoculated biomass.

09 BIOMASS FUELS↗
DOE OSTI2020

Accuracy of HVAC Load Predictions: Validation of EnergyPlus and DOE-2 using FLEXLAB Measurements

The aim of the project reported here was to better understand the level of accuracy of three building energy simulation (BES) engines (‘engines’) — EnergyPlus™, DOE-2.1e, and DOE-2.2 — by identifying and investigating significant deviations between the performance predicted by these engines and actual performance as measured in the FLEXLAB® test facility at Lawrence Berkeley National Laboratory (LBNL). The specific test conditions included some of those prescribed in ANSI/ASHRAE Standard 140 - Standard Method of Test for the Evaluation of Building Energy Analysis Computer Programs. Detailed measurements of FLEXLAB performance, including indoor temperatures and heat fluxes and air-flow and water flow rates and temperatures in the Heating, Ventilating and Air Conditioning (HVAC) system, together with hourly weather data, were recorded and used in analyzing the simulation results from EnergyPlus v8.8, DOE-2.2 v3.65 and DOE-2.1e v127. These engines are commonly used in the United States for building energy code compliance, federal, state, and utility incentives programs, as well as energy efficient design of new buildings and energy retrofit of existing buildings. Seven conventional overhead mixing ventilation scenarios were tested and each engine was found to have a similar level of agreement with the measurements of space-level heating and sensible cooling loads. These results provide useful information regarding the accuracy of these engines in predicting the cooling and heating load elements of whole building energy performance. This information is intended for practitioners who are concerned about transitioning between simulation tools with different engines and for managers of utility programs leveraging these tools for evaluating and/or projecting measure savings to be incentivized under their programs. The results of the comparisons of simulated and measured performance indicate that the predictions from all three engines are not significantly different. The 24-hour average value of the absolute mean bias indicates the likely magnitude of the error in any particular case. The average mean bias is reduced by cancelation of overprediction in one case by underprediction in another. The daytime absolute mean biases, which may be more important for both energy performance and occupant comfort, are ~6%, presumably because of the greater complexity involved in simulating in the presence of solar radiation. EnergyPlus typically overpredicts the cooling load and/or underpredicts the heating load by ~1.5% and the DOE-2 engines typically underpredict the cooling load by approximately the same amount. The Root Mean Square Error is relatively more sensitive to shorter term variations in the difference between predicted and measured loads; the three engines have similar values, ~10%, suggesting that the uncertainties in their predictions of peak loads may also be similar in magnitude. The implication of these results is that users, both designers and program analysts, can use EnergyPlus, DOE-2.1e, or DOE-2.2 to model conventional commercial buildings equipped with overhead mixing ventilation with a similar level of confidence. Further work is required to better understand the variability in the level of agreement between the engine predictions and FLEXLAB measurements, where a particular engine will agree well with FLEXLAB in some cases and not so well in others and another engine will agree or disagree in different cases. As the sources of this variability are identified and eliminated or reduced significantly, it is recommended that the experimental capabilities and methods developed in the study reported here should be applied to validating heating and cooling load calculations for spaces with different types of furniture and miscellaneous loads. These methods should then be applied to low energy space conditioning systems in EnergyPlus including, in particular, radiant slab and radiant ceiling panel cooling and heating systems and ‘mixed mode’ systems that combine mechanical cooling and natural ventilation systems, focusing on controls, including control of thermal mass. The work reported here addresses the conventional method of heating and cooling occupied spaces; other methods, such as the use of radiant heating and cooling systems have the potential to provide equivalent occupant comfort, or better, with lower energy consumption. These systems are addressed more explicitly in EnergyPlus but there is a need for empirical validation to give users the same level of confidence in modeling these systems that they have, or should have, in modeling conventional systems, based on the results presented here.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Open Architecture for Cost Savings in Advanced Nuclear Reactors

Recently, nuclear power plant build projects in the West have run over budget due to high capital costs and schedule overruns. Compared to other sources of energy, nuclear power plants have higher capital costs. Reactors are often different at every site, resulting in a lack of standardization. Nuclear is expected to compete with other low carbon sources of energy which have lower capital costs making it essential for nuclear to develop ways of reducing costs. Strategies such as standardization, learning rates, modularization, and schedule reduction in advanced reactors can reduce nuclear costs by about 40%. Standardization as a way of cutting capital costs has been explored even in large nuclear power plants. Standardization of certain plant components can result in lower component and installation costs and higher learning from experience. Standardization can be achieved by adopting a criterion of key performance indicators and general design principles for a specific system or component such as the balance of plant. Modularization allows the construction of certain components of SMRs in a factory, which saves time, increases productivity, and encourages higher learning rates. Production learning decreases the time and the cost related to an activity. The potential for modularized components of advanced reactors to be manufactured in factories makes it conducive to achieving higher learning rates. Developing large-capacity nuclear programs through sequential builds cultivates a higher learning rate, which in effect may reduce schedule overruns. Open architecture has been identified as a way to drive standardization among advanced reactor designs and result in cost savings. Open architecture (OA) is defined as a design enabling a diverse supply chain by defining and publishing requirements of systems or equipment in functional and/or interface terms, utilizing technical standards in widespread use. Currently, the nuclear industry’s approach is to use closed architecture, making most designs proprietary. However, collaboration between various advanced reactor vendors and suppliers utilizing the concept of open architecture can result in modular and standardized architecture of subsystems or subcomponents of a nuclear power plant. Completely standardizing nuclear power plants may be impossible, however, certain common subsystems amongst the various reactor designs could be standardized and/or access a wider supply chain and leverage existing learning from other sectors. Open architecture will save time and allocate resources to the parts of the plants that have the most unique features. A key advantage of open architecture is its ability to improve production learning across advanced reactors (AR) types in the industry, by providing and utilizing the same kind of component. Sodium fast reactor (SFR), High Temperature Gas Reactor (HTGR) and Molten Salt Reactor (MSR) are the advanced reactors considered for this project. This paper aims to determine the cost savings in advanced reactor programs due to open architecture learning rate. This work is an extension of work done on light water reactor small modular reactors; the cost methodology was utilized to investigate the impact of open architecture on advanced reactors with a particular focus on sodium fast reactors. The cost data on sodium fast reactors used in the model presented the most adequate information required for the analysis.

Advanced Nuclear Reactors↗
DOE OSTI2025

TRAILS Output Files

Overview This data repository contains ZIP files that store compressed versions of the output of running the WaterPaths utility planning and management tool in the DU Re-Evaluation mode (to download the tool, please see this GitHub repository). The tool was used to simulate the six-utility North Carolina Research Triangle problem. Details on the contents of each ZIP file can be seen below. Data details Temporal range: Weekly data for 2,344 weeks from 2015 to 2060 (45 years). Spatial range: Six water utilities in the North Carolina Research Triangle region (0: Chapel Hil/OWASA, 1: Durham, 2: Cary, 3: Raleigh, 4: Pittsboro, and 5: Chatham) File types: CSV and OUT Different solutions available The solution numbers correspond to the different pathway strategies (henceforth referred to as "solutions") discussed in paper's main and supporting text (abstract and link to the paper here). They are as follows: Sol92: The Durham-focused pathway strategy Sol132: The Raleigh-focused pathway strategy Sol140: The regionally-robust pathway strategy Objectives files These files can be accessed by unzipping solXX_objectives_pathways.zip that contains 1,000 Objectives_RDMXX_solsXX_to_XX.csv files. Each CSV file will consist of a row representing all the objective values for that specific solution, while every six columns represents the reliability, restriction frequency, infrastructure net present value ($ mil), peak financial cost, worst-case cost, and unit cost ($ per MG; in that order) for each of the six utilities. There will be 1,000 such files, denoting the performance of the six utilities across the 1,000 deeply uncertain states of the world (DU SOWs). Pathway files These files can be accessed by unzipping solXX_objectives_pathways.zip that contains 1,000 Pathways_sXX_RDMXX.out file. Each OUT corresponds to the set of infrastructure being triggered in a specific DU SOW, and each file will have the name file will consist of four tab-delimited columns that are described as follows: Realization: The realization in which an infrastructure options being triggered utility: The utility currently triggering infrastructure week: The week in which a specific infrastructure option is being triggered infra.: The infrastructure option being triggered If the OUT file contains only the header line, no infrastructure was triggered for that specific DU SOW. Policies files These files can be obtained by unzipping Policies.zip. Each of the 1,000 CSV files within the unzipped folder will contain weekly water use restriction policies for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: 0rest_m: restriction multiplier for utility 0 (values between 0 and 1) 1rest_m: restriction multiplier for utility 1 (values between 0 and 1) 2rest_m: restriction multiplier for utility 2 (values between 0 and 1) 3rest_m: restriction multiplier for utility 3 (values between 0 and 1) 4rest_m: restriction multiplier for utility 4 (values between 0 and 1) 5rest_m: restriction multiplier for utility 5 (values between 0 and 1) 0transf: transfer volume for utility 0 (in MGD) 1transf: transfer volume for utility 1 (in MGD) 2transf: transfer volume for utility 2 (in MGD) 3transf: transfer volume for utility 3 (in MGD) 4transf: transfer volume for utility 4 (in MGD) 5transf: transfer volume for utility 5 (in MGD) Water Sources files These files can be obtained by unzipping WaterSources_subset.zip. Each of the 100 CSV files within the unzipped folder will contain weekly state variables at each water source for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: Xvolume: available water volume from source X (in MGD) Xs_area: surface area of source X (in ACF) Xdemand: demand drawn from a water source from source X (in MGD) Xup_spill: upstream spillage from source X (in MGD) Xww_inflow: wastewater inflow from source X (in MGD) Xcatch_inflow: upstream catchment inflow to source X (in MGD) Xevap: evaporation multiplier for source X (values between 0 and 1) Xds_spill: downstream spillage from source X (in MGD) X_Y_alloc_cap: the allocated capacity from source X to utility Y (values between 0 and 1) X_Y_alloc_dem: the allocated demand from source X to utility Y (values between 0 and 1) Xtrmt_alloc_Y: the allocated treatment capacity from source X to utility Y (values between 0 and 1) Utilities files These files can be obtained by unzipping Utilities_subset.zip. Each of the 100 CSV files within the unzipped folder will contain weekly state variables at each utility for all 1,000 hydroclimatic realizations within a specific DU SOW. The column structure is as follows: Xst_vol: total available storage volume of utility X (in MG) Xcapacity: total storage capacity of utility X (in MG) Xnet_inf: : net inflow for all storage infrastructure for utility X (in MGD) Xst_rof: short term ROF for utility X (values between 0 and 1) Xst_stor_rof: short-term storage ROF for utility X (values between 0 and 1) Xst_trmt_rof: short-term treatment ROF for utility X (values between 0 and 1) Xlt_rof: long-term ROF for utility X (values between 0 and 1) Xlt_stor_rof: long-term storage ROF for utility X (values between 0 and 1) Xlt_trmt_rof: long-term treatment ROF for utility X (values between 0 and 1) Xrest_demand: restricted demand for utility X (in MGD) Xunrest_demand: unrestricted demand for utility X (in MGD) Xunfulf_demand: unfulfilled demand for utility X (in MGD) Xwastewater: wastewater return for utility X (in MGD) Xtreat_capacity: total treatment capacity for utility X (in MG) Xcont_fund: reserve (contingency) fund balance for utility X Xins_pout: insurance payout for utility X (% annual volumetric revenue) Xins_price: insurance price for utility X (% annual volumetric revenue) Xinfra_npv: infrastructure net present value for utility ($mil) Xst_vol: total available storage volume of utility X (in MG) Xdebt_serv: debt service for utility X (usually once per year if the infrastructure is triggered; % annual volumetric revenue) Xstor_vol: total stored volume (in MGD) Xobs_ann_dem: observed annual demand for utility X (in MGD) Xproj_dem: projected annual demand for utility X (in MGD) Xpv_debt_serv: present value of debt service payments for utility X (% annual volumetric revenue) Xgross_rev: gross revenue for utility X ($mil) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

Artificial Intelligence↗
DOE OSTI2025

DIF3D-VARIANT 12.0: Updates and New Features

The DIF3D code has been a workhorse of fast reactor analysis work at Argonne National Laboratory for over 40 years. In 1995, a transport option called VARIANT was added to DIF3D to improve the flux solutions for fast reactor problems which we term DIF3D-VARIANT today. DIF3D-VARIANT performs nodal neutron transport calculations using P N or SP N theory in Cartesian and hexagonal two- and three-dimensional geometries. The limited computing capabilities of the time restricted DIF3D-VARIANT to use at most a 6 th order spatial approximation combined with a P3 flux approximation and P1 scattering kernel for a 33 group structure on most studied reactor problems. Computer capabilities have increased steadily since 1995 and today much larger space-angle-energy approximations are possible. This manuscript serves as an update to the theory section of the original DIF3D-VARIANT manual and details more than twenty years of changes made to DIF3D to make version 12 which was released on November 1 st , 2024. The primary focus of the initial work was to extend the space-angle approximations available in DIF3D-VARIANT such that the error due to transport approximations could be better understood. This work was started and completed in 2002 and marked the official version 10. Unfortunately, those higher order approximations could not be used at that time due to the memory constraints of the BPOINTER part of DIF3D (limited to 2 GB). In version 11, completed in 2012, BPOINTER was circumvented in DIF3D-VARIANT for the largest arrays by introducing a Fortran 90 module called LMA (Large Memory Array). This seamlessly replaces all of the functionality of the BPOINTER concept, but it allows 64 bit addressing for every array such that they can be larger than 2 GB. It is now common for DIF3D-VARIANT jobs to consume 50 GB of memory on modern workstations when using high order space-angle approximations and a large number of groups. Many improvements were made to version 11 from 2012 to 2022 when work to create version 12 started. For version 12, several parts of DIF3D were updated to improve performance and thread parallelism was introduced to further reduce the runtime. Numerous minor bugs were discovered in DIF3D-VARIANT as part of the process of creating the perturbation and sensitivity code PERSENT. All of these algorithmic problems were identified in the transition from version 10 to version 11 which prevented DIF3D-VARIANT from running efficiently and reliably. Firstly, the coarse mesh rebalance scheme would routinely diverge and a study detailed in this report demonstrates how it was also typically not effective. This is not a failure of the coarse mesh rebalance methodology, but a failure of its implementation in DIF3D-VARIANT for hexagonal geometries. The fission source extrapolation algorithm was also found to be unreliable on larger group structure problems, leading to divergence in some cases and a negligible improvement in performance overall. Finally, the “Omega” acceleration applied to the partial current solver routine of DIF3D-VARIANT was found to cause DIF3D-VARIANT to converge to the wrong answer. To resolve these issues, both the coarse mesh rebalance and fission source extrapolation were permanently disabled in version 11. The Tchebychev acceleration was put in as a temporary reliable alternative but it is generally inferior to coarse mesh rebalance or coarse mesh finite difference. For the Omega acceleration, the factor was restricted to guarantee that it would not cause follow-on errors in PERSENT. Due to limited funding to support maintenance and development of DIF3D in the last 10 years, no effort was spent since to resolve the outer iteration acceleration. Except for the threading work, all of the changes discussed in this manuscript refer to changes made between version 10 and version 11. Performance comparisons are done to demonstrate the improvements from version 9 to version 12. As will be demonstrated, the updated versi

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A Comprehensive Economic Coal Transition in South Asia

Many countries are considering accelerating their coal transition. A coal transition refers to an energy sector’s shift from a reliance on coal toward an energy mix largely based on cleaner fuels and renewable energy sources. Such a transition is not just related to greenhouse gas emissions, but also encompasses a range of benefits, recognizing that global energy costs and options are changing. Since 2015, proposed new coal power capacity has dropped by three-quarters globally, leaving only a few countries that develop coal-fired power plants at scale (Littlecott et al., 2021). Historic steps were taken at the United Nations Climate Change 26th Conference of Parties (COP26) in Glasgow, as countries pledged to stop new coal builds, end international coal financing, phase down and phase out unabated coal use, and transition to clean energy. In South Asia, there have been several indicators suggesting that countries may be open to moving toward a coal transition. For example, the number of coal power plants under development across South Asia has decreased by 87% since 2015 (Littlecott et al., 2021). However, the challenges of assuring a just transition are substantial. Because coal plays a critical role in the energy and economic systems in South Asia, especially India, moving away from coal means realizing a broader country-wide economic and social transition. A comprehensive, integrated transition strategy for each state is thus needed urgently. This report briefly reviews the current trends and policies on coal in South Asian countries, develops a framework for a comprehensive economic coal transition, and assesses the opportunities and challenges of the transition in key countries. Several important findings emerge from the analysis. First, a coal transition can support overall economic growth and stability. Financial advantages to a well-planned coal transition include mitigating the risk of stranded assets and taking advantage of low-cost renewables. As a global coal transition proceeds, funds are being diverted from new unabated coal power plants, and utilization rates are declining. The likelihood that coal assets will become stranded is increasing, and the potential for future losses therefore increases as well. Second, coal imports in South Asia are rising. Of the coal consumed in Bangladesh, India, Nepal, and Sri Lanka, 32% is imported; this number increases to 94% when excluding India (International Energy Agency [IEA], 2021d). This illustrates a serious energy security risk. One example is the recent increase in coal prices in South Asia, to be discussed in Section 2.2.1. A diverse energy portfolio that incorporates local renewable energy can provide resilience in the face of changing commodity prices and availability. Third, the social benefits of a coal transition include positive health impacts and broader economic improvements in job creation, although assuring a just transition may be a challenge. Phasing out or phasing down coal can significantly reduce air pollutant emissions and therefore minimize associated premature mortality and improve life expectancy. Additional societal benefits of a coal transition include the high economy-wide potential for job creation, although it creates challenges in terms of reintegration and resettlement for coal miners and their communities.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS2016

Earth Science Data Analytics: Preparing for Extracting Knowledge from Information

Data analytics is the process of examining large amounts of data of a variety of types to uncover hidden patterns, unknown correlations and other useful information. Data analytics is a broad term that includes data analysis, as well as an understanding of the cognitive processes an analyst uses to understand problems and explore data in meaningful ways. Analytics also include data extraction, transformation, and reduction, utilizing specific tools, techniques, and methods. Turning to data science, definitions of data science sound very similar to those of data analytics (which leads to a lot of the confusion between the two). But the skills needed for both, co-analyzing large amounts of heterogeneous data, understanding and utilizing relevant tools and techniques, and subject matter expertise, although similar, serve different purposes. Data Analytics takes on a practitioners approach to applying expertise and skills to solve issues and gain subject knowledge. Data Science, is more theoretical (research in itself) in nature, providing strategic actionable insights and new innovative methodologies. Earth Science Data Analytics (ESDA) is the process of examining, preparing, reducing, and analyzing large amounts of spatial (multi-dimensional), temporal, or spectral data using a variety of data types to uncover patterns, correlations and other information, to better understand our Earth. The large variety of datasets (temporal spatial differences, data types, formats, etc.) invite the need for data analytics skills that understand the science domain, and data preparation, reduction, and analysis techniques, from a practitioners point of view. The application of these skills to ESDA is the focus of this presentation. The Earth Science Information Partners (ESIP) Federation Earth Science Data Analytics (ESDA) Cluster was created in recognition of the practical need to facilitate the co-analysis of large amounts of data and information for Earth science. Thus, from a to advance science point of view: On the continuum of ever evolving data management systems, we need to understand and develop ways that allow for the variety of data relationships to be examined, and information to be manipulated, such that knowledge can be enhanced, to facilitate science. Recognizing the importance and potential impacts of the unlimited ways to co-analyze heterogeneous datasets, now and especially in the future, one of the objectives of the ESDA cluster is to facilitate the preparation of individuals to understand and apply needed skills to Earth science data analytics. Pinpointing and communicating the needed skills and expertise is new, and not easy. Information technology is just beginning to provide the tools for advancing the analysis of heterogeneous datasets in a big way, thus, providing opportunity to discover unobvious scientific relationships, previously invisible to the science eye. And it is not easy It takes individuals, or teams of individuals, with just the right combination of skills to understand the data and develop the methods to glean knowledge out of data and information. In addition, whereas definitions of data science and big data are (more or less) available (summarized in Reference 5), Earth science data analytics is virtually ignored in the literature, (barring a few excellent sources).

data analytics↗
DOE OSTI2024

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗
NASA NTRS2004

SAMZ: Satellite-Derived Management Zones

The term "satellite-derived management zones" (SAMZ) denotes agricultural management zones that are subdivisions of large fields and that are derived from images of the fields acquired by instruments aboard Earth orbiting satellites during approximately the past 15 years. "SAMZ" also denotes the methodology and the software that implements the methodology for creating such zones. The SAMZ approach is one of several products of continuing efforts to realize a concept of precision agriculture, which involves optimal variations in seeding, in application of chemicals, and in irrigation, plus decisions to farm or not to farm certain portions of fields, all in an effort to maximize profitability in view of spatial and temporal variations in the growth and health of crops and in the chemical and physical conditions of soils. As used here, "management zone" signifies, more precisely, a subdivision of a field within which the crop production behavior is regarded as homogeneous. From the perspective of precision agriculture, management zones are the smallest subdivisions between which the seeding, application of chemicals, and other management parameters are to be varied. In the SAMZ approach, the main sources of data are the archives of satellite imagery that have been collected over the years for diverse purposes. One of the main advantages afforded by the SAMZ approach is that the data in these archives can be reused for purposes of precision agriculture at low cost. De facto, these archives contain information on all sources of variability within a field, including weather, crop types, crop management, soil types, and water drainage patterns. The SAMZ methodology involves the establishment of a Web-based interface based on an algorithm that generates management zones automatically and quickly from archival satellite image data in response to requests from farmers. A farmer can make a request by either uploading data describing a field boundary to the Web site or else drawing the boundary on a reference image. Hence, a farmer can start to engage in precision farming shortly after gaining access to the Web site, without need for incurring the high costs of conventional precision-agriculture data-collection practices that include collecting soil samples, mapping electrical conductivity of soil, and compiling multi-year crop-yield data. Given the boundary of a field, a SAMZ server computes the zones within the field in a three-stage process. In the first stage, a vector-valued image of the field is constructed by assembling, from the archives, the equivalent of a stack of the available images of the field (see figure). In the second stage, the vector-valued image is analyzed by use of a wavelet transform that detects spatial variations considered significant for precision farming while suppressing small-scale heterogeneities that are regarded as insignificant. In the third stage, a segmentation algorithm assembles the zones from smaller regions that have been identified in the wavelet analysis.

Source record↗
NASA NTRS2005

Satellite-Derived Management Zones

The term "satellite-derived management zones" (SAMZ) denotes agricultural management zones that are subdivisions of large fields and that are derived from images of the fields acquired by instruments aboard Earth-orbiting satellites during approximately the past 15 years. "SAMZ" also denotes the methodology and the software that implements the methodology for creating such zones. The SAMZ approach is one of several products of continuing efforts to realize a concept of precision agriculture, which involves optimal variations in seeding, in application of chemicals, and in irrigation, plus decisions to farm or not to farm certain portions of fields, all in an effort to maximize profitability in view of spatial and temporal variations in the growth and health of crops, and in the chemical and physical conditions of soils. As used here, "management zone" signifies, more precisely, a subdivision of a field within which the crop-production behavior is regarded as homogeneous. From the perspective of precision agriculture, management zones are the smallest subdivisions between which the seeding, application of chemicals, and other management parameters are to be varied. In the SAMZ approach, the main sources of data are the archives of satellite imagery that have been collected over the years for diverse purposes. One of the main advantages afforded by the SAMZ approach is that the data in these archives can be reused for purposes of precision agriculture at low cost. De facto, these archives contain information on all sources of variability within a field, including weather, crop types, crop management, soil types, and water drainage patterns. The SAMZ methodology involves the establishment of a Web-based interface based on an algorithm that generates management zones automatically and quickly from archival satellite image data in response to requests from farmers. A farmer can make a request by either uploading data describing a field boundary to the Web site or else drawing the boundary on a reference image. Hence, a farmer can start to engage in precision farming shortly after gaining access to the Web site, without the need for incurring the high costs of conventional precision-agriculture data-collection practices that include collecting soil samples, mapping electrical conductivity of soil, and compiling multiyear crop-yield data. Given the boundary of a field, a SAMZ server computes the zones within the field in a three-stage process. In the first stage, a vector-valued image of the field is constructed by assembling, from the archives, the equivalent of a stack of the available images of the field (see figure). In the second stage, the vector-valued image is analyzed by use of a wavelet transform that detects spatial variations considered significant for precision farming while suppressing small-scale heterogeneities that are regarded as insignificant. In the third stage, a segmentation algorithm assembles the zones from smaller regions that have been identified in the wavelet analysis.

Lepoutre, Damien↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗