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,225 records · Page 68

Material Control and Accounting for Liquid-Fueled Molten Salt Reactors: Material Control and Holdup Considerations

The US Nuclear Regulatory Commission (NRC) will likely require license applicants for liquid-fueled molten salt reactors (MSRs) with circulating fuel to submit a nuclear material control and accounting (MC&A) plan or detailed MC&A program description for the facility. In liquid-fueled MSRs with special nuclear material (SNM) in bulk (i.e., not in discrete items) form and rapidly changing quantities due to fuel transmutation and depletion, using traditional nuclear material accounting methods with material balance evaluations is challenging. In reactors with changing inventories, these expected quantities of SNM must be calculated based on operational parameters. Reducing uncertainties on these expected quantities is challenging in the case of MSRs without decades of operational experience to verify and validate predictive computational codes. Moreover, many areas in MSRs are inaccessible because of high-temperature and high-radiation environments, making measurements challenging.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

LDEF geometry/mass model for radiation analyses

A 3-D geometry/mass model of Long Duration Exposure Facility (LDEF) was generated to aid in the interpretation of ionizing radiation experiments relative to the influence of varying shielding distributions around the dosimetry, and to allow more definitive calculations and comparisons with the measured data. This model takes into account the major individual structural members of the LDEF spacecraft, the mass in each experiment tray, and, for selected trays containing ionizing radiation dosimetry, major parts within the tray are modeled. The geometry/mass model, together with a ray tracing algorithms, was programmed for use both as a stand alone code in determining 3-D shielding distributions at dosimetry locations, and as a geometry module that can be interfaced with radiation transport codes.

Colborn, B. L.↗

Technical report series on global modeling and data assimilation. Volume 3: An efficient thermal infrared radiation parameterization for use in general circulation models

A detailed description of a parameterization for thermal infrared radiative transfer designed specifically for use in global climate models is presented. The parameterization includes the effects of the main absorbers of terrestrial radiation: water vapor, carbon dioxide, and ozone. While being computationally efficient, the schemes compute very accurately the clear-sky fluxes and cooling rates from the Earth's surface to 0.01 mb. This combination of accuracy and speed makes the parameterization suitable for both tropospheric and middle atmospheric modeling applications. Since no transmittances are precomputed the atmospheric layers and the vertical distribution of the absorbers may be freely specified. The scheme can also account for any vertical distribution of fractional cloudiness with arbitrary optical thickness. These features make the parameterization very flexible and extremely well suited for use in climate modeling studies. In addition, the numerics and the FORTRAN implementation have been carefully designed to conserve both memory and computer time. This code should be particularly attractive to those contemplating long-term climate simulations, wishing to model the middle atmosphere, or planning to use a large number of levels in the vertical.

Suarex, Max J.↗

On Wave Processes in the Solar Atmosphere

This grant was awarded by NASA/MSFC to The University of Alabama in Huntsville (UAH) to investigate the physical processes responsible for heating and wind acceleration in the solar atmosphere, and to construct theoretical, self-consistent and time-dependent solar wind models based on the momentum deposition by finite amplitude and nonlinear Alfven waves. In summary, there are three main goals of the proposed research: (1) Calculate the wave energy spectra and wave energy fluxes carried by magnetic non- magnetic waves. (2) Find out which mechanism dominates in supplying the wave energy to different parts of the solar atmosphere. (3) Use the results obtained in (1) and (2) to construct theoretical, self-consistent and time- dependent models of the solar wind. We have completed the first goal by calculating the amount of non-radiative energy generated in the solar convection zone as acoustic waves and as magnetic tube waves. To calculate the amount of wave energy carried by acoustic waves, we have used the Lighthill-Stein theory for sound generation modified by Musielak, Rosner, Stein & Ulmschneider (1994). The acoustic wave energy fluxes for stars located in different regions of the Hertzsprung-Russell (H-R) diagram have also been computed. The wave energy fluxes carried by longitudinal and transverse waves along magnetic flux tubes have been calculated by using both analytical and numerical methods. Our analytical approach is based a theory developed by Musielak, Rosner & Ulmschnelder and Musielak, Rosner, Gall & Ulmschneider, which allows computing the wave energy fluxes for linear tube waves. A numerical approach has been developed by Huang, Musielak & Ulmschneider and Ulmschneider & Musielak to compute the energy fluxes for nonlinear tube waves. Both methods have been used to calculate the wave energy fluxes for stars located in different regions of the HR diagram (Musielak, Rosner & Ulmschneider 1998; Ulmschneider, Musielak & Fawzy 1998). Having obtained the wave energy fluxes for acoustic and magnetic tube waves, we have investigated the behavior of these waves in the solar and stellar atmospheres. The results of our extensive studies have been published in many papers and presented at numerous scientific meetings. In these studies we have investigated different aspects of propagation of acoustic and magnetic waves, the efficiency of energy transfer along magnetic structures in the solar atmosphere, and behavior of Alfven waves in stgeady and expanding solar and stellar atmospheres. Recently, we have used some of these results to construct first purely theoretical, two component and time-dependent models of solar and stellar chromospheres. Finally, to address the third goal, we have constructed first fully theoretical, self-consistent and time dependent wind models based on the momentum deposition by non-linear Alfven waves. The full set of single-fluid MHD equations with the background flow has been solved by using a modified version of the ZEUS MHD code. The constructed wind models are radially symmetric with the magnetic field decreasing radially and the initial outflow is described by the standard Parker wind solution. In contrast to previous studies, no assumptions regarding wave linearity, wave damping, and wave-flow interaction are made; the models thus naturally account for the backreaction of the wind on the waves as well as for the nonlinear interaction between different types of MHD waves. The models have been used to explain the origin of fast speed streams in solar coronal holes. The obtained results clearly demonstrate that the momentum deposition by Alfven waves in the solar wind can be sufficient to explain the origin of fast stream components of the solar wind. The range of wave amplitudes required to obtain the desired results seems to be in good agreement with recent observations.

