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 163 records · Page 9

Growth of zinc selenide single crystals by physical vapor transport in microgravity

The goals of this research were the optimization of growth parameters for large (20 mm diameter and length) zinc selenide single crystals with low structural defect density, and the development of a 3-D numerical model for the transport rates to be expected in physical vapor transport under a given set of thermal and geometrical boundary conditions, in order to provide guidance for an advantageous conduct of the growth experiments. In the crystal growth studies, it was decided to exclusively apply the Effusive Ampoule PVT technique (EAPVT) to the growth of ZnSe. In this technique, the accumulation of transport-limiting gaseous components at the growing crystal is suppressed by continuous effusion to vacuum of part of the vapor contents. This is achieved through calibrated leaks in one of the ground joints of the ampoule. Regarding the PVT transport rates, a 3-D spectral code was modified. After introduction of the proper boundary conditions and subroutines for the composition-dependent transport properties, the code reproduced the experimentally determined transport rates for the two cases with strongest convective flux contributions to within the experimental and numerical error.

Rosenberger, Franz↗

The Disk Wind Model of the Broad Line Regions in Active Galactic Nuclei and Cataclysmic Variables

This is the final progress report for our Astrophysics Theory Program (NRA 97-OSS12) grant NAG5-7723. We have made considerable progress on incorporating photoionization calculations with a 2.5D hydrodynamical code to model disk winds in AGNs. Following up on our simultaneous broad band monitoring campaign of the type I Seyfert galaxy NGC 5548, we have investigated the constraints imposed on models of accretion in Seyfert galaxies by their optical, UV, and X-ray spectral energy distributions (SEDs). Using results from thermal Comptonization models that relate the physical properties of the hot inner accretion flow to the thermal reprocessing that occurs in the surrounding colder thin disk, we find that we can constrain the central black hole mass, accretion rate and size scale of the hot central flow. We have applied our model to observations of Seyfert galaxies NGC 3516, NGC 7469 and NGC 5548. Our mass and accretion rate estimates for these objects roughly agree with those found using other methods.

Begelman, Mitchell↗

Convergence Acceleration and Documentation of CFD Codes for Turbomachinery Applications

The development and analysis of turbomachinery components for industrial and aerospace applications has been greatly enhanced in recent years through the advent of computational fluid dynamics (CFD) codes and techniques. Although the use of this technology has greatly reduced the time required to perform analysis and design, there still remains much room for improvement in the process. In particular, there is a steep learning curve associated with most turbomachinery CFD codes, and the computation times need to be reduced in order to facilitate their integration into standard work processes. Two turbomachinery codes have recently been developed by Dr. Daniel Dorney (MSFC) and Dr. Douglas Sondak (Boston University). These codes are entitled Aardvark (for 2-D and quasi 3-D simulations) and Phantom (for 3-D simulations). The codes utilize the General Equation Set (GES), structured grid methodology, and overset O- and H-grids. The codes have been used with success by Drs. Dorney and Sondak, as well as others within the turbomachinery community, to analyze engine components and other geometries. One of the primary objectives of this study was to establish a set of parametric input values which will enhance convergence rates for steady state simulations, as well as reduce the runtime required for unsteady cases. The goal is to reduce the turnaround time for CFD simulations, thus permitting more design parametrics to be run within a given time period. In addition, other code enhancements to reduce runtimes were investigated and implemented. The other primary goal of the study was to develop enhanced users manuals for Aardvark and Phantom. These manuals are intended to answer most questions for new users, as well as provide valuable detailed information for the experienced user. The existence of detailed user s manuals will enable new users to become proficient with the codes, as well as reducing the dependency of new users on the code authors. In order to achieve the objectives listed, the following tasks were accomplished: 1) Parametric Study Of Preconditioning Parameters And Other Code Inputs; 2) Code Modifications To Reduce Runtimes; 3) Investigation Of Compiler Options To Reduce Code Runtime; and 4) Development/Enhancement of Users Manuals for Aardvark and Phantom

Marquart, Jed E.↗

Particle Acceleration and Associated Emission from Relativistic Shocks

