Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source code”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Development of MHD for NekRS

A spectral-element-based formulation of incompressible MHD is presented in the context of the open-source fluid-thermal code, Nek5000/RS. The formulation supports magnetic fields in a solid domain that surrounds the fluid domain. Several steady-state and time-transient model problems are presented as part of the code verification process. Nek5000/RS is designed for large-scale turbulence simulations, which will be the next step with this new MHD capability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

IEA Wind Energy Task 37 - System Engineering - Aerodynamic Optimization Case Study

This paper presents the results from the first IEA Wind Task 37 aerodynamic optimization case study. 8 participants applied their optimization tools to a purely aerodynamic problem and the results were compared. Overall, the different tools produced widely different designs, while there was better agreement in the improvement achieved. This highlights the fact that further investigation is needed to try to understand these differences and develop best practices. There were too many differences between the different analysis and optimization codes to determine the sources of these differences. However, several potential sources of discrepancy were identified for further investigation. One hypothesis is that one potential source of discrepancy is that the design problem itself is relatively flat in design directions of constant loading. This flatness would impact the convergence in the optimization, while at the same time mean that discrepancies in the design itself would be less severe.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2025

Nek5000/RS performance on advanced GPU architectures

The authors explore performance scalability of the open-source thermal-fluids code, NekRS, on the U.S. Department of Energy's leadership computers, Crusher, Frontier, Summit, Perlmutter, and Polaris. Particular attention is given to analyzing performance and time-to-solution at the strong-scale limit for a target efficiency of 80%, which is typical for production runs on the DOE's high-performance computing systems. Several examples of anomalous behavior are also discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Description and evaluation of the Community Ice Sheet Model (CISM) v2.1

Abstract. We describe and evaluate version 2.1 of the Community Ice Sheet Model (CISM). CISM is a parallel, 3-D thermomechanical model, written mainly in Fortran, that solves equations for the momentum balance and the thickness and temperature evolution of ice sheets. CISM's velocity solver incorporates a hierarchy of Stokes flow approximations, including shallow-shelf, depth-integrated higher order, and 3-D higher order. CISM also includes a suite of test cases, links to third-party solver libraries, and parameterizations of physical processes such as basal sliding, iceberg calving, and sub-ice-shelf melting. The model has been verified for standard test problems, including the Ice Sheet Model Intercomparison Project for Higher-Order Models (ISMIP-HOM) experiments, and has participated in the initMIP-Greenland initialization experiment. In multimillennial simulations with modern climate forcing on a 4 km grid, CISM reaches a steady state that is broadly consistent with observed flow patterns of the Greenland ice sheet. CISM has been integrated into version 2.0 of the Community Earth System Model, where it is being used for Greenland simulations under past, present, and future climates. The code is open-source with extensive documentation and remains under active development.

58 GEOSCIENCES↗
DOE OSTI2017

Pececillo

