Error message in fluid samples

[koart81]

Error message appears in fluid samples with the component ScalableTranslationStatistics.Fluid.Model.ScaleTranslationStatistics before simulation: e.g. for the model attached ScalableTranslationStatistics.Fluid.Examples.ScaledNLEquations.NLEquations_5red.zip:

Error in NLEquations_5 (ScalableTranslationStatistics.Fluid.Examples.ScaledNLEquations.NLEquations_5red.mo): Equation system singular (13 variables for 12 equations).
Help: Check.Singular

The structure of the equation system is singular.
Variables: 13
Equations: 12

	Variables		Equations

Real: 13 12 Equations are missing.
Integer: 0 0
Boolean: 0 0
String: 0 0
Variables Equations
continuous-time: 13 12
discrete-time: 0 0
The structure of the equation system is singular.

0: scalableModelicaModel.volume.fluidPortIn.m_flow[2]
1: scalableModelicaModel.simplePressureLossTube.fluidPortOut.m_flow[2]
2: scalableModelicaModel.simplePressureLossTube.fluidPortIn.m_flow[1]
3: scalableModelicaModel.nonlinearTubePumpNetwork.fluidPortOut.m_flow[1]
4: scalableModelicaModel.volume.p[1]
5: scalableModelicaModel.volume.p[2]
6: scalableModelicaModel.volume.p[3]
7: scalableModelicaModel.nonlinearTubePumpNetwork.analytic_pump.fluidPortOut.p[1]
8: scalableModelicaModel.simplePressureLossTube.delta_p[1]
9: scalableModelicaModel.simplePressureLossTube.delta_p[2]
10: der(scalableModelicaModel.volume.m[1])
11: der(scalableModelicaModel.volume.m[2])
12: der(scalableModelicaModel.volume.m[3])

12 Equations [,,,,,]:
0 [1, 0, 0, 0, 1, 0]:
scalableModelicaModel.simplePressureLossTube.fluidPortIn.m_flow[1]=-der(scalableModelicaModel.volume.m[1]);
1 [1, 0, 0, 0, 1, 0]:
scalableModelicaModel.simplePressureLossTube.fluidPortOut.m_flow[2]=-der(scalableModelicaModel.volume.m[3]);
2 [1, 0, 0, 0, 1, 0]: scalableModelicaModel.volume.p[1]=scalableModelicaModel.nonlinearTubePumpNetwork.analytic_pump.fluidPortOut.p[1]+-scalableModelicaModel.nonlinearTubePumpNetwork.analytic_pump.delta_p[1];
3 [1, 0, 0, 0, 1, 0]: scalableModelicaModel.volume.p[1]=scalableModelicaModel.volume.p[2]+scalableModelicaModel.simplePressureLossTube.delta_p[1];
4 [1, 0, 0, 0, 1, 0]: scalableModelicaModel.volume.p[2]=scalableModelicaModel.volume.p[3]+scalableModelicaModel.simplePressureLossTube.delta_p[2];
5 [1, 0, 0, 0, 1, 0]:
der(scalableModelicaModel.volume.m[1])=-scalableModelicaModel.simplePressureLossTube.delta_p[1]* scalableModelicaModel.simplePressureLossTube.delta_p_nom[1]^-1* scalableModelicaModel.simplePressureLossTube.mdot_nom[1];
6 [1, 0, 0, 0, 1, 0]:
der(scalableModelicaModel.volume.m[3])=scalableModelicaModel.simplePressureLossTube.delta_p[2]* scalableModelicaModel.simplePressureLossTube.delta_p_nom[2]^-1* scalableModelicaModel.simplePressureLossTube.mdot_nom[2];
7 [1, 0, 0, 0, 1, 0]:
der(scalableModelicaModel.volume.m[2])=scalableModelicaModel.volume.fluidPortIn.m_flow[2]+-der(scalableModelicaModel.volume.m[3]);
8 [1, 0, 0, 0, 1, 0]:
scalableModelicaModel.volume.p[1]* scalableModelicaModel.volume.V[1]=8.3144626181532395scalableModelicaModel.volume.M[1]^-1 scalableModelicaModel.volume.T[1]* scalableModelicaModel.volume.m[1];
9 [1, 0, 0, 0, 1, 0]:
scalableModelicaModel.volume.p[2]* scalableModelicaModel.volume.V[2]=8.3144626181532395scalableModelicaModel.volume.M[2]^-1 scalableModelicaModel.volume.T[2]* scalableModelicaModel.volume.m[2];
10 [1, 0, 0, 0, 1, 0]:
scalableModelicaModel.volume.p[3]* scalableModelicaModel.volume.V[3]=8.3144626181532395scalableModelicaModel.volume.M[3]^-1 scalableModelicaModel.volume.T[3]* scalableModelicaModel.volume.m[3];
11 [1, 0, 0, 0, 1, 0]: 0=scalableModelicaModel.volume.fluidPortIn.m_flow[2]+scalableModelicaModel.nonlinearTubePumpNetwork.fluidPortOut.m_flow[1]+der(scalableModelicaModel.volume.m[1]);

