From fcfc828dc9f75ef5f213035b64b8c5fe8c838215 Mon Sep 17 00:00:00 2001
From: Wenqing Wang <wenqing.wang@ufz.de>
Date: Fri, 16 Jun 2023 16:49:41 +0200
Subject: [PATCH] [ProcessLib] Update the usage of the vapour diffusion model
---
.../ThermoRichardsFlow/ThermoRichardsFlowFEM-impl.h | 2 ++
.../ConstitutiveCommon/TRMVaporDiffusion.cpp | 13 +++++++------
2 files changed, 9 insertions(+), 6 deletions(-)
diff --git a/ProcessLib/ThermoRichardsFlow/ThermoRichardsFlowFEM-impl.h b/ProcessLib/ThermoRichardsFlow/ThermoRichardsFlowFEM-impl.h
index c3dd6ccb814..f3805644612 100644
--- a/ProcessLib/ThermoRichardsFlow/ThermoRichardsFlowFEM-impl.h
+++ b/ProcessLib/ThermoRichardsFlow/ThermoRichardsFlowFEM-impl.h
@@ -534,6 +534,7 @@ void ThermoRichardsFlowLocalAssembler<ShapeFunction, GlobalDim>::
variables.porosity = phi;
double const D_v =
+ phi * (1.0 - S_L) *
gas_phase->property(MPL::PropertyType::diffusion)
.template value<double>(variables, x_position, t, dt);
@@ -990,6 +991,7 @@ void ThermoRichardsFlowLocalAssembler<ShapeFunction, GlobalDim>::assemble(
variables.porosity = phi;
double const D_v =
+ phi * (1.0 - S_L) *
gas_phase->property(MPL::PropertyType::diffusion)
.template value<double>(variables, x_position, t, dt);
diff --git a/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/TRMVaporDiffusion.cpp b/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/TRMVaporDiffusion.cpp
index 2435e1ea94e..a7c17a5430d 100644
--- a/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/TRMVaporDiffusion.cpp
+++ b/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/TRMVaporDiffusion.cpp
@@ -74,7 +74,10 @@ void TRMVaporDiffusionModel<DisplacementDim>::eval(
double const phi = poro_data.phi;
variables.porosity = phi;
+
+ double const S_g = 1.0 - S_L_data.S_L;
double const D_v =
+ phi * S_g *
gas_phase->property(MPL::PropertyType::diffusion)
.template value<double>(variables, x_t.x, x_t.t, x_t.dt);
@@ -89,21 +92,19 @@ void TRMVaporDiffusionModel<DisplacementDim>::eval(
MaterialPropertyLib::PropertyType::specific_heat_capacity)
.template value<double>(variables, x_t.x, x_t.t, x_t.dt);
- out.M_TT_X_NTN +=
- out.heat_capacity_vapor * rho_wv * (1 - S_L_data.S_L) * phi;
+ out.M_TT_X_NTN += out.heat_capacity_vapor * rho_wv * S_g * phi;
out.storage_coefficient_by_water_vapor =
- phi *
- (rho_wv * dS_L_data.dS_L_dp_cap + (1 - S_L_data.S_L) * drho_wv_dp);
+ phi * (rho_wv * dS_L_data.dS_L_dp_cap + S_g * drho_wv_dp);
- out.M_pT_X_NTN += phi * (1 - S_L_data.S_L) * drho_wv_dT;
+ out.M_pT_X_NTN += phi * S_g * drho_wv_dT;
//
// Latent heat term
//
if (gas_phase->hasProperty(MPL::PropertyType::specific_latent_heat))
{
- double const factor = phi * (1 - S_L_data.S_L) / rho_L_data.rho_LR;
+ double const factor = phi * S_g / rho_L_data.rho_LR;
// The volumetric latent heat of vaporization of liquid water
double const L0 =
gas_phase->property(MPL::PropertyType::specific_latent_heat)
--
GitLab