THE N-BODY CODE - A GENERAL FORTRAN CODE FOR SOLUTION OF PROBLEMS IN SPACE MECHANICS BY NUMERICAL METHODS
Astronomical integration code - fortran solution of space mechanics - computer program
FORTRAN↗
SEARCH · Engineering Papers
Results for “V 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.
Numerical analysis of stiffened shells of revolution. Volume 4: Engineer's program manual for STARS-2S shell theory automated for rotational structures - 2 (statics) digital computer program
The engineering programming information for the digital computer program for analyzing shell structures is presented. The program is designed to permit small changes such as altering the geometry or a table size to fit the specific requirements. Each major subroutine is discussed and the following subjects are included: (1) subroutine description, (2) pertinent engineering symbols and the FORTRAN coded counterparts, (3) subroutine flow chart, and (4) subroutine FORTRAN listing.
Svalbonas, V.↗
Results of a flight experiment on the catalytic efficiency of the Space Shuttle heat shield
An experiment was performed on the second Space Shuttle flight to verify data from arc-heated wind-tunnel tests indicating that the thermal protection tiles used on the Shuttle's surface are noncatalytic to the recombination of dissociated air. The experiment indicates that the tiles are indeed noncatalytic, and explains the lower than expected temperatures observed in the first two Shuttle flights. The data are reasonably well predicted with computations, using a chemically reactive boundary-layer code. The experimental setup and the computational analysis are described.
Rakich, J. V.↗
Deep space optical communication via relay satellite
The application of optical communications for a deep space link via an earth-orbiting relay satellite is discussed. The system uses optical frequencies for the free-space channel and RF links for atmospheric transmission. The relay satellite is in geostationary orbit and contains the optics necessary for data processing and formatting. It returns the data to earth through the RF terrestrial link and also transmits an optical beacon to the satellite for spacecraft return pointing and for the alignment of the transmitting optics. Future work will turn to modulation and coding, pointing and tracking, and optical-RF interfacing.
Dolinar, S.↗
Application of advanced computational codes in the design of an experiment for a supersonic throughflow fan rotor
Increased emphasis on sustained supersonic or hypersonic cruise has revived interest in the supersonic throughflow fan as a possible component in advanced propulsion systems. Use of a fan that can operate with a supersonic inlet axial Mach number is attractive from the standpoint of reducing the inlet losses incurred in diffusing the flow from a supersonic flight Mach number to a subsonic one at the fan face. The design of the experiment using advanced computational codes to calculate the components required is described. The rotor was designed using existing turbomachinery design and analysis codes modified to handle fully supersonic axial flow through the rotor. A two-dimensional axisymmetric throughflow design code plus a blade element code were used to generate fan rotor velocity diagrams and blade shapes. A quasi-three-dimensional, thin shear layer Navier-Stokes code was used to assess the performance of the fan rotor blade shapes. The final design was stacked and checked for three-dimensional effects using a three-dimensional Euler code interactively coupled with a two-dimensional boundary layer code. The nozzle design in the expansion region was analyzed with a three-dimensional parabolized viscous code which corroborated the results from the Euler code. A translating supersonic diffuser was designed using these same codes.
Wood, Jerry R.↗
Application of advanced computational codes in the design of an experiment for a supersonic throughflow fan rotor
Increased emphasis on sustained supersonic or hypersonic cruise has revived interest in the supersonic throughflow fan as a possible component in advanced propulsion systems. Use of a fan that can operate with a supersonic inlet axial Mach number is attractive from the standpoint of reducing the inlet losses incurred in diffusing the flow from a supersonic flight Mach number to a subsonic one at the fan face. The design of the experiment using advanced computational codes to calculate the components required is described. The rotor was designed using existing turbomachinery design and analysis codes modified to handle fully supersonic axial flow through the rotor. A two-dimensional axisymmetric throughflow design code plus a blade element code were used to generate fan rotor velocity diagrams and blade shapes. A quasi-three-dimensional, thin shear layer Navier-Stokes code was used to assess the performance of the fan rotor blade shapes. The final design was stacked and checked for three-dimensional effects using a three-dimensional Euler code interactively coupled with a two-dimensional boundary layer code. The nozzle design in the expansion region was analyzed with a three-dimensional parabolized viscous code which corroborated the results from the Euler code. A translating supersonic diffuser was designed using these same codes.
Wood, Jerry R.↗
Comparison of Coupled Radiative Flow Solutions with Project Fire 2 Flight Data
A nonequilibrium, axisymmetric, Navier-Stokes flow solver with coupled radiation has been developed for use in the design or thermal protection systems for vehicles where radiation effects are important. The present method has been compared with an existing now and radiation solver and with the Project Fire 2 experimental data. Good agreement has been obtained over the entire Fire 2 trajectory with the experimentally determined values of the stagnation radiation intensity in the 0.2-6.2 eV range and with the total stagnation heating. The effects of a number of flow models are examined to determine which combination of physical models produces the best agreement with the experimental data. These models include radiation coupling, multitemperature thermal models, and finite rate chemistry. Finally, the computational efficiency of the present model is evaluated. The radiation properties model developed for this study is shown to offer significant computational savings compared to existing codes.
Olynick, David R.↗
Comparison of Coupled Radiative Navier-Stokes Flow Solutions with the Project Fire II Flight Data
A nonequilibrium, axisymmetric, Navier-Stokes flow solver with coupled radiation has been developed to use in the design of thermal protection systems for vehicles where radiation effects are important. The present method has been compared with an existing flow and radiation solver and with the Project Fire II experimental data. Very good agreement has been obtained over the entire Fire II trajectory with the experimentally determined values of the stagnation radiation intensity in the .2 to 6.2 eV range and with the total stagnation heating. The agreement was significantly better than previous numerical predictions. The effects of a number of flow models are examined to determine which combination of physical models produces the best agreement with the experimental data. These models include radiation coupling, multi-temperature thermal models, finite-rate chemistry, and a quasi-steady-state or Boltzmann assumption for the calculation of the excited electronic states. Finally, the computational efficiency of the present model is evaluated. The radiation properties model developed for this study is shown to offer significant computational savings compared to existing codes.
Olynick, David R.↗
JPSS-1 ATMS Post-launch Active Geolocation Analysis
A NOAA-20 (N20) ATMS active geolocation test was performed around Jan. 2018 with 24 pre-selected coastline crossing scenes. After comprehensive analysis of ATMS stare data and the corresponding VIIRS data, the ATMS pitch, roll and yaw pointing angle errors are found from the nadir perpendicular, the nadir oblique shallow angle, and the off-nadir perpendicular coastline crossing data, respectively. In this study, we first determine the ATMS radiometric coastline crossing time by using the ATMS radiometric count data. Since the coastline can be located anywhere within one ATMS FOV from the passive (regular scanning) geolocation data, depending on the scan starting time, it is not a valid assumption for the inflection point being the same as the coastline location. Consequently, using the passive geolocation data to validate the sensor’s on-orbit pointing angle performance is limited. After finding the ATMS radiometric coastline crossing time from the ATMS data, we compare it with the VIIRS effective time stamp (see main text). Specifically the VIIRS M3, M4 & M5 (true color) & M15 and M16 bands (thermal) data have a much smaller footprint size. Using the differences of the ATMS and VIIRS (effective) coastlines crossing times, the N20 ATMS pitch, roll and yaw pointing angle errors are found to be -0.09o, -0.24o and 0.28o, respectively. To determine ATMS geolocation properly, these on-orbit pointing errors need to be corrected, adding to the ATMS SDR Processing Coefficient Table, and passed on to the operational geolocation processing code.
Active Geolocation↗
Multiphysics Code Coupling for Fusion Blankets
Explore the source record for details and available documents.
Nuclear science and engineering - General studies ↗
Formation and Evolution of Binary Systems Containing Collapsed Stars
This research includes theoretical studies of the formation and evolution of five types of interacting binary systems. Our main focus has been on developing a number of comprehensive population synthesis codes to study the following types of binary systems: (i) cataclysmic variables (#3, #8, #12, #15), (ii) low- and intermediate-mass X-ray binaries (#13, #20, #21), (iii) high-mass X-ray binaries (#14, #17, #22), (iv) recycled binary millisecond pulsars in globular clusters (#5, #10, #ll), and (v) planetary nebulae which form in interacting binaries (#6, #9). The numbers in parentheses refer to papers published or in preparation that are listed in this paper. These codes take a new unified approach to population synthesis studies. The first step involves a Monte Carlo selection of the primordial binaries, including the constituent masses, and orbital separations and eccentricities. Next, a variety of analytic methods are used to evolve the primary star to the point where either a dynamical episode of mass transfer to the secondary occurs (the common envelope phase), or the system evolves down an alternate path. If the residual core of the primary is greater than 2.5 solar mass, it will evolve to Fe core collapse and the production of a neutron star and a supernova explosion. In the case of systems involving neutron stars, a kick velocity is chosen randomly from an appropriate distribution and added to the orbital dynamics which determine the state of the binary system after the supernova explosion. In the third step, all binaries which commence stable mass transfer from the donor star (the original secondary in the binary system) to the compact object, are followed with a detailed binary evolution code. Finally, we include all the relevant dynamics of the binary system. For example, in the case of LMXBs, the binary system, with its recoil velocity from the supernova explosion, is followed in time through its path in the Galactic potential. For our globular cluster study, detailed three-body calculations are used to follow the interactions of the binary with passing field stars.
Rappaport, Saul↗
Magnetic field induced structural changes in magnetite observed by resonant x-ray diffraction and Mössbauer spectroscopy
When a magnetic field is applied to a single crystal of magnetite at 124 K > T > 50 K , the monoclinic c M axis, which is the easy magnetization axis, switches to one of the $\langle$100$\rangle$ cubic directions, nearest to the direction of the magnetic field, and the phenomenon known as an axis switching (AS) occurs. Here, a global symmetry probe, resonant x-ray scattering, and a local probe, Mössbauer spectroscopy, are used to better understand the mechanism of axis switching. The behavior of three subsystems ordered below the Verwey transition temperature T V , i.e., lattice distortion, an orbital, and charge orderings, was observed via resonant x-ray scattering as a function of an external magnetic field. This was preceded by calculation of selected peak intensities using the fdmnes code. The Mössbauer spectroscopy studies confirmed that the magnetic field triggers electronic rearrangements and atomic displacements. The structure observed after the process of axis switching is very similar to the one obtained after cooling below T V with the magnetic field applied along one of the initial $\langle$100$\rangle$ cubic directions and distinct from the cooling in the absence of a magnetic field. From all the experimental observations of the phenomenon done so far, it is clear that AS starts from the fluctuations between octahedral iron orbitals that ultimately lead to the Verwey transition, but also to the higher-temperature trimeron dynamics. Therefore, further observation of the axis switching may be a key point to the understanding of a majority of strongly correlated electronic behavior in magnetite as well as in other transition metal oxides.
36 MATERIALS SCIENCE↗
Implicit highly-coupled single-ion Hall-MHD formulation for hybrid particle-in-cell codes
The rudiments of a particle-based single-fluid two-temperature magnetohydrodynamic (MHD) algorithm have been outlined in Thoma et al. (2013). The extension of this algorithm to include the effect of Hall physics is described in this paper. An implicit leapfrog version of the algorithm, which allows timesteps large compared to the resistive decay time and other relevant timescales, has recently been added to a hybrid particle-in-cell code. In standard MHD the Hall term in the generalized Ohm’s law can often be neglected when the Hall parameter is small. This term must, however, be retained in regimes where it is non-negligible. The retention of displacement current in Maxwell’s equations avoids the numerical difficulties associated with the whistler mode, which are encountered in standard explicit Hall-MHD codes, and allows the algorithm to be incorporated into hybrid particle-in-cell codes, for which particles may migrate from a kinetic to fluid to MHD description based upon local ambient plasma conditions. A highly-coupled implicit Hall-MHD formalism is presented, in which displacement current can either be retained or neglected. Even when displacement current is neglected, the highly-coupled implicit formalism avoids the restrictive timesteps for the whistler mode in explicit Hall-MHD codes. A comparison of numerical and analytic dispersion analysis demonstrates the feasibility of this approach and establishes relevant constraints to assure numerical stability. The implementation of the algorithm is described, and test simulation results in 1D and 2D in both linear and nonlinear regimes are presented.
97 MATHEMATICS AND COMPUTING↗
LSPRAY-V: A Lagrangian Spray Module
LSPRAY-V is a Lagrangian spray solver developed for application with unstructured grids and massively parallel computers. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray encountered over a wide range of operating conditions in modern aircraft engine development. It could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. The manual provides the user with an understanding of various models involved in the spray formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers. With the development of LSPRAY-V, we have advanced the state-of-the-art in spray computations in several important ways.
CFD↗
MOSCATO Development and Integration in Fiscal Year 2022
During FY21, we conducted ongoing development work for the MOSCATO (Molten Salt Chemistry and Transport) solver. The code development work primarily consisted of transitioning capabilities from the original version of the solver, which was written in OpenFOAM, into Nek5000. In doing so, a fast, highly parallelizable solver was created that is capable of complex chemistry and corrosion simulations for engineering-scale molten salt systems. The Nek5000 version of MOSCATO is now fully featured and capable of higher-fidelity simulations than were previously possible. Demonstration cases including a thermal convection loop have been simulated to test these new capabilities. We built upon the work for FY22 and improved the code from several different perspectives. First, we improved the user interface by adding a new component to the official Nek5000 input file (.par). This new part contains documents parameters like, salt properties (density, viscosity, Cp, thermal conductivity), diffusion coefficients, standard potential, etc. Second, we built a conversion script to extract salt properties from the MSTDB-TP salt database and write to MOSCATO input file. Third, we migrated the code to NekRS, which is the GPU branch of Nek5000 and suitable for next generation supercomputers. Verification and Validation (V&V) work was also continued in FY22. Two tasks were performed. The first V&V task involved the validation of the Poisson-Nernst-Planck equation solver and Butler-Volmer electrode kinetics, by comparing with numerical and experimental data about thermoelectric cells. The second task involved the comparisons to corrosion results from a thermal convection loop run during the MSRE era. Satisfactory agreement was obtained from both tasks.
Yuan, Haomin↗
ICME and In-Situ Process Monitoring for Rapid Qualification of Components Made by Laser-based Powder Bed Additive Manufacturing Processes for Nuclear Structural Applications
Additive manufacturing (AM) through selective powder bed melting of successive layers is being considered as disruptive technology for rapid production of low-cost nuclear reactor internal components with complex geometries. However, there is a potential for non-uniform distribution of physical features, such as porosity or microstructural differences, due to variations of temperature across a component build. These heterogeneities make the deployment of traditional material qualification and non-destructive evaluation of AM components difficult. This research project explored the feasibility of using in-situ process monitoring methods and integrated computational materials engineering (ICME) principles as an alternate qualification methodology and approach. This project included six individual tasks: (i) design of artifacts relevant to the nuclear power industry, (ii) evaluation of laser processing and in-situ measurements, (iii) computational modeling, (iv) ex-situ microstructural characterization, (v) evaluation of scaling the methodology for large-scale structures, and (vi) development of a data package to codes and standards organizations. As a part of this research, laser powder bed fusion of metals (316L, Alloy 718, and Ti6Al4V) was explored. This document provides an overview of the research performed over three years and serves as the final U.S. Department of Energy report for this project. In the third year, the project focused on the following activities: (a) in-situ infrared and optical imaging of each layer during processing of representative stainless steel geometries; (b) development of methodologies to analyze the data; and (c) application of ICME methodologies and other experiments to estimate the relevance of defects and microstructure to tensile properties.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
ICME and In-Situ Process Monitoring for Rapid Qualification of Components Made by Laser-based Powder Bed Additive Manufacturing Processes for Nuclear Structural Applications
Additive manufacturing (AM) through selective powder bed melting of successive layers is being considered as disruptive technology for rapid production of low-cost nuclear reactor internal components with complex geometries. However, there is a potential for non-uniform distribution of physical features, such as porosity or microstructural differences, due to variations of temperature across a component build. These heterogeneities make the deployment of traditional material qualification and non-destructive evaluation of AM components difficult. This research project explored the feasibility of using in-situ process monitoring methods and integrated computational materials engineering (ICME) principles as an alternate qualification methodology and approach. This project included six individual tasks: (i) design of artifacts relevant to the nuclear power industry, (ii) evaluation of laser processing and in-situ measurements, (iii) computational modeling, (iv) ex-situ microstructural characterization, (v) evaluation of scaling the methodology for large-scale structures, and (vi) development of a data package to codes and standards organizations. As a part of this research, laser powder bed fusion of metals (316L, Alloy 718, and Ti6Al4V) was explored. This document provides an overview of the research performed over three years and serves as the final U.S. Department of Energy report for this project. In the third year, the project focused on the following activities: (a) in-situ infrared and optical imaging of each layer during processing of representative stainless steel geometries; (b) development of methodologies to analyze the data; and (c) application of ICME methodologies and other experiments to estimate the relevance of defects and microstructure to tensile properties.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
ICME and In-Situ Process Monitoring for Rapid Qualification of Components Made by Laser-based Powder Bed Additive Manufacturing Processes for Nuclear Structural Applications
Additive manufacturing (AM) through selective powder bed melting of successive layers is being considered as disruptive technology for rapid production of low-cost nuclear reactor internal components with complex geometries. However, there is a potential for non-uniform distribution of physical features, such as porosity or microstructural differences, due to variations of temperature across a component build. These heterogeneities make the deployment of traditional material qualification and non-destructive evaluation of AM components difficult. This research project explored the feasibility of using in-situ process monitoring methods and integrated computational materials engineering (ICME) principles as an alternate qualification methodology and approach. This project included six individual tasks: (i) design of artifacts relevant to the nuclear power industry, (ii) evaluation of laser processing and in-situ measurements, (iii) computational modeling, (iv) ex-situ microstructural characterization, (v) evaluation of scaling the methodology for large-scale structures, and (vi) development of a data package to codes and standards organizations. As a part of this research, laser powder bed fusion of metals (316L, Alloy 718, and Ti6Al4V) was explored. This document provides an overview of the research performed over three years and serves as the final U.S. Department of Energy report for this project. In the third year, the project focused on the following activities: (a) in-situ infrared and optical imaging of each layer during processing of representative stainless steel geometries; (b) development of methodologies to analyze the data; and (c) application of ICME methodologies and other experiments to estimate the relevance of defects and microstructure to tensile properties.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