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,117 records · Page 62

Laser-Induced Fluorescence Measurements and Modeling of Nitric Oxide in Counterflow Diffusion Flames

The feasibility of making quantitative nonintrusive NO concentration ([NO]) measurements in nonpremixed flames has been assessed by obtaining laser-induced fluorescence (LIF) measurements of [NO] in counterflow diffusion flames at atmospheric and higher pressures. Comparisons at atmospheric pressure between laser-saturated fluorescence (LSF) and linear LIF measurements in four diluted ethane-air counterflow diffusion flames with strain rates from 5 to 48/s yielded excellent agreement from fuel-lean to moderately fuel-rich conditions, thus indicating the utility of a model-based quenching correction technique, which was then extended to higher pressures. Quantitative LIF measurements of [NO] in three diluted methane-air counterflow diffusion flames with strain rates from 5 to 35/s were compared with OPPDIF model predictions using the GRI (version 2.11) chemical kinetic mechanism. The comparisons revealed that the GRI mechanism underpredicts prompt-NO by 30-50% at atmospheric pressure. Based on these measurements, a modified reaction rate coefficient for the prompt-NO initiation reaction was proposed which causes the predictions to match experimental data. Temperature measurements using thin filament pyrometry (TFP) in conjunction with a new calibration method utilizing a near-adiabatic H2-air Hencken burner gave very good comparisons with model predictions in these counterflow diffusion flames. Quantitative LIF measurements of [NO] were also obtained in four methane-air counterflow partially-premixed flames with fuel-side equivalence ratios (phi(sub B)) of 1.45, 1.6, 1.8 and 2.0. The measurements were in excellent agreement with model predictions when accounting for radiative heat loss. Spatial separation between regions dominated by the prompt and thermal NO mechanisms was observed in the phi(sub B) = 1.45 flame. The modified rate coefficient proposed earlier for the prompt-NO initiation reaction improved agreement between code predictions and measurements in the region where prompt-NO dominates. Finally, LIF measurements of NO were obtained in counterflow diffusion flames at 2 to 5 atm. Comparisons between [NO] measurements and predictions show that the GRI mechanism underpredicts prompt-NO by a factor of two to three at all pressures. In general, the results indicate a need for refinement of the CH chemistry, especially the pressure-dependent CH formation and destruction reactions.

Ravikrishna, Rayavarapu V.↗

Photoionization Modeling of Oxygen K Absorption in the Interstellar Medium: The Chandra Grating Spectra of XTE J1817-330

