Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source Code”

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,009 records · Page 56

Inviscid and Viscous CFD Analysis of Booster Separation for the Space Launch System Vehicle

This paper presents details of Computational Fluid Dynamic (CFD) simulations of the Space Launch System during solid-rocket booster separation using the Cart3D inviscid and Overflow viscous CFD codes. The discussion addresses the use of multiple data sources of computational aerodynamics, experimental aerodynamics, and trajectory simulations for this critical phase of flight. Comparisons are shown between Cart3D simulations and a wind tunnel test performed at NASA Langley Research Center's Unitary Plan Wind Tunnel, and further comparisons are shown between Cart3D and viscous Overflow solutions for the flight vehicle. The Space Launch System (SLS) is a new exploration-class launch vehicle currently in development that includes two Solid Rocket Boosters (SRBs) modified from Space Shuttle hardware. These SRBs must separate from the SLS core during a phase of flight where aerodynamic loads are nontrivial. The main challenges for creating a separation aerodynamic database are the large number of independent variables (including orientation of the core, relative position and orientation of the boosters, and rocket thrust levels) and the complex flow caused by exhaust plumes of the booster separation motors (BSMs), which are small rockets designed to push the boosters away from the core by firing partially in the direction opposite to the motion of the vehicle.

Viscous↗

CFD Simulation of an Inward-Turning Supersonic Inlet Unstart at Flight Mach Number 1.7

Interest in inward-turning supersonic inlets renewed during NASA’s Commercial Supersonic Technology Project, as a means of reducing disturbances caused by propulsion integration. Inward-turning designs have the potential to greatly reduce the strength of shock waves emanating outward from the nacelle, and therefore their contribution to the overall sonic boom signature. For top-mounted, or over-wing installations, disturbances to aircraft aerodynamics are reduced. The nature of the inward-turning geometry results in a boundary-layer thickening effect, opposite that of the familiar thinning of an outward-turning cone boundary-layer. The application of a small amount of boundary-layer bleed at the throat has been shown to alleviate the negative effects of the boundary-layer on subsonic diffuser performance. Boundary-layer bleed, however, comes at the expense of increased complexity, and external disturbances at bleed exits that would partially defeat the purpose of the inward-turning scheme. Vortex generators were proposed as an alternative to bleed, and an initial CFD study of their effectiveness in this application forms the basis of the present work. A source-term model in the Wind-US CFD code was used to model the effects of an array of vortex generators placed at the entrance to the subsonic diffuser. Baseline results with no vortex generators compared favorably to previously reported results with no bleed at back-pressures near the design point, but stability margin in the present case was limited by a curious non-linear phenomena similar to unstart in mixed-compression inlets. While some improvement was noted with vortex generators, their effectiveness was diminished by the limited stability margin. A time-accurate calculation was subsequently performed to better understand this process. Results of this calculation are presented and instantaneous solutions compared with steady-state, subcritical results previously reported. The stability margin limitation is attributed to the use of a constant-pressure outflow boundary condition, which differed from that used in the prior simulations.

Pre-Compressor Cooling, Water Injection, Hypersoni↗

Where IMERG Goes Next: Version 08 and Beyond

With the Version 07 (V07) Integrated Multi-satellitE Retrievals for GPM (IMERG) algorithm finalized and production initiated, the focus turns to enhancements for Version 08. These include innovations not included in V07 due to time constraints, plus issues revealed by the initial V07 products. One high priority is to evaluate and revise the schemes in V07 that rectify temporal artifacts caused by the time interpolation that fills the gaps between the various passive microwave (PMW) sensor overpasses. A second priority is to improve the homogeneity between the TRMM and GPM eras by characterizing differences between the two eras, determining the causes of these differences, and applying corrections as feasible, perhaps by enforcing spatial scale consistency (an overarching issue). Certainly, we must account for GPROF and the Combined Radar-Radiometer Algorithm converting to Machine Learning schemes in V08. Other priority topics include additional automated quality control for artifacts in the IR brightness temperatures and PMW precipitation fields, revisions to the specification algorithm for the probability of liquid precipitation, and accommodating new PMW sensors, which include the next generation of small-sats. We also consider the post-V08 landscape; the final GPM reprocessing will be restricted to fixing known code or algorithmic errors. Nonetheless, there are several data sources on the horizon to consider, including more small-sat PMW radiometers, AVHRR-based precipitation estimates (most useful in high latitudes), and the ISCCP-Next Generation and GEO-Ring projects that could provide easy access to multiple geosynchronous satellite channels and enable significantly improved algorithms compared to GEO-IR alone.

George J. Huffman↗

New features and applications of PRESTO, a computer code for the performance of regenerative, superheated steam turbine cycles