Five talks consist of a research program consisting of numerical simulations and theoretical development designed to provide an understanding of the emission from accelerated particles in relativistic shocks. The goal of this lecture is to discuss the particle acceleration, magnetic field generation, and radiation along with the microphysics of the shock process in a self-consistent manner. The discussion involves the collisionless shocks that produce emission from gamma-ray bursts and their afterglows, and producing emission from supernova remnants and AGN relativistic jets. Recent particle-in-cell simulation studies have shown that the Weibel (mixed mode two-stream filamentation) instability is responsible for particle (electron, positron, and ion) acceleration and magnetic field generation in relativistic collisionless shocks. 3-D RPIC code parallelized with MPI has been used to investigate the dynamics of collisionless shocks in electron-ion and electron-positron plasmas with and without initial ambient magnetic fields. In this lecture we will present brief tutorials of RPIC simulations and RMHD simulations, a brief summary of recent RPIC simulations, mechanisms of particle acceleration in relativistic shocks, and calculation of synchrotron radiation by tracing particles. We will discuss on emission from the collisionless shocks, which will be calculated during the simulation by tracing particle acceleration self-consistently in the inhomogeneous magnetic fields generated in the shocks. In particular, we will discuss the differences between standard synchrotron radiation and the jitter radiation that arises in turbulent magnetic fields.

Nishkawa, Ken-Ichi↗

Assessment of Pole Erosion in a Magnetically Shielded Hall Thruster

Numerical simulations of a 6-kW laboratory Hall thruster called H6 have been performed to quantify the erosion rate at the inner pole. The assessments have been made in two versions of the thruster, namely the unshielded (H6US) and magnetically shielded (H6MS) configurations. The simulations have been performed with the 2-D axisymmetric code Hall2De which employs a new multi-fluid ion algorithm to capture the presence of low-energy ions in the vicinity of the poles. It is found that the maximum computed erosion rate at the inner pole of the H6MS exceeds the measured rate of back-sputtered deposits by ~4.5 times. This explains only part of the surface roughening that was observed after a 150-h wear test, which covered most of the pole area exposed to the plasma. For the majority of the pole surface the computed erosion rates are found to be below the back-sputter rate and comparable to those in the H6US which exhibited little to no sputtering in previous tests. Possible explanations for the discrepancy are discussed.

plasma simulations↗

The Importance of the Cathode Plume and Its Interactions with the Ion Beam in Numerical Simulations of Hall Thrusters

Hall2De is a first-principles, 2-D axisymmetric code that solves the equations of motion for ions, electrons, and neutrals on a magnetic-field-aligned grid. The computational domain downstream of the acceleration channel exit plane is large enough to include self-consistently the cathode boundary. In this paper, we present results from numerical simulations of the H6 laboratory thruster with an internally mounted cathode, with the aim of highlighting the importance of properly accounting for the interactions between the ion beam and cathode plume. The anomalous transport of electrons across magnetic field lines in Hall2De is modelled using an anomalous collision frequency, νanom, yielding νanom approximately equal to omega ce (i.e., the electron cyclotron frequency) in the plume. We first show that restricting the anomalous collision frequency to only regions where the current density of ions is large does not alter the plasma discharge in the Hall thruster as long as the interaction between the ion beam and the cathode plume is captured properly in the computational domain. This implies that the boundary conditions must be placed sufficiently far as to not interfere with the electron transport in this region. These simulation results suggest that electron transport across magnetic field lines occurs largely inside the beam and may be driven by the interactions between beam ions and electrons. A second finding that puts in relevance the importance of including the cathode plume in numerical simulations is on the significance of accounting for the ion acoustic turbulence (IAT), now known to occur in the vicinity of the cathode exit. We have included in the Hall2De simulations a model of the IAT-driven anomalous collision frequency based on Sagdeev's model for saturation of the ion-acoustic instability. This implementation has allowed us to achieve excellent agreement with experimental measurements in the near plume obtained during the operation of the H6 thruster at nominal conditions (300V, 20A) and chamber background pressure of approximately 1.5 x 10(exp -5) Torr. In addition, the numerical results obtained with the latter approach exhibit less sensitivity to background pressure than previous attempts at explaining the features of the plasma properties in the near plume.

