From e0a6142369f19e3bfbd8a1f9607e5d2cb3af12fc Mon Sep 17 00:00:00 2001
From: Thomas Fischer <thomas.fischer@ufz.de>
Date: Tue, 6 Jun 2017 11:37:53 +0200
Subject: [PATCH] [PL] RichardsComponentTransport: Change Bp.
Change the assembly for the pressure related entries of the local right hand side.
---
.../RichardsComponentTransport/RichardsComponentTransportFEM.h | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM.h b/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM.h
index 1fd44aaca6d..f4b1cab4453 100644
--- a/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM.h
+++ b/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM.h
@@ -306,7 +306,7 @@ public:
ip_data.mass_operator;
if (_process_data.has_gravity)
- Bp += w * density * dNdx.transpose() * K_over_mu * b;
+ Bp += w * density * dNdx.transpose() * K * b * (k_rel/mu);
/* with Oberbeck-Boussing assumption density difference only exists
* in buoyancy effects */
}
--
GitLab