Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multidimensional”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

TWO-DIMENSIONAL CORE-COLLAPSE SUPERNOVA SIMULATIONS WITH THE ISOTROPIC DIFFUSION SOURCE APPROXIMATION FOR NEUTRINO TRANSPORT

The neutrino mechanism of core-collapse supernova is investigated via non-relativistic, two-dimensional (2D), neutrino radiation–hydrodynamic simulations. For the transport of electron flavor neutrinos, we use the interaction rates defined by Bruenn and the isotropic diffusion source approximation (IDSA) scheme, which decomposes the transported particles into trapped-particle and streaming-particle components. Heavy neutrinos are described by a leakage scheme. Unlike the “ray-by-ray” approach in some other multidimensional supernova models, we use cylindrical coordinates and solve the trapped-particle component in multiple dimensions, improving the proto-neutron star resolution and the neutrino transport in angular and temporal directions. We provide an IDSA verification by performing one-dimensional (1D) and 2D simulations with 15 and 20 M{sub ⊙} progenitors from Woosley et al. and discuss the difference between our IDSA results and those existing in the literature. Additionally, we perform Newtonian 1D and 2D simulations from prebounce core collapse to several hundred milliseconds postbounce with 11, 15, 21, and 27 M{sub ⊙} progenitors from Woosley et al. with the HS(DD2) equation of state. General-relativistic effects are neglected. We obtain robust explosions with diagnostic energies E{sub dia} ≳ 0.1–0.5 B (1 B ≡ 10{sup 51} erg) for all considered 2D models within approximately 100–300 ms after bounce and find that explosions are mostly dominated by the neutrino-driven convection, although standing accretion shock instabilities are observed as well. We also find that the level of electron deleptonization during collapse dramatically affects the postbounce evolution, e.g., the neglect of neutrino–electron scattering during collapse will lead to a stronger explosion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

PROBING THE EXPANSION HISTORY OF THE UNIVERSE BY MODEL-INDEPENDENT RECONSTRUCTION FROM SUPERNOVAE AND GAMMA-RAY BURST MEASUREMENTS

To probe the late evolution history of the universe, we adopt two kinds of optimal basis systems. One of them is constructed by performing the principle component analysis, and the other is built by taking the multidimensional scaling approach. Cosmological observables such as the luminosity distance can be decomposed into these basis systems. These basis systems are optimized for different kinds of cosmological models that are based on different physical assumptions, even for a mixture model of them. Therefore, the so-called feature space that is projected from the basis systems is cosmological model independent, and it provides a parameterization for studying and reconstructing the Hubble expansion rate from the supernova luminosity distance and even gamma-ray burst (GRB) data with self-calibration. The circular problem when using GRBs as cosmological candles is naturally eliminated in this procedure. By using the Levenberg–Marquardt technique and the Markov Chain Monte Carlo method, we perform an observational constraint on this kind of parameterization. The data we used include the “joint light-curve analysis” data set that consists of 740 Type Ia supernovae and 109 long GRBs with the well-known Amati relation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Type Ia Supernova Explosions From Hybrid Carbon–Oxygen–Neon White Dwarf Progenitors

Motivated by recent results in stellar evolution that predict the existence of hybrid white dwarf (WD) stars with a C–O core inside an O–Ne shell, we simulate thermonuclear (Type Ia) supernovae from these hybrid progenitors. We use the FLASH code to perform multidimensional simulations in the deflagration-to-detonation transition (DDT) explosion paradigm. Our hybrid progenitor models were produced with the MESA stellar evolution code and include the effects of the Urca process, and we map the progenitor model to the FLASH grid. We performed a suite of DDT simulations over a range of ignition conditions consistent with the progenitor's thermal and convective structure assuming multiple ignition points. To compare the results from these hybrid WD stars to previous results from C–O WDs, we construct a set of C–O WD models with similar properties and similarly simulate a suite of explosions. We find that despite significant variability within each suite, trends distinguishing the explosions are apparent in their ${}^{56}\mathrm{Ni}$ yields and the kinetic properties of the ejecta. We compare our results with other recent work that studies explosions from these hybrid progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

An Updated Visual Orbit of the Directly Imaged Exoplanet 51 Eridani b and Prospects for a Dynamical Mass Measurement with Gaia

