From b349a89d03c9c60cd6bcb24a8a4e06d4bc9144ab Mon Sep 17 00:00:00 2001
From: Thomas Fischer <thomas.fischer@ufz.de>
Date: Fri, 17 Jan 2020 12:32:05 +0100
Subject: [PATCH] [PL/ComponentTransport] Check if all
phases/components/properties exists.
---
.../ComponentTransportProcess.cpp | 57 +++++++++++++++++++
1 file changed, 57 insertions(+)
diff --git a/ProcessLib/ComponentTransport/ComponentTransportProcess.cpp b/ProcessLib/ComponentTransport/ComponentTransportProcess.cpp
index d1f4734ce81..38b9b68cad1 100644
--- a/ProcessLib/ComponentTransport/ComponentTransportProcess.cpp
+++ b/ProcessLib/ComponentTransport/ComponentTransportProcess.cpp
@@ -20,6 +20,61 @@ namespace ProcessLib
{
namespace ComponentTransport
{
+
+void checkMPLProperties(MeshLib::Mesh const& mesh,
+ ComponentTransportProcessData const& process_data)
+{
+ DBUG("Check the media properties of ComponentTransport process ...");
+
+ std::array<MaterialPropertyLib::PropertyType, 4> const
+ required_properties_medium = {
+ MaterialPropertyLib::PropertyType::porosity,
+ MaterialPropertyLib::PropertyType::transversal_dispersivity,
+ MaterialPropertyLib::PropertyType::longitudinal_dispersivity,
+ MaterialPropertyLib::PropertyType::permeability};
+
+ std::array<MaterialPropertyLib::PropertyType, 2> const
+ required_properties_liquid_phase = {
+ MaterialPropertyLib::PropertyType::density,
+ MaterialPropertyLib::PropertyType::viscosity};
+
+ std::array<MaterialPropertyLib::PropertyType, 3> const
+ required_properties_components = {
+ MaterialPropertyLib::PropertyType::retardation_factor,
+ MaterialPropertyLib::PropertyType::decay_rate,
+ MaterialPropertyLib::PropertyType::molecular_diffusion};
+
+ for (auto const& element : mesh.getElements())
+ {
+ auto const element_id = element->getID();
+
+ auto const& medium = *process_data.media_map->getMedium(element_id);
+ checkRequiredProperties(medium, required_properties_medium);
+
+ // check if liquid phase definition and the corresponding properties
+ // exists
+ auto const& liquid_phase = medium.phase("AqueousLiquid");
+ checkRequiredProperties(liquid_phase, required_properties_liquid_phase);
+
+ // check if components and the corresponding properties exists
+ auto const number_of_components = liquid_phase.numberOfComponents();
+ for (std::size_t component_id = 0; component_id < number_of_components;
+ ++component_id)
+ {
+ if (!liquid_phase.hasComponent(component_id))
+ {
+ OGS_FATAL(
+ "The component %u in the AqueousLiquid phase isn't "
+ "specified.",
+ component_id);
+ }
+ auto const& component = liquid_phase.component(component_id);
+ checkRequiredProperties(component, required_properties_components);
+ }
+ }
+ DBUG("Media properties verified.");
+}
+
ComponentTransportProcess::ComponentTransportProcess(
std::string name,
MeshLib::Mesh& mesh,
@@ -48,6 +103,8 @@ void ComponentTransportProcess::initializeConcreteProcess(
MeshLib::Mesh const& mesh,
unsigned const integration_order)
{
+ checkMPLProperties(mesh, _process_data);
+
const int process_id = 0;
ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
--
GitLab