The code was designed to analyze performance at valves-wide-open design flow. The code can model conventional steam cycles as well as cycles that include such special features as process steam extraction and induction and feedwater heating by external heat sources. Convenience features and extensions to the special features were incorporated into the PRESTO code. The features are described, and detailed examples illustrating the use of both the original and the special features are given.

Choo, Y. K.↗

Monte Carlo simulation of a new gamma ray telescope

A new Monte Carlo code has been written to simulate the response of the new University of California double scatter gamma ray telescope. This package of modular software routines, written in VAX FORTRAN 77 simulates the detection of 0.1 to 35 MeV gamma rays. The new telescope is flown from high altitude balloons to measure medium energy gamma radiation from astronomical sources. This paper presents (1) the basic physics methods in the code, and (2) the predicted response functions of the telescope. Gamma ray processes include Compton scattering, pair production and photoelectric absorption in plastic scintillator, NaI(Tl) and aluminum. Electron transport processes include ionization energy loss, multiple scattering, production of bremsstrahlung photons and positron annihilation.

Simone, J.↗

Energy Spectrum and Mass Composition of Ultra-high-energy Cosmic Rays Originating from Relativistic Jets of Nearby Radio Galaxies

Relativistic jets of radio galaxies (RGs) are possible sources of ultra-high-energy cosmic rays (UHECRs). Recent studies combining relativistic hydrodynamic simulations with Monte Carlo particle transport have demonstrated that UHECRs can be accelerated to energies beyond 10 20 eV through shocks, turbulence, and relativistic shear in jet-induced flows of Fanaroff–Riley type RGs. The resulting time-asymptotic UHECR spectrum is well modeled by a double power law with an “extended” exponential cutoff, primarily shaped by relativistic shear acceleration. In this study, we adopt this novel source spectrum and simulate the propagation of UHECRs from nearby RGs using the CRPropa code. We focus on Virgo A (Vir A), Centaurus A (Cen A), Fornax A (For A), and Cygnus A (Cyg A), expected to be the most prominent UHECR sources among RGs. We then analyze the energy spectrum and mass composition of UHECRs arriving at Earth. We find that, due to the extended high-energy tail in the source spectrum, UHECRs from Vir A, which has a higher Lorentz factor, exhibit a higher flux at the highest energies and a lighter mass composition at Earth compared to those from Cen A and For A with lower Lorentz factors. Despite Cyg A having an even higher Lorentz factor, the large distance limits its contribution. With a small number of nearby prominent RGs, our findings suggest that if RGs are the major sources of UHECRs, the energy spectrum and mass composition of observed UHECRs would exhibit hemispheric differences between the Northern and Southern skies at the highest energies.

79 ASTRONOMY AND ASTROPHYSICS↗

Use of the MATRIXx Integrated Toolkit on the Microwave Anisotropy Probe Attitude Control System

Recent advances in analytical software tools allow the analysis, simulation, flight code, and documentation of an algorithm to be generated from a single source, all within one integrated analytical design package. NASA's Microwave Anisotropy Probe project has used one such package, Integrated Systems' MATRIXx suite, in the design of the spacecraft's Attitude Control System. The project's experience with the linear analysis, simulation, code generation, and documentation tools will be presented and compared with more traditional development tools. In particular, the quality of the flight software generated will be examined in detail. Finally, lessons learned on each of the tools will be shared.

Ward, David K.↗

Regular Topologies for Gigabit Wide-Area Networks: Congestion Avoidance Testbed Experiments

This document is Volume 3 of the final technical report on the work performed by SRI International (SRI) on SRI Project 8600. The document includes source listings for all software developed by SRI under this effort. Since some of our work involved the use of ST-II and the Sun Microsystems, Inc. (Sun) High-Speed Serial Interface (HSI/S) driver, we have included some of the source developed by LBL and BBN as well. In most cases, our decision to include source developed by other contractors depended on whether it was necessary to modify the original code. If we have modified the software in any way, it is included in this document. In the case of the Traffic Generator (TG), however, we have included all the ST-II software, even though BBN performed the integration, because the ST-II software is part of the standard TG release. It is important to note that all the code developed by other contractors is in the public domain, so that all software developed under this effort can be re-created from the source included here.

Denny, Barbara A.↗

E-Standards For Mass Properties Engineering

A proposal is put forth to promote the concept of a Society of Allied Weight Engineers developed voluntary consensus standard for mass properties engineering. This standard would be an e-standard, and would encompass data, data manipulation, and reporting functionality. The standard would be implemented via an open-source SAWE distribution site with full SAWE member body access. Engineering societies and global standards initiatives are progressing toward modern engineering standards, which become functioning deliverable data sets. These data sets, if properly standardized, will integrate easily between supplier and customer enabling technically precise mass properties data exchange. The concepts of object-oriented programming support all of these requirements, and the use of a JavaTx based open-source development initiative is proposed. Results are reported for activity sponsored by the NASA Langley Research Center Innovation Institute to scope out requirements for developing a mass properties engineering e-standard. An initial software distribution is proposed. Upon completion, an open-source application programming interface will be available to SAWE members for the development of more specific programming requirements that are tailored to company and project requirements. A fully functioning application programming interface will permit code extension via company proprietary techniques, as well as through continued open-source initiatives.

