[cig-commits] r1349 - trunk/aspect/include/aspect/material_model
heien at dealii.org
heien at dealii.org
Tue Nov 6 18:37:48 PST 2012
Author: heien
Date: 2012-11-06 19:37:48 -0700 (Tue, 06 Nov 2012)
New Revision: 1349
Modified:
trunk/aspect/include/aspect/material_model/interface.h
Log:
Added default calculation of thermal diffusivity to material model interface
Modified: trunk/aspect/include/aspect/material_model/interface.h
===================================================================
--- trunk/aspect/include/aspect/material_model/interface.h 2012-11-07 00:44:54 UTC (rev 1348)
+++ trunk/aspect/include/aspect/material_model/interface.h 2012-11-07 02:37:48 UTC (rev 1349)
@@ -182,7 +182,27 @@
const double pressure,
const std::vector<double> &compositional_fields,
const Point<dim> &position) const = 0;
+
/**
+ * Return the thermal diffusivity $\kappa$ of the model as a function of temperature,
+ * pressure and position. This is by default calculated as $\kappa=k/(\rho c_p)$ where
+ * $k$, $\rho$, and $c_p$ are determined by the appropriate material model functions.
+ * $\kappa$ has units $\textrm{m}^2/\textrm{s}$.
+ */
+ virtual double thermal_diffusivity (const double temperature,
+ const double pressure,
+ const std::vector<double> &compositional_fields,
+ const Point<dim> &position) const
+ {
+ double k, rho, c_p;
+
+ k = thermal_conductivity(temperature, pressure, compositional_fields, position);
+ rho = density(temperature, pressure, compositional_fields, position);
+ c_p = specific_heat(temperature, pressure, compositional_fields, position);
+
+ return k/(rho*c_p);
+ };
+ /**
* @}
*/
More information about the CIG-COMMITS
mailing list