SEARCH · Search NASA
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.
Astronomy is an area of applied physics in which unusually beautiful objects challenge the imagination to explain observed phenomena in terms of known laws of physics. It is a field that has stimulated the development of physical laws and of mathematical and computational methods. Current computational applications are discussed in terms of stellar and galactic evolution, galactic dynamics, and particle motions.
The solar nebula, from which the planets in our solar system formed, featured a disk of gas and dust grains in rapid, differential rotation, and at some stage was likely to have been unstable to thermal convection. This situation is suspected by many to lead to significant turbulent Reynolds stress production and angular momentum transport in such systems, and estimates of transport rates have been attempted from unsubstantiated phenomenological models. In order to determine the circumstances and physical conditions under which our own planetary system formed and to explain recent observations of young stellar systems, it is necessary to develop realistic models of heat and angular momentum transport for such flows. Developing an understanding of complicated flows featuring thermal convection, rotation, and shear is also of wide interest in stellar astrophysics and in planetary and terrestrial atmospheric studies. The ultimate objective is to develop workable models based on the numerical simulations for constructing global solar nebula models; viz., relatively simple prescriptions for heat and angular momentum fluxes from given system parameters (e.g., ratios of rotation, shear, and convective lapse rates) are characterized, quantified, and developed. Toward this end, our program has been to attempt to understand the behavior of the direct numerical simulations of Boussinesq convection, which, despite the complexity of the results, is still an overly simplified approximation to the real system and should be more amenable to analysis. These results are also intended to be tested against turbulence models, especially those designed for atmospheric boundary layers, and may provide a basis for subgrid-scale models. In order to make the numerical simulations more realistic with regard to the solar nebula problem, a fully compressible code that will allow incorporation of large density stratifications and realistic thermodynamic and radiative properties is developed. In order to explore the properties of these flows at the very high values of Re found in natural systems and the very low values of Pr found in most astrophysical contexts, we will need to employ large-eddy simulations for which we want to determine the most appropriate subgrid-scale model to incorporate.
The Final Technical Report on a number of related investigations on the physics of high energy emission from active galactic nuclei, such as Seyfert galaxies and quasi-stellar objects is presented. The chief conclusions of the work are briefly described, and citations to the papers supported by this grant and published in the refereed scientific literature are provided. Areas of research included: 'warm' galaxies observed in x rays; x ray/infrared correlations in galaxies; the contribution of active galaxies to the cosmic x ray background radiation; and an unusual x ray emitting starburst galaxy.
We present a rest-frame ultraviolet morphological analysis of 108 z = 2.1 Lyman Alpha Emitters (LAEs) in the Extended Chandra Deep Field South (ECDF-S) and compare it to a similar sample of 171 LAEs at z = 3.1 . Using Hubble Space Telescope (HST) images taken as part of the Galaxy Evolution From Morphology and SEDs survey, Great Observatories Origins Deep Survey, and Hubble Ultradeep Field surveys, we measure the size and photometric component distributions, where photo- metric components are defined as distinct clumps of UV-continuum emission. At both redshifts, the majority of LAEs have observed half-light radii < 2 kpc, but the median half-light radius rises from 0.97 kpc at z = 3.1 to 1.41 kpc at z = 2.1. A similar evolution is seen in the sizes of individual rest-UV components, but there is no evidence for evolution in the number of mUlti-component systems. In the z = 2.1 LAE sample, we see clear correlations between the LAE size and other physical properties derived from its SED. LAEs are found to be larger for galaxies with larger stellar mass, larger star formation rate, and larger dust obscuration, but there is no evidence for a trend between equivalent width and half-light radius at either redshift. The presence of these correlations suggests that a wide range of objects are being selected by LAE surveys at that redshift, including a significant fraction of objects for which a massive and moderately extended population of old stars underlies the young starburst giving rise to the Lya emission.
Approximate analytic expressions are derived for resonance-line wing diagnostics, accounting for frequency redistribution effects, for homogeneous slabs, and slabs with a constant Planck function gradient. Resonance-line emission profiles from a simplified conceptual standpoint are described in order to elucidate the basic physical parameters of the line-forming layers prior to the performance of detailed numerical calculations. An approximate analytic expression is derived for the dependence on stellar surface gravity of the location of the Ca II and Mg II resonance-line profile peaks. An approximate radiative transfer equation using generalized second-order escape probabilities, applicable even in the presence of nearly coherent scattering in the damping wings of resonance lines, is derived. Approximate analytic solutions that can be applied in special regimes and achieve good agreement with accurate numerical results are found.
Dynamic and adaptive mesh refinement is pivotal in high-resolution, multi-physics, multi-model simulations, necessitating precise physics resolution in localized areas across expansive domains. Today’s supercomputers’ extreme heterogeneity presents a significant challenge for dynamically adaptive codes, highlighting the importance of achieving performance portability at scale. Our research focuses on astrophysical simulations, particularly stellar mergers, to elucidate early universe dynamics. Here, we present Octo-Tiger, leveraging Kokkos, HPX, and SIMD for portable performance at scale in complex, massively parallel adaptive multi-physics simulations. Octo-Tiger supports diverse processors, accelerators, and network backends. Experiments demonstrate exceptional scalability across several heterogeneous supercomputers including Perlmutter, Frontier, and Fugaku, encompassing major GPU architectures and x86, ARM, and RISC-V CPUs. Parallel efficiency of 47.59% (110,080 cores and 6880 hybrid A100 GPUs) on a full-system run on Perlmutter (26% HPCG peak performance) and 51.37% (using 32,768 cores and 2048 MI250X) on Frontier are achieved.
Numerical computation of the ideal magnetohydrodynamic (MHD) equilibrium magnetic field is at the base of stellarator optimisation and provides the starting point for solving more sophisticated partial differential equations like transport or turbulence models. Conventional approaches solve for a single stationary point of the ideal MHD equations, which is fully defined by three invariants and the numerical scheme employed by the solver. We present the first numerical approach that can solve for a continuous distribution of equilibria with fixed boundary and rotational transform, varying only the pressure invariant. This approach minimises the force residual by optimising parameters of multilayer perceptrons that map from a scalar pressure multiplier to the Fourier Zernike basis as implemented in the modern stellarator equilibrium solver DESC.
The 12 C + 12 C sub-barrier fusion cross section is calculated within the framework of a time-dependent Hartree-Fock-based classical model using the Feynman path-integral method. The modified astrophysical S* factor is compared to direct and indirect experimental results. A good agreement with the direct data is found. In the lower-energy region where recent analyses of experimental data obtained with the Trojan horse method (THM) lead to contrasting results, the model predicts a nonresonant S* factor half-way between those results. Low-energy resonances revealed in the THM data are added to the calculation, and the relative reaction rate in the Gamow region is calculated. In particular, including 0 + resonances result in some agreement with the THM data. Here, the role of different resonances is discussed in detail, and their influence on the reaction rate at temperatures relevant to stellar evolution is investigated.
We explore the contribution of the Gaia Sausage to the stellar halo of the Milky Way by making use of a Gaussian mixture model (GMM) and applying it to halo star samples of Large Sky Area Multi-Object Fiber Spectroscopic Telescope K giants, Sloan Extension for Galactic Understanding and Exploration K giants, and Sloan Digital Sky Survey blue horizontal branch stars. The GMM divides the stellar halo into two parts, of which one represents a more metal-rich and highly radially biased component associated with an ancient, head-on collision referred to as the Gaia Sausage, and the other one is a more metal-poor and isotropic halo. A symmetric bimodal Gaussian is used to describe the distribution of spherical velocity of the Gaia Sausage, and we find that the mean absolute radial velocity of the two lobes decreases with the Galactocentric radius. We find that the Gaia Sausage contributes about 41%–74% of the inner (Galactocentric radius r gc < 30 kpc) stellar halo. The fraction of stars of the Gaia Sausage starts to decline beyond r gc ~ 25–30 kpc, and the outer halo is found to be significantly less influenced by the Gaia Sausage than the inner halo. After the removal of halo substructures found by integrals of motion, the contribution of the Gaia Sausage falls slightly within r gc ~ 25 kpc but is still as high as 30%–63%. Finally, we select several possible Sausage-related substructures consisting of stars on highly eccentric orbits. The GMM/Sausage component agrees well with the selected substructure stars in their chemodynamical properties, which increases our confidence in the reliability of the GMM fits.
The usage of permanent magnets to shape the confining magnetic field of a stellarator has the potential to reduce or eliminate the need for non-planar coils. As a proof-of-concept for this idea, here we have developed a procedure for designing an array of cubic permanent magnets that works in tandem with a set of toroidal-field (TF) coils to confine a stellarator plasma. All of the magnets in the design are constrained to have identical geometry and one of three polarization types in order to simplify fabrication while still producing sufficient field accuracy. We present some of the key steps leading to the design, including the geometric arrangement of the magnets around the device, the procedure for optimizing the polarizations according to the three allowable magnet types, and the choice of magnet types to be used. We apply these methods to design an array of rare-Earth permanent magnets that can be paired with a set of planar TF coils to confine a quasi-axisymmetric plasma with a toroidal magnetic field strength of about 0.5 T on axis.
We present a 6D map of the Orphan–Chenab (OC) stream by combining the data from Southern Stellar Stream Spectroscopic Survey (S 5 ) and Gaia. We reconstruct the proper motion, radial velocity, distance, on-sky track, and stellar density along the stream with spline models. The stream has a total luminosity of M V = -8.2 and metallicity of [Fe/H] = -1.9, similar to classical Milky Way (MW) satellites like Draco. The stream shows drastic changes in its physical width varying from 200 pc to 1kpc, but a constant line-of-sight velocity dispersion of 5 km s -1 . Despite the large apparent variation in the stellar number density along the stream, the flow rate of stars along the stream is remarkably constant. We model the 6D stream track by a Lagrange-point stripping method with a flexible MW potential in the presence of a moving extended Large Magellanic Cloud (LMC). This allows us to constrain the mass profile of the MW within the distance range 15.6 < r < 55∙5 kpc, with the best measured enclosed mass of (2.85 ± 0.1) X 10 11 M ⊙ within 32.4 kpc. Our stream measurements are highly sensitive to the LMC mass profile with the most precise measurement of its enclosed mass made at 32.8 kpc, 7.02 ± 0.9) X 10 11 M ⊙ . We also detect that the LMC dark matter halo extends to at least 53 kpc. The fitting of the OC stream allows us to constrain the past LMC trajectory and the degree of dynamical friction it experienced. We demonstrate that the stars in the OC stream show large energy and angular momentum spreads caused by LMC perturbation.
The thermal infrared region contains molecular bands of many of the most important species in gaseous astronomical sources. True shapes and frequencies of emission and absorption spectral lines from these constituents of planetary and stellar atmospheres contain unique information on local temperature and abundance distribution, non-thermal effects, composition, local dynamics and winds. Heterodyne spectroscopy in the thermal infrared can remotely measure true line shapes in relatively cool and thin regions and enable the retrieval of detailed information about local physics and chemistry. The concept and techniques for heterodyne detection will be discussed including examples of thermal infrared photomixers and instrumentation used in studies of several astronomical sources. Use of heterodyne detection to study non-LTE phenomena, planetary aurora, minor planetary species and gas velocities (winds) will be discussed. A discussion of future technological developments and relation to space flight missions will be addressed.
The C, N, and Mg isotopic compositions in graphite and SiC grains of carbonaceous chondrites can be explained by nuclear processes in massive O,B stars of second generation passed a stage of WR star with intensive stellar wind, where grain condensation had taken place. The interstellar graphite and SiC grains with anomalous isotopic compositions of C, N, Ne, and Si and other elements of nucleosynthetic origin, found in non-equilibrated chondrites, are most suitable for determination of astro-physical objects, where nucleosynthesis had taken place. Two objects were examined: (1) massive O,B stars of second generation passed a stage of WR star with intensive stellar wind (O,B-WR model) and (2) low-mass stars (1 less than or equal to M/solar mass less than or equal to 3) during thermally pulsing asymptotic giant branch phase (TP-AGB model).
The stability of hybrid stars with first-order phase transitions as determined by calculating fundamental radial oscillation modes is known to differ from the predictions of the widely used Bardeen-Thorne-Meltzer criterion. We consider the effects of out-of-chemical-equilibrium physics on the radial modes and hence stability of these objects. For a barotropic equation of state, this is done by allowing the adiabatic sound speed to differ from the equilibrium sound speed. We show that doing so extends the stable branches of stellar models, allowing stars with rapid phase transitions to support stable higher-order stellar multiplets similarly to stars with multiple slow phase transitions. We also derive a new junction condition to impose on the oscillation modes at the phase transition. Termed the reactive condition, it is physically motivated, consistent with the generalized junction conditions between two phases, and has the common rapid and slow conditions as limiting cases. Unlike the two common cases, it can only be applied to nonbarotropic stars. Here, we apply this junction condition to hybrid stellar models generated using a two-phase equation of state consisting of nuclear matter with unpaired quark matter at high densities joined by a first-order phase transition and show that, like in the slow limiting case, stars that are classically unstable are stabilized by a finite chemical reaction speed.
Coil complexity is a critical consideration in stellarator design. The traditional two-step optimization approach, in which the plasma boundary is optimized for physics properties and the coils are subsequently optimized to be consistent with this boundary, can result in plasma shapes which cannot be produced with sufficiently simple coils. To address this challenge, we propose a method to incorporate considerations of coil complexity in the optimization of the plasma boundary. Coil complexity metrics are computed from the current potential solution obtained with the REGCOIL code (Landreman, Nucl. Fusion , vol. 57, 2017, 046003). While such metrics have previously been included in derivative-free fixed-boundary optimization (Drevlak et al. , Nucl. Fusion , vol. 59, 2018, 016010), we compute the local sensitivity of these metrics with respect to perturbations of the plasma boundary using the shape gradient (Landreman & Paul, Nucl. Fusion , vol. 58, 2018, 076023). We extend REGCOIL to compute derivatives of these metrics with respect to parameters describing the plasma boundary. In keeping with previous research on winding surface optimization (Paul et al. , Nucl. Fusion , vol. 58, 2018, 076015), the shape derivatives are computed with a discrete adjoint method. In contrast with the previous work, derivatives are computed with respect to the plasma surface parameters rather than the winding surface parameters. To further reduce the resolution required to compute the shape gradient, we present a more efficient representation of the plasma surface which uses a single Fourier series to describe the radial distance from a coordinate axis and a spectrally condensed poloidal angle. This representation is advantageous over the standard cylindrical representation used in the VMEC code (Hirshman & Whitson, Phys. Fluids , vol. 26, 1983, pp. 3553–3568), as it provides a uniquely defined poloidal angle, eliminating a null space in the optimization of the plasma surface. In comparison with previous spectral condensation methods (Hirshman & Breslau, Phys. Plasmas , vol. 5, 1998, p. 2664), the modified poloidal angle is obtained algebraically rather than through the solution of a nonlinear optimization problem. The resulting shape gradient highlights features of the plasma boundary that are consistent with simple coils and can be used to couple coil and fixed-boundary optimization.
The formation of the Galactic disc is an enthusiastically debated issue. Numerous studies and models seek to identify the dominant physical process(es) that shaped its observed properties; for example, satellite accretion, starburst, quenching, gas infall, and stellar radial migration. Taking advantage of the improved coverage of the inner Milky Way provided by the SDSS DR16 APOGEE catalogue and of the ages published in the APOGEE-AstroNN Value Added Catalogue (VAC), we examined the radial evolution of the chemical and age properties of the Galactic stellar disc with the aim of better constraining its formation. Using a sample of 199 307 giant stars with precise APOGEE abundances and APOGEE-AstroNN ages, selected in a ±2 kpc layer around the galactic plane, we assessed the dependency as a function of guiding radius of (i) the median metallicity, (ii) the ridge lines of the [Fe/H] – [Mg/Fe] and age–[Mg/Fe] distributions, and (iii) the age distribution function (ADF). The giant star sample allows us to probe the radial behaviour of the Galactic disc from R g = 0 to 14–16 kpc. The thick disc [Fe/H] – [Mg/Fe] ridge lines follow closely grouped parallel paths, supporting the idea that the thick disc did form from a well-mixed medium. However, the ridge lines present a small drift in [Mg/Fe], which decreases with increasing guiding radius. At sub-solar metallicity, the intermediate and outer thin disc [Fe/H] – [Mg/Fe] ridge lines follow parallel sequences shifted to lower metallicity as the guiding radius increases. We interpret this pattern as the signature of a dilution of the interstellar medium from R g ~ 6 kpc to the outskirts of the disc, which occurred before the onset of the thin disc formation. The APOGEE-AstroNN VAC provides stellar ages for statistically significant samples of thin disc stars from the Galactic centre up to R g ~ 14 kpc. An important result provided by this dataset is that the thin disc presents evidence of an inside-out formation up to R g ~ 10 – 12 kpc. Moreover, about ~7 Gyr ago, the [Mg/Fe] ratio in the outer thin disc (R g > 10 kpc) was higher by about ~0.03–0.05 dex than in the more internal regions of the thin disc. This could be the fossil record of a pollution of the outer disc gas reservoir by the thick disc during its starburst phase.
ABSTRACT Primordial black holes may encode the conditions of the early Universe, and may even constitute a significant fraction of cosmological dark matter. Their existence has yet to be established. However, black holes with masses below ${\sim}{1}{\, \mathrm{M}_\odot }$ cannot form as an endpoint of stellar evolution, so the detection of even one such object would be a smoking gun for new physics, and would constitute evidence that at least a fraction of the dark matter consists of primordial black holes. Gravitational wave detectors are capable of making a definitive discovery of this kind by detecting mergers of light black holes. But since the merger rate depends strongly on the shape of the black hole mass function, it is difficult to determine the potential for discovery or constraint as a function of the overall abundance of black holes. Here, we directly maximize and minimize the merger rate to connect observational results to the actual abundance of observable objects. We show that LIGO can discover mergers of light primordial black holes within the next decade even if such black holes constitute only a very small fraction of dark matter. A single merger event involving such an object would (i) provide conclusive evidence of new physics, (ii) establish the nature of some fraction of dark matter, and (iii) probe cosmological history at scales far beyond those observable today.
A generalized formulation is given for treating partial redistribution effects in transfer problems in resonance lines with common upper states. The formulation allows explicitly for the possibility that several spectral lines may arise in transitions from a given upper level to several sharp lower levels, including, for example, the ground state and metastable states. Line profiles for the Ca II H and K lines have been calculated, accounting for the partial frequency coherence of scattered photons. These profiles are compared with calculations made with identical atomic and atmospheric models but assuming complete redistribution. Very significant differences between the profiles obtained using these two different physical descriptions of the scattering process are found, and it is now apparent that the assumption of complete redistribution is a serious oversimplification of the actual physical situation. The results question the validity of equating brightness temperatures observed at K(sub 1) in stellar spectra with minimum temperatures in stellar chromospheres; it appears likely that such a procedure will systematically underestimate the value of T-min.