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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Evaluating the 238 U PFNS Including Chi-Nu Experimental Data

This report documents an evaluation of 238 U prompt fission neutron spectra (PFNS) which is a deliverable for a FY2024 Q4 NCSP (Nuclear Criticality Safety Program) milestone. This evaluation is new; its prior input is based on extended Los Alamos and exciton models implemented in the code CoH. Experimental covariances were estimated for five experimental data sets. One of these data sets that was measured by the Chi-Nu team of LANL and LLNL. It covers the 238 U PFNS for continuous incident-neutron energies of 1–20 MeV and outgoing-neutron energies from 10 keV– 10 MeV with high precision. Contrary to Chi-Nu data, previous data sets were measured in a limited energy range. The resulting evaluated data correspond well to the experimental PFNS taken into account for the evaluation. The evaluated PFNS also produce average mean energies in agreement with associated Chi-Nu data. If one uses the new evaluated data to predict the neutron multiplication factor, k eff , of the Flattop, Flattop-Pu and BigTen ICSBEP critical assemblies (which all have thick reflectors with high percentages of 238 U), the differences of simulated values compared to those using ENDF/B-VIII.1β3 is modest (less than 25 pcm). In addition to that, the new PFNS predict on average 238 U LLNL pulsed-sphere neutron-leakage spectra slightly better than ENDF/BVIIII.0 and ENDF/B-VIII.1β3 PFNS. The differences are, however, well within the experimental uncertainties.

238U↗
NASA NTRS1990

Nonlinear analysis of compressively loaded linked-plate structures

A method is presented for the geometrically nonlinear analysis of the equilibrium behavior of a general class of compressively loaded, prismatic, linked-plate structures. The semi-analytical method accounts for geometric shape imperfections, and allows the use of laminated composite plate properties. Buckling eigenfunctions are used as displacement shape functions, and the associated second-order displacement fields are incorporated. The nonlinear plate equilibrium equations are expressed in terms of the assumed displacement form, and the Galerkin method is used to obtain a set of algebraic equations which are nonlinear in the modal amplitudes and the load parameter. Buckling eigensolutions for the class of structures under consideration can be obtained from the VIPASA computer code; thus, subject to some additional restrictions in loading and plate anisotropy, structural configurations which can be analyzed with regard to buckling response by VIPASA can be analyzed with regard to geometrically nonlinear response using the new method.

Stoll, Frederick↗
NASA NTRS1992

Investigation of the aeroelastic stability of the AFW wind-tunnel model using CAP-TSD

The Computational Aeroelasticity Program - Transonic Small Disturbance (CAP-TSD) code is applied to the Active Flexible Wing (AFW) wind tunnel model for prediction of the model's transonic aeroelastic behavior. A semi-span computational model is used for evaluation of symmetric motions and a full span model is used for evaluation of antisymmetric motions. Static aeroelastic solutions using CAP-TSD are computed. Dynamic flutter analyses are then performed as perturbations about the static aeroelastic deformations and presented as flutter boundaries in terms of Mach number and dynamic pressure. Flutter boundaries that take into account modal refinements, vorticity, and entropy corrections, antisymmetric motions and sensitivity to the modeling of the wing tip ballast stores are also presented and compared with experimental flutter results.

Silva, Walter A.↗
DOE OSTI2016

Advanced Potential Energy Surfaces for Molecular Simulation

