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,693 records · Page 94

Effect of Flight on the Noise from Turbulent Jets in the Generalized Acoustic Analogy

This paper presents an extension of a noise prediction method for static jets based on the Generalized Acoustic Analogy to include a non-zero ambient stream to simulate the effects of flight. Extended formulae for the propagator functions, source terms and the far-field acoustic spectrum are derived. The formulations of the equation and the boundary conditions for numerical solution of the scalar adjoint Green’s function in the presence of a non-zero ambient stream are also given. Wherever possible, the corresponding equations in the original papers for the static case and their extensions to include a flight stream are identified. The extended formulation is implemented into an existing code and used to make predictions of the effect of the flight stream on the observed far-field noise for a series of subsonic, unheated, axisymmetric jets over a range of jet-exit and free-stream velocities using Reynolds-averaged Navier-Stokes flow solutions as input. The predictions are compared with data from experiments conducted at NASA Glenn Research Center which extend a fundamental jet noise database to include a flight stream of various speeds. Descriptions of the experimental setup, data acquisition, post-processing and corrections applied to account for the external shear-layer are given. It is shown that the acoustic analogy-based method can provide reasonably good predictions of the impact of an external flight stream on the spectral characteristics of noise from unheated round jets.

Jet noise↗

Effect of Flight on the Noise from Turbulent Jets in the Generalized Acoustic Analogy

This paper presents an extension of a noise prediction method for static jets based on the Generalized Acoustic Analogy to include a non-zero ambient stream to simulate the effects of flight. Extended formulae for the propagator functions, source terms and the far-field acoustic spectrum are derived. The formulations of the equation and the boundary conditions for numerical solution of the scalar adjoint Green’s function in the presence of a non-zero ambient stream are also given. Wherever possible, the corresponding equations in the original papers for the static case and their extensions to include a flight stream are identified. The extended formulation is implemented into an existing code and used to make predictions of the effect of the flight stream on the observed far-field noise for a series of subsonic, unheated, axisymmetric jets over a range of jet-exit and free-stream velocities using Reynolds-averaged Navier-Stokes flow solutions as input. The predictions are compared with data from experiments conducted at NASA Glenn Research Center which extend a fundamental jet noise database to include a flight stream of various speeds. Descriptions of the experimental setup, data acquisition, post-processing and corrections applied to account for the external shear-layer are given. It is shown that the acoustic analogy-based method can provide reasonably good predictions of the impact of an external flight stream on the spectral characteristics of noise from unheated round jets.

Jet noise↗

Effect of Flight on the Noise from Turbulent Jets in the Generalized Acoustic Analogy

This presents an extension of a noise prediction method for static jets based on the Generalized Acoustic Analogy to include a non-zero ambient stream to simulate the effects of flight. Extended formulae for the propagator functions, source terms and the far-field acoustic spectrum are derived. The formulations of the equation and the boundary conditions for numerical solution of the scalar adjoint Green’s function in the presence of a non-zero ambient stream are also given. Wherever possible, the corresponding equations in the original papers for the static case and their extensions to include a flight stream are identified. The extended formulation is implemented into an existing code and used to make predictions of the effect of the flight stream on the observed far-field noise for a series of subsonic, unheated, axisymmetric jets over a range of jet-exit and free-stream velocities using Reynolds-averaged Navier-Stokes flow solutions as input. The predictions are compared with data from experiments conducted at NASA Glenn Research Center which extend a fundamental jet noise database to include a flight stream of various speeds. Descriptions of the experimental setup, data acquisition, post-processing and corrections applied to account for the external shear-layer are given. It is shown that the acoustic analogy-based method can provide reasonably good predictions of the impact of an external flight stream on the spectral characteristics of noise from unheated round jets.

Jet noise↗

Inertia of Ionospheric Conductance During Electron Precipitation Events in Pulsating Aurora and Polar Arcs

Using coupled SuperThermal Electron Transport (STET) and Super-thermal Proton Electron Atomic Hydrogen – tRansport in the Ionosphere and Thermosphere (SPEAH-RIT) codes, we demonstrate that temporal variability of ionospheric conductance is defined by several time scales: magnetospheric source time variations, starting time of electron precipitation, and termination of the corresponding source of magnetospheric origin. In this case, the time scales are defined by dissipation of energetic electrons and effective recombination processes. The results presented in this paper were applied in the regions of pulsating aurora and polar arcs, demonstrating the fact that ionospheric conductance requires some time to form and decay. These time delays constitute an effective “inertia” in the conductance calculation which is not accounted for in many global models which assume an instantaneous connection between precipitation and conductance. Ionospheric conductance inertia influences the temporal variation in ionospheric and magnetospheric electric fields, and as a result, impacts magnetospheric dynamics and ultimately reconfigures the electron precipitation.

