Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “D CODES”

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 541 records · Page 30

Sensitivity of D meson azimuthal anisotropies to system size and nuclear structure

Recent experimental data at the Large Hadron Collider (LHC) confirmed that the soft sector azimuthal anisotropies in central XeXe collisions are sensitive to deformations in the wave function of the Xenon nucleus. Additionally, the CMS experiment found that D meson flow is slightly suppressed compared to other particle species when considering quark number scaling in small systems of pPb compared to PbPb. In this talk we used the D and B Mesons Modular (DAB-MOD) code coupled to Trento+v-USPhydro to calculate the R AA and v n ’s of D mesons using different energy loss models and Langevin techniques. Comparing D mesons R AA and v n in PbPb to XeXe collisions it is found that in central collisions D meson azimuthal anisotropies are sensitive to details in the nuclear structure wave function of Xenon. Additionally, we find that in mid-central collisions the smaller system size of XeXe suppresses D meson flow.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Code of Record: DOE Standard Canister (DOESC)

The United States Department of Energy (DOE) Spent Nuclear Fuel (SNF) Packaging Demonstration seeks to develop and demonstrate the designs, technology, processes, and regulatory framework for packaging DOE SNF for road-ready dry storage (RRDS); and establish the processes that will be used in a future production facility. The Packaging Demonstration will utilize the DOE Standard Canister (DOESC) for packaging select DOE-managed SNF types for interim storage, transportation and disposal as part of RRDS. Placing DOE-managed SNF into RRDS is part of the strategic framework for SNF on the Idaho National Laboratory (INL) site. To comply with DOE, INL and Nuclear Regulatory Commission (NRC) requirements, this Code of Record establishes the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code (BPVC), Section III, Division 3 as the Code that will govern DOESC and internal support structure constructioni, with certain clarifications. This Code of Record establishes a rationale for proceeding without construction certification of the DOESC (i.e., “N-stamping) and the extent to which a Registered Professional Engineer is required for the DOE Spent Fuel Packaging Demonstration. Given the (i) standard industry practice to pursue independent licensure of commercial storage casks and transportation packages by the NRC in lieu of ASME certification and (ii) guidance provided by 10 CFR 830, DOESC construction activities need not be certified (i.e., “stamped”) to the ASME BPVC. However, additional quality assurance requirements will apply as outlined in this Code of Record.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Small Engine Technology (Set) Task 8 Aeroelastic Prediction Methods

AlliedSignal Engines, in cooperation with NASA LeRC, completed an evaluation of recently developed aeroelastic computer codes using test cases from the AlliedSignal Engines fan blisk database. Test data for this task includes strain gage, light probe, performance, and steady-state pressure information obtained for conditions where synchronous or flutter vibratory conditions were found to occur. Aeroelastic codes evaluated include the quasi 3-D UNSFLO (developed at MIT and modified to include blade motion by AlliedSignal), the 2-D FREPS (developed by NASA LeRC), and the 3-D TURBO-AE (under development at NASA LeRC). Six test cases each where flutter and synchronous vibrations were found to occur were used for evaluation of UNSFLO and FREPS. In addition, one of the flutter cases was evaluated using TURBO-AE. The UNSFLO flutter evaluations were completed for 75 percent radial span and provided good agreement with the experimental test data. Synchronous evaluations were completed for UNSFLO but further enhancement needs to be added to the code before the unsteady pressures can be used to predict forced response vibratory stresses. The FREPS evaluations were hindered as the steady flow solver (SFLOW) was unable to converge to a solution for the transonic flow conditions in the fan blisk. This situation resulted in all FREPS test cases being attempted but no results were obtained during the present program. Currently, AlliedSignal is evaluating integrating FREPS with our existing steady flow solvers to bypass the SFLOW difficulties. ne TURBO-AE steady flow solution provided an excellent match with the AlliedSignal Engines calibrated DAWES 3-D viscous solver. Finally, the TURBO-AE unsteady analyses also matched experimental observations by predicting flutter for the single test case evaluated.

Eick, Chris D.↗

Three-Dimensional Particle-in-Cell Analysis of Mechanical Design Tolerances of the Scorpius Injector With Warp

Here, this article examines the mechanical design tolerances for the Scorpius injector using a 3-D particle-in-cell (PIC) code. In order to meet its objective as a high-fidelity, multipulse flash X-ray source, the accelerator has requirements for the electron beam in order to minimize the spot size of the beam focused onto a foil target to create high fluence X-rays at energies greater than 10 MeV. The requirements for alignment of the cathode, anode beam pipe, focusing solenoids, and transport solenoids specified using Trak and AMBER are independently verified using a 3-D electrostatic/magnetostatic model of the Scorpius injector, including transverse magnetic field errors that arise from mechanical alignment errors. It is predicted that the injector will meet the performance requirements for radiographic operation and has a design margin to correct for potentially unforeseen operating conditions.

