Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “F 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 1,027 records · Page 57

HZEFRG1: An energy-dependent semiempirical nuclear fragmentation model

Methods for calculating cross sections for the breakup of high-energy heavy ions by the combined nuclear and coulomb fields of the interacting nuclei are presented. The nuclear breakup contributions are estimated with an abrasion-ablation model of heavy ion fragmentation that includes an energy-dependent, mean free path. The electromagnetic dissociation contributions arising from the interacting coulomb fields are estimated by using Weizsacker-Williams theory extended to include electric dipole and electric quadrupole contributions. The complete computer code that implements the model is included as an appendix. Extensive comparisons of cross section predictions with available experimental data are made.

Townsend, Lawrence W.↗

Evaluation of candidate working fluid formulations for the electrothermal-chemical wind tunnel

A new hypersonic test facility which can simulate conditions typical of atmospheric flight at Mach numbers up to 20 is currently under study at the NASA/LaRC Hypersonic Propulsion Branch. In the proposed research, it was suggested that a combustion augmented electrothermal wind tunnel concept may be applied to the planned hypersonic testing facility. The purpose of the current investigation is to evaluate some candidate working fluid formulations which may be used in the chemical-electrothermal wind. The efforts in the initial phase of this research were concentrated on acquiring the code used by GASL to model the electrothermal wind tunnel and testing it using the conditions of GASL simulation. The early version of the general chemical kinetics code (GCKP84) was obtained from NASA and the latest updated version of the code (LSENS) was obtained from the author Dr. Bittker. Both codes are installed on a personal computer with a 486 25 MHz processor and 16 Mbyte RAM. Since the available memory was not sufficient to debug LSENS, for the current work GCKP84 was used.

Akyurtlu, Jale F.↗

Coulomb effects in low-energy nuclear fragmentation

Early versions of the Langley nuclear fragmentation code NUCFRAG (and a publicly released version called HZEFRG1) assumed straight-line trajectories throughout the interaction. As a consequence, NUCFRAG and HZEFRG1 give unrealistic cross sections for large mass removal from the projectile and target at low energies. A correction for the distortion of the trajectory by the nuclear Coulomb fields is used to derive fragmentation cross sections. A simple energy-loss term is applied to estimate the energy downshifts that greatly alter the Coulomb trajectory at low energy. The results, which are far more realistic than prior versions of the code, should provide the data base for future transport calculations. The systematic behavior of charge-removal cross sections compares favorably with results from low-energy experiments.

Wilson, John W.↗

Experimental and analytical assessment of the thermal behavior of spiral bevel gears

An experimental and analytical study of spiral bevel gears operating in an aerospace environment has been performed. Tests were conducted within a closed loop test stand at NASA Lewis Research Center. Tests were conducted to 537 kW (720 hp) at 14,400 rpm. The effects of various operating conditions on spiral bevel gear steady state and transient temperature are presented. Also, a three-dimensional analysis of the thermal behavior was conducted using a nonlinear finite element analysis computer code. The analysis was compared to the experimental results attained in this study. The results agreed well with each other for the cases compared and were no more than 10 percent different in magnitude.

Handschuh, Robert F.↗

Analysis of automatic repeat request methods for deep-space downlinks

Automatic repeat request (ARQ) methods cannot increase the capacity of a memoryless channel. However, they can be used to decrease the complexity of the channel-coding system to achieve essentially error-free transmission and to reduce link margins when the channel characteristics are poorly predictable. This article considers ARQ methods on a power-limited channel (e.g., the deep-space channel), where it is important to minimize the total power needed to transmit the data, as opposed to a bandwidth-limited channel (e.g., terrestrial data links), where the spectral efficiency or the total required transmission time is the most relevant performance measure. In the analysis, we compare the performance of three reference concatenated coded systems used in actual deep-space missions to that obtainable by ARQ methods using the same codes, in terms of required power, time to transmit with a given number of retransmissions, and achievable probability of word error. The ultimate limits of ARQ with an arbitrary number of retransmissions are also derived.

Pollara, F.↗

Progressive Failure And Life Prediction of Ceramic and Textile Composites

