From 13cc8508d1a4eeb390a92d40cbfe321b1a81c51e Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <github@naumov.de>
Date: Tue, 28 Jul 2020 18:08:03 +0200
Subject: [PATCH] [PL/LF] Add HydraulicFlow output.
---
ProcessLib/LiquidFlow/LiquidFlowProcess.cpp | 9 +++++++++
ProcessLib/LiquidFlow/LiquidFlowProcess.h | 1 +
2 files changed, 10 insertions(+)
diff --git a/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp b/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
index f8438778edf..cb4fcb4545b 100644
--- a/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
+++ b/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
@@ -18,6 +18,7 @@
#include "MeshLib/PropertyVector.h"
// TODO(TF) used for output of flux, if output classes are ready this has to be changed
#include "MeshLib/IO/writeMeshToFile.h"
+#include "ProcessLib/Utils/ComputeResiduum.h"
#include "ProcessLib/Utils/CreateLocalAssemblers.h"
namespace ProcessLib
@@ -42,6 +43,9 @@ LiquidFlowProcess::LiquidFlowProcess(
_surfaceflux(std::move(surfaceflux))
{
DBUG("Create Liquid flow process.");
+
+ _hydraulic_flow = MeshLib::getOrCreateMeshProperty<double>(
+ mesh, "HydraulicFlow", MeshLib::MeshItemType::Node, 1);
}
void LiquidFlowProcess::initializeConcreteProcess(
@@ -80,6 +84,11 @@ void LiquidFlowProcess::assembleConcreteProcess(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
pv.getActiveElementIDs(), dof_table, t, dt, x, xdot, process_id, M, K,
b, _coupled_solutions);
+
+ auto const residuum = computeResiduum(*x[0], *xdot[0], M, K, b);
+ transformVariableFromGlobalVector(residuum, 0 /*variable id*/,
+ *_local_to_global_index_map,
+ *_hydraulic_flow, std::negate<double>());
}
void LiquidFlowProcess::assembleWithJacobianConcreteProcess(
diff --git a/ProcessLib/LiquidFlow/LiquidFlowProcess.h b/ProcessLib/LiquidFlow/LiquidFlowProcess.h
index 9c74180f54a..09bb2472040 100644
--- a/ProcessLib/LiquidFlow/LiquidFlowProcess.h
+++ b/ProcessLib/LiquidFlow/LiquidFlowProcess.h
@@ -108,6 +108,7 @@ private:
_local_assemblers;
std::unique_ptr<ProcessLib::SurfaceFluxData> _surfaceflux;
+ MeshLib::PropertyVector<double>* _hydraulic_flow = nullptr;
};
} // namespace LiquidFlow
--
GitLab