More accurately simulate the neutron flux profile#3653

Draft

OceanNuclear wants to merge 99 commits intomainfrom

Commits on Oct 28, 2025

First commit - made a mockup of NeutronFluxProfile calculator.
OceanNuclear
committed
Added the effect of neutron multiplication, which is done by increasing the strength of the macroscopic (self and other) scattering matrix.
OceanNuclear
committed
relocated some functions
OceanNuclear
committed
Ruff improvements.
OceanNuclear
committed
Decided to implement a switch case for sinh(x/sqrt(-...))=-sin(x/sqrt(...)), instead of using complex number to cover both +ve and -ve sign in the sqrt.
OceanNuclear
committed
Some further progress on docstrings and equations.
OceanNuclear
committed
WIP - pushing to save progress.
OceanNuclear
committed
Transcribed the one-group constant (manually solved and verified by Desmos). Next step is to write the program for solving multi-group fluxes.
OceanNuclear
committed
Modified expression into sum-of-exponentials form, rather than sum-of-cosh+sinh.
OceanNuclear
committed
Added reaction rates calculations and neutron passage/escape flux.
OceanNuclear
committed
Added a unit test.
OceanNuclear
committed

Commits on Oct 29, 2025

Commits on Oct 30, 2025

Commits on Nov 5, 2025

Created the scattering matrix, and tested it. Next on the to-do list: n2n matrix.
OceanNuclear
committed

Commits on Nov 6, 2025

Commits on Nov 7, 2025

Commits on Nov 9, 2025

Commits on Nov 12, 2025

Commits on Nov 21, 2025

Commits on Nov 23, 2025

Commits on Dec 16, 2025

Commits on Dec 17, 2025

Commits on Feb 1, 2026

Handling duplicate diffusion lengths is simpler than expected!
OceanNuclear
committed