Search NASA⌕ Search

SEARCH · Search NASA

Results for “general relativity equations & solutions”

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 109 records · Page 6

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗

Superluminal matter waves

The Dirac equation has resided among the greatest successes of modern physics since its emergence as the first quantum mechanical theory fully compatible with special relativity. This compatibility ensures that the expectation value of the velocity is less than the vacuum speed of light. Here, we show that the Dirac equation admits free-particle solutions where the peak amplitude of the wave function can travel at any velocity, including those exceeding the vacuum speed of light, despite having a subluminal velocity expectation value. The solutions are constructed by superposing basis functions with correlations in momentum space. These arbitrary velocity wave functions feature a near-constant profile and may impact quantum mechanical processes that are sensitive to the local value of the probability density as opposed to expectation values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Separability of the Dirac equation in the STU black hole spacetime: Pairwise-equal charge case study

We study the separability of the Dirac equation in the background of four dimensional charged rotating asymptotically flat black hole solution of supergravity and the low energy string theory known as STU black hole. In particular, we analyze in detail the separability conditions in the pairwise equal charge STU black hole space-time. While in the latter case the minimally coupled Dirac equation is not separable, the introduction of a specific torsion term ensures the separability. The source of the torsion is the Kalb-Ramond field, which is an integral part of string theory, but further aspects of its properties and coupling to fermionic fields remain to be studied. To derive the torsion, two different approaches are used in conformally related frames, showing that the torsion is not unique. The correspondingly modified Dirac equations in the Einstein and string frames are shown to be separable. Furthermore, the massless radial and angular wave equations are examined; they show close similarity with corresponding equations for the standard Kerr background. A generalization of the Teukolsky equation for the pairwise equal case is conjectured. We also briefly analyze a technically sophisticated radial equation in the massive case.

Dirac equation↗

Perturbations of Gibbons-Maeda black holes in Einstein-Maxwell-dilaton theories

The study of perturbations around black hole backgrounds in general relativity and Einstein-Maxwell theory has a long history, going back to the work of Regge and Wheeler in the 1950s. As part of a broader investigation of perturbations around black holes in supergravity, we describe here our results for the perturbations around the Gibbons-Maeda static charged black holes in a class of Einstein-Maxwell-Dilaton theories. Our analysis follows the general strategy developed by Chandrasekhar and Xanthopoulos for the perturbations of the Reissner-Nordström black hole. Here, the analysis is considerably more involved, because of the presence of the dilaton field, which couples to the other polar modes. We nonetheless find that the problem is completely solvable, in the sense that one can separate variables and eventually describe all the perturbations in terms of diagonalised second-order radial equations. We are able to prove the mode stability of all the Gibbons-Maeda black hole solutions. Published by the American Physical Society 2024

Pope, C. N. (ORCID:0000000284012564)↗

Twistor space origins of the Newman-Penrose map

Recently, we introduced the “Newman-Penrose map”, a novel correspondence between a certain class of solutions of Einstein’s equations and self-dual solutions of the vacuum Maxwell equations, which we showed was closely related to the classical double copy. Here, we give an alternative definition of this correspondence in terms of quantities that are defined naturally on twistor space, and a shear-free null geodesic congruence on Minkowski space whose twistorial character is articulated by the Kerr theorem. The advantage of this reformulation is that it is purely geometrical in nature, being manifestly invariant under both spacetime diffeomorphisms and projective transformations on twistor space. While the original formulation of the map may be more convenient for most explicit calculations, the twistorial formulation we present here may be of greater theoretical utility.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Review—Drag Coefficients of Non-Spherical and Irregularly Shaped Particles

Abstract The knowledge of simple and relatively accurate closure equations for the drag coefficients of nonspherical particles is very important for Eulerian multiphase numerical codes that require such inputs for the solution of the momentum equation. However, the existing correlations for the drag coefficients are laden with high uncertainty, which propagates in the numerical results. This review paper examines critically the state of our knowledge of the drag coefficients of nonspherical particles starting with our understanding of the “size” and “shape” of particles. Lengthscales and other parameters that have been used to characterize sizes of nonspherical particles are presented. Shape factors and shape descriptors that have been used in past correlations are also presented together with an appraisal of their applicability. Several correlations that have been developed for the determination of the drag coefficients of nonspherical particles with regular and irregular shapes are also presented together with available information on their accuracy and applicability. Based on the review of the available correlations general recommendations are offered on the use of the correlations as closure equations and inputs in computatioal fluid dynamics (CFD) codes.

Engineering↗

Extended molecular eigenmodes treatment of dipole–dipole NMR relaxation in real fluids

