void
PorousFlowDensityConstBulk::computeQpProperties()
{
  /// Density and derivatives wrt pressure at the nodes
  _density_nodal[_qp] = density(_porepressure_nodal[_qp][_phase_num]);
  _ddensity_nodal_dp[_qp] = dDensity_dP(_porepressure_nodal[_qp][_phase_num]);

  /// Density and derivatives wrt pressure at the qps
  _density_qp[_qp] = density(_porepressure_qp[_qp][_phase_num]);
  _ddensity_qp_dp[_qp] = dDensity_dP(_porepressure_qp[_qp][_phase_num]);
}
Beispiel #2
0
void
PorousFlowIdealGas::computeQpProperties()
{
    /// Density and derivatives wrt pressure and temperature at the nodes
    _density_nodal[_qp] = density(_porepressure_nodal[_qp][_phase_num], _temperature_nodal[_qp][_phase_num], _molar_mass);
    _ddensity_nodal_dp[_qp] = dDensity_dP(_temperature_nodal[_qp][_phase_num], _molar_mass);
    _ddensity_nodal_dt[_qp] = dDensity_dT(_porepressure_nodal[_qp][_phase_num], _temperature_nodal[_qp][_phase_num], _molar_mass);

    /// Density and derivatives wrt pressure and temperature at the qps
    _density_qp[_qp] = density(_porepressure_qp[_qp][_phase_num], _temperature_qp[_qp][_phase_num], _molar_mass);
    _ddensity_qp_dp[_qp] = dDensity_dP(_temperature_qp[_qp][_phase_num], _molar_mass);
    _ddensity_qp_dt[_qp] = dDensity_dT(_porepressure_qp[_qp][_phase_num], _temperature_qp[_qp][_phase_num], _molar_mass);
}
void
PorousFlowDensityConstBulk::computeQpProperties()
{
  _density[_qp] = density(_porepressure[_qp][_phase_num]);
  _ddensity_dp[_qp] = dDensity_dP(_porepressure[_qp][_phase_num]);
  _ddensity_dt[_qp] = 0.0;
}
Beispiel #4
0
void
PorousFlowMethane::computeQpProperties()
{
  /// Density and derivatives wrt pressure and temperature at the nodes
  _density_nodal[_qp] = density(_porepressure_nodal[_qp][_phase_num], _temperature_nodal[_qp][_phase_num]);
  _ddensity_nodal_dp[_qp] = dDensity_dP(_temperature_nodal[_qp][_phase_num]);
  _ddensity_nodal_dt[_qp] = dDensity_dT(_porepressure_nodal[_qp][_phase_num], _temperature_nodal[_qp][_phase_num]);

  /// Density and derivatives wrt pressure and temperature at the qps
  _density_qp[_qp] = density(_porepressure_qp[_qp][_phase_num], _temperature_qp[_qp][_phase_num]);
  _ddensity_qp_dp[_qp] = dDensity_dP(_temperature_qp[_qp][_phase_num]);
  _ddensity_qp_dt[_qp] = dDensity_dT(_porepressure_qp[_qp][_phase_num], _temperature_qp[_qp][_phase_num]);

  /// Viscosity and derivative wrt temperature at the nodes
  _viscosity_nodal[_qp] = viscosity(_temperature_nodal[_qp][_phase_num]);
  _dviscosity_nodal_dt[_qp] = dViscosity_dT(_temperature_nodal[_qp][_phase_num]);
}