From 701205161723b5da0eae2d590cc6b817c96b91f3 Mon Sep 17 00:00:00 2001
From: Thomas Fetzer <thomas.fetzer@iws.uni-stuttgart.de>
Date: Mon, 30 Jun 2014 13:48:27 +0000
Subject: [PATCH] [multidomain] replaced hourly and diurnal variation by one
general function, which can be adjusted by the paramFile
reviewed by gruenich
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@13001 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
.../2cnistokes2p2cni/stokes2cnisubproblem.hh | 74 ++++++++---------
.../test_2cnistokes2p2cni.input | 17 ++--
.../test_2cnistokes2p2cni_reference.input | 17 ++--
.../2cstokes2p2c/stokes2csubproblem.hh | 80 +++++++++++--------
.../2cstokes2p2c/test_2cstokes2p2c.input | 13 ++-
.../test_2cstokes2p2c_reference.input | 13 ++-
6 files changed, 123 insertions(+), 91 deletions(-)
diff --git a/test/multidomain/2cnistokes2p2cni/stokes2cnisubproblem.hh b/test/multidomain/2cnistokes2p2cni/stokes2cnisubproblem.hh
index e53eebf7a8..cbc66c5dd7 100644
--- a/test/multidomain/2cnistokes2p2cni/stokes2cnisubproblem.hh
+++ b/test/multidomain/2cnistokes2p2cni/stokes2cnisubproblem.hh
@@ -199,22 +199,26 @@ public:
bboxMin_[1] = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Grid, InterfacePos);
bboxMax_[1] = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Grid, YMax);
- refTemperature_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, RefTemperature);
+ refVelocity_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, RefVelocity);
refPressure_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, RefPressure);
refMassfrac_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, RefMassfrac);
- vxMax_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, VxMax);
- bjSlipVel_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, BeaversJosephSlipVel);
+ refTemperature_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, RefTemperature);
- sinusVelVar_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusVelVar);
- sinusPVar_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusPVar);
- sinusTVar_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusTVar);
- sinusXVar_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusXVar);
+ sinusVAmplitude_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusVelAmplitude);
+ sinusVPeriod_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusVelPeriod);
+ sinusPAmplitude_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusPressureAmplitude);
+ sinusPPeriod_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusPressurePeriod);
+ sinusXAmplitude_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusConcentrationAmplitude);
+ sinusXPeriod_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusConcentrationPeriod);
+ sinusTAmplitude_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusTemperatureAmplitude);
+ sinusTPeriod_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusTemperaturePeriod);
+ bjSlipVel_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, BeaversJosephSlipVel);
+ alphaBJ_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, SpatialParams, AlphaBJ);
xMaterialInterface_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, SpatialParams, MaterialInterfaceX);
+
runUpDistanceX_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Problem, RunUpDistanceX); // first part of the interface without coupling
initializationTime_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, TimeManager, InitTime);
-
- alphaBJ_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, SpatialParams, AlphaBJ);
}
catch (Dumux::ParameterException &e) {
std::cerr << e << ". Abort!\n";
@@ -512,17 +516,21 @@ public:
const SpatialParams &spatialParams() const
{ return spatialParams_; }
+ //! \brief Returns the reference velocity.
+ const Scalar refVelocity() const
+ { return refVelocity_ + variation_(sinusVAmplitude_, sinusVPeriod_); }
+
//! \brief Returns the reference pressure.
const Scalar refPressure() const
- { return refPressure_ + diurnalVariation_(sinusPVar_); }
-
- //! \brief Returns the reference temperature.
- const Scalar refTemperature() const
- { return refTemperature_+ diurnalVariation_(sinusTVar_); }
+ { return refPressure_ + variation_(sinusPAmplitude_, sinusPPeriod_); }
//! \brief Returns the reference mass fraction.
const Scalar refMassfrac() const
- { return refMassfrac_ + diurnalVariation_(sinusXVar_); }
+ { return refMassfrac_ + variation_(sinusXAmplitude_, sinusXPeriod_); }
+
+ //! \brief Returns the reference temperature.
+ const Scalar refTemperature() const
+ { return refTemperature_+ variation_(sinusTAmplitude_, sinusTPeriod_); }
private:
/*!
@@ -550,7 +558,7 @@ private:
//! \brief set the profile of the inflow velocity (horizontal direction)
const Scalar xVelocity_(const GlobalPosition &globalPos) const
{
- const Scalar vmax = vxMax_ + hourlyVariation_(sinusVelVar_);
+ const Scalar vmax = refVelocity();
// const Scalar relativeHeight = (globalPos[1]-bboxMin_[1])/height_();
// linear profile
// return vmax*relativeHeight + bjSlipVel_; // BJ slip velocity is added as sqrt(Kxx)
@@ -583,19 +591,9 @@ private:
fluidState.setMoleFraction(phaseIdx, phaseCompIdx, massFraction[phaseCompIdx]*M1/massToMoleDenominator);
}
-
- const Scalar diurnalVariation_(const Scalar value) const
- {
- const Scalar time = this->timeManager().time() + this->timeManager().timeStepSize();
- return sin(2*M_PI*time/86400) * value;
- }
-
-
- const Scalar hourlyVariation_(const Scalar value) const
- {
- const Scalar time = this->timeManager().time();
- return sin(2*M_PI*time/3600) * value;
- }
+ // can be used for the variation of a boundary condition
+ const Scalar variation_(const Scalar amplitude, const Scalar period) const
+ { return sin(2*M_PI*this->timeManager().time()/period) * amplitude; }
bool onLeftBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[0] < bboxMin_[0] + eps_; }
@@ -626,19 +624,23 @@ private:
GlobalPosition bboxMin_;
GlobalPosition bboxMax_;
+ Scalar refVelocity_;
Scalar refPressure_;
- Scalar refTemperature_;
Scalar refMassfrac_;
+ Scalar refTemperature_;
+
+ Scalar sinusVAmplitude_;
+ Scalar sinusVPeriod_;
+ Scalar sinusPAmplitude_;
+ Scalar sinusPPeriod_;
+ Scalar sinusXAmplitude_;
+ Scalar sinusXPeriod_;
+ Scalar sinusTAmplitude_;
+ Scalar sinusTPeriod_;
- Scalar vxMax_;
Scalar bjSlipVel_;
Scalar alphaBJ_;
- Scalar sinusVelVar_;
- Scalar sinusPVar_;
- Scalar sinusTVar_;
- Scalar sinusXVar_;
-
Scalar xMaterialInterface_;
Scalar runUpDistanceX_;
Scalar initializationTime_;
diff --git a/test/multidomain/2cnistokes2p2cni/test_2cnistokes2p2cni.input b/test/multidomain/2cnistokes2p2cni/test_2cnistokes2p2cni.input
index e23bcbb1a2..8a994c82d9 100644
--- a/test/multidomain/2cnistokes2p2cni/test_2cnistokes2p2cni.input
+++ b/test/multidomain/2cnistokes2p2cni/test_2cnistokes2p2cni.input
@@ -78,10 +78,18 @@ StabilizationAlpha = -1.0
#############################################################
[FreeFlow]
#############################################################
+RefVelocity = 3.5
RefPressure = 1e5
+RefMassfrac = 0.008 # 0.01
RefTemperature = 298.15
-RefMassfrac = 0.008
-VxMax = 3.5
+SinusVelAmplitude = 0.0 # 0.1
+SinusVelPeriod = 3600 # in seconds
+SinusPressureAmplitude = 0.0 # 0.5
+SinusPressurePeriod = 3600 # in seconds
+SinusConcentrationAmplitude = 0.0 # 0.5
+SinusConcentrationPeriod = 3600 # in seconds
+SinusTemperatureAmplitude = 0.0 # 0.5
+SinusTemperaturePeriod = 3600 # in seconds
BeaversJosephSlipVel = 0.00134
ExponentMTC = 0.0 # 1./6., Mass transfer coefficient for S^MTC
UseBoundaryLayerModel = 0 # 0 for no boundary layer model, 1 for Blasius, 2 and 3 for turbulent BL, 9 for constant thickness
@@ -89,11 +97,6 @@ BoundaryLayerOffset = 0.0 # for BL model like Blasius, determines a virtual ru
ConstThickness = 0.0016 # for a constant BL thickness, use BL model 9
MassTransferModel = 0 # 0 for none, 1 for power law, 2 for Schluender model
-SinusVelVar = 0.0 # hourly velocity variation
-SinusPVar = 0.0 # diurnal pressure variation
-SinusTVar = 0.0 # diurnal temperature variation
-SinusXVar = 0.0 # diurnal variation of the vapor concentration (massfraction)
-
#############################################################
[PorousMedium]
#############################################################
diff --git a/test/multidomain/2cnistokes2p2cni/test_references/test_2cnistokes2p2cni_reference.input b/test/multidomain/2cnistokes2p2cni/test_references/test_2cnistokes2p2cni_reference.input
index 0f7befcc6f..7d1048ba31 100644
--- a/test/multidomain/2cnistokes2p2cni/test_references/test_2cnistokes2p2cni_reference.input
+++ b/test/multidomain/2cnistokes2p2cni/test_references/test_2cnistokes2p2cni_reference.input
@@ -78,10 +78,18 @@ StabilizationAlpha = -1.0
#############################################################
[FreeFlow]
#############################################################
+RefVelocity = 3.5
RefPressure = 1e5
+RefMassfrac = 0.008 # 0.01
RefTemperature = 298.15
-RefMassfrac = 0.008
-VxMax = 3.5
+SinusVelAmplitude = 0.0 # 0.1
+SinusVelPeriod = 3600 # in seconds
+SinusPressureAmplitude = 0.0 # 0.5
+SinusPressurePeriod = 3600 # in seconds
+SinusConcentrationAmplitude = 0.0 # 0.5
+SinusConcentrationPeriod = 3600 # in seconds
+SinusTemperatureAmplitude = 0.0 # 0.5
+SinusTemperaturePeriod = 3600 # in seconds
BeaversJosephSlipVel = 0.00134
ExponentMTC = 0.0 # 1./6., Mass transfer coefficient for S^MTC
UseBoundaryLayerModel = 0 # 0 for no boundary layer model, 1 for Blasius, 2 and 3 for turbulent BL, 9 for constant thickness
@@ -89,11 +97,6 @@ BoundaryLayerOffset = 0.0 # for BL model like Blasius, determines a virtual ru
ConstThickness = 0.0016 # for a constant BL thickness, use BL model 9
MassTransferModel = 0 # 0 for none, 1 for power law, 2 for Schluender model
-SinusVelVar = 0.0 # hourly velocity variation
-SinusPVar = 0.0 # diurnal pressure variation
-SinusTVar = 0.0 # diurnal temperature variation
-SinusXVar = 0.0 # diurnal variation of the vapor concentration (massfraction)
-
#############################################################
[PorousMedium]
#############################################################
diff --git a/test/multidomain/2cstokes2p2c/stokes2csubproblem.hh b/test/multidomain/2cstokes2p2c/stokes2csubproblem.hh
index cf2f0b581e..f5bfe913c9 100644
--- a/test/multidomain/2cstokes2p2c/stokes2csubproblem.hh
+++ b/test/multidomain/2cstokes2p2c/stokes2csubproblem.hh
@@ -202,18 +202,26 @@ public:
bboxMin_[1] = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Grid, InterfacePos);
bboxMax_[1] = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Grid, YMax);
- refTemperature_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, RefTemperature);
+ refVelocity_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, RefVelocity);
refPressure_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, RefPressure);
refMassfrac_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, RefMassfrac);
- vxMax_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, VxMax);
- bjSlipVel_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, BeaversJosephSlipVel);
- sinusVelVar_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusVelVar);
+ refTemperature_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, RefTemperature);
+
+ sinusVAmplitude_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusVelAmplitude);
+ sinusVPeriod_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusVelPeriod);
+ sinusPAmplitude_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusPressureAmplitude);
+ sinusPPeriod_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusPressurePeriod);
+ sinusXAmplitude_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusConcentrationAmplitude);
+ sinusXPeriod_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusConcentrationPeriod);
+ sinusTAmplitude_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusTemperatureAmplitude);
+ sinusTPeriod_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, SinusTemperaturePeriod);
+ bjSlipVel_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, FreeFlow, BeaversJosephSlipVel);
+ alphaBJ_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, SpatialParams, AlphaBJ);
xMaterialInterface_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, SpatialParams, MaterialInterfaceX);
+
runUpDistanceX_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Problem, RunUpDistanceX); // first part of the interface without coupling
initializationTime_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, TimeManager, InitTime);
-
- alphaBJ_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, SpatialParams, AlphaBJ);
}
catch (Dumux::ParameterException &e) {
std::cerr << e << ". Abort!\n";
@@ -338,9 +346,9 @@ public:
values[velocityXIdx] = xVelocity_(globalPos);
values[velocityYIdx] = 0.0;
- values[pressureIdx] = refPressure_
+ values[pressureIdx] = refPressure()
+ density*this->gravity()[1]*(globalPos[1] - bboxMin_[1]);
- values[massOrMoleFracIdx] = refMassfrac_;
+ values[massOrMoleFracIdx] = refMassfrac();
}
/*!
@@ -379,7 +387,7 @@ public:
&& globalPos[1] > bboxMin_[1] && globalPos[1] < bboxMax_[1])
{
// rho*v*X at inflow
- values[transportEqIdx] = -xVelocity*density*refMassfrac_;
+ values[transportEqIdx] = -xVelocity * density * refMassfrac();
}
}
@@ -512,17 +520,21 @@ public:
const SpatialParams &spatialParams() const
{ return spatialParams_; }
+ //! \brief Returns the reference velocity.
+ const Scalar refVelocity() const
+ { return refVelocity_ + variation_(sinusVAmplitude_, sinusVPeriod_); }
+
//! \brief Returns the reference pressure.
const Scalar refPressure() const
- { return refPressure_; }
-
- //! \brief Returns the reference temperature.
- const Scalar refTemperature() const
- { return refTemperature_; }
+ { return refPressure_ + variation_(sinusPAmplitude_, sinusPPeriod_); }
//! \brief Returns the reference mass fraction.
const Scalar refMassfrac() const
- { return refMassfrac_; }
+ { return refMassfrac_ + variation_(sinusXAmplitude_, sinusXPeriod_); }
+
+ //! \brief Returns the reference temperature.
+ const Scalar refTemperature() const
+ { return refTemperature_+ variation_(sinusTAmplitude_, sinusTPeriod_); }
private:
/*!
@@ -541,15 +553,15 @@ private:
values[velocityXIdx] = xVelocity_(globalPos);
values[velocityYIdx] = 0.;
- values[pressureIdx] = refPressure_
+ values[pressureIdx] = refPressure()
+ density*this->gravity()[1]*(globalPos[1] - bboxMin_[1]);
- values[massOrMoleFracIdx] = refMassfrac_;
+ values[massOrMoleFracIdx] = refMassfrac();
}
//! \brief set the profile of the inflow velocity (horizontal direction)
const Scalar xVelocity_(const GlobalPosition &globalPos) const
{
- const Scalar vmax = vxMax_ + hourlyVariation_(sinusVelVar_);
+ const Scalar vmax = refVelocity();
// parabolic profile
return 4*vmax*(globalPos[1] - bboxMin_[1])*(bboxMax_[1] - globalPos[1])
@@ -569,7 +581,7 @@ private:
void updateFluidStateForBC_(FluidState& fluidState) const
{
fluidState.setTemperature(refTemperature());
- fluidState.setPressure(phaseIdx, refPressure_);
+ fluidState.setPressure(phaseIdx, refPressure());
Scalar massFraction[numComponents];
massFraction[transportCompIdx] = refMassfrac();
@@ -585,19 +597,9 @@ private:
fluidState.setMoleFraction(phaseIdx, phaseCompIdx, massFraction[phaseCompIdx]*M1/massToMoleDenominator);
}
- // can be used for the diurnal variation of a boundary condition
- const Scalar diurnalVariation_(const Scalar value) const
- {
- const Scalar time = this->timeManager().time();
- return sin(2*M_PI*time/86400) * value;
- }
-
- // can be used for the hourly variation of a boundary condition
- const Scalar hourlyVariation_(const Scalar value) const
- {
- const Scalar time = this->timeManager().time();
- return sin(2*M_PI*time/3600) * value;
- }
+ // can be used for the variation of a boundary condition
+ const Scalar variation_(const Scalar amplitude, const Scalar period) const
+ { return sin(2*M_PI*this->timeManager().time()/period) * amplitude; }
bool onLeftBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[0] < bboxMin_[0] + eps_; }
@@ -628,13 +630,21 @@ private:
GlobalPosition bboxMin_;
GlobalPosition bboxMax_;
+ Scalar refVelocity_;
Scalar refPressure_;
- Scalar refTemperature_;
Scalar refMassfrac_;
+ Scalar refTemperature_;
+
+ Scalar sinusVAmplitude_;
+ Scalar sinusVPeriod_;
+ Scalar sinusPAmplitude_;
+ Scalar sinusPPeriod_;
+ Scalar sinusXAmplitude_;
+ Scalar sinusXPeriod_;
+ Scalar sinusTAmplitude_;
+ Scalar sinusTPeriod_;
- Scalar vxMax_;
Scalar bjSlipVel_;
- Scalar sinusVelVar_;
Scalar alphaBJ_;
Scalar xMaterialInterface_;
diff --git a/test/multidomain/2cstokes2p2c/test_2cstokes2p2c.input b/test/multidomain/2cstokes2p2c/test_2cstokes2p2c.input
index c3bf32f300..e2befdcdbd 100644
--- a/test/multidomain/2cstokes2p2c/test_2cstokes2p2c.input
+++ b/test/multidomain/2cstokes2p2c/test_2cstokes2p2c.input
@@ -78,12 +78,19 @@ StabilizationAlpha = -1.0
#############################################################
[FreeFlow]
#############################################################
+RefVelocity = 3.5
RefPressure = 1e5
+RefMassfrac = 0.008 # 0.01
RefTemperature = 298.15
-RefMassfrac = 0.008
-VxMax = 3.5
+SinusVelAmplitude = 0.0 # 0.1
+SinusVelPeriod = 3600 # in seconds
+SinusPressureAmplitude = 0.0 # 0.5
+SinusPressurePeriod = 3600 # in seconds
+SinusConcentrationAmplitude = 0.0 # 0.5
+SinusConcentrationPeriod = 3600 # in seconds
+SinusTemperatureAmplitude = 0.0 # 0.5
+SinusTemperaturePeriod = 3600 # in seconds
BeaversJosephSlipVel = 0.00134
-SinusVelVar = 0.0
ExponentMTC = 0.0 # 1./6., Mass transfer coefficient for S^MTC
UseBoundaryLayerModel = 0 # 0 for no boundary layer model, 1 for Blasius, 2 and 3 for turbulent BL, 9 for constant thickness
BoundaryLayerOffset = 0.0 # for BL model like Blasius, determines a virtual run-up distance for the flow
diff --git a/test/multidomain/2cstokes2p2c/test_references/test_2cstokes2p2c_reference.input b/test/multidomain/2cstokes2p2c/test_references/test_2cstokes2p2c_reference.input
index 63d697d9bb..9e0ee3fe87 100644
--- a/test/multidomain/2cstokes2p2c/test_references/test_2cstokes2p2c_reference.input
+++ b/test/multidomain/2cstokes2p2c/test_references/test_2cstokes2p2c_reference.input
@@ -78,12 +78,19 @@ StabilizationAlpha = -1.0
#############################################################
[FreeFlow]
#############################################################
+RefVelocity = 3.5
RefPressure = 1e5
+RefMassfrac = 0.008 # 0.01
RefTemperature = 298.15
-RefMassfrac = 0.008
-VxMax = 3.5
+SinusVelAmplitude = 0.0 # 0.1
+SinusVelPeriod = 3600 # in seconds
+SinusPressureAmplitude = 0.0 # 0.5
+SinusPressurePeriod = 3600 # in seconds
+SinusConcentrationAmplitude = 0.0 # 0.5
+SinusConcentrationPeriod = 3600 # in seconds
+SinusTemperatureAmplitude = 0.0 # 0.5
+SinusTemperaturePeriod = 3600 # in seconds
BeaversJosephSlipVel = 0.00134
-SinusVelVar = 0.0
ExponentMTC = 0.0 # 1./6., Mass transfer coefficient for S^MTC
UseBoundaryLayerModel = 0 # 0 for no boundary layer model, 1 for Blasius, 2 and 3 for turbulent BL, 9 for constant thickness
BoundaryLayerOffset = 0.0 # for BL model like Blasius, determines a virtual run-up distance for the flow
--
GitLab