From 07b9f1f0ef98ddb22a96e618922f91714261e91d Mon Sep 17 00:00:00 2001
From: HBShaoUFZ <haibing.shao@ufz.de>
Date: Mon, 12 Mar 2018 18:00:10 +0100
Subject: [PATCH] This safety check will fix the nan bug.
---
ProcessLib/HT/HTFEM.h | 23 +++++++++++++++--------
1 file changed, 15 insertions(+), 8 deletions(-)
diff --git a/ProcessLib/HT/HTFEM.h b/ProcessLib/HT/HTFEM.h
index 6d1972a3c11..27e67395b3e 100644
--- a/ProcessLib/HT/HTFEM.h
+++ b/ProcessLib/HT/HTFEM.h
@@ -135,14 +135,21 @@ protected:
_material_properties.thermal_dispersivity_transversal(t, pos)[0];
double const velocity_magnitude = velocity.norm();
- GlobalDimMatrixType const thermal_dispersivity =
- fluid_density * specific_heat_capacity_fluid *
- (thermal_dispersivity_transversal * velocity_magnitude * I +
- (thermal_dispersivity_longitudinal -
- thermal_dispersivity_transversal) /
- velocity_magnitude * velocity * velocity.transpose());
-
- return thermal_conductivity * I + thermal_dispersivity;
+
+ if (velocity_magnitude < std::numeric_limits<double>::epsilon())
+ {
+ return thermal_conductivity * I;
+ }
+ else
+ {
+ GlobalDimMatrixType const thermal_dispersivity =
+ fluid_density * specific_heat_capacity_fluid *
+ (thermal_dispersivity_transversal * velocity_magnitude * I +
+ (thermal_dispersivity_longitudinal -
+ thermal_dispersivity_transversal) /
+ velocity_magnitude * velocity * velocity.transpose());
+ return thermal_conductivity * I + thermal_dispersivity;
+ }
}
std::vector<double> const& getIntPtDarcyVelocityLocal(
--
GitLab