Commit f76510b5 authored by Santiago Ospina's avatar Santiago Ospina

using the reconstructed flux to feed the FV scheme

parent 9936351b
...@@ -14,7 +14,7 @@ private: ...@@ -14,7 +14,7 @@ private:
using RichardsSimulationTraits = Dune::Dorie::RichardsSimulationTraits<BaseTraits,RichardsOrder>; using RichardsSimulationTraits = Dune::Dorie::RichardsSimulationTraits<BaseTraits,RichardsOrder>;
using TransportSimulationTraits = Dune::Dorie::TransportSimulationTraits< using TransportSimulationTraits = Dune::Dorie::TransportSimulationTraits<
BaseTraits, BaseTraits,
typename RichardsSimulationTraits::GFWaterFlux, typename RichardsSimulationTraits::GFWaterFluxReconstruction,
typename RichardsSimulationTraits::GFSaturation typename RichardsSimulationTraits::GFSaturation
>; >;
public: public:
...@@ -41,7 +41,7 @@ private: ...@@ -41,7 +41,7 @@ private:
using GridMapper = typename Traits::GridMapper; using GridMapper = typename Traits::GridMapper;
using GFSaturationPtr = std::shared_ptr<typename RichardsSimulationTraits::GFSaturation>; using GFSaturationPtr = std::shared_ptr<typename RichardsSimulationTraits::GFSaturation>;
using GFWaterFluxPtr = std::shared_ptr<typename RichardsSimulationTraits::GFWaterFlux>; using GFWaterFluxPtr = std::shared_ptr<typename RichardsSimulationTraits::GFWaterFluxReconstruction>;
using TimeField = typename Traits::TimeField; using TimeField = typename Traits::TimeField;
using TimeInterval = typename Traits::TimeInterval; using TimeInterval = typename Traits::TimeInterval;
......
...@@ -360,7 +360,7 @@ public: ...@@ -360,7 +360,7 @@ public:
return satdgf; return satdgf;
} }
std::shared_ptr<GFWaterFlux> get_waterflux() std::shared_ptr<GFWaterFluxReconstruction> get_waterflux()
{ {
update_adapters(); update_adapters();
assert(waterfrgf); assert(waterfrgf);
......
...@@ -6,10 +6,10 @@ _test_command = run ...@@ -6,10 +6,10 @@ _test_command = run
_asset_path = "${CMAKE_CURRENT_LIST_DIR}" _asset_path = "${CMAKE_CURRENT_LIST_DIR}"
grid.gridType = rectangular grid.gridType = rectangular
grid.initialLevel = 0 grid.initialLevel = 2
grid.cells = 50 50 grid.cells = 20 20
grid.mappingFile = "{_asset_path}/maps/cell_ids.h5"
adaptivity.useAdaptivity = false grid.mappingFileDataset = grid_test_2d
#### Richards parameters #### #### Richards parameters ####
...@@ -40,4 +40,4 @@ time.end = 2E5 ...@@ -40,4 +40,4 @@ time.end = 2E5
time.maxTimestep = 1E4 time.maxTimestep = 1E4
time.startTimestep = 1E2 time.startTimestep = 1E2
numerics.timestepMethod = alex2, explicit_euler | expand timescheme numerics.timestepMethod = implicit_euler, explicit_euler | expand timescheme
\ No newline at end of file \ No newline at end of file
...@@ -39,4 +39,4 @@ transport.time.end = 2E5 ...@@ -39,4 +39,4 @@ transport.time.end = 2E5
transport.time.maxTimestep = 1E4 transport.time.maxTimestep = 1E4
transport.time.startTimestep = 1E2 transport.time.startTimestep = 1E2
transport.numerics.timestepMethod = alex2, explicit_euler | expand timescheme transport.numerics.timestepMethod = implicit_euler, explicit_euler | expand timescheme
\ No newline at end of file \ No newline at end of file
...@@ -44,4 +44,4 @@ time.end = 2E5 ...@@ -44,4 +44,4 @@ time.end = 2E5
time.maxTimestep = 1E4 time.maxTimestep = 1E4
time.startTimestep = 1E2 time.startTimestep = 1E2
numerics.timestepMethod = alex2, explicit_euler | expand timescheme numerics.timestepMethod = implicit_euler, explicit_euler | expand timescheme
\ No newline at end of file \ No newline at end of file
...@@ -40,4 +40,4 @@ time.end = 2E5 ...@@ -40,4 +40,4 @@ time.end = 2E5
time.maxTimestep = 1E4 time.maxTimestep = 1E4
time.startTimestep = 1E2 time.startTimestep = 1E2
numerics.timestepMethod = alex2, explicit_euler | expand timescheme numerics.timestepMethod = implicit_euler, explicit_euler | expand timescheme
\ No newline at end of file \ No newline at end of file
...@@ -40,4 +40,4 @@ time.end = 2E5 ...@@ -40,4 +40,4 @@ time.end = 2E5
time.maxTimestep = 1E4 time.maxTimestep = 1E4
time.startTimestep = 1E2 time.startTimestep = 1E2
numerics.timestepMethod = alex2, explicit_euler | expand timescheme numerics.timestepMethod = implicit_euler, explicit_euler | expand timescheme
\ No newline at end of file \ No newline at end of file
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment