Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Relativity”

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 37 records · Page 2

Detectability of QCD phase transitions in binary neutron star mergers: Bayesian inference with the next generation gravitational wave detectors

We study the detectability of postmerger QCD phase transitions in neutron star binaries with next-generation gravitational-wave detectors Cosmic Explorer and Einstein Telescope. We perform numerical relativity simulations of neutron star mergers with equations of state that include a quark deconfinement phase transition through either a Gibbs or Maxwell construction. These are followed by Bayesian parameter estimation of the associated gravitational-wave signals using the nrpmw waveform model, with priors inferred from the analysis of the inspiral signal. We assess the ability of the model to measure the postmerger peak frequency $f$$^{peak}_{2}$ and identify aspects that should be improved in the model. We show that, even at postmerger signal to noise ratios as low as 10, the model can distinguish (at the 90% level) $f$$^{peak}_{2}$ between binaries with and without a phase transition in most cases. Phase-transition induced deviations in the $f$$^{peak}_{2}$ from the predictions of equation-of-state insensitive relations can also be detected if they exceed 1.6⁢σ. Our results suggest that next-generation gravitational wave detectors can measure phase transition effects in binary neutron star mergers. Furthermore, unless the phase transition is “strong,” disentangling it from other hadronic physics uncertainties will require significant theory improvements.

79 ASTRONOMY AND ASTROPHYSICS↗

AT2017gfo: Bayesian inference and model selection of multicomponent kilonovae and constraints on the neutron star equation of state

The joint detection of the gravitational wave GW170817, of the short γ-ray burst GRB170817A and of the kilonova AT2017gfo, generated by the the binary neutron star (NS) merger observed on 2017 August 17, is a milestone in multimessenger astronomy and provides new constraints on the NS equation of state. We perform Bayesian inference and model selection on AT2017gfo using semi-analytical, multicomponents models that also account for non-spherical ejecta. Observational data favour anisotropic geometries to spherically symmetric profiles, with a log-Bayes’ factor of ~104, and favour multicomponent models against single-component ones. The best-fitting model is an anisotropic three-component composed of dynamical ejecta plus neutrino and viscous winds. Using the dynamical ejecta parameters inferred from the best-fitting model and numerical–relativity relations connecting the ejecta properties to the binary properties, we constrain the binary mass ratio to q < 1.54 and the reduced tidal parameter to $120\lt \tilde{\Lambda }\lt 1110$. Finally, we combine the predictions from AT2017gfo with those from GW170817, constraining the radius of a NS of 1.4 M ⊙ to 12.2 ± 0.5 km (1σ level). This prediction could be further strengthened by improving kilonova models with numerical-relativity information.

79 ASTRONOMY AND ASTROPHYSICS↗

AthenaK: A Performance-portable Version of the Athena++ Adaptive Mesh Refinement Framework

We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic, and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test particles (e.g., cosmic rays) are implemented using the framework. In two companion papers, we describe (1) a new solver for the Einstein equations based on the Z4c formalism, and (2) a GRMHD solver in dynamical spacetimes also implemented using the framework, enabling new applications in numerical relativity. By adopting Kokkos, the code can be run on virtually any hardware, including CPUs, GPUs from multiple vendors, and emerging Advanced RISC Machine processors. AthenaK shows excellent performance and weak scaling, achieving over 1 billion cell updates per second for hydrodynamics in three dimensions on a single NVIDIA Grace Hopper processor. It does this with a typical parallel efficiency of 80% on 65,536 AMD GPUs on the OLCF Frontier system. Such performance portability enables AthenaK to leverage modern exascale computing systems for challenging applications in astrophysical fluid dynamics, numerical relativity, and multimessenger astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

Neutrino Flavor Transformation in Neutron Star Mergers

Here, we present the first numerical relativity simulations including neutrino flavor transformations that could result from flavor instabilities, quantum many-body effects, or potential beyond standard model physics in neutron star mergers. We find that neutrino flavor transformations impact the composition and structure of the remnant, potentially leaving an imprint on the postmerger gravitational-wave signal. They also have a significant impact on the composition and nucleosynthesis yields of the ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗

Conservative Numerical Schemes with Optimal Dispersive Wave Relations: Part II. Numerical Evaluations

