Correct issues with concrete moisture transport PDE variables closes i…

…daholab#224

This does two things:
1) Makes the time dependent term be in terms of RH rather than W
2) Removes the dehydration kernel because it's not clear how to cast
   that in terms of RH rather than W. We will add this back in later
   when it is needed, and make sure that is done correctly.

doc/content/source/materials/ConcreteThermalMoisture.md

@@ -350,7 +350,7 @@ A comprehensive set of constitutive models and parameters for moisture diffusion
 The governing equation for moisture diffusion in concrete is formulated by using relative humidity, $H$, as the primary variable:
 
 !equation id=moisture_governing
-\frac{\partial{W}}{\partial{H}} \frac{\partial{H}}{\partial{t}} = \nabla (D_h\nabla H) + \nabla (D_{ht}\nabla T) + \frac{\partial{W_d}}{\partial{t}}
+\frac{\partial{H}}{\partial{t}} = \nabla (D_h\nabla H) + \nabla (D_{ht}\nabla T)
 
 where
 
@@ -360,10 +360,9 @@ $P_{vs}$ = saturate vapor pressure $= P_{atm}\exp\left(4871.3\frac{T-100}{373.15
 $P_{atm}$ = standard atmospheric pressure $= 101.325 Pa$ \\
 $D_h$ = moisture diffusivity (also referred as humidity diffusivity), $m^2/s$\\
 $D_{ht}$= coupled moisture diffusivity under the influence of a temperature gradient, $m^2/s$\\
-$W_d$= total mass of free evaporable water released into the pores by dehydration of the cement paste\\
 $t$  = time, $s$
 
-The term on the left side of [!eqref](moisture_governing) represents time-dependent effects, and is provided by [ConcreteMoistureTimeIntegration](ConcreteMoistureTimeIntegration.md). The first term on the right side of [!eqref](moisture_governing) represents Fickian diffusion, and the second term represents Soret diffusion. These are both provided by [ConcreteMoistureDiffusion](ConcreteMoistureDiffusion.md). The third term on the right hand side of this equation represents a source due to dehydrated water, and is provided by [ConcreteMoistureDehydration](ConcreteMoistureDehydration.md).
+The term on the left side of [!eqref](moisture_governing) represents time-dependent effects, and is provided by [ConcreteMoistureTimeIntegration](ConcreteMoistureTimeIntegration.md). The first term on the right side of [!eqref](moisture_governing) represents Fickian diffusion, and the second term represents Soret diffusion. These are both provided by [ConcreteMoistureDiffusion](ConcreteMoistureDiffusion.md).
 
 Moisture diffusivity $D_h$ depends on the relative humidity, $H$. Thus the moisture diffusion governing equation is highly nonlinear. The following sections describes in detail the constitutive models for moisture diffusivity.
 

include/kernels/ConcreteMoistureDiffusion.h

@@ -60,4 +60,6 @@ class ConcreteMoistureDiffusion : public Diffusion
 
   /// Temperature gradient
   const VariableGradient & _grad_T;
+
+  const MaterialProperty<Real> & _moisture_capacity; //  dW/dH
 };

src/kernels/ConcreteMoistureDiffusion.C

@@ -31,19 +31,22 @@ ConcreteMoistureDiffusion::ConcreteMoistureDiffusion(const InputParameters & par
   : Diffusion(parameters),
     _Dh(getMaterialProperty<Real>("humidity_diffusivity")),
     _Dht(getMaterialProperty<Real>("humidity_diffusivity_thermal")),
-    _grad_T(coupledGradient("temperature"))
+    _grad_T(coupledGradient("temperature")),
+    _moisture_capacity(getMaterialProperty<Real>("moisture_capacity"))
 {
 }
 
 Real
 ConcreteMoistureDiffusion::computeQpResidual()
 {
-  return _Dh[_qp] * Diffusion::computeQpResidual() + _Dht[_qp] * _grad_test[_i][_qp] * _grad_T[_qp];
+  //return _moisture_capacity[_qp] * (_Dh[_qp] * Diffusion::computeQpResidual() + _Dht[_qp] * _grad_test[_i][_qp] * _grad_T[_qp]);
+  return (_Dh[_qp] * Diffusion::computeQpResidual() + _Dht[_qp] * _grad_test[_i][_qp] * _grad_T[_qp]);
 }
 
 Real
 ConcreteMoistureDiffusion::computeQpJacobian()
 {
+  //return _moisture_capacity[_qp] * _Dh[_qp] * Diffusion::computeQpJacobian();
   return _Dh[_qp] * Diffusion::computeQpJacobian();
 }
 

src/kernels/ConcreteMoistureTimeIntegration.C

@@ -34,13 +34,13 @@ Real
 ConcreteMoistureTimeIntegration::computeQpResidual()
 {
   // self accumulation term
-  return _moisture_capacity[_qp] * TimeDerivative::computeQpResidual();
+  return TimeDerivative::computeQpResidual();
 }
 
 Real
 ConcreteMoistureTimeIntegration::computeQpJacobian()
 {
-  return _moisture_capacity[_qp] * TimeDerivative::computeQpJacobian();
+  return TimeDerivative::computeQpJacobian();
 }
 
 Real