[koart81]

@MatthiasBSchaefer After commit BugFix: inherit modifiers for Nin and Nout the same error still appear. Could you provide us for the model attached the equations with variable names (from your GSA) or check/compare it yourself with our equations above?

[MatthiasBSchaefer]

Comparison of the failing equation system:

======================
SIMULATION X

1 scalableModelicaModel.simplePressureLossTube.fluidPortIn.m_flow[1]=-der(scalableModelicaModel.volume.m[1]);
2 scalableModelicaModel.simplePressureLossTube.fluidPortOut.m_flow[2]=-der(scalableModelicaModel.volume.m[3]);
3 scalableModelicaModel.volume.p[1]=scalableModelicaModel.nonlinearTubePumpNetwork.analytic_pump.fluidPortOut.p[1]+-scalableModelicaModel.nonlinearTubePumpNetwork.analytic_pump.delta_p[1]
4 scalableModelicaModel.volume.p[1]=scalableModelicaModel.volume.p[2]+scalableModelicaModel.simplePressureLossTube.delta_p[1]
5 scalableModelicaModel.volume.p[2]=scalableModelicaModel.volume.p[3]+scalableModelicaModel.simplePressureLossTube.delta_p[2]
6 der(scalableModelicaModel.volume.m[1])=-scalableModelicaModel.simplePressureLossTube.delta_p[1]* scalableModelicaModel.simplePressureLossTube.delta_p_nom[1]^-1* scalableModelicaModel.simplePressureLossTube.mdot_nom[1];
7 der(scalableModelicaModel.volume.m[3])=scalableModelicaModel.simplePressureLossTube.delta_p[2]* scalableModelicaModel.simplePressureLossTube.delta_p_nom[2]^-1* scalableModelicaModel.simplePressureLossTube.mdot_nom[2];
8 der(scalableModelicaModel.volume.m[2])=scalableModelicaModel.volume.fluidPortIn.m_flow[2]+-der(scalableModelicaModel.volume.m[3]);
9 scalableModelicaModel.volume.p[1]* scalableModelicaModel.volume.V[1]=8.3144626181532395*scalableModelicaModel.volume.M[1]^-1 scalableModelicaModel.volume.T[1]* scalableModelicaModel.volume.m[1];
10 scalableModelicaModel.volume.p[2]* scalableModelicaModel.volume.V[2]=8.3144626181532395*scalableModelicaModel.volume.M[2]^-1 scalableModelicaModel.volume.T[2]* scalableModelicaModel.volume.m[2];
11 scalableModelicaModel.volume.p[3]* scalableModelicaModel.volume.V[3]=8.3144626181532395*scalableModelicaModel.volume.M[3]^-1 scalableModelicaModel.volume.T[3]* scalableModelicaModel.volume.m[3];
12 0=scalableModelicaModel.volume.fluidPortIn.m_flow[2]+scalableModelicaModel.nonlinearTubePumpNetwork.fluidPortOut.m_flow[1]+der(scalableModelicaModel.volume.m[1]);

==> 12 equations for 13 unkwnons

====================
DYMOLA => After causalization