Musielak, Z. E.↗

High Speed Civil Transport (HSCT) Isolated Nacelle Transonic Boattail Drag Study and Results Using Computational Fluid Dynamics (CFD)

Nozzle boattail drag is significant for the High Speed Civil Transport (HSCT) and can be as high as 25% of the overall propulsion system thrust at transonic conditions. Thus, nozzle boattail drag has the potential to create a thrust-drag pinch and can reduce HSCT aircraft aerodynamic efficiencies at transonic operating conditions. In order to accurately predict HSCT performance, it is imperative that nozzle boattail drag be accurately predicted. Previous methods to predict HSCT nozzle boattail drag were suspect in the transonic regime. In addition, previous prediction methods were unable to account for complex nozzle geometry and were not flexible enough for engine cycle trade studies. A computational fluid dynamics (CFD) effort was conducted by NASA and McDonnell Douglas to evaluate the magnitude and characteristics of HSCT nozzle boattail drag at transonic conditions. A team of engineers used various CFD codes and provided consistent, accurate boattail drag coefficient predictions for a family of HSCT nozzle configurations. The CFD results were incorporated into a nozzle drag database that encompassed the entire HSCT flight regime and provided the basis for an accurate and flexible prediction methodology.

Midea, Anthony C.↗

High Speed Civil Transport (HSCT) Isolated Nacelle Transonic Boattail Drag Study and Results Using Computational Fluid Dynamics (CFD)

Nozzle boattail drag is significant for the High Speed Civil Transport (HSCT) and can be as high as 25 percent of the overall propulsion system thrust at transonic conditions. Thus, nozzle boattail drag has the potential to create a thrust drag pinch and can reduce HSCT aircraft aerodynamic efficiencies at transonic operating conditions. In order to accurately predict HSCT performance, it is imperative that nozzle boattail drag be accurately predicted. Previous methods to predict HSCT nozzle boattail drag were suspect in the transonic regime. In addition, previous prediction methods were unable to account for complex nozzle geometry and were not flexible enough for engine cycle trade studies. A computational fluid dynamics (CFD) effort was conducted by NASA and McDonnell Douglas to evaluate the magnitude and characteristics of HSCT nozzle boattail drag at transonic conditions. A team of engineers used various CFD codes and provided consistent, accurate boattail drag coefficient predictions for a family of HSCT nozzle configurations. The CFD results were incorporated into a nozzle drag database that encompassed the entire HSCT flight regime and provided the basis for an accurate and flexible prediction methodology.

Midea, Anthony C.↗

A Comparison of EAST Shock-Tube Radiation Measurements with a New Air Radiation Model

This paper presents a comparison between the recent EAST shock tube radiation measurements (Grinstead et al., AIAA 2008-1244) and the HARA radiation model. The equilibrium and nonequilibrium radiation measurements are studied for conditions relevant to lunar-return shock-layers; specifically shock velocities ranging from 9 to 11 kilometers per second at initial pressures of 0.1 and 0.3 Torr. The simulated shock-tube flow is assumed one-dimensional and is calculated using the LAURA code, while a detailed nonequilibrium radiation prediction is obtained in an uncoupled manner from the HARA code. The measured and predicted intensities are separated into several spectral ranges to isolate significant spectral features, mainly strong atomic line multiplets. The equations and physical data required for the prediction of these strong atomic lines are reviewed and their uncertainties identified. The 700-1020 nm wavelength range, which accounts for roughly 30% of the radiative flux to a peak-heating lunar return shock-layer, is studied in detail and the measurements and predictions are shown to agree within 15% in equilibrium. The plus or minus 1.5% uncertainty on the measured shock velocity is shown to cause up to a plus or minus 30% difference in the predicted radiation. This band of predictions contains the measured values in almost all cases. For the highly nonequilibrium 0.1 Torr cases, the nonequilibrium radiation peaks are under-predicted by about half. This under-prediction is considered acceptable when compared to the order-of-magnitude over-prediction obtained using a Boltzmann population of electronic states. The reasonable comparison in the nonequilibrium regions provides validation for both the non-Boltzmann modeling in HARA and the thermochemical nonequilibrium modeling in LAURA. The N2 (+)(1-) and N2(2+) molecular band systems are studied in the 290 480 nm wavelength range for both equilibrium and nonequilibrium regimes. The non-Boltzmann rate models for these systems, which have significant uncertainties, are tuned to improve the comparison with measurements.

