Classes
Here you can get necessary information about how the elements in Dyssol (like unit, solver, PSD function, matrices and external solvers) look like on a programming point of view. This is especially important if you want to develope new units and solvers.
All elements mentioned above are defined as classes, which consists of member variables (data fields) and associated functions (methods). Every new independent instance generated under a class is called an object.
For more information on objectoriented programming applied in Dyssol, please refer to internet sources about C++.
Basic unit
Functions in class BaseUnit
are introduced below.
Note
Some definitions:
Compound: a chemical substance defined by particular set of physical properties, calculation methods and data. Examples: water, hydrogen, oxygen.
Phase: stable collection of compounds with a defined amount of substance and a homogeneous composition. It has an associated state of aggregation, e.g. liquid.
State of aggregation: the physical state in which the compounds in that phase occur. Possible state of aggregation: vapor, liquid, solid.
Basic information
std::string GetUnitName()
Returns name of the unit.
std::string GetAuthorName()
Returns name of the unit’s author.
std::string GetUniqueID()
Returns unique identifier of the unit.
std::string GetUnitVersion()
Returns version of the unit.
Ports
unsigned AddPort(std::string PortName, unsigned PortType)
Adds port with PortType
(INPUT_PORT
, OUTPUT_PORT
or UNDEFINED_PORT
) to the unit. Returns index of the port. Should be used in unit’s constructor only. PortName
should be unique within the unit.
unsigned GetPortsNumber()
Returns number of ports in the unit.
unsigned GetPortType(std::string PortName)
Returns type of the port with name PortName. Returning values: INPUT_PORT
, OUTPUT_PORT
. If port with such name has not been defined , UNDEFINED_PORT
will be returned.
CMaterialStream* GetPortStream(std::string PortName)
Returns pointer to the stream which is connected to the port with name PortName
. Returns NULL
if such port has not been defined.
Material streams and holdups
CHoldup* AddHoldup(std::string HoldupName, std::string StreamKey = "")
Adds new holdup with the specified name HoldupName to the unit. HoldupName
should be unique within the unit. The structure of the holdup will be the same as the global stream’s structure (phases, grids, compounds etc.).
Should be used in unit’s constructor; then the holdup will be automatically handled by the simulation system (saved and loaded during the simulation, cleared and removed after use). However, it is allowed to add holdups outside the constructor for temporal purposes, but afterwards you have to save, load or remove this holdup manually by calling functions SaveState()
, LoadState()
and RemoveHoldup()
. Otherwise, all such holdups will be removed at the end of the simulation.
Returns pointer to a created holdup. This pointer should not be used inside the constructor of the unit, since all changes of the holdup made through this pointer will be cancelled during the initialization of the unit.
CHoldup* GetHoldup(std::string HoldupName)
Returns pointer to a holdup with the specified name HoldupName
. This pointer should not be used inside the constructor of the unit, since all changes of the holdup made through this pointer will be cancelled during the initialization of the unit. NULL
will be returned if such holdup has not been defined.
std::vector<CHoldup*> GetHoldups()
Returns pointers to all holdups of the unit. These pointers should not be used inside the constructor of the unit, since all changes of the holdup made through them will be cancelled during the initialization of the unit.
void RemoveHoldup(std::string HoldupName)
Removes holdup with the specified name HoldupName from the unit. Should be used only for those holdups, which have been added to the unit outside the constructor.
CMaterialStream* AddMaterialStream(std::string StreamName, std::string StreamKey = "")
Adds new material stream with the specified name StreamName
for internal temporary use to the unit. StreamName
should be unique within the unit. Structure of the stream will be the same as the global stream’s structure (phases, grids, compounds etc.).
Should be used in unit’s constructor; then the material stream will be automatically handled by the simulation system, this means saved and loaded during the simulation, cleared and removed after use. However, it is allowed to add material outside the constructor for temporal purposes, but you have to save, load or remove this stream manually by calling functions SaveState()
, LoadState()
and RemoveMaterialStream()
. Otherwise, all such material streams will be removed at the end of the simulation.
Returns pointer to a created material stream.
CMaterialStream* GetMaterialStream(std::string StreamName)
Returns pointer to a material stream with specified name StreamName
. NULL
will be returned if such stream has not been defined.
void RemoveMaterialStream(std::string StreamName)
Removes material stream with the specified name StreamName
from the unit. Should be used only for those material streams, which have been added to the unit outside the constructor.
CMaterialStream* AddFeed(std::string FeedName, std::string StreamKey = "")
Adds new feed stream with the name FeedName
to the unit. FeedName
should be unique within the unit. The structure of the stream will be the same as the global stream structure (phases, grids, compounds etc.).
Should be used only in constructor of the unit describing the feed. Returns pointer to a created stream. This pointer should not be used inside the constructor of the unit, since all changes of the stream made through this pointer will be cancelled during the initialization of the unit.
CMaterialStream* GetFeed(std::string FeedName)
Returns pointer to a feed with specified name FeedName
. This pointer should not be used inside the constructor of the unit, since all changes of the stream made through this pointer will be cancelled during the initialization of the unit. NULL
will be returned if such feed has not been defined.
void CopyStreamToStream(CMaterialStream* SrcStream, CMaterialStream DstStream, double Time, bool DeleteDataAfter = true)
Copies all data from SrcStream
to DestStream
for specified time point. If flag DeleteDataAfter
is set to true, all data after the time point Time
in the destination stream will be removed before being copied.
void CopyStreamToStream(CMaterialStream* SrcStream, CMaterialStream DstStream, double StartTime, double EndTime, bool DeleteDataAfter = true)
Copies all data from SrcStream
to DstStream
on a specified time interval. If flag DeleteDataAfter
is set to true, all data after the time point StartTime
in the destination stream will be removed before being copied.
void CopyStreamToPort(CMaterialStream* Stream, std::string PortName, double Time, bool DeleteDataAfter = true)
Copies all data from Stream
to a stream connected to an output port PortName
for specified time point. If flag DeleteDataAfter
is set to true, all data after the time point Time
in the destination stream will be removed before being copied.
void CopyStreamToPort(CMaterialStream* Stream, std::string PortName, double StartTime, double EndTime, bool DeleteDataAfter = true)
Copies all data from Stream
to a stream connected to an output port PortName
for specified time interval. If flag DeleteDataAfter
is set to true, all data after the time point StartTime
in the destination stream will be removed before being copied.
void CopyPortToStream(std::string PortName, CMaterialStream* Stream, double Time, bool DeleteDataAfter = true)
Copies stream’s data of the input port PortName
to another stream Stream
for specified time point. If flag DeleteDataAfter
is set to true, all data after the time point Time in the destination stream will be removed before being copied.
void CopyPortToStream(std::string PortName, CMaterialStream* Stream, double StartTime, double EndTime, bool DeleteDataAfter = true)
Copies stream’s data of the input port PortName
to another stream Stream
for specified time interval. If flag DeleteDataAfter
is set to true, all data after the time point StartTime
in the destination stream will be removed before being copied.
double CalcTemperatureFromProperty(ECompoundTPProperties Property, std::vector<double> CompoundFractions, double Value)
Returns temperature of a generic system of composition CompoundFractions
for a specific value Value
of the property Property
. Possible properties are those defined in material database. For more information, please refer to Thermodynamics on this page.
double CalcPressureFromProperty(ECompoundTPProperties Property, std::vector<double> CompoundFractions, double Value)
Returns pressure of a generic system of composition CompoundFractions
for a specific value Value
of the property Property
. Possible properties are those defined in material database. For more information, please refer to Thermodynamics on this page.
void HeatExchange(CMaterialStream* Stream1, CMaterialStream* Stream2, double Time, double Efficiency);
Performs a heat exchange between material streams Stream1
and Stream2
at specified time point Time
with a specified efficiency Efficiency
ranging between 0 and 1. For more information, please refer to Thermodynamics on this page.
Time points
std::vector<double> GetAllInputTimePoints(double StartTime, double EndTime, bool ForceStartBoundary = false, bool ForceEndBoundary = false)
Returns all time points on which input streams of the unit are defined for specified time interval. Input streams are all streams connected to the input ports. If ForceStartBoundary
and/or ForceEndBoundary
flag is enabled, corresponding boundary points will be forcibly added to the resulting vector.
std::vector<double> GetAllDefinedTimePoints(double StartTime, double EndTime, bool ForceStartBoundary = false, bool ForceEndBoundary = false)
Returns all time points for specified time interval on which input streams and unit parameters are defined. Input streams are all streams connected to the input ports. If ForceStartBoundary
and/or ForceEndBoundary
flag is enabled, corresponding boundary points will be forcibly added to the resulting vector.
std::vector<double> GetAllStreamsTimePoints(std::vector<CMaterialStream*> Srteams, double StartTime, double EndTime)
Returns all time points for specified time interval on which Stream
s are defined.
Unit parameters
unsigned AddConstParameter(std::string Name, double MinValue, double MaxValue, double InitValue, std::string Description = "")
Adds new constant unit parameter to the unit with the name Name
, boundary values MinValue
and MaxValue
, description Description
, and initializes it with the value InitValue
. The name of the parameter should be unique within the unit.
Should be used in unit’s constructor only. Returns index of the parameter.
unsigned AddTDParameter(std::string Name, double MinValue, double MaxValue, double InitValue, std::string Description = "")
Adds new timedependent unit parameter to the unit with the name Name
, boundary values MinValue
and MaxValue
, description Description
, and initializes it in the time point 0 s with the value InitValue
. The name of the parameter should be unique within the unit.
Should be used in unit’s constructor only. Returns index of the parameter.
unsigned AddStringParameter(std::string Name, std::string InitValue = "", std::string Description = "")
Adds new string unit parameter to the unit with the name Name
, description Description
, and initializes it with the value InitValue
. The name of the parameter should be unique within the unit.
Should be used in unit’s constructor only. Returns index of the parameter.
unsigned AddSolverAggregation(std::string Name, std::string Description = "")
Adds new solver parameter of type SOLVER_AGGREGATION_1
with name Name
and description Description
. Allows choosing a specific solver with current type to use it in unit. The name of the parameter should be unique within the unit.
Should be used in unit’s constructor only. Returns index of the parameter.
unsigned GetParametersNumber()
Returns number of unit parameters which have been defined in the unit.
std::string GetParameterName(unsigned Index)
Returns name of the unit parameter with the specified index Index. Empty string is returned if such parameter has not been defined.
double GetParameterMinVal(std::string ParameterName)
Returns minimum allowable value of the timedependent or constant unit parameter with the name ParameterName
. Returns 0
if such parameter has not been defined or this is not a constant or timedependent parameter.
double GetParameterMaxVal(std::string ParameterName)
Returns maximum allowable value of the timedependent or constant unit parameter with the name ParameterName
. Returns 0
if such parameter has not been defined or this is not a constant or timedependent parameter.
double GetConstParameterValue(std::string ParameterName)
Returns value of a constant unit parameter with the name ParameterName
. Returns 0
if such constant parameter has not been defined.
double GetTDParameterValue(std::string ParameterName, double Time)
Returns value of a timedependent unit parameter with the name ParameterName
at the specified time point Time
. If the parameter is not defined at Time
, linear interpolation will be used to obtain the value. Returns 0
if such timedependent parameter has not been defined.
std::string GetStringParameterValue(std::string ParameterName)
Returns value of a string unit parameter. If such string parameter has not been defined, empty string will be returned.
CAggregationSolver* GetSolverAggregation(std::string ParameterName)
Returns pointer to a chosen solver of SOLVER_AGGREGATION_1
type. Returns nullptr
if such unit parameter has not been defined.
void SetParameterValue(std::string ParameterName, double Value, double Time = 0)
Sets Value
of a constant or a time dependent unit parameter in the specified time point Time
. If the time point does not exist, it will be created. If the time point already exists, its value will be overwritten.
void SetParameterValue(std::string ParameterName, std::string Value)
Sets new Value
of a string unit parameter with the specified name ParameterName
.
std::vector<double> GetParameterTimePoints(std::string ParameterName, double TimeStart, double TimeEnd)
Returns all time points for which timedependent parameter is defined within the specified time interval [TimeStart
; TimeEnd
]. Returns empty vector if such parameter has not been defined or this is not a timedependent parameter.
State variables
unsigned AddStateVariable(std::string Name, double InitValue, bool SaveHistory = false)
Adds new state variable and initializes it with InitValue
. Name must be unique within the unit’s state variables. Parameter SaveHistory
specifies if the history of all changes of variable should be saved during calculation for further postprocessing. To save history, function SaveStateVariables()
should be called. All variables which are added with this function will be automatically saved and restored during the simulation. Should be used in Initialize
function of the unit.
Returns index of added variable.
unsigned GetStateVariablesNumber()
Returns number of state variables which have been defined in this unit.
std::string GetStateVariableName(unsigned Index)
Returns the name of the state variable with specified index. Returns empty string if such variable has not been defined.
double GetStateVariable(std::string Name)
Returns value of internal state variable with name Name. Returns 0
if such variable has not been defined.
void SetStateVariable(std::string Name, double Value)
Sets new value Value
of a state variable Name
.
void ClearStateVariables()
Removes all state variables and history of their changes.
void SaveStateVariables(double Time)
Saves values of those internal variables defined as having history at the current time Time
.
Compounds
std::vector<std::string> GetCompoundsList()
Returns unique keys of all compounds defined in the current flowsheet.
std::vector<std::string> GetCompoundsNames()
Returns names of all compounds defined in the current flowsheet.
unsigned GetCompoundsNumber()
Returns number of compounds which are defined in the current flowsheet.
double GetCompoundConstant(std::string CompoundKey, unsigned Constant)
Returns value of constant physical property for specified compound. These properties are stored in the database of materials. Possible constants are listed below:
SOA_AT_NORMAL_CONDITIONS
NORMAL_BOILING_POINT
NORMAL_FREEZING_POINT
CRITICAL_TEMPERATURE
CRITICAL_PRESSURE
MOLAR_MASS
STANDARD_FORMATION_ENTHALPY
HEAT_OF_FUSION_AT_NORMAL_FREEZING_POINT
HEAT_OF_VAPORIZATION_AT_NORMAL_BOILING_POINT
REACTIVITY_TYPE
CONST_PROP_USER_DEFINED_XX
Note
Definition of constant properties:
Definition in class 
Name 
Units 


State of aggregation at normal conditions 
[] 

Normal boiling point 
[K] 

Normal freezing point 
[K] 

Critical temperature 
[K] 

Critical pressure 
[Pa] 

Molar mass 
[kg/mol] 

Standard formation enthalpy 
[J/mol] 

Heat of fusion at normal freezing point 
[J/mol] 

Heat of vaporization at normal boiling point 
[J/mol] 

Reactivity type 
[] 

User defined property 
[] 
double GetCompoundTPDProp(std::string CompoundKey, unsigned Property, double Temperature, double Pressure)
Returns value of temperature/pressuredependent physical properties stored in the material database for compound with specified Temperature
[K] and Pressure
[Pa]. Possible properties are listed below:
DENSITY
HEAT_CAPACITY_CP
ENTHALPY
VAPOR_PRESSURE
VISCOSITY
THERMAL_CONDUCTIVITY
PERMITTIVITY
TP_PROP_USER_DEFINED_XX
Note
Definition of temperaturedependent properties:
Define 
Name 
Units 


Density 
[kg/m \(^3\)] 

Heat 
[J/(kg·K)] 

