Search

While the Boltzmann transport equation can accurately model transport problems with highly forward-peaked scattering, obtaining its solution can become arbitrarily slow due to near-unity spectral radius associated with source iteration. Standard acceleration techniques like diffusion synthetic acceleration and nonlinear diffusion acceleration obtain merely one order of magnitude speedups compared to source iteration due to slowly decaying error moments. Additionally, converging approximations to the Boltzmann equation like Fokker-Planck and Boltzmann Fokker Planck run into similar problems with slow convergence. In this paper we assess the feasibility of using Fourier neural operators to obtain AI-enhanced low order, and single-sweep solutions for the transport equation in slab geometry using a predictor-corrector framework., Comment: Paper draft has been submitted to M&C 2025

The convergence of Boltzmann Fokker Planck solution can become arbitrarily slow with iterative procedures like source iteration. This paper derives and investigates a nonlinear diffusion acceleration scheme for the solution of the Boltzmann Fokker Planck equation in slab geometry. This method is a conventional high order low order scheme with a traditional diffusion-plus-drift low-order system. The method, however, differs from the earlier variants as the definition of the low order equation, which is adjusted according to the zeroth and first moments of the Boltzmann Fokker Planck equation. For the problems considered, we observe that the NDA-accelerated solution follows the unaccelerated well and provides roughly an order of magnitude savings in iteration count and runtime compared to source iteration., Comment: This draft has been submitted to Physor 2024 for review

A Python interface is developed for the GPWR Simulator to automatically simulate cyber-spoofing of different steam generator parameters and plant operation. Specifically, steam generator water level, feedwater flowrate, steam flowrate, valve position, and steam generator controller parameters, including controller gain and time constant, can be directly attacked using command inject, denial of service, and man-in-the-middle type attacks. Plant operation can be initialized to any of the initial conditions provided by the GPWR simulator. Several different diagnostics algorithms have been implemented for anomaly detection, including physics-based diagnostics with Kalman filtering, data-driven diagnostics, noise profiling, and online sensor validation. Industry-standard safety analysis code RELAP5 is also available as a part of the toolkit. Diagnostics algorithms are analyzed based on accuracy and efficiency. Our observations indicate that physics-based diagnostics with Kalman filtering are the most robust. An experimental quantum kernel has been added to the framework for preliminary testing. Our first impressions suggest that while quantum kernels can be accurate, just like any other kernels, their applicability is problem/data dependent, and can be prone to overfitting., Comment: Paper has been submitted to ANS for review

The nonclassical transport equation models particle transport processes in which the particle flux does not decrease as an exponential function of the particle's free-path. Recently, a spectral approach was developed to generate nonclassical spectral S$_N$ equations, which can be numerically solved in a deterministic fashion using classical numerical techniques. This paper introduces a transport synthetic acceleration procedure to speed up the iteration scheme for the solution of the monoenergetic slab-geometry nonclassical spectral S$_N$ equations. We present numerical results that confirm the benefit of the acceleration procedure for this class of problems., Comment: 14 pages, 2 figures

The discrete ordinates method can model forward-peaked transport problems accurately. However, convergence of discrete ordinates solution can become arbitrarily slow upon use of standard iterative procedures like source iteration and GMRES. Standard zeroth and first moment-based acceleration methods like nonlinear diffusion acceleration and diffusion synthetic acceleration are ineffective in accelerating such problems because these methods do not correct higher order Legendre-moments of angular flux. We explore the idea of using Fokker-Planck as a preconditioner to accelerate forward-peaked transport problems in this paper., Comment: This is a draft for our FPSA paper. It will be edited further and submitted to a journal soon

Therapies such as combined-hyperthermia-radiotherapy (CHR) take advantage of excellent radiosensitization properties of hyperthermia and treat of tumors with both radiation and heat. To appropriately model a CHR treatment, features like tumor heating (heat transfer), dosimetry (radiation transport), and tumor dynamics (cell population dynamics) must be considered together. Our first paper on modeling such treatments introduced a one-way coupled model that could only account for isotropic scattering of radiation. In the present work, we extend our transport model to incorporate anisotropic scattering., Comment: 4 pages, 2 figures, submitted to ANS Winter Meeting and Expo. arXiv admin note: substantial text overlap with arXiv:1905.10441

The use of machine learning algorithms to address classification problems is on the rise in many research areas. The current study is aimed at testing the potential of using such algorithms to auto-select the best solvers for transport problems in uniform slabs. Three solvers are used in this work: Richardson, diffusion synthetic acceleration, and nonlinear diffusion acceleration. Three parameters are manipulated to create different transport problem scenarios. Five machine learning algorithms are applied: linear discriminant analysis, K-nearest neighbors, support vector machine, random forest, and neural networks. We present and analyze the results of these algorithms for the test problems, showing that random forest and K-nearest neighbors are potentially the best suited candidates for this type of classification problem., Comment: Accepted -- The International Conference on Mathematics and Computational Methods applied to Nuclear Science and Engineering (M&C 2019), Portland, OR

We develop a multiphysics-based model to predict the response of localized tumors to combined-hyperthermia-radiotherapy (CHR) treatment. This procedure combines hyperthermia (tumor heating) with standard radiotherapy to improve efficacy of the overall treatment. In addition to directly killing tumor cells, tumor heating amends several parameters within the tumor microenvironment. This leads to radiosensitization, which improves the performance of radiotherapy while reducing the side-effects of excess radiation in the surrounding normal tissue. Existing tools to model this kind of treatment consider each of the physics separately. The model presented in this paper accounts for the synergy between hyperthermia and radiotherapy providing a more realistic and holistic approach to simulate CHR treatment. Our model couples radiation transport and heat-transfer with cell population dynamics., Comment: 10 pages, 3 figures, submitted to ANS Topical Meeting