Lopez Ortega, Alejandro↗

Numerical Simulations of the XR-5 Hall Thruster for Life Assessment at Different Operating Conditions

NASA's Jet Propulsion Laboratory has been investigating the applicability of Aerojet Rocketdyne's XR-5 thruster, a 4.5 kW class Hall thruster, for deep-space missions. Major considerations for qualifying the XR-5 for deep-space missions are demonstration of a wide throttling envelope and a usable life capability in excess of 10,000 h. Numerical simulations with the 2-D axisymmetric code Hall2De are employed to inform the qualification process by assessing erosion rates at the thruster surfaces in a wide range of throttling conditions without the need for conducting costly endurance testing. In previous work at JPL by Jorns et al., the anomalous collision frequency distribution for 11 different throttling conditions of the XR-5 spanning 0.3-4.5 kW were identified based on probe measurements of the electron temperature in the near plume region. In this paper, we provide estimates for the erosion rates at the channel walls and pole covers for the same 11 conditions. Uncertainties in the plasma measurements and in the anomalous collision frequency distribution are addressed by determining upper and lower bounds of the erosion rates. Results suggest that erosion of the walls only occurs in the last 5% of the acceleration channel and the rate of such erosion decreases as the geometry of the thruster changes in time due to magnetic shielding. A quasi-zero-erosion state is eventually achieved in all the examined throttling conditions. Examination of the results for pole surface erosion and estimated cathode life indicates that the XR-5 propellant throughput capability will exceed 700 kg, which provides 50% margin over the usable throughput capability of 466 kg as already demonstrated in wear testing.

Lopez Ortega, Alejandro↗

Generalized Bezout's Theorem and its applications in coding theory

This paper presents a generalized Bezout theorem which can be used to determine a tighter lower bound of the number of distinct points of intersection of two or more curves for a large class of plane curves. A new approach to determine a lower bound on the minimum distance (and also the generalized Hamming weights) for algebraic-geometric codes defined from a class of plane curves is introduced, based on the generalized Bezout theorem. Examples of more efficient linear codes are constructed using the generalized Bezout theorem and the new approach. For d = 4, the linear codes constructed by the new construction are better than or equal to the known linear codes. For d greater than 5, these new codes are better than the known codes. The Klein code over GF(2(sup 3)) is also constructed.

Berg, Gene A.↗

FRANC2D: A two-dimensional crack propagation simulator. Version 2.7: User's guide

FRANC 2D (FRacture ANalysis Code, 2 Dimensions) is a menu driven, interactive finite element computer code that performs fracture mechanics analyses of 2-D structures. The code has an automatic mesh generator for triangular and quadrilateral elements. FRANC2D calculates the stress intensity factor using linear elastic fracture mechanics and evaluates crack extension using several methods that may be selected by the user. The code features a mesh refinement and adaptive mesh generation capability that is automatically developed according to the predicted crack extension direction and length. The code also has unique features that permit the analysis of layered structure with load transfer through simulated mechanical fasteners or bonded joints. The code was written for UNIX workstations with X-windows graphics and may be executed on the following computers: DEC DecStation 3000 and 5000 series, IBM RS/6000 series, Hewlitt-Packard 9000/700 series, SUN Sparc stations, and most Silicon Graphics models.

Wawrzynek, Paul↗

Prediction of ice accretion on a swept NACA 0012 airfoil and comparisons to flight test results

In the winter of 1989-90, an icing research flight project was conducted to obtain swept wing ice accretion data. Utilizing the NASA Lewis Research Center's DHC-6 DeHavilland Twin Otter aircraft, research flights were made into known icing conditions in Northeastern Ohio. The icing cloud environment and aircraft flight data were measured and recorded by an onboard data acquisition system. Upon entry into the icing environment, a 24 inch span, 15 inch chord NACA 0012 airfoil was extended from the aircraft and set to the desired sweep angle. After the growth of a well defined ice shape, the airfoil was retracted into the aircraft cabin for ice shape documentation. The ice accretions were recorded by ice tracings and photographs. Ice accretions were mostly of the glaze type and exhibited scalloping. The ice was accreted at sweep angles of 0, 30, and 45 degrees. A 3-D ice accretion prediction code was used to predict ice profiles for five selected flight test runs, which include sweep angle of zero, 30, and 45 degrees. The code's roughness input parameter was adjusted for best agreement. A simple procedure was added to the code to account for 3-D ice scalloping effects. The predicted ice profiles are compared to their respective flight test counterparts. This is the first attempt to predict ice profiles on swept wings with significant scalloped ice formations.

