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.
A study to demonstrate how the dynamics of galaxies may be investigated through the creation of galaxies within a computer model is presented. The numerical technique for simulating galaxies is shown to be both highly efficient and highly robust. Consideration is given to the anatomy of a galaxy, the gravitational N-body problem, numerical approaches to the N-body problem, use of the Poisson equation, and the symplectic integrator.
A large proportion of observed planetary systems contain several planets in a compact orbital configuration, and often harbor at least one close-in object. These systems are then most likely tidally evolving. We investigate how the effects of planet-planet interactions influence the tidal evolution of planets. We introduce for that purpose a new open-source addition to the Mercury N-body code, Mercury-T, which takes into account tides, general relativity and the effect of rotation-induced flattening in order to simulate the dynamical and tidal evolution of multi-planet systems. It uses a standard equilibrium tidal model, the constant time lag model. Besides, the evolution of the radius of several host bodies has been implemented (brown dwarfs, M-dwarfs of mass 0:1 M(sub ⨀), Sun-like stars, Jupiter). We validate the new code by comparing its output for one-planet systems to the secular equations results. We find that this code does respect the conservation of total angular momentum. We applied this new tool to the planetary system Kepler-62. We find that tides influence the stability of the system in some cases. We also show that while the four inner planets of the systems are likely to have slow rotation rates and small obliquities, the fifth planet could have a fast rotation rate and a high obliquity. This means that the two habitable zone planets of this system, Kepler-62e ad f are likely to have very different climate features, and this of course would influence their potential at hosting surface liquid water..
Simulations of shell-forming galaxy mergers are reported in which gas and stars are tracked separately to study the implications of observations that indicate fine structures in elliptical galaxies. The hybrid N-body/hydrodynamics code by Hernquist and Katz (1989) is employed neglecting the self-gravity of galaxies in tidal disruptions by massive primaries of small dwarf galaxies. The evolution of the galactic debris is examined showing that the gas and stars are effectively segregated with the stars forming sharp-edged features and the gas forming disks in the nucleus of the primary. Star formation is theorized to occur in the nuclei of shell galaxies due to shocks in the gas, and the added presence of a black hole can lead to an active galactic nucleus. The simulation is of interest to the interpretation of such observations as those of NGC1275 in which dormant elliptical galaxies are activated by mergers between large galaxies and gas-rich satellites.
Manual revisions are presented for the modified and expanded STEAP series. The STEAP 2 is composed of three independent but related programs: NOMAL for the generation of n-body nominal trajectories performing a number of deterministic guidance events; ERRAN for the linear error analysis and generalized covariance analysis along specific targeted trajectories; and SIMUL for testing the mathematical models used in the navigation and guidance process. The analytic manual provides general problem description, formulation, and solution and the detailed analysis of subroutines. The programmers' manual gives descriptions of the overall structure of the programs as well as the computational flow and analysis of the individual subroutines. The user's manual provides information on the input and output quantities of the programs. These are updates to N69-36472 and N69-36473.
Simulation of a close passage of two galaxies where the self-gravity of the disk stars is taken into account was performed with the aid of a large-scale N-body model (Hohl and Hockney, 1969; Hohl, 1972). The calculation was performed by perturbing an initially stable axisymmetric disk consisting of 100,000 stars by the passage of a companion galaxy having one-fourth the mass of the primary galaxy. Parameters appropriate for M 51 and its companion NGC 5195 were used. The results show the evolution of a weak two arm spiral. The effect of large random velocities on the formation of spiral structure was studied by plotting the spatial distribution of stars for various velocity intervals. The spiral structure is quite pronounced for stars with the lower velocity dispersion, but it can hardly be detected for stars with higher velocities. The tidal interaction causes only slight changes in the azimuthally averaged density distribution.
A series of three computer programs for the mathematical analysis of navigation and guidance of lunar and interplanetary trajectories was developed. All three programs require the integration of n-body trajectories for both interplanetary and lunar missions. The virutal mass technique is used in all three programs. The user's manual contains the information necessary to operate the programs. The input and output quantities of the programs are described. Sample cases are given and discussed.
Detailed maps of most elliptical galaxies reveal that, whereas the greatest part of their luminous mass originates from a smooth distribution with a surface brightness approximated by a de Vaucouleurs law, a small percentage of their light is contributed by low surface brightness distortions termed 'fine structures'. The sharp-edged features called 'shells' are successfully reproduced by merger and infall models involving accretion from less massive companions. In this context, dwarf spheroidal and compact disk galaxies are likely progenitors of these stellar phenomena. However, it is probable that the sources of shell-forming material also contain significant amounts of gas. This component may play an important role in constraining the formation and evolution of shell galaxies. To investigate the effects of the gaseous component, numerical simulations were performed to study the tidal disruption of dwarf galaxies containing both gas and stars by more massive primaries, and the evolution of the ensuing debris. The calculations were performed with a hybrid N-body/hydrodynamics code. Collisionless matter is evolved using a conventional N-body technique and gas is treated using smoothed particle hydrodynamics in which self-gravitating fluid elements are represented as particles evolving according to Lagrangian hydrodynamic equations. An isothermal equation of state is employed so the gas remains at a temperature 104 K. Owing to the large mass ratio between the primary and companion, the primary is modeled as a rigid potential and the self-gravity of both galaxies is neglected.
Observed ring galaxies appear to fall into two major types. The first tends to consist of isolated galaxies which display a smooth, apparently circular ring and a central nucleus. These have been variously classified as R(S) by de Vaucouleurs (1959) and as type O by Few and Madore (1986). The second class of ring galaxy nearly always has a close companion of comparable size (no less than about one tenth that of the ring galaxy). In these objects the ring is knotty in appearance, is usually elliptical, even when deprojected on the sky, and is often open on one side, having a 'horse shoe' or 'banana' shape. The nucleus does not usually appear at the center of the ring and is sometimes apparently absent, giving rise to an 'empty ring' galaxy. deVaucouleurs et al. (1976) designated this second type as RING, while Few and Madore (1986) have classified similar galaxies as P type. These galaxies have elevated far IR emission, bright HII regions, and blue spectral colors. The different environments of the two types or ring galaxy, together with their overall morphological and spectral differences suggest that the R(S)/O type are most probably the result of an instability that occurs in isolated galaxies, whereas the RING/P type appears to be the result of a recent collision between two roughly equal mass objects, at least one of which is a disk galaxy. Theys and Spiegel (1976) studied a sample of this latter type and identified three subclasses: RE: galaxies with crisp, empty rings; RN: galaxies like those of RE but with off-center nuclei; RK: galaxies having single dominant knots or condensations in the rings. A presentation of a preliminary understanding of the connections between these different observed forms in terms of parameters which are intrinsic to the galaxy system, such as time since collision and impact parameter, and in terms of our line of sight view is the purpose of this paper. Here we report results we have obtained from three dimensional computer simulations of collisions between equal mass galaxies, one of which is a rotating, disk galaxy containing both gas and stars and the other is an elliptical containing stars only. We have used a combined n-body/SPH program (see Balsara, 1990) to model fully self consistent models in which the halo mass is 2.5 times that of the disk and gas comprises ten percent of the disk mass.
We leverage the IllustrisTNG cosmological simulations to test how the large-scale cosmic web shapes galaxy morphology and to forecast the early cosmic web structure that the Nancy Grace Roman Space Telescope will reveal. In the hydrodynamic TNG50 and N-body TNG50-Dark runs, we reconstruct the cosmic web at redshifts z = 0, 0.5, 1, 2, 3, and 4 with the Monte Carlo Physarum Machine density estimator and the DisPerSE structure identification framework. We confirm that dark matter halos start out predominantly prolate (elongated), and their shapes are aligned with their nearest filaments; prolate galaxies retain strong shape alignment with their outer halos to later times. At z ≥ 1, the fraction of prolate (spheroidal) halos increases (decreases) toward lower stellar mass, higher redshift, and lower filament density. At z < 1, more spheroidal (oblate) stellar structures preferentially reside in higher-density (lower-density) filaments. We also find that higher-density filaments favor extended rotationally supported disks, whereas lower-density filaments more often host smaller dispersion-supported systems. Then, generating mock galaxy samples from TNG100 and TNG50, we predict the early cosmic web accessible to Roman. We find that the spectroscopic emission-line depth planned for the High-Latitude Wide-Area Survey (HLWAS) yields a highly incomplete galaxy sample that does not accurately trace the z = 1 cosmic web. A survey ≥2.5× deeper over a few square degrees would enable a proper reconstruction and reveal qualitatively correct filament–galaxy morphology relationships. Nevertheless, the planned HLWAS Deep field should still identify most galaxy overdensities; targeted deeper spectroscopy of these regions would efficiently and adequately map the early filamentary structure.
Galactic gas dynamics, spiral structure, and star formation are discussed in relation to N-body computational studies based on a cloud-particle model of the interstellar medium. On the small scale, the interstellar medium is seen as cloud-dominated and supernova-perturbed. It is noted that the cloud-particle model simulates cloud-cloud collisions, the formation of stellar associations, and supernova explosions as dominant local processes. On the large scale, in response to a spiral galactic gravitational field, global density waves and galactic shocks develop having large-scale characteristics similar to those found in continuum gas dynamical studies. Both the system of gas clouds and the system of young stellar associations forming from the clouds figure in the global spiral structure. However, with the attributes of neither assuming a continuum of gas (as in continuum gas dynamical studies) or requiring a prescribed equation of state (such as the isothermal condition), the cloud-particle picture retains much of the detail lost in earlier work. By detail is meant the small-scale features and structures so important in understanding the local, turbulent state of the interstellar medium as well as the degree of raggedness often seen to be superposed on the global spiral structure.
When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature. To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m = 50 cm 2 /g). When simulating these halos, we also vary the mass resolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses. The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10 6 particles per halo), we find that variations of order 10% in collapse time are still present. Here, the results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.
The first major evolution of NASA’s Space Launch System (SLS) will begin its flights starting in the mid-2020’s. This new configuration, called Block 1B, replaces the Interim Cryogenic Propulsion Stage with a larger Exploration Upper Stage (EUS). The additional capability provided by the new upper stage will allow SLS to send heavier payloads into deep space. One destination of interest to the SLS program is called a Near Rectilinear Halo Orbit (NRHO). This is a type of lunar orbit with multiple advantages for deep space exploration. These benefits include Earth/lunar access, low station-keeping requirements, and high communication potential with Earth. Therefore, it is a leading candidate for the proposed Lunar Orbital Platform-Gateway (LOP-G). This paper will provide a detailed assessment of the SLS Block 1B requirements and capabilities for sending payloads to an NRHO. Analysts at Marshall Space Flight Center are producing a multi-year mission availability scan for the SLS Block 1B configuration to a predefined NRHO orbit. The analysis produces an optimized trajectory for each day of the scan window. A maximized payload and minimized propellant requirement are determined for each day. All maneuvers from launch to the end of the Trans-Lunar-Injection (TLI) are being modeled as finite burns. Injections into an NRHO are being modeled as impulsive maneuvers. The payload element is arbitrary, but includes sufficient mass to represent a large habitat or propulsion module. The resulting parameters of payload capability, delta-v requirements, and launch windows length vary over the course of the scan. Many launch days in the scan are eliminated in post-processing as they violate mission constraints such as payload mass and propellant usage to insert into an NRHO. Based off previous one-year scan results for the SLS Block 1B Design Analysis Cycle 2 (DAC-2), it is expected that there will be one to three days per week where the payload is able to insert into an NRHO within SLS constraints objectives. This scan provides results for longer than one-year, allowing analysts to better understand the launch availability and energy requirements of SLS Block 1B over time. The in-space mission design and scans utilize Copernicus, an n-body trajectory optimization tool originally developed out of the University of Texas at Austin with further development at Johnson Space Center in Houston, TX. To seed the in-space trajectory, Copernicus uses a plugin to call a database of SLS ascent trajectories optimized in the Program to Simulate Optimized Trajectories II (POST2). The ascent trajectories are developed using a framework that parametrizes payload mass and LEO inclination, and targets a 100 nmi (nautical mile) altitude circular parking orbit. Though this analysis is specific to the SLS program, it will provide a summary of mission design benefits and constraints associated with generic NRHO access, and may be applied to other programs or concepts that will utilize this orbit.
Theoretical and computational studies were carried out on galactic and planetary disks. With the goal of addressing important open questions centered on galactic structure, the cloudy interstellar medium, giant molecular clouds, and star formation in galactic disks and the collisional dynamics of perturbed particulate matter in planetary disks, focus was largely on the self-gravitational effects, dissipative effects, and collisional dynamics of cloud-particle disks. N-body, 'cloud-particle' computational algorithms were developed for the purpose of isolating the role of gaseous self gravity from the roles of other dominant physical mechanisms and dynamical processes, e.g. the collisional dynamics and dissipative processes. The efforts focused largely on galactic disks show that observational constraints provide stringent tests of the numerical simulation techniques developed. Self gravitational effects of the galactic interstellar medium's gas clouds are included by means of Fourier Transform technique.
Complex mechanisms regulate the cellular distribution of cholesterol, a critical component of eukaryote membranes involved in regulation of membrane protein functions directly and through the physiochemical properties of membranes. StarD4, a member of the steroidogenic acute regulator-related lipid-transfer (StART) domain (StARD)-containing protein family, is a highly efficient sterol-specific transfer protein involved in cholesterol homeostasis. Its mechanism of cargo loading and release remains unknown despite recent insights into the key role of phosphatidylinositol phosphates in modulating its interactions with target membranes. We have used large-scale atomistic Molecular dynamics (MD) simulations to study how the dynamics of cholesterol bound to the StarD4 protein can affect interaction with target membranes, and cargo delivery. We identify the two major cholesterol (CHL) binding modes in the hydrophobic pocket of StarD4, one near S136&S147 (the Ser-mode), and another closer to the putative release gate located near W171, R92&Y117 (the Trp-mode). We show that conformational changes of StarD4 associated directly with the transition between these binding modes facilitate the opening of the gate. To understand the dynamics of this connection we apply a machine-learning algorithm for the detection of rare events in MD trajectories (RED), which reveals the structural motifs involved in the opening of a front gate and a back corridor in the StarD4 structure occurring together with the spontaneous transition of CHL from the Ser-mode of binding to the Trp-mode. Further analysis of MD trajectory data with the information-theory based NbIT method reveals the allosteric network connecting the CHL binding site to the functionally important structural components of the gate and corridor. Mutations of residues in the allosteric network are shown to affect the performance of the allosteric connection. These findings outline an allosteric mechanism which prepares the CHL-bound StarD4 to release and deliver the cargo when it is bound to the target membrane.
An initially stationary and stable axisymmetric disk of stars with a mass spectrum of stars is used to determine the collisional relaxation time. The relaxation time as determined from the rate of energy equipartition was found to be 560 rotation periods for stars with 10 times the mean star mass and was 1700 rotation periods for stars with 0.55 times the mean star mass. These times are in general agreement with theoretical predictions for the relaxation time of the two mass groups. The results show that the model used for the large-scale gravitational N-body calculations is indeed 'collisionless.'
A first attempt is made to calculate the properties of the matter distribution in a universe filled with overlapping bubbles produced by multiple explosions. Each spherical shell follows the cosmological Sedov-Taylor solution until it encounters another shell. Thereafter, mergers are allowed to occur in pairs on the basis of N-body results. At the final epoch, the matrix of overlapping shells is populated with 'galaxies' and the properties of slices through the numerically constructed cube compare well with CfA survey results for specified initial conditions. A statistic is found which measures the distance distribution from uniformly distributed points to the nearest galaxies on the projected plane which appears to provide a good measure of the bubbly character of the galaxy distribution. In a quantitative analysis of the CfA 'slice of the universe', a very good match is found between simulation and the real data for final average bubble radii of (13.5 + or - 1.5)/h Mpc with formal filling factor 1.0-1.5 or actual filling factor of 65-80 percent.
How can one fully harness the power of physics encoded in relativistic N-body phase space? Topologically, phase space is isomorphic to the product space of a simplex and a hypersphere and can be equipped with explicit coordinates and a Riemannian metric. This natural structure that scaffolds the space on which all collider physics events live opens up new directions for machine learning applications and implementation. Here we present a detailed construction of the phase space manifold and its differential line element, identifying particle ordering prescriptions that ensure that the metric satisfies necessary properties. We apply the phase space metric to several binary classification tasks, including discrimination of high-multiplicity resonance decays or boosted hadronic decays of electroweak bosons from QCD processes, and demonstrate powerful performance on simulated data. Our work demonstrates the many benefits of promoting phase space from merely a background on which calculations take place to being geometrically entwined with a theory’s dynamics.
Crater II (CraII), a large and low-density dwarf spheroidal galaxy, has unusual observed properties that are difficult to reproduce in cold dark matter simulations. Ongoing tidal disruption may help explain the discrepancies, as evidenced by the recent discovery of tidal tails. Here we present metallicity-sensitive narrowband photometry of the Ca II H and K lines from the Dark Energy Camera, covering 128 deg$^{2}$ across the center and identified tidal tails of CraII as part of the Mapping the Ancient Galaxy in CaHK survey. Our combined photometric metallicity, color–magnitude, proper motion, and parallax selections identify 162 CraII candidates. Of these, 37 candidates are located in the tidal tails, which extend at least 7° (∼95 kpc) from the center of CraII, suggesting it has lost ≳25% of its initial stellar mass. We confirm low contamination rates with dedicated control fields and highlight the extremely low surface brightness stellar features that can be uncovered with CaHK data, as faint as ∼36 mag arcsec$^{−2}$. We also make the first detection of a metallicity gradient (−0.34 ± 0.17 dex deg$^{−1}$) in the center of the galaxy and infer a stream width of w ∼ 0.°8, roughly 50% larger than the CraII half-light radius. The detection of candidates in the most distant CraII pointings from its center implies that the tidal tails extend beyond our footprint. We compare the CraII stream to N-body models with “cored” and “cuspy” dark matter halo progenitors, determining that CraII’s density profile is still ambiguous and warrants further modeling.