Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code of Accounts”

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 307 records · Page 17

Inverse bremsstrahlung absorption rate for super-Gaussian electron distribution functions including plasma screening

Here we provide analytic expressions for the effective Coulomb logarithm for inverse bremsstrahlung absorption which predict significant corrections to the Langdon effect and overall absorption rate compared to previous estimates. The calculation of the collisional absorption rate of laser energy in a plasma by the inverse bremsstrahlung mechanism usually makes the approximation of a constant Coulomb logarithm. We dispense with this approximation and instead take into account the velocity dependence of the Coulomb logarithm, leading to a more accurate expression for the absorption rate valid in both classical and quantum conditions. In contrast to previous work, the laser intensity enters into the Coulomb logarithm. In most laser-plasma interactions the electron distribution function is super-Gaussian [Langdon, Phys. Rev. Lett. 44, 575 (1980)], and we find the absorption rate under these conditions is increased by as much as ≈ 30% compared to previous estimates at low density. In many cases of interest the correction to Langdon's predicted reduction in absorption is large; for example at Z = 6 and Te = 400 eV the Langdon prediction for the absorption is in error by a factor of ≈ 2. However, we also account for the additional effect of plasma screening, which predicts a reduction in absorption by a similar amount (up to ≈ 30%). These two effects compete to determine the overall absorption, which may be increased or decreased, depending on the conditions. The corrections can be incorporated into radiation-hydrodynamics simulation codes by replacing the familiar Coulomb logarithm with an analytic expression which depends on the super-Gaussian order “M” and the screening length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FY2022 Status Update: A Probabilistic Model for Stress Corrosion Cracking of SNF Dry Storage Canisters

Understanding the potential risk of stress corrosion cracking of spent nuclear fuel dry storage canisters has been identified as a knowledge gap for determining the safety of long-term interim storage of spent nuclear fuel. To address this, the DOE is funding a multi-lab DOE effort to understand the timing, occurrence, and consequences of potential canister SCC. Sandia National Laboratories has developed a probabilistic model for canister penetration by SCC. This model has been continuously updated at SNL since 2014. Model uncertainties are treated using a nested loop structure, where the outer loop accounts for uncertainties due to lack of data and the inner aleatoric loop accounts for uncertainties due to variation in nature. By separating uncertainties into these categories, it is possible to focus future work on reducing the most influential epistemic uncertainties. Several experimental studies have already been performed to improve the modeling approach through expanded process understanding and improved model parameterization. The resulting code is physics-based and intended to inform future work by identifying (1) important modeling assumptions, (2) experimental data needs, and (3) necessary model developments. In this document, several of the sub-models in the probabilistic SCC model have been exercised, and the intermediate results, as the model progresses from one sub-model to the next, are presented. Evaluating the sub-models in this manner provides a better understanding of sub-model outputs and has identified several unintended consequences of model assumptions or parameterizations, requiring updates to the modeling approach. The following updates have been made, and future updates have been identified.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A symbolic framework to obtain mid-fidelity models of flexible multibody systems with application to horizontal-axis wind turbines

Abstract. The article presents a symbolic framework (also called computer algebra program) that is used to obtain, in symbolic mathematical form, the linear and nonlinear equations of motion of a mid-fidelity multibody system including rigid and flexible bodies. Our approach is based on Kane's method and a nonlinear shape function representation for flexible bodies. The shape function approach does not represent the state of the art for flexible multibody dynamics but is an effective trade-off to obtain mid-fidelity models with few degrees of freedom, taking advantage of the separation of space and time. The method yields compact symbolic equations of motion with implicit account of the constraints. The general and automatic framework facilitates the creation and manipulation of models with various levels of complexity by adding or removing degrees of freedom. The symbolic treatment allows for analytical gradients and linearized equations of motion. The linear and nonlinear equations can be exported to Python code or dedicated software. There are multiple applications, such as time domain simulation, stability analyses, frequency domain analyses, advanced controller design, state observers, and digital twins. In this article, we describe the method we used to systematically generate the equations of motion of multibody systems and present the implementation of the framework using the Python package SymPy. We apply the framework to generate illustrative land-based and offshore wind turbine models. We compare our results with OpenFAST simulations and discuss the advantages and limitations of the method. The Python implementation is provided as an open-source project.

Branlard, Emmanuel (ORCID:0000000277506128)↗

