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.
Final report for the DOE Fusion Energy Sciences project titled "A Streamlined, Open Source Neutronics Toolkit for Fusion Reactor Design". This project developed unstructured mesh tracking capabilities in OpenMC as well as two shutdown dose rate methodologies, and performed extensive performance improvements and validation of the open source Monte Carlo code, OpenMC.
Pre-equilibrium de-excitation is a well-established stage of nucleon-induced nuclear reactions but has not previously been incorporated into neutrino event generators. This work implements the Koning-Duijvestijn exciton model within the MARLEY neutrino event generator to simulate energy redistribution and particle emission before compound nucleus equilibrium. The implementation includes modular calculations of particle-hole state densities, internal transition rates, emission rates, and Monte Carlo cascade sampling, and is validated against the TALYS-2.2 nuclear reaction code. A new event record is also proposed to enable future uncertainty reweighting studies. This work establishes the first framework for studying pre-equilibrium effects in neutrino-induced reactions, with the goal of improving predictions of particle multiplicities, γ-ray production, and MeV-scale detector signals relevant to experiments such as DUNE.
Two validation campaigns were examined within the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) Phase 1 project to examine the modeling tools' underprediction of loads and motion of a floating wind semisubmersible (semi) at their surge and pitch natural frequencies. These campaigns were performed at the Maritime Research Institute Netherlands (MARIN) in 2017 and 2018. The load cases (LC) considered include: LC1 – Load measurements across semi under current loading; LC2 - Load measurements across semi under forced surge oscillation; LC3 – Load measurements across semi under wave loading, while held fixed; LC4 – Free-decay motion measurements in surge, pitch, and heave; and LC5 – Motion measurements under wave loading. Details on the results from the OC6 Phase Ia project can be found in the reference, “OC6 Phase 1: Investigating the underprediction of low-frequency hydrodynamic loads and responses of floating wind turbines”, J Phys: Conf Series 1618 032033.
Concern has been raised for the 64-m to 70-m antenna upgrade project that the 70-m system may experience greater S-band beam-pointing perturbations than the 64-m system. The S-band perturbations are due to minor (higher order) mode generation, causing subtle cross-polarization fields affecting beam pointing direction, as described herein. For the antennas in their present configuration (64 m), a slight S-band gain degradation of about 0.05 dB can be attributed to these effects. Therefore, a full physical optics analysis was performed for the present-day 64-m system, as described herein. The results were compared with past analyses and experimental observations in order to verify the algebra and computer code with the intent of deriving a valid analysis method for accurately analyzing the 70-m shaped dual reflector Cassegrainian antenna. The results of the new analysis appear to be in excellent agreement with previous analyses and experimental data.
Until recently, computations of space plasma flow over a spacecraft have been unstable for ratios of spacecraft dimension to Debye length typical of the low Earth orbit environment. Calculations are presented of the spacecraft/environment interaction based on two computer codes, MACH and POLAR. MACH, an inside-out particle tracking code, was developed for the purpose of validating the physics of POLAR in regimes where these are no comprehensive theoretical or experimental results. While the spacecraft which can be treated by MACH are restricted to simple geometries, the methodology is more fundamental than POLAR. MACH generates self-consistent solutions within the context of quasisteady Vlasov plasma flow and achieves Debye ratios previously unobtainable. POLAR uses a three-dimensional finite-element representation of the vehicle in a staggered mesh. The plasma sheath is modeled by outside-in particle tracking. Solutions for the plasma flow, wake and vehicle charging are obtained by Vlasov-Poisson iteration; charge stabilization techniques make the results virtually insensitive to the Debye ratio. POLAR reproduces the Laframboise static plasma solutions for sperical probes and fits the Makita-Kuriki probe data for spheres in a flowing plasma in regions where comparisons are valid. POLAR and MACH solutions for the particle and electrostatic potential structure of the wake of a charged disk in a low-altitude flow are shown for Mach numbers 4, 5, and 8. New features of the solutions include ion focussing in the wake and a definitive determination of the sheath edge in the wake which shows that the sheath is not an equipotential.
Concern has been raised for the 64-m to 70-m antenna upgrade project that the 70-m system may experience greater S-band beam-pointing perturbations than the 64-m system. The S-band perturbations are due to minor (higher order) mode generation, causing subtle cross-polarization fields affecting beam pointing direction, as described herein. For the antennas in their present configuration (64 m), a slight S-band gain degradation of about 0.05 dB can be attributed to these effects. Therefore, a full physical optics analysis was performed for the present-day 64-m system, as described herein. The results were compared with past analyses and experimental observations in order to verify the algebra and computer code with the intent of deriving a valid analysis method for accurately analyzing the 70-m shaped dual reflector Cassegrainian antenna. The results of the new analysis appear to be in excellent agreement with previous analyses and experimental data.
Single and multiple jet interaction in a combustor is simulated numerically to determine the effects of various k-epsilon models and boundary conditions as part of an attempt to understand the physics of the flow field inside the preburner of the Space Shuttle Main Engine. Due to the complexity of the physical and chemical phenomena involved, the model is simplified to two dimensions and a single-component nonreacting fluid. The computer code used is the modified TEACH code, which is a steady state isothermal fluid code using the SIMPLE algorithm. The physical validity and numerical accuracy of the model are established through comparisons with experimentally determined isothermal flow fields in a combustor.
The objective of the research is to develop and validate accurate, user-oriented viscous CFD codes (with inviscid options) for three-dimensional, unsteady aerodynamic flows about arbitrary rotorcraft configurations. Unsteady, three-dimensional Euler and Navier-Stokes codes are developed, adapted, and extended to rotor-body combinations. Flow solvers are coupled with zonal grid topologies, including rotating and nonrotating blocks. Special grid clustering and wave-fitting techniques were developed to capture low-level radiating acoustic waves. Significant progress was made in computing the propagation of acoustic waves due to the interaction of a concentrated vortex and a helicopter airfoil. The need for higher-order schemes was firmly established in relatively inexpensive two-dimensional calculations. In three dimensions, the number of grid points required to capture the low-level acoustic waves becomes very large, so that large supercomputer memory becomes essential. Good agreement was obtained between the numerical results obtained with a thin-layer Navier-Stokes code and experimental data from a model rotor. In addition, several nonrotating configurations that are sometimes proposed to simulate rotor blade tips in conventional wind tunnels were examined, and the complex flow around the radical tip shape of the world's fastest helicopter is under investigation. These studies demonstrate the flexibility and power of CFD to gain physical insight, study novel ideas, and examine various possibilities that might be difficult or impossible to set up in physical experiments. As a prelude to studies of rotor-body aerodynamic interactions, a preliminary grid topology and moving-interface strategy were developed. A new Euler/Navier-Stokes code using these techniques computes the vortical wake directly, rather than modeling it, as in most previous rotorcraft studies. Several hover cases were run for conventional and advanced-geometry blades. Numerical schemes using multi-zones and/or adaptive grids appear to be necessary to simulate the complex vortical flows in rotor wakes.
The Rich-burn/Quick-mix/Lean-burn (RQL) combustor was identified as a potential gas turbine combustor concept to reduce NO(x) emissions in High Speed Civil Transport (HSCT) aircraft. To demonstrate reduced NO(x) levels, cylindrical flametube versions of RQL combustors are being tested at NASA Lewis Research Center. A critical technology needed for the RQL combustor is a method of quickly mixing by-pass combustion air with rich-burn gases. Jet mixing in a cylindrical quick-mix section was numerically analyzed. The quick-mix configuration was five inches in diameter and employed twelve radial-inflow slots. The numerical analyses were performed with an advanced, validated 3-D Computational Fluid Dynamics (CFD) code named REFLEQS. Parametric variation of jet-to-mainstream momentum flux ratio (J) and slot aspect ratio was investigated. Both non-reacting and reacting analyses were performed. Results showed mixing and NO(x) emissions to be highly sensitive to J and slot aspect ratio. Lowest NO(x) emissions occurred when the dilution jet penetrated to approximately mid-radius. The viability of using 3-D CFD analyses for optimizing jet mixing was demonstrated.
Results of the evaluation of the static, buckling, and free vibration analyses capabilities of MHOST for the plate elements are presented. Two large scale, general purpose finite element codes (MARC and MSC/NASTRAN) are used to validate MHOST. Comparisons of MHOST results with those from MARC and MSC/NASTRAN show good agreement and indicate that MHOST can be used with confidence to perform the aforementioned analyses using the plate element.
This paper presents an overview of the primary nonintrusive diagnostic techniques being developed by the NASA Langley Research Center to address the validation needs of Computational Fluid Dynamic (CFD) codes. The techniques include absorption in the UV and IR, Laser Induced Fluorescence, electron beam fluorescence, and a number of scattering techniques including Rayleigh, spontaneous Raman, and several coherent Raman spectroscopies. Most of the techniques are highly specialized, require complex data interpretation, and can satisfy only a few of the CFD needs. For these reasons, the evolving trend in flowfield diagnostics appears to favor a mode in which the diagnostic researcher, the experimental aerodynamicist, and the CFD community jointly define experiments based on the aeronautical requirements and on available diagnostic techniques.
The project goals, plan and accomplishments up to this point are summarized in the viewgraphs. Among the various accomplishments, the most important have been: the development of the prismatic finite element code for doubly curved platforms and its validation with many different antenna configurations; the design and fabrication of a new slot spiral antennas suitable for automobile cellular, GPS and PCs communications; the investigation and development of various mesh truncation schemes, including the perfectly matched absorber and various fast integral equation methods; and the introduction of a frequency domain extrapolation technique (AWE) for predicting broadband responses using only a few samples of the response. This report contains several individual reports most of which have been submitted for publication to referred journals. For a report on the frequency extrapolation technique, the reader is referred to the UM Radiation Laboratory report A total of 14 papers have been published or accepted for publication with the full or partial support of this grant. Several more papers are in preparation.
Determining the baseline aerodynamics of the F/A-18 was one of the major objectives of the High-Angle-of-Attack Technology Program (HATP). This paper will review the key data bases that have contributed to our knowledge of the baseline aerodynamics and the improvements in test techniques that have resulted from the experimental program. Photographs are given highlighting the forebody and leading-edge-extension (LEX) vortices. Other data representing the impact of Mach and Reynolds numbers on the forebody and LEX vortices will also be detailed. The level of agreement between different tunnels and between tunnels and flight will be illustrated using pressures, forces, and moments measured on a 0.06-scale model tested in the Langley 7- by 10-Foot High Speed Tunnel, a 0.16-scale model in the Langley 30- by 60-Foot Tunnel, a full-scale vehicle in the Ames 80- by 120-Foot Wind Tunnel, and the flight F/A-18 High Alpha Research Vehicle (HARV). Next, creative use of wind tunnel resources that accelerated the validation of the computational fluid dynamics (CFD) codes will be described. Lastly, lessons learned, deliverables, and program conclusions are presented.
This report documents the structure and operation of a simulation model of the Benchmark Active Control Technology (BACT) Wind-Tunnel Model. The BACT system was designed, built, and tested at NASA Langley Research Center as part of the Benchmark Models Program and was developed to perform wind-tunnel experiments to obtain benchmark quality data to validate computational fluid dynamics and computational aeroelasticity codes, to verify the accuracy of current aeroservoelasticity design and analysis tools, and to provide an active controls testbed for evaluating new and innovative control algorithms for flutter suppression and gust load alleviation. The BACT system has been especially valuable as a control system testbed.
Performance characteristics of the ejector process within a turbine-based combined-cycle (TBCC) propulsion system are investigated using the NPARC Navier-Stokes code. The TBCC concept integrates a turbine engine with a ramjet into a single propulsion system that may efficiently operate from takeoff to high Mach number cruise. At the operating point considered, corresponding to a flight Mach number of 2.0, an ejector serves to mix flow from the ramjet duct with flow from the turbine engine. The combined flow then passes through a diffuser where it is mixed with hydrogen fuel and burned. Three sets of fully turbulent Navier-Stokes calculations are compared with predictions from a cycle code developed specifically for the TBCC propulsion system. A baseline ejector system is investigated first. The Navier-Stokes calculations indicate that the flow leaving the ejector is not completely mixed, which may adversely affect the overall system performance. Two additional sets of calculations are presented; one set that investigated a longer ejector region (to enhance mixing) and a second set which also utilized the longer ejector but replaced the no-slip surfaces of the ejector with slip (inviscid) walls in order to resolve discrepancies with the cycle code. The three sets of Navier-Stokes calculations and the TBCC cycle code predictions are compared to determine the validity of each of the modeling approaches.
Extensive creep testing of a hot-pressed silicon nitride (NC 132) was performed at 1300 C in air using five different specimen-loading configurations: (1) pure tension, (2) pure compression, (3) four-point uniaxial flexure, (4) ball-on-ring biaxial flexure, and (5) ring-on-ring biaxial flexure. This paper reports experimental results as well as test techniques developed in this work. Nominal creep strain and its rate for a given nominal applied stress were greatest in tension, least in compression, and intermediate in uniaxial and biaxial flexure. Except for the case of compression loading, nominal creep strain generally decreased with time, resulting in a less-defined steady-state condition. Of the four creep formulations-power-law, hyperbolic sine, step, and redistribution--the conventional power-law formulation still provides the most convenient and reasonable estimation of the creep parameters of the NC 132 material. The data base to be obtained will be used to validate the NASA Glenn-developed design code CARES/Creep (ceramics analysis and reliability evaluation of structures and creep).
A design of experiments approach has been implemented using computational hypervelocity impact simulations to determine the most effective place to add mass to an existing metallic Thermal Protection System (TPS) to improve hypervelocity impact protection. Simulations were performed using axisymmetric models in CTH, a shock-physics code developed by Sandia National Laboratories, and validated by comparison with existing test data. The axisymmetric models were then used in a statistical sensitivity analysis to determine the influence of five design parameters on degree of hypervelocity particle dispersion. Several damage metrics were identified and evaluated. Damage metrics related to the extent of substructure damage were seen to produce misleading results, however damage metrics related to the degree of dispersion of the hypervelocity particle produced results that corresponded to physical intuition. Based on analysis of variance results it was concluded that the most effective way to increase hypervelocity impact resistance is to increase the thickness of the outer foil layer. Increasing the spacing between the outer surface and the substructure is also very effective at increasing dispersion.
Computational Fluid Dynamics (CFD) flow solutions for the M2129 diffusing S-duct with and without vane effectors were computed by the Wind-US flow solver. Both structured and unstructured 3-D grids were used. Without vane effectors, the duct exhibited massive flow separation in both experiment and CFD. With vane effectors installed, the flow remained attached and aerodynamic losses were reduced. Total pressure recovery and distortion near the duct outlet were computed from the solutions and compared favorably to experimental values. These calculations are part of a validation effort for the Wind-US code. They also provide an example case to aid engineers in learning to use the Wind-US software.