A hybrid Monte Carlo, discontinuous Galerkin method for linear kinetic transport equations
- Univ. of Tennessee, Knoxville, TN (United States)
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
- University of Notre Dame, IN (United States)
Here we present a hybrid method for time-dependent particle transport problems that combines Monte Carlo (MC) estimation with deterministic solutions based on discrete ordinates. For spatial discretizations, the MC algorithm computes a piecewise constant solution and the discrete ordinates use bilinear discontinuous finite elements. From the hybridization of the problem, the resulting problem solved by Monte Carlo is scattering free, resulting in a simple, efficient solution procedure. Between time steps, we use a projection approach to “relabel” collided particles as uncollided particles. In conclusion, from a series of standard 2-D Cartesian test problems we observe that our hybrid method has improved accuracy and reduction in computational complexity of approximately an order of magnitude relative to standard discrete ordinates solutions.
- Research Organization:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
- Sponsoring Organization:
- USDOE
- Grant/Contract Number:
- AC05-00OR22725
- OSTI ID:
- 2538459
- Journal Information:
- Journal of Computational Physics, Journal Name: Journal of Computational Physics Vol. 514; ISSN 0021-9991
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
A Hybrid Finite-Volume, Discontinuous Galerkin Discretization for the Radiative Transport Equation
A Hybrid Monte Carlo-Deterministic Method for Global Binary Stochastic Medium Transport Problems