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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

A Model for the Sources of the Slow Solar Wind

Models for the origin of the slow solar wind must account for two seemingly contradictory observations: The slow wind has the composition of the closed-field corona, implying that it originates from the continuous opening and closing of flux at the boundary between open and closed field. On the other hand, the slow wind has large angular width, up to approximately 60 degrees, suggesting that its source extends far from the open-closed boundary. We propose a model that can explain both observations. The key idea is that the source of the slow wind at the Sun is a network of narrow (possibly singular) open-field corridors that map to a web of separatrices and quasi-separatrix layers in the heliosphere. We compute analytically the topology of an open-field corridor and show that it produces a quasi-separatrix layer in the heliosphere that extends to angles far front the heliospheric current sheet. We then use an MHD code and MIDI/SOHO observations of the photospheric magnetic field to calculate numerically, with high spatial resolution, the quasi-steady solar wind and magnetic field for a time period preceding the August 1, 2008 total solar eclipse. Our numerical results imply that, at least for this time period, a web of separatrices (which we term an S-web) forms with sufficient density and extent in the heliosphere to account for the observed properties of the slow wind. We discuss the implications of our S-web model for the structure and dynamics of the corona and heliosphere, and propose further tests of the model.

Antiochos, Spiro K.↗
DOE OSTI2016

Technical Note: SPEKTR 3.0—A computational tool for x-ray spectrum modeling and analysis

Purpose: A computational toolkit (SPEKTR 3.0) has been developed to calculate x-ray spectra based on the tungsten anode spectral model using interpolating cubic splines (TASMICS) algorithm, updating previous work based on the tungsten anode spectral model using interpolating polynomials (TASMIP) spectral model. The toolkit includes a MATLAB (The Mathworks, Natick, MA) function library and improved user interface (UI) along with an optimization algorithm to match calculated beam quality with measurements. Methods: The SPEKTR code generates x-ray spectra (photons/mm{sup 2}/mAs at 100 cm from the source) using TASMICS as default (with TASMIP as an option) in 1 keV energy bins over beam energies 20–150 kV, extensible to 640 kV using the TASMICS spectra. An optimization tool was implemented to compute the added filtration (Al and W) that provides a best match between calculated and measured x-ray tube output (mGy/mAs or mR/mAs) for individual x-ray tubes that may differ from that assumed in TASMICS or TASMIP and to account for factors such as anode angle. Results: The median percent difference in photon counts for a TASMICS and TASMIP spectrum was 4.15% for tube potentials in the range 30–140 kV with the largest percentage difference arising in the low and high energy bins due to measurement errors in the empirically based TASMIP model and inaccurate polynomial fitting. The optimization tool reported a close agreement between measured and calculated spectra with a Pearson coefficient of 0.98. Conclusions: The computational toolkit, SPEKTR, has been updated to version 3.0, validated against measurements and existing models, and made available as open source code. Video tutorials for the SPEKTR function library, UI, and optimization tool are available.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Modelling type 1 quasar colours in the era of Rubin and Euclid

We construct a parametric SED model which is able to reproduce the average observed SDSS–UKIDSS–WISE quasar colours to within one-tenth of a magnitude across a wide range of redshift (0 < z < 5) and luminosity (–22 > M i > –29). This model is shown to provide accurate predictions for the colours of known quasars which are less luminous than those used to calibrate the model parameters, and also those at higher redshifts z > 5. Using a single parameter, the model encapsulates an up-to-date understanding of the intra-population variance in the rest-frame ultraviolet and optical emission lines of luminous quasars. At fixed redshift, there are systematic changes in the average quasar colours with apparent i-band magnitude, which we find to be well explained by the contribution from the host galaxy and our parametrization of the emission-line properties. By including redshift as an additional free parameter, the model could be used to provide photometric redshifts for individual objects. For the population as a whole we find that the average emission line and host-galaxy contributions can be well described by simple functions of luminosity which account for the observed changes in the average quasar colours across 18.1 < i AB < 21.5. Finally, we use these trends to provide predictions for quasar colours at the luminosities and redshifts which will be probed by the Rubin Observatory LSST and ESA-Euclid wide survey. The model code is applicable to a wide range of upcoming photometric and spectroscopic surveys, and is made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