Khazanov, George V↗

Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors

Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of 5 subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the 5 pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.

Dim, Odera↗

Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors

Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of five subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the five pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Mechanical behavior simulation of MEMS-based cantilever beam using COMSOL multiphysics

This paper presents the studies of mechanical behavior of MEMS cantilever beam made of poly-silicon material, using the coupling of three application modes (plane strain, electrostatics and the moving mesh) of COMSOL Multi-physics software. The cantilevers playing a key role in Micro Electro-Mechanical Systems (MEMS) devices (switches, resonators, etc) working under potential shock. This is why they require actuation under predetermined conditions, such as electrostatic force or inertial force. In this paper, we present mechanical behavior of a cantilever actuated by an electrostatic force. In addition to the simplification of calculations, the weight of the cantilever was not taken into account. Different parameters like beam displacement, electrostatics force and stress over the beam have been calculated by finite element method after having defining the geometry, the material of the cantilever model (fixed at one of ends but is free to move otherwise) and his operational space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical Study on Crossflow Printed Circuit Heat Exchanger for Advanced Small Modular Reactors

Various fluids such as water, gases (helium), molten salts (FLiNaK, FLiBe) and liquid metal (sodium) are used as a coolant of advanced small modular reactors (SMRs). The printed circuit heat exchanger (PCHE) has been adopted as the intermediate and/or secondary heat exchanger of SMR systems because this heat exchanger is compact and effective. The size and cost of PCHE can be changed by the coolant type of each SMR. In this study, the crossflow PCHE analysis code for advanced small modular reactor has been developed for the thermal design and cost estimation of the heat exchanger. The analytical solution of single pass, both unmixed fluids crossflow heat exchanger model was employed to calculate a two dimensional temperature profile of a crossflow PCHE. The analytical solution of crossflow heat exchanger was simply implemented by using built in function of the MATLAB program. The effect of fluid property uncertainty on the calculation results was evaluated. In addition, the effect of heat transfer correlations on the calculated temperature profile was analyzed by taking into account possible combinations of primary and secondary coolants in the SMR systems. Size and cost of heat exchanger were evaluated for the given temperature requirement of each SMR.

42 ENGINEERING↗

Theory of the Drift-Wave Instability at Arbitrary Collisionality

A numerically efficient framework that takes into account the effect of the Coulomb collision operator at arbitrary collisionalities is introduced. Such a model is based on the expansion of the distribution function on a Hermite-Laguerre polynomial basis to study the effects of collisions on magnetized plasma instabilities at arbitrary mean-free path. Focusing on the drift-wave instability, we show that our framework allows retrieving established collisional and collisionless limits. At the intermediate collisionalities relevant for present and future magnetic nuclear fusion devices, deviations with respect to collision operators used in state-of-the-art turbulence simulation codes show the need for retaining the full Coulomb operator in order to obtain both the correct instability growth rate and eigenmode spectrum, which, for example, may significantly impact quantitative predictions of transport. In conclusion, the exponential convergence of the spectral representation that we propose makes the representation of the velocity space dependence, including the full collision operator, more efficient than standard finite difference methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards Native Execution of Deep Learning on a Leadership-Class HPC System

Large parallel machines generally offer the best parallel performance with "native execution" that is achieved using codes developed with the optimized compilers, communication libraries, and runtimes offered on the machines. In this paper, we report and analyze performance results from native execution of deep learning on a leadership-class high-performance computing (HPC) system. Using our new code called DeepEx, we present a study of the parallel speed up and convergence rates of learning achieved with native parallel execution. In the trade-off between computational parallelism and synchronized convergence, we first focus on maximizing parallelism while still obtaining convergence. Scaling results are reported from execution on up to 15,000 GPUs using two scientific data sets from atom microscopy and protein folding applications, and also using the popular ImageNet data set. In terms of the traditional measure of parallel speed up, excellent scaling is observed up to 12,000 GPUs. Additionally, accounting for convergence rates of deep learning accuracy or error, a deep learning-specific metric called "learning speed up" is also tracked. The performance results indicate the need to evaluate parallel deep learning execution in terms of learning speed up, and point to additional directions for improved exploitation of high-end HPC systems.

Yoginath, Srikanth↗

Transonic Correction Method for Flight Dynamic Stability Analysis of Mach 0.745 Transonic Truss-Braced Wing

