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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,189 records · Page 66

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Transient Vibration Prediction for Rotors on Ball Bearings Using Load-dependent Non-linear Bearing Stiffness

Rolling-element bearing forces vary nonlinearly with bearing deflection. Thus an accurate rotordynamic transient analysis requires bearing forces to be determined at each step of the transient solution. Analyses have been carried out to show the effect of accurate bearing transient forces (accounting for non-linear speed and load dependent bearing stiffness) as compared to conventional use of average rolling-element bearing stiffness. Bearing forces were calculated by COBRA-AHS (Computer Optimized Ball and Roller Bearing Analysis - Advanced High Speed) and supplied to the rotordynamics code ARDS (Analysis of Rotor Dynamic Systems) for accurate simulation of rotor transient behavior. COBRA-AHS is a fast-running 5 degree-of-freedom computer code able to calculate high speed rolling-element bearing load-displacement data for radial and angular contact ball bearings and also for cylindrical and tapered roller beatings. Results show that use of nonlinear bearing characteristics is essential for accurate prediction of rotordynamic behavior.

Fleming, David P.↗

NASA NTRS2005

Code for Analyzing and Designing Spacecraft Power System Radiators

GPHRAD is a computer code for analysis and design of disk or circular-sector heat-rejecting radiators for spacecraft power systems. A specific application is for Stirling-cycle/linear-alternator electric-power systems coupled to radioisotope general-purpose heat sources. GPHRAD affords capabilities and options to account for thermophysical properties (thermal conductivity, density) of either metal-alloy or composite radiator materials.

Juhasz, Albert↗

NASA NTRS2014

Implementing Badhwar-O'Neill Galactic Cosmic Ray Model for the Analysis of Space Radiation Exposure

For the analysis of radiation risks to astronauts and planning exploratory space missions, accurate energy spectrum of galactic cosmic radiation (GCR) is necessary. Characterization of the ionizing radiation environment is challenging because the interplanetary plasma and radiation fields are modulated by solar disturbances and the radiation doses received by astronauts in interplanetary space are likewise influenced. A model of the Badhwar‐O'Neill 2011 (BO11) GCR environment, which is represented by GCR deceleration potential theta, has been derived by utilizing all of the GCR measurements from balloons, satellites, and the newer NASA Advanced Composition Explorer (ACE). In the BO11 model, the solar modulation level is derived from the mean international sunspot numbers with time‐delay, which has been calibrated with actual flight instrument measurements to produce better GCR flux data fit during solar minima. GCR fluxes provided by the BO11 model were compared with various spacecraft measurements at 1 AU, and further comparisons were made for the tissue equivalent proportional counters measurements at low Earth orbits using the high‐charge and energy transport (HZETRN) code and various GCR models. For the comparison of the absorbed dose and dose equivalent calculations with the measurements by Radiation Assessment Detector (RAD) at Gale crater on Mars, the intensities and energies of GCR entering the heliosphere were calculated by using the BO11 model, which accounts for time‐dependent attenuation of the local interstellar spectrum of each element. The BO11 model, which has emphasized for the last 24 solar minima, showed in relatively good agreement with the RAD data for the first 200 sols, but it was resulted in to be less well during near the solar maximum of solar cycle 24 due to subtleties in the changing heliospheric conditions. By performing the error analysis of the BO11 model and the optimization in reducing overall uncertainty, the resultant BO13 model corrects the fit at solar maxima as well as being accurate at solar minima. The BO13 model is implemented to the NASA Space Cancer Risk model for the assessment of radiation risks. Overall cumulative probability distribution of solar modulation parameters represents the percentile rank of the average interplanetary GCR environment, and the probabilistic radiation risks can be assessed for various levels of GCR environment to support mission design and operational planning for future manned space exploration missions.

Kim, Myung-Hee Y.↗

NASA NTRS2021

Magnetosphere-Ionosphere Coupling of Precipitated Electrons in Diffuse Aurora Driven by Time Domain Structures