We present detailed analyses of oxygen K absorption in the interstellar medium (ISM) using four high-resolution Chandra spectra towards the X-ray low-mass binary XTE J1817-330. The 11-25 A broadband is described with a simple absorption model that takes into account the pileup effect and results in an estimate of the hydrogen column density. The oxygen K-edge region (21-25 A) is fitted with the physical warmabs model, which is based on a photoionization model grid generated with the XSTAR code with the most up-to-date atomic database. This approach allows a benchmark of the atomic data which involves wavelength shifts of both the K lines and photoionization cross sections in order to fit the observed spectra accurately. As a result we obtain: a column density of N(sub H) = 1.38 +/- 0.01 x 10(exp 21) cm(exp -2); ionization parameter of log xi = .2.70 +/- 0.023; oxygen abundance of A(sub O) = 0.689(exp +0.015./-0.010); and ionization fractions of O I/O = 0.911, O II/O = 0.077, and O III/O = 0.012 that are in good agreement with previous studies. Since the oxygen abundance in warmabs is given relative to the solar standard of Grevesse and Sauval (1998), a rescaling with the revision by Asplund et al. (2009) yields A(sub O) = 0.952(exp +0.020/-0.013, a value close to solar that reinforces the new standard. We identify several atomic absorption lines.K-alpha , K-beta, and K-gamma in O I and O II; and K-alpha in O III, O VI, and O VII--last two probably residing in the neighborhood of the source rather than in the ISM. This is the first firm detection of oxygen K resonances with principal quantum numbers n greater than 2 associated to ISM cold absorption.

Gatuzz, E.↗

The Safe Guards Analysis Toolbox for Non-Proliferation Modeling

As nuclear technology becomes more widely available outside of nuclear weapon armed nations, the importance of preventing the illicit development of nuclear weapons has increased dramatically. A system of nuclear safeguards prevents proliferation; stopping the diversion of nuclear material from civilian facilities for the construction of nuclear weapons. These safeguards function through a combination of accountancy measurements that confirm the locations and amounts of nuclear material, and physical security at every nuclear facility. Safeguards Analysis (SGA) is a computational tool-kit developed using the MATLAB Simulink software suite that simulates a system of nuclear safeguards. The code allows users to examine current safeguards systems for weaknesses and to pinpoint the best locations for improvement. Developed as a modular framework, SGA does not directly model material flows, and relies on an outside code for this information. Instead, it represents a series of safeguards measurements that represent the values inspectors would see. Furthermore, SGA is designed to be modular in design, allowing the user to construct as complex or simple a model as needed: simulating an entire fuel-cycle comprising multiple facilities to a single measurement. SGA's use of outside codes for its material flow information makes it applicable to any model available to the user, as opposed to previous work that applied only to a single scenario. This paper discusses the initial development of SGA and the testing of these modules. The tests described in this paper demonstrate that the different pieces of SGA currently work as intended. Future work on the toolbox will develop and test SGA against a more realistic and complicated scenario using the various tools available in SGA. SGA currently has the pieces to perform analysis on safeguards systems, but needs to be tested on large, multi-MBA, systems to ensure the code is capable of handling such large models, to see the effect large models have on calculation speed, and if further work is needed to improve code efficiency. (authors)

97 MATHEMATICS AND COMPUTING↗

ESTIMATION OF DISTANCES TO STARS WITH STELLAR PARAMETERS FROM LAMOST

We present a method to estimate distances to stars with spectroscopically derived stellar parameters. The technique is a Bayesian approach with likelihood estimated via comparison of measured parameters to a grid of stellar isochrones, and returns a posterior probability density function for each star’s absolute magnitude. This technique is tailored specifically to data from the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) survey. Because LAMOST obtains roughly 3000 stellar spectra simultaneously within each ∼5° diameter “plate” that is observed, we can use the stellar parameters of the observed stars to account for the stellar luminosity function and target selection effects. This removes biasing assumptions about the underlying populations, both due to predictions of the luminosity function from stellar evolution modeling, and from Galactic models of stellar populations along each line of sight. Using calibration data of stars with known distances and stellar parameters, we show that our method recovers distances for most stars within ∼20%, but with some systematic overestimation of distances to halo giants. We apply our code to the LAMOST database, and show that the current precision of LAMOST stellar parameters permits measurements of distances with ∼40% error bars. This precision should improve as the LAMOST data pipelines continue to be refined.

79 ASTRONOMY AND ASTROPHYSICS↗

Toroidal effect on runaway vortex and avalanche growth rate

The momentum-space dynamics of runaway electrons in a slab geometry, in terms of both the geometry and topological transition of the runaway vortex when synchrotron radiative damping is taking into account, have recently been shown to play a crucial role in runaway mitigation and avoidance. In a tokamak geometry, magnetic trapping arises from parallel motion along the magnetic field that scales as 1/R in strength with R the major radius. Since the transit time for a runaway electron moving along the field is of order 10 -8 sec while the collisional time is of ~ 0.01 sec in ITER-like plasmas, a bounce-averaged formulation can drastically reduce computational cost. Here the LAPS-RFP code’s implementation of a bounce-averaged relativistic Fokker- Planck model, along with the essential physics of synchrotron radiation damping and knock-on collisions, is described. It is found that the magnetic trapping can reduce the volume of the runaway vortex as the momentum-space fluxes are strongly modified inside the trapped-region. As a result, the avalanche growth rate is reduced at off-axis locations. In addition to benchmarking with previous calculations that did not take into account radiation damping, we also clarify how synchrotron radiation damping modifies the avalanche growth rate in a tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three-dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form are iterated in time until a steady solution is reached. Evidence is given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at domain boundaries is proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects are accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. For the first, an investigation on the model behavior in case of multiple boundaries is performed. The flow in a developing S-duct is then solved in the laminar regime at Reynolds number (Re) 790 and in the turbulent regime at Re=40,000 using the Baldwin-Lomax model . The Stanitz elbow is then solved using an inviscid version of the same code at M(sub inlet)=0.4. Grid dependence and convergence rate are investigated showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re=2.5x10(exp 6) is solved with the Baldwin-Lomax and the q-omega models. Both approaches showed satisfactory agreement with experiments, although the q-omega model is slightly more accurate.

