Nonlinear variants of the TR/BDF2 method for thermal radiative diffusion

被引:21
作者
Edwards, Jarrod D. [1 ]
Morel, Jim E. [1 ]
Knoll, Dana A. [2 ]
机构
[1] Texas A&M Univ, Dept Nucl Engn, Zachry Engn Ctr 129, College Stn, TX 77843 USA
[2] Los Alamos Natl Lab, Fluid Dynam & Solid Mech Grp T3, Los Alamos, NM 87545 USA
关键词
Trapezoidal BDF-2; nonlinear radiative diffusion;
D O I
10.1016/j.jcp.2010.10.035
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We apply the Trapezoidal/BDF2 (TR/BDF2) temporal discretization scheme to nonlinear grey radiative diffusion. This is a scheme that is not well-known within the radiation transport community, but we show that it offers many desirable characteristics relative to other second-order schemes. Several nonlinear variants of the TR/BDF2 scheme are defined and computationally compared with the Crank-Nicholson scheme. It is found for our test problems that the most accurate TR/BDF2 schemes are those that are fully iterated to nonlinear convergence, but the most efficient TR/BDF2 scheme is one based upon a single Newton iteration. It is also shown that neglecting the contributions to the Jacobian matrix from the cross-sections, which is often done due to a lack of smooth interpolations for tabular cross-section data, has a significant impact upon efficiency. (C) 2010 Elsevier Inc. All rights reserved.
引用
收藏
页码:1198 / 1214
页数:17
相关论文
共 12 条
[1]  
[Anonymous], 2008, Numerical Methods for Ordinary Differential Equations
[2]   TRANSIENT SIMULATION OF SILICON DEVICES AND CIRCUITS [J].
BANK, RE ;
COUGHRAN, WM ;
FICHTNER, W ;
GROSS, EH ;
ROSE, DJ ;
SMITH, RK .
IEEE TRANSACTIONS ON COMPUTER-AIDED DESIGN OF INTEGRATED CIRCUITS AND SYSTEMS, 1985, 4 (04) :436-451
[3]  
Hairer E., 2002, SPRINGER SERIES COMP
[4]   Numerical analysis of time integration errors for nonequilibrium radiation diffusion [J].
Knoll, D. A. ;
Lowrie, R. B. ;
Morel, J. E. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2007, 226 (02) :1332-1347
[5]   Jacobian-free Newton-Krylov methods: a survey of approaches and applications [J].
Knoll, DA ;
Keyes, DE .
JOURNAL OF COMPUTATIONAL PHYSICS, 2004, 193 (02) :357-397
[6]   ANALYSIS OF A MONTE-CARLO METHOD FOR NONLINEAR RADIATIVE-TRANSFER [J].
LARSEN, EW ;
MERCIER, B .
JOURNAL OF COMPUTATIONAL PHYSICS, 1987, 71 (01) :50-64
[7]  
Lowrie RB, 2004, J COMPUT PHYS, V196, P566, DOI [10.1016/j.jcp.2003.11.016, 10.1016/i.jcp.2003.11.016]
[8]  
Reed WH., 1973, LAUR73479 LOS AL SCI, P20
[9]  
STOER J, 2002, BULIRSCH INTRO NUMER
[10]  
WAREING TA, 1999, P INT C MATH COMP RE, V1, P45