Recent theoretical studies and the spacecraft conjugate observations between the Time History of Events and Macroscale Interactions during Substorms (THEMIS) mission and the low-altitude Enhanced Polar Outflow Probe (e-POP) spacecraft demonstrated the connections between broadband electrostatic fluctuations, the so called time domain structures (TDSs), and electron precipitation in the region of diffuse aurora. In this letter, we used the SuperThermal Electron Transport (STET) code to implement these theoretical and experimental results and analyze magnetosphere-ionosphere energy interplay of the precipitated electrons that are driven by TDSs. To put TDSs electron scattering processes in the context with other wave scattering activities, ECH and whistler chorus waves are taken into account. Integrated electron energy fluxes in the diffuse aurora driven by TDSs are comparable to those driven by whistler waves. We find that including Magnetosphere-Ionosphere coupling processes increases the integrated electron energy fluxes by more than two times than only driven by TDSs.

diffuse aurora↗

DOE OSTI2015

Reassessment of uncertainties for the PU-MET-INTER-002 benchmark - 14432

Uncertainty for the PU-MET-INTER-002 benchmark of ZPR-6/10 was unusually large for a ZPR assembly and was dominated by the uncertainties in matrix tube pitch and composition uncertainties for plutonium, graphite and stainless steel plates. The original uncertainty analysis for PU-MET-INTER-002 was limited by computer resources, software and the nuclear data available in 2002, all of which have been improved since. An as-built Monte Carlo model of ZPR-6/10 was used to quantify components of the transformation from the detailed model to a simplified RZ benchmark in 2002. Compositional uncertainties were estimated based upon TWODANT calculations of the RZ benchmark model (20 group, S{sub 8} - P{sub 1}). MC{sup 2}-2 was used with END/B-V.2 data to generate homogenized cross section data for TWODANT. The MC{sup 2}-2 calculations did not account for lattice cell spatial heterogeneity which was driven by the manual process required to generate ∼80 modifications to the deterministic model and by the software limitations of MC{sup 2}-2. Improvements to DIF3D-VARIANT, development of the MC{sup 2}-3 and BuildZPRModel codes, and better computer resources make accurate 3D transport calculations of reactivity changes practical today. BuildZPRModel imports the as-built geometry and creates a 'drawer' homogenized DIF3D model where each radial DIF3D mesh is the exact size of a single ZPR matrix tube and drawer. A MC{sup 2}-3 input is constructed for each drawer to produce homogenized cross section data where the spatial cell heterogeneity is captured by solving a 1D MOC transport calculation on the ultrafine energy group level. Collapsing 3D drawers into 1D MC{sup 2}-3 drawer models follows conventional procedures and is done by BuildZPRModel. MCNP6 was used to re-compute relevant geometry components of the uncertainty using ENDF/B-VII.0 data with the as-built model. Changes in the geometry components compared with earlier results are insignificant. DIF3D calculations of the compositional uncertainties were performed using 33, 70, 116, 230, 425, and 703 group energy structures and selectively combined with P{sub 3}-P{sub 7} transport approximations and P{sub 1} and P{sub 3} scattering kernels. 3D diffusion theory calculations were also performed for comparison. Many composition uncertainty components were essentially unchanged with the updated methodology. However, the dominant experimental uncertainty for PU-MET-INTER-002 is the composition of stainless steel. The new calculations reduce the stainless steel uncertainty from ∼115 pcm to ∼70 pcm which significantly reduces the total uncertainty for PU-MET-INTER-002. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1994

Large Eddy simulation of turbulence: A subgrid scale model including shear, vorticity, rotation, and buoyancy