Long-run Marginal Emission Rates for Electricity - Workbooks for 2023 Cambium Data

These workbooks contain modeled estimates of long-run marginal emission rates (LRMER) for the contiguous United States' electric sector. A LRMER is an estimate of the rate of emissions that would be either induced or avoided by a change in electric demand, taking into account how the change could influence both the operation as well as the structure of the grid (i.e., the building and retiring of capital assets, such as generators and transmission lines). These workbooks provide data for 18 GEA regions covering the contiguous United States. Mappings of these regions to ZIP codes and counties is given in this workbook in the corresponding tabs. For more data underlying these emissions factors, see the Cambium 2023 project at https://scenarioviewer.nrel.gov/. For more details on input assumptions and methodology see the associated report (Cambium 2023 Scenario Descriptions and Documentation, https://www.nrel.gov/docs/fy24osti/88507.pdf). Users are advised to review section 4 of the report, which discusses limitations and caveats of the data. This data is planned to be updated annually. Information on the latest versions can be found at https://www.nrel.gov/analysis/cambium.html.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Novel Workflow for Coupled Simulation of Hydraulic Stimulation with Simultaneous Injection of Proppant

Enhanced Geothermal Systems (EGS) are located in impermeable host rock. Therefore, to achieve promising circulation, hydraulic stimulation is a key factor. Stimulation involves injecting a combination of fluids and proppant particles into existing fractures under high pressure. The effectiveness of stimulation depends on a variety of different factors, including the pressurization scheme and the optimal choice of carrying fluid and proppant particles to be used. For example, high-density particles can settle out of the injected slurry quickly, whereas lighter particles as well as a dense carrying fluid allow for a long transport into the reservoir. A numerical tool is therefore essential to take all these factors into account and to make robust predictions about the performance of the hydraulic stimulation. We present a novel computational application for coupled numerical simulation of hydraulic stimulation with concurrent injection of proppant in a fractured geothermal reservoir. This application is based on the open-source MOOSE framework that has been enriched with a set of equations to account for the circulation of proppant-laden slurries. It accounts for both the fluid and proppant components of the mixture, as well as particle-driven processes that include gravity settling, particle-particle interactions, and strong density and viscosity contrasts. A particular challenge is the mutual coupling with fracture mechanical processes in a propagating fracture using the FE code Defmod. Aperture and fluid pressure on one side and slurry density and viscosity changes on the other side interact and affect each other. Their basic understanding is therefore essential for the success of the stimulation. We compare our application to common numerical benchmark problems for both slurry flow and hydraulic stimulation.

15 GEOTHERMAL ENERGY↗

Development and application of sorption mass transfer models for fission product transport in TRISO fuel systems using BISON

Fission product mass transfer within and between TRISO particle and compact layers is an important phenomenon. It directly impacts fission product release predictions, which are used as source terms for safety and licensing calculations. Modeling fission product transport within layers and across bonded interfaces is relatively straightforward under the assumptions of isotropic Fickian diffusion, concentration continuity, and flux continuity. Modeling fission product transport across gaps and debonded layers is more difficult. Gaps are known to form between the buffer and inner PyC (IPyC), and debonding may occur between the IPyC and silicon carbide (SiC). A new mass transfer model was developed to provide a more accurate and robust fission product release calculation by accounting for interlayer sorptivity. One importance of the new model is to account for temperature and material changes across the gap based on sorption isotherm. This work presents the development of the sorption behavioral model and a fission product trapping model, and analyses the subsequent fission product diffusion behavior, specifically cesium, through TRISO particles and compact materials using the BISON fuel performance code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Area Localized Coupled Model for Analytical Mean Flow Prediction in Arbitrary Wind Farm Geometries

This work introduces the area localized coupled (ALC) model, which extends the applicability of approaches that couple classical wake superposition models and atmospheric boundary layer models to wind farms with arbitrary layouts. Coupling wake and top–down boundary layer models is particularly challenging since the latter requires averaging over planform areas associated with turbine-specific regions of the flow that need to be specified. The ALC model uses Voronoi tessellation to define this local area around each turbine. A top–down description of a developing internal boundary layer is then applied over Voronoi cells upstream of each turbine to estimate the local mean velocity profile. Coupling between the velocity at hub-height based on this localized top–down model and a wake model is achieved by enforcing a minimum least-square-error in mean velocity in each cell. The wake model in the present implementation takes into account variations in wind farm inflow velocity and represents the wake profile behind each turbine as a super-Gaussian function that smoothly transitions between a top-hat shape in the region immediately following the turbine to a Gaussian profile downstream. Detailed comparisons to large-eddy simulation (LES) data from two different wind farms demonstrate the efficacy of the model in accurately predicting both wind farm power output and local turbine hub-height velocity for different wind farm geometries. These validations using data generated from two different LES codes demonstrate the model's versatility with respect to capturing results from different simulation setups and wind farm configurations.

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

