3.1.1.4. c3po.couplers.DynamicResidualBalanceCoupler module

Contains the class DynamicResidualBalanceCoupler.

class c3po.couplers.DynamicResidualBalanceCoupler.DynamicResidualBalanceCoupler(physics, exchangers, dataManagers)

Bases: Coupler

DynamicResidualBalanceCoupler inherits from Coupler and proposes a dynamic residual balance algorithm. This is a variant of the adaptive residual balance implemented by c3po.couplers.AdaptiveResidualBalanceCoupler.AdaptiveResidualBalanceCoupler.

This algorithm is designed to couple two solvers using an iterative procedure. It controls the accuracy required to each solver in order to limit over-solving and make them converge together.

See R. Delvaux, “Algorithmes de couplage entre neutronique, thermohydraulique et thermique”, PhD Thesis, Institut Polytechnique de Paris, 2022.

DynamicResidualBalanceCoupler works with :

• Two PhysicsDriver, one for each solver. They must implement the iterateTimeStep() method, together with the possibilities to get residual and set target accuracy.

• Four Exchanger: two for exchanges between the PhysicsDriver, and two to get residuals.

• One LocalDataManager (not just a DataManager) which contains the residuals got with the last two exchangers.

Note

Two Exchanger and a DataManager are used to access the residuals in order to support all possible MPI schemes.

The default target accuracies are 1e-4 and the default maximum number of iterations is 100. Use setConvergenceParameters() to change these values.

It may be interesting to use a FixedPointCoupler to add a damping factor and to control the coupling error.

In this case :

• The option setUseIterate(True) of FixedPointCoupler must be used.

• The maximal number of iterations provided to DynamicResidualBalanceCoupler.setConvergenceParameters() is ignored.

• The exchanger ‘2to1’ is also probably redondant with the FixedPointCoupler exchangers and, in this case, can be set to do nothing.

__init__(physics, exchangers, dataManagers)

Build a DynamicResidualBalanceCoupler object.

Parameters:

• physics (list[PhysicsDriver], dict) – List (or dict with keys ['Solver1', 'Solver2']) of two PhysicsDriver. If a list is used, it has to be provided in the same order than the keys here. The provided PhysicsDriver must implement the iterateTimeStep() method (together with solveTimeStep()) and accept new accuracy (for the solveTimeStep() method) through setInputDoubleValue('Accuracy', value).

• exchangers (list[Exchanger], dict) – List (or dict with keys ['1to2', '2to1', 'Residual1', 'Residual2']) of four Exchanger. If a list is used, it has to be provided in the same order than the keys here.

• dataManagers (list[LocalDataManager], dict) – List (or dict with keys ['Residuals']) of one LocalDataManager (not just a DataManager). The residuals must be stored in this DataManager as double values under the names 'Residual1' and 'Residual2'.

getIterateStatus()

See c3po.PhysicsDriver.PhysicsDriver.getSolveStatus().

getSolveStatus()

See c3po.PhysicsDriver.PhysicsDriver.getSolveStatus().

initTimeStep(dt)

See c3po.PhysicsDriver.PhysicsDriver.initTimeStep().

iterateTimeStep()

See c3po.PhysicsDriver.PhysicsDriver.iterateTimeStep().

setConvergenceParameters(targetResidualSolver1, targetResidualSolver2, maxiter)

Set the convergence parameters (target residuals for each solver and maximum number of iterations).

Parameters:

• targetResidualSolver1 – Target residual for solver 1. Default value: 1.E-4.

• targetResidualSolver2 – Target residual for solver 2. Default value: 1.E-4.

• maxiter – The maximal number of iterations. Default value: 100.

setFailureManagement(leaveIfSolvingFailed)

Set if iterations should continue or not in case of solver failure (solveTimeStep() returns False).

Parameters:

leaveIfSolvingFailed (bool) – Set False to continue the iterations, True to stop. Default: False.

setPrintLevel(level)

Set the print level during iterations (0=None, 1 keeps last iteration, 2 prints every iteration).

Parameters:

level (int) – Integer in range [0;2]. Default: 2.

solveTimeStep()

See c3po.PhysicsDriver.PhysicsDriver.solveTimeStep().