diff --git a/.gitignore b/.gitignore
index 4e163b8ed45588bc842a729f41acde13450fb301..c29fadbc9934524e6c2fad7b6d7a9c3106ba7dd6 100644
--- a/.gitignore
+++ b/.gitignore
@@ -1,10 +1,9 @@
-# SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
+# SPDX-FileCopyrightInfo: Copyright © DuMux-Course contributors, see AUTHORS.md in root folder
# SPDX-License-Identifier: CC0-1.0
# build system clutter
build-*
Testing
-!bin/testing
# auto-saved files
*~
@@ -12,14 +11,12 @@ Testing
# hidden files
.cproject
.project
-.vscode
# left overs from git rebase
*.orig
*.rej
# latex clutter
-*.pdf
*.aux
*.blg
*.log
@@ -29,14 +26,6 @@ Testing
*.out
*.tdo
*.toc
-*.synctex.gz
-
-# Python clutter
-*.pyc
-__pycache__
-
-# macOS
-.DS_Store
# always consider files containing source code regardless of their name
!*.cc
@@ -48,9 +37,3 @@ __pycache__
# always consider reference solutions
!*reference.vtu
-
-# ignore files generated during python setup.py sdist
-MANIFEST
-_skbuild/
-dist
-*.egg-info
diff --git a/.patches/exercise-basic/exercise-basic.patch b/.patches/exercise-basic/exercise-basic.patch
index fa829f8d0a6da3c50f346c7d817c44b14d84313b..5b19ef068489b9a70d47d5b3967718a0b7e5343a 100644
--- a/.patches/exercise-basic/exercise-basic.patch
+++ b/.patches/exercise-basic/exercise-basic.patch
@@ -1,5 +1,5 @@
diff -ruN exercises/exercise-basic/2pmain.cc exercises/solution/exercise-basic/2pmain.cc
---- exercises/exercise-basic/2pmain.cc 2024-07-17 14:02:52.088041093 +0200
+--- exercises/exercise-basic/2pmain.cc 2024-07-17 13:58:18.204850249 +0200
+++ exercises/solution/exercise-basic/2pmain.cc 1970-01-01 01:00:00.000000000 +0100
@@ -1,145 +0,0 @@
-// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
@@ -149,7 +149,7 @@ diff -ruN exercises/exercise-basic/2pmain.cc exercises/solution/exercise-basic/2
-} // end main
diff -ruN exercises/exercise-basic/2pnimain.cc exercises/solution/exercise-basic/2pnimain.cc
--- exercises/exercise-basic/2pnimain.cc 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-basic/2pnimain.cc 2024-07-17 14:02:52.092041122 +0200
++++ exercises/solution/exercise-basic/2pnimain.cc 2024-07-17 13:58:18.212850366 +0200
@@ -0,0 +1,138 @@
+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
+// vi: set et ts=4 sw=4 sts=4:
@@ -290,8 +290,8 @@ diff -ruN exercises/exercise-basic/2pnimain.cc exercises/solution/exercise-basic
+ return 0;
+} // end main
diff -ruN exercises/exercise-basic/CMakeLists.txt exercises/solution/exercise-basic/CMakeLists.txt
---- exercises/exercise-basic/CMakeLists.txt 2024-07-17 14:37:31.109952737 +0200
-+++ exercises/solution/exercise-basic/CMakeLists.txt 2024-07-17 14:02:52.092041122 +0200
+--- exercises/exercise-basic/CMakeLists.txt 2024-07-17 14:10:30.239561158 +0200
++++ exercises/solution/exercise-basic/CMakeLists.txt 2024-07-17 14:11:44.100637927 +0200
@@ -1,12 +1,9 @@
# SPDX-FileCopyrightInfo: Copyright © DuMux-Course contributors, see AUTHORS.md in root folder
# SPDX-License-Identifier: GPL-3.0-or-later
@@ -308,242 +308,9 @@ diff -ruN exercises/exercise-basic/CMakeLists.txt exercises/solution/exercise-ba
# add a symlink for each input file
add_input_file_links()
-diff -ruN exercises/exercise-basic/injection2p2cproblem.hh exercises/solution/exercise-basic/injection2p2cproblem.hh
---- exercises/exercise-basic/injection2p2cproblem.hh 2024-07-17 14:37:31.221955214 +0200
-+++ exercises/solution/exercise-basic/injection2p2cproblem.hh 1970-01-01 01:00:00.000000000 +0100
-@@ -1,229 +0,0 @@
--// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
--// vi: set et ts=4 sw=4 sts=4:
--/*****************************************************************************
-- * See the file COPYING for full copying permissions. *
-- * *
-- * This program is free software: you can redistribute it and/or modify *
-- * it under the terms of the GNU General Public License as published by *
-- * the Free Software Foundation, either version 3 of the License, or *
-- * (at your option) any later version. *
-- * *
-- * This program is distributed in the hope that it will be useful, *
-- * but WITHOUT ANY WARRANTY; without even the implied warranty of *
-- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
-- * GNU General Public License for more details. *
-- * *
-- * You should have received a copy of the GNU General Public License *
-- * along with this program. If not, see <http://www.gnu.org/licenses/>. *
-- *****************************************************************************/
--/*!
-- * \file
-- *
-- * \brief The two-phase porousmediumflow problem for exercise-basic
-- */
--
--#ifndef DUMUX_EX_BASIC_PROBLEM_2P2C_HH
--#define DUMUX_EX_BASIC_PROBLEM_2P2C_HH
--
--#include <dumux/common/properties.hh>
--#include <dumux/common/boundarytypes.hh>
--#include <dumux/common/numeqvector.hh>
--#include <dumux/porousmediumflow/problem.hh>
--#include <dumux/material/binarycoefficients/h2o_n2.hh>
--
--namespace Dumux {
--
--/*!
-- * \ingroup TwoPTwoCModel
-- * \ingroup ImplicitTestProblems
-- * \brief Gas injection problem where a gas (here nitrogen) is injected into a fully
-- * water saturated medium. During buoyancy-driven upward migration the gas
-- * passes a high temperature area.
-- *
-- * The domain is sized 60 m times 40 m.
-- *
-- * For the mass conservation equation Neumann boundary conditions are used on
-- * the top, on the bottom and on the right of the domain, while Dirichlet conditions
-- * apply on the left boundary.
-- *
-- * Gas is injected at the right boundary from 7 m to 15 m at a rate of
-- * 0.001 kg/(s m), the remaining Neumann boundaries are no-flow
-- * boundaries.
-- *
-- * At the Dirichlet boundaries a hydrostatic pressure and a gas saturation of zero a
-- *
-- * This problem uses the \ref TwoPModel model.
-- */
--template<class TypeTag>
--class Injection2p2cProblem : public PorousMediumFlowProblem<TypeTag>
--{
-- using ParentType = PorousMediumFlowProblem<TypeTag>;
-- using Scalar = GetPropType<TypeTag, Properties::Scalar>;
-- using Indices = typename GetPropType<TypeTag, Properties::ModelTraits>::Indices;
-- using PrimaryVariables = GetPropType<TypeTag, Properties::PrimaryVariables>;
-- using BoundaryTypes = Dumux::BoundaryTypes<GetPropType<TypeTag, Properties::ModelTraits>::numEq()>;
-- using GridGeometry = GetPropType<TypeTag, Properties::GridGeometry>;
-- using FVElementGeometry = typename GridGeometry::LocalView;
-- using GridView = typename GridGeometry::GridView;
-- using FluidSystem = GetPropType<TypeTag, Properties::FluidSystem>;
-- using NumEqVector = Dumux::NumEqVector<PrimaryVariables>;
--
-- static constexpr int dimWorld = GridView::dimensionworld;
-- using Element = typename GridView::template Codim<0>::Entity;
-- using GlobalPosition = typename Element::Geometry::GlobalCoordinate;
--
--public:
-- Injection2p2cProblem(std::shared_ptr<const GridGeometry> gridGeometry)
-- : ParentType(gridGeometry)
-- {
-- // Initialize the tables of the fluid system
-- FluidSystem::init(/*tempMin=*/273.15,
-- /*tempMax=*/423.15,
-- /*numTemp=*/50,
-- /*pMin=*/0.0,
-- /*pMax=*/30e6,
-- /*numP=*/300);
--
-- // Name of the problem and output file
-- // getParam<TYPE>("GROUPNAME.PARAMNAME") reads and sets parameter PARAMNAME
-- // of type TYPE given in the group GROUPNAME from the input file
-- name_ = getParam<std::string>("Problem.Name");
-- // Depth of the aquifer, unit: m
-- aquiferDepth_ = getParam<Scalar>("Problem.AquiferDepth");
-- // The duration of the injection, unit: seconds
-- injectionDuration_ = getParam<Scalar>("Problem.InjectionDuration");
-- }
--
-- /*!
-- * \brief Returns the problem name
-- *
-- * This is used as a prefix for files generated by the simulation.
-- */
-- std::string name() const
-- { return name_+"-2p2c"; }
--
-- /*!
-- * \brief Specifies which kind of boundary condition should be
-- * used for which equation on a given boundary segment.
-- *
-- * \param globalPos The position for which the bc type should be evaluated
-- */
-- BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const
-- {
-- BoundaryTypes bcTypes;
-- if (globalPos[0] < eps_)
-- bcTypes.setAllDirichlet();
-- else
-- bcTypes.setAllNeumann();
--
-- return bcTypes;
-- }
--
-- /*!
-- * \brief Evaluates the boundary conditions for a Dirichlet
-- * boundary segment
-- *
-- * \param globalPos The global position
-- */
-- PrimaryVariables dirichletAtPos(const GlobalPosition &globalPos) const
-- {
-- return initialAtPos(globalPos);
-- }
--
-- /*!
-- * \brief Evaluate the boundary conditions for a Neumann
-- * boundary segment.
-- *
-- * \param globalPos The position of the integration point of the boundary segment.
-- */
-- NumEqVector neumannAtPos(const GlobalPosition &globalPos) const
-- {
-- // Initialize values to zero, i.e. no-flow Neumann boundary conditions
-- NumEqVector values(0.0);
--
-- // If we are inside the injection zone set inflow Neumann boundary conditions
-- if (injectionActive() && onInjectionBoundary(globalPos))
-- {
-- // TODO: dumux-course-task
-- // Instead of setting -1e-4 here directly use totalAreaSpecificInflow_ in the computation.
--
-- // Inject nitrogen. negative values mean injection.
-- // Convert units from kg/(s*m^2) to mole/(s*m^2)
-- values[Indices::conti0EqIdx + FluidSystem::N2Idx] = -1e-4/FluidSystem::molarMass(FluidSystem::N2Idx);
-- values[Indices::conti0EqIdx + FluidSystem::H2OIdx] = 0.0;
-- }
--
-- return values;
-- }
--
-- /*!
-- * \brief Evaluate the source term for all phases within a given
-- * sub-control-volume.
-- *
-- * \param globalPos The position for which the source term should be evaluated
-- */
-- NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
-- {
-- return NumEqVector(0.0);
-- }
--
-- /*!
-- * \brief Evaluate the initial value for a control volume.
-- *
-- * \param globalPos The position for which the initial condition should be evaluated.
-- */
-- PrimaryVariables initialAtPos(const GlobalPosition &globalPos) const
-- {
-- PrimaryVariables values(0.0);
-- values.setState(Indices::firstPhaseOnly);
-- // Get the water density at atmospheric conditions
-- const Scalar densityW = FluidSystem::H2O::liquidDensity(this->spatialParams().temperatureAtPos(globalPos), 1.0e5);
--
-- // Assume an initially hydrostatic liquid pressure profile
-- // Note: we subtract rho_w*g*h because g is defined negative
-- const Scalar pw = 1.0e5 - densityW*this->spatialParams().gravity(globalPos)[dimWorld-1]*(aquiferDepth_ - globalPos[dimWorld-1]);
--
-- // Initially we have some nitrogen dissolved
-- // Saturation mole fraction would be:
-- // moleFracLiquidN2 = (pw + pc + p_vap^sat)/henry;
-- const Scalar moleFracLiquidN2 = pw*0.95/BinaryCoeff::H2O_N2::henry(this->spatialParams().temperatureAtPos(globalPos));
--
-- // Note that because we start with a single phase system the primary variables
-- // are pl and x^w_N2. This will switch as soon after we start injecting to a two
-- // phase system so the primary variables will be pl and Sn (nonwetting saturation).
-- values[Indices::pressureIdx] = pw;
-- values[Indices::switchIdx] = moleFracLiquidN2;
--
-- return values;
-- }
--
-- //! Set the time for the time dependent boundary conditions (called from main)
-- void setTime(Scalar time)
-- { time_ = time; }
--
-- //! Return true if the injection is currently active
-- bool injectionActive() const
-- { return time_ < injectionDuration_; }
--
-- //! Return true if the given position is in the injection boundary region
-- bool onInjectionBoundary(const GlobalPosition& globalPos) const
-- {
-- return globalPos[1] < 15. + eps_
-- && globalPos[1] > 7. - eps_
-- && globalPos[0] > this->gridGeometry().bBoxMax()[0] - eps_;
-- }
--
--private:
-- static constexpr Scalar eps_ = 1e-6;
-- std::string name_; //! Problem name
-- Scalar aquiferDepth_; //! Depth of the aquifer in m
-- Scalar injectionDuration_; //! Duration of the injection in seconds
-- Scalar time_;
-- //TODO: dumux-course-task
-- //Define the Scalar totalAreaSpecificInflow_ here
--
--};
--
--} //end namespace Dumux
--
--#endif
diff -ruN exercises/exercise-basic/injection2pniproblem.hh exercises/solution/exercise-basic/injection2pniproblem.hh
---- exercises/exercise-basic/injection2pniproblem.hh 2024-07-17 14:37:31.221955214 +0200
-+++ exercises/solution/exercise-basic/injection2pniproblem.hh 2024-07-17 14:37:31.221955214 +0200
+--- exercises/exercise-basic/injection2pniproblem.hh 2024-07-17 13:58:18.208850308 +0200
++++ exercises/solution/exercise-basic/injection2pniproblem.hh 2024-07-17 13:58:18.212850366 +0200
@@ -7,7 +7,7 @@
/*!
* \file
@@ -589,18 +356,14 @@ diff -ruN exercises/exercise-basic/injection2pniproblem.hh exercises/solution/ex
}
/*!
-@@ -142,19 +131,20 @@
- NumEqVector values(0.0);
-
- // If we are inside the injection zone set inflow Neumann boundary conditions
-- if (injectionActive() && onInjectionBoundary(globalPos))
-+ if (injectionActive() && onInjectionBoundary(globalPos))
+@@ -144,17 +133,18 @@
+ // if we are inside the injection zone set inflow Neumann boundary conditions
+ if (injectionActive() && onInjectionBoundary(globalPos))
{
-- // Inject nitrogen. negative values mean injection
+ const Scalar injectionRate = -1e-4;
+
-+ // Inject nitrogen. Negative values mean injection
- // Unit: kg/(s*m^2)
+ // inject nitrogen. Negative values mean injection
+ // unit: kg/(s*m^2)
- values[Indices::conti0EqIdx + FluidSystem::N2Idx] = -1e-4;
+ values[Indices::conti0EqIdx + FluidSystem::N2Idx] = injectionRate;
values[Indices::conti0EqIdx + FluidSystem::H2OIdx] = 0.0;
@@ -608,10 +371,10 @@ diff -ruN exercises/exercise-basic/injection2pniproblem.hh exercises/solution/ex
- /*!
- * TODO:dumux-course-task 4:
- * Set Neumann noflow conditions for the energy equation everywhere else except the left boundary.
-- * Additionally, consider the energy flux at the injection point which is equal to the product of the respective mass flux and the matching enthalpy. Use the function `gasEnthalpy(temperature,pressure)` from the N2 component to access the necessary enthalpy.
+- * Additionally, consider the energy flux at the injection point which is equal to the product of the respective mass flux and the matching enthalpy. Use the function `gasEnthalpy(temperature,pressure)` from the N2 component to access the necessary enthalpy.
- * hint: use `Indices::energyEqIdx` to access the entry belonging to the energy flux.
- */
-+ // Energy fluxes are always mass specific
++ // energy fluxes are always mass specific
+ // unit: W/(m^2)
+ const Scalar temperatureAtInjection = initialAtPos(globalPos)[Indices::temperatureIdx];/*K*/
+ const Scalar pressureAtInjection = initialAtPos(globalPos)[Indices::pressureIdx];/*Pa*/
@@ -619,15 +382,6 @@ diff -ruN exercises/exercise-basic/injection2pniproblem.hh exercises/solution/ex
}
return values;
-@@ -180,7 +170,7 @@
- {
- PrimaryVariables values(0.0);
-
-- // get the water density at atmospheric conditions
-+ // Get the water density at atmospheric conditions
- const Scalar densityW = FluidSystem::H2O::liquidDensity(283.15, 1.0e5);
-
- // Assume an initially hydrostatic liquid pressure profile
@@ -190,13 +180,10 @@
values[Indices::pressureIdx] = pw;
values[Indices::saturationIdx] = 0.0;
@@ -647,7 +401,7 @@ diff -ruN exercises/exercise-basic/injection2pniproblem.hh exercises/solution/ex
}
diff -ruN exercises/exercise-basic/injection2pproblem.hh exercises/solution/exercise-basic/injection2pproblem.hh
---- exercises/exercise-basic/injection2pproblem.hh 2024-07-17 14:37:31.221955214 +0200
+--- exercises/exercise-basic/injection2pproblem.hh 2024-07-17 13:58:18.208850308 +0200
+++ exercises/solution/exercise-basic/injection2pproblem.hh 1970-01-01 01:00:00.000000000 +0100
@@ -1,211 +0,0 @@
-// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
@@ -752,10 +506,10 @@ diff -ruN exercises/exercise-basic/injection2pproblem.hh exercises/solution/exer
- BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const
- {
- BoundaryTypes bcTypes;
-- // Set the left of the domain (with the global position in "0 = x" direction as a Dirichlet boundary
+- // set the left of the domain (with the global position in "0 = x" direction as a Dirichlet boundary
- if (globalPos[0] < eps_)
- bcTypes.setAllDirichlet();
-- // Set all other as Neumann boundaries
+- // set all other as Neumann boundaries
- else
- bcTypes.setAllNeumann();
-
@@ -781,10 +535,10 @@ diff -ruN exercises/exercise-basic/injection2pproblem.hh exercises/solution/exer
- */
- NumEqVector neumannAtPos(const GlobalPosition &globalPos) const
- {
-- // Initialize values to zero, i.e. no-flow Neumann boundary conditions
+- // initialize values to zero, i.e. no-flow Neumann boundary conditions
- NumEqVector values(0.0);
-
-- // If we are inside the injection zone set inflow Neumann boundary conditions
+- // if we are inside the injection zone set inflow Neumann boundary conditions
- // using < boundary + eps_ or > boundary - eps_ is safer for floating point comparisons
- // than using <= or >= as it is robust with regard to imprecision introduced by rounding errors.
- if (injectionActive() && onInjectionBoundary(globalPos))
@@ -819,11 +573,11 @@ diff -ruN exercises/exercise-basic/injection2pproblem.hh exercises/solution/exer
- {
- PrimaryVariables values(0.0);
-
-- // Get the water density at atmospheric conditions
+- // get the water density at atmospheric conditions
- const Scalar densityW = FluidSystem::H2O::liquidDensity(this->spatialParams().temperatureAtPos(globalPos), 1.0e5);
-
-- // Assume an initially hydrostatic liquid pressure profile
-- // Note: we subtract rho_w*g*h because g is defined negative
+- // assume an initially hydrostatic liquid pressure profile
+- // note: we subtract rho_w*g*h because g is defined negative
- const Scalar pw = 1.0e5 - densityW*this->spatialParams().gravity(globalPos)[dimWorld-1]*(aquiferDepth_ - globalPos[dimWorld-1]);
-
- values[Indices::pressureIdx] = pw;
@@ -862,8 +616,8 @@ diff -ruN exercises/exercise-basic/injection2pproblem.hh exercises/solution/exer
-
-#endif
diff -ruN exercises/exercise-basic/params.input exercises/solution/exercise-basic/params.input
---- exercises/exercise-basic/params.input 2024-06-26 11:12:30.497879893 +0200
-+++ exercises/solution/exercise-basic/params.input 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-basic/params.input 2023-10-13 18:08:54.710521000 +0200
++++ exercises/solution/exercise-basic/params.input 2024-07-17 10:58:43.905583103 +0200
@@ -24,7 +24,7 @@
Aquifer.Snr = 0.0
@@ -877,7 +631,7 @@ diff -ruN exercises/exercise-basic/params.input exercises/solution/exercise-basi
+SolidThermalConductivity = 2.8 # solid thermal conducitivity of granite
+SolidHeatCapacity = 790 # solid heat capacity of granite
diff -ruN exercises/exercise-basic/properties2p.hh exercises/solution/exercise-basic/properties2p.hh
---- exercises/exercise-basic/properties2p.hh 2024-07-17 14:02:52.088041093 +0200
+--- exercises/exercise-basic/properties2p.hh 2024-07-17 13:58:18.208850308 +0200
+++ exercises/solution/exercise-basic/properties2p.hh 1970-01-01 01:00:00.000000000 +0100
@@ -1,63 +0,0 @@
-// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
@@ -944,8 +698,8 @@ diff -ruN exercises/exercise-basic/properties2p.hh exercises/solution/exercise-b
-
-#endif
diff -ruN exercises/exercise-basic/properties2pni.hh exercises/solution/exercise-basic/properties2pni.hh
---- exercises/exercise-basic/properties2pni.hh 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-basic/properties2pni.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-basic/properties2pni.hh 2024-07-17 13:58:18.208850308 +0200
++++ exercises/solution/exercise-basic/properties2pni.hh 2024-07-17 13:58:18.212850366 +0200
@@ -24,13 +24,9 @@
namespace Dumux::Properties {
@@ -962,7 +716,7 @@ diff -ruN exercises/exercise-basic/properties2pni.hh exercises/solution/exercise
} // end namespace TTag
diff -ruN exercises/exercise-basic/README.md exercises/solution/exercise-basic/README.md
---- exercises/exercise-basic/README.md 2024-07-17 14:37:31.109952737 +0200
+--- exercises/exercise-basic/README.md 2024-07-17 14:11:07.348102188 +0200
+++ exercises/solution/exercise-basic/README.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,94 +0,0 @@
-# Exercise Basics (DuMuX course)
diff --git a/.patches/exercise-biomineralization/exercise-biomineralization.patch b/.patches/exercise-biomineralization/exercise-biomineralization.patch
index 1f7592f64e9a418b046717bb3637a959bcb5dc17..fb7d5f82bf2d4648ba12177094ef980e789db685 100644
--- a/.patches/exercise-biomineralization/exercise-biomineralization.patch
+++ b/.patches/exercise-biomineralization/exercise-biomineralization.patch
@@ -1,6 +1,6 @@
diff -ruN exercises/exercise-biomineralization/biominproblem.hh exercises/solution/exercise-biomineralization/biominproblem.hh
---- exercises/exercise-biomineralization/biominproblem.hh 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-biomineralization/biominproblem.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-biomineralization/biominproblem.hh 2024-07-17 11:27:28.870349025 +0200
++++ exercises/solution/exercise-biomineralization/biominproblem.hh 2024-07-17 11:27:28.878349199 +0200
@@ -19,10 +19,12 @@
// TODO: dumux-course-task
@@ -70,8 +70,8 @@ diff -ruN exercises/exercise-biomineralization/biominproblem.hh exercises/soluti
{
return Kxx_;
diff -ruN exercises/exercise-biomineralization/biominspatialparams.hh exercises/solution/exercise-biomineralization/biominspatialparams.hh
---- exercises/exercise-biomineralization/biominspatialparams.hh 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-biomineralization/biominspatialparams.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-biomineralization/biominspatialparams.hh 2024-07-17 11:27:28.870349025 +0200
++++ exercises/solution/exercise-biomineralization/biominspatialparams.hh 2024-07-17 11:27:28.878349199 +0200
@@ -18,7 +18,7 @@
#include <dumux/material/fluidmatrixinteractions/porosityprecipitation.hh>
// TODO: dumux-course-task
@@ -91,8 +91,8 @@ diff -ruN exercises/exercise-biomineralization/biominspatialparams.hh exercises/
Scalar initialPorosity_;
diff -ruN exercises/exercise-biomineralization/chemistry/simplebiominreactions.hh exercises/solution/exercise-biomineralization/chemistry/simplebiominreactions.hh
---- exercises/exercise-biomineralization/chemistry/simplebiominreactions.hh 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-biomineralization/chemistry/simplebiominreactions.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-biomineralization/chemistry/simplebiominreactions.hh 2024-07-17 11:27:28.870349025 +0200
++++ exercises/solution/exercise-biomineralization/chemistry/simplebiominreactions.hh 2024-07-17 11:27:28.878349199 +0200
@@ -73,24 +73,26 @@
// define and compute some parameters for convenience:
const Scalar xwCa = volVars.moleFraction(liquidPhaseIdx,CaIdx);
@@ -143,8 +143,8 @@ diff -ruN exercises/exercise-biomineralization/chemistry/simplebiominreactions.h
private:
diff -ruN exercises/exercise-biomineralization/CMakeLists.txt exercises/solution/exercise-biomineralization/CMakeLists.txt
---- exercises/exercise-biomineralization/CMakeLists.txt 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-biomineralization/CMakeLists.txt 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-biomineralization/CMakeLists.txt 2024-07-17 11:27:28.798347462 +0200
++++ exercises/solution/exercise-biomineralization/CMakeLists.txt 2024-07-17 11:27:28.798347462 +0200
@@ -2,9 +2,9 @@
# SPDX-License-Identifier: GPL-3.0-or-later
@@ -158,8 +158,8 @@ diff -ruN exercises/exercise-biomineralization/CMakeLists.txt exercises/solution
# add a symlink for each input file
add_input_file_links()
diff -ruN exercises/exercise-biomineralization/components/biofilm.hh exercises/solution/exercise-biomineralization/components/biofilm.hh
---- exercises/exercise-biomineralization/components/biofilm.hh 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-biomineralization/components/biofilm.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-biomineralization/components/biofilm.hh 2024-07-17 11:27:28.870349025 +0200
++++ exercises/solution/exercise-biomineralization/components/biofilm.hh 2024-07-17 11:27:28.878349199 +0200
@@ -12,9 +12,10 @@
#ifndef DUMUX_BIOFILM_HH
#define DUMUX_BIOFILM_HH
@@ -173,8 +173,8 @@ diff -ruN exercises/exercise-biomineralization/components/biofilm.hh exercises/s
namespace Dumux::Components {
diff -ruN exercises/exercise-biomineralization/fluidmatrixinteractions/permeabilitypowerlaw.hh exercises/solution/exercise-biomineralization/fluidmatrixinteractions/permeabilitypowerlaw.hh
---- exercises/exercise-biomineralization/fluidmatrixinteractions/permeabilitypowerlaw.hh 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-biomineralization/fluidmatrixinteractions/permeabilitypowerlaw.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-biomineralization/fluidmatrixinteractions/permeabilitypowerlaw.hh 2024-07-17 11:27:28.870349025 +0200
++++ exercises/solution/exercise-biomineralization/fluidmatrixinteractions/permeabilitypowerlaw.hh 2024-07-17 11:27:28.878349199 +0200
@@ -41,10 +41,11 @@
using std::pow;
// TODO: dumux-course-task
@@ -189,8 +189,8 @@ diff -ruN exercises/exercise-biomineralization/fluidmatrixinteractions/permeabil
};
diff -ruN exercises/exercise-biomineralization/fluidsystems/biomin.hh exercises/solution/exercise-biomineralization/fluidsystems/biomin.hh
---- exercises/exercise-biomineralization/fluidsystems/biomin.hh 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-biomineralization/fluidsystems/biomin.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-biomineralization/fluidsystems/biomin.hh 2024-07-17 11:27:28.870349025 +0200
++++ exercises/solution/exercise-biomineralization/fluidsystems/biomin.hh 2024-07-17 11:27:28.878349199 +0200
@@ -219,7 +219,7 @@
static void init()
{
@@ -201,8 +201,8 @@ diff -ruN exercises/exercise-biomineralization/fluidsystems/biomin.hh exercises/
}
diff -ruN exercises/exercise-biomineralization/params.input exercises/solution/exercise-biomineralization/params.input
---- exercises/exercise-biomineralization/params.input 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-biomineralization/params.input 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-biomineralization/params.input 2024-07-17 10:58:43.893583024 +0200
++++ exercises/solution/exercise-biomineralization/params.input 2024-07-17 10:58:43.905583103 +0200
@@ -16,10 +16,10 @@
InitBiofilm = 0.05 # [-] initial volumefraction biofilm
@@ -234,8 +234,8 @@ diff -ruN exercises/exercise-biomineralization/params.input exercises/solution/e
+[PowerLaw]
+Exponent = 5.0
diff -ruN exercises/exercise-biomineralization/properties.hh exercises/solution/exercise-biomineralization/properties.hh
---- exercises/exercise-biomineralization/properties.hh 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-biomineralization/properties.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-biomineralization/properties.hh 2024-07-17 11:27:28.870349025 +0200
++++ exercises/solution/exercise-biomineralization/properties.hh 2024-07-17 11:27:28.878349199 +0200
@@ -18,17 +18,20 @@
#include <dumux/discretization/cctpfa.hh>
#include <dumux/porousmediumflow/2pncmin/model.hh>
@@ -269,7 +269,7 @@ diff -ruN exercises/exercise-biomineralization/properties.hh exercises/solution/
public:
using type = FluidSystems::BioMin<Scalar, CO2Impl, H2OType>;
diff -ruN exercises/exercise-biomineralization/README.md exercises/solution/exercise-biomineralization/README.md
---- exercises/exercise-biomineralization/README.md 2024-07-17 14:02:52.088041093 +0200
+--- exercises/exercise-biomineralization/README.md 2024-07-17 10:58:43.889582997 +0200
+++ exercises/solution/exercise-biomineralization/README.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,257 +0,0 @@
-# Exercise Biomineralization (DuMuX Course)
@@ -530,8 +530,8 @@ diff -ruN exercises/exercise-biomineralization/README.md exercises/solution/exer
-output.CellData.append((S_gas_0-S_gas_1),'diffS_gas');
-```
diff -ruN exercises/exercise-biomineralization/solidsystems/biominsolidphase.hh exercises/solution/exercise-biomineralization/solidsystems/biominsolidphase.hh
---- exercises/exercise-biomineralization/solidsystems/biominsolidphase.hh 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-biomineralization/solidsystems/biominsolidphase.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-biomineralization/solidsystems/biominsolidphase.hh 2024-07-17 11:27:28.870349025 +0200
++++ exercises/solution/exercise-biomineralization/solidsystems/biominsolidphase.hh 2024-07-17 11:27:28.878349199 +0200
@@ -15,12 +15,12 @@
#include <string>
#include <dune/common/exceptions.hh>
diff --git a/.patches/exercise-coupling-ff-pm/exercise-coupling-ff-pm.patch b/.patches/exercise-coupling-ff-pm/exercise-coupling-ff-pm.patch
index 34e91e4ceff1033f2102bbe948bd67628ac3d43a..3152274794c413d45b0a8dc17e6f1d2a39b8f484 100644
--- a/.patches/exercise-coupling-ff-pm/exercise-coupling-ff-pm.patch
+++ b/.patches/exercise-coupling-ff-pm/exercise-coupling-ff-pm.patch
@@ -1,6 +1,6 @@
diff -ruN exercises/exercise-coupling-ff-pm/interface/CMakeLists.txt exercises/solution/exercise-coupling-ff-pm/interface/CMakeLists.txt
---- exercises/exercise-coupling-ff-pm/interface/CMakeLists.txt 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/interface/CMakeLists.txt 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-coupling-ff-pm/interface/CMakeLists.txt 2024-07-17 13:58:18.208850308 +0200
++++ exercises/solution/exercise-coupling-ff-pm/interface/CMakeLists.txt 2024-07-17 13:58:18.216850425 +0200
@@ -1,10 +1,30 @@
# SPDX-FileCopyrightInfo: Copyright © DuMux-Course contributors, see AUTHORS.md in root folder
# SPDX-License-Identifier: GPL-3.0-or-later
@@ -36,8 +36,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/interface/CMakeLists.txt exercises/s
# add a symlink for each input file
add_input_file_links()
diff -ruN exercises/exercise-coupling-ff-pm/interface/freeflowsubproblem.hh exercises/solution/exercise-coupling-ff-pm/interface/freeflowsubproblem.hh
---- exercises/exercise-coupling-ff-pm/interface/freeflowsubproblem.hh 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/interface/freeflowsubproblem.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-coupling-ff-pm/interface/freeflowsubproblem.hh 2024-07-17 13:58:18.208850308 +0200
++++ exercises/solution/exercise-coupling-ff-pm/interface/freeflowsubproblem.hh 2024-07-17 13:58:18.216850425 +0200
@@ -71,8 +71,7 @@
const auto& globalPos = scvf.dofPosition();
@@ -135,8 +135,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/interface/freeflowsubproblem.hh exer
}
diff -ruN exercises/exercise-coupling-ff-pm/interface/main.cc exercises/solution/exercise-coupling-ff-pm/interface/main.cc
---- exercises/exercise-coupling-ff-pm/interface/main.cc 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/interface/main.cc 2024-07-17 14:37:31.113952826 +0200
+--- exercises/exercise-coupling-ff-pm/interface/main.cc 2024-07-17 13:58:18.208850308 +0200
++++ exercises/solution/exercise-coupling-ff-pm/interface/main.cc 2024-07-17 13:58:18.252850955 +0200
@@ -52,9 +52,7 @@
using FreeflowTypeTag = Properties::TTag::FreeflowOneP;
using PorousMediumTypeTag = Properties::TTag::PorousMediumFlowOneP;
@@ -263,8 +263,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/interface/main.cc exercises/solution
using PorousMediumSolutionVector = GetPropType<PorousMediumTypeTag, Properties::SolutionVector>;
VtkOutputModule<PorousMediumGridVariables, PorousMediumSolutionVector> porousMediumVtkWriter(*porousMediumGridVariables,
diff -ruN exercises/exercise-coupling-ff-pm/interface/porousmediumsubproblem.hh exercises/solution/exercise-coupling-ff-pm/interface/porousmediumsubproblem.hh
---- exercises/exercise-coupling-ff-pm/interface/porousmediumsubproblem.hh 2024-07-17 14:37:31.109952737 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/interface/porousmediumsubproblem.hh 2024-07-17 14:37:31.113952826 +0200
+--- exercises/exercise-coupling-ff-pm/interface/porousmediumsubproblem.hh 2024-07-17 13:58:18.252850955 +0200
++++ exercises/solution/exercise-coupling-ff-pm/interface/porousmediumsubproblem.hh 2024-07-17 13:58:18.252850955 +0200
@@ -68,9 +68,10 @@
// set Neumann BCs to all boundaries first
values.setAllNeumann();
@@ -278,8 +278,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/interface/porousmediumsubproblem.hh
// set the coupling boundary condition at the interface
if (couplingManager().isCoupledEntity(CouplingManager::darcyIdx, scvf))
diff -ruN exercises/exercise-coupling-ff-pm/interface/properties.hh exercises/solution/exercise-coupling-ff-pm/interface/properties.hh
---- exercises/exercise-coupling-ff-pm/interface/properties.hh 2024-07-17 14:37:31.109952737 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/interface/properties.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-coupling-ff-pm/interface/properties.hh 2024-07-17 13:58:18.252850955 +0200
++++ exercises/solution/exercise-coupling-ff-pm/interface/properties.hh 2024-07-17 13:58:18.216850425 +0200
@@ -16,9 +16,9 @@
#include <dumux/multidomain/staggeredtraits.hh>
#include <dumux/multidomain/boundary/stokesdarcy/couplingmanager.hh>
@@ -336,7 +336,7 @@ diff -ruN exercises/exercise-coupling-ff-pm/interface/properties.hh exercises/so
template<class TypeTag>
diff -ruN exercises/exercise-coupling-ff-pm/interface/readme.md exercises/solution/exercise-coupling-ff-pm/interface/readme.md
---- exercises/exercise-coupling-ff-pm/interface/readme.md 2024-07-17 14:37:31.109952737 +0200
+--- exercises/exercise-coupling-ff-pm/interface/readme.md 2024-07-17 14:00:29.510779420 +0200
+++ exercises/solution/exercise-coupling-ff-pm/interface/readme.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,169 +0,0 @@
-## 1. Changing the interface
@@ -509,8 +509,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/interface/readme.md exercises/soluti
-What changes to the left boundary conditions in the free-flow domain would you make to introduce this? What conditions can be enforced on the right boundary?
-Hint: A relation between velocity and position is used for the vertical velocity component in the original form of the `dirichletAtPos` method.
diff -ruN exercises/exercise-coupling-ff-pm/models/CMakeLists.txt exercises/solution/exercise-coupling-ff-pm/models/CMakeLists.txt
---- exercises/exercise-coupling-ff-pm/models/CMakeLists.txt 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/models/CMakeLists.txt 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-coupling-ff-pm/models/CMakeLists.txt 2024-07-17 13:58:18.208850308 +0200
++++ exercises/solution/exercise-coupling-ff-pm/models/CMakeLists.txt 2024-07-17 13:58:18.216850425 +0200
@@ -1,9 +1,32 @@
# SPDX-FileCopyrightInfo: Copyright © DuMux-Course contributors, see AUTHORS.md in root folder
# SPDX-License-Identifier: GPL-3.0-or-later
@@ -549,8 +549,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/models/CMakeLists.txt exercises/solu
+
+dune_symlink_to_source_files(FILES "params_orig_a.input" "params_b_c.input" plotFluxes.py)
diff -ruN exercises/exercise-coupling-ff-pm/models/main.cc exercises/solution/exercise-coupling-ff-pm/models/main.cc
---- exercises/exercise-coupling-ff-pm/models/main.cc 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/models/main.cc 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-coupling-ff-pm/models/main.cc 2024-07-17 13:58:18.208850308 +0200
++++ exercises/solution/exercise-coupling-ff-pm/models/main.cc 2024-07-17 13:58:18.216850425 +0200
@@ -127,9 +127,15 @@
auto porousMediumGridVariables = std::make_shared<PorousMediumGridVariables>(porousMediumProblem, porousMediumFvGridGeometry);
porousMediumGridVariables->init(sol[porousMediumIdx]);
@@ -572,7 +572,7 @@ diff -ruN exercises/exercise-coupling-ff-pm/models/main.cc exercises/solution/ex
GetPropType<FreeflowTypeTag, Properties::IOFields>::initOutputModule(freeflowVtkWriter);
diff -ruN exercises/exercise-coupling-ff-pm/models/params_b_c.input exercises/solution/exercise-coupling-ff-pm/models/params_b_c.input
--- exercises/exercise-coupling-ff-pm/models/params_b_c.input 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-coupling-ff-pm/models/params_b_c.input 2024-07-16 15:45:21.113991695 +0200
++++ exercises/solution/exercise-coupling-ff-pm/models/params_b_c.input 2024-07-16 16:28:05.559708538 +0200
@@ -0,0 +1,57 @@
+[TimeLoop]
+DtInitial = 100 # s
@@ -632,7 +632,7 @@ diff -ruN exercises/exercise-coupling-ff-pm/models/params_b_c.input exercises/so
+[Assembly]
+NumericDifference.BaseEpsilon = 1e-8
diff -ruN exercises/exercise-coupling-ff-pm/models/params.input exercises/solution/exercise-coupling-ff-pm/models/params.input
---- exercises/exercise-coupling-ff-pm/models/params.input 2024-07-17 14:37:31.109952737 +0200
+--- exercises/exercise-coupling-ff-pm/models/params.input 2024-07-17 13:58:18.252850955 +0200
+++ exercises/solution/exercise-coupling-ff-pm/models/params.input 1970-01-01 01:00:00.000000000 +0100
@@ -1,58 +0,0 @@
-[TimeLoop]
@@ -695,7 +695,7 @@ diff -ruN exercises/exercise-coupling-ff-pm/models/params.input exercises/soluti
-NumericDifference.BaseEpsilon = 1e-8
diff -ruN exercises/exercise-coupling-ff-pm/models/params_orig_a.input exercises/solution/exercise-coupling-ff-pm/models/params_orig_a.input
--- exercises/exercise-coupling-ff-pm/models/params_orig_a.input 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-coupling-ff-pm/models/params_orig_a.input 2024-07-16 15:45:21.113991695 +0200
++++ exercises/solution/exercise-coupling-ff-pm/models/params_orig_a.input 2024-07-16 16:28:05.563708580 +0200
@@ -0,0 +1,56 @@
+[TimeLoop]
+DtInitial = 100 # s
@@ -754,8 +754,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/models/params_orig_a.input exercises
+[Assembly]
+NumericDifference.BaseEpsilon = 1e-8
diff -ruN exercises/exercise-coupling-ff-pm/models/porousmediumsubproblem.hh exercises/solution/exercise-coupling-ff-pm/models/porousmediumsubproblem.hh
---- exercises/exercise-coupling-ff-pm/models/porousmediumsubproblem.hh 2024-07-17 14:37:31.109952737 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/models/porousmediumsubproblem.hh 2024-07-17 14:37:31.113952826 +0200
+--- exercises/exercise-coupling-ff-pm/models/porousmediumsubproblem.hh 2024-07-17 13:58:18.252850955 +0200
++++ exercises/solution/exercise-coupling-ff-pm/models/porousmediumsubproblem.hh 2024-07-17 13:58:18.252850955 +0200
@@ -52,10 +52,15 @@
// primary variable indices
static constexpr int conti0EqIdx = Indices::conti0EqIdx;
@@ -877,8 +877,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/models/porousmediumsubproblem.hh exe
std::shared_ptr<CouplingManager> couplingManager_;
diff -ruN exercises/exercise-coupling-ff-pm/models/properties.hh exercises/solution/exercise-coupling-ff-pm/models/properties.hh
---- exercises/exercise-coupling-ff-pm/models/properties.hh 2024-07-17 14:37:31.109952737 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/models/properties.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-coupling-ff-pm/models/properties.hh 2024-07-17 13:58:18.252850955 +0200
++++ exercises/solution/exercise-coupling-ff-pm/models/properties.hh 2024-07-17 13:58:18.216850425 +0200
@@ -25,12 +25,14 @@
#include <dumux/discretization/cctpfa.hh>
#include <dumux/material/fluidmatrixinteractions/diffusivityconstanttortuosity.hh>
@@ -961,7 +961,7 @@ diff -ruN exercises/exercise-coupling-ff-pm/models/properties.hh exercises/solut
template<class TypeTag>
struct EnableGridGeometryCache<TypeTag, TTag::FreeflowNC> { static constexpr bool value = true; };
diff -ruN exercises/exercise-coupling-ff-pm/models/readme.md exercises/solution/exercise-coupling-ff-pm/models/readme.md
---- exercises/exercise-coupling-ff-pm/models/readme.md 2024-07-17 14:37:31.109952737 +0200
+--- exercises/exercise-coupling-ff-pm/models/readme.md 2024-07-17 14:03:33.957481940 +0200
+++ exercises/solution/exercise-coupling-ff-pm/models/readme.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,118 +0,0 @@
-## 2. Changing the porous medium model
@@ -1083,7 +1083,7 @@ diff -ruN exercises/exercise-coupling-ff-pm/models/readme.md exercises/solution/
-
-Now you are able to simulate a complete drying of the porous medium. Have a look the resulting liquid saturation distribution within the porous medium (with using paraview).
diff -ruN exercises/exercise-coupling-ff-pm/README.md exercises/solution/exercise-coupling-ff-pm/README.md
---- exercises/exercise-coupling-ff-pm/README.md 2024-07-17 14:37:31.109952737 +0200
+--- exercises/exercise-coupling-ff-pm/README.md 2024-07-17 14:00:29.510779420 +0200
+++ exercises/solution/exercise-coupling-ff-pm/README.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,125 +0,0 @@
-# Exercise Coupling free flow/porous medium flow (DuMuX Course)
@@ -1212,8 +1212,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/README.md exercises/solution/exercis
-* [**Exercise 2:** Changing the porous medium model](./models/readme.md)
-* [**Exercise 3:** Introducing a turbulence model in the free flow domain](./turbulence/readme.md)
diff -ruN exercises/exercise-coupling-ff-pm/turbulence/CMakeLists.txt exercises/solution/exercise-coupling-ff-pm/turbulence/CMakeLists.txt
---- exercises/exercise-coupling-ff-pm/turbulence/CMakeLists.txt 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/turbulence/CMakeLists.txt 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-coupling-ff-pm/turbulence/CMakeLists.txt 2024-07-17 13:58:18.208850308 +0200
++++ exercises/solution/exercise-coupling-ff-pm/turbulence/CMakeLists.txt 2024-07-17 13:58:18.216850425 +0200
@@ -1,10 +1,33 @@
# SPDX-FileCopyrightInfo: Copyright © DuMux-Course contributors, see AUTHORS.md in root folder
# SPDX-License-Identifier: GPL-3.0-or-later
@@ -1252,8 +1252,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/turbulence/CMakeLists.txt exercises/
# add a symlink for each input file
add_input_file_links()
diff -ruN exercises/exercise-coupling-ff-pm/turbulence/freeflowsubproblem.hh exercises/solution/exercise-coupling-ff-pm/turbulence/freeflowsubproblem.hh
---- exercises/exercise-coupling-ff-pm/turbulence/freeflowsubproblem.hh 2024-07-17 14:37:31.113952826 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/turbulence/freeflowsubproblem.hh 2024-07-17 14:37:31.113952826 +0200
+--- exercises/exercise-coupling-ff-pm/turbulence/freeflowsubproblem.hh 2024-07-17 13:58:18.252850955 +0200
++++ exercises/solution/exercise-coupling-ff-pm/turbulence/freeflowsubproblem.hh 2024-07-17 13:58:18.252850955 +0200
@@ -8,8 +8,8 @@
* \file
* \brief The free-flow sub problem
@@ -1469,8 +1469,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/turbulence/freeflowsubproblem.hh exe
return values;
}
diff -ruN exercises/exercise-coupling-ff-pm/turbulence/main.cc exercises/solution/exercise-coupling-ff-pm/turbulence/main.cc
---- exercises/exercise-coupling-ff-pm/turbulence/main.cc 2024-07-17 14:02:52.088041093 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/turbulence/main.cc 2024-07-17 14:37:31.113952826 +0200
+--- exercises/exercise-coupling-ff-pm/turbulence/main.cc 2024-07-17 13:58:18.208850308 +0200
++++ exercises/solution/exercise-coupling-ff-pm/turbulence/main.cc 2024-07-17 13:58:18.256851014 +0200
@@ -122,11 +122,12 @@
couplingManager->init(freeflowProblem, porousMediumProblem, sol);
@@ -1501,7 +1501,7 @@ diff -ruN exercises/exercise-coupling-ff-pm/turbulence/main.cc exercises/solutio
porousMediumProblem->postTimeStep(sol[porousMediumIdx], *porousMediumGridVariables, timeLoop->timeStepSize());
diff -ruN exercises/exercise-coupling-ff-pm/turbulence/params_b.input exercises/solution/exercise-coupling-ff-pm/turbulence/params_b.input
--- exercises/exercise-coupling-ff-pm/turbulence/params_b.input 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-coupling-ff-pm/turbulence/params_b.input 2024-07-16 15:45:21.113991695 +0200
++++ exercises/solution/exercise-coupling-ff-pm/turbulence/params_b.input 2024-07-16 16:28:05.563708580 +0200
@@ -0,0 +1,72 @@
+[TimeLoop]
+DtInitial = 1e-1 # [s]
@@ -1577,7 +1577,7 @@ diff -ruN exercises/exercise-coupling-ff-pm/turbulence/params_b.input exercises/
+IsFlatWallBounded = True
diff -ruN exercises/exercise-coupling-ff-pm/turbulence/params_c_d.input exercises/solution/exercise-coupling-ff-pm/turbulence/params_c_d.input
--- exercises/exercise-coupling-ff-pm/turbulence/params_c_d.input 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-coupling-ff-pm/turbulence/params_c_d.input 2024-07-16 15:45:21.113991695 +0200
++++ exercises/solution/exercise-coupling-ff-pm/turbulence/params_c_d.input 2024-07-16 16:28:05.563708580 +0200
@@ -0,0 +1,72 @@
+[TimeLoop]
+DtInitial = 1e-1 # [s]
@@ -1652,7 +1652,7 @@ diff -ruN exercises/exercise-coupling-ff-pm/turbulence/params_c_d.input exercise
+[RANS]
+IsFlatWallBounded = True
diff -ruN exercises/exercise-coupling-ff-pm/turbulence/params.input exercises/solution/exercise-coupling-ff-pm/turbulence/params.input
---- exercises/exercise-coupling-ff-pm/turbulence/params.input 2024-07-16 15:45:21.113991695 +0200
+--- exercises/exercise-coupling-ff-pm/turbulence/params.input 2024-07-16 16:28:05.559708538 +0200
+++ exercises/solution/exercise-coupling-ff-pm/turbulence/params.input 1970-01-01 01:00:00.000000000 +0100
@@ -1,73 +0,0 @@
-[TimeLoop]
@@ -1730,7 +1730,7 @@ diff -ruN exercises/exercise-coupling-ff-pm/turbulence/params.input exercises/so
-SolidHeatCapacity = 790
diff -ruN exercises/exercise-coupling-ff-pm/turbulence/params_orig_a.input exercises/solution/exercise-coupling-ff-pm/turbulence/params_orig_a.input
--- exercises/exercise-coupling-ff-pm/turbulence/params_orig_a.input 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-coupling-ff-pm/turbulence/params_orig_a.input 2024-07-16 15:45:21.113991695 +0200
++++ exercises/solution/exercise-coupling-ff-pm/turbulence/params_orig_a.input 2024-07-16 16:28:05.563708580 +0200
@@ -0,0 +1,72 @@
+[TimeLoop]
+DtInitial = 1e-1 # [s]
@@ -1805,8 +1805,8 @@ diff -ruN exercises/exercise-coupling-ff-pm/turbulence/params_orig_a.input exerc
+[RANS]
+IsFlatWallBounded = True
diff -ruN exercises/exercise-coupling-ff-pm/turbulence/properties.hh exercises/solution/exercise-coupling-ff-pm/turbulence/properties.hh
---- exercises/exercise-coupling-ff-pm/turbulence/properties.hh 2024-07-17 14:37:31.113952826 +0200
-+++ exercises/solution/exercise-coupling-ff-pm/turbulence/properties.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-coupling-ff-pm/turbulence/properties.hh 2024-07-17 13:58:18.252850955 +0200
++++ exercises/solution/exercise-coupling-ff-pm/turbulence/properties.hh 2024-07-17 13:58:18.216850425 +0200
@@ -8,8 +8,8 @@
* \file
* \brief The coupled exercise properties file or the turbulent case.
@@ -1848,7 +1848,7 @@ diff -ruN exercises/exercise-coupling-ff-pm/turbulence/properties.hh exercises/s
// Set the coupling manager
diff -ruN exercises/exercise-coupling-ff-pm/turbulence/readme.md exercises/solution/exercise-coupling-ff-pm/turbulence/readme.md
---- exercises/exercise-coupling-ff-pm/turbulence/readme.md 2024-07-16 15:45:21.113991695 +0200
+--- exercises/exercise-coupling-ff-pm/turbulence/readme.md 2024-07-16 16:28:05.559708538 +0200
+++ exercises/solution/exercise-coupling-ff-pm/turbulence/readme.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,126 +0,0 @@
-## 3. Use a turbulence model in the free flow domain
diff --git a/.patches/exercise-dunemodule/exercise-dunemodule.patch b/.patches/exercise-dunemodule/exercise-dunemodule.patch
index 088dc39d328030ea0f927fbe8398e68b2b304245..381f901fefa706d93c8ce02439e9c8fa56bd227b 100644
--- a/.patches/exercise-dunemodule/exercise-dunemodule.patch
+++ b/.patches/exercise-dunemodule/exercise-dunemodule.patch
@@ -1,5 +1,5 @@
diff -ruN exercises/exercise-dunemodule/README.md exercises/solution/exercise-dunemodule/README.md
---- exercises/exercise-dunemodule/README.md 2024-07-17 14:37:31.113952826 +0200
+--- exercises/exercise-dunemodule/README.md 2024-07-17 14:00:29.510779420 +0200
+++ exercises/solution/exercise-dunemodule/README.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,114 +0,0 @@
-# Exercise New Dune Module (DuMuX course)
diff --git a/.patches/exercise-fluidsystem/exercise-fluidsystem.patch b/.patches/exercise-fluidsystem/exercise-fluidsystem.patch
index 6af0be58a1b56865220b36bc7f59dbe5696e10d7..db24b46fe18b96f2f852e7d78cf9743b6e7284f6 100644
--- a/.patches/exercise-fluidsystem/exercise-fluidsystem.patch
+++ b/.patches/exercise-fluidsystem/exercise-fluidsystem.patch
@@ -1,59 +1,36 @@
-diff -ruN exercises/exercise-fluidsystem/2p2cproblem.hh exercises/solution/exercise-fluidsystem/2p2cproblem.hh
---- exercises/exercise-fluidsystem/2p2cproblem.hh 2024-07-17 14:37:31.333957691 +0200
-+++ exercises/solution/exercise-fluidsystem/2p2cproblem.hh 2024-07-17 14:37:31.333957691 +0200
-@@ -12,7 +12,7 @@
- #ifndef DUMUX_EXERCISE_FLUIDSYSTEM_B_PROBLEM_HH
- #define DUMUX_EXERCISE_FLUIDSYSTEM_B_PROBLEM_HH
-
--// The base porous media box problem.
-+// The base porous media box problem.
- #include <dumux/porousmediumflow/problem.hh>
- #include <dumux/common/boundarytypes.hh>
- #include <dumux/common/properties.hh>
diff -ruN exercises/exercise-fluidsystem/2pproperties.hh exercises/solution/exercise-fluidsystem/2pproperties.hh
---- exercises/exercise-fluidsystem/2pproperties.hh 2024-07-17 14:37:31.357958221 +0200
-+++ exercises/solution/exercise-fluidsystem/2pproperties.hh 2024-07-17 14:37:31.357958221 +0200
-@@ -25,15 +25,13 @@
+--- exercises/exercise-fluidsystem/2pproperties.hh 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fluidsystem/2pproperties.hh 2024-07-17 11:27:28.882349286 +0200
+@@ -25,10 +25,8 @@
#include "spatialparams.hh"
// The components that will be created in this exercise
--// TODO: dumux-course-task 2.2:
--// Include the corresponding header for the task
+-// TODO: dumux-course-task 2.2
+-//****** Include the corresponding header for the task *****//
#include "components/myincompressiblecomponent.hh"
-// #include "components/mycompressiblecomponent.hh"
+#include "components/mycompressiblecomponent.hh"
// We will only have liquid phases here
#include <dumux/material/fluidsystems/1pliquid.hh>
-
--// The two-phase immiscible fluidsystem
-+// The two-phase immiscible fluid system
- #include <dumux/material/fluidsystems/2pimmiscible.hh>
-
- // The interface to create plots during simulation
-@@ -76,12 +74,13 @@
+@@ -77,10 +75,10 @@
using Scalar = GetPropType<TypeTag, Properties::Scalar>;
using TabulatedH2O = Components::TabulatedComponent<Components::H2O<Scalar>>;
using LiquidWaterPhase = typename FluidSystems::OnePLiquid<Scalar, TabulatedH2O>;
-
-- // TODO: dumux-course-task 2.2:
-- // Select the corresponding component for the task
+- // TODO: dumux-course-task 2.2
+- //****** Select the corresponding component for the task *****//
- // Uncomment first line and comment second line for using the compressible component
-- using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyIncompressibleComponent<Scalar> >;
-- // using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyCompressibleComponent<Scalar> >;
+ /*!
-+ * TODO: dumux-course-task 2.2:
+ * Uncomment first line and comment second line for using the incompressible component
+ * Uncomment second line and comment first line for using the compressible component
+ */
-+ using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyIncompressibleComponent<Scalar> >;
-+ //using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyCompressibleComponent<Scalar> >;
+ using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyIncompressibleComponent<Scalar> >;
+ // using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyCompressibleComponent<Scalar> >;
- public:
- using type = typename FluidSystems::TwoPImmiscible<Scalar, LiquidWaterPhase, LiquidMyComponentPhase>;
diff -ruN exercises/exercise-fluidsystem/aparams.input exercises/solution/exercise-fluidsystem/aparams.input
---- exercises/exercise-fluidsystem/aparams.input 2024-07-17 14:37:31.333957691 +0200
-+++ exercises/solution/exercise-fluidsystem/aparams.input 2024-07-17 14:37:31.333957691 +0200
+--- exercises/exercise-fluidsystem/aparams.input 2023-10-13 18:08:54.714521020 +0200
++++ exercises/solution/exercise-fluidsystem/aparams.input 2020-11-12 12:07:04.129814352 +0100
@@ -3,7 +3,7 @@
DtInitial = 10 # initial time step size [s]
@@ -62,7 +39,7 @@ diff -ruN exercises/exercise-fluidsystem/aparams.input exercises/solution/exerci
+Name = exercise-fluidsystem_a_solution # name will be given to e.g. to the vtk result files
[SpatialParams]
- BrooksCoreyPcEntry = 5.0e2 # [Pa]
+ BrooksCoreyPcEntry = 5.0e2 # Pa
@@ -20,4 +20,4 @@
Cells = 60 60 # x-/y-resolution of the grid
@@ -70,8 +47,8 @@ diff -ruN exercises/exercise-fluidsystem/aparams.input exercises/solution/exerci
-PlotDensity = false # plot density over pressure for your component
+PlotDensity = true # plot density over pressure for your component
diff -ruN exercises/exercise-fluidsystem/bparams.input exercises/solution/exercise-fluidsystem/bparams.input
---- exercises/exercise-fluidsystem/bparams.input 2024-07-17 14:37:31.333957691 +0200
-+++ exercises/solution/exercise-fluidsystem/bparams.input 2024-07-17 14:37:31.333957691 +0200
+--- exercises/exercise-fluidsystem/bparams.input 2023-10-13 18:08:54.714521020 +0200
++++ exercises/solution/exercise-fluidsystem/bparams.input 2020-11-12 12:07:04.129814352 +0100
@@ -3,7 +3,7 @@
DtInitial = 10 # initial time step size [s]
@@ -80,10 +57,10 @@ diff -ruN exercises/exercise-fluidsystem/bparams.input exercises/solution/exerci
+Name = exercise-fluidsystem_b_solution # name will be given to e.g. to the vtk result files
[SpatialParams]
- BrooksCoreyPcEntry = 5.0e2 # [Pa]
+ BrooksCoreyPcEntry = 5.0e2 # Pa
diff -ruN exercises/exercise-fluidsystem/CMakeLists.txt exercises/solution/exercise-fluidsystem/CMakeLists.txt
---- exercises/exercise-fluidsystem/CMakeLists.txt 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-fluidsystem/CMakeLists.txt 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-fluidsystem/CMakeLists.txt 2024-07-17 11:27:28.798347462 +0200
++++ exercises/solution/exercise-fluidsystem/CMakeLists.txt 2024-07-17 11:27:28.798347462 +0200
@@ -3,18 +3,16 @@
# executables for exercise part a & b
@@ -110,25 +87,25 @@ diff -ruN exercises/exercise-fluidsystem/CMakeLists.txt exercises/solution/exerc
-# add a symlink for the grids folder
-dune_symlink_to_source_files(FILES grids)
diff -ruN exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh exercises/solution/exercise-fluidsystem/components/mycompressiblecomponent.hh
---- exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh 2024-07-17 14:37:31.265956187 +0200
-+++ exercises/solution/exercise-fluidsystem/components/mycompressiblecomponent.hh 2024-07-17 14:37:31.265956187 +0200
+--- exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fluidsystem/components/mycompressiblecomponent.hh 2024-07-17 11:27:28.882349286 +0200
@@ -38,7 +38,7 @@
static std::string name()
{ return "MyCompressibleComponent"; }
- /*!
+ /*!
- * \brief Returns true if the liquid phase is assumed to be compressible.
+ * \brief Returns true if the liquid phase is assumed to be compressible
*/
static constexpr bool liquidIsCompressible()
@@ -49,11 +49,7 @@
*/
static Scalar molarMass()
{
-- // TODO: dumux-course-task 2.2:
+- // TODO: dumux-course-task 2.2
- // Copy the methods implemented in MyIncompressibleComponent and substitute
- // the density calculation by the expression given in the exercise description.
-- // Replace the line below by a meaningful return statement.
+- // Replace the line below by a meaningful return statement
- DUNE_THROW(Dune::NotImplemented, "Todo: implement molar mass");
+ return 131.39e-3; // [kg/mol]
}
@@ -139,11 +116,11 @@ diff -ruN exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh e
static Scalar liquidDensity(Scalar temperature, Scalar pressure)
{
- // TODO: dumux-course-task 2.2:
-- // Replace the line below by a meaningful return statement.
+- // Replace the line below by a meaningful return statement
- DUNE_THROW(Dune::NotImplemented, "Todo: implement liquid density");
-+ static const Scalar rho_min = 1440.0;
-+ static const Scalar rho_max = 1480.0;
-+ static const Scalar k = 5.0e-7;
++ static const Scalar rho_min = 1440;
++ static const Scalar rho_max = 1480;
++ static const Scalar k = 5e-7;
+
+ using std::exp;
+ return rho_min + (rho_max - rho_min)/(1 + rho_min*exp(-1.0*k*(rho_max - rho_min)*pressure)); // [kg/m^3]
@@ -155,7 +132,7 @@ diff -ruN exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh e
static Scalar liquidViscosity(Scalar temperature, Scalar pressure)
{
- // TODO: dumux-course-task 2.2:
-- // Replace the line below by a meaningful return statement.
+- // Replace the line below by a meaningful return statement
- DUNE_THROW(Dune::NotImplemented, "Todo: implement liquid viscosity");
+ return 5.7e-4;// [Pa*s]
}
@@ -165,23 +142,41 @@ diff -ruN exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh e
*/
static Scalar vaporPressure(Scalar t)
{
-- // TODO: dumux-course-task 3:
-- // Replace the line below by a meaningful return statement.
+- // TODO: dumux-course-task 3
+- // Replace the line below by a meaningful return statement
- DUNE_THROW(Dune::NotImplemented, "Todo: implement vapour pressure");
-+ return 3900.0; // [Pa] (at 20°C)
++ return 3900; // [Pa] (at 20C)
}
};
diff -ruN exercises/exercise-fluidsystem/components/myincompressiblecomponent.hh exercises/solution/exercise-fluidsystem/components/myincompressiblecomponent.hh
---- exercises/exercise-fluidsystem/components/myincompressiblecomponent.hh 2024-07-17 14:37:31.265956187 +0200
-+++ exercises/solution/exercise-fluidsystem/components/myincompressiblecomponent.hh 2024-07-17 14:37:31.265956187 +0200
+--- exercises/exercise-fluidsystem/components/myincompressiblecomponent.hh 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fluidsystem/components/myincompressiblecomponent.hh 2024-07-17 11:27:28.882349286 +0200
+@@ -7,7 +7,7 @@
+ /*!
+ * \file
+ * \ingroup Components
+- * \brief A fictitious component to be implemented in exercise-fluidsystem.
++ * \brief A fictitious component to be implemented in exercise-fluidsystem a.
+ */
+ #ifndef DUMUX_MYINCOMPRESSIBLECOMPONENT_HH
+ #define DUMUX_MYINCOMPRESSIBLECOMPONENT_HH
+@@ -21,7 +21,7 @@
+ {
+ /*!
+ * \ingroup Components
+- * \brief A fictitious component to be implemented in exercise-fluidsystem.
++ * \brief A fictitious component to be implemented in exercise-fluidsystem a.
+ *
+ * \tparam Scalar The type used for scalar values
+ */
@@ -48,10 +48,7 @@
*/
static Scalar molarMass()
{
-- // TODO: dumux-course-task 2.1:
+- // TODO: dumux-course-task 2.1
- // Implement the methods for the component data given in the exercise description.
-- // Replace the line below by a meaningful return statement.
+- // Replace the line below by a meaningful return statement
- DUNE_THROW(Dune::NotImplemented, "Todo: implement molarMass()");
+ return 131.39e-3; // [kg/mol]
}
@@ -191,9 +186,9 @@ diff -ruN exercises/exercise-fluidsystem/components/myincompressiblecomponent.hh
*/
static Scalar liquidDensity(Scalar temperature, Scalar pressure)
{
-- // TODO: dumux-course-task 2.1:
+- // TODO: dumux-course-task 2.1
- // Implement the methods for the component data given in the exercise description.
-- // Replace the line below by a meaningful return statement.
+- // Replace the line below by a meaningful return statement
- DUNE_THROW(Dune::NotImplemented, "Todo: implement liquidDensity()");
+ return 1460.0; // [kg/m^3]
}
@@ -203,9 +198,9 @@ diff -ruN exercises/exercise-fluidsystem/components/myincompressiblecomponent.hh
*/
static Scalar liquidViscosity(Scalar temperature, Scalar pressure)
{
-- // TODO: dumux-course-task 2.1:
+- // TODO: dumux-course-task 2.1
- // Implement the methods for the component data given in the exercise description.
-- // Replace the line below by a meaningful return statement.
+- // Replace the line below by a meaningful return statement
- DUNE_THROW(Dune::NotImplemented, "Todo: implement liquidViscosity()");
+ return 5.7e-4;// [Pa*s]
}
@@ -213,7 +208,7 @@ diff -ruN exercises/exercise-fluidsystem/components/myincompressiblecomponent.hh
diff -ruN exercises/exercise-fluidsystem/components/plotdensityfunction.py exercises/solution/exercise-fluidsystem/components/plotdensityfunction.py
--- exercises/exercise-fluidsystem/components/plotdensityfunction.py 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-fluidsystem/components/plotdensityfunction.py 2024-07-17 14:37:31.265956187 +0200
++++ exercises/solution/exercise-fluidsystem/components/plotdensityfunction.py 2024-07-17 11:27:28.842348417 +0200
@@ -0,0 +1,22 @@
+#!usr/bin/env python
+# SPDX-FileCopyrightInfo: Copyright © DuMux-Course contributors, see AUTHORS.md in root folder
@@ -222,37 +217,34 @@ diff -ruN exercises/exercise-fluidsystem/components/plotdensityfunction.py exerc
+import numpy as np
+import matplotlib.pyplot as plt
+
-+# Function to calculate rho dependent on pressure.
-+rho_min = 1440.0;
-+rho_max = 1480.0;
-+k = 5.0e-7;
++# function to calculate rho dependent on pressure
++rho_min = 1440;
++rho_max = 1480;
++k = 5e-7;
+
+def rho(p):
+ return rho_min + (rho_max - rho_min)/(1 + rho_min*np.exp(-1.0*k*(rho_max - rho_min)*p));
+
-+# Sample pressure in range (1e4, 1e7) and compute corresponding densities.
++# sample pressure in range (1e4, 1e7) and compute corresponding densities
+p = np.logspace(4, 7, 100)
+r = rho(p)
+
-+# Plot density vs. pressure.
++# plot density vs. pressure
+plt.semilogx(p, r)
+plt.show()
diff -ruN exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh exercises/solution/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh
---- exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh 2024-07-17 14:37:31.357958221 +0200
-+++ exercises/solution/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh 2024-07-17 14:37:31.421959636 +0200
-@@ -1,4 +1,3 @@
--// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
- // vi: set et ts=4 sw=4 sts=4:
- //
- // SPDX-FileCopyrightInfo: Copyright © DuMux-Course contributors, see AUTHORS.md in root folder
-@@ -20,14 +19,14 @@
+--- exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh 2024-07-17 11:27:28.882349286 +0200
+@@ -19,15 +19,15 @@
+
#include <dumux/material/fluidsystems/base.hh>
- // The fictitious component that was created in exercise-fluidsystem 2.2
+-// the fictitious component that was created in exercise-fluidsystem a
-#include <exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh>
++// the fictitious component that was created in exercise 3a
+#include <exercises/solution/exercise-fluidsystem/components/mycompressiblecomponent.hh>
- // The binary coefficients corresponding to this fluid system
+ // the binary coefficients corresponding to this fluid system
-#include <exercises/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh>
+#include <exercises/solution/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh>
@@ -262,8 +254,8 @@ diff -ruN exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent
+
/*!
* \brief A compositional fluid consisting of two liquid phases,
- * which are water and a fictitious component from exercise-fluidsystem 2.2.
-@@ -39,7 +38,7 @@
+ * which are water and a fictitious component from exercise-fluidsystem a.
+@@ -39,7 +39,7 @@
{
using ThisType = H2OMyCompressibleComponent<Scalar, H2OType>;
using Base = Dumux::FluidSystems::Base<Scalar, ThisType>;
@@ -272,7 +264,7 @@ diff -ruN exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent
public:
using H2O = H2OType;
using MyCompressibleComponent = Dumux::MyCompressibleComponent<Scalar>;
-@@ -109,11 +108,13 @@
+@@ -109,11 +109,13 @@
assert(0 <= phaseIdx && phaseIdx < numPhases);
return true;
}
@@ -287,16 +279,8 @@ diff -ruN exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent
static constexpr bool isIdealGas(int phaseIdx)
{ return H2O::gasIsIdeal() && MyCompressibleComponent::gasIsIdeal(); }
-@@ -181,7 +182,6 @@
- DUNE_THROW(Dune::InvalidStateException, "Invalid component index " << compIdx);
- }
-
--
- /*!
- * \brief Return the molar mass of a component in [kg/mol].
- */
-@@ -212,16 +212,13 @@
- // see: doctoral thesis of Steffen Ochs 2007
+@@ -211,16 +213,13 @@
+ // See: doctoral thesis of Steffen Ochs 2007
// Steam injection into saturated porous media : process analysis including experimental and numerical investigations
// http://elib.uni-stuttgart.de/bitstream/11682/271/1/Diss_Ochs_OPUS.pdf
+ Scalar rholH2O = H2O::liquidDensity(fluidState.temperature(phaseIdx), fluidState.pressure(phaseIdx));
@@ -309,8 +293,8 @@ diff -ruN exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent
- // Scalar x_H2O = fluidState.moleFraction(phase0Idx, H2OIdx);
- // Scalar x_myComp = fluidState.moleFraction(phase0Idx, NAPLIdx);
-
-- // TODO: dumux-course-task 3:
-- // implement the composition-dependent water density from the exercise sheet.
+- // TODO: dumux-course-task 3
+- // Implement the composition-dependent water density from the exercise sheet.
-
- DUNE_THROW(Dune::NotImplemented, "Todo: implement composition-dependent density");
+ // return composition-dependent water phase density
@@ -318,38 +302,21 @@ diff -ruN exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent
}
else
{
-@@ -397,8 +394,9 @@
- }
-
- /*!
-- * \brief Partial pressures in the gas phase, taken from saturation vapor pressures.
-+ * \brief Partial pressures in the gas phase, taken from saturation vapor pressures.
- */
-+
- template <class FluidState>
- static Scalar partialPressureGas(const FluidState &fluidState, int phaseIdx,
- int compIdx)
-@@ -417,6 +415,7 @@
- /*!
- * \brief Inverse vapor pressures, taken from inverse saturation vapor pressures.
- */
-+
- template <class FluidState>
- static Scalar inverseVaporPressureCurve(const FluidState &fluidState,
- int phaseIdx,
diff -ruN exercises/exercise-fluidsystem/main.cc exercises/solution/exercise-fluidsystem/main.cc
---- exercises/exercise-fluidsystem/main.cc 2024-07-17 14:37:31.333957691 +0200
-+++ exercises/solution/exercise-fluidsystem/main.cc 2024-07-17 14:37:31.333957691 +0200
-@@ -36,7 +36,7 @@
- int main(int argc, char** argv)
- {
- using namespace Dumux;
--
-+
- // TODO: dumux-course-task 3:
- // TYPETAG is set in CMakeLists.txt as compile time definition.
- // Alternatively, you could write `using TypeTag = Properties::TTag::ExerciseFluidsystemTwoP;`
-@@ -80,13 +80,12 @@
+--- exercises/exercise-fluidsystem/main.cc 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fluidsystem/main.cc 2024-07-17 11:27:28.882349286 +0200
+@@ -39,8 +39,8 @@
+
+ // define the type tag for this problem
+ // TYPETAG is set in CMakeLists.txt as compile time definition
+- // alternatively you could write `using TypeTag = Properties::TTag::ExerciseFluidsystemBoxTwoP;`
+- // then, for the 2p2c problem you have to change this line to `using TypeTag = Properties::TTag::ExerciseFluidsystemBoxTwoPTwoC;`
++ // alternatively you could write `using TypeTag = Properties::TTag::ExerciseFluidsystemTwoP;`
++ // then, for the 2p2c problem you have to change this line to `using TypeTag = Properties::TTag::ExerciseFluidsystemTwoPTwoC;`
+ // and recompile the executable
+ using TypeTag = Properties::TTag::TYPETAG;
+
+@@ -80,13 +80,11 @@
auto gridVariables = std::make_shared<GridVariables>(problem, gridGeometry);
gridVariables->init(x);
@@ -360,28 +327,29 @@ diff -ruN exercises/exercise-fluidsystem/main.cc exercises/solution/exercise-flu
const auto maxDt = getParam<Scalar>("TimeLoop.MaxTimeStepSize");
auto dt = getParam<Scalar>("TimeLoop.DtInitial");
-
-+
// initialize the vtk output module
using IOFields = GetPropType<TypeTag, Properties::IOFields>;
VtkOutputModule<GridVariables, SolutionVector> vtkWriter(*gridVariables, x, problem->name());
diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-fluidsystem/README.md
---- exercises/exercise-fluidsystem/README.md 2024-07-17 14:37:31.377958664 +0200
+--- exercises/exercise-fluidsystem/README.md 2023-10-13 18:08:54.714521020 +0200
+++ exercises/solution/exercise-fluidsystem/README.md 1970-01-01 01:00:00.000000000 +0100
-@@ -1,265 +0,0 @@
+@@ -1,266 +0,0 @@
-# Exercise Fluidsystem (DuMuX Course)
-
-The aim of this exercise is to get familiar with the _DuMu<sup>x</sup>_ way of implementing new components and using them in fluid systems (immiscible/mixture). In the scope of this exercise, different two-phase (liquid-liquid) fluid systems are implemented. A new fictitious component is implemented and used to describe a immiscible liquid phase. The second phase will only consist of water (Section 2). Furthermore, the mixture of the two components is implemented as two miscible phases (Section 3).
-
-## Problem set-up
-
--The domain has a size of 60 m x 60 m and contains two low-permeable lenses. Initially, the domain is fully water saturated and the fictitious component is injected through the middle portion of the upper boundary by means of a Neumann boundary condition. The remaining parts of the upper and the entire lower boundary are Neumann no-flow while on the two lateral sides Dirichlet boundary conditions are applied (hydrostatic conditions for the pressure and zero saturation).
+-The domain has a size of 60 x 60 m and contains two low-permeable lenses. Initially, the domain is fully water saturated and the fictitious component is injected through the middle portion of the upper boundary by means of a Neumann boundary condition. The remaining parts of the upper and the entire lower boundary are Neumann no-flow while on the two lateral sides Dirichlet boundary conditions are applied (hydrostatic conditions for the pressure and zero saturation).
-
--<img src="https://git.iws.uni-stuttgart.de/dumux-repositories/dumux-course/raw/master/exercises/extradoc/exercise-fluidsystem_setup.png" width="700">
+-
-
-
--## Task 1: Getting familiar with the code
+-## Preparing the exercise
-
--* Navigate to the directory `dumux-course/exercises/exercise-fluidsystem`.
+-* Navigate to the directory `dumux-course/exercises/exercise-fluidsystem`
+-
+-### 1. Getting familiar with the code
-
-Locate all the files you will need for this exercise
-* The shared __main file__ : `main.cc`
@@ -394,13 +362,13 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
-* The __spatial parameters file__: `spatialparams.hh`
-
-Furthermore you will find the following folders:
--* `binarycoefficients`: Stores headers containing data/methods on binary mixtures.
--* `components`: Stores headers containing data/methods on pure components.
+-* `binarycoefficients`: Stores headers containing data/methods on binary mixtures
+-* `components`: Stores headers containing data/methods on pure components
-* `fluidsystems`: Stores headers containing data/methods on mixtures of pure components. Uses methods from `binarycoefficients`.
-
--To see more components, fluid systems and binarycoefficients implementations, have a look at the folder `dumux/material`.
+-To see more components, fluidsystems and binarycoefficients implementations, have a look at the folder `dumux/material`.
-
--## Task 2: Implement a new component - Part a
+-### 2. Implement a new component
-
-In the following, the basic steps required to set the desired fluid system are outlined. Here, this is done in the __properties file__, i.e. for this part of the exercise the code shown below is taken from the `2pproperties.hh` file.
-
@@ -431,14 +399,14 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
-
-For a cell-centered scheme, you could derive from `CCTpfaModel` or `CCMpfaModel` instead (and, of course, include the right headers).
-
--One of the two phases should only contain the component water. We want to precompute tables on which the properties are then interpolated in order to save computational time. We need the properties of the component water in the properties file and in the problem file. We include the following two headers in our problem file, i.e. `2pproblem.hh` file. By including the problem file in the properties file, it can be accessed through it.
+-One of the two phases should only contain the component water. We want to precompute tables on which the properties are then interpolated in order to save computational time. Since we need this part in both problem file and properties file, we have to include the following two headers in our problem file, i.e. `2pproblem.hh` file, and the properties file has access to them through problem file.
-
-```c++
-// The water component
-#include <dumux/material/components/tabulatedcomponent.hh>
-#include <dumux/material/components/h2o.hh>
-```
--The other phase that will be created only contains our new component, where we want to implement an incompressible (Section 2.1) and a compressible (Section 2.2) variant.
+-The other phase that will be created only contains our new component, where we want to implement an incompressible and a compressible variant.
-The respective headers are prepared, but still incomplete. The compressible variant is still commented so that compilation does not fail when finishing the incompressible variant.
-
-```c++
@@ -460,7 +428,9 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
-#include <dumux/material/fluidsystems/1pliquid.hh>
-```
-
--which creates a _liquid phase_ from a given component. Finally, using all of the included classes, we set the fluid system property. The fluid system consists of two liquids, water (derived from the tabulated water component) and the new incompressible fictitious component. Both components need to be defined as liquid (`OnePLiquid`). The two liquid components are then form an immiscible fluid system (`TwoPImmiscible`).
+-which creates a _liquid phase_ from a given component. Finally, using all of the included classes we set the fluid system property by choosing that the water phase is liquid (`OnePLiquid`) and consists of the tabulated water component, and
+-the other phase is liquid as well and consists of the incompressible fictitious component. Both will make up the immiscible fluid system (`TwoPImmiscible`), consisting of two liquids each consisting of one component:
+-
-
-```c++
-// we use the immiscible fluid system here
@@ -472,8 +442,8 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
- using TabulatedH2O = Components::TabulatedComponent<Components::H2O<Scalar>>;
- using LiquidWaterPhase = typename FluidSystems::OnePLiquid<Scalar, TabulatedH2O>;
-
-- // TODO: dumux-course-task 2.2:
-- // Select the corresponding component for the task
+- // TODO: dumux-course-task 2.2
+- //****** Select the corresponding component for the task *****//
- // Uncomment first line and comment second line for using the compressible component
- using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyIncompressibleComponent<Scalar> >;
- // using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyCompressibleComponent<Scalar> >;
@@ -483,16 +453,16 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
-};
-```
-
--### Task 2.1: Incompressible component
+-### 2.1. Incompressible component
-
--Open the file `myincompressiblecomponent.hh`. A component should always derive from the _Base_ class - thus we include `dumux/material/components/base.hh`-, which defines the interface of a _DuMuX_ component with the possibly required functions to be overloaded by the actual implementation. Additionally, it is required that liquids are derived from the _Liquid_ class (see `dumux/material/components/liquid.hh`), gases from the _Gas_ class (see `dumux/material/components/gas.hh`) and solids from the _Solid_ class (see `dumux/material/components/solid.hh`), with functions specific to liquid, gas or solid.
+-Open the file `myincompressiblecomponent.hh`. A component should always derive from the _Base_ class - thus we include `dumux/material/components/base.hh`-, which defines the interface of a _DuMuX_ component with possibly required functions to be overloaded by the actual implementation. Additionally it is required for liquids to derive from the _Liquid_ class (see `dumux/material/components/liquid.hh`), for gases to derive from the _Gas_ class (see `dumux/material/components/gas.hh`) and for solids to derive from the _Solid_ class (see `dumux/material/components/solid.hh`), with functions specific to liquid, gas or solid.
-
-```c++
-/*!
- * \ingroup Components
-- * \brief A ficitious component to be implemented in exercise-fluidsystem 2.1.
+- * \brief A ficitious component to be implemented in exercise 3.
- *
-- * \tparam Scalar The type used for scalar values.
+- * \tparam Scalar The type used for scalar values
- */
-template <class Scalar>
-class MyIncompressibleComponent
@@ -507,7 +477,7 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
-| Parameter | unit | value |
-| -----| --------| -------- |
-| $`M`$ | $`kg/mol`$ | $`131.39 \cdot 10^{-3}`$ |
--| $`\rho_{liquid}`$ | $`kg/m^3`$ | $`1460.0`$ |
+-| $`\rho_{liquid}`$ | $`kg/m^3`$ | $`1460`$ |
-| $`\mu_{liquid}`$ | $`Pa \cdot s`$ | $`5.7 \cdot 10^{-4}`$ |
-
-In order to do so, have a look at the files `dumux/material/components/base.hh` and `dumux/material/components/liquid.hh` to see how the interfaces are defined and overload them accordingly.
@@ -520,29 +490,30 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
-./exercise_fluidsystem_a aparams.input
-```
-
--The saturation distribution of the nonwetting phase, here S$`_{napl}`$ (the phase consisting of our fictitious incompressible component), at the final simulation time should look like this:
+-The saturation distribution of the nonwetting phase S$`_n`$ (the phase consisting of our fictitious incompressible component) at the final simulation time should look like this:
-
-
-
--### Task 2.2: Compressible component
+-### 2.2. Compressible component
-
-We now want to implement a pressure-dependent density for our component. Open the file `mycompressiblecomponent.hh` and copy in the functions you implemented for the incompressible variant. Now substitute the method that returns the density by the following expression:
-
-$`\displaystyle \rho_{MyComp} = \rho_{min} + \frac{ \rho_{max} - \rho_{min} }{ 1 + \rho_{min}*e^{-1.0*k*(\rho_{max} - \rho_{min})*p} } `$
-
--where $`p`$ is the pressure and $`\rho_{min} = 1440.0 \, kg/m^3`$, $`\rho_{max} = 1480.0 \, kg/m^3`$ and $`k = 5 \cdot 10^{-7} `$. Also, make sure the header is included in the `2pproperties.hh` file by uncommenting the corresponding line. Furthermore, the new component has to be set as a liquid phase in the fluid system. To do so, search for `TODO: dumux-course-task 2.2`. Comment out the corresponding line and uncomment the other. The density distribution of this phase (rho$`_{napl}`$) at the final simulation time should look like this:
+-where $`p`$ is the pressure and $`\rho_{min} = 1440 `$, $`\rho_{max} = 1480 `$ and $`k = 5 \cdot 10^{-7} `$. Also, make sure the header is included in the `2pproperties.hh` file by uncommenting the corresponding line. Furthermore, the new component has to be set as a liquid phase in the fluid system. To do so, search for `TODO: dumux-course-task 2.2`. Comment out the corresponding line and uncomment the other. The density distribution of this phase (rhoN) at the final simulation time should look like this:
-
-
-
-You can plot the density of the phase consisting of your compressible component by setting `PlotDensity` in `aparams.input` to `true` and starting the simulation again.
-Compare the gnuplot output to the following plot of the density function from above:
-
--<img src="https://git.iws.uni-stuttgart.de/dumux-repositories/dumux-course/raw/master/exercises/extradoc/exercise-fluidsystem_a_densityfunction.png" width="500">
+-
-
--## Task 3: Implement a new fluid system - Part b
+-### 3. Implement a new fluid system
-
-The problem file and properties file for this part of the exercise are `2p2cproblem.hh` and `2p2cproperties.hh`, respectively.
--We now want to implement a new fluid system, which still consists of two liquid phases. However, one phase consists mainly of water and the other consists mainly of the previously implemented compressible component. We will now consider compositional effects, which is why we have to derive our _TypeTag_ (`ExerciseFluidsystemTwoPTwoC`) from a _TypeTag_ (`TwoPTwoC`) that holds the miscible two-phase two-component model properties:
+-We now want to implement a new fluid system, which still consists of two liquid phases. However, one phase consists mainly of water and the other consists mainly of the previously implemented compressible component. We will now consider compositional effects, which is why we have to derive our _TypeTag_ (`ExerciseFluidsystemTwoPTwoC`) from a _TypeTag_ (`TwoPTwoC`) that holds the miscible two-phase
+- two-component model properties:
-
-```c++
-// The numerical model
@@ -550,7 +521,7 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
-```
-
-```c++
--// Create new type tags
+-// Create a new type tag for the problem
-namespace TTag {
-struct ExerciseFluidsystemTwoPTwoC { using InheritsFrom = std::tuple<TwoPTwoC, BoxModel>; };
-} // end namespace TTag
@@ -578,10 +549,10 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
-In the `fluidsystems/h2omycompressiblecomponent.hh` file, your implemented compressible component and the binary coefficient files are already included.
-
-```c++
--// The ficitious component that was created in exercise-fluidsystem 2.2
+-// the ficitious component that was created in exercise-fluidsystem a
-#include <exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh>
-
--// The binary coefficients corresponding to this fluid system
+-// the binary coefficients corresponding to this fluid system
-#include <exercises/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh>
-```
-
@@ -599,9 +570,9 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
-./exercise_fluidsystem_b bparams.input
-```
-
--You will observe an error message and an abortion of the program. This is due to the fact that in order for the constraint solver and other mechanisms in the two-phase two-component model to work, an additional functionality in the component has to be implemented: the model has to know the vapour pressure. As in the previous exercise, check the `dumux/material/components/base.hh` file for this function and implement it into `mycompressiblecomponent.hh`. For the vapour pressure, use a value of $`3900.0 \, Pa`$.
+-You will observe an error message and an abortion of the program. This is due to the fact that in order for the constraint solver and other mechanisms in the two-phase two-component model to work, an additional functionality in the component has to be implemented: the model has to know the vapour pressure. As in the previous exercise, check the `dumux/material/components/base.hh` file for this function and implement it into `mycompressiblecomponent.hh`. For the vapour pressure, use a value of $`3900`$ Pa.
-
--## Task 4: Change wettability of the porous medium - Part b
+-### 4. Change wettability of the porous medium
-
-In the `spatialparams.hh` file, we can find the following function, with which we can specify which phase of the fluid system is to be considered as the wetting phase at a given position within the domain:
-
@@ -615,13 +586,9 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
-template<class FluidSystem>
-int wettingPhaseAtPos(const GlobalPosition& globalPos) const
-{
-- // Our fluid system is H2OMyCompressibleComponent.
+- // Our fluid system is H2OMyCompressibleComponent
- // We want to define water as the wetting phase in
-- // the entire domain (see fluid system for the phase indices).
--
-- // TODO: dumux-course-task 4:
-- // Adapt the following line so that the phase of our new component is
-- // the wetting phase, only within the lenses.
+- // the entire domain (see fluid system for the phase indices)
- return FluidSystem::phase0Idx;
-}
-```
@@ -634,37 +601,21 @@ diff -ruN exercises/exercise-fluidsystem/README.md exercises/solution/exercise-f
-./exercise_fluidsystem_b exercise_fluidsystem_b.input
-```
diff -ruN exercises/exercise-fluidsystem/spatialparams.hh exercises/solution/exercise-fluidsystem/spatialparams.hh
---- exercises/exercise-fluidsystem/spatialparams.hh 2024-07-17 14:37:31.333957691 +0200
-+++ exercises/solution/exercise-fluidsystem/spatialparams.hh 2024-07-17 14:37:31.333957691 +0200
-@@ -113,13 +113,17 @@
- template<class FluidSystem>
- int wettingPhaseAtPos(const GlobalPosition& globalPos) const
- {
-- // Our fluid system is H2OMyCompressibleComponent.
-+ // Our fluid system is H2OMyCompressibleComponent.
+--- exercises/exercise-fluidsystem/spatialparams.hh 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fluidsystem/spatialparams.hh 2024-07-17 11:27:28.882349286 +0200
+@@ -117,9 +117,12 @@
// We want to define water as the wetting phase in
-- // the entire domain (see fluid system for the phase indices).
-+ // the entire domain (see fluidsystem for the phase indices).
-+
-+ // TODO: dumux-course-task 4:
-+ // For task 4 (change of wettability), uncomment
-+ // the two commented lines below:
-+
-+ //if (isInLens(globalPos))
-+ // return FluidSystem::phase1Idx;
+ // the entire domain (see fluid system for the phase indices)
-- // TODO: dumux-course-task 4:
+- // TODO: dumux-course-task 4
- // Adapt the following line so that the phase of our new component is
- // the wetting phase, only within the lenses.
++ // For the last task (change of wettability), uncomment
++ // the two commented lines below
++
++ //if (isInLens(globalPos))
++ // return FluidSystem::phase1Idx;
++
return FluidSystem::phase0Idx;
}
-@@ -147,7 +151,7 @@
- Dune::FieldMatrix<Scalar, dim, dim> K_;
- Dune::FieldMatrix<Scalar, dim, dim> KLens_;
-
-- // Object that holds the values/parameters of the selected fluid matrix interaction relationship.
-+ // Object that holds the values/parameters of the selected fluid-matrix interaction relationship.
- const PcKrSwCurve pcKrSwCurve_;
- const PcKrSwCurve lensPcKrSwCurve_;
- };
diff --git a/.patches/exercise-fractures/exercise-fractures.patch b/.patches/exercise-fractures/exercise-fractures.patch
index c2e2d32c591e11814853a8e322f01388381a229f..2fab87df54d02e9c372919c601547d9771046742 100644
--- a/.patches/exercise-fractures/exercise-fractures.patch
+++ b/.patches/exercise-fractures/exercise-fractures.patch
@@ -1,6 +1,6 @@
diff -ruN exercises/exercise-fractures/CMakeLists.txt exercises/solution/exercise-fractures/CMakeLists.txt
---- exercises/exercise-fractures/CMakeLists.txt 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-fractures/CMakeLists.txt 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-fractures/CMakeLists.txt 2024-07-17 11:27:28.798347462 +0200
++++ exercises/solution/exercise-fractures/CMakeLists.txt 2024-07-17 11:27:28.798347462 +0200
@@ -1,9 +1,26 @@
# SPDX-FileCopyrightInfo: Copyright © DuMux-Course contributors, see AUTHORS.md in root folder
# SPDX-License-Identifier: GPL-3.0-or-later
@@ -32,7 +32,7 @@ diff -ruN exercises/exercise-fractures/CMakeLists.txt exercises/solution/exercis
add_input_file_links()
diff -ruN exercises/exercise-fractures/exercise_fractures_solution_a.input exercises/solution/exercise-fractures/exercise_fractures_solution_a.input
--- exercises/exercise-fractures/exercise_fractures_solution_a.input 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-fractures/exercise_fractures_solution_a.input 2024-06-26 11:12:30.529881015 +0200
++++ exercises/solution/exercise-fractures/exercise_fractures_solution_a.input 2020-11-12 12:07:04.129814352 +0100
@@ -0,0 +1,41 @@
+[TimeLoop]
+TEnd = 30000 # [s]
@@ -77,7 +77,7 @@ diff -ruN exercises/exercise-fractures/exercise_fractures_solution_a.input exerc
+SpatialParams.Barrier.Swr = 0.0
diff -ruN exercises/exercise-fractures/exercise_fractures_solution_b.input exercises/solution/exercise-fractures/exercise_fractures_solution_b.input
--- exercises/exercise-fractures/exercise_fractures_solution_b.input 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-fractures/exercise_fractures_solution_b.input 2024-06-26 11:12:30.529881015 +0200
++++ exercises/solution/exercise-fractures/exercise_fractures_solution_b.input 2020-11-12 12:07:04.129814352 +0100
@@ -0,0 +1,41 @@
+[TimeLoop]
+TEnd = 30000 # [s]
@@ -122,7 +122,7 @@ diff -ruN exercises/exercise-fractures/exercise_fractures_solution_b.input exerc
+SpatialParams.Barrier.Swr = 0.0
diff -ruN exercises/exercise-fractures/exercise_fractures_solution_c.input exercises/solution/exercise-fractures/exercise_fractures_solution_c.input
--- exercises/exercise-fractures/exercise_fractures_solution_c.input 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-fractures/exercise_fractures_solution_c.input 2024-06-26 11:12:30.529881015 +0200
++++ exercises/solution/exercise-fractures/exercise_fractures_solution_c.input 2020-11-12 12:07:04.129814352 +0100
@@ -0,0 +1,41 @@
+[TimeLoop]
+TEnd = 30000 # [s]
@@ -166,8 +166,8 @@ diff -ruN exercises/exercise-fractures/exercise_fractures_solution_c.input exerc
+SpatialParams.Barrier.Snr = 0.0
+SpatialParams.Barrier.Swr = 0.0
diff -ruN exercises/exercise-fractures/fractureproblem.hh exercises/solution/exercise-fractures/fractureproblem.hh
---- exercises/exercise-fractures/fractureproblem.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-fractures/fractureproblem.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-fractures/fractureproblem.hh 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fractures/fractureproblem.hh 2024-07-17 11:27:28.882349286 +0200
@@ -59,6 +59,10 @@
std::shared_ptr<typename ParentType::SpatialParams> spatialParams,
const std::string& paramGroup)
@@ -199,8 +199,8 @@ diff -ruN exercises/exercise-fractures/fractureproblem.hh exercises/solution/exe
return values;
}
diff -ruN exercises/exercise-fractures/fracturespatialparams.hh exercises/solution/exercise-fractures/fracturespatialparams.hh
---- exercises/exercise-fractures/fracturespatialparams.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-fractures/fracturespatialparams.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-fractures/fracturespatialparams.hh 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fractures/fracturespatialparams.hh 2024-07-17 11:27:28.886349373 +0200
@@ -41,6 +41,7 @@
using Element = typename GridView::template Codim<0>::Entity;
using GlobalPosition = typename Element::Geometry::GlobalCoordinate;
@@ -322,8 +322,8 @@ diff -ruN exercises/exercise-fractures/fracturespatialparams.hh exercises/soluti
} // end namespace Dumux
diff -ruN exercises/exercise-fractures/main.cc exercises/solution/exercise-fractures/main.cc
---- exercises/exercise-fractures/main.cc 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-fractures/main.cc 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-fractures/main.cc 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fractures/main.cc 2024-07-17 11:27:28.886349373 +0200
@@ -245,4 +245,5 @@
Parameters::print();
@@ -331,8 +331,8 @@ diff -ruN exercises/exercise-fractures/main.cc exercises/solution/exercise-fract
+
}// end main
diff -ruN exercises/exercise-fractures/matrixproblem.hh exercises/solution/exercise-fractures/matrixproblem.hh
---- exercises/exercise-fractures/matrixproblem.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-fractures/matrixproblem.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-fractures/matrixproblem.hh 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fractures/matrixproblem.hh 2024-07-17 11:27:28.886349373 +0200
@@ -65,11 +65,18 @@
: ParentType(gridGeometry, spatialParams, paramGroup)
, boundaryOverPressure_(getParamFromGroup<Scalar>(paramGroup, "Problem.BoundaryOverPressure"))
@@ -456,8 +456,8 @@ diff -ruN exercises/exercise-fractures/matrixproblem.hh exercises/solution/exerc
} // end namespace Dumux
diff -ruN exercises/exercise-fractures/matrixspatialparams.hh exercises/solution/exercise-fractures/matrixspatialparams.hh
---- exercises/exercise-fractures/matrixspatialparams.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-fractures/matrixspatialparams.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-fractures/matrixspatialparams.hh 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fractures/matrixspatialparams.hh 2024-07-17 11:27:28.886349373 +0200
@@ -40,6 +40,7 @@
using Element = typename GridView::template Codim<0>::Entity;
using GlobalPosition = typename Element::Geometry::GlobalCoordinate;
@@ -467,8 +467,8 @@ diff -ruN exercises/exercise-fractures/matrixspatialparams.hh exercises/solution
public:
diff -ruN exercises/exercise-fractures/params.input exercises/solution/exercise-fractures/params.input
---- exercises/exercise-fractures/params.input 2024-06-26 11:12:30.501880032 +0200
-+++ exercises/solution/exercise-fractures/params.input 2024-06-26 11:12:30.529881015 +0200
+--- exercises/exercise-fractures/params.input 2023-10-13 18:08:54.714521020 +0200
++++ exercises/solution/exercise-fractures/params.input 2020-05-08 12:16:18.919855703 +0200
@@ -5,6 +5,9 @@
[Problem]
@@ -510,8 +510,8 @@ diff -ruN exercises/exercise-fractures/params.input exercises/solution/exercise-
+SpatialParams.Swr = 0.0
+SpatialParams.SwrBarrier = 0.0
diff -ruN exercises/exercise-fractures/properties.hh exercises/solution/exercise-fractures/properties.hh
---- exercises/exercise-fractures/properties.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-fractures/properties.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-fractures/properties.hh 2024-07-17 11:27:28.874349112 +0200
++++ exercises/solution/exercise-fractures/properties.hh 2024-07-17 11:27:28.886349373 +0200
@@ -45,6 +45,7 @@
// with tpfa.
#include <dumux/multidomain/facet/cellcentered/tpfa/properties.hh>
@@ -530,7 +530,7 @@ diff -ruN exercises/exercise-fractures/properties.hh exercises/solution/exercise
struct Grid<TypeTag, TTag::MatrixProblem> { using type = Dune::ALUGrid<2, 2, Dune::simplex, Dune::conforming>; };
template<class TypeTag>
diff -ruN exercises/exercise-fractures/README.md exercises/solution/exercise-fractures/README.md
---- exercises/exercise-fractures/README.md 2024-06-26 11:12:30.501880032 +0200
+--- exercises/exercise-fractures/README.md 2023-10-13 18:08:54.714521020 +0200
+++ exercises/solution/exercise-fractures/README.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,216 +0,0 @@
-# Exercise Discrete Fractures (DuMuX Course)
diff --git a/.patches/exercise-grids/exercise-grids.patch b/.patches/exercise-grids/exercise-grids.patch
index 255f321ca57b77705c2d4dc46736738d033febb5..18a57aa57dfb0752bf50daeecc76416694088540 100644
--- a/.patches/exercise-grids/exercise-grids.patch
+++ b/.patches/exercise-grids/exercise-grids.patch
@@ -1,6 +1,6 @@
diff -ruN exercises/exercise-grids/CMakeLists.txt exercises/solution/exercise-grids/CMakeLists.txt
---- exercises/exercise-grids/CMakeLists.txt 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-grids/CMakeLists.txt 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-grids/CMakeLists.txt 2024-07-17 13:58:18.212850366 +0200
++++ exercises/solution/exercise-grids/CMakeLists.txt 2024-07-17 13:58:18.220850484 +0200
@@ -2,7 +2,7 @@
# SPDX-License-Identifier: GPL-3.0-or-later
@@ -11,7 +11,7 @@ diff -ruN exercises/exercise-grids/CMakeLists.txt exercises/solution/exercise-gr
# add a symlink for each input file
diff -ruN exercises/exercise-grids/grids/grid_structured.dgf exercises/solution/exercise-grids/grids/grid_structured.dgf
---- exercises/exercise-grids/grids/grid_structured.dgf 2024-06-26 11:12:30.501880032 +0200
+--- exercises/exercise-grids/grids/grid_structured.dgf 2023-10-13 18:08:54.718521039 +0200
+++ exercises/solution/exercise-grids/grids/grid_structured.dgf 1970-01-01 01:00:00.000000000 +0100
@@ -1,6 +0,0 @@
-DGF
@@ -21,7 +21,7 @@ diff -ruN exercises/exercise-grids/grids/grid_structured.dgf exercises/solution/
-30 20 % cells in x and y direction
-#
diff -ruN exercises/exercise-grids/grids/grid_structured.geo exercises/solution/exercise-grids/grids/grid_structured.geo
---- exercises/exercise-grids/grids/grid_structured.geo 2024-06-26 11:12:30.501880032 +0200
+--- exercises/exercise-grids/grids/grid_structured.geo 2023-10-13 18:08:54.718521039 +0200
+++ exercises/solution/exercise-grids/grids/grid_structured.geo 1970-01-01 01:00:00.000000000 +0100
@@ -1,30 +0,0 @@
-X = 60;
@@ -55,8 +55,8 @@ diff -ruN exercises/exercise-grids/grids/grid_structured.geo exercises/solution/
-Transfinite Volume "*";
-
diff -ruN exercises/exercise-grids/main.cc exercises/solution/exercise-grids/main.cc
---- exercises/exercise-grids/main.cc 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-grids/main.cc 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-grids/main.cc 2024-07-17 13:58:18.212850366 +0200
++++ exercises/solution/exercise-grids/main.cc 2024-07-17 13:58:18.220850484 +0200
@@ -6,7 +6,7 @@
//
/*!
@@ -67,8 +67,8 @@ diff -ruN exercises/exercise-grids/main.cc exercises/solution/exercise-grids/mai
#include <config.h>
diff -ruN exercises/exercise-grids/params.input exercises/solution/exercise-grids/params.input
---- exercises/exercise-grids/params.input 2024-07-17 14:37:31.113952826 +0200
-+++ exercises/solution/exercise-grids/params.input 2024-06-26 11:12:30.529881015 +0200
+--- exercises/exercise-grids/params.input 2024-07-17 14:09:52.163005898 +0200
++++ exercises/solution/exercise-grids/params.input 2024-07-17 14:10:01.563142989 +0200
@@ -3,13 +3,34 @@
TEnd = 3.154e9 # in seconds, i.e ten years
@@ -110,8 +110,8 @@ diff -ruN exercises/exercise-grids/params.input exercises/solution/exercise-grid
[Problem]
Name = grid_exercise
diff -ruN exercises/exercise-grids/properties.hh exercises/solution/exercise-grids/properties.hh
---- exercises/exercise-grids/properties.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-grids/properties.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-grids/properties.hh 2024-07-17 13:58:18.212850366 +0200
++++ exercises/solution/exercise-grids/properties.hh 2024-07-17 13:58:18.220850484 +0200
@@ -5,9 +5,9 @@
// SPDX-License-Identifier: GPL-3.0-or-later
//
@@ -163,7 +163,7 @@ diff -ruN exercises/exercise-grids/properties.hh exercises/solution/exercise-gri
// Set the problem property
template<class TypeTag>
diff -ruN exercises/exercise-grids/README.md exercises/solution/exercise-grids/README.md
---- exercises/exercise-grids/README.md 2024-07-17 14:02:52.092041122 +0200
+--- exercises/exercise-grids/README.md 2024-07-17 10:58:43.901583076 +0200
+++ exercises/solution/exercise-grids/README.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,159 +0,0 @@
-# Exercise Grids (DuMuX course)
diff --git a/.patches/exercise-mainfile/exercise-mainfile.patch b/.patches/exercise-mainfile/exercise-mainfile.patch
index 1a3a3a6b0fb4d61ee9c8c37f557971270f60a2bb..4ef11fa075453c37a192a78f7a20799e39cc158d 100644
--- a/.patches/exercise-mainfile/exercise-mainfile.patch
+++ b/.patches/exercise-mainfile/exercise-mainfile.patch
@@ -1,6 +1,6 @@
diff -ruN exercises/exercise-mainfile/CMakeLists.txt exercises/solution/exercise-mainfile/CMakeLists.txt
---- exercises/exercise-mainfile/CMakeLists.txt 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-mainfile/CMakeLists.txt 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-mainfile/CMakeLists.txt 2024-07-17 14:07:42.177109339 +0200
++++ exercises/solution/exercise-mainfile/CMakeLists.txt 2024-07-17 13:58:18.220850484 +0200
@@ -2,14 +2,9 @@
# SPDX-License-Identifier: GPL-3.0-or-later
@@ -22,7 +22,7 @@ diff -ruN exercises/exercise-mainfile/CMakeLists.txt exercises/solution/exercise
+add_input_file_links()
\ No newline at end of file
diff -ruN exercises/exercise-mainfile/exercise1pamain.cc exercises/solution/exercise-mainfile/exercise1pamain.cc
---- exercises/exercise-mainfile/exercise1pamain.cc 2024-07-17 14:02:52.092041122 +0200
+--- exercises/exercise-mainfile/exercise1pamain.cc 2024-07-17 13:58:18.212850366 +0200
+++ exercises/solution/exercise-mainfile/exercise1pamain.cc 1970-01-01 01:00:00.000000000 +0100
@@ -1,130 +0,0 @@
-// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
@@ -157,7 +157,7 @@ diff -ruN exercises/exercise-mainfile/exercise1pamain.cc exercises/solution/exer
-}// end main
diff -ruN exercises/exercise-mainfile/exercise1pa_solution_main.cc exercises/solution/exercise-mainfile/exercise1pa_solution_main.cc
--- exercises/exercise-mainfile/exercise1pa_solution_main.cc 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-mainfile/exercise1pa_solution_main.cc 2024-07-17 14:02:52.096041149 +0200
++++ exercises/solution/exercise-mainfile/exercise1pa_solution_main.cc 2024-07-17 13:58:18.220850484 +0200
@@ -0,0 +1,130 @@
+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
+// vi: set et ts=4 sw=4 sts=4:
@@ -290,7 +290,7 @@ diff -ruN exercises/exercise-mainfile/exercise1pa_solution_main.cc exercises/sol
+
+}// end main
diff -ruN exercises/exercise-mainfile/exercise1pbmain.cc exercises/solution/exercise-mainfile/exercise1pbmain.cc
---- exercises/exercise-mainfile/exercise1pbmain.cc 2024-07-17 14:02:52.092041122 +0200
+--- exercises/exercise-mainfile/exercise1pbmain.cc 2024-07-17 13:58:18.212850366 +0200
+++ exercises/solution/exercise-mainfile/exercise1pbmain.cc 1970-01-01 01:00:00.000000000 +0100
@@ -1,120 +0,0 @@
-// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
@@ -414,7 +414,7 @@ diff -ruN exercises/exercise-mainfile/exercise1pbmain.cc exercises/solution/exer
-
-}// end main
diff -ruN exercises/exercise-mainfile/exercise1pcmain.cc exercises/solution/exercise-mainfile/exercise1pcmain.cc
---- exercises/exercise-mainfile/exercise1pcmain.cc 2024-07-17 14:02:52.092041122 +0200
+--- exercises/exercise-mainfile/exercise1pcmain.cc 2024-07-17 13:58:18.212850366 +0200
+++ exercises/solution/exercise-mainfile/exercise1pcmain.cc 1970-01-01 01:00:00.000000000 +0100
@@ -1,146 +0,0 @@
-// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
@@ -564,7 +564,7 @@ diff -ruN exercises/exercise-mainfile/exercise1pcmain.cc exercises/solution/exer
-
-}// end main
diff -ruN exercises/exercise-mainfile/exercise_mainfile_a.input exercises/solution/exercise-mainfile/exercise_mainfile_a.input
---- exercises/exercise-mainfile/exercise_mainfile_a.input 2024-06-26 11:12:30.505880174 +0200
+--- exercises/exercise-mainfile/exercise_mainfile_a.input 2023-10-13 18:08:54.718521039 +0200
+++ exercises/solution/exercise-mainfile/exercise_mainfile_a.input 1970-01-01 01:00:00.000000000 +0100
@@ -1,18 +0,0 @@
-[Grid]
@@ -587,7 +587,7 @@ diff -ruN exercises/exercise-mainfile/exercise_mainfile_a.input exercises/soluti
-BaseEpsilon = 1e-10
diff -ruN exercises/exercise-mainfile/exercise_mainfile_a_solution.input exercises/solution/exercise-mainfile/exercise_mainfile_a_solution.input
--- exercises/exercise-mainfile/exercise_mainfile_a_solution.input 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-mainfile/exercise_mainfile_a_solution.input 2024-07-17 14:37:31.113952826 +0200
++++ exercises/solution/exercise-mainfile/exercise_mainfile_a_solution.input 2024-07-17 14:08:11.237533494 +0200
@@ -0,0 +1,18 @@
+[Grid]
+LowerLeft = 0 0
@@ -608,7 +608,7 @@ diff -ruN exercises/exercise-mainfile/exercise_mainfile_a_solution.input exercis
+PriVarMagnitude = 1e5
+BaseEpsilon = 1e-10
diff -ruN exercises/exercise-mainfile/exercise_mainfile_b.input exercises/solution/exercise-mainfile/exercise_mainfile_b.input
---- exercises/exercise-mainfile/exercise_mainfile_b.input 2024-06-26 11:12:30.505880174 +0200
+--- exercises/exercise-mainfile/exercise_mainfile_b.input 2023-10-13 18:08:54.718521039 +0200
+++ exercises/solution/exercise-mainfile/exercise_mainfile_b.input 1970-01-01 01:00:00.000000000 +0100
@@ -1,17 +0,0 @@
-[Grid]
@@ -629,7 +629,7 @@ diff -ruN exercises/exercise-mainfile/exercise_mainfile_b.input exercises/soluti
-[Assembly.NumericDifference]
-PriVarMagnitude = 1e5
diff -ruN exercises/exercise-mainfile/exercise_mainfile_c.input exercises/solution/exercise-mainfile/exercise_mainfile_c.input
---- exercises/exercise-mainfile/exercise_mainfile_c.input 2024-07-17 14:37:31.113952826 +0200
+--- exercises/exercise-mainfile/exercise_mainfile_c.input 2024-07-17 14:08:27.885776444 +0200
+++ exercises/solution/exercise-mainfile/exercise_mainfile_c.input 1970-01-01 01:00:00.000000000 +0100
@@ -1,21 +0,0 @@
-[TimeLoop]
@@ -654,8 +654,8 @@ diff -ruN exercises/exercise-mainfile/exercise_mainfile_c.input exercises/soluti
-[Assembly.NumericDifference]
-PriVarMagnitude = 1e5
diff -ruN exercises/exercise-mainfile/properties.hh exercises/solution/exercise-mainfile/properties.hh
---- exercises/exercise-mainfile/properties.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-mainfile/properties.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-mainfile/properties.hh 2024-07-17 13:58:18.212850366 +0200
++++ exercises/solution/exercise-mainfile/properties.hh 2024-07-17 13:58:18.220850484 +0200
@@ -25,8 +25,8 @@
#include <dumux/porousmediumflow/1p/model.hh>
@@ -679,7 +679,7 @@ diff -ruN exercises/exercise-mainfile/properties.hh exercises/solution/exercise-
// the fluid system for compressible tests
diff -ruN exercises/exercise-mainfile/README.md exercises/solution/exercise-mainfile/README.md
---- exercises/exercise-mainfile/README.md 2024-07-17 14:02:52.092041122 +0200
+--- exercises/exercise-mainfile/README.md 2024-07-17 11:27:28.682344944 +0200
+++ exercises/solution/exercise-mainfile/README.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,196 +0,0 @@
-# Exercise Mainfiles (DuMuX course)
diff --git a/.patches/exercise-model/exercise-model.patch b/.patches/exercise-model/exercise-model.patch
index aa74587fadfc0bbb1739d5dfb5191041cb65c3cc..a7c60cd7ffe4c3e38c93acc6dea9eca2c5e6380d 100644
--- a/.patches/exercise-model/exercise-model.patch
+++ b/.patches/exercise-model/exercise-model.patch
@@ -1,6 +1,6 @@
diff -ruN exercises/exercise-model/CMakeLists.txt exercises/solution/exercise-model/CMakeLists.txt
---- exercises/exercise-model/CMakeLists.txt 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-model/CMakeLists.txt 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-model/CMakeLists.txt 2024-07-17 11:27:28.798347462 +0200
++++ exercises/solution/exercise-model/CMakeLists.txt 2024-07-17 11:27:28.802347549 +0200
@@ -3,5 +3,5 @@
dune_symlink_to_source_files(FILES images params.input)
@@ -9,8 +9,8 @@ diff -ruN exercises/exercise-model/CMakeLists.txt exercises/solution/exercise-mo
+dumux_add_test(NAME exercise_nonlineardiffusion_sol
SOURCES main.cc)
diff -ruN exercises/exercise-model/main.cc exercises/solution/exercise-model/main.cc
---- exercises/exercise-model/main.cc 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-model/main.cc 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-model/main.cc 2024-07-17 11:27:28.878349199 +0200
++++ exercises/solution/exercise-model/main.cc 2024-07-17 11:27:28.886349373 +0200
@@ -29,6 +29,8 @@
#include <dumux/nonlinear/newtonsolver.hh>
#include <dumux/assembly/fvassembler.hh>
@@ -148,7 +148,7 @@ diff -ruN exercises/exercise-model/main.cc exercises/solution/exercise-model/mai
}// end main
diff -ruN exercises/exercise-model/model.hh exercises/solution/exercise-model/model.hh
--- exercises/exercise-model/model.hh 1970-01-01 01:00:00.000000000 +0100
-+++ exercises/solution/exercise-model/model.hh 2024-07-17 14:02:52.096041149 +0200
++++ exercises/solution/exercise-model/model.hh 2024-07-17 11:27:28.886349373 +0200
@@ -0,0 +1,210 @@
+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
+// vi: set et ts=4 sw=4 sts=4:
@@ -361,7 +361,7 @@ diff -ruN exercises/exercise-model/model.hh exercises/solution/exercise-model/mo
+
+#endif
diff -ruN exercises/exercise-model/README.md exercises/solution/exercise-model/README.md
---- exercises/exercise-model/README.md 2024-07-16 15:45:21.113991695 +0200
+--- exercises/exercise-model/README.md 2024-07-16 16:28:05.559708538 +0200
+++ exercises/solution/exercise-model/README.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,195 +0,0 @@
-# Exercise Model (DuMuX course)
diff --git a/.patches/exercise-properties/exercise-properties.patch b/.patches/exercise-properties/exercise-properties.patch
index 51b238414f219128fb4cefa02e10443e07e5bd8d..d181a5253928601f05773c8d0e02f2f964f9e3b5 100644
--- a/.patches/exercise-properties/exercise-properties.patch
+++ b/.patches/exercise-properties/exercise-properties.patch
@@ -1,6 +1,6 @@
diff -ruN exercises/exercise-properties/CMakeLists.txt exercises/solution/exercise-properties/CMakeLists.txt
---- exercises/exercise-properties/CMakeLists.txt 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-properties/CMakeLists.txt 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-properties/CMakeLists.txt 2024-07-17 11:27:28.798347462 +0200
++++ exercises/solution/exercise-properties/CMakeLists.txt 2024-07-17 11:27:28.802347549 +0200
@@ -1,8 +1,8 @@
# SPDX-FileCopyrightInfo: Copyright © DuMux-Course contributors, see AUTHORS.md in root folder
# SPDX-License-Identifier: GPL-3.0-or-later
@@ -13,8 +13,8 @@ diff -ruN exercises/exercise-properties/CMakeLists.txt exercises/solution/exerci
# add a symlink for each input file
diff -ruN exercises/exercise-properties/mylocalresidual.hh exercises/solution/exercise-properties/mylocalresidual.hh
---- exercises/exercise-properties/mylocalresidual.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-properties/mylocalresidual.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-properties/mylocalresidual.hh 2024-07-17 11:27:28.878349199 +0200
++++ exercises/solution/exercise-properties/mylocalresidual.hh 2024-07-17 11:27:28.886349373 +0200
@@ -59,8 +59,6 @@
* \note The volVars can be different to allow computing
* the implicit euler time derivative here
@@ -51,8 +51,8 @@ diff -ruN exercises/exercise-properties/mylocalresidual.hh exercises/solution/ex
auto eqIdx = conti0EqIdx + phaseIdx;
flux[eqIdx] = fluxVars.advectiveFlux(phaseIdx, upwindTerm);
diff -ruN exercises/exercise-properties/problem.hh exercises/solution/exercise-properties/problem.hh
---- exercises/exercise-properties/problem.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-properties/problem.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-properties/problem.hh 2024-07-17 11:27:28.878349199 +0200
++++ exercises/solution/exercise-properties/problem.hh 2024-07-17 11:27:28.886349373 +0200
@@ -11,6 +11,7 @@
#ifndef DUMUX_INCOMPRESSIBLE_TWOP_TEST_PROBLEM_HH
#define DUMUX_INCOMPRESSIBLE_TWOP_TEST_PROBLEM_HH
@@ -77,8 +77,8 @@ diff -ruN exercises/exercise-properties/problem.hh exercises/solution/exercise-p
return values;
}
diff -ruN exercises/exercise-properties/properties.hh exercises/solution/exercise-properties/properties.hh
---- exercises/exercise-properties/properties.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-properties/properties.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-properties/properties.hh 2024-07-17 11:27:28.878349199 +0200
++++ exercises/solution/exercise-properties/properties.hh 2024-07-17 11:27:28.886349373 +0200
@@ -27,9 +27,7 @@
#include "spatialparams.hh"
@@ -103,7 +103,7 @@ diff -ruN exercises/exercise-properties/properties.hh exercises/solution/exercis
// Set the fluid system
template<class TypeTag>
diff -ruN exercises/exercise-properties/README.md exercises/solution/exercise-properties/README.md
---- exercises/exercise-properties/README.md 2024-06-26 11:12:30.505880174 +0200
+--- exercises/exercise-properties/README.md 2023-10-13 18:08:54.718521039 +0200
+++ exercises/solution/exercise-properties/README.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,75 +0,0 @@
-# Exercise Properties (DuMuX course)
diff --git a/.patches/exercise-runtimeparams/exercise-runtimeparams.patch b/.patches/exercise-runtimeparams/exercise-runtimeparams.patch
index e9ba41593573b01d072628cf3e096a4deaf6e775..4de557eb95d314a762991ab043c2542256516e34 100644
--- a/.patches/exercise-runtimeparams/exercise-runtimeparams.patch
+++ b/.patches/exercise-runtimeparams/exercise-runtimeparams.patch
@@ -1,6 +1,6 @@
diff -ruN exercises/exercise-runtimeparams/CMakeLists.txt exercises/solution/exercise-runtimeparams/CMakeLists.txt
---- exercises/exercise-runtimeparams/CMakeLists.txt 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-runtimeparams/CMakeLists.txt 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-runtimeparams/CMakeLists.txt 2024-07-17 13:58:18.212850366 +0200
++++ exercises/solution/exercise-runtimeparams/CMakeLists.txt 2024-07-17 13:58:18.220850484 +0200
@@ -2,7 +2,7 @@
# SPDX-License-Identifier: GPL-3.0-or-later
@@ -11,8 +11,8 @@ diff -ruN exercises/exercise-runtimeparams/CMakeLists.txt exercises/solution/exe
# add a symlink for each input file
diff -ruN exercises/exercise-runtimeparams/params.input exercises/solution/exercise-runtimeparams/params.input
---- exercises/exercise-runtimeparams/params.input 2024-06-26 11:12:30.505880174 +0200
-+++ exercises/solution/exercise-runtimeparams/params.input 2024-07-17 14:37:31.117952914 +0200
+--- exercises/exercise-runtimeparams/params.input 2024-07-17 14:07:04.624561096 +0200
++++ exercises/solution/exercise-runtimeparams/params.input 2024-07-17 14:06:46.236292578 +0200
@@ -12,11 +12,12 @@
AquiferDepth = 2700.0 # m
InjectionDuration = 2.628e6 # in seconds, i.e. one month
@@ -28,8 +28,8 @@ diff -ruN exercises/exercise-runtimeparams/params.input exercises/solution/exerc
Aquitard.Swr = 0.2
Aquitard.Snr = 0.0
diff -ruN exercises/exercise-runtimeparams/problem.hh exercises/solution/exercise-runtimeparams/problem.hh
---- exercises/exercise-runtimeparams/problem.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-runtimeparams/problem.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-runtimeparams/problem.hh 2024-07-17 13:58:18.212850366 +0200
++++ exercises/solution/exercise-runtimeparams/problem.hh 2024-07-17 13:58:18.220850484 +0200
@@ -81,12 +81,16 @@
injectionDuration_ = getParamFromGroup<Scalar>("Problem","InjectionDuration");
// TODO: dumux-course-task 2
@@ -67,8 +67,8 @@ diff -ruN exercises/exercise-runtimeparams/problem.hh exercises/solution/exercis
};
diff -ruN exercises/exercise-runtimeparams/properties.hh exercises/solution/exercise-runtimeparams/properties.hh
---- exercises/exercise-runtimeparams/properties.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-runtimeparams/properties.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-runtimeparams/properties.hh 2024-07-17 13:58:18.212850366 +0200
++++ exercises/solution/exercise-runtimeparams/properties.hh 2024-07-17 13:58:18.220850484 +0200
@@ -58,6 +58,6 @@
FluidSystems::H2ON2DefaultPolicy</*fastButSimplifiedRelations=*/ true> >;
};
@@ -78,7 +78,7 @@ diff -ruN exercises/exercise-runtimeparams/properties.hh exercises/solution/exer
#endif
diff -ruN exercises/exercise-runtimeparams/README.md exercises/solution/exercise-runtimeparams/README.md
---- exercises/exercise-runtimeparams/README.md 2024-07-17 14:02:52.092041122 +0200
+--- exercises/exercise-runtimeparams/README.md 2024-07-17 10:58:43.901583076 +0200
+++ exercises/solution/exercise-runtimeparams/README.md 1970-01-01 01:00:00.000000000 +0100
@@ -1,171 +0,0 @@
-# Exercise Runtime Parameters (DuMuX course)
@@ -253,8 +253,8 @@ diff -ruN exercises/exercise-runtimeparams/README.md exercises/solution/exercise
-
-* > __Task 4__: Using one of the bool `hasParam` functions, place an output in the problem file to alert the user where the parameter value comes from.
diff -ruN exercises/exercise-runtimeparams/spatialparams.hh exercises/solution/exercise-runtimeparams/spatialparams.hh
---- exercises/exercise-runtimeparams/spatialparams.hh 2024-07-17 14:02:52.092041122 +0200
-+++ exercises/solution/exercise-runtimeparams/spatialparams.hh 2024-07-17 14:02:52.096041149 +0200
+--- exercises/exercise-runtimeparams/spatialparams.hh 2024-07-17 13:58:18.212850366 +0200
++++ exercises/solution/exercise-runtimeparams/spatialparams.hh 2024-07-17 13:58:18.220850484 +0200
@@ -12,8 +12,8 @@
* fully implicit model.
*/
diff --git a/exercises/exercise-basic/injection2p2cproblem.hh b/exercises/exercise-basic/injection2p2cproblem.hh
deleted file mode 100644
index 06fa35d86fdcaf0907f00c586891fa1c91350e90..0000000000000000000000000000000000000000
--- a/exercises/exercise-basic/injection2p2cproblem.hh
+++ /dev/null
@@ -1,229 +0,0 @@
-// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
-// vi: set et ts=4 sw=4 sts=4:
-/*****************************************************************************
- * See the file COPYING for full copying permissions. *
- * *
- * This program is free software: you can redistribute it and/or modify *
- * it under the terms of the GNU General Public License as published by *
- * the Free Software Foundation, either version 3 of the License, or *
- * (at your option) any later version. *
- * *
- * This program is distributed in the hope that it will be useful, *
- * but WITHOUT ANY WARRANTY; without even the implied warranty of *
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
- * GNU General Public License for more details. *
- * *
- * You should have received a copy of the GNU General Public License *
- * along with this program. If not, see <http://www.gnu.org/licenses/>. *
- *****************************************************************************/
-/*!
- * \file
- *
- * \brief The two-phase porousmediumflow problem for exercise-basic
- */
-
-#ifndef DUMUX_EX_BASIC_PROBLEM_2P2C_HH
-#define DUMUX_EX_BASIC_PROBLEM_2P2C_HH
-
-#include <dumux/common/properties.hh>
-#include <dumux/common/boundarytypes.hh>
-#include <dumux/common/numeqvector.hh>
-#include <dumux/porousmediumflow/problem.hh>
-#include <dumux/material/binarycoefficients/h2o_n2.hh>
-
-namespace Dumux {
-
-/*!
- * \ingroup TwoPTwoCModel
- * \ingroup ImplicitTestProblems
- * \brief Gas injection problem where a gas (here nitrogen) is injected into a fully
- * water saturated medium. During buoyancy-driven upward migration the gas
- * passes a high temperature area.
- *
- * The domain is sized 60 m times 40 m.
- *
- * For the mass conservation equation Neumann boundary conditions are used on
- * the top, on the bottom and on the right of the domain, while Dirichlet conditions
- * apply on the left boundary.
- *
- * Gas is injected at the right boundary from 7 m to 15 m at a rate of
- * 0.001 kg/(s m), the remaining Neumann boundaries are no-flow
- * boundaries.
- *
- * At the Dirichlet boundaries a hydrostatic pressure and a gas saturation of zero a
- *
- * This problem uses the \ref TwoPModel model.
- */
-template<class TypeTag>
-class Injection2p2cProblem : public PorousMediumFlowProblem<TypeTag>
-{
- using ParentType = PorousMediumFlowProblem<TypeTag>;
- using Scalar = GetPropType<TypeTag, Properties::Scalar>;
- using Indices = typename GetPropType<TypeTag, Properties::ModelTraits>::Indices;
- using PrimaryVariables = GetPropType<TypeTag, Properties::PrimaryVariables>;
- using BoundaryTypes = Dumux::BoundaryTypes<GetPropType<TypeTag, Properties::ModelTraits>::numEq()>;
- using GridGeometry = GetPropType<TypeTag, Properties::GridGeometry>;
- using FVElementGeometry = typename GridGeometry::LocalView;
- using GridView = typename GridGeometry::GridView;
- using FluidSystem = GetPropType<TypeTag, Properties::FluidSystem>;
- using NumEqVector = Dumux::NumEqVector<PrimaryVariables>;
-
- static constexpr int dimWorld = GridView::dimensionworld;
- using Element = typename GridView::template Codim<0>::Entity;
- using GlobalPosition = typename Element::Geometry::GlobalCoordinate;
-
-public:
- Injection2p2cProblem(std::shared_ptr<const GridGeometry> gridGeometry)
- : ParentType(gridGeometry)
- {
- // Initialize the tables of the fluid system
- FluidSystem::init(/*tempMin=*/273.15,
- /*tempMax=*/423.15,
- /*numTemp=*/50,
- /*pMin=*/0.0,
- /*pMax=*/30e6,
- /*numP=*/300);
-
- // Name of the problem and output file
- // getParam<TYPE>("GROUPNAME.PARAMNAME") reads and sets parameter PARAMNAME
- // of type TYPE given in the group GROUPNAME from the input file
- name_ = getParam<std::string>("Problem.Name");
- // Depth of the aquifer, unit: m
- aquiferDepth_ = getParam<Scalar>("Problem.AquiferDepth");
- // The duration of the injection, unit: seconds
- injectionDuration_ = getParam<Scalar>("Problem.InjectionDuration");
- }
-
- /*!
- * \brief Returns the problem name
- *
- * This is used as a prefix for files generated by the simulation.
- */
- std::string name() const
- { return name_+"-2p2c"; }
-
- /*!
- * \brief Specifies which kind of boundary condition should be
- * used for which equation on a given boundary segment.
- *
- * \param globalPos The position for which the bc type should be evaluated
- */
- BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const
- {
- BoundaryTypes bcTypes;
- if (globalPos[0] < eps_)
- bcTypes.setAllDirichlet();
- else
- bcTypes.setAllNeumann();
-
- return bcTypes;
- }
-
- /*!
- * \brief Evaluates the boundary conditions for a Dirichlet
- * boundary segment
- *
- * \param globalPos The global position
- */
- PrimaryVariables dirichletAtPos(const GlobalPosition &globalPos) const
- {
- return initialAtPos(globalPos);
- }
-
- /*!
- * \brief Evaluate the boundary conditions for a Neumann
- * boundary segment.
- *
- * \param globalPos The position of the integration point of the boundary segment.
- */
- NumEqVector neumannAtPos(const GlobalPosition &globalPos) const
- {
- // Initialize values to zero, i.e. no-flow Neumann boundary conditions
- NumEqVector values(0.0);
-
- // If we are inside the injection zone set inflow Neumann boundary conditions
- if (injectionActive() && onInjectionBoundary(globalPos))
- {
- // TODO: dumux-course-task
- // Instead of setting -1e-4 here directly use totalAreaSpecificInflow_ in the computation.
-
- // Inject nitrogen. negative values mean injection.
- // Convert units from kg/(s*m^2) to mole/(s*m^2)
- values[Indices::conti0EqIdx + FluidSystem::N2Idx] = -1e-4/FluidSystem::molarMass(FluidSystem::N2Idx);
- values[Indices::conti0EqIdx + FluidSystem::H2OIdx] = 0.0;
- }
-
- return values;
- }
-
- /*!
- * \brief Evaluate the source term for all phases within a given
- * sub-control-volume.
- *
- * \param globalPos The position for which the source term should be evaluated
- */
- NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
- {
- return NumEqVector(0.0);
- }
-
- /*!
- * \brief Evaluate the initial value for a control volume.
- *
- * \param globalPos The position for which the initial condition should be evaluated.
- */
- PrimaryVariables initialAtPos(const GlobalPosition &globalPos) const
- {
- PrimaryVariables values(0.0);
- values.setState(Indices::firstPhaseOnly);
- // Get the water density at atmospheric conditions
- const Scalar densityW = FluidSystem::H2O::liquidDensity(this->spatialParams().temperatureAtPos(globalPos), 1.0e5);
-
- // Assume an initially hydrostatic liquid pressure profile
- // Note: we subtract rho_w*g*h because g is defined negative
- const Scalar pw = 1.0e5 - densityW*this->spatialParams().gravity(globalPos)[dimWorld-1]*(aquiferDepth_ - globalPos[dimWorld-1]);
-
- // Initially we have some nitrogen dissolved
- // Saturation mole fraction would be:
- // moleFracLiquidN2 = (pw + pc + p_vap^sat)/henry;
- const Scalar moleFracLiquidN2 = pw*0.95/BinaryCoeff::H2O_N2::henry(this->spatialParams().temperatureAtPos(globalPos));
-
- // Note that because we start with a single phase system the primary variables
- // are pl and x^w_N2. This will switch as soon after we start injecting to a two
- // phase system so the primary variables will be pl and Sn (nonwetting saturation).
- values[Indices::pressureIdx] = pw;
- values[Indices::switchIdx] = moleFracLiquidN2;
-
- return values;
- }
-
- //! Set the time for the time dependent boundary conditions (called from main)
- void setTime(Scalar time)
- { time_ = time; }
-
- //! Return true if the injection is currently active
- bool injectionActive() const
- { return time_ < injectionDuration_; }
-
- //! Return true if the given position is in the injection boundary region
- bool onInjectionBoundary(const GlobalPosition& globalPos) const
- {
- return globalPos[1] < 15. + eps_
- && globalPos[1] > 7. - eps_
- && globalPos[0] > this->gridGeometry().bBoxMax()[0] - eps_;
- }
-
-private:
- static constexpr Scalar eps_ = 1e-6;
- std::string name_; //! Problem name
- Scalar aquiferDepth_; //! Depth of the aquifer in m
- Scalar injectionDuration_; //! Duration of the injection in seconds
- Scalar time_;
- //TODO: dumux-course-task
- //Define the Scalar totalAreaSpecificInflow_ here
-
-};
-
-} //end namespace Dumux
-
-#endif
diff --git a/exercises/exercise-basic/injection2pniproblem.hh b/exercises/exercise-basic/injection2pniproblem.hh
index 4206d730affa56b27933fb2d551aa66ac7b362d8..e2b277e76e846f396f5c6425884a585d8d2e9850 100644
--- a/exercises/exercise-basic/injection2pniproblem.hh
+++ b/exercises/exercise-basic/injection2pniproblem.hh
@@ -71,13 +71,13 @@ public:
/*pMax=*/30e6,
/*numP=*/300);
- // Name of the problem and output file
+ // name of the problem and output file
// getParam<TYPE>("GROUPNAME.PARAMNAME") reads and sets parameter PARAMNAME
// of type TYPE given in the group GROUPNAME from the input file
name_ = getParam<std::string>("Problem.Name");
- // Depth of the aquifer, unit: m
+ // depth of the aquifer, unit: m
aquiferDepth_ = getParam<Scalar>("Problem.AquiferDepth");
- // The duration of the injection, unit: seconds
+ // the duration of the injection, unit: seconds
injectionDuration_ = getParam<Scalar>("Problem.InjectionDuration");
}
@@ -138,21 +138,21 @@ public:
*/
NumEqVector neumannAtPos(const GlobalPosition &globalPos) const
{
- // Initialize values to zero, i.e. no-flow Neumann boundary conditions
+ // initialize values to zero, i.e. no-flow Neumann boundary conditions
NumEqVector values(0.0);
- // If we are inside the injection zone set inflow Neumann boundary conditions
+ // if we are inside the injection zone set inflow Neumann boundary conditions
if (injectionActive() && onInjectionBoundary(globalPos))
{
- // Inject nitrogen. negative values mean injection
- // Unit: kg/(s*m^2)
+ // inject nitrogen. Negative values mean injection
+ // unit: kg/(s*m^2)
values[Indices::conti0EqIdx + FluidSystem::N2Idx] = -1e-4;
values[Indices::conti0EqIdx + FluidSystem::H2OIdx] = 0.0;
/*!
* TODO:dumux-course-task 4:
* Set Neumann noflow conditions for the energy equation everywhere else except the left boundary.
- * Additionally, consider the energy flux at the injection point which is equal to the product of the respective mass flux and the matching enthalpy. Use the function `gasEnthalpy(temperature,pressure)` from the N2 component to access the necessary enthalpy.
+ * Additionally, consider the energy flux at the injection point which is equal to the product of the respective mass flux and the matching enthalpy. Use the function `gasEnthalpy(temperature,pressure)` from the N2 component to access the necessary enthalpy.
* hint: use `Indices::energyEqIdx` to access the entry belonging to the energy flux.
*/
}
@@ -183,8 +183,8 @@ public:
// get the water density at atmospheric conditions
const Scalar densityW = FluidSystem::H2O::liquidDensity(283.15, 1.0e5);
- // Assume an initially hydrostatic liquid pressure profile
- // Note: we subtract rho_w*g*h because g is defined negative
+ // assume an initially hydrostatic liquid pressure profile
+ // note: we subtract rho_w*g*h because g is defined negative
const Scalar pw = 1.0e5 - densityW*this->spatialParams().gravity(globalPos)[dimWorld-1]*(aquiferDepth_ - globalPos[dimWorld-1]);
values[Indices::pressureIdx] = pw;
diff --git a/exercises/exercise-basic/injection2pproblem.hh b/exercises/exercise-basic/injection2pproblem.hh
index a171f36f565c3a38319180a558a5155272ff62b0..02e7a99c16ccb8138322d1551d43f11f0e64fba7 100644
--- a/exercises/exercise-basic/injection2pproblem.hh
+++ b/exercises/exercise-basic/injection2pproblem.hh
@@ -100,10 +100,10 @@ public:
BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const
{
BoundaryTypes bcTypes;
- // Set the left of the domain (with the global position in "0 = x" direction as a Dirichlet boundary
+ // set the left of the domain (with the global position in "0 = x" direction as a Dirichlet boundary
if (globalPos[0] < eps_)
bcTypes.setAllDirichlet();
- // Set all other as Neumann boundaries
+ // set all other as Neumann boundaries
else
bcTypes.setAllNeumann();
@@ -129,10 +129,10 @@ public:
*/
NumEqVector neumannAtPos(const GlobalPosition &globalPos) const
{
- // Initialize values to zero, i.e. no-flow Neumann boundary conditions
+ // initialize values to zero, i.e. no-flow Neumann boundary conditions
NumEqVector values(0.0);
- // If we are inside the injection zone set inflow Neumann boundary conditions
+ // if we are inside the injection zone set inflow Neumann boundary conditions
// using < boundary + eps_ or > boundary - eps_ is safer for floating point comparisons
// than using <= or >= as it is robust with regard to imprecision introduced by rounding errors.
if (injectionActive() && onInjectionBoundary(globalPos))
@@ -167,11 +167,11 @@ public:
{
PrimaryVariables values(0.0);
- // Get the water density at atmospheric conditions
+ // get the water density at atmospheric conditions
const Scalar densityW = FluidSystem::H2O::liquidDensity(this->spatialParams().temperatureAtPos(globalPos), 1.0e5);
- // Assume an initially hydrostatic liquid pressure profile
- // Note: we subtract rho_w*g*h because g is defined negative
+ // assume an initially hydrostatic liquid pressure profile
+ // note: we subtract rho_w*g*h because g is defined negative
const Scalar pw = 1.0e5 - densityW*this->spatialParams().gravity(globalPos)[dimWorld-1]*(aquiferDepth_ - globalPos[dimWorld-1]);
values[Indices::pressureIdx] = pw;
diff --git a/exercises/exercise-fluidsystem/2p2cproblem.hh b/exercises/exercise-fluidsystem/2p2cproblem.hh
index 3c78896ff29b1929bf9c620df84822eca7f51ba2..2cf745904c93c067b052bc4b90bd82abe5955a25 100644
--- a/exercises/exercise-fluidsystem/2p2cproblem.hh
+++ b/exercises/exercise-fluidsystem/2p2cproblem.hh
@@ -7,12 +7,12 @@
/*!
* \file
*
- * \brief Tutorial problem for a fully coupled two-phase two-component box model.
+ * \brief Tutorial problem for a fully coupled two phase-two component box model.
*/
#ifndef DUMUX_EXERCISE_FLUIDSYSTEM_B_PROBLEM_HH
#define DUMUX_EXERCISE_FLUIDSYSTEM_B_PROBLEM_HH
-// The base porous media box problem.
+// The base porous media box problem
#include <dumux/porousmediumflow/problem.hh>
#include <dumux/common/boundarytypes.hh>
#include <dumux/common/properties.hh>
@@ -23,7 +23,7 @@ namespace Dumux {
/*!
* \ingroup TwoPBoxModel
*
- * \brief Tutorial problem for a fully coupled two-phase two-component box model.
+ * \brief Tutorial problem for a fully coupled two phase-two component box model.
*/
template <class TypeTag>
class ExerciseFluidsystemProblemTwoPTwoC : public PorousMediumFlowProblem<TypeTag>
@@ -64,7 +64,7 @@ public:
* \brief Specifies which kind of boundary condition should be
* used for which equation on a given boundary segment.
*
- * \param globalPos The position for which the bc type should be evaluated.
+ * \param globalPos The position for which the bc type should be evaluated
*/
BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const
{
@@ -80,9 +80,9 @@ public:
/*!
* \brief Evaluates the boundary conditions for a Dirichlet
- * boundary segment.
+ * boundary segment
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
PrimaryVariables dirichletAtPos(const GlobalPosition &globalPos) const
{
@@ -113,7 +113,7 @@ public:
}
else
{
- // no-flow on the remaining Neumann-boundaries
+ // no-flow on the remaining Neumann-boundaries.
values[Indices::conti0EqIdx + FluidSystem::H2OIdx] = 0;
values[Indices::conti0EqIdx + FluidSystem::NAPLIdx] = 0;
}
@@ -125,7 +125,7 @@ public:
/*!
* \brief Evaluate the initial value for a control volume.
*
- * \param globalPos The position for which the initial condition should be evaluated.
+ * \param globalPos The position for which the initial condition should be evaluated
*
* For this method, the \a values parameter stores primary
* variables.
@@ -149,9 +149,9 @@ public:
/*!
- * \brief Returns the source term.
+ * \brief Returns the source term
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
{
@@ -162,8 +162,8 @@ public:
private:
- Scalar eps_; //! Small epsilon value
- Scalar depthBOR_; //! Depth at the bottom of the reservoir
+ Scalar eps_; //! small epsilon value
+ Scalar depthBOR_; //! depth at the bottom of the reservoir
};
} // end namespace Dumux
diff --git a/exercises/exercise-fluidsystem/2p2cproperties.hh b/exercises/exercise-fluidsystem/2p2cproperties.hh
index fc13a8a68f833c35690ee8fddeba36d6b3cf43ed..5bdc04c72bbc618ddbdd1ef4df3696f56b7c2e2d 100644
--- a/exercises/exercise-fluidsystem/2p2cproperties.hh
+++ b/exercises/exercise-fluidsystem/2p2cproperties.hh
@@ -27,11 +27,12 @@
// The fluid system that is created in this exercise
#include "fluidsystems/h2omycompressiblecomponent.hh"
-// The problem file, where setup-specific boundary and initial conditions are defined
+// The problem file, where setup-specific boundary and initial conditions are defined.
#include"2p2cproblem.hh"
namespace Dumux::Properties {
+// Create a new type tag for the problem
// Create new type tags
namespace TTag {
struct ExerciseFluidsystemTwoPTwoC { using InheritsFrom = std::tuple<TwoPTwoC, BoxModel>; };
diff --git a/exercises/exercise-fluidsystem/2pproblem.hh b/exercises/exercise-fluidsystem/2pproblem.hh
index 4f9d17b2b1396d61ff90dbf7d8050452597f1ae5..8258644f199076f13dd7500364218ebde2ed5051 100644
--- a/exercises/exercise-fluidsystem/2pproblem.hh
+++ b/exercises/exercise-fluidsystem/2pproblem.hh
@@ -7,7 +7,7 @@
/*!
* \file
*
- * \brief Tutorial problem for a fully coupled two-phase box model.
+ * \brief Tutorial problem for a fully coupled twophase box model.
*/
#ifndef DUMUX_EXERCISE_FLUIDSYSTEM_A_PROBLEM_HH
#define DUMUX_EXERCISE_FLUIDSYSTEM_A_PROBLEM_HH
@@ -29,7 +29,7 @@ namespace Dumux {
/*!
* \ingroup TwoPBoxModel
- * \brief Tutorial problem for a fully coupled two-phase box model.
+ * \brief Tutorial problem for a fully coupled twophase box model.
*/
template <class TypeTag>
class ExerciseFluidsystemProblemTwoP : public PorousMediumFlowProblem<TypeTag>
@@ -86,7 +86,7 @@ public:
* \brief Specifies which kind of boundary condition should be
* used for which equation on a given boundary segment.
*
- * \param globalPos The position for which the bc type should be evaluated.
+ * \param globalPos The position for which the bc type should be evaluated
*/
BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const
{
@@ -102,9 +102,9 @@ public:
/*!
* \brief Evaluates the boundary conditions for a Dirichlet
- * boundary segment.
+ * boundary segment
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
PrimaryVariables dirichletAtPos(const GlobalPosition& globalPos) const
{
@@ -136,7 +136,7 @@ public:
}
else
{
- // no-flow on the remaining Neumann-boundaries
+ // no-flow on the remaining Neumann-boundaries.
values[contiWEqIdx] = 0;
values[contiNEqIdx] = 0;
}
@@ -148,7 +148,7 @@ public:
/*!
* \brief Evaluate the initial value for a control volume.
*
- * \param globalPos The position for which the initial condition should be evaluated.
+ * \param globalPos The position for which the initial condition should be evaluated
*
* For this method, the \a values parameter stores primary
* variables.
@@ -167,9 +167,9 @@ public:
/*!
- * \brief Returns the source term.
+ * \brief Returns the source term
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
NumEqVector sourceAtPos(const GlobalPosition& globalPos) const
{
@@ -207,9 +207,9 @@ private:
gnuplot_.plot("YourComponentPhase_density");
}
- Scalar eps_; //! Small epsilon value
- Scalar depthBOR_; //! Depth at the bottom of the reservoir
- Dumux::GnuplotInterface<double> gnuplot_; //! Collects data for plotting
+ Scalar eps_; //! small epsilon value
+ Scalar depthBOR_; //! depth at the bottom of the reservoir
+ Dumux::GnuplotInterface<double> gnuplot_; //! collects data for plotting
};
} // end namespace Dumux
diff --git a/exercises/exercise-fluidsystem/2pproperties.hh b/exercises/exercise-fluidsystem/2pproperties.hh
index 3807ca6453e9288a6ccfd6bfcded6985ebf5bf4a..af59530961b7ef8345375bda8e24605dcd8ecddb 100644
--- a/exercises/exercise-fluidsystem/2pproperties.hh
+++ b/exercises/exercise-fluidsystem/2pproperties.hh
@@ -7,7 +7,7 @@
/*!
* \file
*
- * \brief The properties file for two-phase exercise-fluidsystem.
+ * \brief The properties file for two phase exercise-fluidsystem.
*/
#ifndef DUMUX_EXERCISE_FLUIDSYSTEM_A_PROPERTIES_HH
#define DUMUX_EXERCISE_FLUIDSYSTEM_A_PROPERTIES_HH
@@ -25,25 +25,26 @@
#include "spatialparams.hh"
// The components that will be created in this exercise
-// TODO: dumux-course-task 2.2:
-// Include the corresponding header for the task
+// TODO: dumux-course-task 2.2
+//****** Include the corresponding header for the task *****//
#include "components/myincompressiblecomponent.hh"
// #include "components/mycompressiblecomponent.hh"
// We will only have liquid phases here
#include <dumux/material/fluidsystems/1pliquid.hh>
-// The two-phase immiscible fluidsystem
+// The two-phase immiscible fluid system
#include <dumux/material/fluidsystems/2pimmiscible.hh>
// The interface to create plots during simulation
#include <dumux/io/gnuplotinterface.hh>
-// The problem file, where setup-specific boundary and initial conditions are defined
+// The problem file, where setup-specific boundary and initial conditions are defined.
#include"2pproblem.hh"
namespace Dumux::Properties {
+// Create a new type tag for the problem
// Create new type tags
namespace TTag {
struct ExerciseFluidsystemTwoP { using InheritsFrom = std::tuple<TwoP, BoxModel>; };
@@ -68,7 +69,7 @@ public:
template<class TypeTag>
struct Grid<TypeTag, TTag::ExerciseFluidsystemTwoP> { using type = Dune::YaspGrid<2>; };
-// We use the immiscible fluid system here
+// we use the immiscible fluid system here
template<class TypeTag>
struct FluidSystem<TypeTag, TTag::ExerciseFluidsystemTwoP>
{
@@ -77,8 +78,8 @@ private:
using TabulatedH2O = Components::TabulatedComponent<Components::H2O<Scalar>>;
using LiquidWaterPhase = typename FluidSystems::OnePLiquid<Scalar, TabulatedH2O>;
- // TODO: dumux-course-task 2.2:
- // Select the corresponding component for the task
+ // TODO: dumux-course-task 2.2
+ //****** Select the corresponding component for the task *****//
// Uncomment first line and comment second line for using the compressible component
using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyIncompressibleComponent<Scalar> >;
// using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyCompressibleComponent<Scalar> >;
diff --git a/exercises/exercise-fluidsystem/README.md b/exercises/exercise-fluidsystem/README.md
index 0eb7631c2426f2401ec83e9bbfcf8a45843f5a39..6a5ea99f5de1d8673e5c837b153f2bd8d3958b5a 100644
--- a/exercises/exercise-fluidsystem/README.md
+++ b/exercises/exercise-fluidsystem/README.md
@@ -4,14 +4,16 @@ The aim of this exercise is to get familiar with the _DuMu<sup>x</sup>_ way of i
## Problem set-up
-The domain has a size of 60 m x 60 m and contains two low-permeable lenses. Initially, the domain is fully water saturated and the fictitious component is injected through the middle portion of the upper boundary by means of a Neumann boundary condition. The remaining parts of the upper and the entire lower boundary are Neumann no-flow while on the two lateral sides Dirichlet boundary conditions are applied (hydrostatic conditions for the pressure and zero saturation).
+The domain has a size of 60 x 60 m and contains two low-permeable lenses. Initially, the domain is fully water saturated and the fictitious component is injected through the middle portion of the upper boundary by means of a Neumann boundary condition. The remaining parts of the upper and the entire lower boundary are Neumann no-flow while on the two lateral sides Dirichlet boundary conditions are applied (hydrostatic conditions for the pressure and zero saturation).
-<img src="https://git.iws.uni-stuttgart.de/dumux-repositories/dumux-course/raw/master/exercises/extradoc/exercise-fluidsystem_setup.png" width="700">
+
-## Task 1: Getting familiar with the code
+## Preparing the exercise
-* Navigate to the directory `dumux-course/exercises/exercise-fluidsystem`.
+* Navigate to the directory `dumux-course/exercises/exercise-fluidsystem`
+
+### 1. Getting familiar with the code
Locate all the files you will need for this exercise
* The shared __main file__ : `main.cc`
@@ -24,13 +26,13 @@ Locate all the files you will need for this exercise
* The __spatial parameters file__: `spatialparams.hh`
Furthermore you will find the following folders:
-* `binarycoefficients`: Stores headers containing data/methods on binary mixtures.
-* `components`: Stores headers containing data/methods on pure components.
+* `binarycoefficients`: Stores headers containing data/methods on binary mixtures
+* `components`: Stores headers containing data/methods on pure components
* `fluidsystems`: Stores headers containing data/methods on mixtures of pure components. Uses methods from `binarycoefficients`.
-To see more components, fluid systems and binarycoefficients implementations, have a look at the folder `dumux/material`.
+To see more components, fluidsystems and binarycoefficients implementations, have a look at the folder `dumux/material`.
-## Task 2: Implement a new component - Part a
+### 2. Implement a new component
In the following, the basic steps required to set the desired fluid system are outlined. Here, this is done in the __properties file__, i.e. for this part of the exercise the code shown below is taken from the `2pproperties.hh` file.
@@ -61,14 +63,14 @@ while the `BoxModel` _TypeTag_ can be found in the `discretization/box.hh` heade
For a cell-centered scheme, you could derive from `CCTpfaModel` or `CCMpfaModel` instead (and, of course, include the right headers).
-One of the two phases should only contain the component water. We want to precompute tables on which the properties are then interpolated in order to save computational time. We need the properties of the component water in the properties file and in the problem file. We include the following two headers in our problem file, i.e. `2pproblem.hh` file. By including the problem file in the properties file, it can be accessed through it.
+One of the two phases should only contain the component water. We want to precompute tables on which the properties are then interpolated in order to save computational time. Since we need this part in both problem file and properties file, we have to include the following two headers in our problem file, i.e. `2pproblem.hh` file, and the properties file has access to them through problem file.
```c++
// The water component
#include <dumux/material/components/tabulatedcomponent.hh>
#include <dumux/material/components/h2o.hh>
```
-The other phase that will be created only contains our new component, where we want to implement an incompressible (Section 2.1) and a compressible (Section 2.2) variant.
+The other phase that will be created only contains our new component, where we want to implement an incompressible and a compressible variant.
The respective headers are prepared, but still incomplete. The compressible variant is still commented so that compilation does not fail when finishing the incompressible variant.
```c++
@@ -90,7 +92,9 @@ This fluid system expects __phases__ as input and so far we have only included t
#include <dumux/material/fluidsystems/1pliquid.hh>
```
-which creates a _liquid phase_ from a given component. Finally, using all of the included classes, we set the fluid system property. The fluid system consists of two liquids, water (derived from the tabulated water component) and the new incompressible fictitious component. Both components need to be defined as liquid (`OnePLiquid`). The two liquid components are then form an immiscible fluid system (`TwoPImmiscible`).
+which creates a _liquid phase_ from a given component. Finally, using all of the included classes we set the fluid system property by choosing that the water phase is liquid (`OnePLiquid`) and consists of the tabulated water component, and
+the other phase is liquid as well and consists of the incompressible fictitious component. Both will make up the immiscible fluid system (`TwoPImmiscible`), consisting of two liquids each consisting of one component:
+
```c++
// we use the immiscible fluid system here
@@ -102,8 +106,8 @@ private:
using TabulatedH2O = Components::TabulatedComponent<Components::H2O<Scalar>>;
using LiquidWaterPhase = typename FluidSystems::OnePLiquid<Scalar, TabulatedH2O>;
- // TODO: dumux-course-task 2.2:
- // Select the corresponding component for the task
+ // TODO: dumux-course-task 2.2
+ //****** Select the corresponding component for the task *****//
// Uncomment first line and comment second line for using the compressible component
using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyIncompressibleComponent<Scalar> >;
// using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyCompressibleComponent<Scalar> >;
@@ -113,16 +117,16 @@ public:
};
```
-### Task 2.1: Incompressible component
+### 2.1. Incompressible component
-Open the file `myincompressiblecomponent.hh`. A component should always derive from the _Base_ class - thus we include `dumux/material/components/base.hh`-, which defines the interface of a _DuMuX_ component with the possibly required functions to be overloaded by the actual implementation. Additionally, it is required that liquids are derived from the _Liquid_ class (see `dumux/material/components/liquid.hh`), gases from the _Gas_ class (see `dumux/material/components/gas.hh`) and solids from the _Solid_ class (see `dumux/material/components/solid.hh`), with functions specific to liquid, gas or solid.
+Open the file `myincompressiblecomponent.hh`. A component should always derive from the _Base_ class - thus we include `dumux/material/components/base.hh`-, which defines the interface of a _DuMuX_ component with possibly required functions to be overloaded by the actual implementation. Additionally it is required for liquids to derive from the _Liquid_ class (see `dumux/material/components/liquid.hh`), for gases to derive from the _Gas_ class (see `dumux/material/components/gas.hh`) and for solids to derive from the _Solid_ class (see `dumux/material/components/solid.hh`), with functions specific to liquid, gas or solid.
```c++
/*!
* \ingroup Components
- * \brief A ficitious component to be implemented in exercise-fluidsystem 2.1.
+ * \brief A ficitious component to be implemented in exercise 3.
*
- * \tparam Scalar The type used for scalar values.
+ * \tparam Scalar The type used for scalar values
*/
template <class Scalar>
class MyIncompressibleComponent
@@ -137,7 +141,7 @@ Implement an incompressible component into the file `myincompressiblecomponent.h
| Parameter | unit | value |
| -----| --------| -------- |
| $`M`$ | $`kg/mol`$ | $`131.39 \cdot 10^{-3}`$ |
-| $`\rho_{liquid}`$ | $`kg/m^3`$ | $`1460.0`$ |
+| $`\rho_{liquid}`$ | $`kg/m^3`$ | $`1460`$ |
| $`\mu_{liquid}`$ | $`Pa \cdot s`$ | $`5.7 \cdot 10^{-4}`$ |
In order to do so, have a look at the files `dumux/material/components/base.hh` and `dumux/material/components/liquid.hh` to see how the interfaces are defined and overload them accordingly.
@@ -150,29 +154,30 @@ make exercise_fluidsystem_a
./exercise_fluidsystem_a aparams.input
```
-The saturation distribution of the nonwetting phase, here S$`_{napl}`$ (the phase consisting of our fictitious incompressible component), at the final simulation time should look like this:
+The saturation distribution of the nonwetting phase S$`_n`$ (the phase consisting of our fictitious incompressible component) at the final simulation time should look like this:
-### Task 2.2: Compressible component
+### 2.2. Compressible component
We now want to implement a pressure-dependent density for our component. Open the file `mycompressiblecomponent.hh` and copy in the functions you implemented for the incompressible variant. Now substitute the method that returns the density by the following expression:
$`\displaystyle \rho_{MyComp} = \rho_{min} + \frac{ \rho_{max} - \rho_{min} }{ 1 + \rho_{min}*e^{-1.0*k*(\rho_{max} - \rho_{min})*p} } `$
-where $`p`$ is the pressure and $`\rho_{min} = 1440.0 \, kg/m^3`$, $`\rho_{max} = 1480.0 \, kg/m^3`$ and $`k = 5 \cdot 10^{-7} `$. Also, make sure the header is included in the `2pproperties.hh` file by uncommenting the corresponding line. Furthermore, the new component has to be set as a liquid phase in the fluid system. To do so, search for `TODO: dumux-course-task 2.2`. Comment out the corresponding line and uncomment the other. The density distribution of this phase (rho$`_{napl}`$) at the final simulation time should look like this:
+where $`p`$ is the pressure and $`\rho_{min} = 1440 `$, $`\rho_{max} = 1480 `$ and $`k = 5 \cdot 10^{-7} `$. Also, make sure the header is included in the `2pproperties.hh` file by uncommenting the corresponding line. Furthermore, the new component has to be set as a liquid phase in the fluid system. To do so, search for `TODO: dumux-course-task 2.2`. Comment out the corresponding line and uncomment the other. The density distribution of this phase (rhoN) at the final simulation time should look like this:
You can plot the density of the phase consisting of your compressible component by setting `PlotDensity` in `aparams.input` to `true` and starting the simulation again.
Compare the gnuplot output to the following plot of the density function from above:
-<img src="https://git.iws.uni-stuttgart.de/dumux-repositories/dumux-course/raw/master/exercises/extradoc/exercise-fluidsystem_a_densityfunction.png" width="500">
+
-## Task 3: Implement a new fluid system - Part b
+### 3. Implement a new fluid system
The problem file and properties file for this part of the exercise are `2p2cproblem.hh` and `2p2cproperties.hh`, respectively.
-We now want to implement a new fluid system, which still consists of two liquid phases. However, one phase consists mainly of water and the other consists mainly of the previously implemented compressible component. We will now consider compositional effects, which is why we have to derive our _TypeTag_ (`ExerciseFluidsystemTwoPTwoC`) from a _TypeTag_ (`TwoPTwoC`) that holds the miscible two-phase two-component model properties:
+We now want to implement a new fluid system, which still consists of two liquid phases. However, one phase consists mainly of water and the other consists mainly of the previously implemented compressible component. We will now consider compositional effects, which is why we have to derive our _TypeTag_ (`ExerciseFluidsystemTwoPTwoC`) from a _TypeTag_ (`TwoPTwoC`) that holds the miscible two-phase
+ two-component model properties:
```c++
// The numerical model
@@ -180,7 +185,7 @@ We now want to implement a new fluid system, which still consists of two liquid
```
```c++
-// Create new type tags
+// Create a new type tag for the problem
namespace TTag {
struct ExerciseFluidsystemTwoPTwoC { using InheritsFrom = std::tuple<TwoPTwoC, BoxModel>; };
} // end namespace TTag
@@ -208,10 +213,10 @@ public:
In the `fluidsystems/h2omycompressiblecomponent.hh` file, your implemented compressible component and the binary coefficient files are already included.
```c++
-// The ficitious component that was created in exercise-fluidsystem 2.2
+// the ficitious component that was created in exercise-fluidsystem a
#include <exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh>
-// The binary coefficients corresponding to this fluid system
+// the binary coefficients corresponding to this fluid system
#include <exercises/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh>
```
@@ -229,9 +234,9 @@ make exercise_fluidsystem_b
./exercise_fluidsystem_b bparams.input
```
-You will observe an error message and an abortion of the program. This is due to the fact that in order for the constraint solver and other mechanisms in the two-phase two-component model to work, an additional functionality in the component has to be implemented: the model has to know the vapour pressure. As in the previous exercise, check the `dumux/material/components/base.hh` file for this function and implement it into `mycompressiblecomponent.hh`. For the vapour pressure, use a value of $`3900.0 \, Pa`$.
+You will observe an error message and an abortion of the program. This is due to the fact that in order for the constraint solver and other mechanisms in the two-phase two-component model to work, an additional functionality in the component has to be implemented: the model has to know the vapour pressure. As in the previous exercise, check the `dumux/material/components/base.hh` file for this function and implement it into `mycompressiblecomponent.hh`. For the vapour pressure, use a value of $`3900`$ Pa.
-## Task 4: Change wettability of the porous medium - Part b
+### 4. Change wettability of the porous medium
In the `spatialparams.hh` file, we can find the following function, with which we can specify which phase of the fluid system is to be considered as the wetting phase at a given position within the domain:
@@ -245,13 +250,9 @@ In the `spatialparams.hh` file, we can find the following function, with which w
template<class FluidSystem>
int wettingPhaseAtPos(const GlobalPosition& globalPos) const
{
- // Our fluid system is H2OMyCompressibleComponent.
+ // Our fluid system is H2OMyCompressibleComponent
// We want to define water as the wetting phase in
- // the entire domain (see fluid system for the phase indices).
-
- // TODO: dumux-course-task 4:
- // Adapt the following line so that the phase of our new component is
- // the wetting phase, only within the lenses.
+ // the entire domain (see fluid system for the phase indices)
return FluidSystem::phase0Idx;
}
```
diff --git a/exercises/exercise-fluidsystem/aparams.input b/exercises/exercise-fluidsystem/aparams.input
index 8969d5d5c7a1123d067e3b5f788838467579c108..60897262690e489aea6e05890c9ad50d3dd32dac 100644
--- a/exercises/exercise-fluidsystem/aparams.input
+++ b/exercises/exercise-fluidsystem/aparams.input
@@ -6,11 +6,11 @@ DtInitial = 10 # initial time step size [s]
Name = exercise-fluidsystem_a # name will be given to e.g. to the vtk result files
[SpatialParams]
-BrooksCoreyPcEntry = 5.0e2 # [Pa]
+BrooksCoreyPcEntry = 5.0e2 # Pa
BrooksCoreyLambda = 2.0
Swr = 0.1
Snr = 0.0
-Lens.BrooksCoreyPcEntry = 1e3 # [Pa]
+Lens.BrooksCoreyPcEntry = 1e3 # Pa
Lens.BrooksCoreyLambda = 2.0
Lens.Swr = 0.1
Lens.Snr = 0.0
diff --git a/exercises/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh b/exercises/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh
index 79a4e55a069a1511135cb25d4d44b9caa20c4be9..dbaad539c623257d214dfc075ba2d070273d8060 100644
--- a/exercises/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh
+++ b/exercises/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh
@@ -7,7 +7,7 @@
/*!
* \file
*
- * \brief Binary coefficients for water and a fictitious component implemented in exercise-fluidsystem 2.1.
+ * \brief Binary coefficients for water and a fictitious component implemented in exercise-fluidsystem a.
*/
#ifndef DUMUX_BINARY_COEFF_H2O_MYCOMPRESSIBLECOMPONENT_HH
#define DUMUX_BINARY_COEFF_H2O_MYCOMPRESSIBLECOMPONENT_HH
@@ -18,7 +18,8 @@ namespace BinaryCoeff
{
/*!
- * \brief Binary coefficients for water and a fictitious component implemented in exercise-fluidsystem 2.1.
+ * \brief Binary coefficients for water and a fictitious component implemented in exercise-fluidsystem a
+ * The implementation of the missing methods in this file is part of exercise-fluidsystem b.
*/
class H2O_MyCompressibleComponent
{
@@ -29,7 +30,7 @@ public:
template <class Scalar>
static Scalar henryMyCompressibleComponentInWater(Scalar temperature)
{
- Scalar dumuxH = 1.5e-1 / 101.325; // [(mol/m^3)/Pa]
+ Scalar dumuxH = 1.5e-1 / 101.325; // unit [(mol/m^3)/Pa]
dumuxH *= 18.02e-6; //multiplied by molar volume of reference phase = water
return 1.0/dumuxH; // [Pa]
}
@@ -51,7 +52,7 @@ public:
static Scalar liquidDiffCoeff(Scalar temperature, Scalar pressure)
{
// arbitrary
- return 1.e-9; // [m^2/s]
+ return 1.e-9;
}
};
diff --git a/exercises/exercise-fluidsystem/bparams.input b/exercises/exercise-fluidsystem/bparams.input
index 372f314f8084a41f3c9a8889e1b72e7f1547c16c..2b6595f3ec491dffce95de2385f4b5a41860e294 100644
--- a/exercises/exercise-fluidsystem/bparams.input
+++ b/exercises/exercise-fluidsystem/bparams.input
@@ -6,11 +6,11 @@ DtInitial = 10 # initial time step size [s]
Name = exercise-fluidsystem_b # name will be given to e.g. to the vtk result files
[SpatialParams]
-BrooksCoreyPcEntry = 5.0e2 # [Pa]
+BrooksCoreyPcEntry = 5.0e2 # Pa
BrooksCoreyLambda = 2.0
Swr = 0.1
Snr = 0.0
-Lens.BrooksCoreyPcEntry = 1e3 # [Pa]
+Lens.BrooksCoreyPcEntry = 1e3 # Pa
Lens.BrooksCoreyLambda = 2.0
Lens.Swr = 0.1
Lens.Snr = 0.0
diff --git a/exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh b/exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh
index 9cce5a9843d5636c24bf680cabe6a7664904fca4..308ea171414d88a2f31da9125c06c89f2c76c52e 100644
--- a/exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh
+++ b/exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh
@@ -7,7 +7,7 @@
/*!
* \file
* \ingroup Components
- * \brief A fictitious component to be implemented in exercise 2.2.
+ * \brief A fictitious component to be implemented in exercise 3.
*/
#ifndef DUMUX_MYCOMPRESSIBLECOMPONENT_HH
#define DUMUX_MYCOMPRESSIBLECOMPONENT_HH
@@ -21,7 +21,7 @@ namespace Dumux
{
/*!
* \ingroup Components
- * \brief A fictitious component to be implemented in exercise 2.2.
+ * \brief A fictitious component to be implemented in exercise 3.
*
* \tparam Scalar The type used for scalar values
*/
@@ -39,7 +39,7 @@ public:
{ return "MyCompressibleComponent"; }
/*!
- * \brief Returns true if the liquid phase is assumed to be compressible.
+ * \brief Returns true if the liquid phase is assumed to be compressible
*/
static constexpr bool liquidIsCompressible()
{ return true; }
@@ -49,10 +49,10 @@ public:
*/
static Scalar molarMass()
{
- // TODO: dumux-course-task 2.2:
+ // TODO: dumux-course-task 2.2
// Copy the methods implemented in MyIncompressibleComponent and substitute
// the density calculation by the expression given in the exercise description.
- // Replace the line below by a meaningful return statement.
+ // Replace the line below by a meaningful return statement
DUNE_THROW(Dune::NotImplemented, "Todo: implement molar mass");
}
@@ -66,7 +66,7 @@ public:
static Scalar liquidDensity(Scalar temperature, Scalar pressure)
{
// TODO: dumux-course-task 2.2:
- // Replace the line below by a meaningful return statement.
+ // Replace the line below by a meaningful return statement
DUNE_THROW(Dune::NotImplemented, "Todo: implement liquid density");
}
@@ -92,7 +92,7 @@ public:
static Scalar liquidViscosity(Scalar temperature, Scalar pressure)
{
// TODO: dumux-course-task 2.2:
- // Replace the line below by a meaningful return statement.
+ // Replace the line below by a meaningful return statement
DUNE_THROW(Dune::NotImplemented, "Todo: implement liquid viscosity");
}
@@ -104,8 +104,8 @@ public:
*/
static Scalar vaporPressure(Scalar t)
{
- // TODO: dumux-course-task 3:
- // Replace the line below by a meaningful return statement.
+ // TODO: dumux-course-task 3
+ // Replace the line below by a meaningful return statement
DUNE_THROW(Dune::NotImplemented, "Todo: implement vapour pressure");
}
};
diff --git a/exercises/exercise-fluidsystem/components/myincompressiblecomponent.hh b/exercises/exercise-fluidsystem/components/myincompressiblecomponent.hh
index 3ce865687a0497721b13f165e7487ff43a470dbe..3a5c0f4cd13d32931dd76cbd2ce8aa42e8f4407f 100644
--- a/exercises/exercise-fluidsystem/components/myincompressiblecomponent.hh
+++ b/exercises/exercise-fluidsystem/components/myincompressiblecomponent.hh
@@ -7,7 +7,7 @@
/*!
* \file
* \ingroup Components
- * \brief A fictitious component to be implemented in exercise-fluidsystem 2.1.
+ * \brief A fictitious component to be implemented in exercise-fluidsystem.
*/
#ifndef DUMUX_MYINCOMPRESSIBLECOMPONENT_HH
#define DUMUX_MYINCOMPRESSIBLECOMPONENT_HH
@@ -21,9 +21,9 @@ namespace Dumux
{
/*!
* \ingroup Components
- * \brief A fictitious component to be implemented in exercise-fluidsystem 2.1.
+ * \brief A fictitious component to be implemented in exercise-fluidsystem.
*
- * \tparam Scalar The type used for scalar values.
+ * \tparam Scalar The type used for scalar values
*/
template <class Scalar>
class MyIncompressibleComponent
@@ -38,7 +38,7 @@ public:
{ return "MyIncompressibleComponent"; }
/*!
- * \brief Returns true if the liquid phase is assumed to be compressible.
+ * \brief Returns true if the liquid phase is assumed to be compressible
*/
static constexpr bool liquidIsCompressible()
{ return false; }
@@ -48,9 +48,9 @@ public:
*/
static Scalar molarMass()
{
- // TODO: dumux-course-task 2.1:
+ // TODO: dumux-course-task 2.1
// Implement the methods for the component data given in the exercise description.
- // Replace the line below by a meaningful return statement.
+ // Replace the line below by a meaningful return statement
DUNE_THROW(Dune::NotImplemented, "Todo: implement molarMass()");
}
@@ -63,9 +63,9 @@ public:
*/
static Scalar liquidDensity(Scalar temperature, Scalar pressure)
{
- // TODO: dumux-course-task 2.1:
+ // TODO: dumux-course-task 2.1
// Implement the methods for the component data given in the exercise description.
- // Replace the line below by a meaningful return statement.
+ // Replace the line below by a meaningful return statement
DUNE_THROW(Dune::NotImplemented, "Todo: implement liquidDensity()");
}
@@ -90,9 +90,9 @@ public:
*/
static Scalar liquidViscosity(Scalar temperature, Scalar pressure)
{
- // TODO: dumux-course-task 2.1:
+ // TODO: dumux-course-task 2.1
// Implement the methods for the component data given in the exercise description.
- // Replace the line below by a meaningful return statement.
+ // Replace the line below by a meaningful return statement
DUNE_THROW(Dune::NotImplemented, "Todo: implement liquidViscosity()");
}
};
diff --git a/exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh b/exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh
index 14490adf5a243817386227c4cbb45ba90cb266cd..bbed8af548ce429509948fb3f9a3ea14ec5740ad 100644
--- a/exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh
+++ b/exercises/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh
@@ -8,8 +8,8 @@
* \file
*
* \brief A fluid system with water and a fictitious component, which is to be
- * implemented in exercise-fluidsystem 2.2, as phases and components. This
- * fluid system is to be implemented in exercise-fluidsystem 3.
+ * implemented in exercise-fluidsystem a, as phases and components. This
+ * fluid system is to be implemented in exercise-fluidsystem b.
*/
#ifndef DUMUX_H2O_MYCOMPRESSIBLECOMPONENT_FLUID_SYSTEM_HH
#define DUMUX_H2O_MYCOMPRESSIBLECOMPONENT_FLUID_SYSTEM_HH
@@ -19,10 +19,10 @@
#include <dumux/material/fluidsystems/base.hh>
-// The fictitious component that was created in exercise-fluidsystem 2.2
+// the fictitious component that was created in exercise-fluidsystem a
#include <exercises/exercise-fluidsystem/components/mycompressiblecomponent.hh>
-// The binary coefficients corresponding to this fluid system
+// the binary coefficients corresponding to this fluid system
#include <exercises/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh>
namespace Dumux {
@@ -30,7 +30,7 @@ namespace FluidSystems {
/*!
* \brief A compositional fluid consisting of two liquid phases,
- * which are water and a fictitious component from exercise-fluidsystem 2.2.
+ * which are water and a fictitious component from exercise-fluidsystem a.
*/
template <class Scalar,
class H2OType = Dumux::Components::TabulatedComponent<Dumux::Components::H2O<Scalar> > >
@@ -47,12 +47,12 @@ public:
static constexpr int numPhases = 2;
static constexpr int numComponents = 2;
- static constexpr int phase0Idx = 0; // Index of the first phase
- static constexpr int phase1Idx = 1; // Index of the second phase
+ static constexpr int phase0Idx = 0; // index of the first phase
+ static constexpr int phase1Idx = 1; // index of the second phase
static constexpr int H2OIdx = 0;
static constexpr int NAPLIdx = 1;
- // Export component indices to indicate the main component
+ // export component indices to indicate the main component
// of the corresponding phase at atmospheric pressure 1 bar
// and room temperature 20°C:
static constexpr int comp0Idx = H2OIdx;
@@ -181,7 +181,6 @@ public:
DUNE_THROW(Dune::InvalidStateException, "Invalid component index " << compIdx);
}
-
/*!
* \brief Return the molar mass of a component in [kg/mol].
*/
@@ -209,7 +208,7 @@ public:
if (phaseIdx == phase0Idx)
{
- // see: doctoral thesis of Steffen Ochs 2007
+ // See: doctoral thesis of Steffen Ochs 2007
// Steam injection into saturated porous media : process analysis including experimental and numerical investigations
// http://elib.uni-stuttgart.de/bitstream/11682/271/1/Diss_Ochs_OPUS.pdf
@@ -218,8 +217,8 @@ public:
// Scalar x_H2O = fluidState.moleFraction(phase0Idx, H2OIdx);
// Scalar x_myComp = fluidState.moleFraction(phase0Idx, NAPLIdx);
- // TODO: dumux-course-task 3:
- // implement the composition-dependent water density from the exercise sheet.
+ // TODO: dumux-course-task 3
+ // Implement the composition-dependent water density from the exercise sheet.
DUNE_THROW(Dune::NotImplemented, "Todo: implement composition-dependent density");
}
@@ -238,7 +237,7 @@ public:
* The molar density for the simple relation is defined by the
* mass density \f$\rho_\alpha\f$ and the molar mass of the main component
*
- * The molar density for the complex relation is defined by the
+ * The molar density for the complrex relation is defined by the
* mass density \f$\rho_\alpha\f$ and the mean molar mass \f$\overline M_\alpha\f$:
*
* \f[\rho_{mol,\alpha} = \frac{\rho_\alpha}{\overline M_\alpha} \;.\f]
@@ -396,8 +395,7 @@ public:
DUNE_THROW(Dune::NotImplemented, "FluidSystems::H2OMyCompressibleComponent::kelvinVaporPressure()");
}
- /*!
- * \brief Partial pressures in the gas phase, taken from saturation vapor pressures.
+ /* partial pressures in the gas phase, taken from saturation vapor pressures
*/
template <class FluidState>
static Scalar partialPressureGas(const FluidState &fluidState, int phaseIdx,
@@ -414,8 +412,7 @@ public:
DUNE_THROW(Dune::InvalidStateException, "non-existent component index " << compIdx);
}
- /*!
- * \brief Inverse vapor pressures, taken from inverse saturation vapor pressures.
+ /* inverse vapor pressures, taken from inverse saturation vapor pressures
*/
template <class FluidState>
static Scalar inverseVaporPressureCurve(const FluidState &fluidState,
diff --git a/exercises/exercise-fluidsystem/main.cc b/exercises/exercise-fluidsystem/main.cc
index 5f8394b7c4c15968247664cfecd8fa57c2364765..741f0621acdc27bf102037b6887c441e097e6de4 100644
--- a/exercises/exercise-fluidsystem/main.cc
+++ b/exercises/exercise-fluidsystem/main.cc
@@ -6,7 +6,7 @@
//
/*!
* \file
- * \brief The main file for exercise-fluidsystem. This solves an instationary problem
+ * \brief The main file for exercise 3. This solves an instationary problem
* with an implicit time discretization.
*/
#include <config.h>
@@ -37,11 +37,11 @@ int main(int argc, char** argv)
{
using namespace Dumux;
- // TODO: dumux-course-task 3:
- // TYPETAG is set in CMakeLists.txt as compile time definition.
- // Alternatively, you could write `using TypeTag = Properties::TTag::ExerciseFluidsystemTwoP;`
- // then, for the 2p2c problem you have to change this line to `using TypeTag = Properties::TTag::ExerciseFluidsystemTwoPTwoC;`
- // and recompile the executable.
+ // define the type tag for this problem
+ // TYPETAG is set in CMakeLists.txt as compile time definition
+ // alternatively you could write `using TypeTag = Properties::TTag::ExerciseFluidsystemBoxTwoP;`
+ // then, for the 2p2c problem you have to change this line to `using TypeTag = Properties::TTag::ExerciseFluidsystemBoxTwoPTwoC;`
+ // and recompile the executable
using TypeTag = Properties::TTag::TYPETAG;
// initialize MPI+x, finalize is done automatically on exit
diff --git a/exercises/exercise-fluidsystem/spatialparams.hh b/exercises/exercise-fluidsystem/spatialparams.hh
index 11c74faa4c67390109872cc7e9423da4f96ec233..7601c3af2b1f91f43e07c9228c1e6dc6fb23e602 100644
--- a/exercises/exercise-fluidsystem/spatialparams.hh
+++ b/exercises/exercise-fluidsystem/spatialparams.hh
@@ -13,10 +13,10 @@
#ifndef DUMUX_EXERCISE_FLUIDSYSTEM_SPATIAL_PARAMS_HH
#define DUMUX_EXERCISE_FLUIDSYSTEM_SPATIAL_PARAMS_HH
-// Include parent spatial parameters
+// include parent spatialparameters
#include <dumux/porousmediumflow/fvspatialparamsmp.hh>
-// Include fluid-matrix interaction relationship
+//include fluid matrix interaction relationship
#include <dumux/material/fluidmatrixinteractions/2p/brookscorey.hh>
#include <dumux/material/fluidmatrixinteractions/2p/linearmaterial.hh>
@@ -44,13 +44,13 @@ class ExerciseFluidsystemSpatialParams
using PcKrSwCurve = FluidMatrix::BrooksCoreyDefault<Scalar>;
public:
- // Export permeability type
+ // export permeability type
using PermeabilityType = Dune::FieldMatrix<Scalar, dim, dim>;
/*!
- * \brief The constructor.
+ * \brief The constructor
*
- * \param gridGeometry The finite volume grid geometry.
+ * \param gridGeometry The finite volume grid geometry
*/
ExerciseFluidsystemSpatialParams(std::shared_ptr<const GridGeometry>& gridGeometry)
: ParentType(gridGeometry)
@@ -59,8 +59,8 @@ public:
, pcKrSwCurve_("SpatialParams")
, lensPcKrSwCurve_("SpatialParams.Lens")
{
- // set main diagonal entries of the permeability tensor to a value
- // setting to one value means: isotropic, homogeneous
+ //set main diagonal entries of the permeability tensor to a value
+ //setting to one value means: isotropic, homogeneous
for (int i = 0; i < dim; i++)
{
K_[i][i] = 1e-7;
@@ -71,7 +71,7 @@ public:
/*!
* \brief Define the intrinsic permeability \f$\mathrm{[m^2]}\f$.
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
PermeabilityType permeabilityAtPos(const GlobalPosition& globalPos) const
@@ -84,7 +84,7 @@ public:
/*!
* \brief Define the porosity \f$\mathrm{[-]}\f$.
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
Scalar porosityAtPos(const GlobalPosition& globalPos) const
{
@@ -94,8 +94,8 @@ public:
}
/*!
- * \brief Returns the fluid-matrix interaction law at a given location.
- * \param globalPos The global coordinates for the given location.
+ * \brief Returns the fluid-matrix interaction law at a given location
+ * \param globalPos The global coordinates for the given location
*/
auto fluidMatrixInteractionAtPos(const GlobalPosition& globalPos) const
{
@@ -107,17 +107,17 @@ public:
/*!
* \brief Function for defining which phase is to be considered as the wetting phase.
*
- * \return The wetting phase index.
- * \param globalPos The position of the center of the element.
+ * \return the wetting phase index
+ * \param globalPos The position of the center of the element
*/
template<class FluidSystem>
int wettingPhaseAtPos(const GlobalPosition& globalPos) const
{
- // Our fluid system is H2OMyCompressibleComponent.
+ // Our fluid system is H2OMyCompressibleComponent
// We want to define water as the wetting phase in
- // the entire domain (see fluid system for the phase indices).
+ // the entire domain (see fluid system for the phase indices)
- // TODO: dumux-course-task 4:
+ // TODO: dumux-course-task 4
// Adapt the following line so that the phase of our new component is
// the wetting phase, only within the lenses.
return FluidSystem::phase0Idx;
@@ -133,8 +133,8 @@ public:
}
/*!
- * \brief Returns the temperature at the domain at the given position.
- * \param globalPos The position in global coordinates where the temperature should be specified.
+ * \brief Returns the temperature at the domain at the given position
+ * \param globalPos The position in global coordinates where the temperature should be specified
*/
Scalar temperatureAtPos(const GlobalPosition& globalPos) const
{
@@ -147,7 +147,7 @@ private:
Dune::FieldMatrix<Scalar, dim, dim> K_;
Dune::FieldMatrix<Scalar, dim, dim> KLens_;
- // Object that holds the values/parameters of the selected fluid matrix interaction relationship.
+ // Object that holds the values/parameters of the selected fluid matrix interaction relationship
const PcKrSwCurve pcKrSwCurve_;
const PcKrSwCurve lensPcKrSwCurve_;
};
diff --git a/exercises/extradoc/exercise-fluidsystem_a_solution.png b/exercises/extradoc/exercise-fluidsystem_a_solution.png
index 71a3634f5a5c09646d62ec091ae72ec4e5941cd8..fdf3ef23c81a0bc414554a9bc2e4f1e784e134c4 100644
Binary files a/exercises/extradoc/exercise-fluidsystem_a_solution.png and b/exercises/extradoc/exercise-fluidsystem_a_solution.png differ
diff --git a/exercises/extradoc/exercise-fluidsystem_a_solution2.png b/exercises/extradoc/exercise-fluidsystem_a_solution2.png
index 4004196f000c2de22328d38ed8637bd154616e20..ba90e6fcfa783fcc5a4eaf1c5fbe252dc46ea8d9 100644
Binary files a/exercises/extradoc/exercise-fluidsystem_a_solution2.png and b/exercises/extradoc/exercise-fluidsystem_a_solution2.png differ
diff --git a/exercises/solution/exercise-basic/injection2pniproblem.hh b/exercises/solution/exercise-basic/injection2pniproblem.hh
index a0e7478b5717545f7f4a940bdf83ebbfe8a9f461..aa74091c152dd459ea9ab9fb8c38e618d7f60e4b 100644
--- a/exercises/solution/exercise-basic/injection2pniproblem.hh
+++ b/exercises/solution/exercise-basic/injection2pniproblem.hh
@@ -73,13 +73,13 @@ public:
/*pMax=*/30e6,
/*numP=*/300);
- // Name of the problem and output file
+ // name of the problem and output file
// getParam<TYPE>("GROUPNAME.PARAMNAME") reads and sets parameter PARAMNAME
// of type TYPE given in the group GROUPNAME from the input file
name_ = getParam<std::string>("Problem.Name");
- // Depth of the aquifer, unit: m
+ // depth of the aquifer, unit: m
aquiferDepth_ = getParam<Scalar>("Problem.AquiferDepth");
- // The duration of the injection, unit: seconds
+ // the duration of the injection, unit: seconds
injectionDuration_ = getParam<Scalar>("Problem.InjectionDuration");
}
@@ -127,20 +127,20 @@ public:
*/
NumEqVector neumannAtPos(const GlobalPosition &globalPos) const
{
- // Initialize values to zero, i.e. no-flow Neumann boundary conditions
+ // initialize values to zero, i.e. no-flow Neumann boundary conditions
NumEqVector values(0.0);
- // If we are inside the injection zone set inflow Neumann boundary conditions
- if (injectionActive() && onInjectionBoundary(globalPos))
+ // if we are inside the injection zone set inflow Neumann boundary conditions
+ if (injectionActive() && onInjectionBoundary(globalPos))
{
const Scalar injectionRate = -1e-4;
- // Inject nitrogen. Negative values mean injection
- // Unit: kg/(s*m^2)
+ // inject nitrogen. Negative values mean injection
+ // unit: kg/(s*m^2)
values[Indices::conti0EqIdx + FluidSystem::N2Idx] = injectionRate;
values[Indices::conti0EqIdx + FluidSystem::H2OIdx] = 0.0;
- // Energy fluxes are always mass specific
+ // energy fluxes are always mass specific
// unit: W/(m^2)
const Scalar temperatureAtInjection = initialAtPos(globalPos)[Indices::temperatureIdx];/*K*/
const Scalar pressureAtInjection = initialAtPos(globalPos)[Indices::pressureIdx];/*Pa*/
@@ -170,11 +170,11 @@ public:
{
PrimaryVariables values(0.0);
- // Get the water density at atmospheric conditions
+ // get the water density at atmospheric conditions
const Scalar densityW = FluidSystem::H2O::liquidDensity(283.15, 1.0e5);
- // Assume an initially hydrostatic liquid pressure profile
- // Note: we subtract rho_w*g*h because g is defined negative
+ // assume an initially hydrostatic liquid pressure profile
+ // note: we subtract rho_w*g*h because g is defined negative
const Scalar pw = 1.0e5 - densityW*this->spatialParams().gravity(globalPos)[dimWorld-1]*(aquiferDepth_ - globalPos[dimWorld-1]);
values[Indices::pressureIdx] = pw;
diff --git a/exercises/solution/exercise-fluidsystem/2p2cproblem.hh b/exercises/solution/exercise-fluidsystem/2p2cproblem.hh
index 4f16f28fef60cdaa5cf546e1e536704d1978913e..2cf745904c93c067b052bc4b90bd82abe5955a25 100644
--- a/exercises/solution/exercise-fluidsystem/2p2cproblem.hh
+++ b/exercises/solution/exercise-fluidsystem/2p2cproblem.hh
@@ -7,12 +7,12 @@
/*!
* \file
*
- * \brief Tutorial problem for a fully coupled two-phase two-component box model.
+ * \brief Tutorial problem for a fully coupled two phase-two component box model.
*/
#ifndef DUMUX_EXERCISE_FLUIDSYSTEM_B_PROBLEM_HH
#define DUMUX_EXERCISE_FLUIDSYSTEM_B_PROBLEM_HH
-// The base porous media box problem.
+// The base porous media box problem
#include <dumux/porousmediumflow/problem.hh>
#include <dumux/common/boundarytypes.hh>
#include <dumux/common/properties.hh>
@@ -23,7 +23,7 @@ namespace Dumux {
/*!
* \ingroup TwoPBoxModel
*
- * \brief Tutorial problem for a fully coupled two-phase two-component box model.
+ * \brief Tutorial problem for a fully coupled two phase-two component box model.
*/
template <class TypeTag>
class ExerciseFluidsystemProblemTwoPTwoC : public PorousMediumFlowProblem<TypeTag>
@@ -64,7 +64,7 @@ public:
* \brief Specifies which kind of boundary condition should be
* used for which equation on a given boundary segment.
*
- * \param globalPos The position for which the bc type should be evaluated.
+ * \param globalPos The position for which the bc type should be evaluated
*/
BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const
{
@@ -80,9 +80,9 @@ public:
/*!
* \brief Evaluates the boundary conditions for a Dirichlet
- * boundary segment.
+ * boundary segment
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
PrimaryVariables dirichletAtPos(const GlobalPosition &globalPos) const
{
@@ -113,7 +113,7 @@ public:
}
else
{
- // no-flow on the remaining Neumann-boundaries
+ // no-flow on the remaining Neumann-boundaries.
values[Indices::conti0EqIdx + FluidSystem::H2OIdx] = 0;
values[Indices::conti0EqIdx + FluidSystem::NAPLIdx] = 0;
}
@@ -125,7 +125,7 @@ public:
/*!
* \brief Evaluate the initial value for a control volume.
*
- * \param globalPos The position for which the initial condition should be evaluated.
+ * \param globalPos The position for which the initial condition should be evaluated
*
* For this method, the \a values parameter stores primary
* variables.
@@ -149,9 +149,9 @@ public:
/*!
- * \brief Returns the source term.
+ * \brief Returns the source term
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
{
@@ -162,8 +162,8 @@ public:
private:
- Scalar eps_; //! Small epsilon value
- Scalar depthBOR_; //! Depth at the bottom of the reservoir
+ Scalar eps_; //! small epsilon value
+ Scalar depthBOR_; //! depth at the bottom of the reservoir
};
} // end namespace Dumux
diff --git a/exercises/solution/exercise-fluidsystem/2p2cproperties.hh b/exercises/solution/exercise-fluidsystem/2p2cproperties.hh
index fc13a8a68f833c35690ee8fddeba36d6b3cf43ed..5bdc04c72bbc618ddbdd1ef4df3696f56b7c2e2d 100644
--- a/exercises/solution/exercise-fluidsystem/2p2cproperties.hh
+++ b/exercises/solution/exercise-fluidsystem/2p2cproperties.hh
@@ -27,11 +27,12 @@
// The fluid system that is created in this exercise
#include "fluidsystems/h2omycompressiblecomponent.hh"
-// The problem file, where setup-specific boundary and initial conditions are defined
+// The problem file, where setup-specific boundary and initial conditions are defined.
#include"2p2cproblem.hh"
namespace Dumux::Properties {
+// Create a new type tag for the problem
// Create new type tags
namespace TTag {
struct ExerciseFluidsystemTwoPTwoC { using InheritsFrom = std::tuple<TwoPTwoC, BoxModel>; };
diff --git a/exercises/solution/exercise-fluidsystem/2pproblem.hh b/exercises/solution/exercise-fluidsystem/2pproblem.hh
index 4f9d17b2b1396d61ff90dbf7d8050452597f1ae5..8258644f199076f13dd7500364218ebde2ed5051 100644
--- a/exercises/solution/exercise-fluidsystem/2pproblem.hh
+++ b/exercises/solution/exercise-fluidsystem/2pproblem.hh
@@ -7,7 +7,7 @@
/*!
* \file
*
- * \brief Tutorial problem for a fully coupled two-phase box model.
+ * \brief Tutorial problem for a fully coupled twophase box model.
*/
#ifndef DUMUX_EXERCISE_FLUIDSYSTEM_A_PROBLEM_HH
#define DUMUX_EXERCISE_FLUIDSYSTEM_A_PROBLEM_HH
@@ -29,7 +29,7 @@ namespace Dumux {
/*!
* \ingroup TwoPBoxModel
- * \brief Tutorial problem for a fully coupled two-phase box model.
+ * \brief Tutorial problem for a fully coupled twophase box model.
*/
template <class TypeTag>
class ExerciseFluidsystemProblemTwoP : public PorousMediumFlowProblem<TypeTag>
@@ -86,7 +86,7 @@ public:
* \brief Specifies which kind of boundary condition should be
* used for which equation on a given boundary segment.
*
- * \param globalPos The position for which the bc type should be evaluated.
+ * \param globalPos The position for which the bc type should be evaluated
*/
BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const
{
@@ -102,9 +102,9 @@ public:
/*!
* \brief Evaluates the boundary conditions for a Dirichlet
- * boundary segment.
+ * boundary segment
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
PrimaryVariables dirichletAtPos(const GlobalPosition& globalPos) const
{
@@ -136,7 +136,7 @@ public:
}
else
{
- // no-flow on the remaining Neumann-boundaries
+ // no-flow on the remaining Neumann-boundaries.
values[contiWEqIdx] = 0;
values[contiNEqIdx] = 0;
}
@@ -148,7 +148,7 @@ public:
/*!
* \brief Evaluate the initial value for a control volume.
*
- * \param globalPos The position for which the initial condition should be evaluated.
+ * \param globalPos The position for which the initial condition should be evaluated
*
* For this method, the \a values parameter stores primary
* variables.
@@ -167,9 +167,9 @@ public:
/*!
- * \brief Returns the source term.
+ * \brief Returns the source term
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
NumEqVector sourceAtPos(const GlobalPosition& globalPos) const
{
@@ -207,9 +207,9 @@ private:
gnuplot_.plot("YourComponentPhase_density");
}
- Scalar eps_; //! Small epsilon value
- Scalar depthBOR_; //! Depth at the bottom of the reservoir
- Dumux::GnuplotInterface<double> gnuplot_; //! Collects data for plotting
+ Scalar eps_; //! small epsilon value
+ Scalar depthBOR_; //! depth at the bottom of the reservoir
+ Dumux::GnuplotInterface<double> gnuplot_; //! collects data for plotting
};
} // end namespace Dumux
diff --git a/exercises/solution/exercise-fluidsystem/2pproperties.hh b/exercises/solution/exercise-fluidsystem/2pproperties.hh
index dd901d72a65904907ab2964a3cbe4793a8476406..a89db8fd209e715acb811d91411249cdff067da9 100644
--- a/exercises/solution/exercise-fluidsystem/2pproperties.hh
+++ b/exercises/solution/exercise-fluidsystem/2pproperties.hh
@@ -7,7 +7,7 @@
/*!
* \file
*
- * \brief The properties file for two-phase exercise-fluidsystem.
+ * \brief The properties file for two phase exercise-fluidsystem.
*/
#ifndef DUMUX_EXERCISE_FLUIDSYSTEM_A_PROPERTIES_HH
#define DUMUX_EXERCISE_FLUIDSYSTEM_A_PROPERTIES_HH
@@ -37,11 +37,12 @@
// The interface to create plots during simulation
#include <dumux/io/gnuplotinterface.hh>
-// The problem file, where setup-specific boundary and initial conditions are defined
+// The problem file, where setup-specific boundary and initial conditions are defined.
#include"2pproblem.hh"
namespace Dumux::Properties {
+// Create a new type tag for the problem
// Create new type tags
namespace TTag {
struct ExerciseFluidsystemTwoP { using InheritsFrom = std::tuple<TwoP, BoxModel>; };
@@ -66,7 +67,7 @@ public:
template<class TypeTag>
struct Grid<TypeTag, TTag::ExerciseFluidsystemTwoP> { using type = Dune::YaspGrid<2>; };
-// We use the immiscible fluid system here
+// we use the immiscible fluid system here
template<class TypeTag>
struct FluidSystem<TypeTag, TTag::ExerciseFluidsystemTwoP>
{
@@ -75,12 +76,11 @@ private:
using TabulatedH2O = Components::TabulatedComponent<Components::H2O<Scalar>>;
using LiquidWaterPhase = typename FluidSystems::OnePLiquid<Scalar, TabulatedH2O>;
/*!
- * TODO: dumux-course-task 2.2:
* Uncomment first line and comment second line for using the incompressible component
* Uncomment second line and comment first line for using the compressible component
*/
- using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyIncompressibleComponent<Scalar> >;
- //using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyCompressibleComponent<Scalar> >;
+ using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyIncompressibleComponent<Scalar> >;
+ // using LiquidMyComponentPhase = typename FluidSystems::OnePLiquid<Scalar, MyCompressibleComponent<Scalar> >;
public:
using type = typename FluidSystems::TwoPImmiscible<Scalar, LiquidWaterPhase, LiquidMyComponentPhase>;
diff --git a/exercises/solution/exercise-fluidsystem/aparams.input b/exercises/solution/exercise-fluidsystem/aparams.input
index 5add5c761dcd9881967186533f369d80967b2304..594ad15c015cad4b0d28d64acc4475bec490382f 100644
--- a/exercises/solution/exercise-fluidsystem/aparams.input
+++ b/exercises/solution/exercise-fluidsystem/aparams.input
@@ -6,11 +6,11 @@ DtInitial = 10 # initial time step size [s]
Name = exercise-fluidsystem_a_solution # name will be given to e.g. to the vtk result files
[SpatialParams]
-BrooksCoreyPcEntry = 5.0e2 # [Pa]
+BrooksCoreyPcEntry = 5.0e2 # Pa
BrooksCoreyLambda = 2.0
Swr = 0.1
Snr = 0.0
-Lens.BrooksCoreyPcEntry = 1e3 # [Pa]
+Lens.BrooksCoreyPcEntry = 1e3 # Pa
Lens.BrooksCoreyLambda = 2.0
Lens.Swr = 0.1
Lens.Snr = 0.0
diff --git a/exercises/solution/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh b/exercises/solution/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh
index 79a4e55a069a1511135cb25d4d44b9caa20c4be9..dbaad539c623257d214dfc075ba2d070273d8060 100644
--- a/exercises/solution/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh
+++ b/exercises/solution/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh
@@ -7,7 +7,7 @@
/*!
* \file
*
- * \brief Binary coefficients for water and a fictitious component implemented in exercise-fluidsystem 2.1.
+ * \brief Binary coefficients for water and a fictitious component implemented in exercise-fluidsystem a.
*/
#ifndef DUMUX_BINARY_COEFF_H2O_MYCOMPRESSIBLECOMPONENT_HH
#define DUMUX_BINARY_COEFF_H2O_MYCOMPRESSIBLECOMPONENT_HH
@@ -18,7 +18,8 @@ namespace BinaryCoeff
{
/*!
- * \brief Binary coefficients for water and a fictitious component implemented in exercise-fluidsystem 2.1.
+ * \brief Binary coefficients for water and a fictitious component implemented in exercise-fluidsystem a
+ * The implementation of the missing methods in this file is part of exercise-fluidsystem b.
*/
class H2O_MyCompressibleComponent
{
@@ -29,7 +30,7 @@ public:
template <class Scalar>
static Scalar henryMyCompressibleComponentInWater(Scalar temperature)
{
- Scalar dumuxH = 1.5e-1 / 101.325; // [(mol/m^3)/Pa]
+ Scalar dumuxH = 1.5e-1 / 101.325; // unit [(mol/m^3)/Pa]
dumuxH *= 18.02e-6; //multiplied by molar volume of reference phase = water
return 1.0/dumuxH; // [Pa]
}
@@ -51,7 +52,7 @@ public:
static Scalar liquidDiffCoeff(Scalar temperature, Scalar pressure)
{
// arbitrary
- return 1.e-9; // [m^2/s]
+ return 1.e-9;
}
};
diff --git a/exercises/solution/exercise-fluidsystem/bparams.input b/exercises/solution/exercise-fluidsystem/bparams.input
index 3e73a5ead825b33de6f747b2c690a0d155403631..a421776f1b1bbbdd0976dae7d779f1e716016019 100644
--- a/exercises/solution/exercise-fluidsystem/bparams.input
+++ b/exercises/solution/exercise-fluidsystem/bparams.input
@@ -6,11 +6,11 @@ DtInitial = 10 # initial time step size [s]
Name = exercise-fluidsystem_b_solution # name will be given to e.g. to the vtk result files
[SpatialParams]
-BrooksCoreyPcEntry = 5.0e2 # [Pa]
+BrooksCoreyPcEntry = 5.0e2 # Pa
BrooksCoreyLambda = 2.0
Swr = 0.1
Snr = 0.0
-Lens.BrooksCoreyPcEntry = 1e3 # [Pa]
+Lens.BrooksCoreyPcEntry = 1e3 # Pa
Lens.BrooksCoreyLambda = 2.0
Lens.Swr = 0.1
Lens.Snr = 0.0
diff --git a/exercises/solution/exercise-fluidsystem/components/mycompressiblecomponent.hh b/exercises/solution/exercise-fluidsystem/components/mycompressiblecomponent.hh
index 1c6011b315a34820973c50845f0df2b791086883..2171bd1af543e0567de4ef37f4d846c6445f717f 100644
--- a/exercises/solution/exercise-fluidsystem/components/mycompressiblecomponent.hh
+++ b/exercises/solution/exercise-fluidsystem/components/mycompressiblecomponent.hh
@@ -7,7 +7,7 @@
/*!
* \file
* \ingroup Components
- * \brief A fictitious component to be implemented in exercise 2.2.
+ * \brief A fictitious component to be implemented in exercise 3.
*/
#ifndef DUMUX_MYCOMPRESSIBLECOMPONENT_HH
#define DUMUX_MYCOMPRESSIBLECOMPONENT_HH
@@ -21,7 +21,7 @@ namespace Dumux
{
/*!
* \ingroup Components
- * \brief A fictitious component to be implemented in exercise 2.2.
+ * \brief A fictitious component to be implemented in exercise 3.
*
* \tparam Scalar The type used for scalar values
*/
@@ -39,7 +39,7 @@ public:
{ return "MyCompressibleComponent"; }
/*!
- * \brief Returns true if the liquid phase is assumed to be compressible.
+ * \brief Returns true if the liquid phase is assumed to be compressible
*/
static constexpr bool liquidIsCompressible()
{ return true; }
@@ -61,9 +61,9 @@ public:
*/
static Scalar liquidDensity(Scalar temperature, Scalar pressure)
{
- static const Scalar rho_min = 1440.0;
- static const Scalar rho_max = 1480.0;
- static const Scalar k = 5.0e-7;
+ static const Scalar rho_min = 1440;
+ static const Scalar rho_max = 1480;
+ static const Scalar k = 5e-7;
using std::exp;
return rho_min + (rho_max - rho_min)/(1 + rho_min*exp(-1.0*k*(rho_max - rho_min)*pressure)); // [kg/m^3]
@@ -101,7 +101,7 @@ public:
*/
static Scalar vaporPressure(Scalar t)
{
- return 3900.0; // [Pa] (at 20°C)
+ return 3900; // [Pa] (at 20C)
}
};
diff --git a/exercises/solution/exercise-fluidsystem/components/myincompressiblecomponent.hh b/exercises/solution/exercise-fluidsystem/components/myincompressiblecomponent.hh
index 6bca0349fbd58112010064ea323c3d43744c39cd..40f1989831a61c1ad1aa7960abb2df958a032089 100644
--- a/exercises/solution/exercise-fluidsystem/components/myincompressiblecomponent.hh
+++ b/exercises/solution/exercise-fluidsystem/components/myincompressiblecomponent.hh
@@ -7,7 +7,7 @@
/*!
* \file
* \ingroup Components
- * \brief A fictitious component to be implemented in exercise-fluidsystem 2.1.
+ * \brief A fictitious component to be implemented in exercise-fluidsystem a.
*/
#ifndef DUMUX_MYINCOMPRESSIBLECOMPONENT_HH
#define DUMUX_MYINCOMPRESSIBLECOMPONENT_HH
@@ -21,9 +21,9 @@ namespace Dumux
{
/*!
* \ingroup Components
- * \brief A fictitious component to be implemented in exercise-fluidsystem 2.1.
+ * \brief A fictitious component to be implemented in exercise-fluidsystem a.
*
- * \tparam Scalar The type used for scalar values.
+ * \tparam Scalar The type used for scalar values
*/
template <class Scalar>
class MyIncompressibleComponent
@@ -38,7 +38,7 @@ public:
{ return "MyIncompressibleComponent"; }
/*!
- * \brief Returns true if the liquid phase is assumed to be compressible.
+ * \brief Returns true if the liquid phase is assumed to be compressible
*/
static constexpr bool liquidIsCompressible()
{ return false; }
diff --git a/exercises/solution/exercise-fluidsystem/components/plotdensityfunction.py b/exercises/solution/exercise-fluidsystem/components/plotdensityfunction.py
index 943bfafe5c4d07bd0b78a2dc0e366635ef1671e0..209e191b01dd7bc25a2ff40ff36a78a608b3f9eb 100644
--- a/exercises/solution/exercise-fluidsystem/components/plotdensityfunction.py
+++ b/exercises/solution/exercise-fluidsystem/components/plotdensityfunction.py
@@ -5,18 +5,18 @@
import numpy as np
import matplotlib.pyplot as plt
-# Function to calculate rho dependent on pressure.
-rho_min = 1440.0;
-rho_max = 1480.0;
-k = 5.0e-7;
+# function to calculate rho dependent on pressure
+rho_min = 1440;
+rho_max = 1480;
+k = 5e-7;
def rho(p):
return rho_min + (rho_max - rho_min)/(1 + rho_min*np.exp(-1.0*k*(rho_max - rho_min)*p));
-# Sample pressure in range (1e4, 1e7) and compute corresponding densities.
+# sample pressure in range (1e4, 1e7) and compute corresponding densities
p = np.logspace(4, 7, 100)
r = rho(p)
-# Plot density vs. pressure.
+# plot density vs. pressure
plt.semilogx(p, r)
plt.show()
diff --git a/exercises/solution/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh b/exercises/solution/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh
index a5a4a5f7c5e12cec758fc6e1d33ac31346fed64b..9c47f792166325f8eb9b9c4cb723923d362a4512 100644
--- a/exercises/solution/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh
+++ b/exercises/solution/exercise-fluidsystem/fluidsystems/h2omycompressiblecomponent.hh
@@ -1,3 +1,4 @@
+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
//
// SPDX-FileCopyrightInfo: Copyright © DuMux-Course contributors, see AUTHORS.md in root folder
@@ -7,8 +8,8 @@
* \file
*
* \brief A fluid system with water and a fictitious component, which is to be
- * implemented in exercise-fluidsystem 2.2, as phases and components. This
- * fluid system is to be implemented in exercise-fluidsystem 3.
+ * implemented in exercise-fluidsystem a, as phases and components. This
+ * fluid system is to be implemented in exercise-fluidsystem b.
*/
#ifndef DUMUX_H2O_MYCOMPRESSIBLECOMPONENT_FLUID_SYSTEM_HH
#define DUMUX_H2O_MYCOMPRESSIBLECOMPONENT_FLUID_SYSTEM_HH
@@ -18,10 +19,10 @@
#include <dumux/material/fluidsystems/base.hh>
-// The fictitious component that was created in exercise-fluidsystem 2.2
+// the fictitious component that was created in exercise 3a
#include <exercises/solution/exercise-fluidsystem/components/mycompressiblecomponent.hh>
-// The binary coefficients corresponding to this fluid system
+// the binary coefficients corresponding to this fluid system
#include <exercises/solution/exercise-fluidsystem/binarycoefficients/h2omycompressiblecomponent.hh>
namespace Dumux {
@@ -29,7 +30,7 @@ namespace FluidSystems {
/*!
* \brief A compositional fluid consisting of two liquid phases,
- * which are water and a fictitious component from exercise-fluidsystem 2.2.
+ * which are water and a fictitious component from exercise-fluidsystem a.
*/
template <class Scalar,
class H2OType = Dumux::Components::TabulatedComponent<Dumux::Components::H2O<Scalar> > >
@@ -46,12 +47,12 @@ public:
static constexpr int numPhases = 2;
static constexpr int numComponents = 2;
- static constexpr int phase0Idx = 0; // Index of the first phase
- static constexpr int phase1Idx = 1; // Index of the second phase
+ static constexpr int phase0Idx = 0; // index of the first phase
+ static constexpr int phase1Idx = 1; // index of the second phase
static constexpr int H2OIdx = 0;
static constexpr int NAPLIdx = 1;
- // Export component indices to indicate the main component
+ // export component indices to indicate the main component
// of the corresponding phase at atmospheric pressure 1 bar
// and room temperature 20°C:
static constexpr int comp0Idx = H2OIdx;
@@ -209,7 +210,7 @@ public:
if (phaseIdx == phase0Idx)
{
- // see: doctoral thesis of Steffen Ochs 2007
+ // See: doctoral thesis of Steffen Ochs 2007
// Steam injection into saturated porous media : process analysis including experimental and numerical investigations
// http://elib.uni-stuttgart.de/bitstream/11682/271/1/Diss_Ochs_OPUS.pdf
Scalar rholH2O = H2O::liquidDensity(fluidState.temperature(phaseIdx), fluidState.pressure(phaseIdx));
@@ -235,7 +236,7 @@ public:
* The molar density for the simple relation is defined by the
* mass density \f$\rho_\alpha\f$ and the molar mass of the main component
*
- * The molar density for the complex relation is defined by the
+ * The molar density for the complrex relation is defined by the
* mass density \f$\rho_\alpha\f$ and the mean molar mass \f$\overline M_\alpha\f$:
*
* \f[\rho_{mol,\alpha} = \frac{\rho_\alpha}{\overline M_\alpha} \;.\f]
@@ -393,10 +394,8 @@ public:
DUNE_THROW(Dune::NotImplemented, "FluidSystems::H2OMyCompressibleComponent::kelvinVaporPressure()");
}
- /*!
- * \brief Partial pressures in the gas phase, taken from saturation vapor pressures.
+ /* partial pressures in the gas phase, taken from saturation vapor pressures
*/
-
template <class FluidState>
static Scalar partialPressureGas(const FluidState &fluidState, int phaseIdx,
int compIdx)
@@ -412,10 +411,8 @@ public:
DUNE_THROW(Dune::InvalidStateException, "non-existent component index " << compIdx);
}
- /*!
- * \brief Inverse vapor pressures, taken from inverse saturation vapor pressures.
+ /* inverse vapor pressures, taken from inverse saturation vapor pressures
*/
-
template <class FluidState>
static Scalar inverseVaporPressureCurve(const FluidState &fluidState,
int phaseIdx,
diff --git a/exercises/solution/exercise-fluidsystem/main.cc b/exercises/solution/exercise-fluidsystem/main.cc
index 6b94b931491e4d95c47e4123b281605cf53d2a50..4ad7547d25c844195a1ebcb11610261da40aa118 100644
--- a/exercises/solution/exercise-fluidsystem/main.cc
+++ b/exercises/solution/exercise-fluidsystem/main.cc
@@ -6,7 +6,7 @@
//
/*!
* \file
- * \brief The main file for exercise-fluidsystem. This solves an instationary problem
+ * \brief The main file for exercise 3. This solves an instationary problem
* with an implicit time discretization.
*/
#include <config.h>
@@ -36,12 +36,12 @@
int main(int argc, char** argv)
{
using namespace Dumux;
-
- // TODO: dumux-course-task 3:
- // TYPETAG is set in CMakeLists.txt as compile time definition.
- // Alternatively, you could write `using TypeTag = Properties::TTag::ExerciseFluidsystemTwoP;`
+
+ // define the type tag for this problem
+ // TYPETAG is set in CMakeLists.txt as compile time definition
+ // alternatively you could write `using TypeTag = Properties::TTag::ExerciseFluidsystemTwoP;`
// then, for the 2p2c problem you have to change this line to `using TypeTag = Properties::TTag::ExerciseFluidsystemTwoPTwoC;`
- // and recompile the executable.
+ // and recompile the executable
using TypeTag = Properties::TTag::TYPETAG;
// initialize MPI+x, finalize is done automatically on exit
@@ -85,7 +85,6 @@ int main(int argc, char** argv)
const auto tEnd = getParam<Scalar>("TimeLoop.TEnd");
const auto maxDt = getParam<Scalar>("TimeLoop.MaxTimeStepSize");
auto dt = getParam<Scalar>("TimeLoop.DtInitial");
-
// initialize the vtk output module
using IOFields = GetPropType<TypeTag, Properties::IOFields>;
VtkOutputModule<GridVariables, SolutionVector> vtkWriter(*gridVariables, x, problem->name());
diff --git a/exercises/solution/exercise-fluidsystem/spatialparams.hh b/exercises/solution/exercise-fluidsystem/spatialparams.hh
index 2769e5171a8908e52ca0b59e17094376ef5cf701..a5a51007f79012ee5242bb03b24065a256022215 100644
--- a/exercises/solution/exercise-fluidsystem/spatialparams.hh
+++ b/exercises/solution/exercise-fluidsystem/spatialparams.hh
@@ -13,10 +13,10 @@
#ifndef DUMUX_EXERCISE_FLUIDSYSTEM_SPATIAL_PARAMS_HH
#define DUMUX_EXERCISE_FLUIDSYSTEM_SPATIAL_PARAMS_HH
-// Include parent spatial parameters
+// include parent spatialparameters
#include <dumux/porousmediumflow/fvspatialparamsmp.hh>
-// Include fluid-matrix interaction relationship
+//include fluid matrix interaction relationship
#include <dumux/material/fluidmatrixinteractions/2p/brookscorey.hh>
#include <dumux/material/fluidmatrixinteractions/2p/linearmaterial.hh>
@@ -44,13 +44,13 @@ class ExerciseFluidsystemSpatialParams
using PcKrSwCurve = FluidMatrix::BrooksCoreyDefault<Scalar>;
public:
- // Export permeability type
+ // export permeability type
using PermeabilityType = Dune::FieldMatrix<Scalar, dim, dim>;
/*!
- * \brief The constructor.
+ * \brief The constructor
*
- * \param gridGeometry The finite volume grid geometry.
+ * \param gridGeometry The finite volume grid geometry
*/
ExerciseFluidsystemSpatialParams(std::shared_ptr<const GridGeometry>& gridGeometry)
: ParentType(gridGeometry)
@@ -59,8 +59,8 @@ public:
, pcKrSwCurve_("SpatialParams")
, lensPcKrSwCurve_("SpatialParams.Lens")
{
- // set main diagonal entries of the permeability tensor to a value
- // setting to one value means: isotropic, homogeneous
+ //set main diagonal entries of the permeability tensor to a value
+ //setting to one value means: isotropic, homogeneous
for (int i = 0; i < dim; i++)
{
K_[i][i] = 1e-7;
@@ -71,7 +71,7 @@ public:
/*!
* \brief Define the intrinsic permeability \f$\mathrm{[m^2]}\f$.
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
PermeabilityType permeabilityAtPos(const GlobalPosition& globalPos) const
@@ -84,7 +84,7 @@ public:
/*!
* \brief Define the porosity \f$\mathrm{[-]}\f$.
*
- * \param globalPos The global position.
+ * \param globalPos The global position
*/
Scalar porosityAtPos(const GlobalPosition& globalPos) const
{
@@ -94,8 +94,8 @@ public:
}
/*!
- * \brief Returns the fluid-matrix interaction law at a given location.
- * \param globalPos The global coordinates for the given location.
+ * \brief Returns the fluid-matrix interaction law at a given location
+ * \param globalPos The global coordinates for the given location
*/
auto fluidMatrixInteractionAtPos(const GlobalPosition& globalPos) const
{
@@ -107,19 +107,18 @@ public:
/*!
* \brief Function for defining which phase is to be considered as the wetting phase.
*
- * \return The wetting phase index.
- * \param globalPos The position of the center of the element.
+ * \return the wetting phase index
+ * \param globalPos The position of the center of the element
*/
template<class FluidSystem>
int wettingPhaseAtPos(const GlobalPosition& globalPos) const
{
- // Our fluid system is H2OMyCompressibleComponent.
+ // Our fluid system is H2OMyCompressibleComponent
// We want to define water as the wetting phase in
- // the entire domain (see fluidsystem for the phase indices).
-
- // TODO: dumux-course-task 4:
- // For task 4 (change of wettability), uncomment
- // the two commented lines below:
+ // the entire domain (see fluid system for the phase indices)
+
+ // For the last task (change of wettability), uncomment
+ // the two commented lines below
//if (isInLens(globalPos))
// return FluidSystem::phase1Idx;
@@ -137,8 +136,8 @@ public:
}
/*!
- * \brief Returns the temperature at the domain at the given position.
- * \param globalPos The position in global coordinates where the temperature should be specified.
+ * \brief Returns the temperature at the domain at the given position
+ * \param globalPos The position in global coordinates where the temperature should be specified
*/
Scalar temperatureAtPos(const GlobalPosition& globalPos) const
{
@@ -151,7 +150,7 @@ private:
Dune::FieldMatrix<Scalar, dim, dim> K_;
Dune::FieldMatrix<Scalar, dim, dim> KLens_;
- // Object that holds the values/parameters of the selected fluid-matrix interaction relationship.
+ // Object that holds the values/parameters of the selected fluid matrix interaction relationship
const PcKrSwCurve pcKrSwCurve_;
const PcKrSwCurve lensPcKrSwCurve_;
};