diff --git a/doc/doxygen/extradoc/parameterlist.txt b/doc/doxygen/extradoc/parameterlist.txt
index 16975c5942d89c54560ab7431061cde26d74b7fe..a54116eb3f488c2516572136f1a2ea6070463f08 100644
--- a/doc/doxygen/extradoc/parameterlist.txt
+++ b/doc/doxygen/extradoc/parameterlist.txt
@@ -60,6 +60,21 @@
  * | Component            | SolidThermalConductivity                      | Scalar                   | -               | Thermal conductivity of the component as a solid                                                                                                       |
  * | Component            | TriplePressure                                | Scalar                   | -               | The triple pressure of the component                                                                                                                   |
  * | Component            | TripleTemperature                             | Scalar                   | -               | The triple temperature of the component                                                                                                                |
+ * | \b DualNetwork       | Bi                                            | Scalar                   | 1.0             | The Biot number.                                                                                                                                       |
+ * | DualNetwork          | ConnectionAreaShapeFactor                     | Scalar                   | 0.9             | The factor for the connection are shape.                                                                                                               |
+ * | DualNetwork          | ConvectionLambaExponent                       | Scalar                   | 0.4             | The exponent for the convection lambda.                                                                                                                |
+ * | DualNetwork          | ConvectionLambaFactor                         | Scalar                   | 0.9             | The multiplicative factor for the convection lambda.                                                                                                   |
+ * | DualNetwork          | FixedConvectionLambda                         | Scalar                   | -1.0            | Fixed value for the convecion lambda if used.                                                                                                          |
+ * | DualNetwork          | FluidTemperatureMode                          | std::string              | mean            | The computation mode of the fluid temperature.                                                                                                         |
+ * | DualNetwork          | Nu                                            | Scalar                   | 1.0             | The Nusselt number.                                                                                                                                    |
+ * | DualNetwork          | OverlapFactor                                 | Scalar                   | -               | The overlap factor for the computation of exact pore radius.                                                                                           |
+ * | DualNetwork          | SourceVerboseForDof                           | int                      | -1              | Whether status should be printed.                                                                                                                      |
+ * | DualNetwork          | SphereRadius                                  | Scalar                   | 50e-6           | The radius of a spherical pore.                                                                                                                        |
+ * | DualNetwork          | UseAverageConvectionArea                      | bool                     | false           | Whether to use an averaged convection area.                                                                                                            |
+ * | DualNetwork          | UseExactConnectionAreaSphere                  | bool                     | false           | Whether to use the exact connection area.                                                                                                              |
+ * | DualNetwork          | UseExactConnectionLength                      | bool                     | false           | Whether to use the exact connection length.                                                                                                            |
+ * | DualNetwork          | UseExactPoreRadiusSolid                       | bool                     | false           | Whether to use the exact pore radius for the solid.                                                                                                    |
+ * | DualNetwork          | UseExactPoreRadiusVoid                        | bool                     | false           | Whether to use the exact pore radius for the void space.                                                                                               |
  * | \b ElectroChemistry  | ActivationBarrier                             | Scalar                   | -               | The activation barrier to calculate the exchange current density.                                                                                      |
  * | ElectroChemistry     | CellVoltage                                   | Scalar                   | -               | The voltage of the fuel cell.                                                                                                                          |
  * | ElectroChemistry     | MaxIterations                                 | int                      | -               | The maximum number of iterations in iteatively (Newton solver) calculating the current density.                                                        |
@@ -86,6 +101,7 @@
  * | \b FluxOverSurface   | Verbose                                       | bool                     | false           | For enabling or disabling the console output                                                                                                           |
  * | \b Forchheimer       | MaxIterations                                 | std::size_t              | 30              | The maximum number of Newton iterations for solving the Forchheimer equation                                                                           |
  * | Forchheimer          | NewtonTolerance                               | Scalar                   | 1e-12           | The error tolerance in the Newton method for solving the Forchheimer equation                                                                          |
