A spectral analysis of the domain decomposed Monte Carlo method for linear systems

Slattery, Stuart R.; Evans, Thomas M.; Wilson, Paul P. H.

doi:10.1016/j.nucengdes.2015.07.054

Title: A spectral analysis of the domain decomposed Monte Carlo method for linear systems

Full Record
Other Related Research

Abstract

The domain decomposed behavior of the adjoint Neumann-Ulam Monte Carlo method for solving linear systems is analyzed using the spectral properties of the linear oper- ator. Relationships for the average length of the adjoint random walks, a measure of convergence speed and serial performance, are made with respect to the eigenvalues of the linear operator. In addition, relationships for the effective optical thickness of a domain in the decomposition are presented based on the spectral analysis and diffusion theory. Using the effective optical thickness, the Wigner rational approxi- mation and the mean chord approximation are applied to estimate the leakage frac- tion of random walks from a domain in the decomposition as a measure of parallel performance and potential communication costs. The one-speed, two-dimensional neutron diffusion equation is used as a model problem in numerical experiments to test the models for symmetric operators with spectral qualities similar to light water reactor problems. We find, in general, the derived approximations show good agreement with random walk lengths and leakage fractions computed by the numerical experiments.

Authors:

Slattery, Stuart R. ^[1]; Evans, Thomas M. ^[1]; Wilson, Paul P. H. ^[2]

Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Univ. of Wisconsin, Madison, WI (United States)

Publication Date:: Tue Sep 08 00:00:00 EDT 2015

Research Org.:: Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)

Sponsoring Org.:: USDOE Office of Science (SC)

OSTI Identifier:: 1265476

Alternate Identifier(s):: OSTI ID: 1247411

Grant/Contract Number:: AC05-00OR22725

Resource Type:: Accepted Manuscript

Journal Name:: Nuclear Engineering and Design

Additional Journal Information:: Journal Volume: 295; Journal Issue: C; Journal ID: ISSN 0029-5493

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 97 MATHEMATICS AND COMPUTING

Citation Formats


                    Slattery, Stuart R., Evans, Thomas M., and Wilson, Paul P. H. A spectral analysis of the domain decomposed Monte Carlo method for linear systems.  United States: N. p., 2015. 
Web.  doi:10.1016/j.nucengdes.2015.07.054.

Copy to clipboard


                    Slattery, Stuart R., Evans, Thomas M., & Wilson, Paul P. H. A spectral analysis of the domain decomposed Monte Carlo method for linear systems.  United States.  https://doi.org/10.1016/j.nucengdes.2015.07.054

Copy to clipboard


                    Slattery, Stuart R., Evans, Thomas M., and Wilson, Paul P. H. Tue .  
"A spectral analysis of the domain decomposed Monte Carlo method for linear systems".  United States.  https://doi.org/10.1016/j.nucengdes.2015.07.054.  https://www.osti.gov/servlets/purl/1265476.

