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 469 records · Page 26

Ablative Thermal Response Analysis Using the Finite Element Method

A review of the classic techniques used to solve ablative thermal response problems is presented. The advantages and disadvantages of both the finite element and finite difference methods are described. As a first step in developing a three dimensional finite element based ablative thermal response capability, a one dimensional computer tool has been developed. The finite element method is used to discretize the governing differential equations and Galerkin's method of weighted residuals is used to derive the element equations. A code to code comparison between the current 1-D tool and the 1-D Fully Implicit Ablation and Thermal Response Program (FIAT) has been performed.

Dec John A.↗

Kinetic and Transport Modeling for Entry Flows in Hydrogen-Helium Atmospheres

Several different thermochemical models for H2/He mixtures are used to simulate a 1-D shock in a viscous CFD code at conditions relevant for hypersonic entry into the ice and gas giants. These models range from a 6-species, 1-temperature (1-T) model, to more complex 17 and 25-species state-to-state (StS) models. Corresponding radiance profiles for each of these models are then computed using a radiation code, and results are compared to experimental data from the NASA Ames Electric Arc Shock Tube (EAST) facility. Overall, the 6-species model with a quasi-steady state (QSS) solver for atomic H and modified ionization rates, along with the 17 and 25-species StS models, are found to reproduce the magnitude of the radiance seen in the EAST experiments reasonably accurately. Two additional changes to the modeling of species diffusion fluxes are also investigated: the exclusion of a ∇W term, as well as the inclusion of Soret/ Dufour effects. It is determined that these changes to the transport modeling have a negligible impact on the flowfield and therefore the radiance predictions for the H2/He shock conditions considered in this study.

Outer planets↗

L2-Charged Particle Environment (L2-CPE)Low Energy Radiation Fluence Model

The L2 Charged Particle Environment (L2-CPE) model provides estimates of number flux and fluence for the low energy electron, proton, and alpha particle populations in the near Earth solar wind and the Earth’s distant magnetosheath and magnetotail. The model is an engineering tool for specifying radiation environments over an energy range from a few eV to a few MeV of importance to surface dose and radiation damage to thin space exposed materials and is intended for use in space system design applications. Mission fluences are obtained by simulating a spacecraft flight trajectory through time-dependent bow shock and magnetopause boundaries with dimensions and orientations driven by solar wind parameters. Monte Carlo sampling of flux environments within individual plasma regimes and/or fluence accumulated along a fight trajectory through multiple regimes is used to determine statistical variations (means and extremes) of the differential number flux and fluence environments for each of the three charged particle species. Model output is differential (in energy) number fluence and flux for the three particle species in units of particles/cm2-keV and particles/cm2-sec-keV, respectively. Users can select the energy range of interest but the current version of the code (L2-CPE Version 1.4.2d) is limited to an energy range of 1 eV to 10 MeV. Flux is integrated over angle to give flux and fluence to surfaces in the ±XGSE, ±YGSE, and/or ±ZGSE directions. Figure 1 shows the opening screen from the L2-CPE graphical user interface (GUI) with options for flux, fluence, plotting model output. The flux scene generate is not currently implemented in the GUI.

Joseph I Minow↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A new computer code was developed to solve the 2-D or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating-direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 3 is the Programmer's Reference, and describes the program structure, the FORTRAN variables stored in common blocks, and the details of each subprogram.

Towne, Charles E.↗

CFD Computation of Broadband Fan Interaction Noise

In this study, a 3-D, unsteady, Reynolds Averaged Navier Stokes CFD code coupled to an acoustic calculation is used to predict the contribution of the exit guide vanes to broadband fan noise. The configuration investigated is that corresponding to the NASA Source Diagnostic Test (SDT) 22-in fan rig. Then an acoustic model introduced by Nallasamy which is based on 2-D strip theory is used to compute the broadband rotor-stator interaction noise. One configuration from the SDT matrix is considered here: the fan speed correlating to approach, and outlet guide vane count designed for cut-off of the blade passage frequency. Thus, in the chosen configuration, there are 22 rotor blades and 54 stator blades. The stators are located 2.5 tip chords downstream of the rotor trailing edge. The RANS computations are used to obtain the spectra of the unsteady surface pressure on the exit guide vanes. This surface pressure is then integrated together with the Green's function for and infinite cylindrical duct to obtain the acoustic field. The results from this investigation validate the use of the CFD code along with the acoustic model for broadband fan noise predictions. The validation enables future investigations such as the determination of rotor tip clearance and stator solidity effects on fan rotor-stator interaction noise.

Grace, Sheryl M.↗

Maximum likelihood decoding analysis of accumulate-repeat-accumulate codes