Advanced potential energy surfaces are defined as theoretical models that explicitly include many-body effects that transcend the standard fixed-charge, pairwise-additive paradigm typically used in molecular simulation. However, several factors relating to their software implementation have precluded their widespread use in condensed-phase simulations: the computational cost of the theoretical models, a paucity of approximate models and algorithmic improvements that can ameliorate their cost, underdeveloped interfaces and limited dissemination in computational code bases that are widely used in the computational chemistry community, and software implementations that have not kept pace with modern high-performance computing (HPC) architectures, such as multicore CPUs and modern graphics processing units (GPUs). Here in this Feature Article we review recent progress made in these areas, including well-defined polarization approximations and new multipole electrostatic formulations, novel methods for solving the mutual polarization equations and increasing the MD time step, combining linear-scaling electronic structure methods with new QM/MM methods that account for mutual polarization between the two regions, and the greatly improved software deployment of these models and methods onto GPU and CPU hardware platforms. Finally, we have now approached an era where multipole-based polarizable force fields can be routinely used to obtain computational results comparable to state-of-the-art density functional theory while reaching sampling statistics that are acceptable when compared to that obtained from simpler fixed partial charge force fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Precompound Emission of Energetic Light Fragments in Spallation Reactions [Dissertation]

Emission of light fragments (LF) from nuclear reactions is an open question. Different reaction mechanisms contribute to their production; the relative roles of each, and how they change with incident energy, mass number of the target, and the type and emission energy of the fragments is not completely understood. None of the available models are able to accurately predict emission of LF from arbitrary reactions. However, the ability to describe production of LF (especially at energies ≳ 30 MeV) from many reactions is important for different applications, such as cosmic-ray-induced Single Event Upsets (SEUs), radiation protection, and cancer therapy with proton and heavy-ion beams, to name just a few. The Cascade-Exciton Model (CEM) version 03.03 and the Los Alamos version of the Quark-Gluon String Model (LAQGSM) version 03.03 event generators in Monte Carlo N-Particle Transport Code version 6 (MCNP6) describe quite well the spectra of fragments with sizes up to 4 He across a broad range of target masses and incident energies (up to ~ 5 GeV for CEM and up to ~ 1 TeV/A for LAQGSM). However, they do not predict the high energy tails of LF spectra heavier than 4 He well. Most LF with energies above several tens of MeV are emitted during the precompound stage of a reaction. The current versions of the CEM and LAQGSM event generators do not account for precompound emission of LF larger than 4 He. The aim of our work is to extend the precompound model in them to include such processes, leading to an increase of predictive power of LF-production in MCNP6. This entails upgrading the Modified Exciton Model currently used at the preequilibrium stage in CEM and LAQGSM. It also includes expansion and examination of the coalescence and Fermi break-up models used in the precompound stages of spallation reactions within CEM and LAQGSM. Extending our models to include emission of fragments heavier than 4 He at the precompound stage has indeed provided results that have much better agreement with experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1979

The role of energetic particle precipitation in Jovian magnetospherics. I - Secondary electrons from the ionosphere of Jupiter

Escape of photoelectrons from the Jovian ionosphere produces only a meager source of thermal plasma for the inflated centrifugally unstable magnetosphere and is unable to account for the 100 eV thermal plasma temperatures of the magnetosphere. Since the Jovian magnetosphere is well populated with highly energetic electrons, the creation of secondaries and the energy degradation of the primaries precipitating into the lower ionosphere provide additional sources of magnetospheric thermal plasma as well as of the energy for further elevating the plasma temperatures in the top side ionosphere. The efficiencies with which escaping electrons are created by precipitating electrons with energies up to several MeV are computed using energetic electron transport and thermalization codes. The more energetic incident fluxes are far less efficient in creating escaping electrons than the lower energy fluxes with only 0.001% of the secondaries escaping for a 1 MeV source versus 3% for a 1 keV source. Incident fluxes of the order of 10 to 100 per sq cm per sec per eV between 100 eV and 100 keV are required to produce 50 eV escape fluxes comparable to those generated by solar EUV.

Swartz, W. E.↗
NASA NTRS1999

Multidisciplinary Probabilistic Heat Transfer/Structural Analysis Code Developed: NESTEM