We present a revision to the visual orbit of the young, directly imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric (e=0.53{sub -0.13}{sup +0.09}) orbit at an intermediate inclination (i=136{sub -11}{sup +10} °), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semimajor axis of 11.1{sub -1.3}{sup +4.2} au and a period of 28.1{sub -4.9}{sup +17.2} yr, assuming a mass of 1.75 M{sub ⊙} for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigate the potential of using the absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51 Eri derived from a comparison of the Hipparcos and Gaia catalogs was found to be inconsistent at the 2σ–3σ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs and the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using Gaia astrometry alone for a dynamical mass measurement of the planet by simulating Gaia measurements of the motion of the photocenter of the system over the course of the extended 8 yr mission. We find that such a measurement is only possible (>98% probability) given the most optimistic predictions for the Gaia scan astrometric uncertainties for bright stars and a high mass for the planet (≳3.6 M {sub Jup}).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Significant Increase in Detection of High-resolution Emission Spectra Using a Three-dimensional Atmospheric Model of a Hot Jupiter

High-resolution spectroscopy has opened the way for new, detailed study of exoplanet atmospheres. There is evidence that this technique can be sensitive to the complex, three-dimensional (3D) atmospheric structure of these planets. In this work, we perform cross-correlation analysis of high-resolution (R ∼ 100,000) CRIRES/VLT emission spectra of the hot Jupiter HD 209458b. We generate template emission spectra from a 3D atmospheric circulation model of the planet, accounting for temperature structure and atmospheric motions—winds and planetary rotation—missed by spectra calculated from one-dimensional models. In this first-of-its-kind analysis, we find that using template spectra generated from a 3D model produces a more significant detection (6.9σ) of the planet’s signal than any of the hundreds of one-dimensional models we tested (maximum of 5.1σ). We recover the planet’s thermal emission, its orbital motion, and the presence of CO in its atmosphere at high significance. Additionally, we analyzed the relative influences of 3D temperature and chemical structures in this improved detection, including the contributions from CO and H{sub 2}O, as well as the role of atmospheric Doppler signatures from winds and rotation. This work shows that the hot Jupiter’s 3D atmospheric structure has a first-order influence on its emission spectra at high resolution and motivates the use of multidimensional atmospheric models in high-resolution spectral analysis.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Bayesian Cross-matching of High Proper-motion Stars in Gaia DR2 and Photometric Metallicities for ∼1.7 million K and M Dwarfs

We present a Bayesian method to cross-match 5,827,988 high proper-motion Gaia sources (μ > 40 mas yr{sup −1}) to various photometric surveys: Two Micron All Sky Survey, AllWISE data release from the Wide-field Infrared Explorer (WISE) mission, Galaxy Evolution Explorer, Radial Velocity Experiment, Sloan Digital Sky Survey, and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). To efficiently associate these objects across catalogs, we develop a technique that compares the multidimensional distribution of all sources in the vicinity of each Gaia star to a reference distribution of random field stars obtained by extracting all sources in a region on the sky displaced 2′. This offset preserves the local field stellar density and magnitude distribution, allowing us to characterize the frequency of chance alignments. The resulting catalog with Bayesian probabilities >95% has a marginally higher match rate than current internal Gaia data release 2 (DR2) matches for most catalogs. However, a significant improvement is found with Pan-STARRS, where ∼99.8% of the sample within the Pan-STARRS footprint is recovered, as compared to a low ∼20.8% in Gaia DR2. Using these results, we train a Gaussian process regressor to calibrate two photometric metallicity relationships. For dwarfs of 3500 < T{sub eff} < 5280 K, we use metallicity values of 4378 stars from the Apache Point Observatory Galactic Evolution Experiment and Hejazi et al. to calibrate the relationship, producing results with a 1σ precision of 0.12 dex and few systematic errors. We then indirectly infer the metallicity of 4018 stars with 2850 < T{sub eff} < 3500 K, which are wide companions of primaries whose metallicities are estimated with our first regressor, to produce a relationship with a 1σ precision of 0.21 dex and significant systematic errors. Additional work is needed to better remove unresolved binaries from this second sample to reduce these systematic errors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Constraining the Single-degenerate Channel of Type Ia Supernovae with Stable Iron-group Elements in SNR 3C 397