Cerro, Jeffrey A.↗

Point Spread Function and Transmittance Analyses for Conical and Hexapod Secondary Mirror Support Towers for the Next Generation Space Telescope (NGST)

Different secondary mirror support towers were modeled on the CODE V optical design/analysis program for the NGST Optical Telescope Assembly (OTA) B. The vertices of the NGST OTA B primary and secondary mirrors were separated by close to 9.0 m. One type of tower consisted of a hollow cone 6.0 m long, 2.00 m in diameter at the base, and 0.704 m in diameter at its top. The base of the cone was considered attached to the primary's reaction structure through a hole in the primary. Extending up parallel to the optical axis from the top of this cylinder were eight blades (pyramidal struts) 3.0 m long. A cross section of each these long blades was an isosceles triangle with a base of 0.010 m and a height of 0.100 m with the sharpest part of each triangle pointing inward. The eight struts occurred every 45 deg. The other type of tower was purely a hexapod arrangement and had no blades or cones. The hexapod consisted simply of six, very thin, circular struts, leaving in pairs at 12:00, 4:00, and 8:00 at the primary and traversing to the outer edge of the back of the secondary mount. At this mount, two struts arrived at each of 10:00, 2:00, and 6:00. The struts were attached to the primary mirror in a ring 3.5 m in diameter. They reached the back of the secondary mount, a circle 0.704 m in diameter. Transmittance analyses at two levels were performed on the secondary mirror support towers. Detailed transmittances were accomplished by the use of the CODE V optical design/analysis program and were compared to transmittance calculations that were almost back-of-the-envelope. Point spread function (PSF) calculations, including both diffraction and aberration effects, were performed on CODE V. As one goes out from the center of the blur (for a point source), the two types of support towers showed little difference between their PSF intensities until one reaches about the 3% level. Contours can be delineated on CODE V down to about 10 (exp -8) times the peak intensity, fine detailing the differences in the wings of the PSF. Between the first type of tower and the hexapod tower, below 3% of the peak intensity, the PSF contour maps well delineated significant differences in the outer portions of the PSFS. Such differences in the wings of the PSF significantly impact the resolution of brown dwarfs or extra solar planets from their star. The hexapod tower provided the best PSF.

Wilkerson, Gary W.↗

A liquid xenon imaging telescope for gamma ray astrophysics: Design and expected performance

A high resolution telescope for imaging cosmic x ray sources in the MeV region, with an angular resolution better than 0.5 deg is being developed as balloon-borne payload. The instrument consists of a 3-D liquid xenon TPC as x ray detector, coupled with a coded aperture at a distance of 1 meter. A study of the actual source distribution of the 1.809 MeV line from the decay of Al-26 and the 511 keV positron-electron annihilation line is among the scientific objectives, along with a search for new x ray sources. The telescope design parameters and expected minimum flux sensitivity to line and continuum radiation are presented. The unique capablity of the LXe-TPC as a Compton Polarimeter is also discussed.

Aprile, E.↗

Development of an upwind, finite-volume code with finite-rate chemistry

Under this grant, two numerical algorithms were developed to predict the flow of viscous, hypersonic, chemically reacting gases over three-dimensional bodies. Both algorithms take advantage of the benefits of upwind differencing, total variation diminishing techniques, and a finite-volume framework, but obtain their solution in two separate manners. The first algorithm is a zonal, time-marching scheme, and is generally used to obtain solutions in the subsonic portions of the flow field. The second algorithm is a much less expensive, space-marching scheme and can be used for the computation of the larger, supersonic portion of the flow field. Both codes compute their interface fluxes with a temporal Riemann solver and the resulting schemes are made fully implicit including the chemical source terms and boundary conditions. Strong coupling is used between the fluid dynamic, chemical, and turbulence equations. These codes have been validated on numerous hypersonic test cases and have provided excellent comparison with existing data.

Molvik, Gregory A.↗

Development of an upwind, finite-volume code with finite-rate chemistry

Under this grant, two numerical algorithms were developed to predict the flow of viscous, hypersonic, chemically reacting gases over three-dimensional bodies. Both algorithms take advantage of the benefits of upwind differencing, total variation diminishing techniques and of a finite-volume framework, but obtain their solution in two separate manners. The first algorithm is a zonal, time-marching scheme, and is generally used to obtain solutions in the subsonic portions of the flow field. The second algorithm is a much less expensive, space-marching scheme and can be used for the computation of the larger, supersonic portion of the flow field. Both codes compute their interface fluxes with a temporal Riemann solver and the resulting schemes are made fully implicit including the chemical source terms and boundary conditions. Strong coupling is used between the fluid dynamic, chemical and turbulence equations. These codes have been validated on numerous hypersonic test cases and have provided excellent comparison with existing data. This report summarizes the research that took place from August 1,1994 to January 1, 1995.

