BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ryan G. McClarren

Ryan G. McClarren appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.comp-phComputational Engineering, Finance, and Science

Catalog footprint

What is connected

11works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2024arXiv

Mitigating Spatial Error in the iterative-Quasi-Monte Carlo (iQMC) Method for Neutron Transport Simulations with Linear Discontinuous Source Tilting and Effective Scattering and Fission Rate Tallies

The iterative Quasi-Monte Carlo (iQMC) method is a recently proposed method for multigroup neutron transport simulations. iQMC can be viewed as a hybrid between deterministic iterative techniques, Monte Carlo simulation, and Quasi-Monte Carlo techniques. iQMC holds several algorithmic characteristics that make it desirable for high performance computing environments including a $O(N^{-1})$ convergence scheme, ray tracing transport sweep, and highly parallelizable nature similar to analog Monte Carlo. While there are many potential advantages of using iQMC there are also inherent disadvantages, namely the spatial discretization error introduced from the use of a mesh across the domain. This work introduces two significant modifications to iQMC to help reduce the spatial discretization error. The first is an effective source transport sweep, whereby the source strength is updated on-the-fly via an additional tally. This version of the transport sweep is essentially agnostic to the mesh, material, and geometry. The second is the addition of a history-based linear discontinuous source tilting method. Traditionally, iQMC utilizes a piecewise-constant source in each cell of the mesh. However, through the proposed source tilting technique iQMC can utilize a piecewise-linear source in each cell and reduce spatial error without refining the mesh. Numerical results are presented from the 2D C5G7 and Takeda-1 k-eigenvalue benchmark problems. Results show that the history-based source tilting significantly reduces error in global tallies and the eigenvalue solution in both benchmarks. Through the effective source transport sweep and linear source tilting iQMC was able to converge the eigenvalue from the 2D C5G7 problem to less than $0.04\%$ error on a uniform Cartesian mesh with only $204\times204$ cells.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Accurate solutions to time dependent transport problems with a moving mesh and exact uncollided source treatment

For the purpose of finding benchmark quality solutions to time dependent Sn transport problems, we develop a numerical method in a Discontinuous Galerkin (DG) framework that utilizes time dependent cell edges, which we call a moving mesh, and an uncollided source treatment. The DG method for discretizing space is a powerful solution technique on smooth problems and is robust on non-smooth problems. In order to realize the potential of the DG method to spectrally resolve smooth problems, our moving mesh and uncollided source treatment is devised to circumvent discontinuities in the solution or the first derivative of the solutions that are admitted in transport calculations. The resulting method achieves spectral convergence on smooth problems, like a standard DG implementation. When applied to problems with nonsmooth sources that induce discontinuities, our moving mesh, uncollided source method returns a significantly more accurate solution than the standard DG method. On problems with smooth sources, we observe spectral convergence even in problems with wave fronts. In problems where the angular flux is inherently non-smooth, as in Ganapol's (2001) well known plane pulse benchmark, we do not observe an elevated order of accuracy when compared with static meshes, but there is a reduction in error that is nearly three orders of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Benchmarks for infinite medium, time dependent transport problems with isotropic scattering

The widely used AZURV1 transport benchmarks package provides a suite of solutions to isotropic scattering transport problems with a variety of initial conditions (Ganapol 2001). Most of these solutions have an initial condition that is a Dirac delta function in space; as a result these benchmarks are challenging problems to use for verification tests in computer codes. Nevertheless, approximating a delta function in simulation often leads to low orders of convergence and the inability to test the convergence of high-order numerical methods. While there are examples in the literature of integration of these solutions as Green's functions for the transport operator to produce results for more easily simulated sources, they are limited in scope and briefly explained. For a sampling of initial conditions and sources, we present solutions for the uncollided and collided scalar flux to facilitate accurate testing of source treatment in numerical solvers. The solution for the uncollided scalar flux is found in analytic form for some sources. Since integrating the Green's functions is often nontrivial, discussion of integration difficulty and workarounds to find convergent integrals is included. Additionally, our uncollided solutions can be used as source terms in verification studies, in a similar way to the method of manufactured solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Data Reduction in Deterministic Neutron Transport Calculations Using Machine Learning

Neutron cross section matrices for fission and scattering data are required for each material, temperature, and enrichment level to calculate the neutron transport equation accurately. This information can be a limiting factor when using the multigroup discrete ordinates (SN) method when the number of energy groups is large. Machine Learning (ML) can be used to replace the need for the cross section matrices by reproducing the function that maps the scalar flux to the scattering and fission sources. Through the use of autoencoders and Deep Jointly-Informed Neural Networks (DJINN), the data storage requirements are reduced by 94% of the original data for a 618 group problem. This is accomplished while preserving the scalar flux, maintaining generality, and decreasing wall clock times.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Data-Driven Acceleration of Thermal Radiation Transfer Calculations with the Dynamic Mode Decomposition and a Sequential Singular Value Decomposition