Validation Assessment of Turbulent Reacting Flow Model Using The Area-Validation Metric on Medium-Scale Methanol Pool Fire Results

Accident analysis and ensuring power plant safety are pivotal in the nuclear energy sector. Significant strides have been achieved over the past few decades regarding fire protection and safety, primarily centered on design and regulatory compliance. Yet, after the Fukushima accident a decade ago, the imperative to enhance measures against fire, internal flooding, and power loss has intensified. Hence, a comprehensive, multilayered protection strategy against severe accidents is needed. Consequently, gaining a deeper insight into pool fires and their behavior through extensive validated data can greatly aid in improving these measures using advanced validation techniques. A model validation study was performed at Sandia National Laboratories in which a 30-cm diameter methanol pool fire was modeled using the SIERRA/Fuego turbulent reacting flow code. This validation study used a standard validation experiment to compare model results against, and conclusions have been published. The fire was modeled with a Large Eddy Simulation (LES) turbulence model with subgrid turbulent kinetic energy closure. Combustion was modeled using a strained laminar flamelet library approach. Radiative heat transfer was accounted for with a model utilizing the gray-gas approximation. In the present study, additional validation analysis is performed using the area validation metric (AVM). These activities are done on multiple datasets involving different variables and temporal/spatial ranges and intervals. In conclusion, the results provide insight into the use of the area validation metric on such temporally varying datasets and the importance of physics-aware use of the metric for proper analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Assessment of Microreactor Safety Analysis Challenges and Recommendations for Utilization of the Comprehensive Reactor Analysis Bundle

To enable the broad deployment of microreactors in fundamentally new application regimes (i.e., mobile and autonomous operations), their safety must be indisputable in terms of possessing inherent resistance to severe offsite dose consequences. Therefore, mechanistic beyond-design-basis event source term calculations that demonstrate a sufficiently large margin of safety will be required to accommodate these new application regimes, which have no history of commercial regulation. Even for traditional reactor operation configurations, safety analysis expertise and familiarity with accident phenomena and conditions in microreactors—specifically those with heat pipe primary cooling arrangements—are lacking compared with other advanced reactor concepts and small modular reactors. Recently, modeling and simulation tools to account for unique heat pipe design aspects have been developed by Sandia National Laboratories with MELCOR and by the US Department of Energy’s (DOE’s) Office of Nuclear Energy Advanced Modeling and Simulation Program with BlueCRAB. However, further demonstration and assessment of potential knowledge gaps are needed to support these codes’ broad usage by the microreactor community. Through the DOE Microreactor Program, an initial assessment of these two tools and guidance on how an evaluation model could be constructed was performed and is reported herein. Moreover, a proposed approach for demonstrating an evaluation model using these two tools is outlined.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

User Friendly Web-Based Tool to Assess the Energy Efficiency and Durability of Residential Wall Retrofits

According to the U.S. Department of Energy Windows and Building Envelope Research and Development Roadmap for Emerging Technologies, building envelope wall energy loss in the United States accounts for about 5.9 x 1010 kWh or 2 quads of energy annually, costing homeowners and occupants billions of dollars. Enclosure retrofits targeting these losses can save significant energy, reduce greenhouse gas emissions, and save occupants millions of dollars over time. Older homes, built before 1992 when DOE’s Building Energy Codes Program was established, represent approximately 68 percent of the residential building stock in the country, often having significant air leakage and inadequate insulation. Homes with little to no air sealing or insulation have heating and cooling losses that can represent a substantial portion of utility bills.High-performance building envelope retrofit systems are rarely selected for retrofit applications. Current solutions are expensive and/or unfamiliar to many designers, builders, contractors, and code officials and therefore are perceived as risky. The dominant perceived risk is durability specifically related to condensation and moisture accumulation in the building envelope component.The Building Science Advisor (BSA) is a rule-based expert system web-based tool that was originally developed to assist building professionals in designing energy efficient and durable wall systems for new construction. With the present focus being placed on upgrading the existing building stock, a retrofit module has been developed that, based on the location, existing construction, and planned retrofit strategy, provides recommendations on how to address the retrofit in a manner that will perform in accordance with IECC 2018 building code and not create a durability problem. This paper will describe the development of this tool and demonstrate its features and capabilities.

