modelWrongInitialSolutionSelected

Initialization converges to the wrong solution

Extends from InitialFailure (Initial nonlinear system of equations has no solutions).

Information

This model describes a simple hydraulic system with a pump, followed by a valve, which fills a reservoir.

The operating point of the pumpt is determined by a nonlinear system with five unknowns: w_pump, dp_pump, dp_valve, sqrt_dp, and p1. They can be reduced to one by selecting w_pump as the tearing variable.

At time t=0, this system has two solutions, one with positive w_pump, and the other one with negative w_pump. If the start value of the tearing variable w_pump is chosen incorrectly, the solver will converge to the negative solution, then lock onto it for the rest of the simulation.

When the user sees the negative w_pump, he/she should be able to analyze how this value was found at time t = 0. The debugger should show that w_pump is solved by that nonlinear system, and show the values of the tearing variables and of the torn variables at each iteration. It will then become apparent that the start value of the teaing variable w_pump, which is negative in this test case, causes the solver to converge on the wrong solution, thus suggesting the solution to the problem, i.e., changing the start value of w_pump to a positive value that allows the solver to converge on the desired solution.

Parameters

TypeNameDefaultDescription
SI.Pressurepatm (from InitialFailure)101325Atmospheric pressure
RealKv (from InitialFailure)1e-2Valve coefficient
Realdp_small (from InitialFailure)1Small dp for valve equation
Realdp0 (from InitialFailure)3e5Pump dp @ zero flow
Reala1 (from InitialFailure)1e6Pump coefficient
Reala2 (from InitialFailure)3e2Pump coefficient
Reala3 (from InitialFailure)3e2Pump coefficient
SI.TemperatureT0 (from InitialFailure)20 + 273.15Temperature of incoming fluid
SI.Densityrho (from InitialFailure)995Density of fluid
SI.AreaA (from InitialFailure)0.01Storage tank cross section
SI.MassFlowRatew_extra (from InitialFailure)0Extra mass flow rate into reservoir
SI.Accelerationg (from InitialFailure)Modelica.Constants.g_nAcceleration of gravity
SI.TemperatureTref (from InitialFailure)273.16Reference temperature for specific enthalpy computation
SI.SpecificHeatCapacitycp (from InitialFailure)4186Cp of the fluid

Components

TypeNameDefaultDescription
SI.MassFlowRatew_pump (from InitialFailure)Mass flow rate from the pump
SI.Pressurep1 (from InitialFailure)Pump discharge pressure
SI.Pressurep2 (from InitialFailure)Storage tank inlet pressure
SI.Pressuredp_pump (from InitialFailure)Pump dp
SI.Pressuredp_valve (from InitialFailure)Valve dp
Realsqrt_dp (from InitialFailure)Regularized sqrt(dp)
SI.SpecificEnthalpyh0 (from InitialFailure)Pump inlet specific enthalpy
SI.SpecificEnthalpyh1 (from InitialFailure)Pump discharge specific enthalpy
SI.PowerW (from InitialFailure)
SI.Lengthy (from InitialFailure)Reservoir level
Realeta (from InitialFailure)(p1 - patm)*w_pump/rho/WPump efficiency
SI.TemperatureT1 (from InitialFailure)Pump discharge temperature
SI.Timetau (from InitialFailure)1Time constant of temperature sensor