Add 2D LRMP with DG discretization

cadet · Dec 4, 2024 · ef257ea · ef257ea
1 parent abd22a2
commit ef257ea
Show file tree

Hide file tree

Showing 17 changed files with 7,960 additions and 212 deletions.
diff --git a/src/libcadet/CMakeLists.txt b/src/libcadet/CMakeLists.txt
@@ -216,6 +216,14 @@ if (ENABLE_DG)
 		${CMAKE_SOURCE_DIR}/src/libcadet/model/LumpedRateModelWithoutPoresDG.cpp
 		${CMAKE_SOURCE_DIR}/src/libcadet/model/GeneralRateModelDG.cpp
 	)
+	if (ENABLE_2D_MODELS)
+		list (APPEND LIBCADET_SOURCES
+			${CMAKE_SOURCE_DIR}/src/libcadet/model/parts/TwoDimensionalConvectionDispersionOperatorDG.cpp
+			${CMAKE_SOURCE_DIR}/src/libcadet/model/LumpedRateModelWithPoresDG2D.cpp
+			${CMAKE_SOURCE_DIR}/src/libcadet/model/LumpedRateModelWithPoresDG2D-LinearSolver.cpp
+			${CMAKE_SOURCE_DIR}/src/libcadet/model/LumpedRateModelWithPoresDG2D-InitialConditions.cpp
+		)
+	endif()
 endif()
 
 # Preprocess binding and reaction models

diff --git a/src/libcadet/ModelBuilderImpl.cpp b/src/libcadet/ModelBuilderImpl.cpp
@@ -40,6 +40,9 @@ namespace cadet
 #ifdef ENABLE_2D_MODELS
 		void registerGeneralRateModel2D(std::unordered_map<std::string, std::function<IUnitOperation*(UnitOpIdx, IParameterProvider&)>>& models);
 		void registerMultiChannelTransportModel(std::unordered_map<std::string, std::function<IUnitOperation*(UnitOpIdx, IParameterProvider&)>>& models);
+	#ifdef ENABLE_DG
+		void registerLumpedRateModelWithPoresDG2D(std::unordered_map<std::string, std::function<IUnitOperation* (UnitOpIdx, IParameterProvider&)>>& models);
+	#endif
 #endif
 
 		namespace inlet
@@ -67,7 +70,11 @@ namespace cadet
 #ifdef ENABLE_2D_MODELS
 		model::registerGeneralRateModel2D(_modelCreators);
 		model::registerMultiChannelTransportModel(_modelCreators);
+	#ifdef ENABLE_DG
+			model::registerLumpedRateModelWithPoresDG2D(_modelCreators);
+	#endif
 #endif
+
 		// Register all available inlet profiles
 		model::inlet::registerPiecewiseCubicPoly(_inletCreators);