High-Speed Civil Transport (HSCT) engine combustor liners are subjected to complex thermal environments and have to endure these for thousands of hours with assured reliability. In the past, several deterministic analyses have been performed, including detailed heat transfer analyses to obtain thermal profiles and deterministic stress analyses to identify critical locations of high stresses. Actual rig tests also have been performed for segments by simulating these loading situations as closely as possible. However, it is well known that many uncertainties exist in loading (primarily thermal loads due to heat transfer), boundary conditions (end fixity unknowns), and material properties (moduli, thermal-expansion coefficients, and conductivities). The present in-house effort at the NASA Lewis Research Center is directed toward accounting for these in a formal way to assess the performance of liner components under complex and uncertain loading conditions as well as subject to other geometry- and material-related uncertainties.

Murthy, Pappu L. N.↗
DOE OSTI2017

Neutron and Gamma Correlations using CGM in MCNP 6.2.0

The general-purpose Monte Carlo radiation-transport code MCNP6{sup TM} came from the integration of MCNP5{sup TM} and MCNPX{sup TM}. The most recent release of MCNP, version 6.2.0, is scheduled for release in early 2017, which includes a variety of bug fixes and new features to the code. MCNP traditionally models nuclear reactions using Monte Carlo sampling techniques on measured and modeled cross-sectional data contained within ACE (A Compact Evaluated Nuclear Data File) libraries. In addition to reaction cross-sections, ACE libraries contain information relating to the production of secondary particles, such as secondary neutrons and photons. For example, an incident neutron on a nucleus can produce a nuclear reaction in the form (n, M{sub n}n'), where n is the incident neutron and M{sub n} is the number of secondary neutrons n' produced. Using this notation, different reactions can be described based on their neutron multiplicity, such as neutron capture (M{sub n} = 0), elastic and inelastic scattering (M{sub n} = 1), and reactions with multiple secondary neutrons (M{sub n} > 1). Accounting for secondary γ-rays emitted from the residual nuclei, this reaction can further be generalized as n, Mnn'M to describe a reaction producing M{sub n} neutrons and M -rays. These values M{sub n} and M are referred to as the neutron and multiplicities, respectively. ACE libraries contain the cross-sectional data for each reaction, and their corresponding statistically averaged multiplicities for neutrons and photons, M-bar{sub n} and M-bar{sub γ}. However, there are several limitations to reaction sampling using ACE libraries. First, the data contained within ACE libraries is limited by the accuracy of the model and/or experiment, which out of necessity can greatly simplify the true physics behind a reaction. For example, the value M{sub γ} in a (n, M{sub n}n'M{sub γ}γ) reaction is determined by the statistically averaged multiplicity M-bar{sub γ} instead of a distribution. Consequently, MCNP, unable to produce a continuous distribution, will utilizes a binary sampling around M-bar{sub γ}. While this sampling method provides a statistically average M{sub γ} over a large simulation, it is an inaccurate representation of an individual reaction and its corresponding secondary particles. Secondary neutrons are affected in a similar manner, such that different kinds of secondary neutrons can be sampled from the same interaction (eg., an inelastic neutron with an [n, 2n] neutron). Lastly, the random sampling technique for secondary particles remains completely independent of the sampled reaction, eliminating the ability to correlate secondary particles for a given reaction. To supplement these data libraries, additional codes and models describing certain nuclear reactions are employed. These codes can fill in the information missing from ACE, provide a better model to the true physics distribution, and can provide true correlated secondary particles, which has many useful applications in research and industry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

CFD Simulation of Floating Body Motion with Mooring Dynamics: Coupling MoorDyn with OpenFOAM

