Avoid the tank to ever become completely empty

modelica · Feb 9, 2025 · 64ec590 · 64ec590
1 parent ff95fea
commit 64ec590
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/Modelica/Fluid/Vessels.mo b/Modelica/Fluid/Vessels.mo
@@ -307,7 +307,7 @@ of the modeller. Increase nPorts to add an additional port.
         end for;
         // Check for correct solution
         assert(fluidLevel <= fluidLevel_max, "Vessel is overflowing (fluidLevel > fluidLevel_max = " + String(fluidLevel) + ")");
-        assert(fluidLevel > -1e-6*fluidLevel_max, "Fluid level (= " + String(fluidLevel) + ") is below zero meaning that the solution failed.");
+        assert(fluidLevel > 1e-6*fluidLevel_max, "Fluid level (= " + String(fluidLevel) + ") is too small, the tank cannot be completely empty otherwise the energy balance equation becomes singular.");
 
         // Boundary conditions
 

diff --git a/ModelicaTest/Fluid/TestComponents/Vessels/TestSimpleTank.mo b/ModelicaTest/Fluid/TestComponents/Vessels/TestSimpleTank.mo
@@ -10,10 +10,10 @@ model TestSimpleTank
     height=1,
     nPorts=2,
     portsData={Modelica.Fluid.Vessels.BaseClasses.VesselPortsData(diameter=0.1,
-        height=0),Modelica.Fluid.Vessels.BaseClasses.VesselPortsData(diameter=
+        height=0.001),Modelica.Fluid.Vessels.BaseClasses.VesselPortsData(diameter=
         0.1, height=1)},
     crossArea=1,
-    level_start=0)
+    level_start=0.001)
     annotation (Placement(transformation(extent={{0,0},{40,40}})));
 
   inner Modelica.Fluid.System system(energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial)