THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.
The following task has a new comment added:
FS#172 - Restructure the way heat conduction is implemented in the 2p(2c)ni
boxmodel
User who did this - Klaus Mosthaf (Klaus)
----------
The combination of the heat conductivities of the solid material and the fluids
has to be done at some point. I also thought about putting it somewhere into
the model, as it is done for the diffusion coefficient. However, this averaging
(ldry + sqrt(Sl) * (ldry - lsat)) is just one option and it is possible, that
one might do the calculation of the effective heat conductivity differently. If
the calculation is done in the model, it would be necessary to change/overwrite
the methods there.
After spending some more thoughts and after having talked to Benjamin and
Bernd, we ended up with an even better solution:
The model to compute the effective thermal conductivity is put into a static
class \"Somerton\". It is then treated in the same way as other fluid-matrix
interactions and will be put into dumux/material/fluidmatrixinteractions/2p. A
property, probably called ThermalConductivityModel, will be set in the 2p(2c)ni
model and the default choice is the Somerton model. The property allows to
exchange the method to compute the effective thermal conductivity. The spatial
parameter file requires only the method heatConductivitySolid(...) instead of
*matrixHeatFlux(). A diff of the spatial parameters can be found below.
I would propose to handle the effective diffusion coefficient in the same way;
so far the computation is done within the model.
Index: test/implicit/2p2cni/waterairspatialparams.hh
===================================================================
--- test/implicit/2p2cni/waterairspatialparams.hh (revision 9992)
+++ test/implicit/2p2cni/waterairspatialparams.hh (working copy)
@@ -114,6 +114,9 @@
// porosities
finePorosity_ = 0.3;
coarsePorosity_ = 0.3;
+
+ // heat conductivity of granite
+ lambdaSolid_ = 2.8;
// residual saturations
fineMaterialParams_.setSwr(0.2);
@@ -209,6 +212,18 @@
}
/*!
+ * \\brief Returns the thermal conductivity \\f$[W/m^2]\\f$ of the porous
material.
+ *
+ * \\param pos The global position
+ */
+ Scalar thermalConductivitySolid(const Element &element,
+ const FVElementGeometry &fvGeometry,
+ const int scvIdx) const
+ {
+ return lambdaSolid_;
+ }
+
+ /*!
* \\brief Calculate the heat flux \\f$[W/m^2]\\f$ through the
* rock matrix based on the temperature gradient \\f$[K / m]\\f$
*
@@ -266,6 +281,9 @@
Scalar finePorosity_;
Scalar coarsePorosity_;
+ // heat conductivity of the solid material only
+ Scalar lambdaSolid_;
+
MaterialLawParams fineMaterialParams_;
MaterialLawParams coarseMaterialParams_;
};
----------
More information can be found at the following URL:
http://www.dumux.uni-stuttgart.de/flyspray/index.php?do=details&task_id=172#comment296
You are receiving this message because you have requested it from the Flyspray
bugtracking system. If you did not expect this message or don't want to
receive mails in future, you can change your notification settings at the URL
shown above.
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
