Source author record

Barry D. Ganapol

Barry D. Ganapol appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.comp-ph Biological Physics math-ph math.MP Neural and Evolutionary Computing

Catalog footprint

What is connected

5works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2025arXiv

Evolutionary Discovery of Sequence Acceleration Methods for Slab Geometry Neutron Transport

We present a genetic programming approach to automatically discover convergence acceleration methods for discrete ordinates solutions of neutron transport problems in slab geometry. Classical acceleration methods such as Aitken's delta-squared and Wynn epsilon assume specific convergence patterns and do not generalize well to the broad set of transport problems encountered in practice. We evolved mathematical formulas specifically tailored to SN convergence characteristics in this work. The discovered accelerator, featuring second differences and cross-product terms, achieved over 75 percent success rate in improving convergence compared to raw sequences - almost double that observed for classical techniques for the problem set considered. This work demonstrates the potential for discovering novel numerical methods in computational physics via genetic programming and attempts to honor Prof. Ganapol's legacy of advancing experimental mathematics applied to neutron transport.

preprint2020arXiv

Matrix riccati equation solution of the 1d radiative transfer equation

In recent years, the first author has developed three successful numerical methods to solve the 1D radiative transport equation yielding highly precise benchmarks. The second author has shown a keen interest in novel solution methodologies and an ability for their implementation. Here, we combine talents to generate yet another high precision solution, the Matrix Riccati Equation Method (MREM). MREM features the solution to two of the four matrix Riccati ODEs that arise from the interaction principle of particle transport. Through interaction coefficients, the interaction principle describes how particles reflect from- and transmit through- a single slab. On combination with Taylor series and doubling, a high quality numerical benchmark, to nearly seven places, is established.

preprint2016arXiv

Particle Transport in a 3D duct by adding and doubling

Particle transport through a duct by Lambertian reflection from duct walls is again considered. This popular transport example has been solved by most numerical transport methods except notably one- the method of doubling. We shall show that the method of doubling provides every bit as, or more, accurate reflectances and transmittances as the numerical discrete ordinates (NDO) and analytical discrete ordinates (ADO) methods with less mathematical and numerical effort.

preprint2014arXiv

The Fourier transform solution for the Green's function of monoenergetic neutron transport theory

Nearly 45 years ago, Ken Case published his seminal paper on the singular eigenfunction solution for the Green's function of the monoenergetic neutron transport equation with isotropic scattering. Previously, the solution had been obtained by Fourier transform. While it is apparent the two had to be equivalent, a convincing equivalence proof for general anisotropic scattering remained a challenge until now.

preprint2010arXiv

An Accurate Numerical Solution to the Kinetics of Breakable Filament Assembly

Proteinaceous aggregation occurs through self-assembly-- a process not entirely understood. In a recent article [1], an analytical theory for amyloid fibril growth via secondary rather than primary nucleation was presented. Remarkably, with only a single kinetic parameter, the authors were able to unify growth characteristics for a variety of experimental data. In essence, they seem to have uncovered the underlying allometric laws governing the evolution of filament elongation simply from two coupled non-linear ordinary differential equations (ODEs) stemming from a master equation. While this work adds significantly to our understanding of filament self-assembly, it required an approximate analytical solution representation. Here, we show that the same results are found by purely numerical means once a straightforward and reliable numerical solution to the set of ODEs has been established.