In this paper, the performance of the repeat-accumulate codes with (ML) decoding are analyzed and compared to random codes by very tight bounds. Some simple codes are shown that perform very close to Shannon limit with maximum likelihood decoding.

LDPC codes turbo-like codes maximum likelihood dec↗

Trellis Coding of Non-coherent Multiple Symbol Full Response M-ary CPFSK with Modulation Index 1/M

This paper introduces a trellis coded modulation (TCM) scheme for non-coherent multiple full response M-ary CPFSK with modulation index 1/M. A proper branch metric for the trellis decoder is obtained by employing a simple approximation of the modified Bessel function for large signal to noise ratio (SNR). Pairwise error probability of coded sequences is evaluated by applying a linear approximation to the Rician random variable.

trellis coded modulation TCM↗

Simulations of Solar Wind Plasma Flow Around a Simple Solar Sail

In recent years, a number of solar sail missions of various designs and sizes have been proposed (e.g., Geostorm). Of importance to these missions is the interaction between the ambient solar wind plasma environment and the sail. Assuming a typical 1 AU solar wind environment of 400 km/s velocity, 3.5 cu cm density, ion temperature of approx.10 eV, electron temperature of 40 eV, and an ambient magnetic field strength of 10(exp -4) G, a first order estimate of the plasma interaction with square solar sails on the order of the sizes being considered for a Geostorm mission (50 m x 50 m and 75 m x 75 m corresponding to approx.2 and approx.3 times the Debye length in the plasma) is carried out. First, a crude current balance for the sail surface immersed in the plasma environment and in sunlight was used to estimate the surface potential of the model sails. This gave surface potentials of approx.10 V positive relative to the solar wind plasma. A 3-D, Electrostatic Particle-in-Cell (PIC) code was then used to simulate the solar wind flowing around the solar sail. It is assumed in the code that the solar wind protons can be treated as particles while the electrons follow a Boltzmann distribution. Next, the electric field and particle trajectories are solved self-consistently to give the proton flow field, the electrostatic field around the sail, and the plasma density in 3-D. The model sail was found to be surrounded by a plasma sheath within which the potential is positive compared to the ambient plasma and followed by a separate plasma wake which is negative relative to the plasma. This structure departs dramatically from a negatively charged plate such as might be found in the Earth s ionosphere on the night side where both the plate and its negative wake are contiguous. The implications of these findings are discussed as they apply to the proposed Geostorm solar sail mission.

Garrett, Henry B.↗

Analysis of an advanced ducted propeller subsonic inlet

A time marching Navier-Stokes code called PARC (PARC2D for 2-D/axisymmetric and PARC3D for 3-D flow simulations) was validated for an advanced ducted propeller (ADP) subsonic inlet. The code validation for an advanced ducted propeller (ADP) subsonic inlet. The code validation was implemented for a non-separated flow condition associated with the inlet operating at angles-of-attack of 0 and 25 degrees. The inlet test data were obtained in the 9 x 15 ft Low Speed Wind Tunnel at NASA Lewis Research Center as part of a cooperative study with Pratt and Whitney. The experimental study focused on the ADP inlet performance for take-off and approach conditions. The inlet was tested at a free stream Mach number of 0.2, at angles-of-attack between O and 35 degrees, and at a maximum propeller speed of 12,000 RPM which induced a corrected air flow rate of about 46 lb/sec based on standard day conditions. The computational grid and flow boundary conditions (BC) were based on the actual inlet geometry and the funnel flow conditions. At the propeller face, two types of BC's were applied: a mass flow BC and a fixed flow properties BC. The fixed flow properties BC was based on a combination of data obtained from the experiment and calculations using a potential flow code. Comparison of the computational results with the test data indicates that the PARC code with the propeller face fixed flow properties BC provided a better prediction of the inlet surface static pressures than the predictions when the mass flow BC was used. For an angle-of-attack of 0 degrees, the PARC2D code with the propeller face mass flow BC provided a good prediction of inlet static pressures except in the region of high pressure gradient. With the propeller face fixed flow properties BC, the PARC2D code provided a good prediction of the inlet static pressures. For an angle-of-attack of 25 degrees with the mass flow BC, the PARC3D code predicted statis pressures which deviated significantly from the test data; however, with the fixed flow properties BC, a good comparison with the test data was obtained.

Iek, Chanthy↗

Numerical simulation of viscous vortex rings