An engineering approach to predict the fatigue life and progressive failure of multilayered composite and textile laminates is presented. Analytical models which account for matrix cracking, statistical fiber failures and nonlinear stress-strain behavior have been developed for both composites and textiles. The analysis method is based on a combined micromechanics, fracture mechanics and failure statistics analysis. Experimentally derived empirical coefficients are used to account for the interface of fiber and matrix, fiber strength, and fiber-matrix stiffness reductions. Similar approaches were applied to textiles using Repeating Unit Cells. In composite fatigue analysis, Walker's equation is applied for matrix fatigue cracking and Heywood's formulation is used for fiber strength fatigue degradation. The analysis has been compared with experiment with good agreement. Comparisons were made with Graphite-Epoxy, C/SiC and Nicalon/CAS composite materials. For textile materials, comparisons were made with triaxial braided and plain weave materials under biaxial or uniaxial tension. Fatigue predictions were compared with test data obtained from plain weave C/SiC materials tested at AS&M. Computer codes were developed to perform the analysis. Composite Progressive Failure Analysis for Laminates is contained in the code CPFail. Micromechanics Analysis for Textile Composites is contained in the code MicroTex. Both codes were adapted to run as subroutines for the finite element code ABAQUS and CPFail-ABAQUS and MicroTex-ABAQUS. Graphic user interface (GUI) was developed to connect CPFail and MicroTex with ABAQUS.

Xue, David Y.↗

Effect of a Nonplanar Melt-Solid Interface On Lateral Compositional Distribution During Unidirectional Solidification of a Binary Alloy with a Constant Growth Velocity V: Theory - Pt. 1

Infrared detected materials, such as Hg(1-x)Cd(x)Te, Hg(1-x)Zn(x)Te have energy gaps almost linearly proportional to their composition. Due to the wide separation of liquidus and solidus curves of their phase diagram, there are compositional segregation in both of the axial and radial directions of these crystals grown in the Bridgman system unidirectionally with constant growth rate. It is important to understand the mechanisms, which affect lateral segregation such that large radially uniform composition crystal can be produced. Following Coriel, etc's treatment, we have developed a theory to study the effect of a curved melt-solid interface shape on lateral composition distribution. The model is considered to be a cylindrical system with azimuthal symmetry and a curved melt-solid interface shape which can be expressed as a linear combination of a series of Bessell's functions. The results show that melt-solid interface shape has a dominant effect on the lateral composition distribution of these systems. For small values of beta, the solute concentration at the melt-solid interface scales linearly with interface shape with a proportional constant of the produce of beta and (1 -k), where beta = VR/D, with V as growth velocity, R as the sample radius, D as the diffusion constant and k as the distribution constant. A detailed theory will be presented. A computer code has been developed and simulations have been performed and compared with experimental results. These will be published in another paper.

Wang, Jai-Ching↗

Adventures with Implicit Methods

In this lecture we trace the historical developments of alternating direction implicit methods. In particular, we emphasize contributions originating in the Computational Fluid Dynamics Branch at Ames Research Center in the 1970's and early 1980's. Joe Steger played a seminal role in demonstrating the practicality of using an efficient, vectorized, implicit code for solving the compressible Navier-Stokes equations. Numerous discussions with Joe had a significant impact on our own research and it is a pleasure to dedicate this lecture to honor his memory.

Warming, Robert F.↗

Radiolytic Model for Chemical Composition of Europa's Atmosphere and Surface

The overall objective of the present effort is to produce models for major and selected minor components of Europa s neutral atmosphere in 1-D versus altitude and in 2-D versus altitude and longitude or latitude. A 3-D model versus all three coordinates (alt, long, lat) will be studied but development on this is at present limited by computing facilities available to the investigation team. In this first year we have focused on 1-D modeling with Co-I Valery Shematovich s Direct Simulation Monte Carlo (DSMC) code for water group species (H2O, O2, O, OH) and on 2-D with Co-I Mau Wong's version of a similar code for O2, O, CO, CO2, and Na. Surface source rates of H2O and O2 from sputtering and radiolysis are used in the 1-D model, while observations for CO2 at the Europa surface and Na detected in a neutral cloud ejected from Europa are used, along with the O2 sputtering rate, to constrain source rates in the 2-D version. With these separate approaches we are investigating a range of processes important to eventual implementation of a comprehensive 3-D atmospheric model which could be used to understand present observations and develop science requirements for future observations, e.g. from Earth and in Europa orbit. Within the second year we expect to merge the full water group calculations into the 2-D version of the DSMC code which can then be extended to 3-D, pending availability of computing resources. Another important goal in the second year would be the inclusion of sulk and its more volatile oxides (SO, SO2).