Pececillo is a proxy-app for the open source Truchas metal processing code (LA-CC-15-097). It implements many of the physics models used in Truchas: free-surface, incompressible Navier-Stokes fluid dynamics (e.g., water waves); heat transport, material phase change, view factor thermal radiation; species advection-diffusion; quasi-static, elastic/plastic solid mechanics with contact; electomagnetics (Maxwell's equations). The models are simplified versions that retain the fundamental computational complexity of the Truchas models while omitting many non-essential features and modeling capabilities. The purpose is to expose Truchas algorithms in a greatly simplified context where computer science problems related to parallel performance on advanced architectures can be more easily investigated. While Pececillo is capable of performing simulations representative of typical Truchas metal casting, welding, and additive manufacturing simulations, it lacks many of the modeling capabilites needed for real applications.

Carlson, Neil↗
DOE OSTI2017

Real-Time Digital Simulation of Inertial Response with Hardware-in-the-Loop Implementation on the CART3 Wind Turbine at the National Wind Technology Center

With increasing penetrations of wind power on electric grids, the stability and reliability of interconnected power systems may be impacted. In some countries that have developed renewable energy sources and systems, grid codes have been revised to require wind power plants (WPPs) to provide ancillary services to support the power system frequency in case of severe grid events. To do this, wind turbine generators (WTGs) should be deloaded to reserve a certain amount of active power for primary frequency response; however, deloading curtails annual energy production, and the market for this type of service needs to be further developed. In this report, we focus on the temporary frequency support provided by WTGs through inertial response. WTGs have potential to provide inertial response, but appropriate control methods should be implemented. With the implemented inertial control methods, wind turbines are capable of increasing their active power output by releasing some of their stored kinetic energy when a frequency excursion occurs. Active power can be temporarily boosted above the maximum power points, after which the rotor speed decelerates, and subsequently an active power output reduction restores the kinetic energy. In this report, we develop two types of models for wind power systems: the first is common, based on the wind power aerodynamic equation, and the power coefficient can be regressed using nonlinear functions; the second is much more complicated, wherein the wind turbine system is modeled using the Fatigue, Aerodynamics, Structures, and Turbulence Modeling (FAST) tool with several degrees of freedoms. A nine-bus test power system is built in Simulink and the Real-Time Digital Simulator, respectively, which are used to evaluate the frequency support performance of the WPPs. We implement two distinct types of inertial control methods in the modeled wind turbines: frequency-based inertial control (FBIC) and stepwise inertial control (SIC). We compare the performances of the two methods in terms of their frequency nadirs, rates of change of frequency, and recovery times. We conclude the results under various wind speeds and penetration cases, which provide insight into designing the inertial response of WTGs. Further, we discuss the impact of the parameters on the performance of the inertial control methods. We evaluate both the scaling factors for the FBIC method and the slope values for the TLIC methods. The simulation work shows the characteristics of different inertial responses compared to conventional synchronous generators. Based on the simulation results, we modify, improve, and test the inertial control methods under a more realistic wind turbine model based on FAST. We then validate the inertial responses under highly turbulent wind conditions generated by TurbSim, and we examine their influences on the turbine mechanical components. The extensive simulation proves the effectiveness of the proposed inertial control methods as well as the nine-bus test power system. We then reconsider the parameters. We rebuild the same test power system using Real time Simulator Computer Aided Design (RSCAD), and we implement the inertial control methods in the real Controls Advanced Research Turbine (CART3), which is prepared for the hardware-in-the-loop field-test simulation. After the setups for the hardware and software hybrid simulation platform are complete, the inertial response is further tested on a real wind turbine for the first time, in which CART3 release the controlled inertial response against the emulated frequency excursion, provided by the real-time simulated power system test bed in RTDS.

17 WIND ENERGY↗
DOE OSTI2018

Recent Updates to the System Advisor Model (SAM)

The System Advisor Model (SAM) is a mature suite of techno-economic models for many renewable energy technologies that can be downloaded for free as a desktop application or software development kit. SAM is used for system-level modeling, including generating performance pro the release of the code as an open source project on GitHub. Other additions that will be covered include the ability to download data directly into SAM from the National Solar Radiation Database (NSRDB) and up- dates to a user-interface macro that assists with PV system sizing. A brief update on SAM's battery model and its integration with the detailed photovoltaic model will also be discussed. Finally, an outline of planned work for the next year will be presented, including the addition of a bifacial model, support for multiple MPPT inputs for detailed inverter modeling, and the addition of a model for inverter thermal behavior.

14 SOLAR ENERGY↗
DOE OSTI2021

IEA Wind Energy Task 37 System Engineering Aerodynamic Optimization Case Study: Preprint

This paper presents the results from the first IEA Wind Task 37 aerodynamic optimization case study. Eight participants applied their optimization tools to a purely aerodynamic problem and the results were compared. Overall, the different tools produced widely different designs, while there was better agreement in the improvement achieved. This highlights the fact that further investigation is needed to try to understand these differences and develop best practices. There were too many differences between the different analysis and optimization codes to determine the sources of these differences. However, several potential sources of discrepancy were identified for further investigation. One hypothesis is that one potential source of discrepancy is that the design problem itself is relatively flat in design directions of constant loading. This flatness would impact the convergence in the optimization, and at the same time mean that discrepancies in the design itself would be less severe.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2021

Progress on Optimizing Wind Farms and Rotor Designs Using Adjoints

Modern wind plants are increasingly tasked with multiple performance objectives. In addition to designing plants that maximize power output and minimize the levelized cost of energy (LCOE), the design and operation of wind plants is increasingly influenced by challenges regarding grid integration of variable generation renewables. This places a growing emphasis on making wind plants more controllable and predictable. WindSE is a Reynolds-averaged Navier-Stokes (RANS) model designed around analytical gradient and adjoint methods, with the ability to capture terrain-induced effects, as shown in Figure 1. The recent addition of an unsteady solver with an actuator line method (ALM) and ongoing work to enable massively parallel optimizations gives it a unique niche to explore coupled plant-level controls and design problems. This code is an open source python package built on the FEniCS framework that utilizes fast, parallel PETSc solvers to model fluid flow throughout wind-farm scale domains. Two recent studies performed using WindSE demonstrate the capability to optimize under a wide variety of flow conditions and objective functions. In the first, we present an optimization focused on modifying the layout of a wind farm with a fixed number of turbines for maximum total power output [1]. This study highlights the ability to quickly perform simulations using the steady Navier-Stokes solver combined with rotors represented as actuator disks while also stressing the importance of capturing terrain-induced effects. Gradient-based optimization using the RANS equations is viable due to the inclusion of efficiently computed adjoint derivatives. We interpret the physical results of the optimal layout and also discuss the computational cost of scaling to larger problems. In the second study, we present the capabilities of the unsteady Navier-Stokes solver, where rotor-blade profiles represented by actuator lines are optimized to enhance wake steering effects and overall power production [2]. We quantify the wind plant performance gains obtained from this type of simultaneous control co-design optimization as compared to optimizing the blade design and yaw independently. Figure 2 shows the differences between a baseline two-turbine system and an optimized system where we fine-tune the blade chord profile. Results and challenges from each study are quickly summarized and used to motivate the current development efforts within WindSE. Current and future work is focused on enabling higher-resolution studies with more degrees of freedom through parallelization of both the simulation and optimization algorithms. We present benchmarking results to show that WindSE performs well in both weak- and strong-scaling tests and further demonstrate that the optimizer obtains the same convergence rates in both shared- and distributed-memory environments. Using larger wind farms, we can study deep-array effects within an optimization context, allowing the use of objective functions that have been previously unstudied. As an example, we present ongoing work on a blockage metric which characterizes the loss of available kinetic energy due to wake effects from multiple upstream turbines.

adjoint optimization↗
DOE OSTI2016

Semi-Analytical Benchmarks for MCNP6

Code verification is an extremely important process that involves proving or disproving the validity of code algorithms by comparing them against analytical results of the underlying physics or mathematical theory on which the code is based. Monte Carlo codes such as MCNP6 must undergo verification and testing upon every release to ensure that the codes are properly simulating nature. Specifically, MCNP6 has multiple sets of problems with known analytic solutions that are used for code verification. Monte Carlo codes primarily specify either current boundary sources or a volumetric fixed source, either of which can be very complicated functions of space, energy, direction and time. Thus, most of the challenges with modeling analytic benchmark problems in Monte Carlo codes come from identifying the correct source definition to properly simulate the correct boundary conditions. The problems included in this suite all deal with mono-energetic neutron transport without energy loss, in a homogeneous material. The variables that differ between the problems are source type (isotropic/beam), medium dimensionality (infinite/semi-infinite), etc.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Preliminary Criticality and Radiation Shielding Analysis for the Storage and Transfer of MARVEL Reactor Spent Nuclear Fuel

This report documents the results of preliminary nuclear criticality and radiation shielding assessments during transfer and dry storage of irradiated Microreactor Applications Research Validation and Evaluation (MARVEL) reactor fuel at Idaho National Laboratory (INL). The assessments focus on transfer casks and storage canisters that are currently in use at INL, which may be compatible with the irradiated MARVEL reactor fuel. The criticality assessments were performed using the radiation transport code MCNP6 with 37 MARVEL reactor fuel elements in various configurations and scenarios. All transfer and storage configurations under dry conditions were below the assumed criticality safety limit of 0.93. Some storage and transfer configurations under wet conditions exceeded the criticality safety limit. This suggests that the appropriate administrative and engineering controls, in addition to reducing the number of MARVEL reactor fuel elements per transfer cask or storage canister, can be expected to ensure criticality safety under all scenarios. The radiation shielding assessments were performed by generating conservative neutron and photon source terms using the ORIGEN module in the SCALE suite of codes. These source spectra were used to estimate the dose equivalent rates using the radiation transport code MCNP6, both on contact and 1 m away from the fuel and transfer casks. The maximum estimated dose equivalent rate of 37 unshielded MARVEL reactor fuel elements on contact is approximately 42000 R/hr. The maximum estimated dose equivalent rates on contact to the ATR transfer cask, HFEF-5 transfer cask, and high load charger were approximately 233 mR/hr, 171 mR/hr, and 201 mR/hr, respectively. This suggests that with the appropriate administrative and engineering controls, all three transfer casks analyzed can be expected to provide sufficient radiation shielding to workers during transfer of irradiated MARVEL reactor fuel. These calculations are performed to support the planning and strategy for the MARVEL project and will demonstrate the technical viability of the different configurations discussed and help identify where engineered or administrative controls may be necessary. A complete criticality safety analysis and radiation shielding analysis, including validation and contingency and accident analysis must be completed by licensed and authorized personnel before any transfer or storage of MARVEL reactor nuclear fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Alternative modeling methods for plasma-based Rf ion sources

Rf-driven ion sources for accelerators and many industrial applications benefit from detailed numerical modeling and simulation of plasma characteristics. For instance, modeling of the Spallation Neutron Source (SNS) internal antenna H{sup −} source has indicated that a large plasma velocity is induced near bends in the antenna where structural failures are often observed. This could lead to improved designs and ion source performance based on simulation and modeling. However, there are significant separations of time and spatial scales inherent to Rf-driven plasma ion sources, which makes it difficult to model ion sources with explicit, kinetic Particle-In-Cell (PIC) simulation codes. In particular, if both electron and ion motions are to be explicitly modeled, then the simulation time step must be very small, and total simulation times must be large enough to capture the evolution of the plasma ions, as well as extending over many Rf periods. Additional physics processes such as plasma chemistry and surface effects such as secondary electron emission increase the computational requirements in such a way that even fully parallel explicit PIC models cannot be used. One alternative method is to develop fluid-based codes coupled with electromagnetics in order to model ion sources. Time-domain fluid models can simulate plasma evolution, plasma chemistry, and surface physics models with reasonable computational resources by not explicitly resolving electron motions, which thereby leads to an increase in the time step. This is achieved by solving fluid motions coupled with electromagnetics using reduced-physics models, such as single-temperature magnetohydrodynamics (MHD), extended, gas dynamic, and Hall MHD, and two-fluid MHD models. We show recent results on modeling the internal antenna H{sup −} ion source for the SNS at Oak Ridge National Laboratory using the fluid plasma modeling code USim. We compare demonstrate plasma temperature equilibration in two-temperature MHD models for the SNS source and present simulation results demonstrating plasma evolution over many Rf periods for different plasma temperatures. We perform the calculations in parallel, on unstructured meshes, using finite-volume solvers in order to obtain results in reasonable time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Setting the H-mode pedestal structure: variations of particle source location using gas puff and pellet fueling

Experiments in DIII-D show that a particle source location inside the top of the H-mode pedestal (pellets) maintains a higher pedestal pressure than an edge source (gas fueling) through a widening of the electron temperature pedestal with reduction of the temperature gradient. The effect of these two fueling schemes on the H-mode pedestal structure was examined in DIII-D by comparing controlled pellet-fueled and gas-fueled discharges across a fueling scan up to 40 torr 1 s –1 . High resolution electron profiles reveal that gas fueling lowers the pedestal pressure as the density profile shifts radially outwards and the separatrix density increases, while pellet fueling maintains a constant pedestal pressure. The neutral source locations from pellets and gas are determined with the PELLET and UEDGE codes, respectively, and quantify the particle source localization. Pellets provide significant ionization inside the pedestal top while gas puffing localizes ionization in the scrape-off-layer and pedestal foot, broadly consistent with the density profile structure influenced by the source. ELMs are observed to increase in frequency and reduce impurity content as fueling is increased. Stability analysis with ELITE shows that both conditions are near the type-I ELM corner of the peeling–ballooning stability diagram, which is altered significantly by the introduction of pellets. Since transport mechanisms are not observed to change substantially with particle source location, wider pedestals allow the pellet-fueled discharges to retain higher pedestal temperatures at similar pedestal densities. EPED1 is tested to capture the pedestal pressure, under-predicting the height with pellets and over-predicting the height with gas by ~15%. Overall, these results have important implications for future reactors where pellet fueling will be the primary particle source due to an opaque scrape-off-layer by showing that the neutral source location plays a role in setting the structure of the H-mode pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fixed Source Sensitivity Calculations for Inertial Confinement Fusion Applications

A numerical code library was developed for the radiation transport code MCNP6.3 to calculate generalized response sensitivity coefficients for fixed source neutron transport problems with applications to inertial confinement fusion (ICF) experiments. The new MCNP6.3 dependency is used to generate a novel time convolution response that represents a neutron time-of-flight (nToF) signal. The traditional suite of macroscopic cross-section sensitivities and constrained fixed source probability distribution sensitivities are available for both the standard and the new response tallies in this library. However, novel sensitivity coefficients for the constrained hyperparameters of analytic fixed source probability distributions are emphasized in this work for their connection to ICF neutron transport models. Particularly, advanced Monte Carlo methods are developed for calculating the sensitivity of a nToF signal to perturbations in an ICF plasma’s ion temperature and burn history as well as perturbations in the target liner mass density and the shape parameters of the nToF detector’s impulse response function. Together, these capabilities form an advanced suite of computational tools that can be used to analyze and extract information from any ICF experimental platform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Integrated Decision Methodology for Final Disposal of Radioactive Sources

Radioactive sources are used throughout the world, and their beneficial uses are many. Sources are at the core of many applications in different branches of industry, medicine, agriculture and research. Sources are present in a very wide range of equipment that is used for, inter alia, cancer treatment, killing bacteria in food, sterilizing medical supplies, measuring instruments such as gauges used to measure soil moisture and soil density and many other components, irradiating seeds for enhancing food production, protecting buildings from lightning strikes, mapping underground sources of water, prospecting for oil and gas reserves, measuring density of soil for construction projects, or even detecting smoke. At their “cradle” and when in use, radioactive sources are usually properly managed and controlled. It is when a source has reached the end of its useful life when they are at a higher risk of being under inadequate control, poorly managed, or even becoming orphan. Sources that have reached the end of their useful life must be carefully and safely managed. The “grave” of disused sources needs a proper disposal process to prevent them from posing a potential threat to people or the environment. To protect the public and the environment from the potential hazards of ionizing radiation, and to prevent disused sources from becoming orphan, a “cradle-to-grave” control of radioactive sources is essential, as promoted by the Code of Conduct on the Safety and Security of Radioactive Sources. This approach requires a national policy and strategy, an adequate legal and regulatory framework, and adequate resources and infrastructure. Many IAEAMember States have expressed the need for effective solutions to this challenge. There is also a gap related to the life cycle costs of potential disposal practices. This paper proposes a multi-attribute decision making methodology to compare disposal options using both economic and non-economic criteria leading to a preferred option for disposal and a strategy for final disposal of disused radioactive sources that is specific and optimized for each IAEA Member State.

Beatty, Randy↗
DOE OSTI2021

Introducing PRIMRE's MRE Software Knowledge Hub (February 2021)

This paper focuses on the role of the Marine Renewable Energy (MRE) Software Knowledge Hub on the Portal and Repository for Information on Marine Renewable Energy (PRIMRE). The MRE Software Knowledge Hub provides online services for MRE software users and developers, and seeks to develop assessments and recommendations for improving MRE software in the future. Online software discovery platforms, known as the Code Hub and the Code Catalog, are provided. The Code Hub is a collection of open-source MRE software that includes a landing page with search functionality, linked to files hosted on the MRE Code Hub GitHub organization. The Code Catalog is a searchable online platform for discovery of useful (open-source or commercial) software packages, tools, codes, and other software products. To gather information about the existing MRE software landscape, a software survey is being performed, the preliminary results of which are presented herein. Initially, the data collected in the MRE software survey will be used to populate the MRE Software knowledge hub on PRIMRE, and future work will use data from the survey to perform a gap analysis and develop a vision for future software development. Additionally, as one of PRIMRE's roles is to support development of MRE software within project partners, a silo of knowledge relating to best practices has been gathered. An early draft of new guidance developed from this knowledge is presented.

gap analysis↗
DOE OSTI2019

Time-Encoded Gamma-Ray Imaging Using a 3-D Position-Sensitive CdZnTe Detector Array

The gamma-ray response of thick pixelated CdZnTe detectors varies spatially but is relatively constant as a function of time. Detector irregularities arise from flaws in the material, electrode fabrication, and gaps between the detectors in an array. These nonuniformities generate artifacts in coded aperture images that can obscure a source of interest. This article presents the theory, design, and demonstration of a 3-D CdZnTe-based time-encoded imager. Experimental results show the gamma-ray images of point-like and extended 57 Co and HEU-metal sources formed using a spatially nonuniform 3-D CdZnTe detector array and moving coded aperture mask. Image signal-to-noise ratio (SNR) for a point source was within 10% of its statistical limit at 250k counts.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Demonstration of a code coupling framework for modeling beam-collimator impacts in the advanced photon source

The high-brightness beams being produced in current and future accelerators present new machine protection concerns with the potential for high-energy-density (HED) conditions ( >100 J/mm 3 ) in beam-intercepting components. Simulating HED conditions in accelerators requires utilizing a suite of physics codes for particle dynamics, particle-matter interactions, and hydrodynamics. This paper describes a method of coupling the codes elegant, fluka, and flash to simulate the effects of a rapid beam loss in the advanced photon source storage ring and the resulting interaction of the beam and collimators. This paper expands previous work [J. Dooling et al., Collimator irradiation studies at the advanced photon source, in Proceedings of the IBIC-2023 (2023), pp. 245–249] by introducing a definition of the evolving geometry of the collimator surface as well as providing methods for simulating the absorption of synchrotron radiation and tracking shower particles produced during beam strikes. We demonstrate this framework by simulating machine conditions of the APS ring before and after its recent upgrade. Simulation results are compared with observed damage to collimators and test samples taken from the APS ring.

Accelerator/storage ring control systems↗