pyopensn.solver.PRKSolver

class pyopensn.solver.PRKSolver

General transient solver for point reactor kinetics.

Wrapper of opensn::PRKSolver.

Advance(self: pyopensn.solver.Solver) None

Advance time values function.

Execute(self: pyopensn.solver.Solver) None

Execute the solver.

GetPeriod(self: pyopensn.solver.PRKSolver) float

Get the period computed for the last time step.

GetPopulationNew(self: pyopensn.solver.PRKSolver) float

Get the population at the next time step.

GetPopulationPrev(self: pyopensn.solver.PRKSolver) float

Get the population at the previous time step.

GetSolutionNew(self: pyopensn.solver.PRKSolver) memoryview

Get the solution at the next time step.

GetSolutionPrev(self: pyopensn.solver.PRKSolver) memoryview

Get the solution at the previous time step.

GetTimeNew(self: pyopensn.solver.PRKSolver) float

Get the time computed for the next time step.

GetTimePrev(self: pyopensn.solver.PRKSolver) float

Get the time computed for the last time step.

Initialize(self: pyopensn.solver.Solver) None

Initialize the solver.

SetRho(self: pyopensn.solver.PRKSolver, rho: float) None

Set the reactivity value for the solver in dollars.

Step(self: pyopensn.solver.Solver) None

Step the solver.

__init__(self: pyopensn.solver.PRKSolver, **kwargs) None

Construct a point reactor kinetics solver.

Parameters:
  • precursor_lambdas (List[float], default=[0.0124, 0.0304, 0.111, 0.301, 1.14, 3.01]) – Decay constants for delayed neutron precursors.

  • precursor_betas (List[float], default=[0.00021, 0.00142, 0.00127, 0.00257, 0.00075, 0.00027]) – Fractional delayed neutron fractions.

  • gen_time (float, default=1.0e-5) – Neutron generation time [s].

  • initial_rho (float, default=0.0) – Initial reactivity in dollars.

  • initial_source (float, default=1.0) – Initial source strength [particles/s].

  • initial_population (float, default=1.0) – Initial neutron population.

  • time_integration (str, default='implicit_euler') – One of {‘explicit_euler’,’implicit_euler’,’crank_nicolson’}.