Update: refactor JumpConditionsSolver and shockturbulence subpackage to stabilize RH-slope calculations and enable upstream parametric sweeps

[AlbertoCuadra]

Update JumpConditionsSolver

This PR fixes inconsistent inconsistent values of Gammas1 and Gammas3 in JumpConditionsSolver when using scalar Mach inputs (e.g., M1 = 5). The previous implementation relied on interpolation along the Hugoniot curve, which is not well-defined for a single point and can change when extra points are added.

Gammas1 and Gammas3 now enforce the constant post-shock density constraint using the chain-rule, leveraging derivatives with respect to Mach computed during Gammas2:

$$ \left(\frac{\partial p_2}{\partial \phi}\right)_{\rho_2} = \frac{\partial p_2}{\partial \phi} - \frac{\partial p_2}{\partial M_1} \frac{\left(\partial \rho_2/\partial \phi\right)}{\left(\partial \rho_2/\partial M_1\right)} , \qquad \phi \in {T_1, p_1} $$

Results are now stable and consistent for both scalar and array Mach inputs.

Update ShockTurbulenceSolver

Reorganize and generalize the shockturbulence subpackage to improve maintainability and enable parametric sweeps on upstream parameters beyond the pre-shock Mach number.

• solve() now returns mixArray1 and mixArray2, with LIA outputs stored in mixArray2(i).lia as a structure.
• Refactoring and cleanup of solver/model interfaces for easier extension and reuse.

Bugs

• Invoke clean() / cleanMoles() when updating ChemicalSystem.propertiesMatrix to avoid stale values

Minor changes

• Add: eta, chi, etaVorticity, and lia as properties of Mixture class
• Include print results in Mixtureclass for SHOCKTURBULENCE_* problems
• Add: include missing labels in interpreterLabel utility function