Recent Suzaku X-ray spectra of supernova remnant (SNR) 3C 397 indicate enhanced stable iron group element abundances of Ni, Mn, Cr, and Fe. Seeking to address key questions about the progenitor and explosion mechanism of 3C 397, we compute nucleosynthetic yields from a suite of multidimensional hydrodynamics models in the near-Chandrasekhar-mass, single-degenerate paradigm for Type Ia supernovae (SNe Ia). Varying the progenitor white dwarf (WD) internal structure, composition, ignition, and explosion mechanism, we find that the best match to the observed iron peak elements of 3C 397 are dense (central density ≥6 × 10{sup 9} g cm{sup −3}), low-carbon WDs that undergo a weak, centrally ignited deflagration, followed by a subsequent detonation. The amount of {sup 56}Ni produced is consistent with a normal or bright normal SNe Ia. A pure deflagration of a centrally ignited, low central density (≃2 × 10{sup 9} g cm{sup −3}) progenitor WD, frequently considered in the literature, is also found to produce good agreement with 3C 397 nucleosynthetic yields, but leads to a subluminous SN Ia event, in conflict with X-ray line width data. Additionally, in contrast to prior work that suggested a large supersolar metallicity for the WD progenitor for SNR 3C 397, we find satisfactory agreement for solar- and subsolar-metallicity progenitors. We discuss a range of implications our results have for the single-degenerate channel.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Importance of Resolving the Spectral Support of Beam-plasma Instabilities in Simulations

Many astrophysical plasmas are prone to beam-plasma instabilities. For relativistic and dilute beams, the spectral support of the beam-plasma instabilities is narrow, i.e., the linearly unstable modes that grow with rates comparable to the maximum growth rate occupy a narrow range of wavenumbers. This places stringent requirements on the box-sizes when simulating the evolution of the instabilities. We identify the implied lower limits on the box size imposed by the longitudinal beam plasma instability, i.e., typically the most stringent condition required to correctly capture the linear evolution of the instabilities in multidimensional simulations. We find that sizes many orders of magnitude larger than the resonant wavelength are typically required. Using one-dimensional particle-in-cell simulations, we show that the failure to sufficiently resolve the spectral support of the longitudinal instability yields slower growth and lower levels of saturation, potentially leading to erroneous physical conclusion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Sub-Chandrasekhar-mass White Dwarf Detonations Revisited

The detonation of a sub-Chandrasekhar-mass white dwarf (WD) has emerged as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. Recent studies have suggested that the rapid transfer of a very small amount of helium from one WD to another is sufficient to ignite a helium shell detonation that subsequently triggers a carbon core detonation, yielding a "dynamically driven double-degenerate double-detonation" SN Ia. Because the helium shell that surrounds the core explosion is so minimal, this scenario approaches the limiting case of a bare C/O WD detonation. Motivated by discrepancies in previous literature and by a recent need for detailed nucleosynthetic data, we revisit simulations of naked C/O WD detonations in this paper. We disagree to some extent with the nucleosynthetic results of previous work on sub-Chandrasekhar-mass bare C/O WD detonations; for example, we find that a median-brightness SN Ia is produced by the detonation of a 1.0 ${M}_{\odot }$ WD instead of a more massive and rarer 1.1 ${M}_{\odot }$ WD. The neutron-rich nucleosynthesis in our simulations agrees broadly with some observational constraints, although tensions remain with others. There are also discrepancies related to the velocities of the outer ejecta and light curve shapes, but overall our synthetic light curves and spectra are roughly consistent with observations. We are hopeful that future multidimensional simulations will resolve these issues and further bolster the dynamically driven double-degenerate double-detonation scenario's potential to explain most SNe Ia.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

The Ensemble Photometric Variability of Over 10 5 Quasars in the Dark Energy Camera Legacy Survey and the Sloan Digital Sky Survey

We present the ensemble variability analysis results of quasars using the Dark Energy Camera Legacy Survey (DECaLS) and the Sloan Digital Sky Survey (SDSS) quasar catalogs. Our data set includes 119,305 quasars with redshifts up to 4.89. Combining the two data sets provides a 15 year baseline and permits the analysis of the long timescale variability. Adopting a power-law form for the variability structure function, $V=A(t/$1years)$^γ$, we use the multidimensional parametric fitting to explore the relationships between the quasar variability amplitude and a wide variety of quasar properties, including redshift (positive), bolometric luminosity (negative), rest-frame wavelength (negative), and black hole mass (uncertain). We also find that γ can be also expressed as a function of redshift (negative), bolometric luminosity (positive), rest-frame wavelength (positive), and black hole mass (positive). Tests of the fitting significance with the bootstrap method show that, even with such a large quasar sample, some correlations are marginally significant. The typical value of $γ$ for the entire data set is $≳0.25$, consistent with the results in previous studies on both the quasar ensemble variability and the structure function. A significantly negative correlation between the variability amplitude and the Eddington ratio is found, which may be explained as an effect of accretion disk instability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Exploring Fundamentally Three-dimensional Phenomena in High-fidelity Simulations of Core-collapse Supernovae

