Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “M CODES”

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 685 records · Page 38

Analysis of Radioactive Waste and Contamination using 3D Position-Sensitive CdZnTe Detectors

INTRODUCTION The US has surplus plutonium. This material will be diluted in gloveboxes at SRS [Ref 1 and 2]. This process will create holdup. Holdup is material in an unexpected or unwanted location. Because this holdup material is primarily PuO2, it must be accounted for (measured). This measurement is typically done using a liquid-nitrogen cooled, HPGe detector. To enable more rapid processing, advanced holdup technologies are being explored [Ref 3 and 4]. METHODS For the surplus plutonium dilute and dispose program, the traditional detector system could slow down the optimal operational tempo of the program. A multiple-detector system, providing complete glovebox coverage, based on CdZnTe crystals, 3D-positioning electronics, and a coded aperture mask, is in the process of being developed by ORNL and SRNL [Ref 4].

Whiteside, Tad [Savannah River National Laboratory↗

An independent analysis of bias sources and variability in wind plant pre-construction energy yield estimation methods

The wind resource assessment community has long had the goal of reducing the bias between wind plant pre-construction energy yield assessment (EYA) and the observed annual energy production (AEP). This comparison is typically made between the 50% probability of exceedance (P50) value of the EYA and the long-term corrected operational AEP (hereafter OA P50), and is known as the P50 bias. The industry has critically lacked an independent analysis of bias reduction investigated across multiple consultants to identify the greatest sources of uncertainty and variance in the EYA process and the best opportunities for uncertainty reduction. The present study addresses this gap by benchmarking consultant methodologies against each other and against operational data at a scale not seen before in industry collaborations. We consider data from 10 wind plants and evaluate discrepancies between eight consultancies in the steps taken from estimates of gross to net energy. Consultants tend to overestimate the gross energy produced at the turbines and then compensate by further overestimating downstream losses, leading to a mean P50 bias near zero, still with significant variability among the individual wind plants. Within our data sample, we find that consultant estimates of all loss categories, except environmental losses, tend to reduce the project-to-project variability of the P50 bias. The disagreement between consultants, however, remains flat throughout the addition of losses. Finally, we find that differences in consultants’ estimates of project performance can lead to differences up to $10/MWh in the levelized cost of energy for a wind plant.

Todd, Austin C.↗

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

ASAS J071404+7004.3 – a close, bright nova-like cataclysmic variable with gusty winds

ABSTRACT Despite being bright (V ≃ 11.8) and nearby (d = 212 pc) ASAS J071404+7004.3 has only recently been identified as a nova-like cataclysmic variable. We present time-resolved optical spectroscopy obtained at the Isaac Newton and the Hiltner and McGraw-Hill Telescopes, together with Swift X-ray and ultraviolet observations. We combined these with TESS photometry and find a period of 3.28 h and a mass transfer rate of $4\!-\!9\times 10^{-9}\, {\mathrm{M_{\odot }\, yr}^{-1}}$. Historical photometry shows at least one low state establishing the system as a VY Scl star. Our high-cadence spectroscopy also revealed rapidly changing winds emanating from the accretion disc. We have modelled these using the Monte Carlo python code and shown that all the emission lines could emanate from the wind – which would explain the lack of double-peaked lines in such systems. In passing, we discuss the effect of variability on the position of cataclysmic variables in the Gaia Hertzsprung–Russell diagram.

Inight, K. (ORCID:0000000222002416)↗

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

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

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

42 ENGINEERING↗

Superelement reduction of substructures for sequential load calculations in OpenFAST

This article presents the superelement formulation newly implemented in OpenFAST to simulate fixed-bottom substructures with a reduced-order model similar to a common industry practice. The Guyan and Craig-Bampton methods are used to reduce the number of degrees of freedom and generate a so-called "superelement". The formulation allows manufacturers to exchange such superelements to perform load calculations without revealing sensitive information about the support structure (e.g., foundation, substructure, and/or tower) or turbine. The source code is made publicly available in the OpenFAST repository. Test cases with varying degrees of complexity are presented to validate the technique, and accuracy issues are discussed.

17 WIND ENERGY↗

Improvements to PVWatts for Fixed and One-Axis Tracking Systems

This work presents improvements to the widely used NREL PVWatts photovoltaic system energy model to improve modeling accuracy for typical fixed and one axis system designs. The aim is to calculate losses in the PV system assuming typical modern system design practices, while maintaining simplicity by keeping the required set of input parameters small. These improvements allow users to more credibly and quickly evaluate competing system designs in early stage feasibility. Common submodels for module cover, spectral, snow, tracker, transformer, plant controller, and self-shading losses, in addition to a bifacial gain option, are incorporated into the PVWatts model, and are shown to improve PVWatts' system performance prediction capabilities without major impact to ease of use. We anticipate including these improvements in a future release of NREL's open source PVWatts code, and some of the features may become available in the System Advisor Model (SAM) desktop software as well as the popular PVWatts web application.

14 SOLAR ENERGY↗

PETSc Users Manual (Rev. 3.13)

This manual describes the use of PETSc for the numerical solution of partial differential equations and related problems on high-performance computers. The Portable, Extensible Toolkit for Scientific Computation (PETSc) is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all message-passing communication. PETSc includes an expanding suite of parallel linear solvers, nonlinear solvers, and time integrators that may be used in application codes written in Fortran, C, C++, and Python. PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or dbx, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than \rolling them" yourself.

97 MATHEMATICS AND COMPUTING↗

Oxidation Behavior and Property Degradation of Nuclear Graphites

During its multidecade operation in the core of nuclear reactors, graphite components are subjected to aggressive and continuous exposure to a high field of ionizing and neutron irradiation, high temperature, and various types of present and postulated chemical attacks. High density, high crystallinity polygranular synthetic graphite is unique among other materials for its extraordinary capacity of resisting and adapting to the aggression inflicted by high temperature, high energy neutron bombardment and ionizing gamma radiation. But, as a carbonaceous material, even though of very high purity, graphite is reactive towards common oxidizing agents: oxygen, carbon dioxide, water. Safe operation of HTGRs relies, among other aspects, on engineered safeguard systems for efficient and continuous protection of graphite components against oxidation. Graphite oxidation behavior was, and continues to be, an important direction of theoretical and experimental research, engineering analyses, models and simulations, and design and safety regulations. The avalanche of publications, reports, experimental data, computer codes, and regulatory documents related to oxidation behavior of nuclear graphite is now accelerating to new levels, prompted by the increased interest for nuclear energy as a clean, carbon-free energy source. Even though public’s perception of nuclear energy advantages may still be influenced by the memories of past accidents of nuclear reactors from generations II and III, the community of informed scientists and engineers, regulators and statemen knows that generation IV of nuclear reactors is designed at very high safety standards, doubled by great advances of scientific knowledge and technological progress. One of routes of these recent advances is directed at better understanding of graphite oxidation behavior, its relationship with graphite manufacturing and microstructural properties, along with the effects of various environmental factors and process variables. Together, the recent progress in manufacturing, properties characterization, and modeling of intricated physical and chemical processes that concur to the oxidation behavior led to development of powerful simulation codes able to analyze various scenarios of normal operation and hypothetical off-normal events, and thus to clearly specify the allowable parameters envelopes for the designers, constructors, and operators of current and future modular HTGRs. This review begins with an introduction on manufacturing methods, structure, and properties of nuclear graphite, including basic requirements that this specialty graphite type must satisfy for nuclear use. It continues with a chapter on environmental effects on nuclear graphite, where emphasis is placed less on irradiation and much more on oxidation phenomena, their safety implications, and the basic traits of chronic and acute oxidation by air (oxygen) and water (humidity, steam). Particular attention is placed on the three graphite grades of interest for this document (IG-110, NBG-18, PCEA). A chapter on properties degradation induced by oxidation follows, with focus on density, dimensional, and mechanical properties changes. The next chapter is intended as a brief review of various approaches used for modeling of graphite oxidation behavior. It summarizes the progress of oxidation models, from the early attempts to complex computational approaches interfaced with specialized computer codes designed for nuclear reactor simulations. Last, a list is presented of knowledge gaps where more research is needed. A short summary concludes the review.

36 MATERIALS SCIENCE↗

The waterSHED Model: User Guide

The ideal design and operation of small hydropower plants is a complex optimization problem with economic, social, and environmental objectives. The waterSHED (Water Allocation Tool Enabling Rapid Small Hydropower Environmental Design) model is a user-friendly tool that allows hydropower stakeholders to model the trade-offs among these objectives using the Standard Modular Hydropower (SMH) framework. The SMH framework employs modular technologies that can be represented as blackbox objects and combined within a river to create a hydropower facility. For a given site, the waterSHED model aims to determine which modules should be placed in a facility and how those modules should be operated. This user guide describes how to use the graphical user interface and related functionalities. This document also summarizes the background research and mathematical formulations that are explained indepth in the accompanying doctoral dissertation. This model is an early step toward a new hydropower design process that employs standardization and modularity to reduce costs, development timelines, and challenges regarding social and environmental mitigation measures for low-head, small hydropower development. The waterSHED model is a Python application that will require the ability to download a GitHub repository, import the necessary packages, and run a set of Python script files using an integrated development environment. The script produces a graphical user interface to coordinate inputs, simulate operation, and visualize results, so no coding experience is needed once the script is running. Additionally, the waterSHED Workbook is a Microsoft Excel file that works with the Python script to facilitate data entry

13 HYDRO ENERGY↗

Modernization efforts for the R-Matrix code SAMMY [Slides]

This presentation details modernization efforts for the R-Matrix code SAMMY. It provides an overview of SAMMY, modernization and maintenance goals, modernization accomplishments, covariance data, SAMMY API development, as well as future plans.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NUM-DAT File Format Specification: Used in M-9 Gun Experiment Data Archiving

The M-9 Shock and Detonation Physics group executes experiments on gun and explosive platforms with large numbers of oscilloscopes used for data acquisition. The data acquisition from these oscilloscopes was automated many years ago using a custom piece of software called RunDig . The default save format from this software is a custom structure referred to as "NUM-DAT" format. This file format includes a text ".DAT" file which is a header file used to interpret the binary ".NUM" file which contains the oscilloscope data. The data save format was originally developed by John Vorthman and has been in use by M-9 personnel for over 20 years. This data format has been used for archiving data from experiments performed by M-9 personnel at TA-40, TA-39, and the TA-55 Impact Test Facility. Numerous custom analysis and visualization programs have also been developed, and continue to be used, that utilize this data format. This document describes the NUM-DAT format and provides code examples for reading the format and converting it to other formats.

47 OTHER INSTRUMENTATION↗

Geographically aware estimates of remotely sensed water properties for Chesapeake Bay

This software contains the code for a data-driven model to produce high-resolution (< 1 km2) estimates of temperature, salinity, and turbidity over decadal time scales as required by next-generation Earth Systems Models. Our model fits daily MODIS Aqua reflectance data to surface observations (< 1 m depth) from 2000-2021 in Chesapeake Bay, USA.

Stachelek, Jemma↗

Quasi-Coherent Fluctuations Limiting the Pedestal Growth on Alcator C-Mod: Experiment and Modeling

Performance predictions for future fusion devices rely on an accurate model of the pedestal structure. The leading candidate for predictive pedestal structure is EPED, and it is imperative to test the underlying hypotheses to further gain confidence for ITER projections. Here, we present experimental work testing one of the EPED hypotheses, namely the existence of a soft limit set by microinstabilities such as the kinetic ballooning mode (KBM). This work extends recent work on Alcator C-Mod [Diallo, et al., Phys. Rev. Lett., , (2014), 115001], to include detailed measurements of the edge fluctuations and comparisons of edge simulation codes and experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

1.2.2.405 HydroWIRES Topic D1: Capacity Expansion Model (CEM) Enhancements

Long-term grid planning tools have difficulty representing detailed hydropower operating characteristics, which depend not only on technological specifications but also on water management practices and regulations. As a result, the value of hydropower is incompletely characterized, and the potential role of hydropower in the performance and resiliency of the future electric grid is not fully understood. This work will fill that gap by developing new ways to represent hydropower resource, technology, and operational characteristics in electric sector capacity expansion models and implementing them in the open-source version of the National Renewable Energy Laboratory's Regional Energy Deployment System (ReEDS) model. ReEDS is a well-established national scale grid planning tool used since 2003 by the U.S. Department of Energy and others to explore the evolution of the U.S. electric sector. Improvements will include a comprehensive national resource assessment for pumped storage hydropower and methods for modeling multiple hydropower technology categories characterized by technical, regulatory, and economic characteristics. The project will provide guiding principles and strategies for improving hydropower modeling in capacity expansion models and deliver a first-of-its kind versatile PSH dataset. All data, code, and methods will be publicly available, allowing the industry to better identify the value of hydropower in the future electricity system and make more informed planning decisions.

capacity expansion↗

Constraints on cosmology and baryonic feedback with joint analysis of Dark Energy Survey Year 3 lensing data and ACT DR6 thermal Sunyaev-Zel'dovich effect observations

We present a joint analysis of weak gravitational lensing (shear) data obtained from the first three years of observations by the Dark Energy Survey and thermal Sunyaev-Zel'dovich (tSZ) effect measurements from a combination of Atacama Cosmology Telescope (ACT) and Planck data. A combined analysis of shear (which traces the projected mass) with the tSZ effect (which traces the projected gas pressure) can jointly probe both the distribution of matter and the thermodynamic state of the gas, accounting for the correlated effects of baryonic feedback on both observables. We detect the shear$~\times~$tSZ cross-correlation at a 21$\sigma$ significance, the highest to date, after minimizing the bias from cosmic infrared background leakage in the tSZ map. By jointly modeling the small-scale shear auto-correlation and the shear$~\times~$tSZ cross-correlation, we obtain $S_8 = 0.811^{+0.015}_{-0.012}$ and $\Omega_{\rm m} = 0.263^{+0.023}_{-0.030}$, results consistent with primary CMB analyses from Planck and P-ACT. We find evidence for reduced thermal gas pressure in dark matter halos with masses $M < 10^{14} \, M_{\odot}/h$, supporting predictions of enhanced feedback from active galactic nuclei on gas thermodynamics. A comparison of the inferred matter power suppression reveals a $2-4\sigma$ tension with hydrodynamical simulations that implement mild baryonic feedback, as our constraints prefer a stronger suppression. Finally, we investigate biases from cosmic infrared background leakage in the tSZ-shear cross-correlation measurements, employing mitigation techniques to ensure a robust inference. Our code is publicly available on GitHub.

Pandey, S. [Johns Hopkins U.; Columbia U.] (ORCID:↗