Enthalpy 
[J/kg] 

Vapor 
[Pa] 

Viscosity 
[Pa·s] 

Thermal 
[W/(m·K)] 

Permittivity 
[F/m] 

User 
[] 
double GetCompoundsInteractionProp(std::string CompoundKey1, std::string CompoundKey2, unsigned Property, double Temperature, double Pressure)
Returns the value of the interaction property for selected compounds under the specified Temperature
[K] and Pressure
[Pa]. These properties are stored in the database of materials. Possible properties are listed below:
INTERFACE_TENSION
INT_PROP_USER_DEFINED_XX
bool IsCompoundNameDefined(std::string CompoundName)
Returns true
if compound with specified name has been defined, otherwise returns false
.
bool IsCompoundKeyDefined(std::string CompoundKey)
Returns true
if compound with specified unique key has been defined, otherwise returns false
.
Phases
unsigned GetPhasesNumber()
Returns number of phases which are currently defined in the flowsheet.
std::string GetPhaseName(unsigned PhaseType)
Returns name of the specified phase. Possible values of PhaseType
are:
SOA_SOLID
SOA_LIQUID
SOA_LIQUID2
SOA_VAPOR
Returns empty string if such phase has not been defined.
unsigned GetPhaseSOA(unsigned PhaseIndex)
Returns state of aggregation for the phase with index PhaseIndex
. Returns SOA_UNDEFINED
if the phase with specified index doesn’t exist.
size_t GetPhaseIndex(unsigned PhaseType)
Returns index of the specified phase. Returns 1
if such phase has not been defined. Possible values of PhaseType are:
SOA_SOLID
SOA_LIQUID
SOA_LIQUID2
SOA_VAPOR
bool IsPhaseDefined(unsigned PhaseType)
Returns true
if such phase has been defined in the flowsheet. Possible values of PhaseType are:
SOA_SOLID
SOA_LIQUID
SOA_LIQUID2
SOA_VAPOR
unsigned GetLiquidPhasesNumber()
Returns number of defined liquid phases.
Solid distributed properties
std::vector<EDistrTypes> GetDistributionsTypes()
Returns list of types of solid distributions, which have been defined in the flowsheet. Possible types are:
DISTR_COMPOUNDS
: distribution by compounds.
DISTR_SIZE
: distribution by particle size.
DISTR_PART_POROSITY
: distribution by partcle porosity.
DISTR_FORM_FACTOR
: distribution by particle form factor.
DISTR_COLOR
: distribution by particle color.
DISTR_USER_DEFINED_01
toDISTR_USER_DEFINED_10
: user defined distribution.
std::vector<unsigned> GetDistributionsClasses()
Returns list with number of classes for all defined distributions.
unsigned GetDistributionsNumber()
Returns number of solids distributions, which have been defined in the flowsheet.
EGridTypes GetDistributionGridType(EDistrTypes distrType)
Returns grid’s type, which was defined for specified solid distribution distrType
. Possible values are:
GRID_NUMERIC
GRID_SYMBOLIC
GRID_UNDEFINED
std::vector<double> GetNumericGrid(EDistrTypes distrType)
Returns grid of classes for specified solid distribution for Numeric grid. Returns empty vector if this distribution has Symbolic grid.
std::vector<std::string> GetSymbolicGrid(EDistrTypes distrType)
Returns grid of classes for specified solid distribution for Symbolic grid. Retunrs empty vector if this distribution has Numeric grid.
unsigned GetClassesNumber(EDistrTypes distrType)
Returns number of classes for specified solid distribution. Returns 0
if such distribution has not been defined.
std::vector<double> GetClassesMeans(EDistrTypes distrType)
Returns means of classes for specified solid distribution with Numeric grid. Retunrs empty vector if such distribution has not been defined or has Symbolic grid.
std::vector<double> GetPSDGridDiameters()
Returns size grid for particle diameters. Returns empty vector if DISTR_SIZE
distribution has not been defined.
std::vector<double> GetPSDGridVolumes()
Returns size grid for particle volumes. Returns empty vector if DISTR_SIZE
distribution has not been defined.
std::vector<double> GetPSDMeanDiameters()
Returns mean particle diameters. Returns empty vector if DISTR_SIZE
distribution has not been defined.
std::vector<double> GetPSDMeanSurfaces()
Returns mean particle surfaces. Returns empty vector if DISTR_SIZE
distribution has not been defined.
std::vector<double> GetPSDMeanVolumes()
Returns mean particle volumes. Returns empty vector if DISTR_SIZE
distribution has not been defined.
std::vector<double> GetClassesSizes(EDistrTypes distrType)
Returns sizes of classes for specified solid distribution with Numeric grid. Returns empty vector if such distribution has not been defined or has Symbolic grid.
bool IsDistributionDefined(EDistrTypes distrType)
Returns ture
if such solids distribution has been defined in the flowsheet.
void CalculateTM(EDistrTypes distrType, std::vector<double> InDistr, std::vector<double> OutDistr, CTransformMatrix &outTM)
Calculates transformation matrix for onedimensional distribution with type distrType
according to input and output distributions. Obtained matrix can be applied to the stream instead of direct setting of distribution to retain secondary dimensions in multidimensional distribution.
Following algorithm is applied to setup transformation matrix:
Go through the classes of source and target distributions from left to right.
The mostleft notempty class of the initial distribution proceeds to the mostleft notempty class of the output distribution.
Transition to the next class of the initial distribution is performed if the current class is completely transferred to the output distribution.
Transition to the next class of the output distribution is performed if the current class is already full.
Tolerance
double GetAbsTolerance()
Returns absolute tolerance, which has been defined for the flowsheet.
double GetRelTolerance()
Returns relative tolerance, which has been defined for the flowsheet.
Errors and warnings
void RaiseError(std::string Description = "")
Can be called to indicate that an error occurred. Description will be displayed in the simulation’s log and simulation will be stopped after setting.
void RaiseWarning(std::string Description = "")
Can be called to indicate warning. Description will be displayed in the simulation’s log and simulation will not be stopped.
void ShowInfo(std::string Description)
Can be called to write out messages to the simulation’s log screen during the simulation. Description will be displayed in the simulation’s log.
Plots
int AddPlot(std::string PlotName, std::string XAxisName, std::string YAxisName)
Adds new 2dimensional plot with specified name and axis, returns index of the plot. PlotName
must be unique within the unit’s plots. Returns 1
on error.
int AddPlot(std::string PlotName, std::string XAxisName, std::string YAxisName, std::string ZAxisName)
Adds new 3dimensional plot with specified name and axis, returns index of the plot. PlotName
must be unique within the unit’s plots. Returns 1
on error.
int AddCurveOnPlot(std::string PlotName, std::string CurveName)
Adds new curve with specified name CurveName
on the 2dimensional plot with name PlotName
. Returns index of the curve within specified plot. Returns 1
on error.
int AddCurveOnPlot(std::string PlotName, double ZValue)
Adds new curve with specified ZValue
on the 2 or 3dimensional plot with name PlotName
. Returns index of the curve within specified plot.
ZValue
can be time, …
Returns 1
on error.
void AddPointOnCurve(std::string PlotName, std::string CurveName, double X, double Y)
Adds new point on specified curve CurveName
for 2dimensional plot.
void AddPointOnCurve(std::string PlotName, double ZValue, double X, double Y)
Adds new point on specified curve for 3dimensional plot with specified ZValue
. ZValue
can be time, …
void AddPointOnCurve(std::string PlotName, std::string CurveName, std::vector<double> X, std::vector<double> Y)
Adds new points on specified curve CurveName
for 2dimensional plot named PlotName
.
void AddPointOnCurve(std::string PlotName, double ZValue, std::vector<double> X, std::vector<double> Y)
Adds new points on specified curve ZValue
for 3dimensional plot named PlotName
. ZValue
can be time, …
Virtual functions
Virtual function is declared within the base class, which you can redefine in your derived class.
void Simulate(double Time)
Calculates unit on a specified time point Time
for steadystate units. Is called by the suimulator iteratively for all time points for which this unit should be calculated. All logic of the unit’s model must be implemented here.
void Simulate(double StartTime, double EndTime)
Calculates unit for specified time interval for dynamic units. Is called by the suimulator iteratively for all time points for which this unit should be calculated. All logic of the unit’s model must be implemented here.
void Initialize(double Time)
Initializes unit at the time point Time
. Is called by the simulator only once at the start of the simulation. Here some additional objects can be initialized (for example holdups, material streams or state variables).
void SaveState()
Saves current state of the unit. All time dependent variables which weren’t added to the unit be manually saved here with the help of AddStateVariable, AddMaterialStream or AddHoldup.
For flowsheets containing recycled streams, this function is called when the convergence on the current time interval is reached, this also ensures the return to the previous state of the unit if convergence fails during the calculation.
void LoadState()
Loads last saved state of the unit. All time dependent variables which were previously saved in SaveState()
function should be manually loaded here.
void Finalize()
Finalizes unit. Is called by the simulator only once at the end of the simulation. Here closing and cleaning operations can be performed.
Stream
Functions in classes CMaterialStream
and CHoldup
are introduced below.
Basic stream properties
All functions in this section are for both CMaterialStream
and CHoldup
.
std::string GetStreamName()
Returns the name of the material stream / holdup.
void SetStreamName(std::string Name)
Sets the name of the material stream / holdup.
Time points
All functions in this section are for both CMaterialStream
and CHoldup
.
void AddTimePoint(double Time, double SourceTime = 1)
Adds new time point Time
to the material stream / holdup. Data for this time point is copied from SourceTime
. By default (SourceTime = 1
) data will be copied from the previous time point. If this is the first time point in the material stream / holdup, all data will be set to 0
. If such time point already exists, nothing will be done.
void RemoveTimePoint(double Time)
Removes time point Time
from the material stream / holdup, if such point exists.
void RemoveTimePoints(double Start, double End)
Removes all time points from the specified interval, including boundaries.
void RemoveTimePointsAfter(double Start, bool IncludeStart = false)
Removes all data after the specified time point including (if IncludeStart
is set to true
) or excluding (IncludeStart
is set to false
) point Start
.
std::vector<double> GetAllTimePoints()
Returns all time points which are defined in the material stream / holdup.
std::vector<double> GetTimePointsForInterval(double Start, double End, bool ForceInclBoudaries = false)
Returns the list of time points for the specified time interval (incl. boundary points Start
and End
). If ForceInclBoudaries
is set to true
, resulting vector will contain boundary points even if they have not been defined in the material stream / holdup.
double GetLastTimePoint()
Returns last defined time point in the material stream / holdup. Returns 1
if no time points have been defined.
double GetPreviousTimePoint(double Time)
Returns the nearest time point before Time
. Returns 1
if there is no time points before the specified value.
Overall properties
double GetMassFlow(double Time, unsigned Basis)
Specific function for CMaterialStream
.
Returns mass or mole flow of the material stream at the specified time point Time
. If such time point has not been defined, interpolation of data will be done.
Basis
is a basis of results (BASIS_MASS
in [kg/s] or BASIS_MOLL
in [mol/s]):
BASIS_MASS
: \(\dot m\) in [kg/s], total mass flow of the material stream.
BASIS_MOLL
: \(\sum\limits_i \dfrac{\dot m \cdot w_i}{M_i}\) in [mol/s], with \(w_i\) mass fraction of the phase \(i\), and \(M_i\) molar mass of the phase \(i\).
double GetMass(double Time, unsigned Basis)
Specific function for CHoldup
.
Returns mass or mole of the holdup at the specified time point Time
. If such time point has not been defined, interpolation of data will be done.
Basis
is a basis of results (BASIS_MASS
in [kg] or BASIS_MOLL
in [mol]):
BASIS_MASS
: \(m\) in [kg], total mass of the holdup.
BASIS_MOLL
: \(\sum\limits_i \dfrac{m \cdot w_i}{M_i}\) in [mol], with \(w_i\) mass fraction of the phase \(i\), and \(M_i\) molar mass of the phase \(i\).
void SetMassFlow(double Time, double Value, unsigned Basis)
Specific function for CMaterialStream
.
Sets mass flow of the material stream at the time point Time
. Negative values before setting will be converted to 0
. If the time point Time has not been defined in the material stream, then the value will not be set.
Basis
is a basis of results (BASIS_MASS
in [kg/s] or BASIS_MOLL
in [mol/s]):
BASIS_MASS
: in this case you directly have your input Value
as mass flow: \(\dot{m} =\) Value
in [kg/s]
BASIS_MOLL
: in this case you have Value
as mole flow and this should be converted to mass flow. \(\dot{m} =\) Value
\(\cdot \sum\limits_i M_i \cdot w_i\) in [mol/s], with \(w_i\) mass fraction of the phase \(i\), and \(M_i\) molar mass of the phase \(i\).
void SetMass(double Time, double Value, unsigned Basis)
Specific function for CHoldup
.
Sets mass of the holdup at the time point Time
. Previously set negative values will be converted to 0
. If the time point Time
has not been defined in the holdup, then the value will not be set.
Basis
is a basis of results (BASIS_MASS
in [kg] or BASIS_MOLL
in [mol]):
BASIS_MASS
: in this case you directly have your input Value
as mass value: \(m =\) Value
in [kg].
BASIS_MOLL
: in this case you have Value
as mole amount and this should be converted to mass. \(m =\) Value
\(\cdot \sum\limits_i M_i \cdot w_i\) in [mol], with \(w_i\) mass fraction of the phase \(i\), and \(M_i\) molar mass of the phase \(i\).
double GetTemperature(double Time)
Function for both CMaterialStream
and CHoldup
.
Returns temperature of the material stream / holdup at the specified time point Time
in [K]. If such time point has not been defined, interpolation of data will be done.
void SetTemperature(double Time, double Value)
Function for both CMaterialStream
and CHoldup
.
Sets temperature of the material stream / holdup at the time point Time
in [K]. Negative values before setting will be converted to 0
. If time the point Time
has not been defined in the material stream / holdup, the value will not be set.
double GetPressure(double Time)
Function for both CMaterialStream
and CHoldup
.
Returns pressure of the material stream / holdup at the specified time point Time
in [Pa]. If such time point has not been defined, interpolation of data will be done.
void SetPressure(double Time, double Value)
Function for both CMaterialStream
and CHoldup
.
Sets pressure of the material stream / holdup in the time point Time
in [Pa]. Negative values before setting will be converted to 0
. If the time point Time
has not been defined in the material stream / holdup, the value will not be set.
double GetOverallProperty(double Time, unsigned Property, unsigned Basis)
Returns nonconstant physical property value for the overall mixture at the specified time point Time
. If such time point has not been defined, interpolation of data will be done.
Basis
is a basis of results (BASIS_MASS
or BASIS_MOLL
).
Property
is an identifier of a physical property. Available properties are:
FLOW
andTOTAL_FLOW
for material streamCMaterialStream
: refer to function getMassFlow.
MASS
andTOTAL_MASS
for holdupCHoldup
: refer to function getMass.
TEMPERATURE
: refer to function getTemperature.
PRESSURE
: refer to function getPressure.
MOLAR_MASS
: \(\sum\limits_i M_i \cdot w_i\), with \(M\) molar mass of the total flow, \(w_i\) mass fraction of the phase \(i\), and \(M_i\) molar mass of the phase \(i\).
ENTHALPY
:
Set
Basis
asBASIS_MASS
: \(\sum\limits_i H_i \cdot w_i\), with \(H_i\) the enthalpy of the phase \(i\), and \(w_i\) the mass fraction of the phase \(i\).Set
Basis
asBASIS_MOLL
: \(\sum\limits_i H_i \cdot x_i\), with \(H_i\) the enthalpy of the phase \(i\), and \(x_i\) the mole fraction of the phase \(i\).
Note
Definition of overall mixture properties:
Define 
Name 
Units 


