Implementation of Monte Carlo Transport in the General Relativistic SpEC Code
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
The objective of this report is to summarize the deterministic and probabilistic structural evaluation results of two structures made with advanced ceramic composites (CMC): internally pressurized tube and uniformly loaded flange. The deterministic structural evaluation includes stress, displacement, and buckling analyses. It is carried out using the finite element code MHOST, developed for the 3-D inelastic analysis of structures that are made with advanced materials. The probabilistic evaluation is performed using the integrated probabilistic assessment of composite structures computer code IPACS. The affects of uncertainties in primitive variables related to the material, fabrication process, and loadings on the material property and structural response behavior are quantified. The primitive variables considered are: thermo-mechanical properties of fiber and matrix, fiber and void volume ratios, use temperature, and pressure. The probabilistic structural analysis and probabilistic strength results are used by IPACS to perform reliability and risk evaluation of the two structures. The results will show that the sensitivity information obtained for the two composite structures from the computational simulation can be used to alter the design process to meet desired service requirements. In addition to detailed probabilistic analysis of the two structures, the following were performed specifically on the CMC tube: (1) predicted the failure load and the buckling load, (2) performed coupled non-deterministic multi-disciplinary structural analysis, and (3) demonstrated that probabilistic sensitivities can be used to select a reduced set of design variables for optimization.
To assess engine performance during the testing of Space Shuttle Main Engines (SSMEs), the design of an optimal altitude diffuser is studied for future Space Transportation Systems (STS). For other Space Transportation Systems, rocket propellant using kerosene is also studied. Methane and dodecane have similar reaction schemes as kerosene, and are used to simulate kerosene combustion processes at various temperatures. The equations for the methane combustion mechanism at high temperature are given, and engine combustion is simulated on the General Aerodynamic Simulation Program (GASP). The successful design of an altitude diffuser depends on the study of a sub-scaled diffuser model tested through two-dimensional (2-D) flow-techniques. Subroutines given calculate the static temperature and pressure at each Mach number within the diffuser flow. Implementing these subroutines into program code for the properties of 2-D compressible fluid flow determines all fluid characteristics, and will be used in the development of an optimal diffuser design.
A reliable and efficient iterative method has been developed for designing wing section contours corresponding to a prescribed subcritical pressure distribution. The design process is initialized by using MCAERO (MCAIR 3-D Subsonic Potential Flow Analysis Code) to analyze a baseline configuration. A second program DMCAERO is then used to calculate a matrix containing the partial derivative of potential at each control point with respect to each unknown geometry parameter by applying a first-order expansion to the baseline equations in MCAERO. This matrix is calculated only once but is used in each iteration cycle to calculate the geometry perturbation and to analyze the perturbed geometry. The potential on the new geometry is calculated by linear extrapolation from the baseline solution. This extrapolated potential is converted to velocity by numerical differentiation, and velocity is converted to pressure by using Bernoulli's equation. There is an interactive graphics option which allows the user to graphically display the results of the design process and to interactively change either the geometry or the prescribed pressure distribution.
Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The jitter radiation (Medvedev 2000) from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.
Shock wave acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The 'jitter' radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.
Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.
The Griesmer bound is a classical technique (developed in 1960) for estimating the minimum length n required for a binary linear code with a given dimension k and minimum distance d. In this article, a unified derivation of the Griesmer bound and two new variations on it are presented. The first variation deals with linear codes which contain the all-ones vector; such codes are quite common and are useful in practice because of their 'transparent' properties. The second variation deals with codes that are constrained to contain a word of weight greater than or equal to M. In both cases these constraints (the all-ones word or a word of high weight) can increase the minimum length of a code with given k and d.
A Monte Carlo code simulating neutron transport in infinite cones of water and water-equivalent hydrogen was prepared for an IBM 704 computer. The code was essentially a modification of the point-source, infinite-medium code used in NASA TN D-850. Studies were made of differential neutron number spectra and associated buildup factors for infinite cones having apex half-angles of 15 degrees, 30 degrees, 45 degrees, and 60 degrees. The buildup factors obtained were compared with those for the appropriate infinite medium, which allowed an examination of the effect of solid angle subtended by material on the transport of 6-Mev source neutrons emanating from the cone apex. The variation of number buildup factor with distance for the various cones shows that neutron scattering out of the cones is predominant in the first 30 t o 40 centimeters of material, and that transport beyond this distance is of a similar nature in all the cones.
An experimental and theoretical study was undertaken of the separation delay and dramatic boundary-layer thinning that can occur in vortex-generator installations. Wind tunnel measurements of the dynamic-pressure profile downstream of a vortex generator were found to compare under certain conditions with that downstream of a suction slit, while water-tunnel visualization studies of vortex-generator height and geometry suggested optimum configurations, and only a minor effect of base porosity. A series of progressively more complex inviscid flow models was developed to be applied to a 3-D integral boundary-layer code. This code predicted layer thinning downstream of the suction site of the vortex models, and other observed features. Thin-layer Navier-Stokes equations are now being used with the ultimate goal of clarifying the physical processes involved in vortex generator performance and developing calculational procedures capable of predicting it.
NASA LeRC has identified the Rich burn/Quick mix/Lean burn (RQL) combustor as a potential gas turbine combustor concept to reduce NOx emissions in High Speed Civil Transport (HSCT) aircraft. To demonstrate reduced NOx levels, NASA LeRC soon will test a flametube version of an RQL combustor. The critical technology needed for the RQL combustor is a method of quickly mixing combustion air with rich burn gases. Two concepts were proposed to enhance jet mixing in a circular cross-section: the Asymmetric Jet Penetration (AJP) concept; and the Lobed Mixer (LM) concept. In Phase 1, two preliminary configurations of the AJP concept were compared with a conventional 12-jet radial-inflow slot design. The configurations were screened using an advanced 3-D Computational Fluid Dynamics (CFD) code named REFLEQS. Both non-reacting and reacting analyses were performed. For an objective comparison, the conventional design was optimized by parametric variation of the jet-to-mainstream momentum flux (J) ratio. The optimum J was then employed in the AJP simulations. Results showed that the three-jet AJP configuration was superior in overall mixedness compared to the conventional design. However, in regards to NOx emissions, the AJP configuration was inferior. The higher emission level for AJP was caused by a single hot spot located in the wake of the central jet as it entered the combustor. Ways of maintaining good mixedness while eliminating the hot spot were identified for Phase 2 study. Overall, Phase 1 showed the viability of using CFD analyses to evaluate quick-mix concepts. A high probability exists that advancing mixing concepts will reduce NOx emissions in RQL combustors, and should be explored in Phase 2, by parallel numerical and experimental work.
Radiation dosimetry aboard LDEF, as well as post-flight measurements of the radioactivity induced in numerous LDEF spacecraft components, provide unique data for evaluating the accuracy of current models for predicting both the space radiation environments (trapped proton intensity, spectra, and directionality; cosmic ray fluence) and the radiation environments induced in spacecraft components (fluence, energy spectra, secondary particle, linear energy transfer, etc.). By determining the accuracy of such models using LDEF data, and with model updates where required, improved radiation environment predictions can be made for future missions, which in turn allows improved predictions for specific radiation effects for future spacecraft components (single event upsets of microelectronics, radiation damage to focal plane arrays, noise in sensitive instrumentation, etc.). Herein, the status and results from radiation model predictions and comparisons with LDEF data is given. The calculations are made using radiation transport codes coupled with a 3-D geometry/mass model of LDEF, together with current models of the space radiation environment.
Computational fluid dynamics (CFD) analysis is providing verification of Space Shuttle flight performance details and is being applied to Space Shuttle Main Engine Multiple plume interaction flow field definition. Advancements in real-gas CFD methodology that are described have allowed definition of exhaust plume flow details at Mach 3.5 and 107,000 ft. The specific objective includes the estimate of flow properties at oblique shocks between plumes and plume recirculation into the Space Shuttle Orbiter base so that base heating and base pressure can be modeled accurately. The approach utilizes the Rockwell USA Real Gas 3-D Navier-Stokes (USARG3D) Code for the analysis. The code has multi-zonal capability to detail the geometry of the plumes based region and utilizes finite-rate chemistry to compute the plume expansion angle and relevant flow properties at altitude correctly. Through an improved definition of the base recirculation flow properties, heating, and aerodynamic design environments of the Space Shuttle Vehicle can be further updated.
The flow inside a model exhaust configuration was studied using both experimental and computational techniques. The hardware was tested at the NASA Lewis Research Center's Powered Lift Facility at tailpipe total pressure to ambient static pressure ratios ranging from 1.0 to 5.0. The flow simulations were obtained using the two 3-D Navier-Stokes CFD codes run on the Lewis Cray Y-MP computer. Both codes produced oscillatory solutions due to the inflow boundary condition reflecting acoustic waves. The CFD solutions correctly predicted the flow separation along the inside elbow of the takeoff and also along the walls of the ventral duct. Mass flow rates were overpredicted due to underprediction of the turbulent energy dissipation and subsequent total pressure loss.
As axial and radial turbine designs have been pushed to their aerothermodynamic and mechanical limits, the mixed-flow turbine (MFT) concept has been projected to offer performance and durability improvements, especially when ceramic materials are considered. The objective of this NASA/U.S. Army sponsored mixed-flow turbine (AMFT) program was to determine the level of performance attainable with MFT technology within the mechanical constraints of 1997 projected ceramic material properties. The MFT geometry is similar to a radial turbine, exhibiting a large radius change from inlet to exit, but differing in that the inlet flowpath is not purely radial, nor axial, but mixed; it is the inlet geometry that gives rise to the name 'mixed-flow'. The 'mixed' orientation of the turbine inlet offers several advantages over radial designs by allowing a nonzero inlet blade angle yet maintaining radial-element blades. The oblique inlet not only improves the particle-impact survivability of the design, but improves the aerodynamic performance by reducing the incidence at the blade inlet. The difficulty, however, of using mixed-flow geometry lies in the scarcity of detailed data and documented design experience. This paper reports the design of a MFT stage designed with the intent to maximize aerodynamic performance by optimizing design parameters such as stage reaction, rotor incidence, flowpath shape, blade shape, vane geometry, and airfoil counts using 2-D, 3-D inviscid, and 3-D viscous computational fluid dynamics code. The aerodynamic optimization was accomplished while maintaining mechanical integrity with respect to vibration and stress levels in the rotor. A full-scale cold-flow rig test was performed with metallic hardware fabricated to the specifications of the hot ceramic geometry to evaluate the stage performance.
The hydrostatic bearing promises life and speed characteristics currently unachievable with rolling element bearings alone. In order to achieve the speed and life requirements of the next generation of rocket engines, turbopump manufacturers are proposing hydrostatic bearings to be used in place of, or in series with, rolling element bearings. The design of a hydrostatic bearing is dependent on accurate pressure in the bearing. The stiffness and damping of the hydrostatic bearing is very sensitive to the bearing recess pressure ratio. In the conventional approach, usually ad hoc assumptions were made in determining the bearing pressure of this approach is inherently incorrect. In the present paper, a more elaborate approach to obtain bearing pressure is used. The bearing pressure and complete flow features of the bearing are directly computed by solving the complete 3-D Navier Stokes equation. The code used in the present calculation is a modified version of REACT3D code. Several calculations have been performed for the hydrostatic bearing designed and tested at Texas A&M. Good agreement has been obtained between computed and test results. Detailed flow features in the bearing will also be described and discussed.
Alfven waves are a major mechanism for the transport of electromagnetic energy from the distant part of the magnetosphere to the near-Earth space. This is especially true for the auroral and polar regions of the Earth. However, the mechanisms for their dissipation have remained illusive. One of the mechanisms is the formation of double layers when the current associated with Alfven waves in the inertial regime interact with density cavities, which either are generated nonlinearly by the waves themselves or are a part of the ambient plasma turbulence. Depending on the strength of the cavities, weak and strong double layers could form. Such double layers are transient; their lifetimes depend on that of the cavities. Thus they impulsively accelerate ions and electrons. Another mechanism is the resonant absorption of broadband Alfven- wave noise by the ions at the ion cyclotron frequencies. But this resonant absorption may not be possible for the very low frequency waves, and it may be more suited for electromagnetic ion cyclotron waves. A third mechanism is the excitation of secondary waves by the drifts of electrons and ions in the Alfven wave fields. It is found that under suitable conditions, the relative drifts between different ion species and/or between electrons and ions are large enough to drive lower hybrid waves, which could cause transverse accelerations of ions and parallel accelerations of electrons. This mechanism is being further studied by means of kinetic simulations using 2.5- and 3-D particle-in-cell codes. The ongoing modeling efforts on space weather require quantitative estimates of energy inputs of various kinds, including the electromagnetic energy. Our studies described here contribute to the methods of determining the estimates of the input from ubiquitous Alfven waves.
Pulsed detonation rocket engines (PDREs) have generated considerable research interest in recent years as a chemical propulsion system potentially offering improved performance and reduced complexity compared to conventional rocket engines. The detonative mode of combustion employed by these devices offers a thermodynamic advantage over the constant-pressure deflagrative combustion mode used in conventional rocket engines and gas turbines. However, while this theoretical advantage has spurred a great deal of interest in building PDRE devices, the unsteady blowdown process intrinsic to the PDRE has made realistic estimates of the actual propulsive performance problematic. The recent review article by Kailasanath highlights some of the difficulties in comparing the available experimental measurements with numerical models. In a previous paper by the author, parametric studies of the performance of a single, straight-tube PDRE were reported. A 1-D, unsteady method of characteristics code, employing a constant-gamma assumption behind the detonation front, was developed for that study. Models of this type are computationally inexpensive, and are particularly useful for parametric performance comparisons. For example, a plot showing the specific impulse of various PDRE and steady-state rocket engine (SSRE) configurations as a function of blowdown pressure ratio. The performance curves clearly indicate that a straight-tube PDRE is superior in specific impulse to a SSRE with a sonic nozzle over the entire range of pressure ratios. Note, however, that a straight-tube PDRE in general does not compare favorably to a SSRE fitted with an optimized de Laval supersonic nozzle, particularly at the high pressure ratios typical for boost or in-space rocket applications. However, the calculations also show that if a dynamically optimized, supersonic de Laval nozzle could be could be fitted to a PDRE, then the specific impulse of the device would exceed that of a comparable SSRE. While such a nozzle is a considerable idealization, it is clear that nozzle design and optimization will play a critical role in whether the performance potential of PDREs can be effectively realized in practice. In order to study PDRE nozzle issues with greater accuracy, a quasi-one-dimensional, finite-rate chemistry CFD code has been developed by the author. Comparisons of the code with both the previous MOC model and experimental data from Stanford University are reported. The effect of constant-gamma and finite-rate chemistry assumptions on the flowfield and performance is examined. Parametric studies of the effect of nozzle throat size and expansion ratio, at various blowdown pressure ratios, are reported.