+ * | \b FouriersLaw       | UseVolumeEqualPyramid                         | bool                     | false           | Whether to deduce pyramid base length from volume. The base side length for the frustum of the pyramid is obtained by assuming that the volume of the frustum is the same as half the volume of the pore body. Otherwise the side length is set as the diameter of the pore body. |
  * | \b FreeFlow          | EnableDilatationTerm                          | bool                     | false           | For enabling the turbulent dilation term.                                                                                                              |
  * | FreeFlow             | EnableUnsymmetrizedVelocityGradient           | bool                     | false           | For enabling unsymmetrized velocity gradient. If false consider the shear stress caused by the gradient of the velocities normal to our face of interest. |
  * | FreeFlow             | EnableUnsymmetrizedVelocityGradientForBeaversJoseph | bool                     | false           | For enabling unsymmetrized velocity gradient for the Beavers Joseph coupling condition. If true and if the current scvf is on a boundary and if a Dirichlet BC for the pressure or a BJ condition for the slip velocity is set there, assume a tangential velocity gradient of zero along the lateral face. |
@@ -94,6 +110,7 @@
  * | Grid                 | AllowIntersectingDiagonals                    | bool                     | true            | Whether to allow diagonals to intersect in the context of the generation of a structured-lattice pore-network.                                         |
  * | Grid                 | Angular0/1/2                                  | std::vector<Scalar>      | -               | min/max value for angular coordinate. Cake grids can be created by either specifying Radial,Angular or Axial in all coordinate directions.             |
  * | Grid                 | Axial0/1/2                                    | std::vector<Scalar>      | -               | min/max value for axial coordinate. Cake grids can be created by either specifying Radial,Angular or Axial in all coordinate directions.               |
+ * | Grid                 | BinaryMask                                    | std::string              | -               | Name of the mask file.                                                                                                                                 |
  * | Grid                 | BoundaryPoreLabels                            | std::vector<std::string> | -               | With this, the boundary faces can be set in the format xmin xmax ymin ymax (zmin zmax).                                                                |
  * | Grid                 | BoundarySegments                              | bool                     | false           | For the dune gmsh reader: Whether to insert boundary segments into the grid                                                                            |
  * | Grid                 | CapPoreRadii                                  | bool                     | true            | If true a maximal pore radius is set.                                                                                                                  |
@@ -119,11 +136,12 @@
  * | Grid                 | GmshPhysicalEntityThreshold                   | std::size_t              | 0               | Boundary element index threshold                                                                                                                       |
  * | Grid                 | Grading0/1/2                                  | std::vector<Scalar>      | -               | For a grid with zones, grading factors for the x-zones. 1.0 means all cells within this zone have equal extension in x-direction. Negative factors are possible. |
  * | Grid                 | Image                                         | std::string              | -               | The image file if the sub grid is constructed from a raster image                                                                                      |
+ * | Grid                 | KeepOnlyConnected                             | std::string              | -               | Whether non-connected cells should be removed.                                                                                                         |
  * | Grid                 | KeepPhysicalOverlap                           | bool                     | true            | Whether to keep the physical overlap in physical size or in number of cells upon refinement                                                            |
  * | Grid                 | LeftBoundary                                  | CoordinateType           | 0.0             | The start coordinate of a 1D grid                                                                                                                      |
  * | Grid                 | LowerLeft                                     | GlobalPosition           | GlobalPosition(0.0) | The lowerLeft corner of a structured grid                                                                                                              |
  * | Grid                 | MakeConsistentlyOriented                      | bool                     | true            | Whether to make sure that all element intersections follow the same local orientation and indexing no matter how the elements may be rotated or twisted. |
- * | Grid                 | Marker                                        | bool                     | false           | To customize the subgrid generation.                                                                                                                   |
+ * | Grid                 | Marker                                        | char                     | 0               | To customize the subgrid generation.                                                                                                                   |
  * | Grid                 | MaxPoreInscribedRadius                        | Scalar                   | -               | In the case of a uniform random distribution, this specifies the maximum pore radius.                                                                  |
  * | Grid                 | MeanPoreInscribedRadius                       | Scalar                   | -               | In the case of a lognormal random distribution, this specifies the mean pore radius.                                                                   |
  * | Grid                 | MinPoreInscribedRadius                        | Scalar                   | -               | In the case of a uniform random distribution, this specifies the minimum pore radius.                                                                  |