Mass / mole flow 
[kg/s] or [mol/s] 

Mass / mole 
[kg] or [mol] 

Temperature 
[K] 

Pressure 
[Pa] 

Molar mass 
[kg/mol] 

Enthalpy 
[J/kg/s] or [J/mol/s] for material stream, [J/kg] or [J/mol] for holdup 
double SetOverallProperty(double Time, unsigned Property, double Value, unsigned Basis)
Sets nonconstant physical property value for the overall mixture at the specified time point Time
.
Basis is a basis of the value (BASIS_MASS
or BASIS_MOLL
).
Property
is an identifier of a physical property. Available properties are:
FLOW
andTOTAL_FLOW
for material streamCMaterialStream
: refer to function setMassFlow.
MASS
andTOTAL_MASS
for holdupCHoldup
: refer to function setMass.
TEMPERATURE
: refer to function setTemperature.
PRESSURE
: refer to function setPressure.
Note
Definition of overall mixture properties:
Define 
Name 
Units 


Mass / mole flow 
[kg/s] or [mol/s] 

Mass / mole 
[kg] or [mol] 

Temperature 
[K] 

Pressure 
[Pa] 

Molar mass 
[kg/mol] 

Enthalpy 
[J/kg/s] or [J/mol/s] 
double CalcTemperatureFromProperty(ECompoundTPProperties Property, double Time, double Value)
Function for both CMaterialStream
and CHoldup
.
Returns temperature of the material stream / holdup for a specific value Value
of the property Property
at the time point Time
. Available properties are those defined in material database.
For further information, please refer to Thermodynamics on this page.
double CalcPressureFromProperty(ECompoundTPProperties Property, double Time, double Value)
Function for both CMaterialStream
and CHoldup
.
Returns pressure of the material stream / holdup for a specific value Value
of the property Property
at the time point Time
. Available properties are those defined in material database.
For further information, please refer to Thermodynamics on this page.
Compounds
double GetCompoundFraction(double Time , std::string CompoundKey, unsigned Basis)
Function for both CMaterialStream
and CHoldup
.
Returns total fraction of the compound with key CompoundKey
at the time point Time
. If such time point has not been defined, interpolation of data will be done.
Basis can be BASIS_MASS
or BASIS_MOLL
.
BASIS_MASS
: \(f_i = \sum \limits_i w_i \cdot f_{i}\), with \(f_i\) the mass fraction of compound \(i\), and \(w_i\) the mass fraction of phase \(i\).
BASIS_MOLL
: \(f_i^{mol} = \sum \limits_i w_i \dfrac{f_i}{M_i \cdot \sum\limits_j \frac{f_{i,j}}{M_j}}\), with \(f_i^{mol}\) the mole fraction of compound \(i\), \(f_{i,j}\) the mass fraction of compound \(j\) in phase \(i\), and \(M_j\) the molar mass of compound \(j\).
double GetCompoundPhaseFraction(double Time, std::string CompoundKey, unsigned Phase, unsigned Basis)
Function for both CMaterialStream
and CHoldup
.
Returns fraction of the compound with the key CompoundKey
in the phase Phase
(SOA_SOLID
, SOA_LIQUID
, SOA_VAPOR
) for the time point Time
. If such time point has not been defined, interpolation of data will be done.
Basis can be BASIS_MASS
or BASIS_MOLL
.
BASIS_MASS
: \(f_{i,j}\), mass fraction of compound \(j\) in phase \(i\).
BASIS_MOLL
: \(f_{i,j}^{mol} = \sum \limits_i w_i \dfrac{f_{i,j}}{M_i \cdot \sum\limits_j \frac{f_{i,j}}{M_j}}\), with \(f_{i,j}^{mol}\) the mole fraction of compound \(j\) in phase \(i\), and \(M_j\) the molar mass of compound \(j\).
void SetCompoundPhaseFraction (double Time, std::string CompoundKey, unsigned Phase, double Fraction, unsigned Basis)
Function for both CMaterialStream
and CHoldup
.
Sets fraction of the compound with key CompoundKey
in phase Phase
(SOA_SOLID
, SOA_LIQUID
, SOA_VAPOR
) for the time point Time
. If such time point has not been defined, nothing will be done. Negative values before setting will be converted to 0
.
Basis can be BASIS_MASS
or BASIS_MOLL
.
BASIS_MASS
: in this case you have your input Fraction
directly as mass fraction of compound \(j\) in phase \(i\): \(f_{i,j} =\) Fraction
.
BASIS_MOLL
: \(f_{i,j} =\) Fraction
\(\cdot \dfrac{M_i}{\sum\limits_j \frac{f_{i,j}}{M_j}}\), with \(f_{i,j}^{mol}\) the mole fraction of compound \(j\) in phase \(i\), and \(M_j\) the molar mass of compound \(j\).
double GetCompoundMassFlow(double Time, std::string CompoundKey, unsigned Phase, unsigned Basis)
Specific function for CMaterialStream
.
Returns mass flow of the compound with key CompoundKey in phase Phase (SOA_SOLID
, SOA_LIQUID
, SOA_VAPOR
) for the time point Time
. If such time point has not been defined, interpolation of data will be done.
Basis is a basis of value (BASIS_MASS
in [kg/s] or BASIS_MOLL
in [mol/s]).
BASIS_MASS
: \(\dot{m}_{i,j} = w_i \cdot f_{i,j} \cdot \dot{m}\), with \(\dot m_{i,j}\) the mass flow of compound \(j\) in phase \(i\), \(w_i\) the mass fraction of phase \(i\), and \(f_{i,j}\) the mass fraction of compound \(j\) in phase \(i\).
BASIS_MOLL
: \(\dot{m}_{i,j} = w_i \cdot f_{i,j} \cdot \sum\limits_k \dfrac{\dot{m} \cdot w_k}{M_k}\), with \(\dot m\) the total mass flow of the material stream, and \(M_k\) the molar mass of phase \(k\).
double GetCompoundMass(double Time, std::string CompoundKey, unsigned Phase, unsigned Basis)
Specific function for CHoldup
.
Returns mass of the compound with key CompoundKey in phase Phase (SOA_SOLID
, SOA_LIQUID
, SOA_VAPOR
) for the time point Time
. If such time point has not been defined, interpolation of data will be done.
Basis is a basis of value (BASIS_MASS
in [kg] or BASIS_MOLL
in [mol]).
BASIS_MASS
: \(m_{i,j} = w_i \cdot f_{i,j} \cdot m\), with \(m_{i,j}\) the mass of compound \(j\) in phase \(i\), \(w_i\) the mass fraction of phase \(i\), and \(f_{i,j}\) the mass fraction of compound \(j\) in phase \(i\).
BASIS_MOLL
: \(m_{i,j} = w_i \cdot f_{i,j} \cdot \sum\limits_k \dfrac{m \cdot w_k}{M_k}\), with \(m\) the total mass of the holdup, and \(M_k\) the molar mass of phase \(k\).
double GetCompoundConstant(std::string CompoundKey, ECompoundConstProperties ConstProperty)
Function for both CMaterialStream
and CHoldup
.
Returns value of the constant physical property ConstProperty
for the specified compound. These properties are stored in material database. Available constants are:
SOA_AT_NORMAL_CONDITIONS
NORMAL_BOILING_POINT
NORMAL_FREEZING_POINT
CRITICAL_TEMPERATURE
CRITICAL_PRESSURE
MOLAR_MASS
STANDARD_FORMATION_ENTHALPY
HEAT_OF_FUSION_AT_NORMAL_FREEZING_POINT
HEAT_OF_VAPORIZATION_AT_NORMAL_BOILING_POINT
REACTIVITY_TYPE
CONST_PROP_USER_DEFINED_XX
Note
Definition of constant properties for pure compounds:
Define 
Name 
Unit 


State of aggregation at normal conditions 
[] 

Normal boiling point 
[K] 

Normal freezing point 
[K] 

Critical temperature 
[K] 

Critical pressure 
[Pa] 

Molar mass 
[kg/mol] 

Standard formation enthalpy 
[J/mol] 

Heat of fusion at normal freezing point 
[J/mol] 

Heat of vaporization at normal boiling point 
[J/mol] 

Reactivity type 
[] 

User defined property 
[] 
double GetCompoundTPDProp(std::string CompoundKey, unsigned Property, double Temperature, double Pressure)
Function for both CMaterialStream
and CHoldup
.
Returns value of the temperature / pressuredependent physical Property (which are stored in the database of materials) for the compound with the specified Temperature
in [K] and Pressure
in [Pa]. Available properties are:
DENSITY
HEAT_CAPACITY_CP
VAPOR_PRESSURE
VISCOSITY
THERMAL_CONDUCTIVITY
PERMITTIVITY
ENTHALPY
TP_PROP_USER_DEFINED_XX
Note
Definition of temperaturedependent compound properties:
Define 
Name 
Unit 


Density 
[kg/m \(^3\)] 

Heat capacity \(C_p\) 
[J/(kg·K)] 

Vapor pressure 
[Pa] 

Viscosity 
[Pa·s] 

Thermal conductivity 
[W/(m·K)] 

Permittivity 
[F/m] 

Enthalpy 
[J/kg/s] or [J/mol/s] for material stream, [J/kg] or [J/mol] for holdup 

User defined property 
[] 
double GetCompoundTPDProp(double Time, std::string CompoundKey, unsigned Property)
Function for both CMaterialStream
and CHoldup
.
Returns value of the temperature / pressuredependent physical Property
(which are stored in the database of materials) for the compound with the current temperature and pressure. Available properties are:
DENSITY
HEAT_CAPACITY_CP
VAPOR_PRESSURE
VISCOSITY
THERMAL_CONDUCTIVITY
PERMITTIVITY
ENTHALPY
TP_PROP_USER_DEFINED_XX
Note
Definition of temperaturedependent compound properties:
Define 
Name 
Unit 


Density 
[kg/m \(^3\)] 

Heat capacity \(C_p\) 
[J/(kg·K)] 

Vapor pressure 
[Pa] 

Viscosity 
[Pa·s] 

Thermal conductivity 
[W/(m·K)] 

Permittivity 
[F/m] 

Enthalpy 
[J/kg/s] or [J/mol/s] for material stream, [J/kg] or [J/mol] for holdup 

User defined property 
[] 
double GetCompoundInteractionProp(std::string CompoundKey1, std::string CompoundKey2, unsigned Property, double Temperature, double Pressure)
Function for both CMaterialStream
and CHoldup
.
Returns the value of the interaction property Property
for the selected compounds under the specified Temperature
in [K] and Pressure
in [Pa]. These properties are stored in the material database. Available properties are:
INTERFACE_TENSION
INT_PROP_USER_DEFINED_XX
Note
Definition of interaction properties between two pure compounds:
Define 
Name 
Unit 


Interface tension 
[N/m] 

User defined property 
[] 
double GetCompoundInteractionProp(double Time, std::string CompoundKey1, std::string CompoundKey2, unsigned Property)
Function for both CMaterialStream
and CHoldup
.
Returns the value of the interaction property Property
for the selected compounds under the current temperature and pressure. These properties are stored in the material database. Available properties are:
INTERFACE_TENSION
INT_PROP_USER_DEFINED_XX
Note
Definition of interaction properties between two pure compounds:
Define 
Name 
Unit 


Interface tension 
[N/m] 