Michelassi, V.↗

Efficient construction of mock catalogs for baryon acoustic oscillation surveys

Precision measurements of the large scale structure of the Universe require large numbers of high fidelity mock catalogs to accurately assess, and account for, the presence of systematic effects. We introduce and test a scheme for generating mock catalogs rapidly using suitably derated N-body simulations. Our aim is to reproduce the large scale structure and the gross properties of dark matter halos with high accuracy, while sacrificing the details of the halo's internal structure. By adjusting global and local time-steps in an N-body code, we demonstrate that we recover halo masses to better than 0.5% and the power spectrum to better than 1% both in real and redshift space for k =1 h Mpc{sup −1}, while requiring a factor of 4 less CPU time. We also calibrate the redshift spacing of outputs required to generate simulated light cones. We find that outputs separated by Δ z =0.05 allow us to interpolate particle positions and velocities to reproduce the real and redshift space power spectra to better than 1% (out to k =1 h Mpc{sup −1}). We apply these ideas to generate a suite of simulations spanning a range of cosmologies, motivated by the Baryon Oscillation Spectroscopic Survey (BOSS) but broadly applicable to future large scale structure surveys including eBOSS and DESI. As an initial demonstration of the utility of such simulations, we calibrate the shift in the baryonic acoustic oscillation peak position as a function of galaxy bias with higher precision than has been possible so far. This paper also serves to document the simulations, which we make publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Computer simulation of nonstationary thermal fields in design and operation of northern oil and gas fields

A mathematical model, numerical algorithm and program code for simulation and long-term forecasting of changes in permafrost as a result of operation of a multiple well pad of northern oil and gas field are presented. In the model the most significant climatic and physical factors are taken into account such as solar radiation, determined by specific geographical location, heterogeneous structure of frozen soil, thermal stabilization of soil, possible insulation of the objects, seasonal fluctuations in air temperature, and freezing and thawing of the upper soil layer. Results of computing are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulation and Experimental Studies of Illumination Effects on the Current Transport of Nitridated GaAs Schottky Diode

In this paper, we studied the electrical characteristic of Schottky diodes based on gold contact on nitridated GaAs substrates. The used (100) GaAs substrate is n-type with concentration of N{sub d} = 4.9 × 10{sup 15} cm{sup –3}. Nitridation process was performed using a N2 glow discharge source (GDS) creating N atomic species. A ultra-thin film with a thickness of 2.2 nm GaN is performed on GaAs surface. In order to study the electric characteristics under illumination, we use of a He–Ne laser of 1 mW power and 632.8 nm wavelength. The current–voltage (I–V) of the Au/GaN/GaAs structures was investigated at room temperature. The saturation current I{sub S}, the series resistance R{sub S} and the mean ideality factor n are, respectively, equal to 4.46 × 10{sup –07} A, 172 Ohm, 1.4 in the dark and to 5.64 × 10{sup –07} A, 148 Ohm, 1.21 under illumination. To analyze these results, a 1D-simulation code of forward and reverse current–voltage characteristics versus the critical parameters of a nitridated GaAs Schottky diode is implemented. The algorithm is based on the solution of the system composed by Poisson’s and continuities equations. In this calculation, we take into account the existence of the GaN layer and the localization of traps states in the perturbed interface. We have considered the W-shaped and U-shaped distribution of traps states in the band gap. The effects of the doping concentration of GaAs, the traps states density, light intensities and the work function of gold ϕ{sub m} are investigated. By fitting the experimental curve, we can deduce the values of the traps states and the exact value of the work function of gold, and consequently we validate the developed model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Implementation and Validation of Photon Transport in OpenMC