The details of the physical mechanism that drives core-collapse supernovae (CCSNe) remain uncertain. While there is an emerging consensus on the qualitative outcome of detailed CCSN mechanism simulations in 2D, only recently have high-fidelity 3D simulations become possible. In this work, we present the results of an extensive set of 3D CCSN simulations using high-fidelity multidimensional neutrino transport, high-resolution hydrodynamics, and approximate general relativistic gravity. We employ a state-of-the-art 20 M ⊙ progenitor generated using Modules for Experiments in Stellar Astrophysics, and the SFHo equation of state. While none of our 3D CCSN simulations explode within ~500 ms after core bounce, we find that the presence of large-scale aspherical motion in the Si and O shells aid shock expansion and bring the models closer to the threshold of explosion. We also find some dependence on resolution and geometry (octant versus full 4π). As has been noted in other recent works, we find that the post-shock turbulence plays an important role in determining the overall dynamical evolution of our simulations. We find a strong standing accretion shock instability (SASI) that develops at late times. The SASI produces transient shock expansions, but these do not result in any explosions. We also report that for a subset of our simulations, we find conclusive evidence for the lepton-number emission self-sustained asymmetry, which until now has not been confirmed by independent simulation codes. Finally, both the progenitor asphericities and the SASI-induced transient shock expansion phases generate transient gravitational waves and neutrino signal modulations via perturbations of the protoneutron star by turbulent motions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Relativistic Hydrodynamics with Wavelets

Methods to solve the relativistic hydrodynamic equations are important in a large number of astrophysical simulations, which may be very dynamic and involve multiscale features. This requires computational methods that are highly adaptive and capable of automatically resolving numerous localized features and instabilities that emerge across the computational domain and over many temporal scales. While this has been historically accomplished with adaptive-mesh-refinement-based methods, alternatives using wavelet bases and the wavelet transformation have recently achieved significant success in adaptive representation for advanced engineering applications. The current work presents a new method, extending the wavelet adaptive multiresolution representation method, for the integration of the relativistic hydrodynamic equations using iterated interpolating wavelets and introduces a highly adaptive implementation for multidimensional simulation. We note the wavelet coefficients provide a direct measure of the local approximation error for the solution and place collocation points that naturally adapt to the fluid flow while providing good conservation of fluid quantities. The resulting implementation, OAHU , is applied to a series of demanding 1D and 2D problems that explore high Lorentz factor outflows and the formation of several instabilities, including the Kelvin–Helmholtz instability and the Rayleigh–Taylor instability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Numerical Stability of Detonations in White Dwarf Simulations

Some simulations of SNe Ia feature self-consistent thermonuclear detonations. However, these detonations are not meaningful if the simulations are not resolved, so it is important to establish the requirements for achieving a numerically converged detonation. In this study we examine a test detonation problem inspired by collisions of white dwarfs. This test problem demonstrates that achieving a converged thermonuclear ignition requires a spatial resolution much finer than 1 km in the burning region. Current computational resource constraints place this stringent resolution requirement out of reach for multidimensional supernova simulations. Consequently, contemporary simulations that self-consistently demonstrate detonations are possibly not converged and should be treated with caution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Toward Rate Estimation for Transient Surveys. I. Assessing Transient Detectability and Volume Sensitivity for iPTF

The past couple of decades have seen an emergence of transient detection facilities in various avenues of time-domain astronomy that have provided us with a rich data set of transients. The rates of these transients have implications in star formation, progenitor models, evolution channels, and cosmology measurements. The crucial component of any rate calculation is the detectability and spacetime volume sensitivity of a survey to a particular transient type as a function of many intrinsic and extrinsic parameters. Fully sampling that multidimensional parameter space is challenging. Instead, we present a scheme to assess the detectability of transients using supervised machine learning. The data product is a classifier that determines the detection likelihood of sources resulting from an image subtraction pipeline associated with time-domain survey telescopes, taking into consideration the intrinsic properties of the transients and the observing conditions. We apply our method to assess the spacetime volume sensitivity of type Ia supernovae (SNe Ia) in the intermediate Palomar Transient Factory (iPTF) and obtain the result, &lang; VT &rang; Ia = (2.93 ± 0.21) × 10 -2 Gpc 3 yr. With rate estimates in the literature, this volume sensitivity gives a count of 680-1160 SNe Ia detectable by iPTF, which is consistent with the archival data. With a view toward wider applicability of this technique we do a preliminary computation for long-duration type IIp supernovae (SNe IIp) and find &lang; VT &rang; IIa = (7.80 ± 0.76) × 10 -4 Gpc 3 yr. This classifier can be used for computationally fast spacetime volume sensitivity calculation of any generic transient type using their light-curve properties. Hence, it can be used as a tool to facilitate calculation of transient rates in a range of time-domain surveys, given suitable training sets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Simulating Turbulence-aided Neutrino-driven Core-collapse Supernova Explosions in One Dimension