User defined property 
[] 
Phases
double GetPhaseMassFlow(double Time, unsigned Phase, unsigned Basis = BASIS_MASS)
Specific function for CMaterialStream
.
Returns mass flow of the specified phase Phase
(SOA_SOLID
, SOA_LIQUID
, SOA_VAPOR
) in the material stream for the time point Time
. If such time point has not been defined, the value will be interpolated.
Basis is a basis of value (BASIS_MASS
in [kg/s] or BASIS_MOLL
in [mol/s]).
BASIS_MASS
: \(\dot{m}_i = \dot{m} \cdot w_i\), with \(\dot{m}_i\) the mass flow of phase \(i\), \(w_i\) the mass fraction of phase \(i\), and \(\dot{m}\) the total mass flow of the material stream.
BASIS_MOLL
: \(\dot{n}_i = \dfrac{\dot{m} \cdot w_i}{M_i}\), with \(\dot{n}_i\) the mole flow of phase \(i\), \(w_i\) the mass fraction of phase \(i\), \(\dot{m}\) the total mass flow of the material stream, and \(M_i\) the molar mass of phase \(i\).
double GetPhaseMass(double Time, unsigned Phase, unsigned Basis = BASIS_MASS)
Specific function for CHoldup
.
Returns mass of the specified phase Phase
(SOA_SOLID
, SOA_LIQUID
, SOA_VAPOR
) in the holdup for the time point Time
. If such time point has not been defined, the value will be interpolated.
Basis is a basis of value (BASIS_MASS
in [kg] or BASIS_MOLL
in [mol]).
BASIS_MASS
: \(m_i = m \cdot w_i\), with \(m_i\) the mass of phase \(i\), \(w_i\) the mass fraction of phase \(i\), and \(m\) the total mass flow of the material stream.
BASIS_MOLL
: \(n_i = \dfrac{m \cdot w_i}{M_i}\), with \(n_i\) the mole of phase \(i\), \(w_i\) the mass fraction of phase \(i\), \(m\) the total mass flow of the material stream, and \(M_i\) the molar mass of phase \(i\).
void SetPhaseMassFlow(double Time, unsigned Phase, double Value, unsigned Basis)
Specific function for CMaterialStream
.
Sets mass flow of the specified phase Phase
(SOA_SOLID
, SOA_LIQUID
, SOA_VAPOR
) in the material stream for the time point Time
.
Is performed by calculation and setting of a new total mass flow of the material stream and new phase fractions (according to the new mass flow of the specified phase). Negative values before setting will be converted to 0
. If there is no specified time point or phase in the material stream, the value will not be set.
Basis is a basis of value (BASIS_MASS
in [kg/s] or BASIS_MOLL
in [mol/s]).
BASIS_MASS
: in this case you have your input Value
as mass flow of one defined phase: \(\dot m_i =\) Value
and \(w_i = \dot m_i / \dot m\). Meanwhile, the total mass flow \(\dot m\) changes due to assignment for \(\dot m_i\): \(\dot m = \dot m_{old} + (\) Value
\( \dot{m}_{i,old})\). Here \(\dot{m}_i\) stands for the mass flow of phase \(i\), \(w_i\) for the mass fraction of phase \(i\), and \(\dot{m}\) for the total mass flow of the material stream.
BASIS_MOLL
: in this case you have your input Value
as mole flow of one defined phase: \(\dot m_i =\) Value
\(\cdot M_i ` and :math:`w_i = \dot m_i / \dot m\). Meanwhile, the total mass flow \(\dot m\) changes due to assignment for \(\dot m_i\): \(\dot m = \dot m_{old} + (\) Value
\(\cdot M_i  \dot{m}_{i,old})\). Here \(w_i\) stands for the mass fraction of phase \(i\), \(\dot{m}\) for the total mass flow of the material stream, and \(M_i\) for the molar mass of phase \(i\).
void SetPhaseMass(double Time, unsigned Phase, double Value, unsigned Basis)
Specific function for CHoldup
.
Sets mass of the specified phase Phase
(SOA_SOLID
, SOA_LIQUID
, SOA_VAPOR
) in the holdup for the time point Time
.
Is performed by calculation and setting of a new total mass flow of the holdup and new phase fractions (according to the new mass of the specified phase). Negative values before setting will be converted to 0
. If there is no specified time point or phase in the holdup, the value will not be set.
Basis is a basis of value (BASIS_MASS
in [kg] or BASIS_MOLL
in [mol]).
BASIS_MASS
: in this case you have your input Value
as the mass of one defined phase: \(m_i =\) Value
and \(w_i = m_i / m\). Meanwhile, the total mass \(m\) changes due to assignment for \(m_i\): \(\dot m = \dot m_{old} + (\) Value
\( \dot{m}_{i,old})\). Here \(m_i\) stands for the mass of phase \(i\), \(w_i\) for the mass fraction of phase \(i\), and \(m\) for the total mass of the holdup.
BASIS_MOLL
: in this case you have your input Value
as mole flow of one defined phase: \(m_i = ` ``Value`\) \(\cdot M_i\) , \(w_i = m_i / m\). Meanwhile, the total mass \(m\) changes due to assignment for \(m_i\): \(m = m_{old} + (\) Value
\(\cdot M_i  m_{i,old})\). Here \(m_i\) stands for the mass of phase \(i\), \(w_i\) for the mass fraction of phase \(i\), \(m\) for the total mass of the holdup, and \(M_i\) for the molar mass of phase \(i\).
double GetSinglePhaseProp(double Time, unsigned Property, unsigned Phase, unsigned Basis)
Function for both CMaterialStream
and CHoldup
.
Returns nonconstant physical property value for the phase mixture Phase
(SOA_SOLID
, SOA_LIQUID
, SOA_VAPOR
) for the specified time point. If such time point has not been defined, interpolation of data will be done.
Basis is a basis of results (BASIS_MASS
or BASIS_MOLL
).
Property
is an identifier of a physical property. Available properties are:
FLOW
: only for classCMaterialStream
. Refer to function getMassFlow.
MASS
: only for classCHoldup
. Refer to function getMass.
TEMPERATURE
: refer to function getTemperature.
PRESSURE
: refer to function getPressure.
PHASE_FRACTION
,FRACTION
:
BASIS_MASS
: function returns \(w_i\), mass fraction of phase \(i\).
BASIS_MOLL
: function returns result of \(\left ( \dfrac{w_i}{M_i \cdot \sum\limits_j \frac{w_j}{M_j}} \right )\), with \(w_i\) the mass fraction of phase \(i\), and \(M_i\) the molar mass of phase \(i\).
MOLAR_MASS
: calculate the molar mass of the phase \(M\) by \(\left ( \frac{1}{M} = \sum\limits_i \frac{w_i}{M_i} \right )\), with \(M_i\) the molar mass of phase \(i\), and \(w_i\) the mass fraction of phase \(i\).
DENSITY
: refer to function getPhaseTPDProp.
HEAT_CAPACITY_CP
: refer to function getPhaseTPDProp.
THERMAL_CONDUCTIVITY
: refer to function getPhaseTPDProp.
VISCOSITY
: refer to function getPhaseTPDProp.
VAPOR_PRESSURE
: refer to function getPhaseTPDProp.
PERMITTIVITY
: refer to function getPhaseTPDProp.
TP_PROP_USER_DEFINED_XX
: refer to function getPhaseTPDProp.
ENTHALPY
:For solid and liquid phase: \(h = h_0 + C_p \cdot \Delta T + \frac{M}{\rho} (P  P_0)\)
BASIS_MASS
: \(H = \sum\limits_i \frac{h_i \cdot f_i}{M_i}\)
BASIS_MOLL
: \(H = \sum\limits_i h_i \cdot f_i\)For vapor phase: \(h = h_0 + C_p \cdot \Delta T\)
BASIS_MASS
: \(H = \sum\limits_i \frac{h_i \cdot f_i}{M_i}\)
BASIS_MOLL
: \(H = \sum\limits_i h_i \cdot f_i\)Note
Notations for enthalpy:
\(H\) – enthalpy of the phase. [J/kg/s] or [J/mol/s] for material stream; [J/kg] or [J/mol] for holdup
\(h_i\) – enthalpy of the compound \(i\) [J/mol]
\(f_i\) – mass fraction of the compound \(i\) in phase
\(M_i\) – molar mass of the compound \(i\)
\(h_0\) – formation enthalpy [J/mol]
\(C_p\) – heat capacity for constant pressure of the compound
\(\Delta T\) – difference between the temperature at normal conditions (298.15 K) and current temperature
\(P\) – current pressure
\(P_0\) – pressure at normal conditions (101325 Pa)
Note
Definition of singlephase mixture properties:
Define 
Name 
Unit 


Mass flow 
[kg/s] or [mol/s] 

Mass 
[kg] or [mol] 

Temperature 
[K] 

Pressure 
[Pa] 

Phase fraction 
[] 

Molar mass 
[kg/mol] 

Density 
[kg/m \(^3\)] 

Heat capacity \(C_p\) 
[J/(kg·K)] 

Thermal conductivity 
[W/(m·K)] 

Viscosity 
[Pa·s] 

Vapor pressure 
[Pa] 

Enthalpy 
[J/kg/s]or [J/mol/s] for material stream, [J/kg] or [J/mol] for holdup 

Permittivity 
[F/m] 

User defined property 
[] 
void SetSinglePhaseProp(double Time, unsigned Property, unsigned Phase, double Value, unsigned Basis)
Function for both CMaterialStream
and CHoldup
.
Sets nonconstant physical property value for phase mixture Phase
(SOA_SOLID
, SOA_LIQUID
, SOA_VAPOR
) for the specified time point Time
. If there is no specified time point or phase in the material stream or holdup, the value will not be set.
Property
is an identifier of a physical property. Available properties are:
FLOW
: only for classCMaterialStream
. Refer to function setPhaseMassFlow.
MASS
: only for classCHoldup
. Refer to function setPhaseMass.
FRACTION
: mass fraction of thePhase
is set toValue
.
Basis is a basis of value (BASIS_MASS
or BASIS_MOLL
).
double GetPhaseTPDProp(double Time, unsigned Property, unsigned Phase)
Function for both CMaterialStream
and CHoldup
.
Returns value of temperature / pressuredependent physical property for specified phase (SOA_SOLID
, SOA_LIQUID
, SOA_VAPOR
) for the time point Time
. If such time point has not been defined, interpolation of data will be done.
Available properties are:
DENSITY
:
For solid phase: is calculated by \(\rho = \sum\limits_{i,j} \rho_i \, (1  \varepsilon_j)\,f_{i,j}\), with \(\varepsilon_j\) the porosity in interval \(j\), and \(f_{i,j}\) the mass fraction of compound \(i\) with porosity \(j\).
For liquid and vapor phase: is calculated by \(\frac{1}{\rho} = \sum\limits_i \frac{w_i}{\rho_i}\), with \(w_i\) the mass fraction of compound \(i\) in
Phase
.
HEAT_CAPACITY_CP
: is calculated by \(C_p = \sum\limits_i w_i \cdot C_{p,i}\), with \(C_{p,i}\) the heat capacity of compound \(i\), and \(w_i\) the mass fraction of compound \(i\) inPhase
.
VAPOR_PRESSURE
: is calculated by \(P_v = \min\limits_{i} (P_v)_i\), with \((P_v)_i\) vapor pressure of compound \(i\).
VISCOSITY
:
For solid phase: is calculated by \(\eta = \sum\limits_i w_i\, \eta_i\), with \(\eta_i\) the viscosity of compound \(i\), and \(w_i\) the mass fraction of compound \(i\).
For liquid phase: is calculated by \(\ln \eta = \dfrac{\sum\limits_i w_i\,\ln \eta_i}{\sum\limits_i x_i\,\sqrt{M_i}}\), with \(\eta_i\) the viscosity of compound \(i\), \(w_i\) the mass fraction of compound \(i\) in Phase and \(x_i\) the mole fraction of compound \(i\) in
Phase
.For vapor phase: \(\eta = \dfrac{\sum\limits_i x_i\,\sqrt{M_i}\,\eta_i}{\sum\limits_i x_i\,\sqrt{M_i}}\), with \(\eta_i\) the viscosity of compound \(i\), \(w_i\) the mass fraction of compound \(i\) in Phase, and \(x_i\) the mole fraction of compound \(i\) in
Phase
.
THERMAL_CONDUCTIVITY
:
For solid phase: is calculated by \(\lambda = \sum\limits_i w_i \, \lambda_i\), with \(\lambda_i\) the thermal conductivity of compound \(i\).
For liquid phase: is calculated by \(\lambda = \dfrac{1}{\sqrt{\sum\limits_i x_i \, \lambda_i^{2}}}\), with \(\lambda_i\) the thermal conductivity of compound \(i\).
For vapor phase: is calculated by \(\lambda = \sum\limits_i \dfrac{x_i\,\lambda_i}{\sum\limits_j x_j\, F_{i,j}}\), \(F_{i,j} = \frac{(1 + \sqrt{\lambda_i^4 / \lambda_j} \sqrt{M_j / M_i})^2}{\sqrt{8(1 + M_i / M_j)}}\). With \(M_i\) the molar mass of compound \(i\).
PERMITTIVITY
: is calculated by \(\varepsilon = \sum\limits_i w_i\,\varepsilon_i\), with \(\varepsilon_i\) the permittivity of compound \(i\), and \(w_i\) the mass fraction of compound \(i\) inPhase
.
ENTHALPY
: is calculated by \(H = \sum\limits_i w_i\,H_i\), with \(H_i\) the enthalpy of compound \(i\), and \(w_i\) the mass fraction of compound \(i\) inPhase
.
TP_PROP_USER_DEFINED_XX
: is calculated by \(Y = \sum\limits_i w_i\,Y_i\), with \(Y_i\) the property value of compound \(i\), and \(w_i\) the mass fraction of compound \(i\) inPhase
.
Note
Definition of temperaturedependent compound properties:
Define 
Name 
Unit 


Density 
[kg/m \(^3\)] 

Heat capacity \(C_p\) 
[J/(kg·K)] 

Vapor pressure 
[Pa] 

Viscosity 
[Pa·s] 

Thermal conductivity 
[W/(m·K)] 

Permittivity 
[F/m] 

Enthalpy 
[J/kg/s] or [J/mol/s] for material stream, [J/kg] or [J/mol] for holdup 