This paper presents a transonic correction method for obtaining dynamic stability derivatives for flight dynamic stability analysis. The method provides a transonic correction to the Theodorsen's theory of unsteady aerodynamics using FUN3D CFD solver of unsteady Reynolds-averaged Navier-Stokes equations (RANS) for a series of wing sections of the Mach 0.745 Transonic Truss-Braced Wing in pitch and plunge oscillations. Unsteady lift and pitching moment coefficients are obtained and used to develop the correction terms in the Theodorsen's theory to account for transonic aerodynamics. The unsteady lift and pitching moment derivatives with respect to the unsteady angle of attack are obtained as frequency response functions of the reduced frequency. These frequency response functions are used to compute the dynamic stability derivatives of lift and pitching moment due to the angle of attack and pitch rate and the dynamic stability derivatives for the rolling moment and yawing moment with respect to the roll rate and yaw rate. A transonic correction is applied to steady-state stability derivatives computed by VSPAERO solver using transonic small disturbance code TSFOIL coupled to an integral boundary method. A dynamic stability analysis is conducted for longitudinal and lateral-directional motions. Without transonic corrections and dynamic stability derivatives, the analysis indicates an unstable phugoid mode. The transonic correction applied to the steady-state stability derivatives computed by VSPAERO shows a stable phugoid mode. This is due to the increase of the drag stability derivatives as a result of the additional wave drag contribution in transonic flow. The effect of the transonic dynamic stability derivatives is observed to be a significant contributor to the increase in the damping values of all the flight dynamic modes of the Mach 0.745 Transonic-Truss Braced Wing.

Transonic↗

Inertia of Ionospheric Conductance During Electron Precipitation Events

Using the SuperThermal Electron Transport (STET) model coupled with the Super‐thermal Proton Electron Atomic Hydrogen—tRansport in the Ionosphere and Thermosphere (SPEAH‐RIT) codes, we demonstrate that temporal variability of ionospheric conductance is defined by several time scales: the temporal variations of the magnetospheric source, the start time of electron precipitation, and the termination of the corresponding source. In these cases, the time scales are defined by dissipation of energetic electrons and effective recombination processes. The results presented in this paper were applied in the regions of pulsating aurora and polararcs, demonstrating the fact that ionospheric conductance requires some time to form and decay. These time delays constitute an effective “inertia” in the conductance calculation which is- not accounted for in many global models which assume an instantaneous connection between precipitation and conductance. Ionospheric conductance inertia influences the temporal variation in ionospheric and magnetospheric electric fields, and as a result, impacts magnetosphere‐ionosphere (MI) dynamics on a global scale.

George V Khazanov↗

Modeling solvation effects in real-space and real-time within density functional approaches

The Polarizable Continuum Model (PCM) can be used in conjunction with Density Functional Theory (DFT) and its time-dependent extension (TDDFT) to simulate the electronic and optical properties of molecules and nanoparticles immersed in a dielectric environment, typically liquid solvents. In this contribution, we develop a methodology to account for solvation effects in real-space (and real-time) (TD)DFT calculations. The boundary elements method is used to calculate the solvent reaction potential in terms of the apparent charges that spread over the van der Waals solute surface. In a real-space representation, this potential may exhibit a Coulomb singularity at grid points that are close to the cavity surface. We propose a simple approach to regularize such singularity by using a set of spherical Gaussian functions to distribute the apparent charges. We have implemented the proposed method in the OCTOPUS code and present results for the solvation free energies and solvatochromic shifts for a representative set of organic molecules in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prediction of tissue optical properties using the Monte Carlo modeling of photon transport in turbid media and integrating spheres

Monte Carlo methods are an established technique for simulating light transport in biological tissue. Integrating spheres make experimental measurements of the reflectance and transmittance of a sample straightforward and inexpensive. This work presents an extension to existing Monte Carlo photon transport methods to simulate integrating sphere experiments. Crosstalk between spheres in dual-sphere experiments is accounted for in the method. Analytical models, previous works on Monte Carlo photon transport, and experimental measurements of a synthetic tissue phantom validate this method. We present two approaches for using this method to back-calculate the optical properties of samples. Experimental and simulation uncertainties are propagated through both methods. Both back-calculation methods find the optical properties of a sample accurately and precisely. Our model is implemented in standard Python 3 and CUDA C++ [J. Nickolls, I. Buck, M. Garland, and K. Skadron, ACM Queue 6 , 40 ( 2008 ) ] and is publicly available in Code 1.

Cook, Patrick D. (ORCID:0000000279345428)↗

Characterization of Vacuum Facility Background Gas Through Simulation and Considerations for Electric Propulsion Ground Testing

The background gas in a vacuum facility for electric propulsion ground testing is examined in detail through a series of cold flow simulations using a direct simulation Monte Carlo (DSMC) code. The focus here is on the background gas itself, its structure and characteristics, rather than assessing its interaction and impact on thruster operation. The background gas, which is often incorrectly characterized as uniform, is found to have a notable velocity within a test facility. The gas velocity has an impact on the proper measurement of pressure and the calculation of ingestion flux to a thruster. There are also considerations for best practices for tests that involve the introduction of supplemental gas flows to artificially increase the background pressure. All of these effects need to be accounted for to properly characterize the operation of electric propulsion thrusters across different ground test vacuum facilities.

