American Nuclear Society
Home

Home / Store / Journals / Electronic Articles / Nuclear Science and Engineering / Volume 155 / Number 2 / Pages 290-299

A Quasilinear Implicit Riemann Solver for the Time-Dependent Pn Equations

Ryan G. McClarren, James Paul Holloway, Thomas A. Brunner, Thomas A. Mehlhorn

Nuclear Science and Engineering / Volume 155 / Number 2 / Pages 290-299

February 2007

Format:

Price:$30.00
Member Price:$27.00
Member Savings:$3.00

An implicit Riemann solver for the one- and two-dimensional time-dependent spherical harmonics approximation (Pn) to the linear transport equation is presented. This spatial discretization scheme is based on cell-averaged quantities and uses a monotonicity-preserving high resolution method to achieve second-order accuracy (away from extreme points in the solution). Such a spatial scheme requires a nonlinear method of reconstructing the slope within a spatial cell. We have devised a means of creating an implicit (in time) method without the necessity of a nonlinear solver. This is done by computing a time step using a first-order scheme and then, based on that solution, reconstructing the slope in each cell, an implementation that we justify by analyzing the model equation for the method. This quasilinear approach produces smaller errors in less time than both a first-order scheme and a method that solves the full nonlinear system using a Newton-Krylov method.

 
 
 
Questions or comments about the site? Contact the ANS Webmaster.
advertisement