User defined property 
[] 
Solid distributed properties
All functions in this section are for both CMaterialStream
and CHoldup
.
double GetFraction(double Time, std::vector<unsigned> Coords)
Returns solid mass fraction by specified coordinates according to all defined distributions. If such time point has not been defined, interpolation of data will be done.
void SetFraction(double Time, std::vector<unsigned> Coords, double Value)
Sets solid mass fraction by specified coordinates according to all defined distributions. If such time point has not been defined in the material stream / holdup, nothing will be done.
Direct setting of fractions to the material stream / holdup leads to a change of all dependent distributions. Approach with transformation matrix should be used to avoid this.
bool GetDistribution(double Time, EDistrTypes Dim, std::vector<double>& Result)
Returns vector of distributed property for specified time point Time
and dimension Dim
. If such time point has not been defined in the material stream / holdup, then linear interpolation will be used to obtain data.
Returns false
on error.
bool GetDistribution(double Time, EDistrTypes Dim1, EDistrTypes Dim2, CDense2DMatrix& Result)
Returns matrix of two distributed dependent properties Dim1
and Dim2
for the specified time point Time
.
If such time point has not been defined in the material stream / holdup, then linear interpolation will be used to obtain data. Rows of resulting matrix will correspond to Dim1
, columns – to Dim2
.
Returns false
on error.
bool GetDistribution(double Time, std::vector<EDistrTypes> Dims, CDenseMDMatrix& Result)
Returns multidimensional matrix of distributed dependent properties for specified time point Time
and dimensions Dims
. If such time point has not been defined in the material stream / holdup, then linear interpolation will be used to obtain data.
Returns false
on error.
bool GetDistribution(double Time, EDistrTypes Dim, std::string Compound, std::vector<double>& Result)
Returns vector of distributed property for specified time point Time
, dimension Dim
and compound Compound
.
Input dimensions should not include distribution by compounds (DISTR_COMPOUNDS
). If specified compound has not been defined in the material stream / holdup, nothing will be done. If specified time point has not been defined, then linear interpolation will be used to obtain data.
Returns false
on error.
bool GetDistribution(double Time, EDistrTypes Dim1, EDistrTypes Dim2, std::string Compound, CDense2DMatrix& 2DResult)
Returns matrix of two distributed dependent properties Dim1
and Dim2
for specified compound Compound
and time point Time
.
Input dimensions should not include distribution by compounds (DISTR_COMPOUNDS
). If specified compound has not been defined in the material stream / holdup, nothing will be done. If specified time point has not been defined, then linear interpolation will be used to obtain data. Rows of resulting matrix will correspond to Dim1
, columns to Dim2
.
Returns false
on error.
bool GetDistribution(double Time, std::vector<EDistrTypes> Dims, std::string Compound, CDenseMDMatrix& MDResult)
Returns multidimensional matrix of distributed dependent properties for specified time point Time
, dimensions Dims
and compound Compound
.
Input dimensions should not include distribution by compounds (DISTR_COMPOUNDS
). If specified compound has not been defined in the material stream / holdup, nothing will be done. If specified time point has not been defined, then linear interpolation will be used to obtain data.
Returns false
on error.
bool SetDistribution(double Time, EDistrTypes Dim, std::vector<double> Distr)
Sets distributed property Distr
of type Dim
for specified time point Time
.
If such time point or dimension doesn’t exist, nothing will be done. Returns false
on error.
Direct setting of distribution to the material stream / holdup leads to a change of all dependent distributions. Approach with transformation matrix should be used to avoid this.
bool SetDistribution(double Time, EDistrTypes Dim1, EDistrTypes Dim2, CDense2DMatrixDistr)
Sets matrix Distr
of two dependent distributed properties of types Dim1
and Dim2
for specified time point Time
. If such time point or dimensions don’t exist nothing will be done.
Returns false
on error.
Direct setting of distribution to the material stream / holdup leads to a change of all dependent distributions. Approach with transformation matrix should be used to avoid this.
bool SetDistribution(double Time, CDenseMDMatrix Distr)
Sets multidimensional matrix Distr
of dependent distributed properties for specified time point Time
. If such time point or dimensions, which are specified in Distr
, don’t exist, nothing will be done.
Returns false
on error.
Direct setting of distribution to the material stream / holdup leads to a change of all dependent distributions. Approach with transformation matrix should be used to avoid this.
bool SetDistribution(double Time, EDistrTypes Dim, std::string Compound, std::vector<double> Distr)
Sets distributed property Distr
of type Dim
for specified compound Compound
and time point Time
. If such time point, compound or dimension doesn’t exist, nothing will be done. Input dimensions should not include distribution by compounds (DISTR_COMPOUNDS
).
Returns false
on error.
Direct setting of distribution to the holdup leads to a change of all dependent distributions. Approach with transformation matrix should be used to avoid this.
bool SetDistribution(double Time, EDistrTypes Dim1, EDistrTypes Dim2, std::string Compound, CDense2DMatrix 2DDistr)
Sets matrix 2DDistr
of two dependent distributed properties of types Dim1
and Dim2
for specified compound Compound
and time point Time
. If such time point, compound or dimensions don’t exist, nothing will be done. Input dimensions should not include distribution by compounds (DISTR_COMPOUNDS
).
Returns false
on error.
Direct setting of distribution to the holdup leads to a change of all dependent distributions. Approach with transformation matrix should be used to avoid this.
bool SetDistribution(double Time, std::string Compound, CDenseMDMatrix MDDistr)
Sets multidimensional matrix MDDistr
of dependent distributed properties for specified compound Compound
and time point Time
. If such time point, compound or dimensions, which are specified in MDDistr
, don’t exist, nothing will be done. Input dimensions should not include distribution by compounds (DISTR_COMPOUNDS
).
Returns false
on error.
Direct setting of distribution to the holdup leads to a change of all dependent distributions. Approach with transformation matrix should be used to avoid this.
bool ApplyTM(double Time, CTransformMatrix Transformation)
Transforms matrix of distributed parameters of solids for time point Time
by applying a movement matrix Transformation
. Returns true
if the transformation was successful.
bool ApplyTM (double Time, std::string Compound, CTransformMatrix Transformation)
Transforms matrix of distributed parameters of solids for specified compound Compound
and time point Time
by applying a movement matrix Transformation
. Dimensions of transformation matrix should not include distribution by compounds (DISTR_COMPOUNDS
). Returns true
if the transformation was successful.
void NormalizeDistribution(double Time)
Normalizes data in solid distribution matrix at the specified time point Time
. If Time
has not been defined, nothing will be done.
void NormalizeDistribution(double Start, double End)
Normalizes data in solid distribution matrix in each time point from interval [Start
; End
].
void NormalizeDistribution()
Normalizes data in solid distribution matrix in all defined time points.
Praticle size distribution
All functions in this section are for both CMaterialStream
and CHoldup
.
std::vector<double> GetPSD(double Time, EPSDType PSDType, EPSDGridType PSDGridType)
Returns particle size distribution of the total mixture of the solid phase at the time point Time
.
PSDGridType
defines grid units if needed: EPSDGridType::DIAMETER
for diameter in [m]; or EPSDGridType::VOLUME
for volume in [m \(^3\)].
PSDType
is a type of distribution. Available types are: PSD_q0
, PSD_Q0
, PSD_q2
, PSD_Q2
, PSD_q3
, PSD_Q3
, PSD_MassFrac
, PSD_Number
.
PSD data is originally stored in a form of mass fractions and all transformations are performed by the following equations:
PSD_q0
: numberrelated distribution of particles: \(q_{0,i} = \dfrac{N_i}{N_{tot} \cdot \Delta d_i}\);PSD_Q0
: \(Q_{0,i} = Q_{0,i1} + q_{0,i} \cdot \Delta d_i\)PSD_q2
: surfacearearelated distribution of particles: \(q_{2,i} = \dfrac{Q_{2,i}  Q_{2,i1}}{\Delta d_i}\)PSD_Q2
: \(Q_{3,i} = \dfrac{\sum\limits_{j=0}^i N_j\,\pi\,d_j^2}{\sum\limits_j N_j \,\pi d_j^2}\)PSD_q3
: \(q_{3,i} = w_i / \Delta d_i\)PSD_Q3
: \(Q_{3,0} = w_0\), \(Q_{3,i} = Q_{3,i1} + w_i\)PSD_MassFrac
: returns the size distribution in the form of mass fractions with the total sum of 1.PSD_Number
: obtains numberrelated distribution of particles depends on several conditions. Three cases of calculation can be distinguished:If only one compound is specified: \(N_i = \dfrac{m_i}{\rho \, \frac{\pi}{6}\, d_i^3}\).
For several compounds: \(N_i = \sum\limits_j \dfrac{M_{tot} \cdot w_{i,j}}{\frac{\pi \cdot d_i^3}{6} \cdot \rho_j}\).
3. If distribution by particle porosity has been defined: \(N_i = \sum\limits_j N_{i,j}\), with \(N_{i,j} = \sum\limits_k \dfrac{M_{tot} \cdot w_{i,j,k}}{\frac{\pi \cdot d_i^3}{6} \cdot \rho_j \cdot (1  \varepsilon_k)}\).
Note
Notations:
\(i\) – index of size classes
\(j\) – index of compounds
\(k\) – index of porosities
\(d_i\) – particle diameter of class \(i\)
\(\Delta d_i\) – size of the class \(i\)
\(m_i\) – mass of particles of class \(i\)
\(M_{tot}\) – total mass of particles
\(N_i\) – number of particles of class \(i\)
\(N_{i,j}\) – number of particles of compound \(j\) with size class \(i\)
\(N_{tot}\) – total number of particles
\(w_i\) – mass fraction of particles of class \(i\)
\(w_{i,j}\) – mass fraction of particles of compound \(j\) with size class \(i\)
\(w_{i,j,k}\) – mass fraction of particles of compound \(j\) with size class \(i\) and porosity \(k\)
\(\rho_j\) – density of compound \(j\)
\(\varepsilon_k\) – porosity of class \(k\)
\(q_0\) – numberrelated density distribution
\(Q_0\) – numberrelated cumulative distribution
\(q_2\) – surfacearearelated density distribution
\(Q_2\) – surfacearearelated cumulative distribution
\(q_3\) – massrelated density distribution
\(Q_3\) – massrelated cumulative distribution
std::vector<double> GetPSD(double Time, EPSDType PSDType, std::string Compound, EPSDGridType PSDGridType)
Returns particle size distribution of compound Compound
of the solid phase of the material stream / holdup at the time point Time
.
PSDType
is a type of distribution, please refer to function getPSD(Time, PSDType, PSDGridType) for detailed information.
PSDGridType
defines grid units if needed: EPSDGridType::DIAMETER
for diameter in [m]; or EPSDGridType::VOLUME
for volume in [m \(^3\)].
void SetPSD(double Time, EPSDType PSDType, std::vector<double> PSD, EPSDGridType PSDGridType)
Sets particle size distribution with type ··PSDType·· to the solid phase of the material stream / holdup for time point Time
.
Direct setting of PSD to the material stream / holdup leads to a change of all dependent distributions. Approach with transformation matrix should be used to avoid this.
Available PSDType
are: PSD_q0
, PSD_Q0
, PSD_q2
, PSD_Q2
, PSD_q3
, PSD_Q3
, PSD_MassFrac
, PSD_Number
.
Please note that when using PSD_Number
, if the total mass of the particles given in the number distribution differs from the material stream / holdup’s total particle mass, the number will be normalized to match the material stream / holdup’s particle mass.
As mass fractions are used to store data, PSD will be converted using functions Convertq0ToMassFractions()
, ConvertQ0ToMassFractions()
, Convertq3ToMassFractions()
, ConvertQ3ToMassFractions()
, ConvertNumbersToMassFractions()
, please refer to section Particle size distribution for detailed explanation about those functions.
PSDGridType
defines grid units if needed: EPSDGridType::DIAMETER
for diameter in [m]; or EPSDGridType::VOLUME
for volume in [m \(^3\)].
void SetPSD(double Time, EPSDType PSDType, std::string Compound, std::vector<double> PSD, EPSDGridType PSDGridType)
Sets PSD with type PSDType
for the specific compound Compound
to the solid phase of the material stream / holdup for time point Time.
Direct setting of PSD to the material stream / holdup leads to a change of all dependent distributions. Approach with transformation matrix should be used to avoid this.
Available PSDType
are: PSD_q0
, PSD_Q0
, PSD_q2
, PSD_Q2
, PSD_q3
, PSD_Q3
, PSD_MassFrac
, PSD_Number
.
Please note that when using PSD_Number, If the total mass of the particles given in the number distribution differs from the material stream / holdup’s total particle mass, the number will be normalized to match the material stream/holdup’s particle mass.
As mass fractions are used to store data, PSD will be converted using functions Convertq0ToMassFractions()
, ConvertQ0ToMassFractions()
, Convertq3ToMassFractions()
, ConvertQ3ToMassFractions()
, ConvertNumbersToMassFractions()
, please refer to section Particle size distribution for detailed explanation about those functions.
PSDGridType
defines grid units if needed: EPSDGridType::DIAMETER
for diameter in [m]; or EPSDGridType::VOLUME
for volume in [m \(^3\)].
Lookup tables
All functions in this section are for both CMaterialStream
and CHoldup
.
CLookupTable* GetLookupTable(ECompoundTPProperties Property, EDependencyTypes DependencyType, double Time)
Creates (if not yet exists), fills with compounds fractions and returns a corresponding lookup table for the specified Property
(see table at the end of this section), DependencyType
(DEPENDENCE_TEMP
or DEPENDENCE_PRES
) and Time
.
For more information, please refer to section Thermodynamics.
double CalcTemperatureFromProperty(ECompoundTPProperties Property, double Time, double Value)
Reads the temperature from the corresponding lookup table for a specific Value
of the selected Property
(see table at the end of this section) at the corresponding time point Time
.
For more information, please refer to section Thermodynamics.
double CalcPressureFromProperty(ECompoundTPProperties Property, double Time, double Value)
Reads the pressure from the corresponding lookup table for a specific Value
of the selected Property
(see table at the end of this section) at the corresponding time point Time
.
For more information, please refer to section Thermodynamics.
double CalcPropertyFromTemperature(ECompoundTPProperties Property, double Time, double T)
Reads the value of the specified Property
(see table at the end of this section) at Time
from the corresponding lookup table for the given temperature T
.
For more information, please refer to section Thermodynamics.
double CalcPropertyFromPressure(ECompoundTPProperties Property, double Time, double P)
Reads the value of the specified Property
(see table at the end of this section) at Time
from the corresponding lookup table for the given pressure P
.
For more information, please refer to section Thermodynamics.
Note
Definition of temperaturedependent compound properties:
Define 
Name 
Unit 


Density 
[kg/m \(^3\)] 

Heat capacity \(C_p\) 
[J/(kg·K)] 

Vapor pressure 
[Pa] 

Viscosity 
[Pa·s] 

Thermal conductivity 
[W/(m·K)] 

Permittivity 
[F/m] 

Enthalpy 
[J/kg/s] or [J/mol/s] for material stream, [J/kg] or [J/mol] for holdup 

