1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
|
/* Copyright (c) 2020, Dyssol Development Team. All rights reserved. This file is part of Dyssol. See LICENSE file for license information. */
#define _USE_MATH_DEFINES
#include "Holdup.h"
#include "Stream.h"
#include "TimeDependentValue.h"
#include "Phase.h"
#include "ContainerFunctions.h"
#include <cmath>
CHoldup::CHoldup(const std::string& _key) :
CBaseStream{ _key }
{
m_overall[EOverall::OVERALL_MASS]->SetName("Mass");
m_overall[EOverall::OVERALL_MASS]->SetUnits("kg");
}
CHoldup::CHoldup(const std::string& _key, const CMaterialsDatabase* _materialsDB, const CMultidimensionalGrid& _grid,
const std::vector<SOverallDescriptor>* _overall, const std::vector<SPhaseDescriptor>* _phases,
const SCacheSettings* _cache, const SToleranceSettings* _tolerance, const SThermodynamicsSettings* _thermodynamics) :
CBaseStream{ _key, _materialsDB, _grid, _overall, _phases, _cache, _tolerance, _thermodynamics }
{
m_overall[EOverall::OVERALL_MASS]->SetName("Mass");
m_overall[EOverall::OVERALL_MASS]->SetUnits("kg");
}
CHoldup::CHoldup(const CBaseStream& _other) :
CBaseStream{ _other }
{
m_overall[EOverall::OVERALL_MASS]->SetName("Mass");
m_overall[EOverall::OVERALL_MASS]->SetUnits("kg");
}
void CHoldup::CopyFromHoldup(double _time, const CHoldup* _source)
{
Copy(_time, *_source);
}
void CHoldup::CopyFromHoldup(double _timeBeg, double _timeEnd, const CHoldup* _source)
{
Copy(_timeBeg, _timeEnd, *_source);
}
void CHoldup::CopyFromHoldup(double _timeDst, const CHoldup* _source, double _timeSrc)
{
Copy(_timeDst, *_source, _timeSrc);
}
/* All time points after _timeEnd are removed, since after mixing they make no sense anymore.
Nevertheless, _timeEnd should stay to enable successive adding of several streams to the holdup.
If holdup has no time points after _timeBeg, adds the stream to the values in _timeBeg.
If there are other time points after _timeBeg, adds the stream to the values in _timeEnd (possibly extrapolated from the last defined time point)*/
void CHoldup::AddStream(double _timeBeg, double _timeEnd, const CStream* _source)
{
if (_timeBeg > _timeEnd) return;
if (!HaveSameStructure(*this, *_source)) return;
// remove discarded time points
RemoveTimePointsAfter(_timeEnd);
// special handling if there are no time points in the holdup
if (m_timePoints.empty())
{
Copy(_timeBeg, *_source);
SetMass(_timeBeg, 0.0);
}
// ensure the stream has only these two time points for proper interpolation
const CStream* stream = _source;
if (std::fabs(_source->GetPreviousTimePoint(_timeEnd) - _timeBeg) > m_epsilon) // there are some points in-between - create a copy with only two points
{
auto* temp = new CStream();
temp->SetupStructure(_source);
temp->SetCacheSettings(SCacheSettings{ false, 0, {} });
temp->Copy(_timeBeg, *_source);
temp->Copy(_timeEnd, *_source);
stream = temp;
}
// calculate time point in the stream
const double timeDelta = (_timeEnd - _timeBeg) / M_SQRT2;
const bool timeGrowth = stream->GetMassFlow(_timeBeg) < stream->GetMassFlow(_timeEnd);
const double timeSrc = timeGrowth ? _timeBeg + timeDelta : _timeEnd - timeDelta;
// get masses
const double massSrc = stream->GetMassFlow((_timeBeg + _timeEnd) / 2) * (_timeEnd - _timeBeg);
const double massDst = GetLastTimePoint() <= _timeBeg ? GetMass(_timeBeg) : GetMass(_timeEnd);
// calculate mixture
const mix_type mix = CalculateMix(timeSrc, *stream, massSrc, _timeEnd, *this, massDst);
// set obtained mixture to the temporary holdup
SetMix(_timeEnd, mix);
// remove time points between begin and end, since they are not consistent anymore
RemoveTimePoints(_timeBeg, _timeEnd, false);
// clean up temporary stream if necessary
if (stream != _source)
delete stream;
}
void CHoldup::AddHoldup(double _time, const CHoldup* _source)
{
Add(_time, *_source);
}
void CHoldup::AddHoldup(double _timeBeg, double _timeEnd, const CHoldup* _source)
{
Add(_timeBeg, _timeEnd, *_source);
}
// TODO: move it somewhere
////////////////////////////////////////////////////////////////////////////////
/// Deprecated functions
void CHoldup::CopyFromHoldup(const CHoldup* _source, double _time)
{
CopyFromHoldup(_time, _source);
}
void CHoldup::AddHoldup(const CHoldup* _source, double _time)
{
AddHoldup(_time, _source);
}
void CHoldup::AddStream(const CStream* _source, double _timeBeg, double _timeEnd)
{
AddStream(_timeBeg, _timeEnd, _source);
}
void CHoldup::AddStream2(const CStream* _source, double _timeBeg, double _timeEnd)
{
AddStream(_timeBeg, _timeEnd, _source);
}
double CHoldup::GetMass(double _time, unsigned _basis) const
{
if (_basis == 0)
return GetMass(_time);
else
return GetMol(_time);
}
void CHoldup::SetMass(double _time, double _value, unsigned _basis)
{
if (_basis == 0)
SetMass(_time, _value);
else
SetMol(_time, _value);
}
double CHoldup::GetCompoundMass(double _time, const std::string& _compoundKey, unsigned _soa, unsigned _basis) const
{
if (_basis == 0)
return GetCompoundMass(_time, _compoundKey, SOA2EPhase(_soa));
else
return GetCompoundMol(_time, _compoundKey, SOA2EPhase(_soa));
}
double CHoldup::GetPhaseMass(double _time, unsigned _soa, unsigned _basis) const
{
if (_basis == 0)
return GetPhaseMass(_time, SOA2EPhase(_soa));
else
return GetPhaseMol(_time, SOA2EPhase(_soa));
}
void CHoldup::SetPhaseMass(double _time, unsigned _soa, double _value, unsigned _basis)
{
if (_basis == 0)
SetPhaseMass(_time, SOA2EPhase(_soa), _value);
else
SetPhaseMol(_time, SOA2EPhase(_soa), _value);
}
|