OpenMC is a Monte Carlo particle transport code focused on reactor physics calculations. It stochastically simulates neutrons moving through 3D models using constructive solid geometry and continuous-energy cross sections. Recently, a photon transport capability was added to OpenMC. In reactor simulations, photon transport is needed, for example, to account for photon heating: although the majority of the energy deposited at a fission site comes from the kinetic energy of the fission fragments, a nontrivial portion can be carried away by photons. Therefore, photons produced in neutron reactions must be simulated in order to accurately model energy deposition. OpenMC can now simulate neutron-induced photon production as well as perform pure photon calculations. The four basic photon interactions with matter — coherent (Rayleigh) scattering, incoherent (Compton) scattering, the photoelectric effect, and pair production — are modeled. Secondary processes that can create new photons (i.e., atomic relaxation, electron-positron annihilation, and bremsstrahlung) are also included. The methods implemented here were validated against MCNP6 for a variety of materials and energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Integrated Approach to Fluoride High Temperature Reactor Technology and Licensing Challenges (FHR-IRP)

This document constitutes progress made on this project by the Georgia Tech led Integrated Approach to Fluoride High Temperature Reactor Technology and Licensing Challenges (FHR-IRP). This is the first draft of the project final report. Chapter 3 is focused on Tritium management. Tritium management in Fluoride-salt High-temperature Reactors (FHRs) presents a major operational concern and a significant challenge. Tritium is generated via multiple pathways including various interactions in fuel and structural materials, as well as coolant (primarily interaction with lithium and beryllium), leading to a complex evolution and spatial distribution of tritium generation rate. FHR generates orders of magnitude more tritium than in LWRs, and has potential for a significant release, via permeation through structural materials (due to high temperature), and through heat exchanger(s). Therefore, accurate estimate is necessary. For that purpose relevant tritium production pathways were identified. A methodology to obtain space-time dependent tritium generation rate was developed and implemented, employing SCALE and MCNP codes. Detailed models of the reference Advanced High Temperature Reactor (AHTR) design were developed. A multi-parametric metrics was introduced to capture relevant tritium production indicators, including the instantaneous tritium production rate, total tritium production (neglecting subsequent decay), and accumulated tritium production (accounting for the decay). The obtained results will inform tritium management, including the sizing of the components, and to gain insight on the impact of the core design parameters on tritium generation rates. Research effort on Liquid Salt Impurity Removal, Redox and Corrosion Control is discussed in Chapter 4. Section 4.1 is an introduction chapter. Based on the current available knowledge and operation experience on previous MSRs and other types of nuclear reactors, the report identifies the major impurities that are related to either material compatibility, salt properties or radioactive source terms, and the potential contamination sources of these impurities, their transport in the primary coolant loop and available cleanup methodologies for a particular group of impurity. In section 4.2, the available purification methods were discussed, including the liquid/liquid extraction using liquid Bi-Li which has been applied in MSBR and electrochemical purification which is developed based on Pyroprocessing. In this section, the new multi-stage and multi-salt-component was developed and discussed. Section 4.3 discussed the redox potential control for salt chemistry and corrosion mitigation, including the control theory and the soluble/soluble buffer control method. The section also developed method on how to analyze soluble/soluble redox reaction in a molten salt system and designed an electrochemical sensor for redox potential control. Corrosion models were also reported in this section. Section 4.4 focused on how to measure the fundamental properties in a molten salt. As an example, the section focused on La fluoride in the molten salt. The section developed methods on how analyze electrochemical measurement data to obtain the fundamental properties including the redox potential, diffusion coefficient and exchange current density. All these fundamental information is important for salt purification system development. Section 4.5 focused on measurements of corrosion products in the salts. Because of the different behaviors between Cr and La in the salt, new analysis method for Cr fluorides in the salt was developed. Some results for Fe fluoride in the salt and corrosion measurements were also reported. Section 4.6 demonstrated the separation of oxide from the salt. This section also measured the solubility of La 2 O 3 in the salt as function of temperature. Section 4.7 focused on the separation of La and Ce from the molten salt. These two elements were selected because La was found to be separated by using inert electrode, while for separation of Ce, reactive electrode must be selected. Section 4.8 demonstrated the separation iodine ions from the molten salt. Two kinds of working electrode materials (metal W and graphite) were investigated. Section 4.9 focused on the fundamental data of fission products and actinides in liquid Bi. The chapter reviewed the available data, then based a CALPHAD model, the available data were broadened to the concentration-dependent cases.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Laying the Foundation for a Predictive Theory of Electroweak Nuclear Phenomena