The core-collapse supernova (CCSN) mechanism is fundamentally 3D, with instabilities, convection, and turbulence playing crucial roles in aiding neutrino-driven explosions. Simulations of CCNSe including accurate treatments of neutrino transport and sufficient resolution to capture key instabilities remain among the most expensive numerical simulations in astrophysics, prohibiting large parameter studies in 2D and 3D. Studies spanning a large swath of the incredibly varied initial conditions of CCSNe are possible in 1D, though such simulations must be artificially driven to explode. We present a new method for including the most important effects of convection and turbulence in 1D simulations of neutrino-driven CCSNe, called Supernova Turbulence In Reduced-dimensionality, or STIR. Our new approach includes crucial terms resulting from the turbulent and convective motions of the flow. We estimate the strength of convection and turbulence using a modified mixinglength theory approach, introducing a few free parameters to the model that are fit to the results of 3D simulations. For sufficiently large values of the mixing-length parameter, turbulence-aided neutrino-driven explosions are obtained. We compare the results of STIR to high-fidelity 3D simulations and perform a parameter study of CCSN explosion using 200 solar-metallicity progenitor models from 9 to 120 M ⊙ . We find that STIR is a better predictor of which models will explode in multidimensional simulations than other methods of driving explosions in 1D. We also present a preliminary investigation of predicted observable characteristics of the CCSN population from STIR, such as the distributions of explosion energies and remnant masses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Dynamics of Laterally Propagating Flames in X-Ray Bursts. I. Burning Front Structure

We investigate the structure of laterally propagating flames through the highly stratified burning layer in an X-ray burst. Two-dimensional hydrodynamics simulations of flame propagation are performed through a rotating plane-parallel atmosphere, exploring the structure of the flame. We discuss the approximations needed to capture the length and timescales at play in an X-ray burst and describe the flame acceleration observed. Finally, our studies complement other multidimensional studies of burning in X-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Carbon–Oxygen Classical Novae Are Galactic 7 Li Producers as well as Potential Supernova Ia Progenitors

Here, we report on studies of classical nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on carbon–oxygen (CO) white dwarfs (WDs). We vary both the mass of the WD (from 0.6 M ⊙ to 1.35 M ⊙ ) and the composition of the accreted material. Our simulations are guided by the results of multidimensional studies of TNRs in WDs, which find that sufficient mixing with WD core material occurs after the TNR is well underway, and levels of enrichment are reached that agree with observations of CN ejecta abundances. We use NOVA (our one-dimensional hydrodynamic code) to accrete solar matter until the TNR is ongoing and then switch to a mixed composition (either 25% WD material and 75% solar or 50% WD material and 50% solar). Because the amount of accreted material is inversely proportional to the initial 12 C abundance, by first accreting solar matter the amount of material taking part in the outburst is larger than in those simulations where we assume a mixed composition from the beginning. Our results show large enrichments of 7 Be in the ejected gases, implying that CO CNe may be responsible for a significant fraction (~100 M ⊙ ) of the 7 Li in the galaxy (~1000 M ⊙ ). Although the ejected gases are enriched in WD material, the WDs in these simulations eject less material than they accrete. We predict that the WD is growing in mass as a consequence of the accretion–outburst–accretion cycle, and CO CNe may be an important channel for SN Ia progenitors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Markov Chain Monte Carlo Predictions of Neutron-rich Lanthanide Properties as a Probe of r -process Dynamics

Lanthanide element signatures are key to understanding many astrophysical observables, from merger kilonova light curves to stellar and solar abundances. To learn about the lanthanide element synthesis that enriched our solar system, we apply the statistical method of Markov Chain Monte Carlo to examine the nuclear masses capable of forming the r-process rare-earth abundance peak. We describe the physical constraints we implement with this statistical approach and demonstrate the use of the parallel chains method to explore the multidimensional parameter space. We apply our procedure to three moderately neutron-rich astrophysical outflows with distinct types of r-process dynamics. We show that the mass solutions found are dependent on outflow conditions and are related to the r-process path. We describe in detail the mechanism behind peak formation in each case. We then compare our mass predictions for neutron-rich neodymium and samarium isotopes to the latest experimental data from the CPT at CARIBU. Here, we find our mass predictions given outflows that undergo an extended (n,γ)⇌(γ,n) equilibrium to be those most compatible with both observational solar abundances and neutron-rich mass measurements.

79 ASTRONOMY AND ASTROPHYSICS↗