@@ -194,6 +212,7 @@
  * | KEpsilon             | YPlusThreshold                                | Scalar                   | 30.0            | yPlus below this value is considered as near-wall region                                                                                               |
  * | \b KOmega            | EnableDissipationLimiter                      | bool                     | true            | Whether to enable the dissipation limiter                                                                                                              |
  * | KOmega               | EnableProductionLimiter                       | bool                     | false           | Whether to enable the production limiter                                                                                                               |
+ * | \b LinearPDESolver   | EnableDynamicOutput                           | bool                     | true            | Whether dynamic output should be enabled.                                                                                                              |
  * | \b LinearSolver      | GMResRestart                                  | int                      | 10              | cycles before restarting                                                                                                                               |
  * | LinearSolver         | MaxIterations                                 | int                      | 250             | The maximum iterations of the linear solver                                                                                                            |
  * | LinearSolver         | MaxOrthogonalizationVectors                   | int                      | 10              | Maximal number of previous vectors which are orthogonalized against the new search direction                                                           |
@@ -230,7 +249,6 @@
  * | LinearSolver         | ResidualReduction                             | double                   | 1e-13(linear solver),1e-6(nonlinear) | The residual reduction threshold, i.e. stopping criterion                                                                                              |
  * | LinearSolver         | Restart                                       | int                      | 10              | cycles before restarting                                                                                                                               |
  * | LinearSolver         | Type                                          | std::string              | -               | The type of linear solver, e.g. restartedflexiblegmressolver or uzawa                                                                                  |
- * | LinearSolver         | UMFPackOrdering                               | int                      | 1               | You can choose from one of the following ordering strategies: 0: UMFPACK_ORDERING_CHOLMOD, 1: UMFPACK_ORDERING_AMD (default), 2: UMFPACK_ORDERING_GIVEN, 3: UMFPACK_ORDERING_METIS, 4: UMFPACK_ORDERING_BEST, 5: UMFPACK_ORDERING_NONE, 6: UMFPACK_ORDERING_USER. See https://fossies.org/linux/SuiteSparse/UMFPACK/Doc/UMFPACK_UserGuide.pdf page 17 for details. |
  * | LinearSolver         | Verbosity                                     | int                      | 0               | The verbosity level of the linear solver                                                                                                               |
  * | \b LoadSolution      | CellCenterPriVarNames                         | std::vector<std::string> | -               | Names of cell-centered primary variables of a model with staggered grid discretization                                                                 |
  * | LoadSolution         | FacePriVarNames                               | std::vector<std::string> | -               | Names of primary variables on the cell faces of a model with staggered grid discretization                                                             |
@@ -281,8 +299,10 @@
  * | \b Problem           | EnableGravity                                 | bool                     | -               | Whether to enable the gravity term                                                                                                                     |
  * | Problem              | EnableInertiaTerms                            | bool                     | -               | Whether to enable the inertia terms                                                                                                                    |
  * | Problem              | Name                                          | std::string              | -               | Set a name for a problem                                                                                                                               |
+ * | Problem              | PoissonRatio                                  | Scalar                   | -               | The poisson ratio.                                                                                                                                     |
  * | Problem              | SandGrainRoughness                            | Scalar                   | -               | The sand grain roughness                                                                                                                               |
  * | Problem              | UsePrimaryVariableSwitch                      | bool                     | -               | Whether to perform variable switch at a degree of freedom location                                                                                     |