User defined property 
[] 
Other streams
Following functions are for class CMaterialStream
.
void CopyFromStream(CMaterialStream *SrcStream, double Time, bool DeleteDataAfter = true)
Copies all stream data from SrcStream
for specified time point Time
to the current material stream.
If flag DeleteDataAfter
is set to true
, all data after the time point Time
in the destination stream will be removed before being copied.
void CopyFromStream(CMaterialStream *SrcStream, double Start, double End, bool DeleteDataAfter = true)
Copies all stream data from SrcStream
on the certain time interval to the current material stream. Boundary points Start
and End
are included into this interval.
If flag DeleteDataAfter
is set to true
, all data after the time point Time
in the destination stream will be removed before being copied.
void CopyFromStream(double TimeDst, CMaterialStream *SrcStream, double TimeSrc, bool DeleteDataAfter = true)
Copies all stream data from time point TimeSrc
of material stream SrcStream
to the time point TimeDst
of this material stream.
If flag DeleteDataAfter
is set to true
, all data after the time point Time
in the destination stream will be removed before being copied.
void CopyFromHoldup(CHoldup *SrcHoldup, double Time, double MassFlow, bool DeleteDataAfter = true)
Copies all data from SrcHoldup
for the specified time point to the current material stream and sets mass flow MassFlow
.
If flag DeleteDataAfter
is set to true
, all data after the time point Time
in the destination stream will be removed before being copied.
void CopyFromHoldup(double TimeDst, CHoldup *SrcHoldup, double TimeSrc, double MassFlow, bool DeleteDataAfter = true)
Copies all stream data from time point TimeSrc
of holdup SrcHoldup
to the time point TimeDst
of this material stream with setting of new mass flow MassFlow
.
If flag DeleteDataAfter
is set to true
, all data after the time point Time
in the destination stream will be removed before being copied.
void AddStream (CMaterialStream *Stream, double Time)
Performs a mixing of this material stream with material stream Stream
for the specified time point Time
.
void AddStream (CMaterialStream *Stream, double Start, double End, unsigned TPType = BOTH_TP)
Performs a mixing of this material stream with material stream Stream
for the specified time interval. Boundary points Start
and End
are included into this interval.
Parameter TPType
specifies which time points will be present in the resulting stream: combining points from two streams (BOTH_TP
), only from the first stream (DST_TP
), or only from the second stream (SRC_TP
).
Data for nonexistent points are obtained by linear interpolation.
Following functions are for class CHoldup
.
void CopyFromHoldup(CHoldup *SrcHoldup, double Time, bool DeleteDataAfter = true)
Copies all holdup data from SrcHoldup
for the specified time point Time
to the current holdup.
If flag DeleteDataAfter
is set to true
, all data after the time point Time
in the destination holdup will be removed before being copied.
void CopyFromHoldup(CHoldup *SrcHoldup, double Start, double End, bool DeleteDataAfter = true)
Copies all holdup data from SrcHoldup
on the certain time interval to the current holdup. Boundary points Start
and End
are included into this interval.
If flag DeleteDataAfter
is set to true
, all data after the time point Time
in the destination holdup will be removed before being copied.
void CopyFromHoldup(double TimeDst, CHoldup *SrcHoldup, double TimeSrc, bool DeleteDataAfter = true)
Copies all holdup data from time point TimeSrc
of holdup SrcHoldup
to the time point TimeDst
of this holdup.
If flag DeleteDataAfter
is set to true
, all data after the time point Time
in the destination holdup will be removed before being copied.
void AddHoldup(CHoldup *Holdup, double Time)
Performs a mixing of this holdup with holdup Holdup
for the specified time point Time
.
void AddHoldup(CHoldup *Holdup, double Start, double End, unsigned TPType)
Performs a mixing of this holdup with Holdup
for the specified time interval. Boundary points Start
and End
are included into this interval.
Parameter TPType
specifies which time points will be present in the resulting holdup: combining points from two holdups (BOTH_TP
), only from the first holdup (DST_TP
), or only from the second holdup (SRC_TP
).
Data for nonexistent points are obtained by linear approximation.
void AddStream(CMaterialStream *Stream, double Start, double End)
Performs a mixing of this holdup with material stream Stream
for the specified time interval. Boundary points Start
and End
are included into this interval.
Data for nonexistent points are obtained by linear interpolation.
Particle size distribution
Several global functions are defined to work with particle size distributions. These functions can be called from any place of the code.
All functions receive grid (Grid
) as the input parameter. The grid can be previously obtained with the help of the function GetNumericGrid
, for more information please refer to getNumericGrid in section Basic unit.
Note
Notations:
\(d_i\) – diameter of particle in class \(i\)
\(\Delta d_i\) – size of the class \(i\)
\(M_k\) – \(k\)th moment
\(q\) – density distribution
\(q_0\) – number related density distribution
\(Q_0\) – number related cumulative distribution
\(q_2\) – surfacearearelated density distribution
\(Q_2\) – surfacearearelated cumulative distribution
\(q_3\) – massrelated density distribution
\(Q_3\) – massrelated cumulative distribution
\(w_i\) – mass fraction of particles of class \(i\)
\(N_i\) – number of particles of class \(i\)
\(N_{tot}\) – total number of particles
double GetMMoment(Moment, Grid, InDistr)
Calculates moment of the density distribution by \(M_k = \sum\limits_i d_i^k \, q_i \, \Delta d_i\).
vector<double> ConvertQ0Toq0(Grid, InDistr)
Performs conversion from \(Q_0\) to \(q_0\) distributions using information about the size grid: \(q_{0,0} = \dfrac{Q_{0,0}}{\Delta d_i}\) and \(q_{0,i} = \dfrac{Q_{0,i}  Q_{0,i1}}{\Delta d_i}\).
vector<double> Convertq0ToQ0(Grid, InDistr)
Performs conversion from \(q_0\) to \(Q_0\) distributions using information about the size grid: \(Q_{0,i} = \sum\limits_i q_{0,i} \, \Delta d_i = Q_{0,i1} + q_{0,i} \, \Delta d_i\).
vector<double> ConvertQ2Toq2(Grid, InDistr)
Performs conversion from \(Q_2\) to \(q_2\) distributions using information about the size grid: \(q_{2,0} = \dfrac{Q_{2,0}}{\Delta d_i}\) and \(q_{2,i} = \dfrac{Q_{2,i}  Q_{2,i1}}{\Delta d_i}\).
vector<double> Convertq2ToQ2(Grid, InDistr)
Performs conversion from \(q_2\) to \(Q_2\) distributions using information about the size grid: \(Q_{2,i} = \sum\limits_i q_{2,i} \, \Delta d_i = Q_{2,i1} + q_{2,i} \, \Delta d_i\).
vector<double> ConvertQ3Toq3(Grid, InDistr)
Performs conversion from \(Q_3\) to \(q_3\) distributions using information about the size grid: \(q_{3,0} = \dfrac{Q_{3,0}}{\Delta d_i}\) and \(q_{3,i} = \dfrac{Q_{3,i}  Q_{3,i1}}{\Delta d_i}\).
vector<double> Convertq3ToQ3(Grid, InDistr)
Performs conversion from \(q_3\) to \(Q_3\) distributions using information about the size grid: \(Q_{3,i} = \sum\limits_i q_{3,i} \, \Delta d_i = Q_{3,i1} + q_{3,i} \, \Delta d_i\).
vector<double> Convertq0Toq2(Grid, InDistr)
Performs conversion from \(q_0\) to \(q_2\) distributions using information about the size grid by \(q_{2,i} = \dfrac{d_i^2 \, q_{0,i}}{M_2(q_0)}\).
vector<double> Convertq0Toq3(Grid, InDistr)
Performs conversion from \(q_0\) to \(q_3\) distributions using information about the size grid by \(q_{3,i} = \dfrac{d_i^3 \, q_{0,i}}{M_3(q_0)}\).
vector<double> Convertq2Toq0(Grid, InDistr)
Performs conversion from \(q_2\) to \(q_0\) distributions using information about the size grid by \(q_{0,i} = \dfrac{d_i^{2} \, q_{2,i}}{M_{2}(q_2)}\).
vector<double> Convertq2Toq3(Grid, InDistr)
Performs conversion from \(q_2\) to \(q_3\) distributions using information about the size grid by \(q_{3,i} = \dfrac{d_i \, q_{2,i}}{M_1(q_2)}\).
vector<double> Convertq3Toq0(Grid, InDistr)
Performs conversion from \(q_3\) to \(q_0\) distributions using information about the size grid by \(q_{0,i} = \dfrac{d_i^{3} \, q_{3,i}}{M_{3}(q_3)}\).
vector<double> Convertq3Toq2(Grid, InDistr)
Performs conversion from \(q_3\) to \(q_2\) distributions using information about the size grid by \(q_{2,i} = \dfrac{d_i^{1} \, q_{3,i}}{M_{1}(q_3)}\).
vector<double> ConvertMassFractionsToq3(Grid, InDistr)
Calculates \(q_3\) distribution using the size grid and the distribution of mass fractions by \(q_3 = w_i / \Delta d_i\).
vector<double> ConvertMassFractionsToQ3(InDistr)
Calculates \(Q_3\) distribution using the distribution of mass fractions: \(Q_{3,0} = w_i\) and \(Q_{3,i} = Q_{3,i1} + w_i\).
vector<double> ConvertMassFractionsToq0(Grid, InDistr)
Calculates \(q_0\) distribution using the functions ConvertMassFractionsToq3 and Convertq3Toq0.
vector<double> ConvertMassFractionsToQ0(Grid, InDistr)
Calculates \(Q_0\) distribution using the functions ConvertMassFractionsToq0 and Convertq0ToQ0.
vector<double> ConvertMassFractionsToq2(Grid, InDistr)
Calculates \(q_0\) distribution using the functions ConvertMassFractionsToq3 and Convertq3Toq2.
vector<double> ConvertMassFractionsToQ2(Grid, InDistr)
Calculates \(Q_0\) distribution using the functions ConvertMassFractionsToq2 and Convertq2ToQ2.
vector<double> Convertq3ToMassFractions(Grid, InDistr)
Calculates mass fractions from \(q_3\) distribution using the size grid by \(w_i = q_{3,i}\cdot \Delta d_i\).
vector<double> ConvertQ3ToMassFractions(InDistr)
Calculates mass fractions from \(Q_3\) distribution using the size grid: \(w_0 = Q_{3,0}\) and \(w_i = Q_{3,i}  Q_{3,i1}\).
vector<double> Convertq0ToMassFractions(Grid, InDistr)
Calculates mass fractions from \(q_0\) distribution using the functions Convertq0Toq3 and Convertq3ToMassFractions.
vector<double> ConvertQ0ToMassFractions(Grid, InDistr)
Calculates mass fractions from \(Q_0\) distribution using the functions ConvertQ0Toq0 and Convertq0ToMassFractions.
vector<double> Convertq2ToMassFractions(Grid, InDistr)
Calculates mass fractions from \(q_2\) distribution using the functions Convertq2Toq3 and Convertq3ToMassFractions.
vector<double> ConvertQ2ToMassFractions(Grid, InDistr)
Calculates mass fractions from \(Q_2\) distribution using the functions ConvertQ2Toq2 and Convertq2ToMassFractions.
vector<double> ConvertNumbersToq0(Grid, InDistr)
Calculates \(q_0\) distribution using the number distribution and the size grid by \(q_{0,i} = \dfrac{N_i}{\Delta d_i \, N_{tot}}\).
vector<double> ConvertNumbersToQ2(Grid, InDistr)
Calculates \(Q_2\) distribution using the number distribution and the size grid by \(Q_{2,i} = \dfrac{\sum\limits_{j=0}^i N_j \, \pi\,d_j^2}{\sum\limits_j N_j\, \pi\,d_j^2}\).
vector<double> ConvertNumbersToQ0(Grid, InDistr)
Calculates \(Q_0\) distribution using the number distribution and the functions ConvertNumbersToq0 and Convertq0ToQ0.
vector<double> ConvertNumbersToq2(Grid, InDistr)
Calculates \(q_2\) distribution using the number distribution and the functions ConvertNumbersToQ2 and ConvertQ2Toq2.
vector<double> ConvertNumbersToq3(Grid, InDistr)
Calculates \(q_3\) distribution using the number distribution and the functions ConvertNumbersToq0 and Convertq0Toq3.
vector<double> ConvertNumbersToQ3(Grid, InDistr)
Calculates \(Q_3\) distribution using the number distribution and the functions ConvertNumbersToq3 and Convertq3ToQ3.
vector<double> ConvertNumbersToMassFractions(Grid, InDistr)
Calculates mass fractions from the number distribution using the functions ConvertNumberToq0 and Convertq0ToMassFractions.
vector<double> Convertq0Toq0(OldGrid, OldDistr, NewGrid)
Converts \(q_0\) distribution to the same distribution on the modified size grid.
vector<double> Convertq2Toq2(OldGrid, OldDistr, NewGrid)
Converts \(q_2\) distributions to the same distribution on the modified size grid.
vector<double> Convertq3Toq3(OldGrid, OldDistr, NewGrid)
Converts \(q_3\) distributions to the same distribution on the modified size grid.
NormalizeDensityDistribution(Grid, qiDistr)
Normalizes density distribution \(q_0\) or \(q_3\) by \(q_i = \dfrac{q_i}{\sum\limits_j q_j\,\Delta d_j}\).
double GetDistributionMedian(Grid, QxDistr)
Returns median in [m] of \(Q_0\) or \(Q_3\) distribution. Median is a diameter, which corresponds to a value of distribution equal to 0.5.
double GetDistributionValue (Grid, QxDistr, Val)
Returns diameter in [m], which corresponds to a specified value of cumulative distribution \(Q_0\) or \(Q_3\). Input value Val
should range between 0 and 1.
double GetDistributionMode (Grid, qxDistr)
Returns diameter in [m], which corresponds to a maximum value of density distribution.
double GetAverageDiameter (Grid, qxDistr)
Returns average diameter in [m] of the distribution \(q_0\) or \(q_3\).
double GetSauterDiameter (Grid, q3Distr)
Calculates Sauter diameter (\(d_{32}\)) of \(q_3\) distribution in [m].
double GetSpecificSurface (Grid, q3Distr)
Calculates specific surface of \(q_3\) distribution in [m \(^2\)].
Thermodynamics
In Dyssol, thermodynamic functions are applied to calculate a material stream / holdup temperature \(T\) from a given enthalpy value \(H\). In most cases, the enthalpy of a compound correlates with temperature nonlinearly. Thus, for backward calculation of the temperature from an enthalpy value, you need to solve a nonlinear equation \(T = f^{1}(H)\).
This solving step is replaced by a lookup table functionality, i.e. codependent temperature and enthalpy are precalculated for a certain range of values and then stored in a table. This table can be used to determine temperatureenthalpy value pairs by interpolation between two neighbor points.
To add more functionality, these lookup tables can be called for every property defined in the materials database.
The temperature or pressure (e.g. pressure from saturation pressure or different properties) is then calculated in the following steps:
Create tables of value pairs for each compound that is present in the system, i.e. temperature / pressure and respective dependent property value over a certain range of the temperature / pressure.
Combine the compoundlookup tables by massweighted summation of the tables.
Readout of temperature / pressure from the resulting lookuptable at a certain point of the property by interpolation between neighbor value pairs.
The functions in base unit are less timeconsuming than the material stream / holdupfunction and therefore are implemented for usage during solving process (i.e. function CalculateResiduals), if the composition of the material stream / holdup varies in each iteration.
In Unit library, Heater
and HeatExchanger
apply thermodynamic calculations.
For material stream and holdup
double CalcTemperatureFromProperty(ECompoundTPProperties Property, double Time, double Value)
Returns temperature of the material stream / holdup for a specific value Value
of the property Property
at the time point Time
. Possible properties are those defined in Material database.
double CalcPressureFromProperty(ECompoundTPProperties Property, double Time, double Value)
Returns pressure of the material stream / holdup for a specific value Value
of the property Property
at the time point Time
. Possible properties are those defined in Material database.
For base unit
double CalcTemperatureFromProperty(ECompoundTPProperties Property,vector<double>& CompoundFractions, double Value)
Returns temperature of a generic system of composition CompoundFractions
for a specific value Value
of the property Property
. Possible properties are those defined in Material database.
double CalcPressureFromProperty(ECompoundTPProperties Property,vector<double>& CompoundFractions, double Value)
Returns pressure of a generic system of composition CompoundFractions
for a specific value Value
of the property Property
. Possible properties are those defined in Material database.
void HeatExchange(CMaterialStream* Stream1, CMaterialStream* Stream2, double Time, double Efficiency);
Performs a heat exchange between material streams Stream1
and Stream2
at specified time point Time
with a specified efficiency (0 ≤ Efficiency
≤ 1).
\(\varepsilon\) stands for efficiency and \(\theta\) for temperature as an integration variable.
External solver
As mentioned in section Solver development, you can generate your own solvers in Dyssol under class CExternalSolver
.
In Dyssol, external solvers are connected with opensource programm packages. Thus you can only do some basic operations for these solvers, such as generate new solver, name it, get ID and version, initialize and finalize, etc.
In this section, the functions and variables applied in external solver are introduced.
CExternalSolver()
Basic constructor of the solver. Creates an empty external solver. Called only once when solver is added to the unit.
Internal variables are:
m_solverName
: name of the solver that will be displayed in user interface of Dyssol.m_authorName
: solver’s authorm_solverUniqueKey
: unique identificator of the solver. Simulation environment distinguishes different solvers with the help of this identificator. You must ensure that ID of your solver is unique. This ID can be created manually or using GUIDgenerator of Visual Studio (Tools → GUID Genarator).m_solverVersion
: current version of the solver.
Basic information
std::string GetName()
Returns name of the solver.
std::string GetUniqueID()
Returns string key, which is unique among all solvers.
std::string GetAuthorName()
Returns name of the solver’s author.
unsigned GetVersion()
Returns version of the solver.
Virtual functions
virtual void Initialize()
Solver‘s initialization. This function is called only once for each simulation during the initialization of unit. Implementation of this function is not obligatory and can be skipped.
virtual void Finalize()
Unit‘s finalization. This function is called only once for each simulation during the finalization of unit. Implementation of this function is not obligatory and can be skipped.
DAE Systems
In Dyssol, you can solve systems of DAE automatically. In this case, the unit should contain one or several additional objects of class CDAEModel
. This class is used to describe DAE systems and can be automatically solved with class CDAESolver
.
DAE model
CDAEModel()
Basic constructor. Creates an empty DAE model.
Variables
unsigned AddDAEVariable(bool _bIsDifferentiable, double _dVariableInit, double _dDerivativeInit, double _dConstraint)
Adds new algebraic (_bIsDifferentiable = false
) or differential (_bIsDifferentiable = true
) variable with initial values _dVariableInit
and _dDerivativeInit
.
Should be called in function Initialize(Time)
of the unit. Returns unique index of added variable.
_dConstraint
sets the constraint for the variable, different values of _dConstraint
are defined as follows:
0.0: no constraint
1.0: Variable ≥ 0.0
−1.0: Variable ≤ 0.0
2.0: Variable > 0.0
−2.0: Variable < 0.0
unsigned GetVariablesNumber()
Returns current number of defined variables.
void ClearVariables()
Removes all defined variables from the model.
Tolerance
void SetTolerance(double _dRTol, double _dATol)
Sets values of relative and absolute tolerances for all variables.
void SetTolerance(double _dRTol, std::vector<double> &_vATol)
Sets value of relative tolerance for all variables and absolute tolerances for each variable separately.
double GetRTol()
Returns current relative tolerance.
double GetATol(unsigned _dIndex)
Returns current absolute tolerance for specified variable.
Virtual functions
virtual void CalculateResiduals(double _dTime, double *_pVars, double *_pDerivs, double *_pRes, void *_pUserData)
Computes the problem residual for given values of the independent variable _dTime
, state vector _pVars
, and derivatives _pDerivs
. Here the DAE system should be specified in implicit form. Function will be called by solver automatically.
_dTime
: current value of the independent variable \(t\).
_pVars
: pointer to an array of the dependent variables, \(y(t)\).
_pDerivs
: pointer to an array of derivatives \(y'(t)\).
_pRes
: output residual vector \(F(t, y, y')\).
_pUserData
: pointer to user’s data. Is used to provide access from this function to unit’s data.
virtual void ResultsHandler(double _dTime, double *_pVars, double *_pDerivs, void *_pUserData)
Processing the results returned by the solver at each calculated step. Called by solver every time when the solution in new time point is ready.
_dTime
: current value of the independent variable \(t\).
_pVars
: current values of the dependent variables, \(y(t)\).
_pDerivs
: current values of derivatives \(y'(t)\).
_pUserData
: pointer to user’s data. Is used to provide access from this function to unit’s data.
Other functions
void Clear()
Removes all data from the model.
void SetUserData(void *_pUserData)
Set pointer to user’s data. This data will be returned in overloaded functions CalculateResiduals and ResultsHandler. Usually is used to provide access from these functions to unit’s data.
DAE solver
CDAESolver()
Basic constructor. Creates an empty solver.
Model
bool SetModel(CDAEModel *_pModel)
Sets model to a solver. Should be called in function Initialize of the unit.
Returns false
on error.
bool SetMaxStep()
Sets maximum time step for solver. Should be used in Initialize before the function SetModel.
bool Calculate(double _dTime, unsigned _nModel = KIN_NONE)
Solves problem on a given time point _dTime
. Should be called in function Simulate of the unit.
Default value of _nModel
is KIN_NONE
, a Newtonbased iteration method. Other available models can be found in declaration for KINSOL parameters in the *.cpp
file, including: KIN_FP
(fixed point), KIN_LINESEARCH
and KIN_PICARD
.
Returns false
on error.
bool Calculate(double _dStartTime, double _dEndTime)
Solves problem on a given time interval. Should be called in function Simulate of the unit. Returns false
on error.
Other functions
void SaveState()
Saves current state of the solver. Should be called in function SaveState of the unit.
void LoadState()
Loads last saved state of the solver. Should be called in function LoadState of the unit.
std::string GetError()
Returns error’s description. Can be called if function SetModel or Calculate returns false
.
Application example
Assume that you are going to solve a dynamic DAE system,
where \(x\), \(y\), \(z\) are fractions of solid, liquid and vapor phases of the output stream of the unit.
Now you want to develope a new unit for for automatic calculation of this DAE by using builtin solvers of Dyssol, just do the following steps:
Add a new template unit to your solution and name it
DynamicUnitWithSolver
, please refer to Unit development.In file
Unit.h
, there is already a class of DAE modelCMyDAEModel
defined with:two functions, which must be overridden (CalculateResiduals and ResultsHandler) and
a class of unit
CUnit
with two additional variables (for DAE modelCMyDAEModel m_Model
, for DAE solverCDAESolver m_Solver
).
Add three variables in class
CMyDAEModel
to store indices for \(x\), \(y\) and \(z\) variables, name themm_nS
,m_nL
andm_nV
respectively.After adding description of class
CMyDAEModel
, your code should look like this:class CMyDAEModel : public CDAEModel { public: unsigned m_nS; //solid fraction unsigned m_nL; //liquid fraction unsigned m_nV; //vapor fraction public: void CalculateResiduals(double _dTime, double* _pVars, double* _pDers, double* _pRes, void* _pUserData); void ResultsHandler(double _dTime, double* _pVars, double* _pDers, void *_pUserData); };
Setup unit’s basic info (name, author’s name, unique ID, ports) as described in Unit development.Then provide model with pointer to your unit, in order to have an access to the unit from functions of the model. Now your unit’s constructor should look like this:
CUnit::CUnit() { // Unit basic information m_sUnitName = "DummyUnit4"; m_sAuthorName = "Author"; m_sUniqueID = "344BCC0048AA4c3a9117F20A9F8AF9A8"; //an example ID // Add ports for in and outlets AddPort("InPort", INPUT_PORT); AddPort("OutPort", OUTPUT_PORT); // Set this unit as user data of the model applied m_Model.SetUserData(this); }
Implement function
Initialize
of the unit:Since function
Initialize
is called every time when simulation starts, all variables must be previously removed from the model by callingClearVariables
.m_Model.ClearVariables();
Now you can add 3 variables with specified initial conditions to the model (using function AddDAEVariable) according to the equation system. Use phase fractions of the input stream as initials.
Set tolerances to the model using function SetTolerance. As tolerances for the model, global tolerances of the system can be used.
m_Model.SetTolerance( GetRelTolerance(), GetAbsTolerance() );
Now, you can connect your model to the solver by calling function SetModel. To receive errors from solver, connect it to the global errors handling procedure.
if( !m_Solver.SetModel(&m_Model) ) RaiseError(m_Solver.GetError());
Your code should look like this after following all steps above:
void CUnit::Initialize(double _dTime) { // Get pointer to inlet stream CMaterialStream* pInpStream = GetPortStream("InPort"); // Clear previous state variables in model m_Model.ClearVariables(); // Add state variables to model double dSFr = pInpStream>GetSinglePhaseProp(_dTime, FRACTION, SOA_SOLID); double dLFr = pInpStream>GetSinglePhaseProp(_dTime, FRACTION, SOA_LIQUID); double dVFr = pInpStream>GetSinglePhaseProp(_dTime, FRACTION, SOA_VAPOR); m_Model.m_nS = m_Model.AddDAEVariable(true, dSFr, 0.04, 1.0); m_Model.m_nL = m_Model.AddDAEVariable(true, dLFr, 0.04, 1.0); m_Model.m_nV = m_Model.AddDAEVariable(false, dVFr, 0, 1.0); // Set tolerance to model m_Model.SetTolerance( GetRelTolerance(), GetAbsTolerance() ); // Set model to solver if( !m_Solver.SetModel(&m_Model) ) RaiseError(m_Solver.GetError()); }
Connect solver to a system saving / loading procedure in functions
SaveState()
andLoadState()
of the unit:void CUnit::SaveState() { m_Solver.SaveState(); } void CUnit::LoadState() { m_Solver.LoadState(); }
In function
Simulate
of the unit, calculation procedure should be run by calling functionCalculate
of the solver. Additionally, solver can be connected to the system’s errors handling procedure to receive possible errors during the calculation. Unit’sSimulate
function after that must look as follows:void CUnit::Simulate(double _dStartTime, double _dEndTime) { // Get pointers to inlet and outlet streams CMaterialStream *pInputStream = GetPortStream("InPort"); CMaterialStream *pOutputStream = GetPortStream("OutPort"); pOutputStream>RemoveTimePointsAfter(_dStartTime, true); pOutputStream>CopyFromStream(pInputStream, _dStartTime); // Copy the inlet stream information to outlet stream // Run solver and check if errors take place if( !m_Solver.Calculate(_dStartTime, _dEndTime) ) RaiseError( m_Solver.GetError() ); }
In function
CalculateResiduals
, DAE system in implicit form should be described. According to the given equation system\[ \begin{align}\begin{aligned}\begin{cases} \dfrac{dx}{dt} = 0.04 \cdot x + 10^4 \cdot y \cdot z &x(0) = 0.04\\ \dfrac{dy}{dt} = 0.04\cdot x  10^4 \cdot y \cdot z  3 \cdot 10^7 \cdot y^2 &y(0) = 0.04\\ x + y + z = 1 \end{cases}\end{aligned}\end{align} \]and the definition of residual
_pRes
(difference between new and old values from last calculation), the equaitons are expressed below. Here the values are arranged in vector_pVars
according to sequence \([x, y, z]\), in_pDers
according to \([dx, dy]\).void CMyDAEModel::CalculateResiduals(double _dTime, double* _pVars, double* _pDers, double* _pRes, void* _pUserData) { _pRes[0] = _pDers[0]  (0.04*_pVars[0] + 1.0e4*_pVars[1]*_pVars[2]); _pRes[1] = _pDers[1]  ( 0.04*_pVars[0]  1.0e4*_pVars[1]*_pVars[2]  3.0e7*_pVars[1]*_pVars[1] ); _pRes[2] = _pVars[0] + _pVars[1] + _pVars[2]  1; }
Last step is handling of results from the solver (
CMyDAEModel::ResultsHandler
). Calculated fractions can be set here to the output stream of the unit. To access to the unit’s data, use the pointer_pUserData
defined previously:void CMyDAEModel::ResultsHandler(double _dTime, double* _pVars, double* _pDerivs, void *_pUserData) { // Get pointers to streams CUnit *unit = static_cast<CUnit*>(_pUserData); CMaterialStream *pStream = static_cast<CMaterialStream*>(unit>GetPortStream("OutPort")); // Add time point to outlet stream pStream>AddTimePoint(_dTime); // Set calculated results to corresponding phases pStream>SetSinglePhaseProp(_dTime, FRACTION, SOA_SOLID, _pVars[0]); pStream>SetSinglePhaseProp(_dTime, FRACTION, SOA_LIQUID, _pVars[1]); pStream>SetSinglePhaseProp(_dTime, FRACTION, SOA_VAPOR, _pVars[2]); }
Matrices
Several types of matrix classes, including the following types, are introduced in this section.
Transform matrix:
CTransformMatrix
.Dense 2dimensional matrix:
CDense2DMatrix
.Dense multidimensional matrix:
CDenseMDMatrix
.
Transformation matrix
Transformation matrices (class CTransformMatrix
) are applied to describe laws of changes for multidimensional distributions. Each cell of matrix describes how much of material will be transferred from one class of multidimensional distribution to another.
CTransformMatrix()
Basic constructor. Creates an empty matrix.
CTransformMatrix (unsigned _nType, unsigned _nClasses)
Creates matrix to transform onedimensional distribution with type _nType
and _nClasses
classes. _nType
is one of the predefined types of solid distributions.
All values in matrix will be set to 0.
CTransformMatrix(unsigned _nType1, unsigned _nClasses1, unsigned _nType2, unsigned _nClasses2)
Creates matrix to transform twodimensional distribution with types _nType1
and _nType2
and classes _nClasses1
and _nClasses2
. _nType1
and _nType2
are types from predefined types of solid distributions.
All values in matrix will be set to 0.
CTransformMatrix(const std::vector<unsigned> &_vTypes, const std::vector<unsigned> &_vClasses)
Creates transformation matrix for distribution with specified types and classes. _vTypes
and _vClasses
must have the same length. _vTypes
is the vector of types from predefined types of solid distributions.
All values in matrix will be set to 0.
Dimensions
bool SetDimensions(unsigned _nType, unsigned _nClasses)
Sets new dimensions set to the matrix in order to transform onedimensional distribution with type _nType
and _nClasses
classes. _nType
is one of the predefined types of solid distributions.
Old data will be erased and matrix will be initialized with zeroes. Returns false
on error.
bool SetDimensions(unsigned _nType1, unsigned _nClasses1, unsigned _nType2, unsigned _nClasses2)
Sets new dimensions set to the matrix in order to transform twodimensional distribution.
_nType1
and _nType2
are types of the predefined types of solid distributions. Types must be unique. _nClasses1
and _nClasses2
are number of classes in corresponding distributions.
Old data will be erased and matrix will be initialized with zeroes. Returns false
on error.
bool SetDimensions(unsigned _nType1, unsigned _nClasses1, unsigned _nType2, unsigned _nClasses2, unsigned _nType3, unsigned _nClasses3)
Sets new dimensions set to the matrix in order to transform threedimensional distribution.
_nType1
, _nType2
and _nType3
are one of the predefined types of solid distributions. Types must be unique. _nClasses1
, _nClasses2
and _nClasses3
are number of classes in corresponding distributions.
Old data will be erased and matrix will be initialized with zeroes. Returns false
on error.
bool SetDimensions(const std::vector<unsigned> &_vTypes, const std::vector<unsigned> &_vClasses)
Sets new dimensions set with types _vTypes
and numbers of classes _vClasses
. _vTypes
is the vector of predefined types of solid distributions.
All old data will be erased and matrix will be initialized with zeroes. Sizes of vectors _vTypes
and _vClasses
must be equal. Returns false
on error.
std::vector<unsigned> GetDimensions()
Returns vector with all current defined dimensions types.
std::vector<unsigned> GetClasses()
Returns vector with current numbers of classes.
unsigned GetDimensionsNumber()
Returns current number of dimensions.
Get data
double GetValue(unsigned _nCoordSrc, unsigned _nCoordDst)
Returns value by specified coordinates according to all defined dimensions in transformation matrix for onedimensional distribution. _nCoordSrc
is coordinate of a source class, _nCoordDst
is coordinate of a destination class.
Returning value is a mass fraction, which will be transferred from the source class to the destination class. Works with onedimensional distribution only. Returns 1
on error.
double GetValue(unsigned _nCoordSrc1, unsigned _nCoordSrc2, unsigned _nCoordDst1, unsigned _nCoordDst2)
Returns value by specified coordinates according to all defined dimensions in transformation matrix for twodimensional distribution. _nCoordSrc1
and _nCoordSrc2
are coordinates of a source class, _nCoordDst1
and _nCoordDst2
are coordinate of a destination class.
Returning value is a mass fraction, which will be transferred from the source class to the destination class. Works with twodimensional distribution only. Returns 1
on error.
double GetValue(const std::vector<unsigned> &_vCoordsSrc, const std::vector<unsigned> &_vCoordsDst)
Returns value by specified coordinates according to all defined dimensions. _vCoordsSrc
are coordinates of a source class, _vCoordsDst
are coordinates of a destination class. Sizes of vectors _vCoordsSrc
and _vCoordsDst
must be equal and must correspond to the number of currently defined dimensions.
Returning value is a mass fraction, which will be transferred from the source class to the destination class. Returns 1
on error.
double GetValue(const std::vector<unsigned> &_vDimsSrc, const std::vector<unsigned> &_vCoordsSrc, const std::vector<unsigned> &_vDimsDst, const std::vector<unsigned> &_vCoordsDst)
Returns value according to specified coordinates and dimensions. Number of dimensions must be the same as defined in the transformation matrix, but their sequence can be different. Sizes of all vectors must be equal.
Returning value is a mass fraction, which will be transferred from the source class to the destination class. Returns 1
on error.
bool GetVectorValue(const std::vector<unsigned> &_vCoordsSrc, const std::vector<unsigned> &_vCoordsDst, std::vector<double> &_vResult)
Returns vector value by specified coordinates according to all defined dimensions. Size of one vector of coordinates must be equal to the number of dimensions in transformation matrix; size of the second one must be one less.
Returning value _vResult
is a vector of mass fractions, which will be transferred from the source to the destination. Returns false
on error.
bool GetVectorValue(const std::vector<unsigned> &_vDimsSrc, const std::vector<unsigned> &_vCoordsSrc, const std::vector<unsigned> &_vDimsDst, const std::vector<unsigned> &_vCoordsDst, std::vector<double> &_vResult)
Returns vector of values according to specified coordinates and dimensions sequence. Number of dimensions must be the same as defined in the transformation matrix, but their sequence can be different. Size of one vector of coordinates must be equal to the number of dimensions in transformation matrix; size of the second one must be one less.
Returning value _vResult
is a vector of mass fractions, which will be transferred from the source to the destination. Returns false
on error.
Set data
bool SetValue(unsigned _nCoordSrc, unsigned _nCoordDst, double _dValue)
Sets value by specified coordinates for onedimensional distribution. _nCoordSrc
is a coordinate of the source class; _nCoordDst
is a coordinate of the destination class. _dValue
is a mass fraction, which will be transferred from the source class to the destination class.
Returns false
on error.
bool SetValue (unsigned _nCoordSrc1, unsigned _nCoordSrc2, unsigned _nCoordDst1, unsigned _nCoordDst2, double _dValue)
Sets value by specified coordinates for twodimensional distribution.
_nCoordSrc1
and _nCoordSrc2
are coordinate of the source class; _nCoordDst1
and _nCoordDst2
are coordinate of the destination class. _dValue
is a mass fraction, which will be transferred from the source class to the destination class.
Returns false
on error.
bool SetValue (const std::vector<unsigned> &_vCoordsSrc, const std::vector<unsigned> &_vCoordsDst, double _dValue)
Sets value by specified coordinates and full dimensions set. _vCoordsSrc
are coordinates of the source class, _vCoordsDst
are coordinates of the destination class. Sizes of vectors _vCoordsSrc
and _vCoordsDst
must be equal. _dValue
is a mass fraction, which will be transferred from the source class to the destination class.
Returns false
on error.
bool SetValue (const std::vector<unsigned> &_vDimsSrc, const std::vector<unsigned> &_vCoordsSrc, const std::vector<unsigned> &_vDimsDst, const std::vector<unsigned> &_vCoordsDst, double _dValue)
Sets value according to specified coordinates and dimensions. Number of dimensions must be the same as defined in the transformation matrix, but their sequence can be different. Sizes of all vectors must be equal. _dValue
is a mass fraction, which will be transferred from the source class to the destination class.
Returns false
on error.
bool SetVectorValue (const std::vector<unsigned> &_vCoordsSrc, const std::vector<unsigned> &_vCoordsDst, const std::vector<double> &_vValue)
Sets vector of values by specified coordinates according to all defined dimensions. Size of one vector of coordinates must be equal to the number of dimensions in transformation matrix; size of the second one must be one less. _vValue
is a vector of mass fractions, which will be transferred from the source to the destination.
Returns false
on error.
bool SetVectorValue (const std::vector<unsigned> &_vDimsSrc, const std::vector<unsigned> &_vCoordsSrc, const std::vector<unsigned> &_vDimsDst, const std::vector<unsigned> &_vCoordsDst, const std::vector<double> &_vValue)
Sets vector of values according to specified coordinates and dimensions sequence. Number of dimensions must be the same as defined in the transformation matrix, but their sequence can be different. Size of one vector of coordinates must be equal to the number of dimensions in transformation matrix; size of the second one must be one less. _vValue
is a vector of mass fractions, which will be transferred from the source to the destination.
Returns false
on error.
Other functions
void Normalize()
Normalizes data in matrix: sets sum of material which transfers from each single class to 1.
void ClearData()
Sets all data in matrix to 0.
void Clear()
Removes all data and information about dimensions from the matrix.
Twodimensional matrix
Basic information and functions of class CDense2DMatrix
are introduced below.
CDense2DMatrix()
Basic constructor. Creates empty matrix with zero in all rows and columns.
CDense2DMatrix(CDense2DMatrix &_matrix)
Copy constructor. Creates matrix with the same dimensions as in _matrix
and copies all data.
CDense2DMatrix(unsigned _nRows, unsigned _nColumns)
Creates new matrix with specified number of rows and columns. All data will be set to 0.
Dimensions
void SetDimensions(unsigned _nRows, unsigned _nColumns)
Sets new dimensions to the matrix. Old data is removed and all entries are set to zero.
unsigned GetRowsNumber()
Returns number of rows in the matrix.
unsigned GetColumnsNumber()
Returns number of columns in the matrix.
Get data
double GetValue(unsigned _nRow, unsigned _nColumn)
Returns data by the specified indexes. Returns 0
if such indexes do not exist.
std::vector<double> GetRow(unsigned _nRow)
Returns vector of data for specified row. Returns empty vector if such row does not exist.
std::vector<double> GetColumn(unsigned _nColumn)
Returns vector of data for specified column. Returns empty vector if such column does not exist.
std::vector<std::vector<double>> GetMatrix()
Returns all data in form of vectorofvectors.
Set data
void SetValue(unsigned _nRow, unsigned _ nColumn, double _dValue)
Sets data _dValue
by the specified indexes.
void SetRow(unsigned _nRow, const std::vector<double>& _Values)
Sets data _Values
to a specified row.
void SetColumn(unsigned _nColumn, const std::vector<double>& _Values)
Sets data _Values
to a specified column.
void SetMatrix(const std::vector<std::vector<double>>& _matr)
Sets all values in form vectorofvectors _matr
to matrix. _matr
must have the same dimensions as the matrix itself.
Overloaded operators
CDense2DMatrix operator+(const CDense2DMatrix& _Matrix1, const CDense2DMatrix &_Matrix2)
Performs addition of two matrices with the same dimensions. Returns an empty matrix in case of different dimensions.
CDense2DMatrix operator(const CDense2DMatrix &_Matrix1, const CDense2DMatrix &_Matrix2)
Performs subtraction of two matrices with the same dimensions. Returns an empty matrix in case of different dimensions.
CDense2DMatrix operator*(double _dMultiplier)
Performs multiplication of the matrix with a coefficient _dMultiplier
.
CDense2DMatrix& operator=(const CDense2DMatrix &_matrix)
Sets dimenions and data from the _matrix
to a left matrix.
Other functions
void Normalize()
Normalizes the matrix so that the sum of all elements equals to 1.
void ClearData()
Sets all elements in matrix to 0.
void Clear()
Removes all data and sets number of rows and columns equal to 0.
Multidimensional matrix
Basic information and functions of class CDenseMDMatrix
are introduced below.
CDenseMDMatrix()
Basic constructor. Creates an empty matrix.
CDenseMDMatrix(const CDenseMDMatrix &_source)
Copy constructor. Creates matrix with the same dimensions as in _source
and copies all data from there.
Dimensions
bool SetDimensions(unsigned _nType, unsigned _nClasses)
Sets onedimensional distribution of type _nType and numbers of classes _nClasses
.
_nType
is one of predefined types of solid distributions. All old data will be erased. Matrix will be initialized with zero values.
Returns false
on error.
bool SetDimensions(unsigned _nType1, unsigned _nClasses1, unsigned _nType2, unsigned _nClasses2)
Sets twodimensional distribution of types _nType1
and _nType2
and numbers of classes _nClasses1
and _nClasses2
.
_nType1
and _nType2
are from predefined types of solid distributions. All types must be unique. All old data will be erased. Matrix will be initialized with zero values.
Returns false
on error.
bool SetDimensions(unsigned _nType1, unsigned _nClasses1, unsigned _nType2, unsigned _nClasses2, unsigned _nType3, unsigned _nClasses3)
Sets threedimensional distribution of types _nType1
, nType2
and _nType3
and numbers of classes _nClasses1
, _nClasses2
and _nClasses3
.
_nType1
, nType2
and _nType3
are from predefined types of solid distributions. All types must be unique. All old data will be erased. Matrix will be initialized with zero values.
Returns false
on error.
bool SetDimensions(const std::vector<unsigned> &_vTypes, const std::vector<unsigned> &_vClasses)
Sets types _vTypes
of dimensions and numbers of classes _vClasses
.
_vTypes
is the vector of predefined types of solid distributions. All types must be unique. All old data will be erased. Matrix will be initialized with zero values.
Returns false
on error.
std::vector<unsigned> GetDimensions()
Returns vector with all currently defined dimensions types.
std::vector<unsigned> GetClasses()
Returns vector with current numbers of classes.
unsigned GetDimensionsNumber()
Returns current number of dimensions.
Get data
double GetValue(unsigned _nDim, unsigned _nCoord)
Returns value with specified coordinate _nCoord
of specified dimension _nDim
. It is possible to use this function if matrix has more than one dimension: if number of dimensions does not conform to the matrix, the sum of values by remaining dimensions will be returned.
Returns 1
on error.
double GetValue(unsigned _nDim1, unsigned _nCoord1, unsigned _nDim2, unsigned _nCoord2)
Returns value according to specified coordinates and dimensions. It is possible to use this function if matrix has more than two dimensions: if number of dimensions does not conform to the matrix, the sum of values by remaining dimensions will be returned. Sequence of dimensions may not match the sequence, which was defined in the matrix.
Returns 1
on error.
double GetValue(unsigned _nDim1, unsigned _nCoord1, unsigned _nDim2, unsigned _nCoord2, unsigned _nDim3, unsigned _nCoord3)
Returns value according to specified coordinates and dimensions. It is possible to use this function if matrix has more than three dimensions: if number of dimensions does not conform to the matrix, the sum of values by remaining dimensions will be returned. Sequence of dimensions may not match the sequence, which was defined in the matrix.
Returns 1
on error.
double GetValue(const std::vector<unsigned> &_vDims, const std::vector<unsigned> &_vCoords)
Returns value according to specified coordinates and dimensions. It is possible to use this function if matrix has more dimensions than was defined in _vDims
: if number of dimensions does not conform to the matrix, the sum of values by remaining dimensions will be returned. Sequence of dimensions may not match the sequence, which was defined in the matrix. Number of dimensions _vDims
and coordinates _vCoords
must be the same.
Returns 1
on error.
double GetValue(const std::vector<unsigned> &_vCoords)
Returns value by specified coordinates according to the full defined set of dimensions. Returns 1
on error.
bool GetVectorValue(unsigned _nDim, std::vector<double> &_vResult)
Returns vector _vResult
according to specified dimension. If number of dimensions in the matrix is more than one, then the sum of values by remaining dimensions will be returned.
Returns false
on error.
bool GetVectorValue(unsigned _nDim1, unsigned _nCoord1, unsigned _nDim2, std::vector<double> &_vResult)
Returns vector of values _vResult
according to specified dimensions and coordinate. If number of dimensions in the matrix is more than two, then the sum of values by remaining dimensions will be returned. Sequence of dimensions may not match the sequence, which was defined in the matrix.
Returns false
on error.
bool GetVectorValue(unsigned _nDim1, unsigned _nCoord1, unsigned _nDim2, unsigned _nCoord2, unsigned _nDim3, std::vector<double> &_vResult)
Returns vector of values _vResult
according to specified dimensions and coordinates. If number of dimensions in the matrix is more than three, then the sum of values by remaining dimensions will be returned. Sequence of dimensions may not match the sequence, which was defined in the matrix.
Returns false
on error.
bool GetVectorValue(const std::vector<unsigned> &_vDims, const std::vector<unsigned> &_vCoords, std::vector<double> &_vResult)
Returns vector of values _vResult
according to specified dimensions and coordinates. If number of dimensions in the matrix is more than it was specified in _vDims
, then the sum of values by remaining dimensions will be returned. Sequence of dimensions may not match the sequence, which was defined in the matrix. Number of coordinates _vCoords
must be one less than the number of dimensions _vDims
.
Returns false
on error.
bool GetVectorValue(const std::vector<unsigned> &_vCoords, std::vector<double> &_vResult)
Returns vector of values _vResult
by specified coordinates according to the full defined set of dimensions. Returns false
on error.
Set data
bool SetValue(unsigned _nCoord, double _dValue)
Sets value _dValue
with coordinate _nCoord
into onedimensional matrix. Sets the value only if the matrix has one dimension.
Returns false
on error.
bool SetValue(unsigned _nDim1, unsigned _nCoord1, unsigned _nCoord2, double _dValue)
Sets value _dValue
according to specified coordinates and dimensions into twodimensional matrix. Sets the value only if the matrix has two dimensions. Sequence of dimensions may not match the sequence, which was defined in the matrix.
Returns false on error
.
bool SetValue(unsigned _nDim1, unsigned _nCoord1, unsigned _nDim2, unsigned _nCoord2, unsigned _nCoord3, double _dValue)
Sets value _dValue
according to specified coordinates and dimensions into threedimensional matrix. Sets the value only if the matrix has three dimensions. Sequence of dimensions may not match the sequence, which was defined in the matrix.
Returns false
on error.
bool SetValue(const std::vector<unsigned> &_vDims, const std::vector<unsigned> &_vCoords, double _dValue)
Sets value _dValue
according to specified coordinates and dimensions. Sets the value only if the number of dimensions is the same as in the matrix. Number of dimensions _vDims
and coordinates _vCoords
must be the same. Sequence of dimensions may not match the sequence, which was defined in the matrix.
Returns false
on error.
bool SetValue(const std::vector<unsigned> &_vCoords, double _dValue)
Sets value _dValue
according to specified coordinates for full set of dimensions. Sets the value only if the number of coordinates is the same as the number of dimensions in the matrix. Number of coordinates _vCoords
must be equal to a number of dimensions in matrix.
Returns false
on error.
bool SetVectorValue(const std::vector<double> &_vValue)
Sets vector of values _vValue
in onedimensional matrix. Sets the values only if the matrix has one dimension.
Returns false
on error.
bool SetVectorValue(unsigned _nDim, unsigned _nCoord, const std::vector<double> &_vValue)
Sets vector of values _vValue
according to specified dimension and coordinate in twodimensional matrix. Sets the values only if the matrix has two dimensions. Sequence of dimensions may not match the sequence, which was defined in the matrix.
Returns false
on error.
bool SetVectorValue(unsigned _nDim1, unsigned _nCoord1, unsigned _nDim2, unsigned _nCoord2, const std::vector<double> &_vValue)
Sets vector of values _vValue
according to specified dimensions and coordinates in threedimensional matrix. Sets the values only if the matrix has three dimensions. Sequence of dimensions may not match the sequence, which was defined in the matrix.
Returns false
on error.
bool SetVectorValue(const std::vector<unsigned> &_vDims, const std::vector<unsigned> &_vCoords, const std::vector<double> &_vValue)
Sets vector of values _vValue
according to specified dimensions and coordinates. Sets values only if the number of dimensions in _vDims
is one less than in the matrix. Number of dimensions _vDims
and coordinates _vCoords
must be equal.
Returns false
on error.
bool SetVectorValue (const std::vector<unsigned> &_vCoords, const std::vector<double> &_vValue)
Sets vector of values _vValue
according to specified coordinates for full set of dimensions, which were defined in the matrix. Sets values only if the number of coordinates is one less than number of dimensions in the matrix.
Returns false
on error.
Overloaded operators
CDenseMDMatrix operator+(const CDenseMDMatrix &_matrix)
Performs addition of the matrix with the same dimensions. Returns an empty matrix with 0 defined dimensions in case of different dimenions.
CDenseMDMatrix operator(const CDenseMDMatrix &_matrix)
Performs subtraction of matrices with the same dimensions. Returns an empty matrix with 0 defined dimensions in case of different dimenions.
CDenseMDMatrix operator*(double _dFactor)
Performs multiplication of the matrix by a coefficient _dFactor
.
CDenseMDMatrix& operator=(const CDenseMDMatrix &_matrix)
Sets dimenions and data from the _matrix
to a left matrix.
Other functions
void Normalize()
Normalizes the matrix so that the sum of all elements equals to 1.
bool IsNormalized()
Returns true
if the matrix is normalized.
void ClearData()
Sets all elements in matrix equal to 0.
void Clear()
Removes all data and information about dimensions from the matrix.
List of universal constants
In Dyssol, some universal constants have been defined, which you may use directly in your codes.
Define 
Typical value 
Unit 


\(6.02214199(47) \cdot 10^{23}\) 
[1/mol] 

\(1.3806503(24) \cdot 10^{23}\) 
[J/K] 

\(101\,325\) 
[Pa] 

\(8.314472(15)\) 
[J/mol/K] 

\(2.99792458(1) \cdot 10^8\) 
[m/s] 

\(9.80665\) 
[m/s \(^2\)] 

\(298.15\) 
[K] 

\(101\,325\) 
[Pa] 

\(3.14159265358979323846\) 
[] 