We present a method for accelerating discrete ordinates radiative transfer calculations for radiative transfer. Our method works with nonlinear positivity fixes, in contrast to most acceleration schemes. The method is based on the dynamic mode decomposition (DMD) and using a sequence of rank-one updates to compute the singular value decomposition needed for DMD. Using a sequential method allows us to automatically determine the number of solution vectors to include in the DMD acceleration. We present results for slab geometry discrete ordinates calculations with the standard temperature linearization. Compared with positive source iteration, our results demonstrate that our acceleration method reduces the number of transport sweeps required to solve the problem by a factor of about 3 on a standard diffusive Marshak wave problem, a factor of several thousand on a cooling problem where the effective scattering ratio approaches unity, and a factor of 20 improvement in a realistic, multimaterial radiating shock problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Frequency-Dependent Material Motion Benchmarks for Radiative Transfer

We present a general solution for the radiation intensity in front of a purely absorbing slab moving toward an observer at constant speed and with a constant temperature. The solution is obtained by integrating the lab-frame radiation transport equation through the slab to the observer. We present comparisons between our benchmark and results from the Kull simulation code for an aluminum slab moving toward the observer at 2% the speed-of-light. We demonstrate that ignoring certain material motion correction terms in the transport equation can lead to 20-80% errors with the error magnitude growing as the frequency resolution is improved. Our results also indicate that our benchmark can identify potential errors in the implementation of material motion corrections.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Mathematical and numerical validation of the simplified spherical harmonics approach for time-dependent anisotropic-scattering transport problems in homogeneous media

In this work, we extend the solid harmonics derivation, which was used by Ackroyd et al to derive the steady-state SP$_N$ equations, to transient problems. The derivation expands the angular flux in ordinary surface harmonics but uses harmonic polynomials to generate additional surface spherical harmonic terms to be used in Galerkin projection. The derivation shows the equivalence between the SP$_N$ and the P$_N$ approximation. Also, we use the line source problem and McClarren's "box" problem to demonstrate such equivalence numerically. Both problems were initially proposed for isotropic scattering, but here we add higher-order scattering moments to them. Results show that the difference between the SP$_N$ and P$_N$ scalar flux solution is at the roundoff level.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effective Non-oscillatory Regularized L$_1$ Finite Elements for Particle Transport Simulations

In this work, we present a novel regularized L$_1$ (RL$_1$) finite element spatial discretization scheme for radiation transport problems. We review the recently developed least-squares finite element method in nuclear applications. We then derive an L$_1$ finite element by minimizing the L$_1$ norm of the transport residual. To ensure the stability on incident boundary, we newly develop a consistent L$_1$ boundary condition (BC). The numerical tests demonstrate such a method effectively prevents the oscillations which would occur to least-squares finite element when discontinuity exists such as void and absorber. Further, the RL$_1$ method is accurate in problems with scattering.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Globally Conservative, Hybrid Self-Adjoint Angular Flux and Least-Squares Method Compatible with Void

In this paper, we derive a method for the second-order form of the transport equation that is both globally conservative and compatible with voids, using Continuous Finite Element Methods (CFEM). The main idea is to use the Least-Squares (LS) form of the transport equation in the void regions and the Self-Adjoint Angular Flux (SAAF) form elsewhere. While the SAAF formulation is globally conservative, the LS formulation need a correction in void. The price to pay for this fix is the loss of symmetry of the bilinear form. We first derive this Conservative LS (CLS) formulation in void. Second we combine the SAAF and CLS forms and end up with an hybrid SAAF-CLS method, having the desired properties. We show that extending the theory to near-void regions is a minor complication and can be done without affecting the global conservation of the scheme. Being angular discretization agnostic, this method can be applied to both discrete ordinates (SN) and spherical harmonics (PN) methods. However, since a globally conservative and void compatible second-order form already exists for SN (Wang et al. 2014), but is believed to be new for PN, we focus most of our attention on that latter angular discretization. We implement and test our method in Rattlesnake within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework. Results comparing it to other methods are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Implicit Filtered PN for High-Energy Density Thermal Radiation Transport using Discontinuous Galerkin Finite Elements

In this work, we provide a fully-implicit implementation of the time-dependent, filtered spherical harmonics (FPN) equations for non-linear, thermal radiative transfer. We investigate local filtering strategies and analyze the effect of the filter on the conditioning of the system, showing in particular that the filter improves the convergence properties of the iterative solver. We also investigate numerically the rigorous error estimates derived in the linear setting, to determine whether they hold also for the non-linear case. Finally, we simulate a standard test problem on an unstructured mesh and make comparisons with implicit Monte-Carlo (IMC) calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Moment Closures Based on Minimizing the Residual of the P$_N$ Angular Expansion in Radiation Transport

In this work we present two new closures for the spherical harmonics (PN) method in slab geometry transport problems. Our approach begins with an analysis of the squared-residual of the transport equation where we show that the standard truncation and diffusive closures do not minimize the residual of the PN expansion. Based on this analysis we derive two models, a moment-limited diffusive MLDN closure and a transient PN (TPN) closure that attempt to address shortcomings of common closures. The form of these closures is similar to flux-limiters for diffusion with the addition of a time-derivative in the definition of the closure. Numerical results on a pulsed plane source problem, the Gordian knot of slab-geometry transport problems, indicate that our new closure outperforms existing linear closures. Additionally, on a deep penetration problem we demonstrate that the TPN closure does not suffer from the artificial shocks that can arise in the MN entropy-based closure. Finally, results for Reed's problem demonstrate that the TPN solution is as accurate as the PN+3 solution. We further extend the TPN closure to 2D Cartesian geometry. The line source test problem demonstrates the model effectively damps oscillations and negative densities

0 citations0 reviewsNo code linkedOpen access