Reliable model-independent predictions of nuclear processes such as beta-decay and neutrinonucleus scattering are important for many fundamental science and mission-related investigations; For example, accurate descriptions of as beta-decay are required for modeling the anti-neutrino spectrum of nuclear reactor, whether for trying to detect new physics such as sterile neutrinos, or remotely identifying the purpose of that reactor. For these activities, we need accurate calculations, as well as careful and thorough analyses of theoretical error and uncertainty in the underlying models and the many-body methods and codes used in solving that model. Unfortunately, modern theoretical calculations of these processes are inadequate for this task. At the heart of this inadequacy is the interaction between nucleons and other particles known as nuclear currents. Historically, descriptions of nuclear currents took into account only the simplest of particle interactions. Recently, the development of an effective field theory known as chiral EFT has been used to derive large corrections to these particle interactions. However, the complex nature of these corrections inhibited their use in ab initio methods of predicting nuclear processes. We developed a novel method to address these complexities and applied it in a series of benchmark and demonstrative calculations using state-of-the-art quantum few-nucleon solvers. Specifically, we studied the role of two-nucleon currents in electroweak processes using light nuclei (i.e., nuclei containing fewer than five nucleons) as a workbench to compute these processes and quantify their uncertainties using Bayesian extensions of the EFT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A New CIERRA Gridded Product Climatology for LIS/OTD

The GLM-CIERRA processing generates two products: A level-2 cluster feature dataset that merges flashes that should have been a single flash following the clustering algorithm flash definition (for GLM, this is due to the LCFA thresholds); A set of level-3 gridded products including Flash Extent Density, Convective Probability, etc. The CIERRA Level-3 products are generated on a 0.1 degree grid. This is suboptimal given the spatial variations in GLM pixel size AND the size of gridpoints in the output grid. As a stepping-stone towards developing improved GLM-CIERRA grids, I am testing pixel matching techniques on the LIS/OTD data. The goal is to maintain accuracy while vastly improving computational efficiency. The CIERRA reclustering codes are also being applied to LIS / OTD to merge flashes split by the “first fit” clustering technique. The end result will be standardized grids at the nominal resolution of each instrument (5 km for LIS, 10 km for OTD) that take into account the pixel geometry for each event that comprises a given reclustered flash.

54 ENVIRONMENTAL SCIENCES↗

The use of multigrid techniques in the solution of the Elrod algorithm for a dynamically loaded journal bearing

A numerical solution to a theoretical model of vapor cavitation in a dynamically loaded journal bearing is developed, utilizing a multigrid iterative technique. The code is compared with a presently existing direct solution in terms of computational time and accuracy. The model is based on the Elrod algorithm, a control volume approach to the Reynolds equation which mimics the Jakobssen-Floberg and Olsson cavitation theory. Besides accounting for a moving cavitation boundary and conservation of mass at the boundary, it also conserves mass within the cavitated region via liquid striations. The mixed nature of the equations (elliptic in the full film zone and nonelliptic in the cavitated zone) coupled with the dynamic aspects of the problem create interesting difficulties for the present solution approach. Emphasis is placed on the methods found to eliminate solution instabilities. Excellent results are obtained for both accuracy and reduction of computational time.