+ * | Problem              | YoungsModulus                                 | Scalar                   | -               | The value of Young's modulus.                                                                                                                          |
  * | \b RANS              | EddyViscosityModel                            | std::string              | "vanDriest"     | Choose the eddy viscosity model                                                                                                                        |
  * | RANS                 | FlowDirectionAxis                             | int                      | 0               | The flow direction axis                                                                                                                                |
  * | RANS                 | IsFlatWallBounded                             | bool                     | false           | Set to true, if geometry consists of flat walls                                                                                                        |
diff --git a/doc/doxygen/extradoc/parameters.json b/doc/doxygen/extradoc/parameters.json
index b8038d4337617dc8d70b71013324d6a7f71ec0bc..36ea55025d5ad2fb9b819740679d78f1edd47931 100644
--- a/doc/doxygen/extradoc/parameters.json
+++ b/doc/doxygen/extradoc/parameters.json
@@ -541,6 +541,165 @@
             "The triple temperature of the component"
         ]
     },
+    "DualNetwork.Bi": {
+        "explanation": [
+            "The Biot number."
+        ],
+        "type": [
+            "Scalar"
+        ],
+        "defaultValue": [
+            "1.0"
+        ]
+    },
+    "DualNetwork.ConnectionAreaShapeFactor": {
+        "explanation": [
+            "The factor for the connection are shape."
+        ],
+        "type": [
+            "Scalar"
+        ],
+        "defaultValue": [
+            "0.9"
+        ]
+    },
+    "DualNetwork.ConvectionLambaExponent": {
+        "explanation": [
+            "The exponent for the convection lambda."
+        ],
+        "type": [
+            "Scalar"
+        ],
+        "defaultValue": [
+            "0.4"
+        ]
+    },
+    "DualNetwork.ConvectionLambaFactor": {
+        "explanation": [
+            "The multiplicative factor for the convection lambda."
+        ],
+        "type": [
+            "Scalar"
+        ],
+        "defaultValue": [
+            "0.9"
+        ]
+    },
+    "DualNetwork.FixedConvectionLambda": {
+        "explanation": [
+            "Fixed value for the convecion lambda if used."
+        ],
+        "type": [
+            "Scalar"
+        ],
+        "defaultValue": [
+            "-1.0"
+        ]
+    },
+    "DualNetwork.FluidTemperatureMode": {
+        "explanation": [
+            "The computation mode of the fluid temperature."
+        ],
+        "type": [
+            "std::string"
+        ],
+        "defaultValue": [
+            "mean"
+        ]
+    },
+    "DualNetwork.Nu": {
+        "explanation": [
+            "The Nusselt number."
+        ],
+        "type": [
+            "Scalar"
+        ],
+        "defaultValue": [
+            "1.0"
+        ]
+    },
+    "DualNetwork.OverlapFactor": {
+        "explanation": [
+            "The overlap factor for the computation of exact pore radius."
+        ]
+    },
+    "DualNetwork.SourceVerboseForDof": {
+        "explanation": [
+            "Whether status should be printed."
+        ],
+        "type": [
+            "int"
+        ],
+        "defaultValue": [
+            "-1"
+        ]
+    },
+    "DualNetwork.SphereRadius": {
+        "explanation": [
+            "The radius of a spherical pore."
+        ],
+        "type": [
+            "Scalar"
+        ],
+        "defaultValue": [
+            "50e-6"
+        ]
+    },
+    "DualNetwork.UseAverageConvectionArea": {
+        "explanation": [
+            "Whether to use an averaged convection area."
+        ],
+        "type": [
+            "bool"
+        ],
+        "defaultValue": [
+            "false"
+        ]
+    },
+    "DualNetwork.UseExactConnectionAreaSphere": {
+        "explanation": [
+            "Whether to use the exact connection area."
+        ],
+        "type": [
+            "bool"
+        ],
+        "defaultValue": [
+            "false"
+        ]
+    },
+    "DualNetwork.UseExactConnectionLength": {
+        "explanation": [
+            "Whether to use the exact connection length."
+        ],
+        "type": [
+            "bool"
+        ],
+        "defaultValue": [
+            "false"
+        ]
+    },
+    "DualNetwork.UseExactPoreRadiusSolid": {
+        "explanation": [
+            "Whether to use the exact pore radius for the solid."
+        ],
+        "type": [
+            "bool"
+        ],
+        "defaultValue": [
+            "false"
+        ]
+    },
+    "DualNetwork.UseExactPoreRadiusVoid": {
+        "explanation": [
+            "Whether to use the exact pore radius for the void space."
+        ],
+        "type": [
+            "bool"
+        ],
+        "defaultValue": [
+            "false"
+        ]
+    },
     "ElectroChemistry.ActivationBarrier": {
         "explanation": [
             "The activation barrier to calculate the exchange current density."
@@ -671,6 +830,17 @@
             "The error tolerance in the Newton method for solving the Forchheimer equation"
         ]
     },