42 ENGINEERING↗

Validation of a wall-interference assessment/correction procedure for airfoil tests in the Langley 0.3-meter Transonic Cryogenic Tunnel

Validation of a wall-interference assessment/correction (WIAC) procedure for a given facility requires its application to several sets of test data. It is necessary therefore to create a somewhat automated method for processing data through the various steps in the procedure. An automated procedure is also welcomed by the test engineer or eventual user in order to reduce the required effort and opportunity for error. Such a procedure has been developed for the Langley 0.3-m Transonic Cryogenic Tunnel using the TWINTN4 WIAC code. This code provides a four-wall, 2-D, transonic correction; that is, it accounts for sidewall boundary-layer effects, as well as for top and bottom wall effects on the airfoil tests. The TWINTN4 code utilizes measured pressure data at the tunnel walls and on the airfoil model; thus, classical homogeneous-wall boundary conditions are not used in the correction procedure.

Gumbert, C. R.↗

An improved version of NCOREL: A computer program for 3-D nonlinear supersonic potential flow computations

A computer code called NCOREL (for Nonconical Relaxation) has been developed to solve for supersonic full potential flows over complex geometries. The method first solves for the conical at the apex and then marches downstream in a spherical coordinate system. Implicit relaxation techniques are used to numerically solve the full potential equation at each subsequent crossflow plane. Many improvements have been made to the original code including more reliable numerics for computing wing-body flows with multiple embedded shocks, inlet flow through simulation, wake model and entropy corrections. Line relaxation or approximate factorization schemes are optionally available. Improved internal grid generation using analytic conformal mappings, supported by a simple geometric Harris wave drag input that was originally developed for panel methods and internal geometry package are some of the new features.

Siclari, Michael J.↗

Coupled 2-dimensional cascade theory for noise an d unsteady aerodynamics of blade row interaction in turbofans. Volume 2: Documentation for computer code CUP2D

A two dimensional linear aeroacoustic theory for rotor/stator interaction with unsteady coupling was derived and explored in Volume 1 of this report. Computer program CUP2D has been written in FORTRAN embodying the theoretical equations. This volume (Volume 2) describes the structure of the code, installation and running, preparation of the input file, and interpretation of the output. A sample case is provided with printouts of the input and output. The source code is included with comments linking it closely to the theoretical equations in Volume 1.

Hanson, Donald B.↗

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗

Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions

A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.

plasma↗

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗

Entropy Generation/Availability Energy Loss Analysis Inside MIT Gas Spring and "Two Space" Test Rigs

The results of the entropy generation and availability energy loss analysis under conditions of oscillating pressure and oscillating helium gas flow in two Massachusetts Institute of Technology (MIT) test rigs piston-cylinder and piston-cylinder-heat exchanger are presented. Two solution domains, the gas spring (single-space) in the piston-cylinder test rig and the gas spring + heat exchanger (two-space) in the piston-cylinder-heat exchanger test rig are of interest. Sage and CFD-ACE+ commercial numerical codes are used to obtain 1-D and 2-D computer models, respectively, of each of the two solution domains and to simulate the oscillating gas flow and heat transfer effects in these domains. Second law analysis is used to characterize the entropy generation and availability energy losses inside the two solution domains. Internal and external entropy generation and availability energy loss results predicted by Sage and CFD-ACE+ are compared. Thermodynamic loss analysis of simple systems such as the MIT test rigs are often useful to understand some important features of complex pattern forming processes in more complex systems like the Stirling engine. This study is aimed at improving numerical codes for the prediction of thermodynamic losses via the development of a loss post-processor. The incorporation of loss post-processors in Stirling engine numerical codes will facilitate Stirling engine performance optimization. Loss analysis using entropy-generation rates due to heat and fluid flow is a relatively new technique for assessing component performance. It offers a deep insight into the flow phenomena, allows a more exact calculation of losses than is possible with traditional means involving the application of loss correlations and provides an effective tool for improving component and overall system performance.

Ebiana, Asuquo B.↗

The 3-D viscous flow CFD analysis of the propeller effect on an advanced ducted propeller subsonic inlet

A time marching Navier-Stokes code called PARC3D was used to study the 3-D viscous flow associated with an advanced ducted propeller (ADP) subsonic inlet at take-off operating conditions. At a free stream Mach number of 0.2, experimental data for the inlet-with-propeller test model indicated that the airflow was attached on the cowl windward lip at an angle of attack of 25 degrees became unstable at 29 degrees, and separated at 30 degrees. An experimental study with a similar inlet and with no propeller (through-flow) indicated that flow separation occurred at an angle of attack a few degrees below the value observed when the inlet was tested with the propeller. This tends to indicate that the propeller exerts a favorable effect on the inlet performance. During the through-flow experiment a stationary blockage device was used to successfully simulate the propeller effect on the inlet flow field at angles of attack. In the present numerical study, this flow blockage was modeled via a PARC3D computational boundary condition (BC) called the screen BC. The principle formulation of this BC was based on the one-and-half dimension actuator disk theory. This screen BC was applied at the inlet propeller face station of the computational grid. Numerical results were obtained with and without the screen BC. The application of the screen BC in this numerical study provided results which are similar to the results of past experimental efforts in which either the blockage device or the propeller was used.

