mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-02-25 18:55:30 -06:00
rename "bbox" to "boundingBox"
This commit is contained in:
parent
3278f66be5
commit
29926a7d06
@ -577,7 +577,7 @@ private:
|
|||||||
{ return pos[0] < eps_; }
|
{ return pos[0] < eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onInlet_(const GlobalPosition &pos) const
|
bool onInlet_(const GlobalPosition &pos) const
|
||||||
{ return onRightBoundary_(pos) && (5 < pos[1]) && (pos[1] < 15); }
|
{ return onRightBoundary_(pos) && (5 < pos[1]) && (pos[1] < 15); }
|
||||||
|
@ -460,13 +460,13 @@ private:
|
|||||||
{ return pos[0] < eps_; }
|
{ return pos[0] < eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onLowerBoundary_(const GlobalPosition &pos) const
|
bool onLowerBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] < eps_; }
|
{ return pos[1] < eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &pos) const
|
bool onUpperBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] > this->bboxMax()[1] - eps_; }
|
{ return pos[1] > this->boundingBoxMax()[1] - eps_; }
|
||||||
|
|
||||||
bool isContaminated_(const GlobalPosition &pos) const
|
bool isContaminated_(const GlobalPosition &pos) const
|
||||||
{
|
{
|
||||||
|
@ -294,7 +294,7 @@ public:
|
|||||||
|
|
||||||
private:
|
private:
|
||||||
bool onLeftSide_(const GlobalPosition &pos) const
|
bool onLeftSide_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] < (this->bboxMin()[0] + this->bboxMax()[0]) / 2; }
|
{ return pos[0] < (this->boundingBoxMin()[0] + this->boundingBoxMax()[0]) / 2; }
|
||||||
|
|
||||||
void setupInitialFluidStates_()
|
void setupInitialFluidStates_()
|
||||||
{
|
{
|
||||||
|
@ -431,21 +431,21 @@ public:
|
|||||||
|
|
||||||
private:
|
private:
|
||||||
bool onLeftBoundary_(const GlobalPosition &pos) const
|
bool onLeftBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] < this->bboxMin()[0] + eps_; }
|
{ return pos[0] < this->boundingBoxMin()[0] + eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onLowerBoundary_(const GlobalPosition &pos) const
|
bool onLowerBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] < this->bboxMin()[1] + eps_; }
|
{ return pos[1] < this->boundingBoxMin()[1] + eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &pos) const
|
bool onUpperBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] > this->bboxMax()[1] - eps_; }
|
{ return pos[1] > this->boundingBoxMax()[1] - eps_; }
|
||||||
|
|
||||||
bool onInlet_(const GlobalPosition &pos) const
|
bool onInlet_(const GlobalPosition &pos) const
|
||||||
{
|
{
|
||||||
Scalar width = this->bboxMax()[0] - this->bboxMin()[0];
|
Scalar width = this->boundingBoxMax()[0] - this->boundingBoxMin()[0];
|
||||||
Scalar lambda = (this->bboxMax()[0] - pos[0]) / width;
|
Scalar lambda = (this->boundingBoxMax()[0] - pos[0]) / width;
|
||||||
|
|
||||||
if (!onUpperBoundary_(pos))
|
if (!onUpperBoundary_(pos))
|
||||||
return false;
|
return false;
|
||||||
|
@ -536,16 +536,16 @@ public:
|
|||||||
|
|
||||||
private:
|
private:
|
||||||
bool onLeftBoundary_(const GlobalPosition &pos) const
|
bool onLeftBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] < this->bboxMin()[0] + eps_; }
|
{ return pos[0] < this->boundingBoxMin()[0] + eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onLowerBoundary_(const GlobalPosition &pos) const
|
bool onLowerBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] < this->bboxMin()[1] + eps_; }
|
{ return pos[1] < this->boundingBoxMin()[1] + eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &pos) const
|
bool onUpperBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] > this->bboxMax()[1] - eps_; }
|
{ return pos[1] > this->boundingBoxMax()[1] - eps_; }
|
||||||
|
|
||||||
void computeHeatCondParams_(HeatConductionLawParams ¶ms, Scalar poro)
|
void computeHeatCondParams_(HeatConductionLawParams ¶ms, Scalar poro)
|
||||||
{
|
{
|
||||||
|
@ -334,7 +334,7 @@ private:
|
|||||||
{ return pos[dim - 1] < eps_; }
|
{ return pos[dim - 1] < eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &pos) const
|
bool onUpperBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[dim - 1] > this->bboxMax()[dim - 1] - eps_; }
|
{ return pos[dim - 1] > this->boundingBoxMax()[dim - 1] - eps_; }
|
||||||
|
|
||||||
bool isInLens_(const GlobalPosition &pos) const
|
bool isInLens_(const GlobalPosition &pos) const
|
||||||
{
|
{
|
||||||
|
@ -387,13 +387,13 @@ private:
|
|||||||
{ return pos[0] < eps_; }
|
{ return pos[0] < eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onLowerBoundary_(const GlobalPosition &pos) const
|
bool onLowerBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] < eps_; }
|
{ return pos[1] < eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &pos) const
|
bool onUpperBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] > this->bboxMax()[1] - eps_; }
|
{ return pos[1] > this->boundingBoxMax()[1] - eps_; }
|
||||||
|
|
||||||
bool onInlet_(const GlobalPosition &pos) const
|
bool onInlet_(const GlobalPosition &pos) const
|
||||||
{ return onUpperBoundary_(pos) && 50 < pos[0] && pos[0] < 75; }
|
{ return onUpperBoundary_(pos) && 50 < pos[0] && pos[0] < 75; }
|
||||||
|
@ -409,8 +409,8 @@ public:
|
|||||||
|
|
||||||
// set wetting phase pressure and saturation
|
// set wetting phase pressure and saturation
|
||||||
if (onLeftBoundary_(pos)) {
|
if (onLeftBoundary_(pos)) {
|
||||||
Scalar height = this->bboxMax()[1] - this->bboxMin()[1];
|
Scalar height = this->boundingBoxMax()[1] - this->boundingBoxMin()[1];
|
||||||
Scalar depth = this->bboxMax()[1] - pos[1];
|
Scalar depth = this->boundingBoxMax()[1] - pos[1];
|
||||||
Scalar alpha = (1 + 1.5 / height);
|
Scalar alpha = (1 + 1.5 / height);
|
||||||
|
|
||||||
// hydrostatic pressure scaled by alpha
|
// hydrostatic pressure scaled by alpha
|
||||||
@ -418,7 +418,7 @@ public:
|
|||||||
Sw = 1.0;
|
Sw = 1.0;
|
||||||
}
|
}
|
||||||
else {
|
else {
|
||||||
Scalar depth = this->bboxMax()[1] - pos[1];
|
Scalar depth = this->boundingBoxMax()[1] - pos[1];
|
||||||
|
|
||||||
// hydrostatic pressure
|
// hydrostatic pressure
|
||||||
pw = 1e5 - densityW * this->gravity()[1] * depth;
|
pw = 1e5 - densityW * this->gravity()[1] * depth;
|
||||||
@ -472,7 +472,7 @@ public:
|
|||||||
int timeIdx) const
|
int timeIdx) const
|
||||||
{
|
{
|
||||||
const GlobalPosition &pos = context.pos(spaceIdx, timeIdx);
|
const GlobalPosition &pos = context.pos(spaceIdx, timeIdx);
|
||||||
Scalar depth = this->bboxMax()[1] - pos[1];
|
Scalar depth = this->boundingBoxMax()[1] - pos[1];
|
||||||
|
|
||||||
Opm::ImmiscibleFluidState<Scalar, FluidSystem> fs;
|
Opm::ImmiscibleFluidState<Scalar, FluidSystem> fs;
|
||||||
fs.setPressure(wPhaseIdx, /*pressure=*/1e5);
|
fs.setPressure(wPhaseIdx, /*pressure=*/1e5);
|
||||||
@ -529,21 +529,21 @@ private:
|
|||||||
}
|
}
|
||||||
|
|
||||||
bool onLeftBoundary_(const GlobalPosition &pos) const
|
bool onLeftBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] < this->bboxMin()[0] + eps_; }
|
{ return pos[0] < this->boundingBoxMin()[0] + eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onLowerBoundary_(const GlobalPosition &pos) const
|
bool onLowerBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] < this->bboxMin()[1] + eps_; }
|
{ return pos[1] < this->boundingBoxMin()[1] + eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &pos) const
|
bool onUpperBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] > this->bboxMax()[1] - eps_; }
|
{ return pos[1] > this->boundingBoxMax()[1] - eps_; }
|
||||||
|
|
||||||
bool onInlet_(const GlobalPosition &pos) const
|
bool onInlet_(const GlobalPosition &pos) const
|
||||||
{
|
{
|
||||||
Scalar width = this->bboxMax()[0] - this->bboxMin()[0];
|
Scalar width = this->boundingBoxMax()[0] - this->boundingBoxMin()[0];
|
||||||
Scalar lambda = (this->bboxMax()[0] - pos[0]) / width;
|
Scalar lambda = (this->boundingBoxMax()[0] - pos[0]) / width;
|
||||||
return onUpperBoundary_(pos) && 0.5 < lambda && lambda < 2.0 / 3.0;
|
return onUpperBoundary_(pos) && 0.5 < lambda && lambda < 2.0 / 3.0;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
@ -244,16 +244,16 @@ private:
|
|||||||
}
|
}
|
||||||
|
|
||||||
bool onLeftBoundary_(const GlobalPosition &globalPos) const
|
bool onLeftBoundary_(const GlobalPosition &globalPos) const
|
||||||
{ return globalPos[0] < this->bboxMin()[0] + eps_; }
|
{ return globalPos[0] < this->boundingBoxMin()[0] + eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &globalPos) const
|
bool onRightBoundary_(const GlobalPosition &globalPos) const
|
||||||
{ return globalPos[0] > this->bboxMax()[0] - eps_; }
|
{ return globalPos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onLowerBoundary_(const GlobalPosition &globalPos) const
|
bool onLowerBoundary_(const GlobalPosition &globalPos) const
|
||||||
{ return globalPos[1] < this->bboxMin()[1] + eps_; }
|
{ return globalPos[1] < this->boundingBoxMin()[1] + eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &globalPos) const
|
bool onUpperBoundary_(const GlobalPosition &globalPos) const
|
||||||
{ return globalPos[1] > this->bboxMax()[1] - eps_; }
|
{ return globalPos[1] > this->boundingBoxMax()[1] - eps_; }
|
||||||
|
|
||||||
Scalar eps_;
|
Scalar eps_;
|
||||||
};
|
};
|
||||||
|
@ -291,7 +291,7 @@ private:
|
|||||||
{ return pos[0] < eps_; }
|
{ return pos[0] < eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
template <class FluidState, class Context>
|
template <class FluidState, class Context>
|
||||||
void initialFluidState_(FluidState &fs, const Context &context,
|
void initialFluidState_(FluidState &fs, const Context &context,
|
||||||
@ -300,9 +300,9 @@ private:
|
|||||||
Scalar T = temperature(context, spaceIdx, timeIdx);
|
Scalar T = temperature(context, spaceIdx, timeIdx);
|
||||||
// Scalar rho = FluidSystem::H2O::liquidDensity(T, /*pressure=*/1.5e5);
|
// Scalar rho = FluidSystem::H2O::liquidDensity(T, /*pressure=*/1.5e5);
|
||||||
// Scalar z = context.pos(spaceIdx, timeIdx)[dim - 1] -
|
// Scalar z = context.pos(spaceIdx, timeIdx)[dim - 1] -
|
||||||
// this->bboxMax()[dim - 1];
|
// this->boundingBoxMax()[dim - 1];
|
||||||
// Scalar z = context.pos(spaceIdx, timeIdx)[dim - 1] -
|
// Scalar z = context.pos(spaceIdx, timeIdx)[dim - 1] -
|
||||||
// this->bboxMax()[dim - 1];
|
// this->boundingBoxMax()[dim - 1];
|
||||||
|
|
||||||
fs.setSaturation(/*phaseIdx=*/0, 1.0);
|
fs.setSaturation(/*phaseIdx=*/0, 1.0);
|
||||||
fs.setPressure(/*phaseIdx=*/0, 1e5 /* + rho*z */);
|
fs.setPressure(/*phaseIdx=*/0, 1e5 /* + rho*z */);
|
||||||
|
@ -361,10 +361,10 @@ public:
|
|||||||
|
|
||||||
private:
|
private:
|
||||||
bool onLeftBoundary_(const GlobalPosition &pos) const
|
bool onLeftBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] < this->bboxMin()[0] + eps_; }
|
{ return pos[0] < this->boundingBoxMin()[0] + eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
void setupInitialFluidState_()
|
void setupInitialFluidState_()
|
||||||
{
|
{
|
||||||
|
@ -424,10 +424,10 @@ public:
|
|||||||
int spaceIdx, int timeIdx) const
|
int spaceIdx, int timeIdx) const
|
||||||
{
|
{
|
||||||
const auto &pos = context.pos(spaceIdx, timeIdx);
|
const auto &pos = context.pos(spaceIdx, timeIdx);
|
||||||
Scalar x = pos[0] - this->bboxMin()[0];
|
Scalar x = pos[0] - this->boundingBoxMin()[0];
|
||||||
Scalar y = pos[dim - 1] - this->bboxMin()[dim - 1];
|
Scalar y = pos[dim - 1] - this->boundingBoxMin()[dim - 1];
|
||||||
Scalar height = this->bboxMax()[dim - 1] - this->bboxMin()[dim - 1];
|
Scalar height = this->boundingBoxMax()[dim - 1] - this->boundingBoxMin()[dim - 1];
|
||||||
Scalar width = this->bboxMax()[0] - this->bboxMin()[0];
|
Scalar width = this->boundingBoxMax()[0] - this->boundingBoxMin()[0];
|
||||||
if ((onLeftBoundary_(pos) || onRightBoundary_(pos)) && y < height / 2) {
|
if ((onLeftBoundary_(pos) || onRightBoundary_(pos)) && y < height / 2) {
|
||||||
// injectors
|
// injectors
|
||||||
auto fs = initialFluidState_;
|
auto fs = initialFluidState_;
|
||||||
@ -566,7 +566,7 @@ private:
|
|||||||
{ return pos[0] < eps_; }
|
{ return pos[0] < eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onInlet_(const GlobalPosition &pos) const
|
bool onInlet_(const GlobalPosition &pos) const
|
||||||
{ return onRightBoundary_(pos) && (5 < pos[1]) && (pos[1] < 15); }
|
{ return onRightBoundary_(pos) && (5 < pos[1]) && (pos[1] < 15); }
|
||||||
|
@ -371,21 +371,21 @@ public:
|
|||||||
|
|
||||||
private:
|
private:
|
||||||
bool onLeftBoundary_(const GlobalPosition &pos) const
|
bool onLeftBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] < this->bboxMin()[0] + eps_; }
|
{ return pos[0] < this->boundingBoxMin()[0] + eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onLowerBoundary_(const GlobalPosition &pos) const
|
bool onLowerBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] < this->bboxMin()[1] + eps_; }
|
{ return pos[1] < this->boundingBoxMin()[1] + eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &pos) const
|
bool onUpperBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] > this->bboxMax()[1] - eps_; }
|
{ return pos[1] > this->boundingBoxMax()[1] - eps_; }
|
||||||
|
|
||||||
bool onInlet_(const GlobalPosition &pos) const
|
bool onInlet_(const GlobalPosition &pos) const
|
||||||
{
|
{
|
||||||
Scalar width = this->bboxMax()[0] - this->bboxMin()[0];
|
Scalar width = this->boundingBoxMax()[0] - this->boundingBoxMin()[0];
|
||||||
Scalar lambda = (this->bboxMax()[0] - pos[0]) / width;
|
Scalar lambda = (this->boundingBoxMax()[0] - pos[0]) / width;
|
||||||
return onUpperBoundary_(pos) && 0.5 < lambda && lambda < 2.0 / 3.0;
|
return onUpperBoundary_(pos) && 0.5 < lambda && lambda < 2.0 / 3.0;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
@ -214,10 +214,10 @@ public:
|
|||||||
// parabolic profile
|
// parabolic profile
|
||||||
const Scalar v1 = 1.0;
|
const Scalar v1 = 1.0;
|
||||||
values[velocity0Idx + 1]
|
values[velocity0Idx + 1]
|
||||||
= -v1 * (globalPos[0] - this->bboxMin()[0])
|
= -v1 * (globalPos[0] - this->boundingBoxMin()[0])
|
||||||
* (this->bboxMax()[0] - globalPos[0])
|
* (this->boundingBoxMax()[0] - globalPos[0])
|
||||||
/ (0.25 * (this->bboxMax()[0] - this->bboxMin()[0])
|
/ (0.25 * (this->boundingBoxMax()[0] - this->boundingBoxMin()[0])
|
||||||
* (this->bboxMax()[0] - this->bboxMin()[0]));
|
* (this->boundingBoxMax()[0] - this->boundingBoxMin()[0]));
|
||||||
|
|
||||||
Scalar moleFrac[numComponents];
|
Scalar moleFrac[numComponents];
|
||||||
if (onUpperBoundary_(globalPos))
|
if (onUpperBoundary_(globalPos))
|
||||||
@ -273,16 +273,16 @@ public:
|
|||||||
|
|
||||||
private:
|
private:
|
||||||
bool onLeftBoundary_(const GlobalPosition &globalPos) const
|
bool onLeftBoundary_(const GlobalPosition &globalPos) const
|
||||||
{ return globalPos[0] < this->bboxMin()[0] + eps_; }
|
{ return globalPos[0] < this->boundingBoxMin()[0] + eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &globalPos) const
|
bool onRightBoundary_(const GlobalPosition &globalPos) const
|
||||||
{ return globalPos[0] > this->bboxMax()[0] - eps_; }
|
{ return globalPos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onLowerBoundary_(const GlobalPosition &globalPos) const
|
bool onLowerBoundary_(const GlobalPosition &globalPos) const
|
||||||
{ return globalPos[1] < this->bboxMin()[1] + eps_; }
|
{ return globalPos[1] < this->boundingBoxMin()[1] + eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &globalPos) const
|
bool onUpperBoundary_(const GlobalPosition &globalPos) const
|
||||||
{ return globalPos[1] > this->bboxMax()[1] - eps_; }
|
{ return globalPos[1] > this->boundingBoxMax()[1] - eps_; }
|
||||||
|
|
||||||
Scalar eps_;
|
Scalar eps_;
|
||||||
};
|
};
|
||||||
|
@ -213,9 +213,9 @@ public:
|
|||||||
moleFrac[AirIdx] = 1 - moleFrac[H2OIdx];
|
moleFrac[AirIdx] = 1 - moleFrac[H2OIdx];
|
||||||
|
|
||||||
// parabolic velocity profile
|
// parabolic velocity profile
|
||||||
Scalar y = this->bboxMax()[1] - pos[1];
|
Scalar y = this->boundingBoxMax()[1] - pos[1];
|
||||||
Scalar x = pos[0] - this->bboxMin()[0];
|
Scalar x = pos[0] - this->boundingBoxMin()[0];
|
||||||
Scalar width = this->bboxMax()[0] - this->bboxMin()[0];
|
Scalar width = this->boundingBoxMax()[0] - this->boundingBoxMin()[0];
|
||||||
|
|
||||||
// parabolic velocity profile
|
// parabolic velocity profile
|
||||||
const Scalar maxVelocity = 1.0;
|
const Scalar maxVelocity = 1.0;
|
||||||
@ -285,16 +285,16 @@ public:
|
|||||||
|
|
||||||
private:
|
private:
|
||||||
bool onLeftBoundary_(const GlobalPosition &pos) const
|
bool onLeftBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] < this->bboxMin()[0] + eps_; }
|
{ return pos[0] < this->boundingBoxMin()[0] + eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onLowerBoundary_(const GlobalPosition &pos) const
|
bool onLowerBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] < this->bboxMin()[1] + eps_; }
|
{ return pos[1] < this->boundingBoxMin()[1] + eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &pos) const
|
bool onUpperBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] > this->bboxMax()[1] - eps_; }
|
{ return pos[1] > this->boundingBoxMax()[1] - eps_; }
|
||||||
|
|
||||||
bool onBoundary_(const GlobalPosition &pos) const
|
bool onBoundary_(const GlobalPosition &pos) const
|
||||||
{
|
{
|
||||||
|
@ -180,8 +180,8 @@ public:
|
|||||||
{
|
{
|
||||||
const GlobalPosition &pos = context.pos(spaceIdx, timeIdx);
|
const GlobalPosition &pos = context.pos(spaceIdx, timeIdx);
|
||||||
|
|
||||||
Scalar y = pos[1] - this->bboxMin()[1];
|
Scalar y = pos[1] - this->boundingBoxMin()[1];
|
||||||
Scalar height = this->bboxMax()[1] - this->bboxMin()[1];
|
Scalar height = this->boundingBoxMax()[1] - this->boundingBoxMin()[1];
|
||||||
|
|
||||||
// parabolic velocity profile
|
// parabolic velocity profile
|
||||||
const Scalar maxVelocity = 1.0;
|
const Scalar maxVelocity = 1.0;
|
||||||
@ -221,8 +221,8 @@ public:
|
|||||||
{
|
{
|
||||||
const auto &pos = context.pos(spaceIdx, timeIdx);
|
const auto &pos = context.pos(spaceIdx, timeIdx);
|
||||||
|
|
||||||
Scalar y = pos[1] - this->bboxMin()[1];
|
Scalar y = pos[1] - this->boundingBoxMin()[1];
|
||||||
Scalar height = this->bboxMax()[1] - this->bboxMin()[1];
|
Scalar height = this->boundingBoxMax()[1] - this->boundingBoxMin()[1];
|
||||||
|
|
||||||
// parabolic velocity profile on boundaries
|
// parabolic velocity profile on boundaries
|
||||||
const Scalar maxVelocity = 1.0;
|
const Scalar maxVelocity = 1.0;
|
||||||
@ -280,16 +280,16 @@ public:
|
|||||||
|
|
||||||
private:
|
private:
|
||||||
bool onLeftBoundary_(const GlobalPosition &pos) const
|
bool onLeftBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] < this->bboxMin()[0] + eps_; }
|
{ return pos[0] < this->boundingBoxMin()[0] + eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onLowerBoundary_(const GlobalPosition &pos) const
|
bool onLowerBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] < this->bboxMin()[1] + eps_; }
|
{ return pos[1] < this->boundingBoxMin()[1] + eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &pos) const
|
bool onUpperBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] > this->bboxMax()[1] - eps_; }
|
{ return pos[1] > this->boundingBoxMax()[1] - eps_; }
|
||||||
|
|
||||||
bool onBoundary_(const GlobalPosition &pos) const
|
bool onBoundary_(const GlobalPosition &pos) const
|
||||||
{
|
{
|
||||||
|
@ -444,13 +444,13 @@ private:
|
|||||||
{ return pos[0] < eps_; }
|
{ return pos[0] < eps_; }
|
||||||
|
|
||||||
bool onRightBoundary_(const GlobalPosition &pos) const
|
bool onRightBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[0] > this->bboxMax()[0] - eps_; }
|
{ return pos[0] > this->boundingBoxMax()[0] - eps_; }
|
||||||
|
|
||||||
bool onLowerBoundary_(const GlobalPosition &pos) const
|
bool onLowerBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] < eps_; }
|
{ return pos[1] < eps_; }
|
||||||
|
|
||||||
bool onUpperBoundary_(const GlobalPosition &pos) const
|
bool onUpperBoundary_(const GlobalPosition &pos) const
|
||||||
{ return pos[1] > this->bboxMax()[1] - eps_; }
|
{ return pos[1] > this->boundingBoxMax()[1] - eps_; }
|
||||||
|
|
||||||
bool onInlet_(const GlobalPosition &pos) const
|
bool onInlet_(const GlobalPosition &pos) const
|
||||||
{ return onLowerBoundary_(pos) && (15.0 < pos[0]) && (pos[0] < 25.0); }
|
{ return onLowerBoundary_(pos) && (15.0 < pos[0]) && (pos[0] < 25.0); }
|
||||||
|
@ -244,7 +244,7 @@ public:
|
|||||||
|
|
||||||
values.setFreeFlow(context, spaceIdx, timeIdx, fs);
|
values.setFreeFlow(context, spaceIdx, timeIdx, fs);
|
||||||
}
|
}
|
||||||
else if (pos[0] > this->bboxMax()[0] - eps_) {
|
else if (pos[0] > this->boundingBoxMax()[0] - eps_) {
|
||||||
// forced outflow at the right boundary
|
// forced outflow at the right boundary
|
||||||
RateVector massRate(0.0);
|
RateVector massRate(0.0);
|
||||||
|
|
||||||
|
Loading…
Reference in New Issue
Block a user