+    "FouriersLaw.UseVolumeEqualPyramid":{
+        "explanation": [
+            "Whether to deduce pyramid base length from volume. The base side length for the frustum of the pyramid is obtained by assuming that the volume of the frustum is the same as half the volume of the pore body. Otherwise the side length is set as the diameter of the pore body."
+        ],
+        "type": [
+            "bool"
+        ],
+        "defaultValue": [
+            "false"
+        ]
+    },
     "FreeFlow.EnableDilatationTerm": {
         "explanation": [
             "For enabling the turbulent dilation term."
@@ -725,6 +895,15 @@
             "std::vector<Scalar>"
         ]
     },
+    "Grid.BinaryMask": {
+        "explanation": [
+            "Name of the mask file."
+        ],
+        "type": [
+            "std::string"
+        ]
+    },
+
     "Grid.BoundaryPoreLabels": {
         "explanation": [
             "With this, the boundary faces can be set in the format xmin xmax ymin ymax (zmin zmax)."
@@ -908,6 +1087,11 @@
             "The image file if the sub grid is constructed from a raster image"
         ]
     },
+    "Grid.KeepOnlyConnected": {
+        "explanation": [
+            "Whether non-connected cells should be removed."
+        ]
+    },
     "Grid.KeepPhysicalOverlap": {
         "explanation": [
             "Whether to keep the physical overlap in physical size or in number of cells upon refinement"
@@ -1535,6 +1719,11 @@
             "Whether to enable the production limiter"
         ]
     },
+    "LinearPDESolver.EnableDynamicOutput": {
+        "explanation": [
+            "Whether dynamic output should be enabled."
+        ]
+    },
     "LinearSolver.DirectSolverForVelocity": {
         "mode": "ignore"
     },
@@ -1976,11 +2165,6 @@
             "std::string"
         ]
     },
-    "LinearSolver.UMFPackOrdering": {
-        "explanation": [
-            "You can choose from one of the following ordering strategies: 0: UMFPACK_ORDERING_CHOLMOD, 1: UMFPACK_ORDERING_AMD (default), 2: UMFPACK_ORDERING_GIVEN, 3: UMFPACK_ORDERING_METIS, 4: UMFPACK_ORDERING_BEST, 5: UMFPACK_ORDERING_NONE, 6: UMFPACK_ORDERING_USER. See https://fossies.org/linux/SuiteSparse/UMFPACK/Doc/UMFPACK_UserGuide.pdf page 17 for details."
-        ]
-    },
     "LinearSolver.Verbosity": {
         "explanation": [
             "The verbosity level of the linear solver"
@@ -2254,6 +2438,11 @@
             "Set a name for a problem"
         ]
     },
+    "Problem.PoissonRatio": {
+        "explanation": [
+            "The poisson ratio."
+        ]
+    },
     "Problem.SandGrainRoughness": {
         "explanation": [
             "The sand grain roughness"
@@ -2264,6 +2453,14 @@
             "Whether to perform variable switch at a degree of freedom location"
         ]
     },
+    "Problem.YoungsModulus": {
+        "explanation": [
+            "The value of Young's modulus."
+        ],
+        "type": [
+            "Scalar"
+        ]
+    },
     "RANS.EddyViscosityModel": {
         "explanation": [
             "Choose the eddy viscosity model"