Reehorst, Andrew L.↗

Prediction of ice accretion on a swept NACA 0012 airfoil and comparisons to flight test results

In the winter of 1989-90, an icing research flight project was conducted to obtain swept wing ice accretion data. Utilizing the NASA Lewis Research Center's DHC-6 DeHavilland Twin Otter aircraft, research flights were made into known icing conditions in Northeastern Ohio. The icing cloud environment and aircraft flight data were measured and recorded by an onboard data acquisition system. Upon entry into the icing environment, a 24 inch span, 15 inch chord NACA 0012 airfoil was extended from the aircraft and set to the desired sweep angle. After the growth of a well defined ice shape, the airfoil was retracted into the aircraft cabin for ice shape documentation. The ice accretions were recorded by ice tracings and photographs. Ice accretions were mostly of the glaze type and exhibited scalloping. The ice was accreted at sweep angles of 0, 30, and 45 degrees. A 3-D ice accretion prediction code was used to predict ice profiles for five selected flight test runs, which include sweep angle of zero, 30, and 45 degrees. The code's roughness input parameter was adjusted for best agreement. A simple procedure was added to the code to account for 3-D ice scalloping effects. The predicted ice profiles are compared to their respective flight test counterparts. This is the first attempt to predict ice profiles on swept wings with significant scalloped ice formations.

Reehorst, Andrew L.↗

Prediction of effects of wing contour modifications on low-speed maximum lift and transonic performance for the EA-6B aircraft

Computational predictions of the effects of wing contour modifications on maximum lift and transonic performance were made and verified against low speed and transonic wind tunnel data. This effort was part of a program to improve the maneuvering capability of the EA-6B electronics countermeasures aircraft, which evolved from the A-6 attack aircraft. The predictions were based on results from three computer codes which all include viscous effects: MCARF, a 2-D subsonic panel code; TAWFIVE, a transonic full potential code; and WBPPW, a transonic small disturbance potential flow code. The modifications were previously designed with the aid of these and other codes. The wing modifications consists of contour changes to the leading edge slats and trailing edge flaps and were designed for increased maximum lift with minimum effect on transonic performance. The prediction of the effects of the modifications are presented, with emphasis on verification through comparisons with wind tunnel data from the National Transonic Facility. Attention is focused on increments in low speed maximum lift and increments in transonic lift, pitching moment, and drag resulting from the contour modifications.

Allison, Dennis O.↗

Numerical, analytical, experimental study of fluid dynamic forces in seals

NASA/Lewis Research Center is sponsoring a program for providing computer codes for analyzing and designing turbomachinery seals for future aerospace and engine systems. The program is made up of three principal components: (1) the development of advanced three dimensional (3-D) computational fluid dynamics codes, (2) the production of simpler two dimensional (2-D) industrial codes, and (3) the development of a knowledge based system (KBS) that contains an expert system to assist in seal selection and design. The first task has been to concentrate on cylindrical geometries with straight, tapered, and stepped bores. Improvements have been made by adoption of a colocated grid formulation, incorporation of higher order, time accurate schemes for transient analysis and high order discretization schemes for spatial derivatives. This report describes the mathematical formulations and presents a variety of 2-D results, including labyrinth and brush seal flows. Extensions of 3-D are presently in progress.

Shapiro, William↗

SAGE: A 2-D self-adaptive grid evolution code and its application in computational fluid dynamics