Johnston, Christopher O.↗

Rossi and Feynman-alpha formulas including prompt neutron decay and detector die away time constants - 25396

A simple Rossi-alpha formula and the corresponding Feynman-alpha formula have been derived that take into account two time constants; one characterizing the prompt neutron decay in the fissile source domain and the other the detector die-away time. The derivation of the formula is an extension of the classic heuristic method. A test problem is defined and simulated using the Monte Carlo MCNPX code in order to illustrate the applicability of the formulas. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

SU-G-IeP3-12: Preliminary Report On the Experience of Patient Radiation Dose Monitoring and Tracking Systems; PEMNET, Radimetrics and DoseWatch

Purpose: Three patient radiation dose monitoring and tracking (PRDMT) systems have been in operation at this institution for the past 6 months. There are useful information that should be disseminated to those who are considering installation of PRDMT programs. In addition, there are “problems” uncovered in the process of estimating fluoroscopic “peak” skin dose (PSD), especially, for those patients who received interventional angiographic studies and in conjunction with surgical procedures. Methods: Upon exporting the PRDMT data to Microsoft Excel program, the peak skin dose can be estimated by applying various correction factors including; attenuation due to the tabletop and examination mattress, table height, tabletop translation, backscatter, etc. A procedure was established to screen and divide the PRDMT reported radiation dose and estimated PSD to three different levels of threshold to assess the potential skin injuries, to assist patient follow-up, risk management and provide radiation dosimetry information in case of “Sentinel Event”. Results: The Radiation Dose Structured Report (RDSR) was found to be the prerequisite for the PRDMT systems to work seamlessly. And, the geometrical parameters (gantry and table orientation) displayed by the equipment are not necessarily implemented in the “patient centric” manner which could result in a large error in the PSD estimation. Since, the PRDMT systems obtain their pertinent data from the DICOM tags including the polarity (+ and − signs), the geometrical parameters need to be verified. Conclusion: PRDMT systems provide a more accurate PSD estimation than previously possible as the air-kerma-area dose meter become widely implemented. However, care should be exercised to correctly apply the geometrical parameters in estimating the patient dose. In addition, further refinement is necessary for these software programs to account for all geometrical parameters such as the tabletop translation in the z-direction in particular.

60 APPLIED LIFE SCIENCES↗

Influence of laser induced hot electrons on the threshold for shock ignition of fusion reactions

The effects of Hot Electrons (HEs) generated by the nonlinear Laser-Plasma Interaction (LPI) on the dynamics of Shock Ignition Inertial Confinement Fusion targets are investigated. The coupling between the laser beam, plasma dynamics and hot electron generation and propagation is described with a radiative hydrodynamics code using an inline model based on Paraxial Complex Geometrical Optics [Colaïtis et al., Phys. Rev. E 92, 041101 (2015)]. Two targets are considered: the pure-DT HiPER target and a CH-DT design with baseline spike powers of the order of 200–300 TW. In both cases, accounting for the LPI-generated HEs leads to non-igniting targets when using the baseline spike powers. While HEs are found to increase the ignitor shock pressure, they also preheat the bulk of the imploding shell, notably causing its expansion and contamination of the hotspot with the dense shell material before the time of shock convergence. The associated increase in hotspot mass (i) increases the ignitor shock pressure required to ignite the fusion reactions and (ii) significantly increases the power losses through Bremsstrahlung X-ray radiation, thus rapidly cooling the hotspot. These effects are less prominent for the CH-DT target where the plastic ablator shields the lower energy LPI-HE spectrum. Simulations using higher laser spike powers of 500 TW suggest that the CH-DT capsule marginally ignites, with an ignition window width significantly smaller than without LPI-HEs, and with three quarters of the baseline target yield. The latter effect arises from the relation between the shock launching time and the shell areal density, which becomes relevant in presence of a LPI-HE preheating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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↗

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↗

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↗

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↗

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.↗