Vortex interactions and their role in turbulent flow are examined. The objectives are twofold. First, to use the existing axisymmetric code to study the annihilation process of colliding vortex rings and determine the relevance of this problem to similar 3-D phenomena. The second objective is to extend the code to three dimensions. The code under development is unique in that it can compute flows in a truly infinite domain (i.e., without periodic boundary conditions or approximations from truncating the domain). Because of this, the far field sound can be computed, and therefore, contribute to improved models of turbulence generated noise for this class of flows. Issues which can be addressed by the code include: effects of viscosity on mode selection in azimuthal breakdown of vortex rings (i.e., the Widnall instability); reconnection, the associated production of small scales, and the time scale of the process.

Stanaway, Sharon Kay↗

NAS Experiences of Porting CM Fortran Codes to HPF on IBM SP2 and SGI Power Challenge

Current Connection Machine (CM) Fortran codes developed for the CM-2 and the CM-5 represent an important class of parallel applications. Several users have employed CM Fortran codes in production mode on the CM-2 and the CM-5 for the last five to six years, constituting a heavy investment in terms of cost and time. With Thinking Machines Corporation's decision to withdraw from the hardware business and with the decommissioning of many CM-2 and CM-5 machines, the best way to protect the substantial investment in CM Fortran codes is to port the codes to High Performance Fortran (HPF) on highly parallel systems. HPF is very similar to CM Fortran and thus represents a natural transition. Conversion issues involved in porting CM Fortran codes on the CM-5 to HPF are presented. In particular, the differences between data distribution directives and the CM Fortran Utility Routines Library, as well as the equivalent functionality in the HPF Library are discussed. Several CM Fortran codes (Cannon algorithm for matrix-matrix multiplication, Linear solver Ax=b, 1-D convolution for 2-D datasets, Laplace's Equation solver, and Direct Simulation Monte Carlo (DSMC) codes have been ported to Subset HPF on the IBM SP2 and the SGI Power Challenge. Speedup ratios versus number of processors for the Linear solver and DSMC code are presented.

Saini, Subhash↗

CFD Study of Turbo-Ramjet Interactions in Hypersonic Airbreathing Propulsion System

Advanced airbreathing propulsion systems used in Mach 4-6 mission scenarios, usually involve turbo-ramjet configurations. As the engines transition from turbojet to ramjet, there is an operational envelope where both engines operate simultaneously. In the first phase of our study, an over/under nozzle configuration was analyzed. The two plumes from the turbojet and ramjet interact at the end of a common 2-D cowl, where they both reach an approximate Mach 3.0 condition and then jointly expand to Mach 3.6 at the common nozzle exit plane. For the problem analyzed, the turbojet engine operates at a higher nozzle pressure ratio than the ramjet, causes the turbojet plume overpowers the ramjet plume, deflecting it approximately 12 degrees downward and in turn the turbojet plume is deflected 6 degrees upward. In the process, shocks were formed at the deflections and a shear layer formed at the confluence of the two jets. This particular case was experimentally tested and the data were used to compare with a computational fluid dynamics (CFD) study using the PARC2D code. The CFD results were in good agreement with both static pressure distributions on the cowl separator and on nozzle walls. The thrust coefficients were also in reasonable agreement. In addition, inviscid relationships were developed around the confluence point, where the two exhaust jets meet, and these results compared favorably with the CFD results. In the second phase of our study, a 3-D CFD solution was generated to compare with the 2-D solution. The major difference between the 2-D and 3-D solutions was the interaction of the shock waves, generated by the plume interactions, on the sidewall. When a shock wave interacts with a sidewall and sidewall boundary layer, it is called a glancing shock sidewall interaction. These interactions entrain boundary layer flow down the shockline into a vortical flow pattern. The 3-D plots show the streamlines being entrained down the shockline. The pressure of the flow also decreases slightly as the sidewall is approached. Other difference between the 2-D and 3-D solutions were a lowering of the nozzle thrust coefficient value from 0.9850 (2-D) to 0.9807 (3-D), where the experimental value was 0.9790. In the third phase of our study, a different turbo-ramjet configuration was analyzed. The confluence of a supersonic turbojet and a subsonic ramjet in the turbine based combined-cycle (TBCC) propulsion system was studied by a 2-D CFD code. In the analysis, Mach 1.4 primary turbojet was mixed with the subsonic ramjet secondary flow in an ejector mode operation. Reasonable agreements were obtained with the supplied I-D TBCC solutions. For low downstream backpressure, the Fabri choke condition (Break-Point condition) was observed in the secondary flow within mixing zone. For sufficient high downstream backpressure, the Fabri choke no longer exist, the ramjet flow was reduced and the ejector flow became backpressure dependent. Highly non-uniform flow at ejector exit were observed, indicated that for smooth downstream combustion, the mixing of the two streams probably required some physical devices.

Chang, Ing↗

The Atmospheric Effects of HSCT Emissions Simulated by a 3-Wave Interactive Model