PORFLOW 6.42.9 Testing and Verification Document

This report is a continuation of the series of PORFLOW QA documents including (Aleman, 2007; Whiteside, 2007, 2008, 2010a,c,b, 2016b,a, 2017a,b). In this report we have updated the test cases where necessary to account for new features in PORFLOW and to ensure those features perform as described in the documentation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Requirements Description of DASSH-F

This report reviews the modeling and simulation capabilities of Argonne National Laboratory’s DASSH code that is used in present reactor analysis activities. These capabilities will be used to establish the set of verification tasks necessary to verify DASSH for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The DASSH program is a thermal analysis code designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. DASSH is a follow-on development to the SE2-ANL software and SUPERENERGY-2 software that it is based upon. DASSH was designed to account for both neutron and gamma heating and is inherently connected to the GAMSOR part of the ARC suite of fast reactor analysis software. SE2-ANL is a developed piece of software from the 1980s while DASSH is a modern implementation with notable improvements in geometry handling. The most important upgrade of DASSH relative to SE2-ANL is that it can analyze multiple time points in a single run where SE2-ANL can only treat a single time point. This allows the user to understand the impact of and search the flow distribution for the entire operational period of a reactor design considering pressure drop, peak coolant and fuel temperatures, and thermal striping. DASSH has three input paths that have to be verified. The first input path builds the geometry and power distribution based upon the DIF3D model but ignores the gamma heating aspects of the problem. The second input path also builds the geometry from the DIF3D model but it takes the neutron and gamma heating distributions from GAMSOR. The third input path is to take the geometry and power distribution directly from user input (i.e. not coupled to DIF3D or GAMSOR). DASSH also has many built in correlations for material properties along with a user defined specification of the fuel, structure, and coolant properties. There are correlations for flow split, mixing, pressure drop, and heat transfer coefficients (subchannel rather than a direct methodology). In total, verification of DASSH will require an extensive testing to cover all possible user features of the software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

THE ROLE OF COSMIC-RAY PRESSURE IN ACCELERATING GALACTIC OUTFLOWS

We study the formation of galactic outflows from supernova (SN) explosions with the moving-mesh code AREPO in a stratified column of gas with a surface density similar to the Milky Way disk at the solar circle. We compare different simulation models for SN placement and energy feedback, including cosmic rays (CRs), and find that models that place SNe in dense gas and account for CR diffusion are able to drive outflows with similar mass loading as obtained from a random placement of SNe with no CRs. Despite this similarity, CR-driven outflows differ in several other key properties including their overall clumpiness and velocity. Moreover, the forces driving these outflows originate in different sources of pressure, with the CR diffusion model relying on non-thermal pressure gradients to create an outflow driven by internal pressure and the random-placement model depending on kinetic pressure gradients to propel a ballistic outflow. CRs therefore appear to be non-negligible physics in the formation of outflows from the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Dependence of Gamma-Ray Burst Jet Collimation on Black Hole Spin

Gamma-ray bursts (GRBs) are the most luminous events in the Universe and are excellent laboratories to study extreme physical phenomena in the cosmos. Despite a long trajectory of progress in understanding these highly energetic events, there are still many observed features that are yet to be fully explained. Observations of the jet opening angle of long gamma-ray bursts (LGRBs) suggest that LGRB jets are narrower for those GRBs at higher redshift. This phenomenon has been explained in the context of collimation by the stellar envelope, with denser (lower metallicity) stars at higher redshifts able to collimate the jet more effectively. However, until now, the dependence of the jet opening angle on the properties of the central engine has not been explored. We investigate the effect of black hole spin on the jet collimation angle for a magnetically launched jet, using the general relativistic radiation magnetohydrodynamical code ν bhlight . We present 3D results for a range of spin values. The simulations show that higher-spinning black holes tend to create narrower jets. If indeed LGRB progenitors in the early Universe are able to produce black hole central engines with higher spin, this could account for at least some of the observed jet opening angle-redshift correlation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Analyzing high dimensional data

