From 8d1348c42391a81195d94525db080ffab54b9e5c Mon Sep 17 00:00:00 2001
From: Wenqing Wang <wenqing.wang@ufz.de>
Date: Thu, 18 Jan 2024 13:29:47 +0100
Subject: [PATCH] [LIE/HM] Removed the old output of nodal aperture and jump
---
.../HydroMechanics/HydroMechanicsProcess.cpp | 54 -------------------
.../HydroMechanicsProcessData.h | 2 -
2 files changed, 56 deletions(-)
diff --git a/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcess.cpp b/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcess.cpp
index 19046161d0e..c1ecac885c1 100644
--- a/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcess.cpp
+++ b/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcess.cpp
@@ -340,18 +340,6 @@ void HydroMechanicsProcess<GlobalDim>::initializeConcreteProcess(
_process_data.mesh_prop_fracture_shear_failure =
mesh_prop_fracture_shear_failure;
- auto mesh_prop_nodal_w = MeshLib::getOrCreateMeshProperty<double>(
- const_cast<MeshLib::Mesh&>(mesh), "nodal_w",
- MeshLib::MeshItemType::Node, GlobalDim);
- mesh_prop_nodal_w->resize(mesh.getNumberOfNodes() * GlobalDim);
- _process_data.mesh_prop_nodal_w = mesh_prop_nodal_w;
-
- auto mesh_prop_nodal_b = MeshLib::getOrCreateMeshProperty<double>(
- const_cast<MeshLib::Mesh&>(mesh), "nodal_aperture",
- MeshLib::MeshItemType::Node, 1);
- mesh_prop_nodal_b->resize(mesh.getNumberOfNodes());
- _process_data.mesh_prop_nodal_b = mesh_prop_nodal_b;
-
auto mesh_prop_nodal_p = MeshLib::getOrCreateMeshProperty<double>(
const_cast<MeshLib::Mesh&>(mesh), "pressure_interpolated",
MeshLib::MeshItemType::Node, 1);
@@ -454,48 +442,6 @@ void HydroMechanicsProcess<GlobalDim>::postTimestepConcreteProcess(
(*x[process_id])[global_index];
}
}
-
- // compute nodal w and aperture
- auto const& R = _process_data.fracture_property->R;
- auto const& b0 = _process_data.fracture_property->aperture0;
- MeshLib::PropertyVector<double>& vec_w = *_process_data.mesh_prop_nodal_w;
- MeshLib::PropertyVector<double>& vec_b = *_process_data.mesh_prop_nodal_b;
-
- auto compute_nodal_aperture =
- [&](std::size_t const node_id, double const w_n)
- {
- // skip aperture computation for element-wise defined b0 because there
- // are jumps on the nodes between the element's values.
- if (dynamic_cast<ParameterLib::MeshElementParameter<double> const*>(
- &b0))
- {
- return std::numeric_limits<double>::quiet_NaN();
- }
-
- ParameterLib::SpatialPosition x;
- x.setNodeID(node_id);
- return w_n + b0(/*time independent*/ 0, x)[0];
- };
-
- Eigen::VectorXd g(GlobalDim);
- Eigen::VectorXd w(GlobalDim);
- for (MeshLib::Node const* node : _vec_fracture_nodes)
- {
- auto const node_id = node->getID();
- g.setZero();
- for (int k = 0; k < GlobalDim; k++)
- {
- g[k] = mesh_prop_g[node_id * GlobalDim + k];
- }
-
- w.noalias() = R * g;
- for (int k = 0; k < GlobalDim; k++)
- {
- vec_w[node_id * GlobalDim + k] = w[k];
- }
-
- vec_b[node_id] = compute_nodal_aperture(node_id, w[GlobalDim - 1]);
- }
}
template <int GlobalDim>
diff --git a/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcessData.h b/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcessData.h
index ec6d974b10e..14055231579 100644
--- a/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcessData.h
+++ b/ProcessLib/LIE/HydroMechanics/HydroMechanicsProcessData.h
@@ -119,8 +119,6 @@ struct HydroMechanicsProcessData
MeshLib::PropertyVector<double>* mesh_prop_b = nullptr;
MeshLib::PropertyVector<double>* mesh_prop_k_f = nullptr;
MeshLib::PropertyVector<double>* mesh_prop_fracture_shear_failure = nullptr;
- MeshLib::PropertyVector<double>* mesh_prop_nodal_w = nullptr;
- MeshLib::PropertyVector<double>* mesh_prop_nodal_b = nullptr;
MeshLib::PropertyVector<double>* mesh_prop_nodal_p = nullptr;
MeshLib::PropertyVector<double>* mesh_prop_nodal_forces = nullptr;
--
GitLab