Showing 1–6 of 6 results for author: Warsa, J S

Consistent Second Moment Methods with Scalable Linear Solvers for Radiation Transport

Authors: Samuel Olivier, Ben S. Southworth, James S. Warsa, HyeongKae Park

Abstract: 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 correc… ▽ More 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. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Reduced-Memory Methods for Linear Discontinuous Discretization of the Time-Dependent Boltzmann Transport Equation

Authors: Rylan C. Paye, Dmitriy Y. Anistratov, Jim E. Morel, James S. Warsa

Abstract: In this paper, new implicit methods with reduced memory are developed for solving the time-dependent Boltzmann transport equation (BTE). One-group transport problems in 1D slab geometry are considered. The reduced-memory methods are formulated for the BTE discretized with the linear-discontinuous scheme in space and backward-Euler time integration method. Numerical results are presented to demonst… ▽ More In this paper, new implicit methods with reduced memory are developed for solving the time-dependent Boltzmann transport equation (BTE). One-group transport problems in 1D slab geometry are considered. The reduced-memory methods are formulated for the BTE discretized with the linear-discontinuous scheme in space and backward-Euler time integration method. Numerical results are presented to demonstrate performance of the proposed numerical methods. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: 11 pages, 5 figures

Report number: LA-UR-23-25062

Multilevel Second-Moment Methods with Group Decomposition for Multigroup Transport Problems

Authors: Dmitriy Y. Anistratov, Joseph M. Coale, James S. Warsa, Jae H. Chang

Abstract: This paper presents multilevel iterative schemes for solving the multigroup Boltzmann transport equations (BTEs) with parallel calculation of group equations. They are formulated with multigroup and grey low-order equations of the Second-Moment (SM) method. The group high-order BTEs and low-order SM (LOSM) equations are solved in parallel. To further improve convergence and increase computational… ▽ More This paper presents multilevel iterative schemes for solving the multigroup Boltzmann transport equations (BTEs) with parallel calculation of group equations. They are formulated with multigroup and grey low-order equations of the Second-Moment (SM) method. The group high-order BTEs and low-order SM (LOSM) equations are solved in parallel. To further improve convergence and increase computational efficiency of algorithms Anderson acceleration is applied to inner iterations for solving the system of multigroup LOSM equations. Numerical results are presented to demonstrate performance of the multilevel iterative methods. △ Less

Submitted 17 February, 2021; originally announced February 2021.

The Role of Inhomogeneities in Supernova Shock Breakout Emission

Authors: Chris L. Fryer, Christopher J. Fontes, James S. Warsa, Pete W. A. Roming, Shane X. Coffing, Suzannah R. Wood

Abstract: The breakout of a supernova blast wave from its progenitor star provides strong constraints on the star and its immediate surroundings. These surroundings are shaped by mass loss from the star and can include a wide variety of inhomogeneities. Here we present results of multi-dimensional radiation-hydrodynamics calculations of the interactions of the supernova blast wave with inhomogeneities in th… ▽ More The breakout of a supernova blast wave from its progenitor star provides strong constraints on the star and its immediate surroundings. These surroundings are shaped by mass loss from the star and can include a wide variety of inhomogeneities. Here we present results of multi-dimensional radiation-hydrodynamics calculations of the interactions of the supernova blast wave with inhomogeneities in the immediate surroundings of a massive Wolf-Rayet star, calculating the effect these interactions have on the shock breakout signal from supernovae. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 14 pages, 12 figures, submitted to ApJ

Acceleration of the Sn Equations with Highly Anisotropic Scattering using the Fokker-Planck Equation

Authors: Japan K. Patel, James S. Warsa, Anil K. Prinja

Abstract: The discrete ordinates method can model forward-peaked transport problems accurately. However, convergence of discrete ordinates solution can become arbitrarily slow upon use of standard iterative procedures like source iteration and GMRES. Standard zeroth and first moment-based acceleration methods like nonlinear diffusion acceleration and diffusion synthetic acceleration are ineffective in accel… ▽ More The discrete ordinates method can model forward-peaked transport problems accurately. However, convergence of discrete ordinates solution can become arbitrarily slow upon use of standard iterative procedures like source iteration and GMRES. Standard zeroth and first moment-based acceleration methods like nonlinear diffusion acceleration and diffusion synthetic acceleration are ineffective in accelerating such problems because these methods do not correct higher order Legendre-moments of angular flux. We explore the idea of using Fokker-Planck as a preconditioner to accelerate forward-peaked transport problems in this paper. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: This is a draft for our FPSA paper. It will be edited further and submitted to a journal soon

Nonlinear Krylov Acceleration Applied to a Discrete Ordinates Formulation of the k-Eigenvalue Problem

Authors: Matthew T. Calef, Erin D. Fichtl, James S. Warsa, Markus Berndt, Neil N. Carlson

Abstract: We compare variants of Anderson Mixing with the Jacobian-Free Newton-Krylov and Broyden methods applied to an instance of the k-eigenvalue formulation of the linear Boltzmann transport equation. We present evidence that one variant of Anderson Mixing finds solutions in the fewest number of iterations. We examine and strengthen theoretical results of Anderson Mixing applied to linear problems. We compare variants of Anderson Mixing with the Jacobian-Free Newton-Krylov and Broyden methods applied to an instance of the k-eigenvalue formulation of the linear Boltzmann transport equation. We present evidence that one variant of Anderson Mixing finds solutions in the fewest number of iterations. We examine and strengthen theoretical results of Anderson Mixing applied to linear problems. △ Less

Submitted 17 January, 2013; v1 submitted 15 December, 2011; originally announced December 2011.

Comments: This final revision includes results of the C5G7-MOX problem; Nonlinear Krylov Acceleration Applied to a Discrete Ordinates Formulation of the k-Eigenvalue Problem, Accepted by the Journal of Computational Physics December 2012