Desjarlais, Andre Omer↗

Modeling and simulation of oxygen transport in high burnup LWR fuel

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. Here, the transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Demonstration of μCHP in Light Commercial Hot Water Applications

Internal combustion engine driven combined heat and power (CHP) systems produce power while the waste heat is recovered and used for another purpose. For µCHP systems (less than 50kWe output) overall efficiencies can reach over 90% for products currently on the market throughout the world. On a fuel basis this offers an advantage over using grid-based power and another fuel for heat. Additionally, the cost of natural gas on an energy unit basis is significantly less expensive than electric power. This difference can lead to substantial savings for the end user of a CHP system. In some parts of the world, particularly Germany and Japan, µCHP has become a common part of the heat and distributed generation equipment mix. Market analysis conducted in 2012 showed that in North America there is potential of a 5,300 units per year market for µCHP, resulting in an annual energy savings of approximately 13.5 TBtu and more than $84 million in annual customer savings. However, µCHP still has not seen much market penetration in North America. This study helps to identify some of the reasons why, address roadblocks that still need to be overcome and provide experience to established best practices. Nine commercial sites were initially identified (fitness club, surgical center, senior community, resort and spa, commercial laundry, two multi-family housing, restaurant and manufacturing facility). Systems from two different manufacturers were installed at five sites total, but due to unforeseen circumstances, only three of those sites were finally commissioned: test data were collected by Oak Ridge National Laboratory. While not all systems were installed as intended, the project resulted in the following key findings: Micro combined heat and power is an unfamiliar technology for building owners and plumbing and electrical contractors that serve the buildings most suited for µCHP. Manufacturers will need to be aggressive in providing training and providing broad support to successfully bring the technology to market.; The proper integration of µCHP into other building systems is the most critical part of a successful installation. Manufacturers will need to be able to provide the expertise to train, advise and support the mechanical contractors during the design and installation of the heat integration system.; Regulations governing µCHP varies across the country. There are a few UL standards (UL 2200 and UL 1741) that apply nationally. Emissions are governed by EPA except in California where the California Air Resources Board (CARB) has established their own standard. Requirements that dictate installation and interconnection are dependent on local building codes and local utility processes and requirements. Some out of date building codes do not offer clear guidance for these requirements, leading to confusion. This makes developing training and determining broadly applicable best practices difficult.; Installation costs of µCHP are a key factor in making an installation successful. It is critical that these costs and all national, state and local regulations be understood and accounted for by all parties involved in the installation before a project is undertaken.; A robust customer service organization is important to meet the needs of the customers. These systems require annual maintenance and service and local capability is needed to keep costs inline.

03 NATURAL GAS↗

Slender-body approach for computing second-order wave loads in the frequency domain

This work presents a slender-body approach to evaluate the second-order wave loads acting on a floating structure in the frequency domain. The approach is in the same spirit as the common use of Morison’s equation to approximate the wave loads without solving the radiation/diffraction problem. To do so, we employ Rainey’s equation, which can be seen as an extension of the inertial part of Morison’s equation to include nonlinear effects. We introduce modifications to Rainey’s formulation in order to evaluate wave kinematics at the mean body position instead of the original approach of considering instantaneous displacements. We also propose a simple approximation to partially account for wave scattering effects on the second-order loads based on the analytical solution of a surface-piercing bottom-mounted vertical circular cylinder. Though limited to structures composed of cylinders, this slender-body approach is orders of magnitude faster than computing second-order wave coefficients with a radiation/diffraction code. We implemented this approach for difference-frequency (slow drift) loads in an open-source frequency-domain floating wind turbine model. We present comparisons against results obtained with radiation/diffraction theory for three reference floating wind turbine designs: the OC3-Hywind spar, the OC4-DeepCwind semisubmersible, and the VolturnUS-S semisubmersible. In general, the results show that the proposed slender-body approach with the correction to approximate wave scattering effects provides useful estimations of the difference-frequency wave loads and the resulting motions of the floater.

17 WIND ENERGY↗