Copy to clipboard


                    
@article{osti_1265476,

  title        = {A spectral analysis of the domain decomposed Monte Carlo method for linear systems},

  author       = {Slattery, Stuart R. and Evans, Thomas M. and Wilson, Paul P. H.},

  abstractNote = {The domain decomposed behavior of the adjoint Neumann-Ulam Monte Carlo method for solving linear systems is analyzed using the spectral properties of the linear oper- ator. Relationships for the average length of the adjoint random walks, a measure of convergence speed and serial performance, are made with respect to the eigenvalues of the linear operator. In addition, relationships for the effective optical thickness of a domain in the decomposition are presented based on the spectral analysis and diffusion theory. Using the effective optical thickness, the Wigner rational approxi- mation and the mean chord approximation are applied to estimate the leakage frac- tion of random walks from a domain in the decomposition as a measure of parallel performance and potential communication costs. The one-speed, two-dimensional neutron diffusion equation is used as a model problem in numerical experiments to test the models for symmetric operators with spectral qualities similar to light water reactor problems. We find, in general, the derived approximations show good agreement with random walk lengths and leakage fractions computed by the numerical experiments.},

  doi          = {10.1016/j.nucengdes.2015.07.054},

  journal      = {Nuclear Engineering and Design},

  number       = C,

  volume       = 295,

  place        = {United States},

  year         = {Tue Sep 08 00:00:00 EDT 2015},

  month        = {Tue Sep 08 00:00:00 EDT 2015}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1016/j.nucengdes.2015.07.054

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 1 work

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Similar Records in DOE PAGES and OSTI.GOV collections:

A spectral analysis of the domain decomposed Monte Carlo method for linear systems

Conference Slattery, S. R. ; Wilson, P. P. H. ; Evans, T. M.

The domain decomposed behavior of the adjoint Neumann-Ulam Monte Carlo method for solving linear systems is analyzed using the spectral properties of the linear operator. Relationships for the average length of the adjoint random walks, a measure of convergence speed and serial performance, are made with respect to the eigenvalues of the linear operator. In addition, relationships for the effective optical thickness of a domain in the decomposition are presented based on the spectral analysis and diffusion theory. Using the effective optical thickness, the Wigner rational approximation and the mean chord approximation are applied to estimate the leakage fraction ofmore »« less
A Multiple-Set Overlapping-Domain Decomposed Monte Carlo Synthetic Acceleration Method for Linear Systems

Conference Slattery, Stuart R. ; Evans, Thomas M. ; Wilson, Paul P.H. ; ... - SNA + MC 2013 - JOINT INTERNATIONAL CONFERENCE ON SUPERCOMPUTING IN NUCLEAR APPLICATIONS + MONTE CARLO

We present a multiple-set overlapping-domain decomposed strategy for parallelizing the Monte Carlo Synthetic Acceleration method. Monte Carlo Synthetic Acceleration methods use the Neumann-Ulam class of Monte Carlo solvers for linear systems to accelerate a fixed-point iteration sequence. Effective parallel algorithms for these methods require the parallelization of the underlying Neumann-Ulam solvers. To do this in a domain decomposed environment, we borrow strategies traditionally implemented in Monte Carlo particle transport to parallelize the problem. The parallel Neumann-Ulam and multiple-set overlapping-domain decomposition algorithms are presented along with parallel scaling data for the resulting implementation using the Titan Cray XK7 machine at Oakmore » Ridge National Laboratory.« less
https://doi.org/10.1051/snamc/201404211
Transport in Stochastic Media with Random Chord Length Distributions

Technical Report Prinja, Anil K. ; Skinner, Corey

Thermal radiation transport computations in binary Markovian random mixtures rely almost exclusively on the Levermore-Pomraning (LP) model which is obtained by applying a heuristic closure to the ensemble averaged random medium transport equation. The validity of this model has been extensively tested by comparing numerical results over a broad parameter range (material types and mixing parameters) against benchmark solutions in planar geometry. The conditions under which the LP-model provides useful results and when it breaks down are now well established, but work to date has been largely restricted to homogeneous mixing statistics, i.e., the mean chord lengths of both materialsmore »« less
https://doi.org/10.2172/1820548

Full Text Available
Fast multiscale contrast independent preconditioners for linear elastic topology optimization problems

Journal Article Zambrano, Miguel ; Serrano, Sintya ; Lazarov, Boyan S. ; ... - Journal of Computational and Applied Mathematics

The goal of this work is to present a fast and viable approach for the numerical solution of the high-contrast state problems arising in topology optimization. The optimization process is iterative, and the gradients are obtained by an adjoint analysis, which requires the numerical solution of large high-contrast linear elastic problems with features spanning several length scales. The size of the discretized problems forces the utilization of iterative linear solvers with solution time dependent on the quality of the preconditioner. The lack of clear separation between the scales, as well as the high-contrast, imposes severe challenges on the standard preconditioningmore »« less
https://doi.org/10.1016/j.cam.2020.113366

Full Text Available
RADIANT-ENERGY TRANSFER IN GASEOUS FLOWS. Paper 2 of GENERAL RESEARCH IN FLIGHT SCIENCES, JANUARY 1959-JANUARY 1960. VOLUME 1. FLUID MECHANICS. PART 1

Technical Report Tellep, D M ; Edwards, D K

The basic integral expressions for radiant energy transfer near a diffuse reflecting surface with uniform temperature are reviewed. The radiation diffusion approxi-mation is also reviewed. Conservation of energy transferred by radiation, diffusion, conduction, and convection is formulated as an integro- partial differential equation for a gray chemically reacting gas in laminar boundary-layer flow over a black wall. A similarity transform is introduced to reduce to flow equations a pair of simultaneous, ordinary integro-differential equations, and an iterative method of solution is indicated. It is also shown that when the radiation mean free path is small compared with the thermal boundarymore »« less

Similar Records