Molvik, Gregory A.↗

The Seismic Signatures of Recently Formed Impact Craters on Mars

We investigated the seismic signatures of recent impact crater clusters on Mars that would be recorded by the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight)seismometers. We used a database of 77 measured and dated impact sites, with craters with diameters between 2.1 and 33.8 m, along with inferred impact angle, bolide trajectory, and varying target material properties to empirically scale for the momentum, expected seismic source function, and radiation pattern of impacts. The impact source is simulated in a local 3‐Dfinite difference wave propagation code and coupled to teleseismic distances by scaling the spectra of 1‐D global synthetic seismograms. We use the InSight seismometer noise floors to estimate detectability of impact(s) across azimuth and distance. Our experiments reveal that impact clusters have a higher peak corner frequency resulting from energy contributed by smaller craters to the power spectrum. We also find that the time separation between individual impacts in a cluster is small (<10–15 ms) and requires a seismometer closely situated to the source (<10 km) and a high sampling rate (>100 Hz) to resolve individual impacts within the cluster. Two of the clusters in our database(>20‐m effective diameter) would have been detectable by InSight, with the assumptions that the Martian background noise and seismic attenuation are both low. Joint detection of surface changes from newly formed crater(s) in images and by InSight will provide precise source locations that are crucial for constraining the internal structure of Mars

N. C. Schmerra↗

Chasing ghosts: characterization of artifact generation in coded aperture decoding due to experimental implementation

Coded aperture imaging is a form of lensless aperture imaging that projects multiple overlapping images of the source onto the detector, enhancing signal strength, which is advantageous for low-flux sources or high-resolution imaging. This technique requires decoding of the detector signal to reconstruct the original source, which involves convolving the detector data with the aperture pattern. When the signal is from a centered point source, the reconstructed source image is known as the point spread function (PSF). A clean PSF without artifacts is a Dirac delta function [Appl. Opt. 20, 1858 (1981)]. This paper examines the robustness of the decoding process against variations in experimental tolerances by analyzing artifact growth in the reconstructed PSF. We illustrate the effects of incorrect magnification, rotation, and detector size and find that aperture–detector rotational misalignment about the imaging axis is the most sensitive parameter, with significant artifact generation occurring with angular offsets of less than one degree. We discuss compensation methods for imperfect aperture placement, finding that small detector sizes produce uncompensatable artifact generation, and compare theoretical predictions with experimental PSF measurements of a rank , 6.8 mm thick (less than one mean free path) coded aperture with a 3.5 mm cell size, conducted at the MegaJOuLe Neutron Imaging Radiography dense plasma focus [IEEE Trans. Plasma Sci. 49, 3299 (2021)] using a 2.45 MeV neutron source. Based on our findings, we recommend using magnified coded apertures in the under-sampled regime, which allows for the inclusion of fiducial markers to characterize aperture–detector rotational offsets and the addition of mechanical coupling, where possible, to constrain rotational and magnification offsets.

Selwood, M. P. [Lawrence Livermore National Labora↗

User's manual for airfoil flow field computer code SRAIR

A two dimensional unsteady Navier-Stokes calculation procedure with specific application to the isolated airfoil problem is presented. The procedure solves the full, ensemble averaged Navier-Stokes equations with turbulence represented by a mixing length model. The equations are solved in a general nonorthogonal coordinate system which is obtained via an external source. Specific Cartesian locations of grid points are required as input for this code. The method of solution is based upon the Briley-McDonald LBI procedure. The manual discusses the analysis, flow of the program, control steam, input and output.

Shamroth, S. J.↗

A Python Script to Read MCNP6.3 Surface-Source Files

This report provides a Python script to read an MCNP ® surface source file created with the SSW card with SYM = 0 (i.e., the default symmetry treatment). For background: the general format of an MCNP surface-source file is described in; however, that document did not provide coding and/or a tool to interrogate such files. The current format will not be given in this document other than through the record-read statements necessary for the script to function. The reader capability in this report is augmented with the ability to directly write a couple demonstrative outputs: 1. A comma-separated value (CSV) file containing particle phase-space state information and 2. A Matplotlib histogram of the energy distribution of the particles. This report also describes accompanying verification work that shows the script performing as required with MCNP6.2, MCNP6.3, and (expected) MCNP6.4 surface-source files. However, users of the enclosed script must still verify that the script is behaving correctly for their own work.

97 MATHEMATICS AND COMPUTING↗

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