Problems encountered in analyzing high dimensional data are discussed and possible solutions are proposed. The increased importance of second-order statistics in analyzing high dimensional data and the shortcoming of the minimum distance classifier in high dimensional data are recognized. By investigating characteristics of high dimensional data, it is shown that second-order statistics must be taken into account in high dimensional data. There is a need to represent second order statistics effectively. As the data dimensionality increases, it becomes more difficult to perceive and compare information present in statistics derived from data. In order to overcome this problem, a method to visualize statistics using color code is proposed. By representing statistics using a color code, the first and the second statistics can be more readily compared.

Lee, Chulhee↗
DOE OSTI2022

Evaluation of the Effect of Burnup Acceleration on UO 2 Microstructure Evolution

Accelerated fuel qualification (AFQ) is a methodology by which new nuclear fuels are developed in an accelerated time frame compared with historical fuel qualification approaches. AFQ generally relies on high-fidelity physics-based modeling and simulation tools to adequately describe fuel performance as well as on revolutionary methods to accelerate burnup accumulation and collect relevant data more quickly. This report summarizes the use of advanced fuel modeling and simulation tools to evaluate microstructures from commercially irradiated fuel and microstructures from proposed MiniFuel irradiations, in which burnup accumulation is accelerated while prototypic temperature conditions are maintained. In this milestone, we used the mesoscale fuel performance code MARMOT to model the evolution of irradiated UO 2 microstructures and their potential restructuring at high burnup. The simulation conditions were informed by BISON models of both commercially irradiated fuel and MiniFuel. A first set of simulations investigated the recrystallization behavior of fully dense microstructures and showed full recrystallization at burnups as low as 52 MWd/kgU at 950°C. However, these simulations did not account for the presence of fission gas bubbles (predicted by BISON). Therefore, a second set of simulations including fission gas bubbles was performed and indicated that at the lowest temperature considered (650°C), the porous UO 2 microstructures have the highest total Gibbs free energies and are likely to recrystallize earlier than higher temperature cases (800 and 950°C), which agrees with high-burnup fuel characterization data. The results also showed that at lower temperature (650°C), the total free energies of the PWR fuel and MiniFuel microstructures are not significantly different. However, at the highest temperature (950°C), MiniFuel microstructures have a lower free energy than that of the PWR fuel microstructure. The competing effects between the temperature-dependent grain nucleation rate and the reduction of the free energy of the microstructure at higher temperature as a result of diffusion indicated that restructuring may occur at even higher temperatures than those considered in this study. In addition to the microstructure evolution modeling efforts, the burnup gradient across a single fuel specimen was also considered. This evaluation was for the VXF-15 position of the High Flux Isotope Reactor (HFIR) using the code suite HFIRCON, which was developed to automate the workflow for evaluating targets and fuel as they are irradiated in HFIR. The burnup gradient evaluation showed a dependence on both the axial and radial locations within the specimen, with a maximum difference of 1.7 between the inner and outermost radial layers. This relationship was further supported by considering the fission product speciation with respect to location within the specimen, which showed a higher concentration of 239 Pu, 240 Pu, and 241 Pu on the outer radial locations of the specimen than the center. The findings of the burnup and speciation evaluation show that some amount of self-shielding is occurring in the specimen when irradiated in the high-flux environment of HFIR; however, this impact is more pronounced for natural uranium when compared to 6% enrichment due to the higher ratio of 238 U in the specimen. Further analyses are required to understand the sensitivity of this gradient to spatial mesh and enrichment of the specimen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Quantum effects on dislocation motion from ring-polymer molecular dynamics

Quantum motion of atoms known as zero-point vibration was recently proposed to explain a long-standing discrepancy between theoretically computed and experimentally measured low-temperature plastic strength of iron and possibly other metals with high atomic masses. This finding challenges the traditional notion that quantum motion of atoms is relatively unimportant in solids comprised of heavy atoms. Here we report quantum dynamic simulations of quantum effects on dislocation motion within the exact formalism of Ring-Polymer Molecular Dynamics (RPMD). To extend the reach of quantum atomistic simulations to length and time scales relevant for extended defects in materials, we implemented RPMD in the open-source code LAMMPS thus making the RPMD method widely available to the community. We use our RPMD/LAMMPS approach for direct calculations of dislocation mobility and its effects on the yield strength of α-iron. Our simulation results establish that quantum effects are noticeable at temperatures below 50 K but account for only a modest (≈13% at T = 0 K) overall reduction in the Peierls barrier, at variance with the factor of two reduction predicted earlier based on the more approximate framework of harmonic transition state theory. Our results confirm that zero-point vibrations provide ample additional agitation for atomic motion that increases with decreasing temperature, however its enhancing effect on dislocation mobility is largely offset by an increase in the effective atom size, an effect known as quantum dispersion that has not been accounted for in the previous calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Implementation of ATLAS Distributed Computing monitoring dashboards using InfluxDB and Grafana

