modelLiquidPhaseSplitC3Model

Liquid phase split with C = 3, p. 22 in preprint

Parameters

TypeNameDefaultDescription
BooleanstartAtSolfalse
Real[:]xstar1{0.6859023, 0.002836175, 0.9548265, 0.3112616, 0.4171931, 0.0006027087}Solution 1
Real[:]xstar2{0.0006027087, 0.9548265, 0.002836175, 0.04457077, 0.5828069, 0.6859023}Solution 2
BooleanincludeTrivialSolutionMetricMessagefalseInclude a warning message when approaching the unwanted trivial solution
Realz10.2Feed composition component 1
Realz20.4Feed composition component 2
Realz30.4Feed composition component 3
RealT298.15Temperature
Realy2_startif startAtSol then xstar1[1] else 0.1
Realy3_startif startAtSol then xstar1[2] else 0.1
Realx3_startif startAtSol then xstar1[3] else 0.1
Realy1_startif startAtSol then xstar1[4] else 0.1
Reallambda_startif startAtSol then xstar1[5] else 0.1
Realx2_startif startAtSol then xstar1[6] else 0.1

Components

TypeNameDefaultDescription
Real[:]x{y2, y3, x3, y1, lambda, x2}
Realerror1Modelica.Math.Vectors.norm(x - xstar1, Modelica.Constants.inf)
Realerror2Modelica.Math.Vectors.norm(x - xstar2, Modelica.Constants.inf)
Realx1Mole fraction of component 1 in phase I
Realx2Mole fraction of component 2 in phase I
Realx3Mole fraction of component 3 in phase I
Realy1Mole fraction of component 1 in phase II
Realy2Mole fraction of component 2 in phase II
Realy3Mole fraction of component 3 in phase II
ReallambdaRelative amount of phase I
Realmetric(x1 - y1)^2 + (x2 - y2)^2Metric to identify trivial solutions, < 1e-4 for probably trivial solutions with z_i = x_i = y_i