Add multiple reactions in particles

[AntoniaBerger]

This pull request introduces a new reaction interface to support multiple reactions per particle type.
The changes include adding a new reaction model (MassActionLawReactionCrossPhase), updating the reaction model factory, and modifying the GeneralRateModel and LumpedRateModelWithPores to handle the new interface. Below are the key changes grouped by theme:

New Reaction Model

• Added MassActionLawReactionCrossPhase to the build system in CMakeLists.txt and registered it in the ReactionModelFactory. Due to backwards compatibility, the old MAL don't loos any feature.

GeneralRateModel 1D Enhancements

• Introduced support for multiple reactions per particle type in GeneralRateModel. This includes a new _oldReactionInterface flag to distinguish between the old and new interfaces, and _numReactionsPerParticle to track the number of reactions per particle type. [1] [2]
• Updated configureModelDiscretization to initialize and configure multiple reactions for each particle type.
• Enhanced configure to handle the new reaction interface, ensuring proper configuration of multiple reactions. [1] [2]
• Adjusted memory allocation in threadLocalMemorySize to account for multiple reactions. [1] [2]
• Modified residualParticle to compute residuals and Jacobians for multiple reactions per particle type.

LumpedRateModelWithPores 1D Enhancements

• Similar updates as GeneralRateModel to support multiple reactions per particle type, including _oldReactionInterface and _numReactionsPerParticle.
• Updated configureModelDiscretization and configure methods to handle the new reaction interface.
• Adjusted threadLocalMemorySize to allocate memory for multiple reactions.

CSTR Enhancements

• Similar updates as GeneralRateModel to support multiple reactions per particle type, including _oldReactionInterface and _numReactionsPerParticle.
• Updated configureModelDiscretization and configure methods to handle the new reaction interface.
• Adjusted threadLocalMemorySize to allocate memory for multiple reactions.

These changes enhance the flexibility of the reaction models, allowing for more complex scenarios with multiple reactions per particle type.

[AntoniaBerger]

The new suggested user interface for cross phase reactions will be similar to the interface for bulk reactions:

model.root.input.model.unit_000.particle_type_000.n_cross_phase_reaction = 1
model.root.input.model.unit_000.particle_type_000.cross_phase_reaction_000.type = "MASS_ACTION_LAW_CROSS_PHASE"
model.root.input.model.unit_000.particle_type_000.cross_phase_reaction_000.mal_stoich_liquid = [-1, 0]
model.root.input.model.unit_000.particle_type_000.cross_phase_reaction_000.mal_stoich_solid = [0, 1]
model.root.input.model.unit_000.particle_type_000.cross_phase_reaction_000.mal_kfwd_liquid = [1.0]
model.root.input.model.unit_000.particle_type_000.cross_phase_reaction_000.mal_kbwd_solid = [0.5]

Up to discussion is still, how we want to continue with parameters like MAL_EXPONENTS_LIQUID_FWD_MODSOLID.

And for reactions located only in the pore (liquid) phase I purpose:

model.root.input.model.unit_000.particle_type_000.n_reaction_pore= 1
model.root.input.model.unit_000.particle_type_000.pore_reaction_000.type = "MASS_ACTION_LAW"
model.root.input.model.unit_000.particle_type_000.pore_reaction_000.mal_stoich = [-1, 0]
model.root.input.model.unit_000.particle_type_000.pore_reaction_000.mal_kfwd = [1.0]
model.root.input.model.unit_000.particle_type_000.pore_reaction_000.mal_kbwd = [0.5]

To use the same reaction modul for different phases, I gave up the _liquid / _solid description.
I also prefer reaction_cross_phase, reaction_bulk and reaction_solid. Since this is more general and also is valid for CSTR and LRM where there is no particle phase.
For more information, see #92