It is increasingly popular to use Computational Fluid Dynamics (CFD) models to study floating structures subjected to ocean waves, especially when it comes to applications of floating offshore wind turbines and Wave Energy Converters (WECs). Mooring dynamics are currently lacking in most of these applications. This paper presents a coupled simulation study of moored floating body motion by coupling two open-source libraries: a finite volume CFD toolbox, OpenFOAM, and a lumped-mass mooring model, MoorDyn. The instantaneous floating body position and velocity are passed from the body motion solver in the CFD model to the mooring model to calculate the fairlead kinematics. The mooring reaction forces, which are calculated by MoorDyn after updating the mooring system states, are then returned to the body motion solver to update the floating body motion. Both mesh deformation and overset mesh methods are used as the mesh motion solver in the CFD model to account for the floating body motion. The coupled model was validated against experimental measurements for a floating box moored with four catenary lines under the action of regular waves, which came from a preliminary test campaign for WECs. Apart from the lumped-mass mooring model, the present work also coupled a quasi-static mooring model and a finite element model with the floating body motion solver in OpenFOAM. The mooring line tensions predicted by these models were compared. The coupled model equipped with three mooring line codes may be further used to carry out survivability studies of FOWTs and WECs subject to severe sea states.

floating body↗
DOE OSTI2014

Triphoton production at hadron colliders

We present next-to-leading order predictions for the production of triphoton final states at the LHC and the Tevatron. Our results include the effect of photon fragmentation for the first time and we are able to quantify the impact of different isolation prescriptions. We find that calculations accounting for fragmentation effects at leading order, and those employing a smooth cone isolation where no fragmentation contribution is required, are in reasonable agreement with one another. However, larger differences in the predicted rates arise when higher order corrections to the fragmentation functions are included. In addition we present full analytic results for the $\gamma\gamma\gamma$ and $\gamma\gamma+$jet one-loop amplitudes. These amplitudes, which are particularly compact, may be useful to future higher-order calculations. Our results are available in the Monte Carlo code MCFM.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1981

The STD/MHD codes - Comparison of analyses with experiments at AEDC/HPDE, Reynolds Metal Co., and Hercules, Inc

Data from sufficiently well-instrumented, short-duration experiments at AEDC/HPDE, Reynolds Metal Co., and Hercules, Inc., are compared to analyses with multidimensional and time-dependent simulations with the STD/MHD computer codes. These analyses reveal detailed features of major transient events, severe loss mechanisms, and anomalous MHD behavior. In particular, these analyses predicted higher-than-design voltage drops, Hall voltage overshoots, and asymmetric voltage drops before the experimental data were available. The predictions obtained with these analyses are in excellent agreement with the experimental data and the failure predictions are consistent with the experiments. The design of large, high-interaction or advanced MHD experiments will require application of sophisticated, detailed and comprehensive computational procedures in order to account for the critical mechanisms which led to the observed behavior in these experiments.

Vetter, A. A.↗
NASA NTRSDate not supplied

Automation of the Uncertainty Quantification Process Based on Probability Boxes with DAKOTA

To date, while the use of CFD is prevalent, very few efforts have been undertaken that truly attempt to document all (or even most) of the sources of uncertainty in the simulations. Instead, the current state-of-the-art relies heavily on the experience of the CFD practitioner to estimate the uncertainty associated with their simulations through simple sensitivity studies or subject matter expertise. This practice will have to be replaced with a formal uncertainty quantification (UQ) process if CFD is to play an expanded role in the design research and engineering community, test and evaluation community, and ultimately certification for flight. This is especially true for hypersonic air-breathing propulsion systems due to the environment, scale, and duration limitations of ground test facilities. Accounting for uncertainties in a formal manner is a tedious process. Moreover, the typical CFD practitioner is not likely to be familiar with formal UQ methods. Hence, a major obstacle that has prevented the adoption of UQ methods for engineering design and development work is the lack of a tool set to automate most (if not all) of the UQ workflow. Towards this end, the SANDIA package DAKOTA (which has been developed to drive both UQ and optimization processes) will be tightly wrapped around the VULCAN-CFD code to automate the uncertainty quantification process. The automated process will be applied to an isolator turbulence model validation exercise that has previously been documented using a manual approach to the UQ process. Hence, the focus of this paper will be documenting the level to which automation can hide the UQ process details from the CFD practitioner rather than the UQ method itself.

CFD↗
DOE OSTI2016

SU-F-T-516: Effects of Inter-Fraction Organ Displacement/deformation On the Delivered Doses to the Heart, Esophagus, and Lungs in Patients Receiving Thoracic Radiotherapy

