Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power engineering computing”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Models implemented in the methodological approach to design the initial STEP first wall contour

The official Spherical Tokamak for Energy Production mission aims to demonstrate the ability to generate net electricity from fusion with the STEP Prototype Power plant. One of the key technological and engineering challenges in fusion power plants is managing the loads on the first wall within acceptable limits. Therefore, the conceptual design development of the STEP Prototype Power plant needs to be based on load estimates derived using legitimate plasma physics assumptions through dynamic and flexible tools. The current design foresees the STEP main chamber first wall to withstand steady-state heat loads of up to ~1 MW/m 2 , excluding critical regions expected to receive higher heat loads such as the baffle regions approaching the divertors. These critical areas will require ad hoc assessments and will be designed with the presence of limiters. This article focuses on the models and methodology adopted for designing the 2-D poloidal contour of the STEP first wall, based on the anticipated charged particle and radiation heat loads during normal operation. Firstly, the models adopted for calculating the charged particle and radiation heat loads are introduced. The first model is validated through benchmarking against the particle tracing code SMARDDA, while the second model is verified by comparing it with data from the MAST-U experiment. Secondly, the model used to design the 2-D first wall contour according to the heat loads is explained. We acknowledge that this preliminary design stage assumes certain simplifications, notably an axisymmetric geometry, for computational efficiency and clarity in presentation. It is understood that subsequent design phases will address the complexities of real-world engineering, including non-axisymmetric effects, transient plasma scenarios, and the impact of disruptions on the first wall design. Finally, an automatic procedure based on these models is presented for defining the 2-D poloidal contour of the STEP first wall to minimize heat loads, taking into account the need to radiate most of the alpha-particle and auxiliary heating power. Here, by providing an overview of the models, methodology, and an automatic procedure, this paper contributes to the design process of the STEP first wall, addressing the engineering challenges associated with fusion power plant development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An Adaptive Virtual Oscillator Control Structure for Grid-Forming Inverters

The electrical grid is facing unprecedented challenges due to the increasing penetration of inverter-based resources. Grid forming inverters (GFMIs) are a promising technology to address these challenges. Recently, the virtual oscillator based GFMI control is attracting more attention due to its superior performance over other control strategies. In this article, an adaptive control strategy is proposed to provide flexible operation and transition between grid-connected and islanded modes. The controller adapts the virtual oscillator's parameter values depending on the operation mode. It also provides a feedback signal using a measured frequency to account for any steady-state errors and to allow a seamless transition from grid-connected to islanded mode. Finally, to show the feasibility of the proposed controller, this article discusses the simulation results from the implementation of the controller on a single inverter system and on a group of inverters on a large practical utility feeder, the IEEE 13 node feeder, using the DIgSILENT simulation environment.

adaptive control↗
DOE OSTI2023

Visualizing Fault Induced Traveling Waves In Medium Voltage Systems

Traveling waves are induced in power systems during most transient events in the grid. These waves travel close to the speed of light in overhead lines and 50% to 60% the speed of light in underground cables. Even though traveling wave-based protection schemes for transmission systems are available commercially, traveling waves in medium voltage distribution networks are still in research space. Compared to transmission system, medium voltage distribution systems contain more reflections and refractions. Thus, visualization is challenging and is critical in locating faults in distribution network. To address this visualization challenge, this paper presents an open-source tool to visualize the traveling waves using Bewley lattice approach. The developed visualization tool will be useful for the protection engineers to detect and triangulate fault locations in the medium voltage systems and isolate the faults.

Bewley Lattice↗
DOE OSTI2025

Growth and characterization of high-quality Zr doped AlN epilayers

AlN stands out for its remarkable figures of merit for electronic and photonic devices, attributed to its ultrawide bandgap of ∼6.1 eV and an exceptionally high critical field of ∼15 MV/cm. More recently, zirconium (Zr) doped AlN (AlN:Zr) has also been identified as a promising material platform for the exploration of solid-state qubits for quantum information and technology, high performance piezoelectric acoustic wave resonators, and optically triggered ultrafast power switching devices facilitated by optically activating Zr related impurities. Despite the significant potential, the ability for producing AlN:Zr epitaxial structures has yet to be established. In this study, we have achieved AlN:Zr epilayers with a high Zr doping level [NZr] of up to 1020 cm−3 using industrial standard metal-organic chemical vapor deposition growth technique. High crystalline quality of AlN:Zr was confirmed by x-ray diffraction, revealing a narrow full width at half maximum of the (002) rocking curve at 216 arcsec for 1.8 μm thick epilayers deposited on sapphire at [NZr]=1020 cm−3. Zr doping was observed to slightly increase the c-lattice constant to 4.992 Å for AlN:Zr (at [NZr]=1020 cm−3) compared to 4.980 Å for undoped AlN. X-ray photoelectron spectroscopy measurement results verified the substitution of Zr at the Al site (ZrAl). The formation of (ZrAl–VN) complexes, which are predicted to possess all the desired properties required by quantum qubits, was confirmed through optical absorption studies. The realization of high-quality AlN:Zr epilayers significantly broadens the scope of technologically significant device applications for AlN.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electrical circuit control in power systems

