[FEATURE] Add homotopy in pressure loss equation

[pierre-eliep-met]

What is your precise need ?

Improve complex model convergence

What solution do you see ?

Add homotopy in Stodola equation in turbines : Replace

DP_f = - (1+ fouling/100)*Kfr*Q*abs(Q)/rho_in;

By

DP_f = homotopy(- (1+ fouling/100)*Kfr*Q*abs(Q)/rho_in, DP_0/Q_0*Q);

Simplified equation : DP_0/Q_0*Q We should test that it does not break convergence on some complex projects, and add the possibility to use homotopy or not, with an outer parameter