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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Detailed modeling of electron emission for transpiration cooling of hypersonic vehicles

Electron transpiration cooling (ETC) is a recently proposed approach to manage the high heating loads experienced at the sharp leading edges of hypersonic vehicles. Computational fluid dynamics (CFD) can be used to investigate the feasibility of ETC in a hypersonic environment. In this paper, a modeling approach is presented for ETC, which includes developing the boundary conditions for electron emission from the surface, accounting for the space-charge limit effects of the near-wall plasma sheath. The space-charge limit models are assessed using 1D direct-kinetic plasma sheath simulations, taking into account the thermionically emitted electrons from the surface. The simulations agree well with the space-charge limit theory proposed by Takamura et al. for emitted electrons with a finite temperature, especially at low values of wall bias, which validates the use of the theoretical model for the hypersonic CFD code. The CFD code with the analytical sheath models is then used for a test case typical of a leading edge radius in a hypersonic flight environment. The CFD results show that ETC can lower the surface temperature of sharp leading edges of hypersonic vehicles, especially at higher velocities, due to the increase in ionized species enabling higher electron heat extraction from the surface. Finally, the CFD results also show that space-charge limit effects can limit the ETC reduction of surface temperatures, in comparison to thermionic emission assuming no effects of the electric field within the sheath.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

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↗
DOE OSTI2017

Toward a complete theory for predicting inclusive deuteron breakup away from stability

We present an account of the current status of the theoretical treatment of inclusive (d, p) reactions in the breakup-fusion formalism, pointing to some applications and making the connection with current experimental capabilities. Three independent implementations of the reaction formalism have been recently developed, making use of different numerical strategies. The codes also originally relied on two different but equivalent representations, namely the prior (Udagawa-Tamura, UT) and the post (Ichimura-Austern-Vincent, IAV) representations. The different implementations have been benchmarked for the first time, and then applied to the Ca isotopic chain. The neutron-Ca propagator is described in the Dispersive Optical Model (DOM) framework, and the interplay between elastic breakup (EB) and non-elastic breakup (NEB) is studied for three Ca isotopes at two different bombarding energies. The accuracy of the description of different reaction observables is assessed by comparing with experimental data of (d, p) on 40,48Ca. We discuss the predictions of the model for the extreme case of an isotope (60Ca) currently unavailable experimentally, though possibly available in future facilities (nominally within production reach at FRIB). We explore the use of (d, p) reactions as surrogates for (n,γ) processes, by using the formalism to describe the compound nucleus formation in a (d,pγ) reaction as a function of excitation energy, spin, and parity. The subsequent decay is then computed within a Hauser-Feshbach formalism. Comparisons between the (d,pγ) and (n,γ) induced gamma decay spectra are discussed to inform efforts to infer neutron captures from (d,pγ) reactions. Finally, we identify areas of opportunity for future developments, and discuss a possible path toward a predictive reaction theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Input-Output Model for MACCS Nuclear Accident Imp Acts Estimation

Since the original economic model for MACCS was developed, better quality economic data (as well as the tools to gather and process it) and better computational capabilities have become available. The update of the economic impacts component of the MACCS legacy model will provide improved estimates of business disruptions through the use of Input-Output based economic impact estimation. This paper presents an updated MACCS model, bases on Input-Output methodology, in which economic impacts are calculated using the Regional Economic Accounting analysis tool (REAcct) created at Sandia National Laboratories. This new GDP-based model allows quick and consistent estimation of gross domestic product (GDP) losses due to nuclear power plant accidents. This paper outlines the steps taken to combine the REAcct Input-Output-based model with the MACCS code, describes the GDP loss calculation, and discusses the parameters and modeling assumptions necessary for the estimation of long-term effects of nuclear power plant accidents.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

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↗
DOE OSTI2022

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)↗
DOE OSTI2024

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↗
NASA NTRS1982

Speed-dependent collision effects on radar back-scattering from the ionosphere

