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
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.
Astronomical integration code - fortran solution of space mechanics - 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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Explore the source record for details and available documents.
Fatigue is a key consideration in the design of a long-life Wind Turbine Generator (WTG) system. This paper discusses the fatigue aspects of the large Mod-1 horizontal-axis WTG design starting with the characterization of the environment and proceeding through the design. Major sources of fatigue loading are discussed and methods of limiting fatigue loading are described. NASTRAN finite element models are used to determine dynamic loading and internal cyclic stresses. Recent developments in determining the allowable fatigue stress consistent with present construction codes are discussed relative to their application to WTG structural design.
An experiment to evaluate the possibility of disseminating a time code via the SIRIO-1 satellite, by using the onboard VHF repeater is described. The precision in the synchronization of remote clocks was expected to be of the order of 0.1 to 1 ms. The RF carrier was in the VHF band, so that low cost receivers could be used and then a broader class of users could be served. An already existing repeater, even if not designed specifically for communications could be utilized; the operation of this repeater was not intended to affect any other function of the spacecraft (both the SHF repeater and the VHF telemetry link were active during the time code dissemination via the VHF transponder).
A local area network is provided for a plurality of autonomous computers which operate at different rates and under different protocols coupled by network bus adapters to a global bus. A host computer (HC) divides a message file to be transmitted into blocks, each with a header that includes a data type identifier and a trailer. The associated network bus adapter (NBA) then divides the data into packets, each with a header to which a transport header and trailer is added with frame type code which specifies one of three modes of addressing in the transmission of data, namely a physical address mode for computer to computer transmission using two bytes for source and destination addresses, a logical address mode and a data type mode. In the logical address mode, one of the two addressing bytes contains a logical channel number (LCN) established between the transmitting and one or more receiving computers. In the data type mode, one of the addressing bytes contains a code identifying the type of data.
An experimental study of low Reynolds number nozzle flow was performed. A brief comparison was made between some of the experimental performance data and performance predicted by a viscous flow code. The performance of 15, 20, and 25 deg conical nozzles, bell nozzles, and trumpet nozzles was evaluated with unheated nitrogen and hydrogen. The numerical analysis was applied to the conical nozzles only, using an existing viscous flow code that was based on a slender-channel approximation. Although the trumpet and 25 deg conical nozzles had slightly better performance at lower Reynolds numbers, it is unclear which nozzle is superior as all fell within the experimental error band. The numerical rssults were found to agree with experimental results for nitrogen and for some of the hydrogen data. Some code modification is recommended to improve confidence in the performance prediction.
ATIC(Advanced Thin Ionization Calorimeter) is a balloon borne experiment designed to measure the cosmic ray composition for elements from hydrogen to iron and their energy spectra from approx.50 GeV to near 100 TeV. It consists of a Si-matrix detector to determine the charge of a CR particle, a scintillator hodoscope for tracking, carbon interaction targets and a fully active BGO calorimeter. ATIC had its first 16-day flight from McMurdo, Antarctica from 28/12/2000 to 13/01/2000. The ATIC flight collected approximately 25 million events. To measure charge of primary particle in presence of radiation scattered back from the interaction and subsequent shower development in the calorimeter a charge detector must be a mosaic of small detector pads so that the pad containing the signal from the incident particle has no additional signal from albedo particles. Therefore the silicon matrix was built of 4480 individual silicon pads each 2 cm x 1.5 cm. The matrix consists of four planes of detectors and the active detector area, in these planes are partially overlapped to completely cover the aperture. The lateral and amplitude distributions of albedo signals in Si-matrix are analyzed for different primary nuclei and different energy deposits in BGO calorimeter. The greater part of albedo signals has Q near 1, where Q = square root of Amplitude(MIP). The albedo distribution exponentially decreases up to Q near 8. These high values are produced by slow protons and plans. There are also a small number of signals of Q > 8, mainly for heavy nucleus primaries. These signals are apparently generated by neutrons. The comparison of the experimental data and simulations with GEANT 3-21 code using QGSM generator for nucleus-nucleus interactions is presented.
We review evidence supporting the idea that the DNA sequence in genes containing non-coding regions is correlated, and that the correlation is remarkably long range--indeed, nucleotides thousands of base pairs distant are correlated. We do not find such a long-range correlation in the coding regions of the gene. We resolve the problem of the "non-stationarity" feature of the sequence of base pairs by applying a new algorithm called detrended fluctuation analysis (DFA). We address the claim of Voss that there is no difference in the statistical properties of coding and non-coding regions of DNA by systematically applying the DFA algorithm, as well as standard FFT analysis, to every DNA sequence (33301 coding and 29453 non-coding) in the entire GenBank database. Finally, we describe briefly some recent work showing that the non-coding sequences have certain statistical features in common with natural and artificial languages. Specifically, we adapt to DNA the Zipf approach to analyzing linguistic texts. These statistical properties of non-coding sequences support the possibility that non-coding regions of DNA may carry biological information.