Traditional models of NMR relaxation fail to account for the complex, multi-exponential behavior of the autocorrelation function in realistic systems characterized by soft-interactions and molecules that are chemically and physically complex. Here, in this study, we describe the relative diffusion of the spin dipoles by means of a Fokker–Planck equation that includes an interaction potential of mean force to account for the response of the physical/chemical environment around the dipoles. By numerically solving the Fokker–Planck equation for the diffusion propagator, we estimate dipole–dipole NMR relaxation for like- and unlike-spin systems via its eigenmode solution. We test the model against molecular simulations of diffusing dipoles with harmonic potentials and also validate using experimental longitudinal relaxation data from real systems, including Gd(III)–aqua and Gd(III)–DO3A–butrol complexes, the latter being an important MRI contrast agent. Using this novel approach, we predict both the inner- and outer-shell contributions to the relaxivity rates with excellent accuracy at frequencies relevant to MRI. We also show that, under the appropriate assumptions, our framework naturally recovers the Bloembergen–Purcell–Pound, the Solomon–Bloembergen–Morgan, and the Hwang–Freed models. Our implementation is general and publicly available for application to a broad range of systems.

Pinheiro dos Santos, Thiago J. [Rice Univ., Housto↗

Redshifting the Cosmological Constant in Unimodular Gravity via Nonlinear Quantum Mechanics

The cosmological constant problem represents a profound conflict between quantum field theory and general relativity. Unimodular gravity offers a compelling starting point by de-gravitating the vacuum energy of the Standard Model, but this framework traditionally trades the problem of vacuum energy for a fine-tuning of initial conditions, which manifest as a ``shadow" cosmological constant. In this paper, we resolve this initial conditions problem by proposing a novel modification to gravity based on nonlinear quantum mechanics. We introduce specific state-dependent terms to the Hamiltonian, constructed from expectation values of the metric such as the average Ricci scalar. These terms alter the dynamical equations of gravity such that the shadow energy density associated with unconstrained initial conditions redshifts away with cosmic expansion, rendering it negligible at late times. The resulting cosmology is naturally dominated by matter and radiation without fine-tuning. We demonstrate that this significant infrared modification of gravity is consistent with local and cosmological tests of gravity. We comment on the possibility of testing this solution in cosmological measurements of Newton's constant.

Kaplan, David E. [Johns Hopkins U.; Tokyo U., IPMU↗

Uniform-density Bose-Einstein condensates of the Gross-Pitaevskii equation found by solving the inverse problem for the confining potential

Here, in this work, we consider a “reverse-engineering” approach to construct confining potentials that support exact, constant density kovaton solutions to the classical Gross-Pitaevskii equation (GPE) also known as the nonlinear Schr¨odinger equation (NLSE). In the one-dimensional case, the exact solution is the sum of stationary kink and anti-kink solutions, i.e. a kovaton, and in the overlapping region, the density is constant. In higher dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed out finite square well with minima also at the edges. When self-interactions are added, a term proportional to ±gψ*ψ gets added to the confining potential and ±gM, where M is the norm, gets added to the total energy. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of M (related to the number of particles in the BEC). Comparing the stability criteria from Derrick’s theorem and Bogoliubov-de Gennes analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass M for attractive interactions. However, the numerical analysis gives a much lower critical mass. The numerical analysis shows further that the initial instabilities violate the symmetry x → -x assumed by Derrick’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A High-Order Eulerian–Lagrangian Runge–Kutta Finite Volume (EL–RK–FV) Method for Scalar Nonlinear Conservation Laws

Abstract We present a class of high-order Eulerian–Lagrangian Runge–Kutta finite volume methods that can numerically solve Burgers’ equation with shock formations, which could be extended to general scalar conservation laws. Eulerian–Lagrangian (EL) and semi-Lagrangian (SL) methods have recently seen increased development and have become a staple for allowing large time-stepping sizes. Yet, maintaining relatively large time-stepping sizes post shock formation remains quite challenging. Our proposed scheme integrates the partial differential equation on a space-time region partitioned by linear approximations to the characteristics determined by the Rankine–Hugoniot jump condition. We trace the characteristics forward in time and present a merging procedure for the mesh cells to handle intersecting characteristics due to shocks. Following this partitioning, we write the equation in a time-differential form and evolve with Runge–Kutta methods in a method-of-lines fashion. High-resolution methods such as ENO and WENO-AO schemes are used for spatial reconstruction. Extension to higher dimensions is done via dimensional splitting. Numerical experiments demonstrate our scheme’s high-order accuracy and ability to sharply capture post-shock solutions with large time-stepping sizes.

Chen, Jiajie↗

Nonlocal Kernel Network (NKN): a Stable and Resolution-Independent Deep Neural Network.

Neural operators have recently become popular tools for designing solution maps between function spaces in the form of neural networks. Differently from classical scientific machine learning approaches that learn parameters of a known partial differential equation (PDE) for a single instance of the input parameters at a fixed resolution, neural operators approximate the solution map of a family of PDEs [6, 7]. Despite their success, the uses of neural operators are so far restricted to relatively shallow neural networks and confined to learning hidden governing laws. In this work, we propose a novel nonlocal neural operator, which we refer to as nonlocal kernel network (NKN), that is resolution independent, characterized by deep neural networks, and capable of handling a variety of tasks such as learning governing equations and classifying images. Our NKN stems from the interpretation of the neural network as a discrete nonlocal diffusion reaction equation that, in the limit of infinite layers, is equivalent to a parabolic nonlocal equation, whose stability is analyzed via nonlocal vector calculus. The resemblance with integral forms of neural operators allows NKNs to capture long-range dependencies in the feature space, while the continuous treatment of node-to-node interactions makes NKNs resolution independent. The resemblance with neural ODEs, reinterpreted in a nonlocal sense, and the stable network dynamics between layers allow for generalization of NKN’s optimal parameters from shallow to deep networks. This fact enables the use of shallow-to-deep initialization techniques [8]. Our tests show that NKNs outperform baseline methods in both learning governing equations and image classification tasks and generalize well to different resolutions and depths.

97 MATHEMATICS AND COMPUTING↗

Revisiting matrix-based inversion of scanning mobility particle sizer (SMPS) and humidified tandem differential mobility analyzer (HTDMA) data

Abstract. Tikhonov regularization is a tool for reducing noise amplification during data inversion. This work introduces RegularizationTools.jl, a general-purpose software package for applying Tikhonov regularization to data. The package implements well-established numerical algorithms and is suitable for systems of up to ∼ 1000 equations. Included is an abstraction to systematically categorize specific inversion configurations and their associated hyperparameters. A generic interface translates arbitrary linear forward models defined by a computer function into the corresponding design matrix. This obviates the need to explicitly write out and discretize the Fredholm integral equation, thus facilitating fast prototyping of new regularization schemes associated with measurement techniques. Example applications include the inversion involving data from scanning mobility particle sizers (SMPSs) and humidified tandem differential mobility analyzers (HTDMAs). Inversion of SMPS size distributions reported in this work builds upon the freely available software DifferentialMobilityAnalyzers.jl. The speed of inversion is improved by a factor of ∼ 200, now requiring between 2 and 5 ms per SMPS scan when using 120 size bins. Previously reported occasional failure to converge to a valid solution is reduced by switching from the L-curve method to generalized cross-validation as the metric to search for the optimal regularization parameter. Higher-order inversions resulting in smooth, denoised reconstructions of size distributions are now included in DifferentialMobilityAnalyzers.jl. This work also demonstrates that an SMPS-style matrix-based inversion can be applied to find the growth factor frequency distribution from raw HTDMA data while also accounting for multiply charged particles. The outcome of the aerosol-related inversion methods is showcased by inverting multi-week SMPS and HTDMA datasets from ground-based observations, including SMPS data obtained at Bodega Marine Laboratory during the CalWater 2/ACAPEX campaign and co-located SMPS and HTDMA data collected at the US Department of Energy observatory located at the Southern Great Plains site in Oklahoma, USA. Results show that the proposed approaches are suitable for unsupervised, nonparametric inversion of large-scale datasets as well as inversion in real time during data acquisition on low-cost reduced-instruction-set architectures used in single-board computers. The included software implementation of Tikhonov regularization is freely available, general, and domain-independent and thus can be applied to many other inverse problems arising in atmospheric measurement techniques and beyond.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the role of the lattice in metal–insulator phase transitions

Many materials exhibit phase transitions at which both the electronic properties and the crystal structure change. Some authors have argued that the change in electronic order is primary, with the lattice distortion a relatively minor side-effect, and others have argued that the lattice distortions play an essential role in the energetics of the transition. In this paper, we introduce a formalism that resolves this long-standing problem. The methodology works with any electronic structure method that produces solutions of the equation of state determining the electronic order parameter as a function of lattice distortion. We use the formalism to settle the question of the physics of the metal–insulator transitions in the rare-earth perovskite nickelates (RNiO 3 ) and Ruddlesden–Popper calcium ruthenates (Ca 2 RuO 4 ) in bulk, heterostructure, and epitaxially strained thin film forms, finding that electron-lattice coupling is key to stabilizing the insulating state in both classes of materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Current-conserving relativistic linear response for collisional plasmas

Here we investigate the response of a relativistic plasma to electromagnetic fields in the framework of the Boltzmann equation incorporating a collision term in the relaxation rate approximation selected in a form assuring current conservation. We obtain an explicit solution for the linearized perturbation of the Fermi–Dirac equilibrium distribution in terms of the average relaxation rate k. We study the resulting covariant, gauge invariant, and current conserving form of the polarization tensor in the ultrarelativistic and non-relativistic limits. We evaluate the susceptibility in the ultrarelativistic limit and explore their dependence on k. Finally, we study the dispersion relations for the longitudinal and transverse poles of the propagator. We show that for k >2ω p , where ω p is the plasma frequency, the plasma wave modes are overdamped. In the opposite case, k$\ll$ ω p , the propagating plasma modes are weakly damped.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

New heavenly double copies

The double copy relates scattering amplitudes and classical solutions in Yang-Mills theory, gravity, and related field theories. Previous work has shown that this has an explicit realisation in self-dual YM theory, where the equation of motion can be written in a form that maps directly to Plebański’s heavenly equation for self-dual gravity. The self-dual YM equation involves an area-preserving diffeomorphism algebra, two copies of which appear in the heavenly equation. In this paper, we show that this construction is a special case of a wider family of heavenly-type examples, by (i) performing Moyal deformations, and (ii) replacing the area-preserving diffeomorphisms with a less restricted algebra. As a result, we obtain a double-copy interpretation for hyper-Hermitian manifolds, extending the previously known hyper-Kähler case. We also introduce a double-Moyal deformation of the heavenly equation. The examples where the construction of Lax pairs is possible are manifestly consistent with Ward’s conjecture, and suggest that the classical integrability of the gravity-type theory may be guaranteed in general by the integrability of at least one of two gauge-theory-type single copies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Technology Strategy Assessment: Findings from Storage Innovations 2030 Bidirectional Hydrogen Storage

Hydrogen is the most common element in the universe, comprising nearly 75% of all normal matter, and it has been used by scientists for centuries, but it was not fully recognized as an element until 1766, when it was isolated by Henry Cavendish. Early work focused on the generation of hydrogen through the oxidation of metals in water, which released hydrogen gas. Hydrogen’s lighter-than-air and flammable properties were immediately used in engines, zeppelins, and as feedstock for a wide variety of chemical reactions. Several approaches were developed for the production of hydrogen with the most common being associated with the production and conversion of hydrocarbon-based fuels. Coal gasification, steam methane reforming, and other reformation processes provide the majority of current hydrogen production due to the relatively low cost of hydrogen produced through these processes. More than 95% of hydrogen production is used for industrial processes rather than energy storage. To facilitate affordable decarbonization of these industrial processes and to advance the use of hydrogen as a fuel in transportation, DOE launched the Hydrogen Shot as part of the Energy Earthshots Initiative. The goal of the Hydrogen Shot is to reduce the cost of clean hydrogen by 80% to $1/kg of clean hydrogen production within one decade (known as the “1 1 1” goal). This is distinct from the Long-Duration Storage Shot, which is the primary focus of this report; however, it is intrinsically linked to bidirectional hydrogen storage. Several important chemical synthesis processes are dependent upon hydrogen, and the production and use of hydrogen is generally driven by its connection to one of these markets. For example, ammonia is one of the most highly produced chemicals in the world and it depends chiefly on hydrogen. Ammonia is primarily used for agricultural fertilizer and is considered to be largely responsible for a doubling of agricultural production per unit of land over the last century. Another one of hydrogen’s primary uses is as a catalyst in petroleum refining during the desulfurization process. Beyond chemical production, hydrogen is used as a reductant in the production of steel and has been demonstrated as a substitute for metallurgical coal in the production of raw iron. It is even used in the hydrogenation reaction for food products to create more shelf-stable semi-solid fats. However, while hydrogen is produced on the order of 100 million metric tons/year globally to feed these industries, more than 95% of hydrogen is produced from hydrocarbons that emit CO2 during the process. Conversely, electrolysis is a process by which electricity is used to separate hydrogen and oxygen in water molecules, usually across a membrane. Hydrogen production via electrolysis lowers the carbon intensity of produced hydrogen when coupled with low-carbon electricity. Currently, global electrolysis capacity is on the order of 1 GW, which equates to about 500 metric tons/day of hydrogen production. To support large-scale industrial decarbonization, capacity will likely need to increase by two to three orders of magnitude. Electrolysis technology is broadly separated into groups that are defined by the electrolyte used, with further subdivision based on the operating characteristics. The majority of commercial electrolyzer systems are based around three main technology groups: liquid alkaline, proton exchange membrane, and solid oxide. Liquid Alkaline (LA) electrolysis is the oldest, most mature, least expensive, and most common commercial technology, with 400 plants in operation by 1902. Its hydrogen output is low relative to the size of the system due to a low current density. LA electrolysis utilizes a liquid potassium hydroxide solution as the electrolyte. Proton exchange membrane (PEM) electrolysis (also known as polymer electrolyte membrane electrolysis), described in 1960, relies on an acid-impregnated polymer membrane as the electrolyte and typically offers three to six times higher hydrogen production per unit cell area than LA electrolysis. Solid oxide electrolysis, or high-temperature electrolysis, utilizes a ceramic cell as the electrolyte and operates on steam rather than liquid water, enabling electrical efficiencies of more than 90%, which is up from 60% with PEM. Two pre-commercial electrolyzer technologies to note are alkaline exchange membrane (AEM) and proton-conducting solid oxide electrolysis cell (SOEC). AEM potentially has the advantages of both LA and PEM technologies in that it is able to use low-cost materials like LA but with the ability to operate at higher output pressures with a smaller footprint like PEM. Proton-conducting SOEC is similar to commercial SOEC, which uses an oxide-conducting ceramic; however, it uses a proton-conducting ceramic that has the potential to operate at lower temperatures and has lower capital costs. Each of these technologies is experiencing a rapid improvement in performance and a reduction in installed cost, and each appears to be well suited to specific applications. Besides differences in the type of electrolyzer used, the main difference in the architecture of bidirectional hydrogen systems is how the hydrogen is stored. Currently, the most cost-effective way to store large amounts of hydrogen gas is underground, such as in large salt caverns that have been hollowed out. These salt caverns are geographically concentrated in small portions of the United States and are not generally near large metropolitan areas; however, other subsurface architectures are being investigated to expand this reach. A more widely deployable option is aboveground pressurized tanks. These systems are about 10 times as expensive because of the materials and safety margins required to hold hydrogen at high pressures. A third option is using materials-based storage, such as liquid organic hydrogen carriers. By reversibly attaching the produced hydrogen to other molecules, it can be stored at near atmospheric pressure and room temperature. This has the potential to reduce the material cost of storage but may result in a reduction in the efficiency of the process because there are both hydrogen uptake and release processes. While materials-based storage has not been used extensively for large-scale hydrogen storage in the past, there is currently significant activity regarding developing materials and processes for use in large-scale hydrogen storage applications. Electrolysis-produced hydrogen offers an unusual opportunity for energy storage applications. Unlike more conventional energy storage approaches, such as batteries, which operate entirely within electrical markets, hydrogen is a valuable product beyond the electric market and can be directed to the most lucrative use. Hydrogen also can be directly converted back to electricity using either a fuel cell or turbine, or it can be sold to other markets, such as chemical synthesis, steel production, or even export. In this way, excess electricity can be upgraded to the most valuable product. Finally, its use can be actively managed between multiple off-takers; for example, local hydrogen storage can provide a specific amount of stored electricity and any excess can be exported to ammonia production. This flexibility is amplified by the fact that hydrogen storage has fully decoupled power and energy components, which allows for affordable scaling options. Together, this allows a substantial amount of creativity to enable the economic utilization of variable power resources while supporting decarbonization of the industry.

08 HYDROGEN↗

The r -modes of slowly rotating, stratified neutron stars

ABSTRACT The only r-modes that exist in a globally barotropic, rotating, Newtonian star are the fundamental l = |m| solutions, where l and m are the indices of the spherical harmonic $Y_l^m$ that describe the mode’s angular dependence. This is in stark contrast to a stellar model that is non-barotropic throughout its interior, which hosts all the l ≥ |m| perturbations including radial overtones. In reality, neutron stars are stratified with locally barotropic regions. Therefore, we explore how stratification alters a star’s ability to support r-modes. We consider the globally stratified case and examine the behaviour of the modes as the star gets close to barotropicity. In this limit, we find that all but the fundamental l = |m| perturbations change character and become generic inertial modes. Restricting the analysis to l = |m| perturbations, we develop the r-mode equations in order to consider stellar models that exhibit local barotropicity. Our results for such models show that the r-mode overtones diverge and join the inertial modes. In order to see which r-modes persist and retain their character in realistic neutron stars, these calculations will need to be brought into full general relativity.

79 ASTRONOMY AND ASTROPHYSICS↗

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.

79 ASTRONOMY AND ASTROPHYSICS↗