Woods, Claudia M.↗

A model for the progressive failure of laminated composite structural components

Laminated continuous fiber polymeric composites are capable of sustaining substantial load induced microstructural damage prior to component failure. Because this damage eventually leads to catastrophic failure, it is essential to capture the mechanics of progressive damage in any cogent life prediction model. For the past several years the authors have been developing one solution approach to this problem. In this approach the mechanics of matrix cracking and delamination are accounted for via locally averaged internal variables which account for the kinematics of microcracking. Damage progression is predicted by using phenomenologically based damage evolution laws which depend on the load history. The result is a nonlinear and path dependent constitutive model which has previously been implemented to a finite element computer code for analysis of structural components. Using an appropriate failure model, this algorithm can be used to predict component life. In this paper the model will be utilized to demonstrate the ability to predict the load path dependence of the damage and stresses in plates subjected to fatigue loading.

Allen, D. H.↗

Production of Pions in pA-collisions

Accurate knowledge of pion production cross section in PA-collisions is of interest for astrophysics, CR physics, and space radiation studies. Meanwhile, pion production in pA-reactions is often accounted for by simple scaling of that for pp-collisions, which is not enough for many real applications. We evaluate the quality of existing parameterizations using the data and simulations with the Los Alamos version of the Quark-Gluon String Model code LAQGSM and the improved Cascade-Exciton Model code CEM2k. The LAQGSM and CEM2k models have been shown to reproduce well nuclear reactions and hadronic data in the range 0.01-800 GeV/nucleon.

Moskalenko, I. V.↗

Fluid Distribution for In-space Cryogenic Propulsion

The ultimate goal of this task is to enable the use of a single supply of cryogenic propellants for three distinct spacecraft propulsion missions: main propulsion, orbital maneuvering, and attitude control. A fluid distribution system is sought which allows large propellant flows during the first two missions while still allowing control of small propellant flows during attitude control. Existing research has identified the probable benefits of a combined thermal management/power/fluid distribution system based on the Solar Integrated Thermal Management and Power (SITMAP) cycle. Both a numerical model and an experimental model are constructed in order to predict the performance of such an integrated thermal management/propulsion system. This research task provides a numerical model and an experimental apparatus which will simulate an integrated thermal/power/fluid management system based on the SITMAP cycle, and assess its feasibility for various space missions. Various modifications are done to the cycle, such as the addition of a regeneration process that allows heat to be transferred into the working fluid prior to the solar collector, thereby reducing the collector size and weight. Fabri choking analysis was also accounted for. Finally the cycle is to be optimized for various space missions based on a mass based figure of merit, namely the System Mass Ratio (SMR). -. 1 he theoretical and experimental results from these models are be used to develop a design code (JETSIT code) which is able to provide design parameters for such a system, over a range of cooling loads, power generation, and attitude control thrust levels. The performance gains and mass savings will be compared to those of existing spacecraft systems.

Lear, William↗

The Variation of Slat Noise with Mach and Reynolds Numbers

The slat noise from the 30P30N high-lift system has been computed using a computational fluid dynamics code in conjunction with a Ffowcs Williams-Hawkings solver. By varying the Mach number from 0.13 to 0.25, the noise was found to vary roughly with the 5th power of the speed. Slight changes in the behavior with directivity angle could easily account for the different speed dependencies reported in the literature. Varying the Reynolds number from 1.4 to 2.4 million resulted in almost no differences, and primarily served to demonstrate the repeatability of the results. However, changing the underlying hybrid Reynolds-averaged-Navier-Stokes/Large-Eddy-Simulation turbulence model significantly altered the mean flow because of changes in the flap separation. However, the general trends observed in both the acoustics and near-field fluctuations were similar for both models.

Lockard, David P.↗