Purpose: To quantitatively assess the effects of inter-fraction changes in organ shape and location on the delivered dose distribution to the organs at risk (OAR) in lung cancer patients. Methods: This study analyzes treatment data of 10 patients, who were treated to 60Gy in 30 fractions. In each fraction a cone beam CT (CBCT) was acquired. Each CBCT was registered with the planning CT using deformable registration tools within MIM Software. The daily setup shifts were used to translate the planned dose distribution on the deformed planning CT. The structures of lungs, esophagus and heart were re-delineated by a physician on each CBCT. The doses delivered to each OAR, reflecting changes in the position and shape variations, were recomputed. Resultant daily dose volume histograms (DVHs) for OARs were computed and compared to those from the planning CT. Results: Based on the findings of two patients and 24 CBCTs analyzed so far, higher doses are delivered to the lungs and esophagus compared to the treatment plan. The dose differences per fraction between the delivered doses and those in the treatment plan are: for patient 1, lung mean dose = 5.3±1.3cGy and esophagus mean dose = 3.4±3.5cGy. For patient 2, lung mean dose = 12.0±3.9cGy and esophagus mean dose = 34.2±7.5cGy. Regarding the maximum dose to heart, the results varied (−18.9±22.0cGy for patient1 and 53.0±62.2cGy for patient2). Conclusion: The dosimetric effects of inter-fractional anatomical variations could be estimated using deformable image registration and manual organ segmentation for each CBCT. A considerable dose distribution variation between fractions was observed for the OARs. These changes are currently not taken into account while treating the patients and these may explain cases with severe side effects even when the treatment plan looks satisfactory. These results suggest the need for automated daily dose tracking and accumulation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Production of Energetic Light Fragments in the Cascade Exciton Model (CEM)

Emission of light fragments (LF) from nuclear reactions is an open question because different reaction mechanisms contribute to their production, and the relative roles of each, and how they change with incident energy, mass number of the target, and the type and emission energy of the fragments is not completely understood. None of the available models are able to accurately predict emission of LF from arbitrary reactions. However, the ability to describe production of LF (especially at energies ≳ 30 MeV) from many reactions is important for different applications, such as cosmic-ray-induced Single Event Upsets (SEUs), radiation protection, and cancer therapy with proton and heavy ion beams, to name just a few. The Cascade-Exciton Model (CEM) version 03.03 and the Los Alamos version of the Quark-Gluon String Model (LAQGSM) version 03.03 event generators in Monte Carlo N-Particle Transport Code version 6 (MCNP6) describe well the spectra of fragments with sizes up to 4 He across a broad range of target masses and incident energies (up to ∼ 5 GeV for CEM and up to ∼ 1 TeV/A for LAQGSM). However, they do not predict the high-energy tails of LF spectra heavier than 4 He well because most LF with energies above several tens of MeV are emitted during the precompound stage of a reaction, and the current versions of the CEM and LAQGSM event generators do not account for pre compound emission of LF larger than 4 He. The aim of our work is to extend the precompound model in them to include emission of light fragments up to 28 Mg, leading to an increase of predictive power of LF-production in MCNP6. This entails upgrading the Modified Exciton Model currently used at the preequilibrium stage in CEM and LAQGSM. It will also include investigation and possible expansion of the coalescence and Fermi break-up models used in the precompound stages of spallation reactions within CEM and LAQGSM. Extending our models to include emission of fragments heavier than 4 He at the precompound stage has already provided preliminary results that have much better agreement with experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1992

Further investigations of the aeroelastic behavior of the AFW wind-tunnel model using transonic small disturbance theory

