Researcher profile

Ryan G. McClarren

Ryan G. McClarren contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 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.

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.

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.

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.

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.

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.