Search

Higher order finite difference Weighted Essentially Non-Oscillatory (WENO) schemes have been constructed for conservation laws. For multidimensional problems, they offer high order accuracy at a fraction of the cost of a finite volume WENO or DG scheme of comparable accuracy. This makes them quite attractive for several science and engineering applications. But, to the best of our knowledge, such schemes have not been extended to non-linear hyperbolic systems with non-conservative products. In this paper, we perform such an extension which improves the domain of applicability of such schemes. The extension is carried out by writing the scheme in fluctuation form. We use the HLLI Riemann solver of Dumbser and Balsara (2016) as a building block for carrying out this extension. Because of the use of an HLL building block, the resulting scheme has a proper supersonic limit. The use of anti-diffusive fluxes ensures that stationary discontinuities can be preserved by the scheme, thus expanding its domain of applicability. Our new finite difference WENO formulation uses the same WENO reconstruction that was used in classical versions, making it very easy for users to transition over to the present formulation. For conservation laws, the new finite difference WENO is shown to perform as well as the classical version of finite difference WENO, with two major advantages:- 1) It can capture jumps in stationary linearly degenerate wave families exactly. 2) It only requires the reconstruction to be applied once. Several examples from hyperbolic PDE systems with non-conservative products are shown which indicate that the scheme works and achieves its design order of accuracy for smooth multidimensional flows. Stringent Riemann ... *Abstract truncated, see PDF*, Comment: Accepted in Communications on Applied Mathematics and Computation

GPU computing is expected to play an integral part in all modern Exascale supercomputers. It is also expected that higher order Godunov schemes will make up about a significant fraction of the application mix on such supercomputers. It is, therefore, very important to prepare the community of users of higher order schemes for hyperbolic PDEs for this emerging opportunity. We focus on three broad and high-impact areas where higher order Godunov schemes are used. The first area is computational fluid dynamics (CFD). The second is computational magnetohydrodynamics (MHD) which has an involution constraint that has to be mimetically preserved. The third is computational electrodynamics (CED) which has involution constraints and also extremely stiff source terms. Together, these three diverse uses of higher order Godunov methodology, cover many of the most important applications areas. In all three cases, we show that the optimal use of algorithms, techniques and tricks, along with the use of OpenACC, yields superlative speedups on GPUs! As a bonus, we find a most remarkable and desirable result: some higher order schemes, with their larger operations count per zone, show better speedup than lower order schemes on GPUs. In other words, the GPU is an optimal stratagem for overcoming the higher computational complexities of higher order schemes! Several avenues for future improvement have also been identified. A scalability study is presented for a real-world application using GPUs and comparable numbers of high-end multicore CPUs. It is found that GPUs offer a substantial performance benefit over comparable number of CPUs, especially when all the methods designed in this paper are used., Comment: 49 pages

Recent observations of Type Ia supernovae (SNe) by SH0ES collaboration (R11 and R16) diverge from the value reported by recent CMBR observations utilising the Planck satellite and application of the $\Lambda CDM$ cosmological model by at least $3 \sigma$. It is among the most challenging problems in contemporary cosmology and is known as the Hubble tension. The SNe Ia in R11 and R16 were calibrated through cepheid variables in three distinct galaxies: Milky Way, LMC, and NGC4258. Carnegie Hubble Program (CHP) observations of type Ia SNe calibrated using the tip of the red giant approach yielded a somewhat different estimate for the Hubble constant. This decreased the Hubble tension from over 3$\sigma$ to below 2$\sigma$. It is a legitimate question to answer whether there are any issues with SNe Ia calibration and to investigate whether the Hubble tension is real or not. We use statistical techniques namely, ANOVA, K-S test, and t-test to examine whether the cepheid calibration is host-dependent. Our analysis shows that (i) both R11 and R16 data suffer from non-Gaussian systematic effects, (ii) $H_0$ values in the sub-samples (different anchor-based) in both R11 and R16 groups are significantly different at a 99\% confidence level, and (iii) neglecting the metal-rich MW sample does not reduce the $H_0$ value significantly, and thus Hubble tension persists. A small reduction in the Hubble constant could be linked to the differences in the host environment. Hence instead of using a single universal relation environment based slope and zero point should be preferred., Comment: 15 pages (single column), 7 figures