SAGE is a user-friendly, highly efficient, two-dimensional self-adaptive grid code based on Nakahashi and Deiwert's variational principles method. Grid points are redistributed into regions of high flowfield gradients while maintaining smoothness and orthogonality of the grid. Efficiency is obtained by splitting the adaption into 2 directions and applying one-sided torsion control, thus producing a 1-D elliptic system that can be solved as a set of tridiagonal equations.

Davies, Carol B.↗

Noisy quantum trees: infinite protection without correction

We study quantum networks with tree structures, in which information propagates from a root to leaves. At each node in the network, the received qubit unitarily interacts with fresh ancilla qubits, after which each qubit is sent through a noisy channel to a different node in the next level. Therefore, as the tree depth grows, there is a competition between the irreversible effect of noise and the protection against such noise achieved by the delocalization of information. In the classical setting, where each node simply copies the input bit into multiple output bits, this model has been studied as the broadcasting or reconstruction problem on trees, which has broad applications. In this work, we study the quantum version of this problem. We consider a Clifford encoder at each node that encodes the input qubit in a stabilizer code, along with a single qubit Pauli noise channel at each edge. Such noisy quantum trees describe a scenario in which one has access to a stream of fresh (low-entropy) ancilla qubits, but cannot perform error correction. Therefore, they provide a different perspective on quantum fault tolerance. Furthermore, they provide a useful model for describing the effect of noise within the encoders of concatenated codes. We prove that above certain noise thresholds, which depend on the properties of the code such as its distance, as well as the properties of the encoder, information decays exponentially with the depth of the tree. On the other hand, by studying certain efficient decoders, we prove that for codes with distance d ≥ 2 and for sufficiently small (but non-zero) noise, classical information and entanglement propagate over a noisy tree with infinite depth. Indeed, we find that this remains true even for binary trees with certain 2-qubit encoders at each node, which encodes the received qubit in the binary repetition code with distance d = 1.

Quantum information↗

Quantitative proton radiography and shadowgraphy for arbitrary intensities

Charged-particle radiography and shadowgraphy data can be directly inverted to obtain a line-integrated transverse Lorentz force or a line-integrated transverse refractive index gradient if intensity modulations due to scattering and absorption are negligible, and angular deflections are small. We develop a new direct-inversion algorithm based on plasma physics and compare it to a new Monge–Ampère code and an existing power diagram code. The measured or source intensity is represented by electrons subject to drag, and the other intensity by fixed ions. The decrease in kinetic plus electrostatic energy determines convergence. The displacement of the electrons from their initial to their equilibrium positions determines the line-integrated force or refractive index gradient. We have implemented two approaches: PIC (particle in cell) and Lagrangian fluid, in 1-D and 2-D. The PIC code works for arbitrary intensities, can work efficiently in parallel, and can make use of existing codes. The Lagrangian code requires less memory and is faster than the PIC code without massively parallel processing, but fails in 2-D for large intensity modulations. The Monge–Ampère code is by far the fastest in 2-D, without massively parallel processing, but fails for intensities with large voids, high contrast ratios and large deflections across the boundaries, and could not obtain the degree of convergence possible with the PIC code. As a result, the power diagram code was by far the slowest and most memory intensive, and failed for large peaks in the measured intensity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE ART Graphite Program Intro Status 2023

General description of the ART Graphite Program with Status for 2023. This includes list of staff members, and brief description of behavior models, licensing & code, graphite R&D, As-Fab'd Properties, Machanism and Analysis, Irradiation. Also detail FY23 Activites, with AGC Update, Oxidation Activiteis, Baseline STatus, ASME Code Development, VIC Project Update, Molten salt studies. Slide included to highlight each employee contributions for work and papers and NEUPs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Parallelized real-time physics codes for plasma control on DIII-D

A real-time safe multi-threading library was developed on the DIII-D plasma control system to optimize the real-time TORBEAM and real-time STRIDE physics codes. These physics codes are crucial for future fusion power plant operation as they provide information about electron cyclotron wave propagation and heating as well as inform about ideal plasma stability limits. The real-time TORBEAM code executed consistently in under 20 ms while the real-time STRIDE code computes in 100 ms. The multi-threading library developed in this work can be applied to other real-time physics-based codes that will be crucial for the next generation of fusion devices.

DIII-D↗