Commit 861b8bca authored by Katharina Heck's avatar Katharina Heck
Browse files

update files for review

parent b4cb9f26
add_subdirectory(example1)
add_subdirectory(example2)
add_subdirectory(example2_turbulent)
add_subdirectory(example3)
add_subdirectory(example4)
......@@ -3,5 +3,3 @@ add_input_file_links()
dune_add_test(NAME test_stokes1p2cdarcy2p2c_radiation
SOURCES main.cc
CMAKE_GUARD HAVE_UMFPACK)
set(CMAKE_BUILD_TYPE Debug)
......@@ -30,7 +30,7 @@
#include <dumux/porousmediumflow/2p2c/model.hh>
#include <dumux/porousmediumflow/problem.hh>
#include <dumux/material/fluidsystems/h2oair.hh>
#include <dumux/h2oair.hh>
#include <dumux/io/gnuplotinterface.hh>
......@@ -54,6 +54,7 @@ template<class TypeTag>
struct Problem<TypeTag, TTag::DarcyTwoPTwoCTypeTag> { using type = Dumux::DarcySubProblem<TypeTag>; };
// the fluid system
// for silt the kelvin equation has to be switched to true
template<class TypeTag>
struct FluidSystem<TypeTag, TTag::DarcyTwoPTwoCTypeTag> { using type = FluidSystems::H2OAir<GetPropType<TypeTag, Properties::Scalar>>; };
......@@ -135,6 +136,17 @@ public:
time_ = 0.0;
initializationTime_ = getParamFromGroup<Scalar>(this->paramGroup(), "Problem.InitializationTime");
problemName_ = getParam<std::string>("Vtk.OutputName") + "_" + getParamFromGroup<std::string>(this->paramGroup(), "Problem.Name");
}
/*!
* \brief The problem name.
*/
const std::string& name() const
{
return problemName_;
}
void setTime(Scalar time)
......@@ -163,9 +175,9 @@ public:
Scalar energyFlux = 0.0;
Scalar netRadiation = 0.0;
// bulk elements
for (const auto& element : elements(this->fvGridGeometry().gridView()))
for (const auto& element : elements(this->gridGeometry().gridView()))
{
auto fvGeometry = localView(this->fvGridGeometry());
auto fvGeometry = localView(this->gridGeometry());
fvGeometry.bindElement(element);
auto elemVolVars = localView(gridVariables.curGridVolVars());
......@@ -219,8 +231,7 @@ public:
gnuplot_.setYRange(0, 1e-5);
gnuplot_.setXlabel("time [s]");
gnuplot_.setYlabel("kg/s");
gnuplot_.addDataSetToPlot(x_, y_, "evaporation");
gnuplot_.plot("evaporation");
gnuplot_.addDataSetToPlot(x_, y_, name()+"evaporation");
//do a gnuplot
y2_.push_back(massWater);
......@@ -230,8 +241,8 @@ public:
gnuplot2_.setYRange(0, 20);
gnuplot2_.setXlabel("time [s]");
gnuplot2_.setYlabel("kg");
gnuplot2_.addDataSetToPlot(x_, y2_, "water mass");
gnuplot2_.plot("watermass");
gnuplot2_.addDataSetToPlot(x_, y2_, name()+"water mass");
//do a gnuplot
y3_.push_back(averageTemperature);
......@@ -241,8 +252,8 @@ public:
gnuplot3_.setYRange(280,300);
gnuplot3_.setXlabel("time [s]");
gnuplot3_.setYlabel("K");
gnuplot3_.addDataSetToPlot(x_, y3_, "temperatureAverage");
gnuplot3_.plot("temperature");
gnuplot3_.addDataSetToPlot(x_, y3_, name()+"temperatureAverage");
//do a gnuplot
y4_.push_back(netRadiation);
......@@ -251,8 +262,8 @@ public:
gnuplot4_.setYRange(0,600);
gnuplot4_.setXlabel("time [s]");
gnuplot4_.setYlabel("W/m^2");
gnuplot4_.addDataSetToPlot(x_, y4_, "net radiation");
gnuplot4_.plot("radiation");
gnuplot4_.addDataSetToPlot(x_, y4_, name()+"net radiation");
}
}
......@@ -287,13 +298,10 @@ public:
BoundaryTypes boundaryTypes(const Element &element, const SubControlVolumeFace &scvf) const
{
BoundaryTypes values;
const auto& globalPos = scvf.center();
values.setAllNeumann();
if (couplingManager().isCoupledEntity(CouplingManager::darcyIdx, scvf))
values.setAllCouplingNeumann();
// else if (onLowerBoundary_(globalPos))
// values.setAllDirichlet();
else
values.setAllNeumann();
return values;
......@@ -311,7 +319,6 @@ public:
{
PrimaryVariables values(0.0);
values = initialAtPos(scvf.center());
// values[switchIdx]= 0.8;
return values;
}
......@@ -398,7 +405,7 @@ public:
PrimaryVariables values(0.0);
values.setState(initialPhasePresence_);
values[pressureIdx] = pressure_ + 1. * this->spatialParams().gravity(globalPos)[1] * (globalPos[1] - this->fvGridGeometry().bBoxMax()[1]);
values[pressureIdx] = pressure_ + 1. * this->spatialParams().gravity(globalPos)[1] * (globalPos[1] - this->gridGeometry().bBoxMax()[1]);
values[switchIdx] = initialSw_;
values[energyEqIdx] = temperature_; //20
return values;
......@@ -420,16 +427,16 @@ public:
private:
bool onLeftBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[0] < this->fvGridGeometry().bBoxMin()[0] + eps_; }
{ return globalPos[0] < this->gridGeometry().bBoxMin()[0] + eps_; }
bool onRightBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[0] > this->fvGridGeometry().bBoxMax()[0] - eps_; }
{ return globalPos[0] > this->gridGeometry().bBoxMax()[0] - eps_; }
bool onLowerBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[1] < this->fvGridGeometry().bBoxMin()[1] + eps_; }
{ return globalPos[1] < this->gridGeometry().bBoxMin()[1] + eps_; }
bool onUpperBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[1] > this->fvGridGeometry().bBoxMax()[1] - eps_; }
{ return globalPos[1] > this->gridGeometry().bBoxMax()[1] - eps_; }
Scalar pressure_;
Scalar initialSw_;
......@@ -450,6 +457,8 @@ private:
Dumux::GnuplotInterface<double> gnuplot2_;
Dumux::GnuplotInterface<double> gnuplot3_;
Dumux::GnuplotInterface<double> gnuplot4_;
std::string problemName_;
};
} //end namespace
......
......@@ -4,31 +4,31 @@ MaxTimeStepSize = 3600 # [s] (12 hours)
TEnd = 864000 # [s] (10 days)
[Stokes.Grid]
Positions0 = -0.25 0.5
Positions0 = -0.25 0.0 0.25 1.0
Positions1 = 0.25 0.5
Grading0 = 1.0
Grading1 = 1.15
Cells0 = 75
Cells1 = 25
Grading0 = 1.0 1.0 1.0
Grading1 = 1.1
Cells0 = 25 25 50
Cells1 = 30
[Darcy.Grid]
Positions0 = 0.0 0.25
Positions1 = 0.0 0.25
Cells0 = 25
Cells1 = 25
Cells1 = 30
Grading0 = 1.0
Grading1 = -1.1
[Stokes.Problem]
Name = stokes
RefVelocity = 1 # [m/s]
Name = stokes_vel2
RefVelocity = 2 # [m/s]
RefPressure = 1e5 # [Pa]
refMoleFrac = 0.001 # [-]
refMoleFrac = 0.008 # [-]
RefTemperature = 283.15 # [K]
InitializationTime = 0.01
[Darcy.Problem]
Name = darcy
Name = darcy_vel2
Pressure = 1.0e5
Saturation = 0.7 # initial Sw
Temperature = 283.15 # [K]
......@@ -48,12 +48,13 @@ VgN = 8.0
[Problem]
Name = test_stokes1p2cdarcy2p2chorizontal
Name = example1
EnableGravity = true
[Vtk]
AddVelocity = true
WriteFaceData = false
OutputName = example1
[Component]
SolidDensity = 2700
......@@ -63,6 +64,3 @@ SolidHeatCapacity = 790
[RANS]
TurbulentSchmidtNumber = 0.7
TurbulentPrandtlNumber = 0.85
[Newton]
MaxRelativeShift = 1e-5
[TimeLoop]
DtInitial = 0.1 # [s]
MaxTimeStepSize = 3600 # [s] (12 hours)
TEnd = 864000 # [s] (10 days)
[Stokes.Grid]
Positions0 = -0.25 0.0 0.25 1.0
Positions1 = 0.25 0.5
Grading0 = 1.0 1.0 1.0
Grading1 = 1.1
Cells0 = 25 25 50
Cells1 = 30
[Darcy.Grid]
Positions0 = 0.0 0.25
Positions1 = 0.0 0.25
Cells0 = 25
Cells1 = 30
Grading0 = 1.0
Grading1 = -1.1
[Stokes.Problem]
Name = stokes_vel1_silt
RefVelocity = 1 # [m/s]
RefPressure = 1e5 # [Pa]
refMoleFrac = 0.008 # [-]
RefTemperature = 283.15 # [K]
InitializationTime = 0.01
[Darcy.Problem]
Name = darcy_vel1_silt
Pressure = 1.0e5
Saturation = 0.7 # initial Sw
Temperature = 283.15 # [K]
InitPhasePresence = 3 # bothPhases
InitializationTime = 0.01
PlexiglassThermalConductivity = 0.184 # [W/(m*K)] 0.116
PlexiglassThickness = 0.005 # [m]
[SpatialParams]
Porosity = 0.35
Permeability = 1.08e-12
AlphaBJ = 1.0
Swr = 0.057
Snr = 0.029
VgAlpha = 4.28e-5
VgN = 1.5
[Problem]
Name = example1
EnableGravity = true
[Vtk]
AddVelocity = true
WriteFaceData = false
OutputName = example1
[Component]
SolidDensity = 2700
SolidThermalConductivity = 2.8
SolidHeatCapacity = 790
[RANS]
TurbulentSchmidtNumber = 0.7
TurbulentPrandtlNumber = 0.85
......@@ -50,10 +50,10 @@ class TwoPTwoCSpatialParams
using Problem = GetPropType<TypeTag, Properties::Problem>;
using GridView = GetPropType<TypeTag, Properties::GridView>;
using Element = typename GridView::template Codim<0>::Entity;
using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>;
using FVElementGeometry = typename FVGridGeometry::LocalView;
using GridGeometry = GetPropType<TypeTag, Properties::GridGeometry>;
using FVElementGeometry = typename GridGeometry::LocalView;
using SubControlVolume = typename FVElementGeometry::SubControlVolume;
using ParentType = FVSpatialParams<FVGridGeometry, Scalar, TwoPTwoCSpatialParams<TypeTag>>;
using ParentType = FVSpatialParams<GridGeometry, Scalar, TwoPTwoCSpatialParams<TypeTag>>;
using GlobalPosition = typename Element::Geometry::GlobalCoordinate;
using EffectiveLaw = RegularizedVanGenuchten<Scalar>;
......@@ -64,8 +64,8 @@ public:
using PermeabilityType = Scalar;
TwoPTwoCSpatialParams(std::shared_ptr<const FVGridGeometry> fvGridGeometry)
: ParentType(fvGridGeometry)
TwoPTwoCSpatialParams(std::shared_ptr<const GridGeometry> gridGeometry)
: ParentType(gridGeometry)
{
permeability_ = getParam<Scalar>("SpatialParams.Permeability");
porosity_ = getParam<Scalar>("SpatialParams.Porosity");
......
......@@ -27,7 +27,7 @@
#include <dune/grid/yaspgrid.hh>
#include <dumux/material/fluidsystems/1padapter.hh>
#include <dumux/material/fluidsystems/h2oair.hh>
#include <dumux/h2oair.hh>
#include <dumux/freeflow/navierstokes/problem.hh>
#include <dumux/discretization/staggered/freeflow/properties.hh>
......@@ -74,7 +74,7 @@ template<class TypeTag>
struct Problem<TypeTag, TTag::StokesOnePTwoCTypeTag> { using type = Dumux::StokesSubProblem<TypeTag> ; };
template<class TypeTag>
struct EnableFVGridGeometryCache<TypeTag, TTag::StokesOnePTwoCTypeTag> { static constexpr bool value = true; };
struct EnableGridGeometryCache<TypeTag, TTag::StokesOnePTwoCTypeTag> { static constexpr bool value = true; };
template<class TypeTag>
struct EnableGridFluxVariablesCache<TypeTag, TTag::StokesOnePTwoCTypeTag> { static constexpr bool value = true; };
template<class TypeTag>
......@@ -142,7 +142,7 @@ public:
fluidState.setMoleFraction(phaseIdx, 0/*FluidSystem::H2OIdx*/, 0.0);
Scalar density = FluidSystem::density(fluidState, 0);
Scalar kinematicViscosity = FluidSystem::viscosity(fluidState, phaseIdx) / density;
Scalar diameter = this->fvGridGeometry().bBoxMax()[1] - this->fvGridGeometry().bBoxMin()[1];
Scalar diameter = this->gridGeometry().bBoxMax()[1] - this->gridGeometry().bBoxMin()[1];
turbulentKineticEnergy_ = turbulenceProperties.turbulentKineticEnergy(refVelocity_, diameter, kinematicViscosity);
dissipation_ = turbulenceProperties.dissipationRate(refVelocity_, diameter, kinematicViscosity);
......@@ -168,9 +168,9 @@ public:
Scalar evaporation = 0.0;
Scalar energyFlux = 0.0;
// bulk elements
for (const auto& element : elements(this->fvGridGeometry().gridView()))
for (const auto& element : elements(this->gridGeometry().gridView()))
{
auto fvGeometry = localView(this->fvGridGeometry());
auto fvGeometry = localView(this->gridGeometry());
fvGeometry.bindElement(element);
auto elemVolVars = localView(gridVariables.curGridVolVars());
......@@ -292,7 +292,7 @@ public:
const auto globalPos = scv.center();
PrimaryVariables values(initialAtPos(globalPos));
using std::pow;
unsigned int elementIdx = this->fvGridGeometry().elementMapper().index(element);
unsigned int elementIdx = this->gridGeometry().elementMapper().index(element);
const auto wallDistance = ParentType::wallDistance_[elementIdx];
values[Indices::dissipationEqIdx] = 6.0 * ParentType::kinematicViscosity_[elementIdx]
/ (ParentType::betaOmega() * pow(wallDistance, 2));
......@@ -374,11 +374,12 @@ public:
FluidState fluidState;
fluidState.setPressure(0, 1e5);
fluidState.setTemperature(refTemperature());
fluidState.setMoleFraction(0, 0, 1);
fluidState.setMoleFraction(0, 0, 1-refMoleFrac());
fluidState.setMoleFraction(0, 1, refMoleFrac());
Scalar density = FluidSystem::density(fluidState, 0);
PrimaryVariables values(0.0);
values[Indices::pressureIdx] = refPressure() + density*this->gravity()[1]*(globalPos[1] - this->fvGridGeometry().bBoxMin()[1]);
values[Indices::pressureIdx] = refPressure() + density*this->gravity()[1]*(globalPos[1] - this->gridGeometry().bBoxMin()[1]);
values[Indices::conti0EqIdx + 1] = refMoleFrac();
values[Indices::velocityXIdx] = refVelocity();
......@@ -455,20 +456,20 @@ public:
private:
bool onLeftBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[0] < this->fvGridGeometry().bBoxMin()[0] + eps_; }
{ return globalPos[0] < this->gridGeometry().bBoxMin()[0] + eps_; }
bool onRightBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[0] > this->fvGridGeometry().bBoxMax()[0] - eps_; }
{ return globalPos[0] > this->gridGeometry().bBoxMax()[0] - eps_; }
bool onLowerBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[1] < this->fvGridGeometry().bBoxMin()[1] + eps_; }
{ return globalPos[1] < this->gridGeometry().bBoxMin()[1] + eps_; }
bool onUpperBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[1] > this->fvGridGeometry().bBoxMax()[1] - eps_; }
{ return globalPos[1] > this->gridGeometry().bBoxMax()[1] - eps_; }
// the height of the free-flow domain
const Scalar height_() const
{ return this->fvGridGeometry().bBoxMax()[1] - this->fvGridGeometry().bBoxMin()[1]; }
{ return this->gridGeometry().bBoxMax()[1] - this->gridGeometry().bBoxMin()[1]; }
Scalar eps_;
......
add_input_file_links()
dune_add_test(NAME test_stokes1p2cdarcy2p2c_radiation_laminar_fixeddiffusion
SOURCES main.cc
CMAKE_GUARD HAVE_UMFPACK
......@@ -20,4 +18,4 @@ dune_add_test(NAME test_stokes1p2cdarcy2p2c_radiation_laminar_allvary
CMAKE_GUARD HAVE_UMFPACK
COMPILE_DEFINITIONS FLUIDSYSTEM=FluidSystems::H2OAir)
set(CMAKE_BUILD_TYPE Debug)
add_input_file_links()
......@@ -137,6 +137,16 @@ public:
time_ = 0.0;
initializationTime_ = getParamFromGroup<Scalar>(this->paramGroup(), "Problem.InitializationTime");
problemName_ = getParam<std::string>("Vtk.OutputName") + "_" + getParamFromGroup<std::string>(this->paramGroup(), "Problem.Name");
}
/*!
* \brief The problem name.
*/
const std::string& name() const
{
return problemName_;
}
void setTime(Scalar time)
......@@ -165,9 +175,9 @@ public:
Scalar energyFlux = 0.0;
Scalar netRadiation = 0.0;
// bulk elements
for (const auto& element : elements(this->fvGridGeometry().gridView()))
for (const auto& element : elements(this->gridGeometry().gridView()))
{
auto fvGeometry = localView(this->fvGridGeometry());
auto fvGeometry = localView(this->gridGeometry());
fvGeometry.bindElement(element);
auto elemVolVars = localView(gridVariables.curGridVolVars());
......@@ -221,8 +231,7 @@ public:
gnuplot_.setYRange(0, 1e-5);
gnuplot_.setXlabel("time [s]");
gnuplot_.setYlabel("kg/s");
gnuplot_.addDataSetToPlot(x_, y_, "evaporation");
gnuplot_.plot("evaporation");
gnuplot_.addDataSetToPlot(x_, y_, name()+"evaporation");
//do a gnuplot
y2_.push_back(massWater);
......@@ -232,8 +241,7 @@ public:
gnuplot2_.setYRange(0, 20);
gnuplot2_.setXlabel("time [s]");
gnuplot2_.setYlabel("kg");
gnuplot2_.addDataSetToPlot(x_, y2_, "water mass");
gnuplot2_.plot("watermass");
gnuplot2_.addDataSetToPlot(x_, y2_, name()+"water mass");
//do a gnuplot
y3_.push_back(averageTemperature);
......@@ -243,8 +251,7 @@ public:
gnuplot3_.setYRange(280,300);
gnuplot3_.setXlabel("time [s]");
gnuplot3_.setYlabel("K");
gnuplot3_.addDataSetToPlot(x_, y3_, "temperatureAverage");
gnuplot3_.plot("temperature");
gnuplot3_.addDataSetToPlot(x_, y3_, name()+"temperatureAverage");
//do a gnuplot
y4_.push_back(netRadiation);
......@@ -253,8 +260,7 @@ public:
gnuplot4_.setYRange(0,600);
gnuplot4_.setXlabel("time [s]");
gnuplot4_.setYlabel("W/m^2");
gnuplot4_.addDataSetToPlot(x_, y4_, "net radiation");
gnuplot4_.plot("radiation");
gnuplot4_.addDataSetToPlot(x_, y4_, name()+"net radiation");
}
}
......@@ -289,13 +295,10 @@ public:
BoundaryTypes boundaryTypes(const Element &element, const SubControlVolumeFace &scvf) const
{
BoundaryTypes values;
const auto& globalPos = scvf.center();
values.setAllNeumann();
if (couplingManager().isCoupledEntity(CouplingManager::darcyIdx, scvf))
values.setAllCouplingNeumann();
// else if (onLowerBoundary_(globalPos))
// values.setAllDirichlet();
else
values.setAllNeumann();
return values;
......@@ -400,7 +403,7 @@ public:
PrimaryVariables values(0.0);
values.setState(initialPhasePresence_);
values[pressureIdx] = pressure_ + 1. * this->spatialParams().gravity(globalPos)[1] * (globalPos[1] - this->fvGridGeometry().bBoxMax()[1]);
values[pressureIdx] = pressure_ + 1. * this->spatialParams().gravity(globalPos)[1] * (globalPos[1] - this->gridGeometry().bBoxMax()[1]);
values[switchIdx] = initialSw_;
values[energyEqIdx] = temperature_; //20
return values;
......@@ -422,16 +425,16 @@ public:
private:
bool onLeftBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[0] < this->fvGridGeometry().bBoxMin()[0] + eps_; }
{ return globalPos[0] < this->gridGeometry().bBoxMin()[0] + eps_; }
bool onRightBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[0] > this->fvGridGeometry().bBoxMax()[0] - eps_; }
{ return globalPos[0] > this->gridGeometry().bBoxMax()[0] - eps_; }
bool onLowerBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[1] < this->fvGridGeometry().bBoxMin()[1] + eps_; }
{ return globalPos[1] < this->gridGeometry().bBoxMin()[1] + eps_; }
bool onUpperBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[1] > this->fvGridGeometry().bBoxMax()[1] - eps_; }
{ return globalPos[1] > this->gridGeometry().bBoxMax()[1] - eps_; }
Scalar pressure_;
Scalar initialSw_;
......@@ -452,6 +455,8 @@ private:
Dumux::GnuplotInterface<double> gnuplot2_;
Dumux::GnuplotInterface<double> gnuplot3_;
Dumux::GnuplotInterface<double> gnuplot4_;
std::string problemName_;
};
} //end namespace
......
......@@ -69,7 +69,7 @@ struct H2OAirDefaultPolicy
template <class Scalar,
class H2Otype = Components::TabulatedComponent<Components::H2O<Scalar> >,
class Policy = H2OAirDefaultPolicy<>,
bool useKelvinVaporPressure = false>
bool useKelvinVaporPressure = true>
class H2OAir
: public Base<Scalar, H2OAir<Scalar, H2Otype, Policy> >
{
......@@ -272,7 +272,7 @@ public:
? fluidState.pressure(AirIdx)-fluidState.pressure(H2OIdx)
: fluidState.pressure(H2OIdx)-fluidState.pressure(AirIdx);
return H2O::vaporPressure(t)*exp( -pc * molarMass(H2OIdx)
/ density(fluidState, H2OIdx)
/ H2O::liquidDensity(t, fluidState.pressure(H2OIdx))
/ (Dumux::Constants<Scalar>::R*t) );
}
}
......
......@@ -51,7 +51,7 @@ template<bool fastButSimplifiedRelations = false>
struct H2OAirFixedDensityDefaultPolicy
{
static constexpr bool useH2ODensityAsLiquidMixtureDensity() { return fastButSimplifiedRelations; }
static constexpr bool useIdealGasDensity() { return fastButSimplifiedRelations; }
static constexpr bool useIdealGasDensity() { return true; }
static constexpr bool useAirViscosityAsGasMixtureViscosity() { return fastButSimplifiedRelations; }
};
......
[TimeLoop]
DtInitial = 0.1 # [s]
MaxTimeStepSize = 1800 # [s] (12 hours)
TEnd = 86400 # [s] (15 hours)
TEnd = 864000 # [s] (15 hours)
EpisodeLength = 1800
[Stokes.Grid]
Positions0 = -0.5 0.0 0.5
Positions0 = -0.25 0.0 0.25 1.0
Positions1 = 0.25 0.5
Grading0 = 1.0 1.0
Grading0 = 1.0 1.0 1.0
Grading1 = 1.1
Cells0 = 25 50
Cells1 = 40
Cells0 = 25 25 50
Cells1 = 30
[Darcy.Grid]
Positions0 = 0.0 0.25
......@@ -18,18 +18,18 @@ Positions1 = 0.0 0.25
Cells0 = 25
Cells1 = 30
Grading0 = 1.0
Grading1 = -1.2
Grading1 = -1.1
[Stokes.Problem]
Name = stokes_fixedviscosity_noGrav
Name = stokes_fixeddensity
RefVelocity = 0.1 # [m/s]
RefPressure = 1e5 # [Pa]
refMoleFrac = 0.001 # [-]
refMoleFrac = 0.008 # [-]
RefTemperature = 283.15 # [K]
InitializationTime = 0.01
[Darcy.Problem]
Name = darcy_fixedviscosity_noGrav
Name = darcy_fixeddensity
Pressure = 1.0e5
Saturation = 0.7 # initial Sw
Temperature = 283.15 # [K]
......@@ -49,11 +49,12 @@ VgN = 8.0
[Problem]
Name = test_stokes1p2cdarcy2p2chorizontal
EnableGravity = false
EnableGravity = true
[Vtk]
AddVelocity = true