Spin-orbit coupling in heavy 5$d$ metal oxides, in particular, iridates have received tremendous interest in recent years due to the realization of exotic electronic and magnetic phases. Here, we report the synthesis, structural, magnetic, thermodynamic, and optical properties of the ternary iridate Pr$_3$IrO$_7$. Single crystals of Pr$_3$IrO$_7$ have been grown by the KF flux method. Structural analysis shows that Pr$_3$IrO$_7$ crystallizes in an orthorhombic phase with $Cmcm$ symmetry. The electron energy loss spectroscopy study indicates that Pr is in a 3+ valence state, which implies a 5+ oxidation state of Ir. Magnetization data measured at high and low magnetic fields do not exhibit any bifurcation between $M_{ZFC}$ and $M_{FC}$, however, a weak hump in $M(T)$ is observed at $T^*$$\sim$10.4~K. The specific heat data reveal two maxima at $\sim$253 K and $\sim$4.8 K. The optical conductivity $\sigma_1(\omega)$ spectrum shows 24 infrared-active phonon modes and reveals an insulating behavior with an optical gap $\Delta_{OP}$ of size $\sim$500~meV. During cooling down, the temperature-dependent reflectivity spectrum reveals eight extra phonon modes below the structural phase transition ($\sim$ 253 K). An anomaly is observed at around $T^*$ in the temperature evolution of infrared-active mode frequencies suggesting the presence of significant spin-phonon coupling in the system., Comment: 18 pages, 8 figures, to appear in J. Phys.: Condens. Matter

In this article, we propose high-order finite-difference entropy stable schemes for the two-fluid relativistic plasma flow equations. This is achieved by exploiting the structure of the equations, which consists of three independent flux components. The first two components describe the ion and electron flows, which are modeled using the relativistic hydrodynamics equation. The third component is Maxwell's equations, which are linear systems. The coupling of the ion and electron flows, and electromagnetic fields is via source terms only. Furthermore, we also show that the source terms do not affect the entropy evolution. To design semi-discrete entropy stable schemes, we extend the RHD entropy stable schemes in Bhoriya et al. to three dimensions. This is then coupled with entropy stable discretization of the Maxwell's equations. Finally, we use SSP-RK schemes to discretize in time. We also propose ARK-IMEX schemes to treat the stiff source terms; the resulting nonlinear set of algebraic equations is local (at each discretization point). These equations are solved using the Newton's Method, which results in an efficient method. The proposed schemes are then tested using various test problems to demonstrate their stability, accuracy and efficiency.

Understanding the explosion mechanism of type Ia supernova is among the most challenging issues in astrophysics. Accretion of matter on a carbon-oxygen white dwarf from a companion star is one of the most important keys in this regard. Our aim is to study the effects of WD composition on various parameters during the accretion of helium rich matter at a slow rate. We have used the computer simulation code Modules for Experiments in Stellar Astrophysics (MESA) to understand the variations in the properties such as specific heat (\textit{$C_P$}) and degeneracy parameter (\textit{$\eta$}). The profile of specific heat shows a discontinuity and that of the degeneracy parameter shows a dip near the ignition region. As expected, the size of WD decreases and \textit{g} increases during the accretion. However, a red-giant-like expansion is observed after the rapid ignition towards the end. Our study explains the reason behind the delay in onset of helium ignition due to the difference in carbon abundance in a CO-WD. We find that white dwarfs of the lower abundance of carbon accrete slightly longer before the onset of helium ignition., Comment: 5 pages, 4 figures, 3 tables. arXiv admin note: substantial text overlap with arXiv:2111.13141