An interactive model which couples a semi-spectral dynamical model, a radiative transfer code and a two-dimensional chemistry transport model (2-D CTM), is used to assess the atmospheric effects of the High-Speed Civil Transport (HSCT) engine emissions. The residual mean meridional circulation, the zonal-mean temperature and the eddy diffusion coefficients are calculated using zonal means and three longest zonal waves of dynamical variables integrated in the semi-spectral dynamical model. They are used in the 2-D CTM to simulate the distribution of trace gases in the atmosphere. The simulated ozone is sent to the radiative transfer code to calculate the heating rates, which drive the dynamics. This radiative coupling connects the dynamical and photochemical processes and creates feedback when the atmosphere is perturbed. It is found that in most areas the ozone depletion caused by HSCT emissions calculated using the 3-wave model has the features similar to, but with significantly larger magnitude than that calculated by the AER 2-D CTM with prescribed transport parameters and temperature. The difference is mostly due to the differences in the circulation in the two models. The radiative feedback effects are investigated by comparing the ozone depletion calculated with the baseline dynamics and with the dynamics perturbed by the HSCT emissions. The feedback through changes in the residual mean meridional circulation and the eddy diffusion coefficients has moderate effects on the simulated ozone depletion. It reduced the ozone depletion by 20-30% in northern mid and high-latitudes. However, the feedback through changes in the zonal- mean temperature is negligible.

Shia, R.-L.↗

The Navy/NASA Engine Program (NNEP89): A user's manual

An engine simulation computer code called NNEP89 was written to perform 1-D steady state thermodynamic analysis of turbine engine cycles. By using a very flexible method of input, a set of standard components are connected at execution time to simulate almost any turbine engine configuration that the user could imagine. The code was used to simulate a wide range of engine cycles from turboshafts and turboprops to air turborockets and supersonic cruise variable cycle engines. Off design performance is calculated through the use of component performance maps. A chemical equilibrium model is incorporated to adequately predict chemical dissociation as well as model virtually any fuel. NNEP89 is written in standard FORTRAN77 with clear structured programming and extensive internal documentation. The standard FORTRAN77 programming allows it to be installed onto most mainframe computers and workstations without modification. The NNEP89 code was derived from the Navy/NASA Engine program (NNEP). NNEP89 provides many improvements and enhancements to the original NNEP code and incorporates features which make it easier to use for the novice user. This is a comprehensive user's guide for the NNEP89 code.

Plencner, Robert M.↗

High-Energy Space Propulsion Based on Magnetized Target Fusion

Magnetized target fusion is an approach in which a magnetized target plasma is compressed inertially by an imploding material wall. A high energy plasma liner may be used to produce the required implosion. The plasma liner is formed by the merging of a number of high momentum plasma jets converging towards the center of a sphere where two compact toroids have been introduced. Preliminary 3-D hydrodynamics modeling results using the SPHINX code of Los Alamos National Laboratory have been very encouraging and confirm earlier theoretical expectations. The concept appears ready for experimental exploration and plans for doing so are being pursued. In this talk, we explore conceptually how this innovative fusion approach could be packaged for space propulsion for interplanetary travel. We discuss the generally generic components of a baseline propulsion concept including the fusion engine, high velocity plasma accelerators, generators of compact toroids using conical theta pinches, magnetic nozzle, neutron absorption blanket, tritium reprocessing system, shock absorber, magnetohydrodynamic generator, capacitor pulsed power system, thermal management system, and micrometeorite shields.

Thio, Y. C. F.↗

Particle Acceleration, Magnetic Field Generation and Emission from Relativistic Jets and Supernova Remnants

We performed numerical simulations of particle acceleration, magnetic field generation, and emission from shocks in order to understand the observed emission from relativistic jets and supernova remnants. The investigation involves the study of collisionless shocks, where the Weibel instability is responsible for particle acceleration as well as magnetic field generation. A 3-D relativistic particle-in-cell (RPIC) code has been used to investigate the shock processes in electron-positron plasmas. The evolution of theWeibe1 instability and its associated magnetic field generation and particle acceleration are studied with two different jet velocities (0 = 2,5 - slow, fast) corresponding to either outflows in supernova remnants or relativistic jets, such as those found in AGNs and microquasars. Slow jets have intrinsically different structures in both the generated magnetic fields and the accelerated particle spectrum. In particular, the jet head has a very weak magnetic field and the ambient electrons are strongly accelerated and dragged by the jet particles. The simulation results exhibit jitter radiation from inhomogeneous magnetic fields, generated by the Weibel instability, which has different spectral properties than standard synchrotron emission in a homogeneous magnetic field.

Nishikawa, K.-I.↗