A new energy and enstrophy conserving scheme (EEC) for the shallow water equations is proposed and evaluated using a suite of test cases over the global spherical or bounded domain. The evaluation is organized around a set of pre-defined properties: accuracy of individual operators, accuracy of the whole scheme, conservation of key quantities, control of the divergence variable, representation of the energy and enstrophy spectra, and simulation of nonlinear dynamics. The results confirm that the scheme is between the first and second order accurate, and conserves the total energy and potential enstrophy up to the time truncation errors. Here, the scheme is capable of producing more physically realistic energy and enstrophy spectra, indicating that it can help prevent the unphysical energy cascade towards the finest resolvable scales. With an optimal representation of the dispersive wave relations, the scheme is able to keep the flow close to being non-divergent, and maintain the geostrophically balanced structures with large-scale geophysical flows over long-term simulations.

54 ENVIRONMENTAL SCIENCES↗

Oscillon preheating in full general relativity

Abstract Oscillons are dense nonperturbative objects that may be copiously produced in the preheating period after inflation. Oscillon preheatings are usually simulated with nonlinear matter interactions but in a rigid FLRW background, without taking into account the gravitational backreactions. We investigate the oscillon preheating scenario in full general relativity with a full numerical relativity scheme, and compare the general-relativistic simulations with the traditional ones done in an FLRW background. We find that in certain parameter regions the general-relativistic corrections are significant, producing more and denser oscillons, and can be strong enough to collapse the oscillons to black holes.

Astronomy & Astrophysics↗

Second release of the CoRe database of binary neutron star merger waveforms

Abstract We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity ( CoRe ) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to ℓ = m = 4 . They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Signatures of deconfined quark phases in binary neutron star mergers

We investigate the quark deconfinement phase transition in the context of binary neutron star (BNS) mergers. We employ a new finite-temperature composition-dependent equation of state (EOS) with a first order phase transition between hadrons and deconfined quarks to perform numerical relativity simulations of BNS mergers. Here, the softening of the EOS due to the phase transition causes the merger remnants to be more compact and to collapse to a black hole (BH) at earlier times. The phase transition is imprinted on the postmerger gravitational wave (GW) signal duration, amplitude, and peak frequency. However, this imprint is only detectable for binaries with sufficiently long-lived remnants. Moreover, the phase transition does not result in significant deviations from quasi-universal relations for the postmerger GW peak frequency. We also study the impact of the phase transition on dynamical ejecta, remnant accretion disk masses, r-process nucleosynthetic yields and associated electromagnetic (EM) counterparts. While there are differences in the EM counterparts and nucleosynthesis yields between the purely hadronic models and the models with phase transitions, these can be primarily ascribed to the difference in remnant collapse time between the two. An exception is the non-thermal afterglow caused by the interaction of the fastest component of the dynamical ejecta and the interstellar medium, which is systematically boosted in the binaries with phase transition as a consequence of the more violent merger they experience.

79 ASTRONOMY AND ASTROPHYSICS↗

Reverse phase transitions in binary neutron-star systems with exotic-matter cores

Multimessenger observations of binary neutron-star mergers provide a unique opportunity to constrain the dense-matter equation of state. Although it is known from quantum chromodynamics that hadronic matter will undergo a phase transition to exotic forms of matter—e.g., quark matter—the onset density of such a phase transition cannot be computed from first principles. Hence, it remains an open question if such phase transitions occur inside isolated neutron stars or during binary neutron-star mergers, or if they appear at even higher densities that are not realized in the cosmos. Here, in this article, we perform numerical relativity simulations of neutron-star mergers and investigate scenarios in which the onset density of such a phase transition is exceeded in at least one inspiraling binary component. Our simulations reveal that shortly before the merger, it is possible that such stars undergo a “reverse phase transition”—i.e., densities decrease and the quark core inside the star disappears, leaving a purely hadronic star at merger. After the merger, when densities increase once more, the phase transition occurs again and leads, in the cases considered in this work, to the rapid formation of a black hole. We compute the gravitational-wave signal and the mass ejection for our simulations of such scenarios and find clear signatures that are related to the postmerger phase transition—e.g., smaller ejecta masses due to the softening of the equation of state through the quark core formation. Unfortunately, we do not find measurable imprints of the reverse phase transition.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark magnetohydrodynamics: Black hole accretion in superradiant dark photon clouds