Test Facilities↗

Oil Spill Thickness Determination from L-band Synthetic Aperture Radar

Efficient and effective oil spill response requires accurate information regarding slick location, extent, and thickness to direct remediation activities. Of these three, the most challenging to determine is thickness. Ideally, the needed information would be provided by remote sensing instruments, particularly those operating from space. In this study we consider the capability of L-band synthetic aperture radar (SAR) for identifying oil layer thickness for slicks in open water given the range of oil properties and environmental conditions typical in this setting. The goal is to determine layer thickness with precision equivalent to that of the Bonn Agreement Oil Appearance Code. Here we report preliminary results of an ongoing study to determine whether either relative or absolute slick thickness can be determined from L-band SAR. The study has an experimental component, which uses low noise airborne SAR data acquired over slicks to evaluate the capability to determine relative thickness, i.e., to differentiate sheen from zones of varying thicker oil within a slick. The component of the study to evaluate whether absolute thickness can be determined from SAR uses backscatter simulations within a physics- based model of electromagnetic scattering from slicked and unslicked ocean surfaces accounting for oil properties, meteorological conditions, and sea state. As part of the theoretical component of the study, we evaluate the backscatter sensitivity to oil characteristics to determine which properties most influenced the SAR measurements. These results are used to determine whether the oil-to-water ratio or the oil thickness can be quantified with reasonable accuracy based upon SAR backscatter intensities alone or whether it requires calibration to go from relative to absolute thickness. The ratio of the backscatter contrast between clean and slicked ocean surfaces is shown to be sensitive to variations within slicks that well correlate with the oil layer thickness. Determination of absolute thickness is much more challenging given the variation of oil properties as the oil weathers on the sea surface.

Holt, Benjamin↗

USE OF SMALL SPECIMENS FOR FRACTURE TOUGHNESS EVALUATION OF RPV STEELS

Small specimens are playing the key role in evaluating properties of irradiated materials. The use of small specimens provides several advantages. Typically, only a small volume of material can be irradiated in a reactor at desirable conditions in terms of temperature, neutron flux, and neutron dose. A small volume of irradiated material may also allow for easier handling of specimens. Smaller specimens reduce the amount of radioactive material, minimizing personnel exposures and waste disposal. However, use of small specimens imposes a variety of challenges as well. These challenges are associated with proper accounting for size effects and transferability of small specimen data to the real structures of interest. Any fracture toughness specimen that can be made out of standard Charpy specimen or its broken half may have exceptional utility for evaluation of reactor pressure vessels (RPVs) since it would allow one to determine and monitor directly actual fracture toughness instead of requiring indirect predictions using correlations established with impact data. The Charpy V-notch specimen is the most commonly used specimen geometry in surveillance programs and most likely to be used in advanced reactors as per ASME code.

Sokolov, Mikhail A↗

A progress update on selected AREVA NP advanced BWR methodologies

This article presents an overview of select advances in AREVA NP methods concerning 3 topics: Analytical fuel rod model, Transient and accident simulation, and ATWSi post dryout treatment. Concerning the first topic, the analytical fuel rod model relies heavily on first-principles closed-form solutions instead of ever finer nodalization and higher dimensions of standard numerical methods such as finite differences or finite elements. This analytical fuel rod model takes full account of the complications that made numerical analysis necessary, e.g. variation of pellet thermal conductivity with temperature and burnup, and the radial dependence of fission power deposition in pellets. The model also calculates thermo-mechanical parameters such as clad stress and strain, fission gas release, pellet sintering and swelling, and pellet-clad gap size and thermal resistance. The model is suited for innovative fuel designs such as doped pellets. Concerning the second topic, the BWR evaluation model, AURORA-B, is an advanced transient and accident simulator including SRELAP5, a two fluid thermal hydraulic code comprehensively reviewed by international safety authorities for BWR and PWR applications; RODEX4, a best estimate fuel performance code used to evaluate the local thermal-mechanical behavior of fuel rods during postulated transients and accidents; and MB2-K, a neutron kinetics extension of the steady state core simulator MICROBURN-B2. The analytical domains will include transients as well as control rod drop and loss of coolant accidents. Extension of RODEX4 have been made to accommodate analysis of chromium-doped fuel. Concerning the third issue,The new AREVA post-dryout models can accurately predict dryout/re-wet followed by an ultimate failure to re-wet under oscillatory ATWSi (Anticipated Transient Without SCRAM with core instability) conditions. These models, while not a boiling curve based method, can qualitatively reproduce the characteristics of the boiling curve.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