The Computational Aeroelasticity Program-Transonic Small Disturbance (CAP-TSD) code, developed at LaRC, is applied to the active flexible wing wind-tunnel model for prediction of transonic aeroelastic behavior. A semi-span computational model is used for evaluation of symmetric motions, and a full-span model is used for evaluation of antisymmetric motions, and a full-span model is used for evaluation of antisymmetric motions. Static aeroelastic solutions using CAP-TSD are computed. Dynamic deformations are presented as flutter boundaries in terms of Mach number and dynamic pressure. Flutter boundaries that take into account modal refinements, vorticity and entropy corrections, antisymmetric motion, and sensitivity to the modeling of the wing tip ballast stores are also presented with experimental flutter results.

Silva, Walter A.↗
DOE OSTI2018

Airtightness of Commercial Building Envelopes: Where are we and where could we go?

Air leakage in commercial buildings accounts for about one quad of energy annually; therefore, efforts have been put in place to decrease this energy penalty. Data from blower door tests indicate that the airtightness of new building envelopes has been improving in certain niches because of requirements set by local governments and the US Army Corps of Engineers (USACE), as well as the availability of air barrier consultants. Advances in air barrier technologies and construction practices have also contributed to these improvements. Even with these promising results, efforts have been slow to revise building codes so that they establish a maximum air leakage rate requirement for new buildings and envelope retrofits. This paper summarizes current techniques used to conduct whole building air leakage tests, where these tests are mandatory, how the air leakage rate of buildings appears to be decreasing, and what stakeholders think about blower door tests. The paper also includes the nationwide source energy and financial savings potential for standalone retail buildings that were calculated using a newly developed simulation tool that more accurately estimates the effect of air leakage on energy use.

Hun, Diana E.↗
DOE OSTI2024

Subhourly Clipping Correction Model Comparison

This work will compare the Allen method and Walker method of accounting for subhourly inverter clipping power losses in hourly PV performance models. The Allen method uses a matrix lookup based on DNI clearness and clipping potential to assign a clipping correction loss at each simulation timestep. The Walker method models the PV DC power input to the inverter as a distribution over the hourly timestep and uses integration over the timestep to determine the amount of clipping that occurs within the timestep. Both these models have been recently implemented in the System Advisor Model's (SAM) open-source code, and will be applied to hourly SURFRAD datasets to analyze the subhourly clipping loss predicted by each model for different system designs and inverter loading conditions. Both models will be compared to "true" 1-minute SURFRAD data simulations to see their accuracy against more accurate 1-minute clipping correction loss predictions. This model comparisons will be investigated in more detail at the PVSC conference in Seattle, Washington June 2024.

ENERGY PLANNING, POLICY, AND ECONOMY,MATHEMATICS A↗
DOE OSTI2022

Computational investigation of incident ion angles and material erosion at rough graphite and silicon carbide divertor surfaces

We present a computational investigation of the dependence of material erosion on the incident ion angle at rough graphite and silicon carbide divertor surfaces. Ion angle distributions (IADs) for D plasmas at NSTX-U and DIII-D divertors were calculated by an equation-of-motion model that traces the ion trajectories in the sheath. Then, the effective sputtering yields and ion shadowed area fractions were calculated by a Monte Carlo micro-patterning and roughness code that applied the calculated IADs to surface topographic data that were obtained from optical confocal microscopy of rough graphite and SiC divertor surfaces from NSTX-U and DIII-D experiments. The calculations found that the effective sputtering yields, the sputtering pattern, and the shadowed area are determined by the detailed surface topology rather than the root mean square roughness R RMS , which represents deviations from a flat surface. The suppression of the effective sputtering yields for rough surfaces compared to the yield for a smooth surface was accounted for by the change of the mean local incident ion angle (LIIA) (θ'). The mean surface inclination angle distribution (SIAD) (δ) was found to be a useful parameter to estimate the LIIA from the calculated IADs. We report global empirical formulas for the mean LIIA and fraction of the area shadowed from the main ions for D plasmas for rough surfaces with B-field incident angles α = 85°–89° as a function of the mean SIAD (δ). We propose the use of the mean LIIA (θ') to estimate the sputtering yield for rough surfaces from the angular dependence of the sputtering yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