Cooper, John F.↗

Digital Equivalent Data System for XRF Labeling of Objects

A digital equivalent data system (DEDS) is a system for identifying objects by means of the x-ray fluorescence (XRF) spectra of labeling elements that are encased in or deposited on the objects. As such, a DEDS is a revolutionary new major subsystem of an XRF system. A DEDS embodies the means for converting the spectral data output of an XRF scanner to an ASCII alphanumeric or barcode label that can be used to identify (or verify the assumed or apparent identity of) an XRF-scanned object. A typical XRF spectrum of interest contains peaks at photon energies associated with specific elements on the Periodic Table (see figure). The height of each spectral peak above the local background spectral intensity is proportional to the relative abundance of the corresponding element. Alphanumeric values are assigned to the relative abundances of the elements. Hence, if an object contained labeling elements in suitably chosen proportions, an alphanumeric representation of the object could be extracted from its XRF spectrum. The mixture of labeling elements and for reading the XRF spectrum would be compatible with one of the labeling conventions now used for bar codes and binary matrix patterns (essentially, two-dimensional bar codes that resemble checkerboards). A further benefit of such compatibility is that it would enable the conversion of the XRF spectral output to a bar or matrix-coded label, if needed. In short, a process previously used only for material composition analysis has been reapplied to the world of identification. This new level of verification is now being used for "authentication."

Schramm, Harry F.↗

Covariance generation for the prompt neutron multiplicity of 239 Pu including the (n, γ f) process

Fission cross section of 239 Pu can be seen as a sum of the “immediate" fission and “two-step" (n,γf) reactions. In the Resolved Resonance Range of the reaction cross sections, the contribution of the (n,γf) process has an impact on the determination of the partial widths magnitude involved in the Reich-Moore approximation of the R-matrix theory. The present work aims to investigate this impact by using the CONRAD code and the partial width Γ γf for the (n,γf) reaction calculated by Lynn et al. A special attention will be paid to the covariance matrix obtained on ν p .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Experimental observation of the geodesic acoustic frequency limit for the NBI-driven Alfvén eigenmodes in TJ-II

We study Alfvén Eigenmodes (AEs) in the TJ-II heliac in hydrogen plasmas heated by hydrogen cofield NBI. Taking advantage of the unique TJ-II flexibility in varying plasma current, we have observed strong variation of the AE frequency from f AE ~ 30 to ~ 220 kHz for selected modes. An advanced heavy ion beam probe diagnostic determines the spatial location and internal amplitudes of the modes. The modes satisfy a local AE dispersion relation including the geodesic acoustic frequency that represents the lowest frequency of the mode. Linear MHD-modeling with STELLGAP and FAR3D codes shows that the calculated temporal evolution of the mode frequency reproduces the observed maxima and minima at the same time intervals with a similar frequency range and the radial profile peaks near the outer edge of the observed one.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Initial calculations for source term of Molten Salt Reactors

This paper provides an overview of the current MSR design space and lists unique features of the various designs under consideration. Some general considerations for source terms calculation for Molten Salt Reactors (MSRs) are explained. Applicability and limitations of terminology currently defined for legacy light water reactor (LWR) systems are discussed in the view of MSRs and the need for updated terminology is discussed. Calculations carried out for the Molten Salt Reactor Experiment (MSRE) are discussed with a qualitative comparison to the designs presented. The nature of the fission products (FPs) and actinides for Low enriched uranium, thorium and fast U/Pu fuel cycles employed in representative molten salt reactor systems are discussed. Computational results are obtained from a code (Serpent 2) with online reprocessing. Divergence in source terms when fission product bubbling is demonstrated. The source release for each molten salt reactor during postulated accidents is also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Classical-quantum simulation of non-equilibrium Marshak waves

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Benchmark Exercise for the Control Rod Swelling Evaluation

