SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Abstract not provided.
In this report, the activities and achievements made by Argonne National Laboratory for the Nuclear Energy Advanced Modeling and Simulation (NEAMS) BISON code metallic fuel validation and verification project in FY2022 are summarized. The FIPD-BISON integration powered metallic fuel low-burnup swelling framework has been enhanced to cover both radial and axial swelling analyses. Axial-dependent as-fabricated fuel radius profiles have been used to improve the accuracy in evaluating radial swelling strain. The simulations of the IFR experiment X447, which was focused on FCCI/CCCI induced cladding degradation and failure, have been converted into a BISON metallic fuel assessment case with detailed documentation, which is also the first FIPD-BISON integration powered assessment case in BISON. A FIPD-BISON integration repository has been established to be an optional submodule of BISON to support FIPD-BISON integration powered assessment cases. Preliminary trials to simulate the in-cell out-of-pile transient experiment have been made after implementing some essential models such as cladding plasticity models and liquid phase penetration models. The initial results are promising and help identify gaps that need to be filled in FY2023 and beyond, to eventually achieve OPTD-BISON integration powered assessment cases.
Increasing renewable energy requires improving the electricity grid flexibility. Existing measures include power plant cycling and grid-level energy storage, but they incur high operational and investment costs. Using a systems modeling and optimization framework, we study the integration of electrochemical energy storage with individual power plants at various renewable penetration levels. Furthermore, our techno-economic analysis includes both Li-ion and NaS batteries to encompass different technology maturity levels. A California case-study indicates localized integration to be cost-effective for greater grid flexibility. Li-ion batteries can mitigate the residual demand fluctuations of small to medium-sized plants, while NaS batteries would be best-suited for larger storage with higher renewable penetration. Overall, the battery-enabled renewable integration could reduce the unmet grid demand by 75%, the renewable curtailment by 58%, and the CO 2 emission intensity by 16% while including the life cycle emissions of the battery and the renewable farm. Our scenario-based analysis also indicates that rather than replacing all fossil power plants, it is more economical to combine batteries and renewables with individual fossil plants to achieve a clean energy grid.
Nanophotonic resonators are central to numerous applications, from efficient spin–photon interfaces to laser oscillators and precision sensing. A leading approach consists of photonic crystal (PhC) cavities, which have been realized in a wide range of dielectric materials. However, translating proof-of-concept devices into a functional system entails a number of additional challenges, inspiring new approaches that combine resonators with wavelength-scale confinement and high quality factors; scalable integration with integrated circuits and photonic circuits; electrical or mechanical cavity tuning; and, in many cases, a need for heterogeneous integration with functional materials such as III–V semiconductors or diamond color centers for spin–photon interfaces. Here we introduce a concept that generates a finely tunable PhC cavity at a selected wavelength between two heterogeneous optical materials whose properties satisfy the above requirements. The cavity is formed by stamping a hard-to-process material with simple waveguide geometries on top of an easy-to-process material consisting of dielectric grating mirrors and active tuning capability. We simulate our concept for the particularly challenging design problem of multiplexed quantum repeaters based on arrays of cavity-coupled diamond color centers, achieving theoretically calculated unloaded quality factors of 106, mode volumes as small as 1.2(λ/neff)3, and maintaining >60% total on-chip collection efficiency of fluorescent photons. We further introduce a method of low-power piezoelectric tuning of these hybrid diamond cavities, simulating optical resonance shifts up to ∼760 GHz and color center fluorescence tuning of 5 GHz independent of cavity tuning. These results will motivate integrated photonic cavities toward larger scale systems-compatible designs.
Integrative modeling computes a model based on varied types of input information, be it from experiments or prior models. Often, a type of input information will be best handled by a specific modeling software package. In such a case, we desire to integrate our integrative modeling software package, Integrative Modeling Platform (IMP), with software specialized to the computational demands of the modeling problem at hand. After several attempts, however, we have concluded that even in collaboration with the software's developers, integration is either impractical or impossible. The reasons for the intractability of integration include software incompatibilities, differing modeling logic, the costs of collaboration, and academic incentives. In the integrative modeling software ecosystem, several large modeling packages exist with often redundant tools. So we reason, therefore, that the other development groups have similarly concluded that the benefit of integration does not justify the cost. As a result, modelers are often restricted to the set of tools within a single software package. The inability to integrate tools from distinct software negatively impacts the quality of the models and the efficiency of the modeling. As the complexity of modeling problems grows, we seek to galvanize developers and modelers to consider the long-term benefit that software interoperability yields. In this article, we formulate a demonstrative set of software standards for implementing a model search using tools from independent software packages and discuss our efforts to integrate IMP and the crystallography suite Phenix within the Bayesian modeling framework.
To mitigate the variation in demand on the electric grid, thermal energy storage (TES) is an alternative to electric batteries or installing new peaking power plants. Stakeholders and policy makers across the United States have expressed interests in promoting TES, as demonstrated by the US Department of Energy’s Grid-Interactive Efficient Buildings program and the efforts of various state legislatures. However, the cost value provided by TES are unclear. If reliable cost benefits were determined, stakeholders would have a clearer picture of the financial returns that can be gained from their investment in TES. The study in this report is conducted by ORNL with collaboration with Emerson the Helix Innovation Center. In the first part of this report, EnergyPlus was used to perform whole-building simulations for two residential buildings in Indianapolis and Atlanta. The HVAC system in both buildings were equipped with phase change material TES. The TES tank was charged in off-peak hours and discharged in peak hours to perform load shifting. First, the economic value implied by existing time-of-use (TOU) rates offered by utility companies was analyzed via whole-building simulation. Second, existing demand reduction (DR) incentives sourced from 3 different electrical grid administrators (i.e., California, Texas, and New England region) were surveyed to determine their implied value. The study suggests that the traditional value analysis that focuses on ROI for the building owner significantly undervalues TES technology making economic viability difficult. A more comprehensive value analysis that includes peak demand management and deferred capital for peaking power plants shows that TES should be economically viable but here the value is greater for the utility and requires large market penetration and aggregation to fully realize the benefits. Therefore, to facilitate commercialization, new business models are needed that include a broader range of stakeholders and distribute the value of TES proportionally. In the 2 nd part of this project, the benefits of a novel phase change material (PCM) integrated heat pump configuration were evaluated via detailed component based simulation. A one-dimensional PCM heat exchanger model which discretizes the PCM tank and refrigerant tubes into small control volumes is developed. Each control volume can have different PCM temperatures, PCM properties, and heat transfer coefficients. The PCM tank is charged by a wrapped tank condenser and discharged by an internal refrigerant coil. The PCM heat exchanger model is integrated into DOE/ORNL Heat Pump Design Model for heat pump system simulation. To demonstrate the performance of the PCM integrated heat pump, a case study in Chicago was performed. A Time-of-Use utility structure-based control strategy is implemented to schedule the PCM tank charging and discharging mode switching. Compared with a conventional electric heat pump, the PCM integrated heat pump shows superior performance on load shifting and utility cost reduction. As a result, the proposed system demonstrates 24.6% utility saving for cooling application and 25.8% utility saving for heating application.
To mitigate the variation in demand on the electric grid, thermal energy storage (TES) is an alternative to electric batteries and constructing new peaking power plants. The aggressive carbon reduction goals of the United States are accelerating the electrification of building equipment. If space heating and cooling in the US were electrified using conventional heat pump technologies, the winter electric grid peak would approximately double. Integrating TES with heat pump (HP) enables electrification of space cooling and heating devices without overtaxing the grid.This paper evaluates the benefits of a novel phase change material (PCM) integrated heat pump configuration via simulation. A one-dimensional PCM heat exchanger model which discretizes the PCM tank and refrigerant tubes into small control volumes is developed. Each control volume can have different PCM temperatures, PCM properties, and heat transfer coefficients. The PCM tank is charged by a wrapped tank condenser and discharged by an internal refrigerant coil. The PCM heat exchanger model is integrated into DOE/ORNL Heat Pump Design Model for heat pump system simulation.To demonstrate the performance of the PCM integrated heat pump, a case study in Chicago was performed. A Time-of-Use utility structure-based control strategy is implemented to schedule the PCM tank charging and discharging mode switching. Compared with a conventional electric heat pump, the PCM integrated heat pump shows superior performance on load shifting and utility cost reduction. As a result, the proposed system demonstrates 24.6% utility saving for cooling application and 25.8% utility saving for heating application.
The Resource Forecast and Ramp Visualization for Situational Awareness (RAVIS) is an open-source tool for visualizing variable renewable resource forecasts and ramp alerts for significant up/down ramps in renewable resource and the consequent net-load. The modular dashboard of RAVIS contains configurable panes for viewing- probabilistic time series forecasts, ramp event alerts on the look-ahead timeline, spatially resolved resource sites and forecasts, and system simulation and market clearing data such as transmission lines utilization, nodal prices and available generation flexibility. RAVIS uses a technology suite that is assembled to provide optimum visualization facility while maintaining a wide pool of potential deployment and client environments. The tool is designed to take advantage of web application technologies, open source visualization libraries and tooling. Utilizing this technology will enable deployment in any environment, using any operating system, and is scalable to much higher spatial and temporal scales of visualization. As a prototype of the tool and demonstrating a use case of variable renewable integration, RAVIS currently integrates site-specific solar power forecasts in the California Independent System Operator (CAISO) and Mid-continent ISO (MISO) footprint from the IBM WattSun forecasting platform, and also superimposes market simulation data for CAISO footprint from as in-house NREL market clearing tool. The tool has the ability to alert the viewer for excessive up or down ramps for both individual solar sites as well as regionally aggregated net-load ramps, and alerts can also be qualified with respect to available flexible generation. This report will provide an introduction to the RAVIS tool, and summarize the above mentioned capabilities, typical use cases and possible extensions of the tool. The RAVIS development team believes there are likely to be high economic and reliability benefits of integrating probabilistic forecasts of variable renewables into control center visualizations and improved ramp events situational awareness for system operators and forecasting teams in the ISOs and electric utilities in comparison with their business as usual practices.
Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.
The main contribution of this work is to increase the coding productivity of GPU programming by using the concept of Static Graphs. GPU capabilities have been increasing significantly in terms of performance and memory capacity. However, there are still some problems in terms of scalability and limitations to the amount of work that a GPU can perform at a time. To minimize the overhead associated with the launch of GPU kernels, as well as to maximize the use of GPU capacity, we have combined the new CUDA Graph API with the CUDA programming model (including CUDA math libraries) and the OpenACC programming model. We use as test cases two different, well-known and widely used problems in HPC and AI: the Conjugate Gradient method and the Particle Swarm Optimization. In the first test case (Conjugate Gradient) we focus on the integration of Static Graphs with CUDA. In this case, we are able to significantly outperform the NVIDIA reference code, reaching an acceleration of up to 11x thanks to a better implementation, which can benefit from the new CUDA Graph capabilities. In the second test case (Particle Swarm Optimization), we complement the OpenACC functionality with the use of CUDA Graph, achieving again accelerations of up to one order of magnitude, with average speedups ranging from 2x to 4x, and performance very close to a reference and optimized CUDA code. Our main target is to achieve a higher coding productivity model for GPU programming by using Static Graphs, which provides, in a very transparent way, a better exploitation of the GPU capacity. The combination of using Static Graphs with two of the current most important GPU programming models (CUDA and OpenACC) is able to reduce considerably the execution time w.r.t. the use of CUDA and OpenACC only, achieving accelerations of up to more than one order of magnitude. Finally, we propose an interface to incorporate the concept of Static Graphs into the OpenACC Specifications.
Resiliently designed and constructed integrated gas-electric distribution networks (GEDNs) against natural disasters are crucial to social welfare. In this study, a two-stage robust optimisation-based co-expansion planning model is proposed to attain an integrated GEDN with a given resilience level, by optimising the investment strategies of hardening and selective expansion of power distribution feeders and natural gas pipelines, as well as the location and capacity of natural-gas-fired distributed generation. In the first stage, the overall annual investment and operation cost is minimised under normal operation conditions while in the second stage, the feasibility of the investment decisions under the identified worst-case natural disaster scenario is checked with an adjustable load shedding cost criterion. The proposed model is formulated as a mixed integer second-order cone programming problem with the column and constraint generation algorithm employed to seek the optimal solution. Case studies on two integrated GEDNs demonstrate the performance of the proposed methodology.
In recent years, there has been considerable effort to modernize existing and new nuclear power plants with digital instrumentation and control systems (DI&C). However, there has also been considerable concern both by industry and regulatory bodies for the risk and consequence analysis of these systems. Of particular concern are digital common cause failures (CCFs) specifically related to software defects. These “misbehaviors” by the software can occur in both the control and monitoring of a system. While many new methods have been proposed to identify potential software failure modes, such as Systems-theoretic Process Analysis (STPA), Hazard and Consequence Analysis for Digital Systems (HAZCADS), etc., these methods are focused primarily on the control action pathway of a system. In contrast, the information feedback pathway lacks unsafe control actions (UCAs), which are typically related to software basic events; thus, assessment of software basic events in such systems is unclear. In this work, we present the idea of intermediate processors and unsafe information flow (UIF) to help safety analysts trace failure mechanisms in the feedback pathway and how they can be integrated into a fault tree for improved assessment capability. The concepts presented are demonstrated in two comprehensive case studies, a smart sensor integrated platform for unmanned autonomous vehicles and another on a representative advanced human system interface (HSI) for safety critical plant monitoring. The qualitative software basic events are identified, and a fault tree analysis is conducted based on a modified Redundancy-guided Systems-theoretic Hazard Analysis (RESHA) methodology. The case studies demonstrate the use of UIF and intermediate processors in the fault tree to improve traceability of software failures in highly complex digital instrumentation feedback. The improved method can also clarify fault tree construction when multiple component dependencies are present in the system.
In recent years, there has been considerable effort to modernize existing and new nuclear power plants with digital instrumentation and control systems (DI&C). However, there has also been considerable concern both by industry and regulatory bodies for the risk and consequence analysis of these systems. Of particular concern are digital common cause failures (CCFs) specifically related to software defects. These “misbehaviors” by the software can occur in both the control and monitoring of a system. While many new methods have been proposed to identify potential software failure modes, such as Systems-theoretic Process Analysis (STPA), Hazard and Consequence Analysis for Digital Systems (HAZCADS), etc., these methods are focused primarily on the control action pathway of a system. In contrast, the information feedback pathway lacks unsafe control actions (UCAs), which are typically related to software basic events; thus, assessment of software basic events in such systems is unclear. In this work, we present the idea of intermediate processors and unsafe information flow (UIF) to help safety analysts trace failure mechanisms in the feedback pathway and how they can be integrated into a fault tree for improved assessment capability. The concepts presented are demonstrated in two comprehensive case studies, a smart sensor integrated platform for unmanned autonomous vehicles and another on a representative advanced human system interface (HSI) for safety critical plant monitoring. The qualitative software basic events are identified, and a fault tree analysis is conducted based on a modified Redundancy-guided Systems-theoretic Hazard Analysis (RESHA) methodology. The case studies demonstrate the use of UIF and intermediate processors in the fault tree to improve traceability of software failures in highly complex digital instrumentation feedback. The improved method can also clarify fault tree construction when multiple component dependencies are present in the system.
The adopted level densities (LD) for the nuclei produced through different reaction mechanisms significantly impact the accurate calculation of cross sections for the different reaction channels. Many common LD models make simplified assumptions regarding the overall behavior of the total LD and the intrinsic spin and parity distributions of the excited states. However, very few experimental constraints are taken into account in these models: LD at neutron separation energy coming from average spacings of s- and p-wave resonances ( D 0 and D 1 , respectively) whenever they have been previously measured, and the sometimes subjective extrapolation of discrete levels. These, however, constrain the LD only in very specific regions of excitation energy, and for specific spins and parities. This work aims to establish additional experimental constraints on LD through quantitative correlations between cross sections and LD. This allows for fitting and the determination of detailed structures in LD. For this we use the microscopic Hartree-Fock-Bogoliubov (HFB) LD as a starting point as the HFB LD provide a more realistic spin and parity distributions than phenomenological models such as Gilbert-Cameron (GC). We then associate variations predicted by the HFB model with the structure observed in double-differential cross sections at low outgoing neutron energy, a region that is dominated by the LD input. We also use (n, p) on 56 Fe , as an example case where angle-integrated cross sections are extremely sensitive to LD. For comparison purposes we also perform calculations with the GC model. With this approach we are able to perform fits of the LD based on actual experimental data, constraining the model and ensuring its consistency. This approach can be particularly useful in extrapolating the LD to nuclei for which high-excited discrete levels and/or values of D 0 or D 1 are unknown. Finally, it also predicts neutron-induced inelastic γ cross sections that in some cases can differ significantly from more standard phenomenological LD models such as GC.
This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse vertical seismic profiling (VSP) integration. The selected study site is a region around the injector well located within Farnsworth field unit (FWU), where there is an ongoing CO 2 -enhanced oil recovery (EOR) operation. In our study, a site-specific rock physics model was created from extensive geological, geophysical, and geomechanical characterization through 3D seismic data, well logs, and core assessed as part of the 1D MEM conducted on the characterization well within the study area. Here, the Biot-Gassmann workflow was utilized to combine the rock physics and reservoir simulation outputs to determine the seismic velocity change due to fluid substitution. Modeled seismic velocities attributed to mean effective stress were determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from ultrasonic seismic velocity measurements. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer (PSO) was utilized to minimize a penalty function created between the modeled seismic velocities and the observed time-lapse VSP dataset. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.
The increase in the spectral width of an initially monodisperse population of cloud droplets in homogeneous isotropic turbulence is investigated by applying a finite-difference fluid flow model combined with either Eulerian bin microphysics or a Lagrangian particle-based scheme. The turbulence is forced applying a variant of the so-called linear forcing method that maintains the mean turbulent kinetic energy (TKE) and the TKE partitioning between velocity components. The latter is important for maintaining the quasi-steady forcing of the supersaturation fluctuations that drive the increase in the spectral width. We apply a large computational domain (64 3 m 3 ), one of the domains considered in Thomas et al. (2020). The simulations apply 1 m grid length and are in the spirit of the implicit large eddy simulation (ILES), that is, with small-scale dissipation provided by the model numerics. This is in contrast to the scaled-up direct numerical simulation (DNS) applied in Thomas et al. (2020). Two TKE intensities and three different droplet concentrations are considered. Analytic solutions derived in Sardina et al. (2015), valid for the case when the turbulence integral timescale is much larger than the droplet phase relaxation timescale, are used to guide the comparison between the two microphysics simulation techniques. The Lagrangian approach reproduces the scalings relatively well. Representing the spectral width increase in time is more challenging for the bin microphysics because appropriately high resolution in the bin space is needed. The bin width of 0.5 µm is only sufficient for the lowest droplet concentration (26 cm -3 ). For the highest droplet concentration (650 cm -3 ), an order of magnitude smaller bin size is barely sufficient. The scalings are not expected to be valid for the lowest droplet concentration and the high-TKE case, and the two microphysics schemes represent similar departures. Finally, because the fluid flow is the same for all simulations featuring either low or high TKE, one can compare point-by-point simulation results. Such a comparison shows very close temperature and water vapor point-by-point values across the computational domain and larger differences between simulated mean droplet radii and spectral width. The latter are explained by fundamental differences in the two simulation methodologies, numerical diffusion in the Eulerian bin approach and a relatively small number of Lagrangian particles that are used in the particle-based microphysics.
Creep crack growth is a phenomenon which arises in damaged metallic structures under combined primary and secondary loads in the creep regime. The High Temperature Flaw Evaluation Code Committee of the American Society of Mechanical Engineers Boiler and Pressure Vessel Code (BPVC) is evaluating methods in extending Code Case N-934 to capture transient creep crack growth. Here, this paper provides context and basic examples on the selected approach to transient creep crack growth methods based on the C(t)-integral. The basis for selected C(t)-integral solution as analytical method is established, including an overview of its derivation. Practical assessments of a crack growing under creep conditions in a realistic component are conducted to illustrate the analytical approach. Additional considerations in the application of the analytical methods and limitations are discussed.