For the last 10 years, the ATLAS Distributed Computing project has based its monitoring infrastructure on a set of custom designed dashboards provided by CERN. This system functioned very well for LHC Runs 1 and 2, but its maintenance has progressively become more difficult and the conditions for Run 3, starting in 2021, will be even more demanding; hence a more standard code base and more automatic operations are needed. A new infrastructure has been provided by CERN, based on InfluxDB as the data store and Grafana as the display environment. ATLAS has adapted and further developed its monitoring tools to use this infrastructure for data and workflow management monitoring and accounting dashboards, expanding the range of previous possibilities with the aim to achieve a single, simpler, environment for all monitoring applications. This document describes these tools and the data flows for monitoring and accounting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Technical Note: Effect of explicit M and N-shell atomic transitions on a low-energy x-ray source

Purpose: In EGSnrc, atomic transitions to and from the M and N-shells are treated in an average way by default. This approach is justified in which the energy difference between explicit and average M and N-shell binding energies is less than 1 keV, and for most applications can be considered negligible. However, for simulations of low energy x-ray sources on thin, high-Z targets, characteristic x-rays can make up a significant portion of the source spectra. As of release V4-2.4.0, EGSnrc has included an option to enable a more complete algorithm of all atomic transitions available in the EADL compilation. In this paper, the effect of M and N-shell averaging on the calculation of half-value layer (HVL) and relative depth dose (RDD) curve of a 50 kVp intraoperative x-ray tube with a thin gold target was investigated. Methods: A 50 kVp miniature x-ray source with a gold target (The INTRABEAM System, Carl Zeiss, Germany) was modeled with the EGSnrc user code cavity, both with and without M and N-shell averaging. From photon fluence spectra simulations, the source HVLs were determined analytically. The same source model was then used with egs-chamber to calculate RDD curves in water. Results: A 4% increase of HVL was reported when accounting for explicit M and N-shell transitions, and up to a 9% decrease in local relative dose for normalization at 3 mm depth in water. Conclusions: The EGSnrc default of using averaged M and N-shell binding energies has an observable effect on the HVL and RDD of a low energy x-ray source with high-Z target. For accurate modeling of this class of devices, explicit atomic transitions should be included.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Studies of Deteriorated Heat Transfer in Prismatic Cores Stemming from Irradiation-Induced Geometry Distortion

A reference design for the Next Generation Nuclear Plant (NGNP) is to use General Atomics Modular High Temperature Gas-cooled Reactor (MHTGR). For such a configuration in normal operation, the helium coolant flow proceeds from the upper plenum to the lower plenum principally through the core coolant channels and the interstitial gaps (bypass flow) that separate the prismatic blocks from one another. Only the core prismatic blocks have coolant channels. The interstitial gaps are present throughout the core, the inner reflector region, and the out reflector region. The bypass flows in a prismatic gas-cooled reactor (GCR) are of potential concern because they reduce the desired flow rates in the coolant channels and, thereby, can increase outlet gas temperatures and maximum fuel temperatures. Consequently, it is appropriate to account for bypass flows in reactor thermal gas dynamic analyses. The objectives of this project include the following: fundamentally understand bypass flow and heat transfer at scaled, undistorted conditions and with geometry distortions; develop improved estimates of associated loss coefficients, surface friction and heat transfer for systems and network codes; and obtain related data for validation of CFD (computational fluid dynamic) or system (e.g., RELAP5) codes which can be employed in predictions for a GCR for normal power, reduced power, and residual heat removal operations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

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

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

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

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

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