The VTR core has six reactivity control assemblies and three safety assemblies. The control assemblies or primary control rods are adjusted during the normal operation to balance the core reactivity and to control the reactor power. A typical control assembly radial layout is presented in Figure 1. The figure shows the swelled absorber (B 4 C) rod. Initially, helium gas fills the gap between the pin and the cladding before irradiation swelling takes place. For VTR, HT9 steel was selected as the cladding and duct material. The main neutron absorbing material used in the VTR is B 4 C. When residing in the core, the neutronics, thermophysical, and mechanical properties of the materials used in a control assembly will degrade due to accumulated neutron damage. Material degradation limits how long a control assembly can reside in the core. Many phenomena affect the control assembly lifetime, such as the loss of reactivity worth due to B 4 C depletion, the mechanical interaction of the absorber rod and the cladding due to B 4 C swelling, the helium gas buildup in the pin due to B-10 capture, etc. B 4 C swelling, which causes closure of the gap between the absorber rod and the cladding, is usually considered as the main limiting factor from past experience. An initial study was conducted at PNNL to evaluate the irradiation behavior of a VTR control assembly. The evaluation was performed using the CNRD2 code that was initially developed for the FFTF. The study also included an assessment of the VTR control assembly and focused on a 61-pin control assembly design, which is different from that used (37-pin design) in the core design study. The study conducted by PNNL was reviewed independently by ANL. A Python script referred to as the Control Assembly Evaluation Script (CAES) was developed for the independent review and additional assessment of 37-pin control assembly design. The script has focused on the assessment of the absorber rod swelling for its importance in determining the control assembly lifetime. CAES uses geometry, neutronics, materials data as input to predict the swelling of the absorber rod during its residence in the reactor core. The results from CAES showed some non-negligible differences against the PNNL results. Some of the differences can be attributed to the different interpretation of the control rod assembly dimensions. To resolve this issue, a benchmark exercise was proposed. The benchmark specification was developed by PNNL. The benchmark exercise was performed independently at PNNL and ANL using different codes/scripts (CRND2 and CAES). This memo documents the results calculated using the different codes. However, this report is limited to presenting the results obtained. Further investigation of the cause of the observed difference will be performed as part of future activities, pending continuation of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Aerodynamic Characteristics of a Supersonic Fighter Aircraft Model at Mach 0.40 to 2.47

The aerodynamic characteristics of an advanced twin-engine fighter aircraft designed for supersonic cruise have been studied in the Langley 16-Foot Transonic Tunnel and the Lewis 10- by 10-Foot Supersonic Tunnel. The objective of this investigation was to establish an aerodynamic data base for the configuration with flow-through nacelles and representative inlets. The use of a canard for trim and the effects of fairing over the inlets were assessed. Comparisons between experimental and theoretical results were also made. The theoretical results were determined by using a potential vortex lift code for subsonic speeds and a linear aerodynamic code for supersonic speeds. This investigation was conducted at Mach numbers from 0.40 to 2.47, at angles of attack from 0 deg to about 20 deg, and at inlet capture ratios of about 0.5 to 1.4.

Capone, F. J.↗

An experimental and computational study of turbulent boundary layers in oscillating flows

An experimental-computational study of the behavior of turbulent boundary layers for oscillating air flows over a plane surface with a small favorable mean pressure gradient is described. Experimental data were obtained for boundary layers generated on the test section wall of a new type of facility that produces a mean flow velocity with a superposed sinusoidal component over a wide frequency range. Experimental results for unsteady velocity profiles and turbulence intensity obtained from hot wire measurements at a fixed stream-wise position on the surface are presented for frequencies ranging from 5 to 21 Hz. Numerical predictions were obtained using an unsteady boundary layer computational code and the Cebeci-Smith and Glushko turbulence models. Predicted quantities related to unsteady velocity profiles exhibit fair agreement with experiment when the Cebeci-Smith model is used.

Cook, W. J.↗

Relaxation algorithms for the Euler equations

The alternating direction implicit central-difference scheme algorithms may be replaced by successive-line-relaxation (SLR) procedures that are stable in three dimensions. Several Beamand Warming-type codes are presently modified on the basis of SLR solution procedures in order to compute the flows over a cylinder, a NACA-0012 airfoil, and a shock wave reflecting off a plate. Since the codes are written in a 'delta form', the implementation of relaxation procedures that rely only on the corrections is relatively easy; the steady state residual computation and any artificial viscosity terms may be left unchanged.

Wornom, S. F.↗