The Reynolds numbers that characterize geophysical and astrophysical turbulence (Re approximately equals 10(exp 8) for the planetary boundary layer and Re approximately equals 10(exp 14) for the Sun's interior) are too large to allow a direct numerical simulation (DNS) of the fundamental Navier-Stokes and temperature equations. In fact, the spatial number of grid points N approximately Re(exp 9/4) exceeds the computational capability of today's supercomputers. Alternative treatments are the ensemble-time average approach, and/or the volume average approach. Since the first method (Reynolds stress approach) is largely analytical, the resulting turbulence equations entail manageable computational requirements and can thus be linked to a stellar evolutionary code or, in the geophysical case, to general circulation models. In the volume average approach, one carries out a large eddy simulation (LES) which resolves numerically the largest scales, while the unresolved scales must be treated theoretically with a subgrid scale model (SGS). Contrary to the ensemble average approach, the LES+SGS approach has considerable computational requirements. Even if this prevents (for the time being) a LES+SGS model to be linked to stellar or geophysical codes, it is still of the greatest relevance as an 'experimental tool' to be used, inter alia, to improve the parameterizations needed in the ensemble average approach. Such a methodology has been successfully adopted in studies of the convective planetary boundary layer. Experienc e with the LES+SGS approach from different fields has shown that its reliability depends on the healthiness of the SGS model for numerical stability as well as for physical completeness. At present, the most widely used SGS model, the Smagorinsky model, accounts for the effect of the shear induced by the large resolved scales on the unresolved scales but does not account for the effects of buoyancy, anisotropy, rotation, and stable stratification. The latter phenomenon, which affects both geophysical and astrophysical turbulence (e.g., oceanic structure and convective overshooting in stars), has been singularly difficult to account for in turbulence modeling. For example, the widely used model of Deardorff has not been confirmed by recent LES results. As of today, there is no SGS model capable of incorporating buoyancy, rotation, shear, anistropy, and stable stratification (gravity waves). In this paper, we construct such a model which we call CM (complete model). We also present a hierarchy of simpler algebraic models (called AM) of varying complexity. Finally, we present a set of models which are simplified even further (called SM), the simplest of which is the Smagorinsky-Lilly model. The incorporation of these models into the presently available LES codes should begin with the SM, to be followed by the AM and finally by the CM.

Canuto, V. M.↗

DOE OSTI2019

Adaptive mesh refinement for steady flows in Nek5000

Adaptive mesh refinement is performed in the framework of the spectral element method augmented by approaches to error estimation and control. The h-refinement technique is used for adapting the mesh, where selected grid elements are split by a quadtree (2D) or octree (3D) structure. Continuity between parent-child elements is enforced by high-order interpolation of the solution across the common faces. Parallel mesh partitioning and grid management respectively, are taken care of by the external libraries ParMETIS and p4est. Two methods are considered for estimating and controlling the error of the solution. The first error estimate is local and based on the spectral properties of the solution on each element. This method gives a local measure of the L-2-norm of the solution over the entire computational domain. The second error estimate uses the dual-weighted residuals method - it is based on and takes into account both the local properties of the solution and the global dependence of the error in the solution via an adjoint problem. The objective of this second approach is to optimize the computation of a given functional of physical interest. The simulations are performed by using the code Nek5000 and three steady-state test cases are studied: a two-dimensional lid-driven cavity at Re = 7, 500, a two-dimensional flow past a cylinder at Re = 40, and a three-dimensional lid-driven cavity at Re = 2,000 with a moving lid tilted by an angle of 30 degrees. Finally, the efficiency of both error estimators is compared in terms of refinement patterns and accuracy on the functional of interest. In the case of the adjoint error estimators, the trend on the error of the functional is shown to be correctly represented up to a multiplicative constant.

42 ENGINEERING↗

DOE OSTI2016

TH-CD-BRA-03: Direct Measurement of Magnetic Field Correction Factors, KQB, for Application in Future Codes of Practice for Reference Dosimetry

Purpose: With an MR-linac, radiation is delivered in the presence of a magnetic field. Modifications in the codes of practice (CoPs) for reference dosimetry are required to incorporate the effect of the magnetic field. Methods: In most CoPs the absorbed dose is determined using the well-known kQ formalism as the product of the calibration coefficient, the corrected electrometer reading and kQ, to account for the difference in beam quality. To keep a similar formalism a single correction factor is introduced which replaces kQ, and which corrects for beam quality and B-field, kQ,B. In this study we propose a method to determine kQ,B under reference conditions in the MRLinac without using a primary standard, as the product of:- the ratio between detector readings without and with B-field (kB),- the ratio between doses in the point of measurement with and without B-field (rho),- kQ in the absence of the B-field in the MRLinac beam (kQmrl0,Q0),The ratio of the readings, which covers the change in detector reading due to the different electron trajectories in the detector, was measured with a waterproof ionization chamber (IBA-FC65g) in a water phantom in the MRLinac without and with B-field. The change in dose-to-water in the point of measurement due to the B-field was determined with a Monte Carlo based TPS. Results: For the presented approach, the measured ratio of readings is 0.956, the calculated ratio of doses in the point of measurement is 0.995. Based on TPR20,10 measurements kQ was calculated as 0.989 using NCS-18. This yields a value of 0.9408 for kQ,B. Conclusion: The presented approach to determine kQ,B agrees with a method based on primary standards within 0.4% with an uncertainty of 1% (1 std.uncert). It differs from a similar approach using a PMMA-phantom and an NE2571 chamber with 1.3%.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

A New Equivalence Theory Method for Treating Doubly Heterogeneous Fuel - I. Theory

A new methodology has been developed to treat resonance self-shielding in doubly heterogeneous very high temperature gas-cooled reactor systems in which the fuel compact region of a reactor lattice consists of small fuel grains dispersed in a graphite matrix. This new method first homogenizes the fuel grain and matrix materials using an analytically derived disadvantage factor from a two-region problem with equivalence theory and intermediate resonance method. This disadvantage factor accounts for spatial self-shielding effects inside each grain within the framework of an infinite array of grains. Then the homogenized fuel compact is self-shielded using a Bondarenko method to account for interactions between the fuel compact regions in the fuel lattice. In the final form of the equations for actual implementations, the double-heterogeneity effects are accounted for by simply using a modified definition of a background cross section, which includes geometry parameters and cross sections for both the grain and fuel compact regions. With the new method, the doubly heterogeneous resonance self-shielding effect can be treated easily even with legacy codes programmed only for a singly heterogeneous system by simple modifications in the background cross section for resonance integral interpolations. This paper presents a detailed derivation of the new method and a sensitivity study of double-heterogeneity parameters introduced during the derivation. The implementation of the method and verification results for various test cases are presented in the companion paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Enskog Transport Calculations with Magpie

he revised Enskog theory for mixtures has been implemented, as a library, in the LANL institutional thermochemical code Magpie. This report provides the necessary physical background and several representative examples of transport calculations with this library. The particular attention is paid to how the contributions of non-translational degrees of freedom (molecular vibrations, rotations, electronic excitations) to the energy transport can be consistently accounted for via the so called Eucken correction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Energy Deposition Analysis of VERA Progression Problems using MCNP6

The current energy deposition model in the CASL neutronics code MPACT assumes all the fission energy is deposited locally. This approximation limits the accurate represen tation of the heating source distribution in the reactor core, which is essential for the thermal/hydraulic coupling calculations. There has been ongoing work to develop im proved models in MPACT to account for the energy deposition distribution of gamma rays, including direct fission gamma (prompt and delayed) and the gamma rays from neutron capture. Also, the neutron/photon coupling transport calculation has been considered as a potential capability to be added in MPACT to provide high-fidelity heating calculations. In order to provide a reference solution for these gamma heating models, MCNP has been used in this work to develop a set of energy deposition bench mark problems, based on the elective VERA Progression Problems, including single pin cells, 2-D assemblies and 2-D quarter core cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

A Flexible Tool for Beam Induced Background Simulations at a Muon Collider

A Muon Collider represents a very interesting possibility for a future machine to explore the energy frontier in particle physics.However, to reach the needed luminosity, beam intensities of the order of $10^{9-11}$ muons per bunch are needed. In this context, the Beam Induced Background must be taken into account for its effects on the magnets and on the detector. Several optimizations can be conceived with the aim to mitigate them. In this view, it is of crucial importance to develop a flexible tool that allows to easily reconstruct the machine geometry in a Monte Carlo code, allowing to simulate in detail the interaction of muon decay products in the machine, while being able to change the machine optics itself to find the best configuration.In this contribution, a possible approach to such a purpose is presented, based on FLUKA for the Monte Carlo simulation and on LineBuilder for the geometry reconstruction.First results for the 1.5TeV machine optics developed by the MAP collaboration are discussed in this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1990

Analysis of coolant flow and heat transfer in the SSME HPOTP Number 4 bearing assembly

The PHOENICS code has been applied to simulate the flow of liquid oxygen through the number 4 ball bearing assembly of the Space Shuttle Main Engine High-Pressure Oxidizer Turbopump. The Body-Fitted Coordinate system capability of PHOENICS was utilized to create a geometrically accurate model. Effects accounted for by the model include the rotation rate of the calculation domain, viscous heating in the liquid oxygen and two-phase effects due to LOX boiling. A separation program was used to account for conduction within a ball bearing element which was coupled with the PHOENICS flow and heat transfer analysis. Solutions have been obtained for the velocity and temperature fields within the LOX coolant and the temperature within the ball element. The predicted ball temperatures indicate that the coolant is boiling in the region near the ball surface.

Owens, S. F.↗

NASA NTRS2021

Assessment and Improvement of RANS-based Transition Models based on Experimental Data of the Common Research Model with Natural Laminar Flow

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework often rely upon the correlations that were derived from a limited number of low-speed experiments and do not account for all of the transition mechanisms and/or their variation with the significant flow parameters. Available data from a recent experiment in the National Transonic Facility at the NASA Langley Research Center are used to assess the current transition modeling capability in NASA's OVERFLOW 2.2o code for a swept wing configuration at transonic cruise conditions. Specifically, the OVERFLOW solutions are used together with detailed stability analysis of the boundary layer flow over the new Common Research Model with Natural Laminar Flow (CRM-NLF) to evaluate the accuracy and the robustness of the transport-equation-based transition models, with the goal of proposing improvements that would help to strengthen the physical basis of these models for the important class of flows involving the combined effects of crossflow and flow compressibility. Results highlight the significant underprediction of the laminar flow extent within the inboard region of the wing, wherein the onset of transition may be attributed to a gradual amplification of Tollmien-Schlichting instabilities.

Boundary layer transition↗

NASA NTRSDate not supplied

Application of the NASA Multiscale Analysis Tool: Multiscale Integration and Interoperability

In order to demonstrate NASMAT’s multiscale operability, a series of illustrative examples will be presented that focus on the application of NASMAT to practical problems. First, the multiscale integration and data recursion is demonstrated by performing a multiscale analysis using only built-in micromechanics methods. NASMAT’s integration is then highlighted by running a multiscale analysis where an external finite element software calls NASMAT. In this case, at each integration point within the finite element model, a local NASMAT analysis is performed to account for failure behavior at the constituent scale. In a similar example, an external program is called from within NASMAT. This case would be relevant for a user wanting to implement an outside micromechanics technique. A combination of these examples is then presented to further illustrate the code’s flexibility when interfacing with outside codes in a multiscale framework. For all examples, data is presented using a custom-developed visualization tool. Additional potential use cases are also addressed. Finally, the plan for upcoming features and added capabilities is discussed.

NASMAT↗

DOE OSTI2017

Fragmentation Prediction for a Brittle Ceramic: An Integrated Two-Scale Model Approach

A new implementation of a multi-scale fragmentation model is introduced that combines existing independent macro and micro-scale models into a single material model within the Library of Advanced Materials for Engineering (LAME) that can be used directly by the finite element structural mechanics code SIERRA/SM, developed at Sandia National Laboratories. The implementation not only joins existing independent fragmentation models into a single framework, but also enables further development by virtue of being integrated into SIERRA/SM. Integration with a finite element structural mechanics code for a true two-way coupling of the fragmentation models that were previously independent. Hereafter, the prior implementation will be referred to as the sequentially applied fragmentation model. In contrast, the newly developed model will be called the integrated two-scale model. The sequentially applied fragmentation model consists of two independent models applied sequentially to account for the macro-scale and microscale fragmentation. The macro-scale fragmentation behavior is calculated using an elastic gradient damage explicit material model developed to model fracture in brittle materials. The micro-scale fragmentation model is an implementation of the 1D cohesive-zone model developed by Zhou, et al. The micro-scale model is implemented by post-processing the macro-scale data that provides initial conditions for the 1D analysis. The micro-scale model is necessary to predict the resulting fragment size distribution at length scales inaccessible to the finite element based macro-scale model. The full range of fragmentation behavior is then predicted by combining the results from the two independent models. The sequentially applied fragmentation model was used to predict the fragments that result from impact loading of a brittle material, which in that case was a ceramic cylinder of commercial UO{sub 2} nuclear fuel. Understanding the fragmentation behavior of such materials is important for evaluating the risk of airborne particulate release under impact loading scenarios. This work will verify the implementation of the integrated two-scale model using the same experiment modeled by the existing sequentially applied model. The purpose of this manuscript is to introduce the new integrated two-scale fragmentation model and validate the implementation against the existing sequentially applied model for a specific case of impact loading of a brittle material. This validation case will not fully utilize the benefits of having the model integrated into the finite element code as a material model. Rather, this study will demonstrate consistency between the two fragmentation models to ensure that the success of the sequentially applied model will be duplicated by the integrated two-scale approach. Once validated, further development of the model can be made by utilizing the integration into the finite element code to create a fully two-way coupled fragmentation model that will be extended to include higher fidelity models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2011

A Model fot 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 also has large angular width, up to approx.60deg, 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 from the heliospheric current sheet. We then use an MHD code and MDI/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 2008 August 1 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. Key words: solar wind - Sun: corona - Sun: magnetic topology

Antiochos, S. K.↗

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