Black holes threaded by massive vector fields can be subject to a superradiant instability, growing a cloud of massive vector particles around it. In this work, we consider what happens if such a dark matter candidate field mimicking a dark photon interacts with an accretion flow onto the black hole. By including a kinetic mixing term with the standard model photon, we extend the commonly used equations of general-relativistic magnetohydrodynamics to a dark photon constituent. The coupling to the dark photon then appears as an effective dynamo term together with a dark Lorentz force acting on the accreting matter. We numerically study the interactions between the superradiant dark photon cloud and the inner accretion flow by solving the coupled system in full numerical relativity. By parameterically varying the mixing parameter between the dark and standard model sector, we provide a first investigation of how the accretion flow could be modified. In conclusion, depending on the coupling strength, our solutions exhibit increased wind launching, as well as oscillation modes in the disk.

79 ASTRONOMY AND ASTROPHYSICS↗

Impact of neutrino flavor conversions on neutron star merger dynamics, ejecta, nucleosynthesis, and multimessenger signals

We present numerical relativity simulations of binary neutron star mergers incorporating neutrino flavor transformations triggered by fast flavor instability, quantum many-body effects, or potential beyond standard model physics. In both long-lived and short-lived remnant scenarios, neutrino flavor conversions modify species-dependent neutrino luminosities and mean energies, and drive the matter towards more neutron rich conditions. They produce up to more neutron rich ejecta and significantly boost the r-process yields, especially in low-density, near-equatorial outflows. We identify regions unstable to fast flavor instabilities and find that these instabilities persist despite flavor conversions. We further test the sensitivity to the equilibration timescale of the flavor conversions, finding that slower flavor conversions can interact with thermodynamic equilibration, and increase the neutron richness of the ejecta. Flavor conversions may also contribute to stronger gravitational wave and neutrino emissions, pointing to a correlation between neutrino transport and merger dynamics. Furthermore, these results highlight the potential impact of flavor conversions while motivating future work to improve on theoretical understanding of flavor instabilities in global simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Do black holes remember what they are made of?

We study the ringdown signal of black holes formed in prompt-collapse binary neutron star mergers. We analyze data from 47 numerical relativity simulations. We show that the ($l$ = 2, $m$ = 2) and ($l$ = 2, $m$ = 1) multipoles of the gravitational wave signal are well fitted by decaying damped exponentials, as predicted by black-hole perturbation theory. We show that the ratio of the amplitude in the two modes depends on the progenitor binary mass ratio q and reduced tidal parameter $\overline{Λ}$. Unfortunately, the numerical uncertainty in our data is too large to fully quantify this dependency. If confirmed, these results will enable novel tests of general relativity in the presence of matter with next-generation gravitational-wave observatories.

79 ASTRONOMY AND ASTROPHYSICS↗

SERS and the scientific career of Richard P. Van Duyne (1945–2019)

This article describes the scientific career of our dear friend Richard P. Van Duyne. Rick was the discoverer of surface-enhanced Raman spectroscopy (SERS), and over the period from the mid-1970s until his death in 2019, he served as both the intellectual leader of the field, and the technical genius who pushed the field forward through his numerous inventions and discoveries. Rick's work went through many phases, and even involved a significant period of inactivity that Rick was able to break out of, but especially after 2000, there was a period of incredible creativity that ultimately led to a detailed understanding of what the SERS enhancement factor is all about and also led to fabrication techniques for making optimum SERS substrates such that applications of SERS to a wide variety of real-world problems in sensing became possible. In addition, Rick's SERS work spawned numerous related technologies, including index of refraction sensing based on localized surface plasmon resonance (LSPR) measurements, and the extension of SERS to its tip-enhanced cousin TERS, to electrochemical versions of SERS and TERS, and to ultrafast nonlinear techniques. The practical applications he pursued included diagnostics for diseases such as cancer, hand-held glucose sensing, and SERS studies related to art conservation. Furthermore, Rick's legacy will probably be best known through the numerous students, postdocs, and collaborators that passed through his lab and who are now active researchers in areas related to SERS at universities and industrial labs around the world.

47 OTHER INSTRUMENTATION↗

From Data to Knowledge: A Graph-Based Reliability Approach to Assess System Health

