[cig-commits] commit 2054 by buerg to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 2054

Remove pressure and strain rate from material models.

U   trunk/aspire/include/aspect/material_model/interface.h
U   trunk/aspire/include/aspect/simulator.h
U   trunk/aspire/source/material_model/interface.cc
U   trunk/aspire/source/mesh_refinement/density.cc
U   trunk/aspire/source/mesh_refinement/thermal_energy_density.cc
U   trunk/aspire/source/postprocess/heat_flux_statistics.cc
U   trunk/aspire/source/postprocess/visualization/compressibility.cc
U   trunk/aspire/source/postprocess/visualization/density.cc
U   trunk/aspire/source/postprocess/visualization/diffusivity.cc
U   trunk/aspire/source/postprocess/visualization/specific_heat.cc
U   trunk/aspire/source/postprocess/visualization/thermal_conductivity.cc
U   trunk/aspire/source/postprocess/visualization/thermal_expansion_coefficient.cc
U   trunk/aspire/source/postprocess/visualization/viscosity.cc
U   trunk/aspire/source/simulator/assembly.cc


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2054&peg=2054

Diff:
Modified: trunk/aspire/include/aspect/material_model/interface.h
===================================================================
--- trunk/aspire/include/aspect/material_model/interface.h	2013-11-25 14:36:21 UTC (rev 2053)
+++ trunk/aspire/include/aspect/material_model/interface.h	2013-11-25 17:55:08 UTC (rev 2054)
@@ -498,29 +498,11 @@
             * Temperature values at the points given in the #position vector.
             */
            std::vector<double> temperature;
-           /**
-            * Pressure values at the points given in the #position vector.
-            */
-           std::vector<double> pressure;
             /**
              * Values of the compositional fields at the points given in the #position vector:
              composition[i][c] is the compositional field c at point i.
              */
             std::vector<std::vector<double> > composition;
-            /**
-             * Strain rate at the points given in the #position vector. Only
-             * the viscosity may depend on these values. This std::vector
-             * can be set to size 0 if the viscosity is not needed.
-             *
-             * @note The strain rate
-             * is computed as $arepsilon(\mathbf u)=rac 12 (
abla \mathbf u +
-             * 
abla \mathbf u^T)$, regardless of whether the model is
-             * compressible or not. This is relevant since in some other contexts,
-             * the strain rate in the compressible case is computed as
-             * $arepsilon(\mathbf u)=rac 12 (
abla \mathbf u +
-             * 
abla \mathbf u^T) - rac 13 
abla