A Monte Carlo synthetic acceleration method for the non-linear, time-dependent diffusion equation

T. M. Evans, S. W. Mosher

Research output: Chapter in Book/Report/Conference proceedingConference contributionpeer-review

3 Scopus citations

Abstract

We present a Monte Carlo synthetic-acceleration method for solving the time-dependent, nonlinear, equilibrium diffusion equation in three-dimensions. The new scheme uses the adjoint Monte Carlo method as a relaxation step that accelerates standard Jacobi iteration. Results show that this method is 40% faster than regular Jacobi iteration and is competitive with Jacobi-preconditioned Conjugate Gradient methods. Furthermore, the new method is not limited to symmetric, positive-definite systems, and therefore, it can be used for non-symmetric systems. Such systems arise in fully coupled, nonlinear-consistent (Newton) solvers.

Original languageEnglish
Title of host publicationAmerican Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009
Pages1361-1370
Number of pages10
StatePublished - 2009
EventInternational Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009 - Saratoga Springs, NY, United States
Duration: May 3 2009May 7 2009

Publication series

NameAmerican Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009
Volume2

Conference

ConferenceInternational Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009
Country/TerritoryUnited States
CitySaratoga Springs, NY
Period05/3/0905/7/09

Keywords

  • Linear solvers
  • Monte Carlo
  • Radiation diffusion
  • Synthetic acceleration

Fingerprint

Dive into the research topics of 'A Monte Carlo synthetic acceleration method for the non-linear, time-dependent diffusion equation'. Together they form a unique fingerprint.

Cite this