With the goal of maximizing plant reliability and availability, complex systems such as nuclear power plants continuously monitor and record the performance and the health status of many components, assets, and systems. Such data may take the form of online monitoring data, condition reports, and maintenance reports and it carries the potential to provide system engineers with insights into anomalous behaviors or degradation trends as well as the possible causes behind them and to predict their direct consequences. The analysis of such data poses however few challenges. While some of these challenges are technical in nature (i.e., data are often distributed over several physical servers or databases), others are conceptual in nature (i.e., data elements come in different formats, numeric or textual), and measured values have different scales (e.g., vibration spectra and oil temperature). This paper directly tackles these challenges, and it focuses on the integration of all these data elements in order to assist plant system engineers in analyzing component, assets, and systems performances and optimize maintenance activities. This is performed by 1) extracting knowledge from textual data via technical language processing methods, and 2) quantifying system, asset, and component health from numeric condition-based data. We rely on model-based system engineering (MBSE) models of systems and assets to identify their architecture and functional (i.e., cause and effect) relations. Numeric and textual data elements are then associated with an MBSE graph element, based on their nature. This bonding of MBSE models and data elements constitutes a first-of-its-kind knowledge graph of a nuclear power plants system, with data elements being organized in a structured manner that enables system engineers to identify cause-effect trends in data elements and carry out appropriate actions in response.

97 MATHEMATICS AND COMPUTING↗

Robustness of slow contraction to cosmic initial conditions

We present numerical relativity simulations of cosmological scenarios in which the universe is smoothed and flattened by undergoing a phase of slow contraction and test their sensitivity to a wide range of initial conditions. Our numerical scheme enables the variation of all freely specifiable physical quantities that characterize the initial spatial hypersurface, such as the initial shear and spatial curvature contributions as well as the initial field and velocity distributions of the scalar that drives the cosmological evolution. In particular, we include initial conditions that are far outside the perturbative regime of the well-known attractor scaling solution. We complement our numerical results by analytically performing a complete dynamical systems analysis and show that the two approaches yield consistent results.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1) – Part 2: A semi-discrete error analysis framework for assessing coupling schemes

Abstract. Part 1 (Wan et al., 2024) of this study discusses the motivation and empirical evaluation of a revision to the aerosol-related numerical process coupling in the atmosphere component of the Energy Exascale Earth System Model version 1 (EAMv1) to address the previously reported issue of strong sensitivity of the simulated dust aerosol lifetime and dry removal rate to the model's vertical resolution. This paper complements that empirical justification of the revised scheme with a mathematical justification leveraging a semi-discrete analysis framework for assessing the splitting error of process coupling methods. The framework distinguishes the error due to numerical splitting from the error due to the time integration method(s) used for each individual process. Such a distinction results in a framework that provides an intuitive understanding of the causes of the splitting error. The application of this framework to the dust life cycle in EAMv1 confirms (i) that the original EAMv1 scheme artificially strengthens the effect of dry removal processes and (ii) that the revised splitting reduces that artificial strengthening. While the error analysis framework is presented in the context of the dust life cycle in EAMv1, the framework can be broadly leveraged to evaluate process coupling schemes, both in other physical problems and for any number of processes. This framework will be particularly powerful when the various process implementations support a variety of time integration approaches. Whereas traditional local truncation error approaches require separate consideration of each combination of time integration methods, this framework enables evaluation of coupling schemes independent of particular time integration approaches for each process while still allowing for the incorporation of these specific time integration errors if so desired. The framework also explains how the splitting error terms result from (i) the integration of individual processes in isolation from other processes and (ii) the choices of input state and time step size for the isolated integration of processes. Such a perspective has the potential for the rapid development of alternative coupling approaches that utilize knowledge both about the desired accuracy and about the computational costs of individual processes.

58 GEOSCIENCES↗

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spontaneous Creation of Circularly Polarized Photons in Chiral Astrophysical Systems

This work establishes a relation between chiral anomalies in curved spacetimes and the radiative content of the gravitational field. In particular, we show that a flux of circularly polarized gravitational waves triggers the spontaneous creation of photons with net circular polarization from the quantum vacuum. Using waveform catalogs, we identify precessing binary black holes as astrophysical configurations that emit such gravitational radiation and then solve the fully nonlinear Einstein’s equations with numerical relativity to evaluate the net effect. The quantum amplitude for a merger is comparable to the Hawking emission rate of the final black hole and small to be directly observed. However, the implications for the inspiral of binary neutron stars could be more prominent, as argued on symmetry grounds.

79 ASTRONOMY AND ASTROPHYSICS↗