The question whether the differences between fluctuation spectra for linearly speed-dependent and speed-independent collision frequencies could account for disagreements between rocket and incoherent scatter estimate was addressed. The basic theory used for computing the fluctuation spectrum is outlined. The speed-dependence of the charge-neutral collision frequency is discussed, with special emphasis on its derivation from the mobility measurements. Various developments of the computer code used for the computation of the fluctuation spectrum are described. The range of values of input parameters typical to the collision-dominated ionosphere are briefly described. The computational results are presented, and the significance and limitation of these results and the future scope of the research are discussed.

Behl, Y. K.↗
NASA NTRS1992

LDEF: Dosimetric measurement results (AO 138-7 experiment)

One of the objectives of the AO 138-7 experiment on board the LDEF was a total dose measurement with Thermo Luminescent Detectors (TLD 100). Two identical cases, both including 5 TLDs inside various aluminum shields, are exposed to the space environment in order to obtain the absorbed dose profile induced. Radiation fluence received during the total mission length was computed, taking into account the trapped particles (solar maximum and solar minimum periods) and the cosmic rays; due to the magnetospheric shielding, the solar proton fluences are negligible on the LDEF orbit. The total dose induced by these radiations inside a semi-infinite plane shield of Al are computed with radiation transport codes. TLD reading are performed after flight; due to the mission duration increase, a post-flight calibration was necessary in order to cover the range of the flight induced dose. The results obtained, similar (+ or - 30 pct.) in both cases, are compared with the dose profile computation. In practice, these LDEF results, with less than a factor 1.4 between measurements and forecasts, reinforce the validity of the computation methods and models used for the long term evaluation of space radiation intensity on low inclination Earth orbits.

Bourrieau, J.↗
NASA NTRS1997

Study of Magnetic Motions in the Solar Photosphere and their Implications for Heating the Solar Atmosphere

We continued our program of CO observations with the McMath-Pierce facility at Kitt Peak National Solar Observatory. Uitenbroek has developed a two-and three dimensional radiative transfer code that now includes chemical equilibrium calculations. This code allows us to compute a CO spectrum from for instance a snapshot of a solar granulation simulation (e.g. Stein & Nordlund 1989, Apj 342, L95) and compare these theoretical spectra with our spatially resolved CO spectroscopy. Van Ballegooijen and Uitenbroek have started calculations of two-dimensional fluxtube models that account consistently for hydrogen ionization in the calculation of the electron density. To this end we solve radiative transfer for hydrogen (bound-bound and bound-free transitions) in the two-dimensional models, including the effect of partial frequency redistribution (PRD) in the Lyman (alpha) and (beta) lines. From our internally consistent models we will calculate emergent spectra and the way these vary with location of some well-known spectral diagnostics and compare our results with observed line profiles. We can readily compare theoretical CO profiles from out models with our spatially resolved CO observations. Also we can compare with spatially resolved Ca II (ground based) and Mg II (we have observations done with the UVSP/SMM instrument), and Lyman (alpha) observations that should be available from SUMER/SOHO.

Noyes, Robert W.↗
NASA NTRS2002

Impact of MIE-Resonances on the Atmospheric Absorption of Water Clouds

Clouds strongly modulate radiative transfer processes in the Earth's atmosphere. Studies, which simulate bulk properties of clouds, such as absorption, require methods that accurately account for multiple scattering among individual cloud particles. Multiple scattering processes are well described by MIE-theory, if interacting particles have a spherical shape. This is a good assumption for water droplets. Thus, simulations for water clouds (especially for interactions with solar radiation) usually apply readily available MIE-codes. The presence of different drop-sizes, however, necessitates repetitive calculations for many sizes. The usual representation by a few sizes is likely to miss contributions from densely distributed, sharp resonances. Despite their usually narrow width, integrated over the entire size-spectrum of a cloud droplet distribution, the impact of missed resonances could add up. The consideration of these resonances tends to increase cloud extinction and cloud absorption. This mechanism for a larger (than by MIE-methods predicted) solar absorption has the potential to explain observational evidence of larger than predicted cloud absorption at solar wavelengths. The presentation will address the absorption impact of added resonances for typical properties of water clouds (e.g. drop size distributions, drop concentrations and cloud geometry). Special attention will be given to scenarios with observational evidence of law than simulated solar absorption; particularly if simultaneous measurements of cloud micro- and macrophysical properties are available.