Iek, Chanthy↗

Controlling the size of non-axisymmetric magnetic footprints using resonant magnetic perturbations

The structure of the non-axisymmetric heat load distribution at the divertor plates is determined not only by the toroidal but also from the poloidal spectrum of non-axisymmetric field perturbations. Whether they are intrinsic, like error fields, or they are applied through 3D coils, the non-axisymmetric fields produce complex 3D edge magnetic topologies (footprints) that alter the properties of the heat and particle flux distributions on the divertor target plates. In this manuscript, a study of the impact of applied 3D field poloidal spectrum on the footprint size and structure is done for the DIII-D tokamak using the resistive MHD code M3D-C1 coupled with the field line tracing code TRIP3D. To resolve the impact of the poloidal spectrum of the magnetic perturbation, the relative phase of the two rows of in-vessel 3D coils used to produce both a n = 2 and a n = 3 perturbation is varied, where n is the toroidal harmonic of the magnetic perturbation. This shows that the largest footprint is predicted when the relative phase of the two rows is close to zero, which is also where the resonant coupling with the plasma is maximized. These results suggest that it will be challenging to decouple the footprint size from the requisite resonant coupling for RMP–ELM control. The correlation between the measured heat load and particle flux distributions at the outer divertor plates in DIII-D and the magnetic measurements is in good agreement with the predicted dependence of the magnetic footprint size on the amplitude of the resonant component of the plasma response.

3D fields↗

Application of Aeroelastic Solvers Based on Navier-Stokes Equations

A pre-release version of the Navier-Stokes solver (TURBO) was obtained from MSU. Along with Dr. Milind Bakhle of the University of Toledo, subroutines for aeroelastic analysis were developed and added to the TURBO code to develop versions 1 and 2 of the TURBO-AE code. For specified mode shape, frequency and inter-blade phase angle the code calculates the work done by the fluid on the rotor for a prescribed sinusoidal motion. Positive work on the rotor indicates instability of the rotor. The version 1 of the code calculates the work for in-phase blade motions only. In version 2 of the code, the capability for analyzing all possible inter-blade phase angles, was added. The version 2 of TURBO-AE code was validated and delivered to NASA and the industry partners of the AST project. The capabilities and the features of the code are summarized in Refs. [1] & [2]. To release the version 2 of TURBO-AE, a workshop was organized at NASA Lewis, by Dr. Srivastava and Dr. M. A. Bakhle, both of the University of Toledo, in October of 1996 for the industry partners of NASA Lewis. The workshop provided the potential users of TURBO-AE, all the relevant information required in preparing the input data, executing the code, interpreting the results and bench marking the code on their computer systems. After the code was delivered to the industry partners, user support was also provided. A new version of the Navier-Stokes solver (TURBO) was later released by MSU. This version had significant changes and upgrades over the previous version. This new version was merged with the TURBO-AE code. Also, new boundary conditions for 3-D unsteady non-reflecting boundaries, were developed by researchers from UTRC, Ref. [3]. Time was spent on understanding, familiarizing, executing and implementing the new boundary conditions into the TURBO-AE code. Work was started on the phase lagged (time-shifted) boundary condition version (version 4) of the code. This will allow the users to calculate non-zero interblade phase angles using, only one blade passage for analysis.

Keith, Theo G., Jr.↗

Analytical benchmark solution for 1-D neutron transport coupled with thermal conduction and material expansion

In this work we present fully analytical solutions for a class of finite, homogeneous, 1-D slab benchmark problems with nonlinear temperature feedback effects. The proposed class of benchmarks include multiplicative 1-D neutron transport (limited to quasistatic S{sub 2} with μ = ±1) coupled with thermal conduction, convection, Doppler broadening, and expansion effects along the length of the slab. This class of benchmark models, along with the corresponding analytical solutions, are valuable for validating multiphysics analysis frameworks that support coupled neutronics/thermal/structural calculations. Analytical solutions for the benchmarks are obtained by introducing an ansatz that the equilibrium flux and temperature distributions in the slab have the same shape. Specific values for the total microscopic cross section, σ{sub t,0}, and conductive heat transfer coefficient, h, that satisfy the assumed ansatz are then determined. A discussion of the procedure for generating benchmark models and analytical solutions is provided, along with numerical results for an example set of model parameters and thoughts on practical applications of the benchmarks for multiphysics code validation. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