Evaluation of Gamma Ray Production - Quarterly Report FY22Q2

The LLNL team has performed the first systematic study of the effect of two-body weak currents in β decay rates across the entire mass table. In light and medium-mass nuclei, two-body weak currents from chiral effective field theory account for a significant portion of the phenomenological quenching of Gamow-Teller transition matrix elements. We examined the systematic effects of two-body axial currents on Gamow-Teller strength and β-decay rates in heavy nuclei within energy-density functional theory. We extended the charge-changing finiteamplitude pnFAM code to incorporate the contributions of two-body currents to the usual onebody linear response in the Gamow-Teller channel. We found that the two-body currents, as expected, usually quench both summed Gamow-Teller strength and decay rates, but by an amount that decreases as the neutron excess grows; see Fig. 1. In addition, they can enhance individual low-lying transitions, leading to decay rates that are quite different from those that an energy-independent quenching would produce, particularly in neutron-rich nuclei. This new physics capability is unique and provides the highest-fidelity description of β decay on the market. The underlying computational implementation can be reused to calculate γ-strength functions for finite-range nuclear potentials such as the Gogny force, which will allow us to benchmark some of our future developments against published results. The paper is in print in Phys. Rev. C.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ion-temperature- and density-gradient-driven instabilities and turbulence in Wendelstein 7-X close to the stability threshold

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the δƒ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Surrogate modeling of Cellular-Potts agent-based models as a segmentation task using the U-Net neural network architecture

The Cellular-Potts model is a powerful and ubiquitous framework for developing computational models for simulating complex multicellular biological systems. Cellular-Potts models (CPMs) are often computationally expensive due to the explicit modeling of interactions among large numbers of individual model agents and diffusive fields described by partial differential equations (PDEs). In this work, we develop a convolutional neural network (CNN) surrogate model using a U-Net architecture that accounts for periodic boundary conditions. We use this model to accelerate the evaluation of a mechanistic CPM previously used to investigate in vitro vasculogenesis. The surrogate model was trained to predict 100 computational steps ahead (Monte-Carlo steps, MCS), accelerating simulation evaluations by a factor of 562 times compared to single-core CPM code execution on CPU. Over short timescales of up to 3 recursive evaluations, or 300 MCS, our model captures the emergent behaviors demonstrated by the original Cellular-Potts model such as vessel sprouting, extension and anastomosis, and contraction of vascular lacunae. This approach demonstrates the potential for deep learning to serve as a step toward efficient surrogate models for CPM simulations, enabling faster evaluation of computationally expensive CPM simulations of biological processes.

97 MATHEMATICS AND COMPUTING↗

New Conductive Opacities for White Dwarf Envelopes

Thanks to their continuous cooling and relative simplicity, white dwarf stars are routinely used to measure the ages of stellar populations. The usefulness of white dwarfs as cosmochronometers depends on the availability of accurate cooling models. A key ingredient of those models is conductive opacities, which largely govern the cooling rate. In this work, we present improved conductive opacities for the regime of moderate coupling and moderate degeneracy that characterizes an important portion of the envelopes of DA and DB white dwarfs. Here, we find differences of up to a factor of 3 between our calculations and the commonly used opacities of Cassisi et al., which we attribute to an improved account of electron–electron scattering. The cooling models are strongly affected by those changes in the conductive opacities: the age of a 4000 K white dwarf can be reduced by as much as 2 Gyr. We provide analytical fits to our new opacities to facilitate the implementation of this important effect in white dwarf evolution codes.

79 ASTRONOMY AND ASTROPHYSICS↗

The 2021 Annual CSEWG Meeting at Nuclear Data Week(s)

From November 15-19, 2021, CSEWG held its annual meeting during Brookhaven National Laboratory’s Nuclear Data Week(s)1. The combined event was held via zoom (due to COVID travel restrictions) and there were 175 registered attendees. This was CSEWG’s 3rd fully virtual meeting and, by all accounts, it was a success. The timetable and talks are available at https://indico.bnl.gov/event/13121/overview for registered participants. This meeting began with the welcoming of the new ENDF Library Manager, Gustavo Nobre, the announcement that ENDF/B-VIII.1’s release will be delayed until February 2024 due to COVID induced delays, and the adoption of a draft Code of Conduct. This last point is notable since CSEWG is known for its lively physics and engineering discussions. However, as professionals, we must endeavor to keep our discussions civil, collegial, and inclusive.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