Wiscombe, W.↗
NASA NTRS2006

A Thermodynamically Consistent Damage Model for Advanced Composites

A continuum damage model for the prediction of damage onset and structural collapse of structures manufactured in fiber-reinforced plastic laminates is proposed. The principal damage mechanisms occurring in the longitudinal and transverse directions of a ply are represented by a damage tensor that is fixed in space. Crack closure under load reversal effects are taken into account using damage variables established as a function of the sign of the components of the stress tensor. Damage activation functions based on the LaRC04 failure criteria are used to predict the different damage mechanisms occurring at the ply level. The constitutive damage model is implemented in a finite element code. The objectivity of the numerical model is assured by regularizing the dissipated energy at a material point using Bazant's Crack Band Model. To verify the accuracy of the approach, analyses of coupon specimens were performed, and the numerical predictions were compared with experimental data.

Maimi, Pere↗
DOE OSTI2024

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↗
DOE OSTI2023

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↗
DOE OSTI2016

Experimental and Monte Carlo studies of fluence corrections for graphite calorimetry in low- and high-energy clinical proton beams

Purpose: The aim of this study was to determine fluence corrections necessary to convert absorbed dose to graphite, measured by graphite calorimetry, to absorbed dose to water. Fluence corrections were obtained from experiments and Monte Carlo simulations in low- and high-energy proton beams. Methods: Fluence corrections were calculated to account for the difference in fluence between water and graphite at equivalent depths. Measurements were performed with narrow proton beams. Plane-parallel-plate ionization chambers with a large collecting area compared to the beam diameter were used to intercept the whole beam. High- and low-energy proton beams were provided by a scanning and double scattering delivery system, respectively. A mathematical formalism was established to relate fluence corrections derived from Monte Carlo simulations, using the FLUKA code [A. Ferrari et al., “FLUKA: A multi-particle transport code,” in CERN 2005-10, INFN/TC 05/11, SLAC-R-773 (2005) and T. T. Böhlen et al., “The FLUKA Code: Developments and challenges for high energy and medical applications,” Nucl. Data Sheets 120, 211–214 (2014)], to partial fluence corrections measured experimentally. Results: A good agreement was found between the partial fluence corrections derived by Monte Carlo simulations and those determined experimentally. For a high-energy beam of 180 MeV, the fluence corrections from Monte Carlo simulations were found to increase from 0.99 to 1.04 with depth. In the case of a low-energy beam of 60 MeV, the magnitude of fluence corrections was approximately 0.99 at all depths when calculated in the sensitive area of the chamber used in the experiments. Fluence correction calculations were also performed for a larger area and found to increase from 0.99 at the surface to 1.01 at greater depths. Conclusions: Fluence corrections obtained experimentally are partial fluence corrections because they account for differences in the primary and part of the secondary particle fluence. A correction factor, F(d), has been established to relate fluence corrections defined theoretically to partial fluence corrections derived experimentally. The findings presented here are also relevant to water and tissue-equivalent-plastic materials given their carbon content.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Screening sensitivity analysis of a PWR fuel assembly FEM structural model

Within the present project, computational methods to describe the static mechanical response of PWR fuel assemblies (FA) in the reactor core are developed. For this purpose a structural model of a PWR fuel assembly has been implemented with the Finite Element Method (FEM) using the ANSYS APDL Mechanical code. The purpose of the model is to investigate the long-term creep deformation of the fuel assemblies in the reactor core, commonly known as fuel assembly bow. FA bow is considered to be a complex process with a large number of influencing mechanisms and many only roughly known boundary conditions. Uncertainty and sensitivity analyses are a common way to investigate the predictive power and reliability of such complex models. Hence, as a first approach to understand the importance of the many model parameters, a qualitative sensitivity analysis of the model was conducted using the screening method proposed by Morris. In the first part of the paper, the FA structural model and its principal parameters are presented. The second part presents first the different simulation experiments used for the screening analyses. Based on the model description, relevant input parameters are then identified and a typical range of values is determined. The sensitivity analysis is performed in two subsequent steps. First, the model is screened for influential stiffness parameters with the objective of stiffening the FA structure. Second, the influence of the uncertain boundary conditions and material models on the FA creep deformation is tested. For the considered FA design it is concluded that the influence of the involved uncertainties is comparable to the influence of the parameters which can increase the stiffness. An integral approach accounting for these uncertainties is hence necessary to predict FA bow reliably in order to take additional measures for a further reduction of FA bow amplitudes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Dynamic Succession of Groundwater Sulfate-Reducing Communities during Prolonged Reduction of Uranium in a Contaminated Aquifer