Electrical circuit control techniques in power systems are disclosed herein. In one embodiment, a supervisory computer in the power system can be configured to fit phasor measurement data from phasor measurement units into a Gaussian distribution with a corresponding Gaussian confidence level. When the Gaussian confidence level of the fitted Gaussian distribution is above a Gaussian confidence threshold, the supervisory computer can be configured to perform an ambient analysis on the received phasor measurement data to determine an operating characteristic of the power system. The supervisory computer can then automatically applying at least one electrical circuit control action to the power system in response to the determined operating characteristic.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Identifying limitations in screening high-throughput photocatalytic bimetallic nanoparticles with machine-learned hydrogen adsorptions

The Sabatier principle is of fundamental importance to computational catalyst discovery, saving researchers time and expense by predicting catalytic activity in silico at scale. However, as polycrystalline and nanoscale catalysts increasingly dominate industry, computational screening tools must be adapted to these uses. In this work, we demonstrate the effectiveness of computational adsorption energy screening in nanocatalysis by comparing a multisite adsorption energy prediction workflow against a large experimental dataset of hydrogen evolution activities over bimetallic nanoparticles. Comparing 16 million hydrogen adsorption energy predictions with the hydrogen evolution activity of 5300 experiments across 84 monometallic and bimetallic systems, we discover that favorable adsorption energies are a necessary condition for experimental activity, but other factors often determine trends in practice. Finally, about half of the bimetallic search space can be excluded from experimental screens using hydrogen adsorption predictions, but this method may become significantly more powerful when combined with other screening tools.

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

SEAS Communication Engine: An Extensible, Flexible Wrapper for Co-Simulation Agents

When modeling and analyzing the power grid and other large scale systems, researchers often express scenarios as optimization problems and feed them into advanced software solvers. In order to allow multiple solvers to communicate with each other and share data from different domains, the National Renewable Energy Laboratory (NREL) and associated Department of Energy (DOE) labs have developed a software framework called the Hierarchical Engine for Large-scale Infrastructure Co-Simulation (HELICS). HELICS allows cosimulation via a collection of client libraries for different languages that can be called from the appropriate optimization software. However, these client libraries do not provide a higher level of abstraction beyond reading and writing data off of the shared HELICS bus. In this paper, we describe a new software library called the SEAS Communication Engine that exposes a higher-level API for running cosimulation problems. The SEAS Engine provides a class-based abstraction on top of the Python HELICS client, in order to allow users to implement their domain-specific cosimulations without needing to interact with core HELICS primitives. This will make adoption of HELICS and cosimulation in general easier, by exposing a simpler API. In the second part of the paper, we validate our library on a collection of different simulation examples, including the canonical IEEE 13 Bus Feeder. Lastly, we demonstrate using the SEAS Engine to directly call domain-specific code written in the Julia programming language. Our hope is that this will serve as a template for easily calling software in different programming languages via the SEAS Engine, thereby avoiding code duplication and complexity.

co-simulation↗
DOE OSTI2025

Hexagonal Boron Nitride: Physical Properties, Hydride Vapor‐Phase Epitaxy Growth of Large‐Diameter Quasi‐Bulk Wafers and Applications

