Search

Second Moment Methods (SMMs) are developed that are consistent with the Discontinuous Galerkin (DG) spatial discretization of the discrete ordinates (or \Sn) transport equations. The low-order (LO) diffusion system of equations is discretized with fully consistent \Pone, Local Discontinuous Galerkin (LDG), and Interior Penalty (IP) methods. A discrete residual approach is used to derive SMM correction terms that make each of the LO systems consistent with the high-order (HO) discretization. We show that the consistent methods are more accurate and have better solution quality than independently discretized LO systems, that they preserve the diffusion limit, and that the LDG and IP consistent SMMs can be scalably solved in parallel on a challenging, multi-material benchmark problem.

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic S$_N$ setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. Results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used., Comment: To appears in JCP

We present high-order, finite element-based Second Moment Methods (SMMs) for solving radiation transport problems in two spatial dimensions. We leverage the close connection between the Variable Eddington Factor (VEF) method and SMM to convert existing discretizations of the VEF moment system to discretizations of the SMM moment system. The moment discretizations are coupled to a high-order Discontinuous Galerkin discretization of the Discrete Ordinates transport equations. We show that the resulting methods achieve high-order accuracy on high-order (curved) meshes, preserve the thick diffusion limit, and are effective on a challenging multi-material problem both in outer fixed-point iterations and in inner preconditioned iterative solver iterations for the discrete moment systems. We also present parallel scaling results and provide direct comparisons to the VEF algorithms the SMM algorithms were derived from.

We apply high-order mixed finite element discretization techniques and their associated preconditioned iterative solvers to the Variable Eddington Factor (VEF) equations in two spatial dimensions. The mixed finite element VEF discretizations are coupled to a high-order Discontinuous Galerkin (DG) discretization of the Discrete Ordinates transport equation to form effective linear transport algorithms that are compatible with high-order (curved) meshes. This combination of VEF and transport discretizations is motivated by the use of high-order mixed finite element methods in hydrodynamics calculations at the Lawrence Livermore National Laboratory. Due to the mathematical structure of the VEF equations, the standard Raviart Thomas (RT) mixed finite elements cannot be used to approximate the vector variable in the VEF equations. Instead, we investigate three alternatives based on the use of continuous finite elements for each vector component, a non-conforming RT approach where DG-like techniques are used, and a hybridized RT method. We present numerical results that demonstrate high-order accuracy, compatibility with curved meshes, and robust and efficient convergence in iteratively solving the coupled transport-VEF system and in the preconditioned linear solvers used to invert the discretized VEF equations.

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with a high-order Discontinuous Galerkin transport discretization to form an effective high-order, linear transport method. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. In addition, a parallel weak scaling study of the interior penalty VEF discretization demonstrates the scalability of the method out to 1152 processors.

We present a new approach for solving high-order thermal radiative transfer (TRT) using the Variable Eddington Factor (VEF) method (also known as quasidiffusion). Our approach leverages the VEF equations, which consist of the first and second moments of the $S_N$ transport equation, to more efficiently compute the TRT solution for each time step. The scheme consists of two loops - an outer loop to converge the Eddington tensor and an inner loop to converge the iteration between the temperature equation and the VEF system. By converging the outer iteration, one obtains the fully implicit TRT solution for the given time step with a relatively low number of transport sweeps. However, one could choose to perform exactly one outer iteration (and therefore exactly one sweep) per time step, resulting in a semi-implicit scheme that is both highly efficient and robust. Our results indicate that the error between the one-sweep and fully implicit variants of our scheme may be small enough for consideration in many problems of interest., Comment: Submitted to the American Nuclear Society Mathematics & Computation 2021 Conference (The International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, Raleigh, North Carolina, October 3-7, 2021)

For 2x2 block matrices, it is well-known that block-triangular or block-LDU preconditioners with an exact Schur complement (inverse) converge in at most two iterations for fixed-point or minimal-residual methods. Similarly, for saddle-point matrices with a zero (2,2)-block, block-diagonal preconditioners converge in at most three iterations for minimal-residual methods, although they may diverge for fixed-point iterations. But, what happens for non-saddle-point matrices and block-diagonal preconditioners with an exact Schur complement? This note proves that minimal-residual methods applied to general 2x2 block matrices, preconditioned with a block-diagonal preconditioner, including an exact Schur complement, do not (necessarily) converge in a fixed number of iterations. Furthermore, examples are constructed where (i) block-diagonal preconditioning with an exact Schur complement converges no faster than block-diagonal preconditioning using diagonal blocks of the matrix, and (ii) block-diagonal preconditioning with an approximate Schur complement converges as fast as the corresponding block-triangular preconditioning. The paper concludes by discussing some practical applications in neutral-particle transport, introducing one algorithm where block-triangular or block-LDU preconditioning are superior to block-diagonal, and a second algorithm where block-diagonal preconditioning is superior both in speed and simplicity., Comment: 13 pages

We present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin discretizations of the SN transport equation. This approach is sweep-compatible: as the transport sweep is performed, a local quadratic programming (QP) problem is solved in each spatial cell to ensure that the solution satisfies certain physical constraints, including local particle balance. The constraints can be chosen in two ways, leading to two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). The coefficients of the solution are constrained to be nonnegative in QPZ, and they are constrained by an approximate discrete maximum principle in QPMP. There are two primary takeaways in this paper. First, it is shown that the QPMP method, when used with the positive Bernstein basis, eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The second takeaway is that the Variable Eddington Factor (VEF) method can be used to accelerate the convergence of source iteration with fixup for problems with optically thick regions. Regardless of whether a fixup is applied, source iteration converges slowly when optically thick regions are present and acceleration is needed. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle., Comment: Revised draft submitted to Journal of Computational Physics

We present results from general relativistic calculations of the tidal disruption of white dwarf stars from near encounters with intermediate mass black holes. We follow the evolution of 0.2 and $0.6 M_\odot$ stars on parabolic trajectories that approach $10^3$ - $10^4 M_\odot$ black holes as close as a few Schwarzschild radii at periapsis, paying particular attention to the effect tidal disruption has on thermonuclear reactions and the synthesis of intermediate to heavy ion elements. These encounters create diverse thermonuclear environments characteristic of Type I supernovae and capable of producing both intermediate and heavy mass elements in arbitrary ratios, depending on the strength (or proximity) of the interaction. Nuclear ignition is triggered in all of our calculations, even at weak tidal strengths $\beta \sim 2.6$ and large periapsis radius $R_P \sim 28$ Schwarzschild radii. A strong inverse correlation exists between the mass ratio of calcium to iron group elements and tidal strength, with $\beta \lesssim 5$ producing predominately calcium-rich debris. At these moderate to weak interactions, nucleosynthesis is not especially efficient, limiting the total mass and outflows of calcium group elements to $< 15$\% of available nuclear fuel. Iron group elements however continue to be produced in greater quantity and ratio with increasing tidal strength, peaking at $\sim 60$\% mass conversion efficiency in our closest encounter cases. These events generate short bursts of gravitational waves with characteristic frequencies 0.1-0.7 Hz and strain amplitudes $0.5\times10^{-22}$ - $3.5\times10^{-22}$ at 10 Mpc source distance., Comment: 24 pages, 17 figures, accepted for publication in The Astrophysical Journal