To further understand the diversity and dynamics of SRB in response to substrate amendment, we sequenced genes coding for the dissimilatory sulfite reductase (dsrA) in groundwater samples collected after an emulsified vegetable oil (EVO) amendment, which sustained U(VI)-reducing conditions for one year in a fast-flowing aquifer. EVO amendment significantly altered the composition of groundwater SRB communities. Sequences having no closely related-described species dominated (80%) the indigenous SRB communities in nonamended wells. After EVO amendment, Desulfococcus , Desulfobacterium , and Desulfovibrio , known for long-chain-fatty-acid, short-chain-fatty-acid and H 2 oxidation and U(VI) reduction, became dominant accounting for 7 ± 2%, 21 ± 8%, and 55 ± 8% of the SRB communities, respectively. Succession of these SRB at different bioactivity stages based on redox substrates/products (acetate, SO 4 –2 , U(VI), NO 3 – , Fe(II), and Mn(II)) was observed. Desulfovibrio and Desulfococcus dominated SRB communities at 4–31 days, whereas Desulfobacterium became dominant at 80–140 days. By the end of the experiment (day 269), the abundance of these SRB decreased but the overall diversity of groundwater SRB was still higher than non-EVO controls. Up to 62% of the SRB community changes could be explained by groundwater geochemical variables, including those redox substrates/products. A significant ( P < 0.001) correlation was observed between groundwater U(VI) concentrations and Desulfovibrio abundance. Lastly, our results showed that the members of SRB and their dynamics were correlated significantly with slow EVO biodegradation, electron donor production and maintenance of U(VI)-reducing conditions in the aquifer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Expanding the capability of reaction-diffusion codes using pseudo traps and temperature partitioning: Applied to hydrogen uptake and release from tungsten

Simulating the implantation and thermal desorption evolution in a reaction-diffusion model requires solving a set of coupled differential equations that describe the trapping and release of atomic species in Plasma Facing Materials (PFMs). These fundamental equations are well outlined by the Tritium Migration Analysis Program (TMAP) which can model systems with no more than three active traps per atomic species. To overcome this limitation, we have developed a Pseudo Trap and Temperature Partition (PTTP) scheme allowing us to lump multiple inactive traps into one pseudo trap, simplifying the system of equations to be solved. For all temperatures, we show the trapping of atoms from solute is exactly accounted for when using a pseudo trap. However, a single effective pseudo trap energy can not well replicate the release from multiple traps, each with its own detrapping energy. However, atoms held in a high energy trap will remain trapped at relatively low temperatures, and thus there is a temperature range in which release from high energy traps is effectively inactive. By partitioning the temperature range into segments, a pseudo trap can be defined for each segment to account for multiple high energy traps that are actively trapping but are effectively not releasing atoms. With increasing temperature, as in controlled thermal desorption, the lowest energy trap is nearly emptied and can be removed from the set of coupled equations, while the next higher energy trap becomes an actively releasing trap. Each segment is thus calculated sequentially, with the last time step of a given segment solution being used as an initial input for the next segment as only the pseudo and actively releasing traps are modeled. This PTTP scheme is then applied to experimental thermal desorption data for tungsten (W) samples damaged with heavy ions, which display six distinct release peaks during thermal desorption. Without modifying the TMAP7 source code the PTTP scheme is shown to successfully model the D retention in all six traps. In conclusion, we demonstrate the full reconstruction from the plasma implantation phase through the controlled thermal desorption phase with detrapping energies near 0.9, 1.1, 1.4, 1.7, 1.9 and 2.1 eV for a W sample damaged at room temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