Hexagonal boron nitride ( h ‐BN), with its ultrawide bandgap and 2D structure, holds an immense promise for advanced semiconductor applications. Scaling bulk crystals to large‐diameter wafers, crucial for complex device fabrication, remains a challenge with high temperature, high pressure, and metal flux solution methods. To address this, recent efforts have focused on hydride vapor‐phase epitaxy (HVPE) for producing large diameter thick h ‐BN quasi‐bulk wafers, reaching hundreds of micrometers. These HVPE‐grown quasi‐bulk crystals exhibit excellent c ‐axis long‐range order. Notably, the in‐plane mobility‐lifetime products for both electrons and holes surpass 10 −4 cm 2 V −1 , which are two orders of magnitude greater than the out‐of‐plane (vertical) values, highlighting the potential for high‐performance devices leveraging superior lateral transport. Lateral detectors fabricated from 100 μm thick B‐10 enriched h‐ BN wafers have achieved a record 60% thermal neutron detection efficiency. Based on its physical properties, h‐ BN appears to be an outstanding material of choice for light‐triggered electronic power switches capable of supporting high ‐ voltage and high ‐ power operations. These recent advancements in large‐diameter h ‐BN quasi‐bulk crystal growth, enabled by HVPE, pave the way for applications spanning deep UV photonics, high‐power electronics, high‐efficiency neutron detection, and quantum information technologies, establishing h ‐BN as both a versatile active semiconductor and an ideal substrate.