10 scalableModelicaModel.volume[2].p = scalableModelicaModel.volume[2].m*8.31446261815324*scalableModelicaModel.volume[2].T/scalableModelicaModel.volume[2].M/scalableModelicaModel.volume[2].V;
11 scalableModelicaModel.volume[3].p = scalableModelicaModel.volume[3].m*8.31446261815324*scalableModelicaModel.volume[3].T/scalableModelicaModel.volume[3].M/scalableModelicaModel.volume[3].V;
5 scalableModelicaModel.simplePressureLossTube[2].delta_p = scalableModelicaModel.simplePressureLossTube[2].fluidPortIn.p-scalableModelicaModel.simplePressureLossTube[2].fluidPortOut.p; (using aliases)
A scalableModelicaModel.volume[2].fluidPortOut.m_flow = - scalableModelicaModel.simplePressureLossTube[2].mdot_nom*scalableModelicaModel.simplePressureLossTube[2].delta_p/scalableModelicaModel.simplePressureLossTube[2].delta_p_nom;
9 scalableModelicaModel.volume[1].p = scalableModelicaModel.volume[1].m*8.31446261815324*scalableModelicaModel.volume[1].T/scalableModelicaModel.volume[1].M/scalableModelicaModel.volume[1].V;
""""" Nonlinear system of equations:
B scalableModelicaModel.nonlinearTubePumpNetwork[1] : analytic_pump.delta_p := (analytic_pump.n_const/ analytic_pump.n_nom)^2* (1-(fluidPortOut.m_flow*analytic_pump.n_nom/(analytic_pump.mdot_nom*analytic_pump.n_const))^2)*analytic_pump.delta_p_nom;
CscalableModelicaModel.nonlinearTubePumpNetwork[1].parallelResistors.pressureLossTube[1].delta_p:= scalableModelicaModel.volume[1].p+scalableModelicaModel.nonlinearTubePumpNetwork[1].analytic_pump.delta_p-scalableModelicaModel.volume[2].p;
D scalableModelicaModel.nonlinearTubePumpNetwork[1].fluidPortOut.m_flow := - (scalableModelicaModel.nonlinearTubePumpNetwork[1].parallelResistors.pressureLossTube[1].delta_p/scalableModelicaModel.nonlinearTubePumpNetwork[1].parallelResistors.pressureLossTube[1].delta_p_nom )^(1/3)* scalableModelicaModel.nonlinearTubePumpNetwork[1].parallelResistors.pressureLossTube[1].mdot_nom
"""" End nonlinear system of equations
EscalableModelicaModel.nonlinearTubePumpNetwork[1].analytic_pump.fluidPortOut.p = scalableModelicaModel.nonlinearTubePumpNetwork[1].analytic_pump.fluidPortIn.p+scalableModelicaModel.nonlinearTubePumpNetwork[1].analytic_pump.delta_p;
4 scalableModelicaModel.simplePressureLossTube[1].delta_p = scalableModelicaModel.simplePressureLossTube[1].fluidPortIn.p-scalableModelicaModel.simplePressureLossTube[1].fluidPortOut.p; (using aliases)
FscalableModelicaModel.simplePressureLossTube[1].mdot = scalableModelicaModel.simplePressureLossTube[1].mdot_nom*scalableModelicaModel.simplePressureLossTube[1].delta_p/scalableModelicaModel.simplePressureLossTube[1].delta_p_nom;
12 scalableModelicaModel.volume[2].fluidPortIn.m_flow = ( -scalableModelicaModel.simplePressureLossTube[1].fluidPortOut.m_flow)-scalableModelicaModel.nonlinearTubePumpNetwork[1].fluidPortOut.m_flow; (using aliases)
8 der(scalableModelicaModel.volume[2].m) = scalableModelicaModel.volume[2].fluidPortIn.m_flow-( -scalableModelicaModel.volume[2].fluidPortOut.m_flow);
6 der(scalableModelicaModel.volume[1].m) = ( -scalableModelicaModel.nonlinearTubePumpNetwork[1].fluidPortIn.m_flow)-scalableModelicaModel.simplePressureLossTube[1].fluidPortIn.m_flow;

==> equations 1,2,3 and 7 from SimulationX do not appear in Dymola
==> equations A-F from Dymola do not appear in Simulation X
==> 14 equations for 14 unknwons

===============================================
@koart81 There seems to be larger differences in the causalization process between the tools
I'll send you the whole flattened and causalized model code via mail
Is the problem specific for this parametrization of the model or does it appear more often ?

[koart81]

@MatthiasBSchaefer the problem appears in all models of fluid examples: ScalableTranslationStatistics.Fluid.Examples