Jiang, Hongxing [Department of Electrical and Comp↗
DOE OSTI2023

Snapshot multispectral imaging using a diffractive optical network

Abstract Multispectral imaging has been used for numerous applications in e.g., environmental monitoring, aerospace, defense, and biomedicine. Here, we present a diffractive optical network-based multispectral imaging system trained using deep learning to create a virtual spectral filter array at the output image field-of-view. This diffractive multispectral imager performs spatially-coherent imaging over a large spectrum, and at the same time, routes a pre-determined set of spectral channels onto an array of pixels at the output plane, converting a monochrome focal-plane array or image sensor into a multispectral imaging device without any spectral filters or image recovery algorithms. Furthermore, the spectral responsivity of this diffractive multispectral imager is not sensitive to input polarization states. Through numerical simulations, we present different diffractive network designs that achieve snapshot multispectral imaging with 4, 9 and 16 unique spectral bands within the visible spectrum, based on passive spatially-structured diffractive surfaces, with a compact design that axially spans ~72 λ m , where λ m is the mean wavelength of the spectral band of interest. Moreover, we experimentally demonstrate a diffractive multispectral imager based on a 3D-printed diffractive network that creates at its output image plane a spatially repeating virtual spectral filter array with 2 × 2 = 4 unique bands at terahertz spectrum. Due to their compact form factor and computation-free, power-efficient and polarization-insensitive forward operation, diffractive multispectral imagers can be transformative for various imaging and sensing applications and be used at different parts of the electromagnetic spectrum where high-density and wide-area multispectral pixel arrays are not widely available.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tools Assessing Performance (FY2022 Quarter 4 Report)

We were finally successful in acquiring TTU’s 200-meter tower data for the duration of the Isolated Building Wake Experiment (IBWE). This data will provide an independent measurement of the undisturbed wind profile including heights far above those that are typical of meteorological sensor networks. These measurements will characterize the mean and turbulence throughout the atmospheric surface layer, which will be particularly useful for validation of highfidelity computational fluid dynamics models. This dataset can be shared freely among the TAP team but not outside of it as TTU sells this data to the wind power industry. The journal article overview has been submitted to LANL’s publication review and will be submitted to the Journal of Wind Engineering and Industrial Aerodynamics once the review is complete. We are nearing completion on another two articles. The first discusses analyses of the evolution of the measurements of mean wind speed along the centerline of the building as a function of atmospheric stability. Rather than using individual time steps, all of the wind speed measurements are normalized using the upwind measurement at 2 building heights (H) and filtered by wind direction and atmospheric stability. The second discusses the development of the diffusive wake model and uses the averaged profiles from the first article to validate the behavior of the model. The portion of the section below on the diffusive wake model discussing comparisons with the averaged profiles, will be included in the validation section of this paper.

17 WIND ENERGY↗
DOE OSTI2020

MULTI-DIMENSIONAL MODELING OF THE CFR ENGINE FOR THE INVESTIGATION OF SI NATURAL GAS COMBUSTION AND CONTROLLED END-GAS AUTOIGNITION

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by low emission limits and the wide range of chemical reactivity in pipeline quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its methane number (analogous to the octane number for liquid SI fuels) can span from 65 - 95. In order to realize diesel-like efficiencies, SI natural gas engines must be designed to operate at high BMEP near knock limits over a wide range of fuel reactivity. This requires a deep understanding regarding the combustion-engine interactions pertaining to flame propagation and end-gas autoignition (EGAI). However, EGAI, if controlled, provides an opportunity to increase SI natural gas engine efficiency by increasing combustion rate and the total burned fuel, mitigating the effects of the slow flame speeds of natural gas fuels which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, in order to study EGAI phenomenon, the present work highlights multi-dimensional computational fluid dynamics (CFD) models of the Cooperative Fuel Research (CFR) engine. The CFR engine models are used to investigate fuel-engine interactions that lead to EGAI with natural gas, including effects of fuel reactivity, engine operating parameters, and exhaust gas recirculation (EGR). A Three-Pressure Analysis, performed with GT-Power, was used to estimate initial and boundary conditions for the three-dimensional CFD model. CONVERGE CFD v2.4 was used for the three-dimensional CFD modeling where the level set G-Equation model and SAGE detailed chemical kinetics solver were used. An assessment of the different modeling approaches is also provided to evaluate their limitations, advantages and disadvantages, and for which situations they are most applicable. Model validation was performed with experimental data taken with a CFR engine over varying compression ratio, CA50, EGR fraction, and IMEP and shows good agreement in Peak Cylinder Pressure (PCP), PCP crank angle, and the location of the 10%, 50%, and 90% mass fraction burned (CA10, CA50, and CA90, respectively). The models can predict the onset crank angle and pressure rise rate for light, medium, and heavy EGAI under a variety of fuel reactivities and engine operating conditions.

Bestel, Diego↗
DOE OSTI2025

Fabrication and Testing of Solid-Solution Strengthened Corrosion Resistant Alloys For Service in Molten Fluoride Environments

The demand for higher system thermal efficiencies requires the operation of power generation cycles and heat conversion systems at progressively higher temperatures. As the system operating temperature increases, existing materials may not provide adequate mechanical properties or environmental compatibility or both. There is an increasing commercial interest in the development and deployment of liquid-fueled Molten Salt Reactors (MSRs). Hastelloy®N, the highest performing candidate MSR structural alloy, is not capable of operations at temperatures above 700°C, thus limiting the performance of these systems. Using an Integrated Computational Materials Engineering (ICME)-approach and small laboratory scale heats, ORNL developed a class of patented alloys covered by U.S. Patent 9, 435, 011 B2, “Creep-resistant, Cobalt-free alloys for high temperature, liquid-salt heat exchanger systems,” similar to Hastelloy®N in that they are primarily solid solution strengthened. In contrast to precipitation strengthened alloys, the microstructure of solid solution alloys and hence the high temperature mechanical properties are stable for extended periods of time allowing long reactor operating life. The new alloys have shown to possess good resistance to liquid fluorides at temperatures up to 850°C and have significantly improved creep properties when compared to Hastelloy®N. The purpose of the CRADA project was for ORNL to collaborate with Haynes International- a materials producer, MetalTek International- a foundry, and Kairos Power – an advanced reactor developer – to scale-up selected alloys, evaluate their properties, and identify one solid solution strengthened alloy that can meet the property requirements for the reactor being developed by Kairos Power and other similar liquid fluoride-salt cooled reactors. As part of the project, eight alloys were down-selected and fabricated in larger industrial scale heats by Haynes International. Resistance to molten salt was evaluated in flowing FLiNaK and FLiBe by Kairos Power using their Rotating Cage Loop (RCL) system. Accounting for iron deposition during these tests, the new alloys displayed very low net mass change showing excellent corrosion performance in molten salt. Creep properties evaluated at ORNL were found to be better than that of Hastelloy®N and 316 stainless steel. Long-term stabilities of the alloys evaluated by Haynes International showed that these alloys have excellent thermal stability in the temperature range 704.4-815.6°C, with the change in strength and ductility being less than 10-15% after a 4000-hour exposure at 815.6°C. Autogenously Gas Tungsten Arc Welding (GTAW) welded samples showed less than 10% change in yield strength / ultimate tensile strength / total elongation compared to the basemetal, indicating that the alloys have excellent weldability. Three parts were successfully investment-cast using one alloy with very little voiding showing feasibility of fabricating parts using the casting process. This project enabled extensive interaction between the material producer Haynes International, casting supplier MetalTek, and reactor developer Kairos Power. This facilitated testing of materials and components produced using the newly developed alloys by the end-user. This allowed the generation of critical dataset required for down-selection of a few promising alloys for further development. This data is also currently being shared with other reactor designers for them to evaluate the suitability of this alloy for their reactor design. The availability of this alloy will ultimately enable the design and development and deployment of MSRs with increased temperature of operation and thus, improved efficiencies.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Fabrication and Testing of Solid-Solution Strengthened Corrosion Resistant Alloys For Service in Molten Fluoride Environments

The demand for higher system thermal efficiencies requires the operation of power generation cycles and heat conversion systems at progressively higher temperatures. As the system operating temperature increases, existing materials may not provide adequate mechanical properties or environmental compatibility or both. There is an increasing commercial interest in the development and deployment of liquid-fueled Molten Salt Reactors (MSRs). Hastelloy®N, the highest performing candidate MSR structural alloy, is not capable of operations at temperatures above 700°C, thus limiting the performance of these systems. Using an Integrated Computational Materials Engineering (ICME)-approach and small laboratory scale heats, ORNL developed a class of patented alloys covered by U.S. Patent 9,435,011 B2, “Creep-resistant, Cobalt-free alloys for high temperature, liquid-salt heat exchanger systems,” similar to Hastelloy®N in that they are primarily solid solution strengthened. In contrast to precipitation strengthened alloys, the microstructure of solid solution alloys and hence the high temperature mechanical properties are stable for extended periods of time allowing long reactor operating life. The new alloys have shown to possess good resistance to liquid fluorides at temperatures up to 850°C and have significantly improved creep properties when compared to Hastelloy®N. The purpose of the CRADA project was for ORNL to collaborate with Haynes International- a materials producer, MetalTek International- a foundry, and Kairos Power – an advanced reactor developer – to scale-up selected alloys, evaluate their properties, and identify one solid solution strengthened alloy that can meet the property requirements for the reactor being developed by Kairos Power and other similar liquid fluoride-salt cooled reactors. As part of the project, eight alloys were down-selected and fabricated in larger industrial scale heats by Haynes International. Resistance to molten salt was evaluated in flowing FLiNaK and FLiBe by Kairos Power using their Rotating Cage Loop (RCL) system. Accounting for iron deposition during these tests, the new alloys displayed very low net mass change showing excellent corrosion performance in molten salt. Creep properties evaluated at ORNL were found to be better than that of Hastelloy®N and 316 stainless steel. Long-term stabilities of the alloys evaluated by Haynes International showed that these alloys have excellent thermal stability in the temperature range 704.4-815.6°C, with the change in strength and ductility being less than 10-15% after a 4000-hour exposure at 815.6°C. Autogenously Gas Tungsten Arc Welding (GTAW) welded samples showed less than 10% change in yield strength / ultimate tensile strength / total elongation compared to the basemetal, indicating that the alloys have excellent weldability. Three parts were successfully investment-cast using one alloy with very little voiding showing feasibility of fabricating parts using the casting process. This project enabled extensive interaction between the material producer Haynes International, casting supplier MetalTek, and reactor developer Kairos Power. This facilitated testing of materials and components produced using the newly developed alloys by the end-user. This allowed the generation of critical dataset required for down-selection of a few promising alloys for further development. This data is also currently being shared with other reactor designers for them to evaluate the suitability of this alloy for their reactor design. The availability of this alloy will ultimately enable the design and development and deployment of MSRs with increased temperature of operation and thus, improved efficiencies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Techno-Economic Analysis of CSP Incorporating sCO2 Brayton Power Cycles: Trade-Off Between Cost and Performance

Concentrating solar power (CSP) plants, thanks to the implementation of cost-competitive thermal energy storage, represent a dispatchable zero-emission alternative to traditional fossil fuel power plants. Next generation solar towers are expected to adopt high temperature receivers (>700 degrees C) coupled to sCO2-based power blocks, which optimal design is generally pushed towards the maximum cycle efficiency, often neglecting the economic impact with the justification that the main share of the capital cost is represented by the heliostat field. As result, the scientific literature lacks in comprehensive studies on techno-economic evaluation of CSP+sCO2 power plants addressing the important correlation that exists between system cost and performances. This work provides a preliminary techno-economic analysis of a solar power tower comparing four different cycle configurations for the sCO2 power block. Results have been reported on a Pareto front, highlighting the tradeoff between the plant investment cost and the solar-to-electricity plant efficiency. The trends of the optimization variables and cycle results have been reported to give useful insights about proper assumptions for the sCO2 power block design. The recompressed cycle with intercooling resulted as the most promising configuration and it has been further analyzed through a comparison of different solutions on the Pareto front. The cost breakdown of the sCO2 power block has been reported to highlight which components have the greatest impact on the overall plant cost and how they vary along the optimal solutions front. Eventually, the optimization has been repeated introducing a correlation to compute the turbomachinery isentropic efficiencies, to investigate their effect on the techno- economic analysis.

concentrated solar power↗
DOE OSTI2020

MASTODON: An Open-Source Software for Seismic Analysis and Risk Assessment of Critical Infrastructure

Seismic analysis and risk assessment of safety-critical infrastructure like hospitals, nuclear power plants, dams, and facilities handling radioactive materials involve computationally intensive numerical models and coupled multiphysics scenarios. They are also performed in a strict regulatory environment that requires high software quality assurance standards, and in the case of safety-related nuclear facilities, a conformance to the American Society of Mechanical Engineers Nuclear Quality Assurance (NQA-1) standard. This paper introduces the open-source finite-element software, MASTODON (Multi-hazard Analysis of Stochastic Time-Domain Phenomena), which implements state-of-the-art seismic analysis and risk assessment tools in a quality-controlled environment. MASTODON is built on MOOSE (Multi-physics Object-Oriented Simulation Environment), which is a highly parallelizable, NQA-1 conforming, coupled multiphysics, finite-element framework developed at Idaho National Laboratory. MASTODON is capable of fault rupture and source-to-site wave propagation using the domain reduction method, nonlinear site response, and soil-structure interaction analysis, implicit and explicit time integration, automated stochastic simulations, and seismic probabilistic risk assessment. When coupled with other MOOSE applications, MASTODON can also solve strongly and weakly coupled multiphysics problems. This paper presents a summary of the capabilities of MASTODON and some demonstrative examples.

42 ENGINEERING↗
DOE OSTI2024

Estimating the carbon footprint of computational fluid dynamics

In this study, computational resources have grown exponentially in the past few decades. These machines make possible research and design in fields as diverse as medicine, astronomy, and engineering. Despite ever-increasing computational capabilities, direct simulation of complex systems has remained challenging owing to the degrees of freedom involved. At the cusp of exascale computing, high-resolution simulation of practical problems with minimal model assumptions may soon experience a renaissance. However, growing reliance on modern computers comes at the cost of a growing carbon footprint. To illustrate this, we examine historic computations in fluid dynamics where larger computers have afforded the opportunity to simulate flows at increasingly relevant Reynolds numbers. Under a variety of flow configurations, the carbon footprint of such simulations is found to scale roughly with the fourth power of Reynolds number. This is primarily explained by the computation cost in core-hours, which is also described by similar scaling, though regional differences in renewable energy use also play a role. Using the established correlation, we examine a large database of simulations to develop estimates for the carbon footprint of computational fluid dynamics in a given year. Collectively, the analysis provides an additional benchmark for new computations where, in addition to balancing considerations of model fidelity, carbon footprint should also be considered.

42 ENGINEERING↗
DOE OSTI2024

ParaStell: parametric modeling and neutronics support for stellarator fusion power plants

The three-dimensional variation inherent to stellarator geometries and fusion sources motivates three-dimensional modeling to obtain accurate results from computational modeling in support of design and analysis of first wall, blanket, and shield (FWBS) systems. Manually constructing stellarator fusion power plant geometries in computer-aided design (CAD) and defining the corresponding fusion source can be cumbersome and challenging. The open-source parametric modeling toolset ParaStell has been developed to automate construction of such geometries in low-fidelity. Low-fidelity modeling is useful during the conceptual phase of engineering design as a means of rapidly exploring the design space of a given device. The modeling capability of ParaStell includes in-vessel components and magnets, for any given stellarator configuration, using a parametric definition and plasma equilibrium data. Furthermore, the toolset automates the generation of detailed, tetrahedral neutron source definitions and DAGMC geometries for use in neutronics modeling. ParaStell assists rapid design iteration, parametric study, and design optimization of stellarator fusion cores. As a demonstration of the design iteration capability, the effect of the three-dimensional parameter space on tritium breeding and magnet shielding is investigated, using the WISTELL-D configuration as a design basis. Blanket and shield thicknesses are varied in three dimensions, using the space available between the plasma edge and magnet coils as a constraint. The corresponding effects on tritium breeding ratio and magnet heating are tallied using the open-source Monte Carlo particle transport code OpenMC. The inclusion of additional and higher-fidelity modeling capabilities is planned for ParaStell’s future, as well as its implementation in machine-driven optimization.

Moreno, Connor↗