diff --git a/test/freeflow/CMakeLists.txt b/test/freeflow/CMakeLists.txt index d041b5ec5e3c26146b12fdb969645164f1702997..c0a5176e7a88498302dda88012b67fb87dd4d851 100644 --- a/test/freeflow/CMakeLists.txt +++ b/test/freeflow/CMakeLists.txt @@ -1,7 +1,6 @@ add_subdirectory("navierstokes") add_subdirectory("staggered") add_subdirectory("staggerednc") -add_subdirectory("staggeredni") add_subdirectory("stokes") add_subdirectory("stokes2c") add_subdirectory("stokes2cni") diff --git a/test/freeflow/staggeredni/CMakeLists.txt b/test/freeflow/staggeredni/CMakeLists.txt deleted file mode 100644 index 0428efb58716ef5cbf3f3b7eaab5ad4fc7441156..0000000000000000000000000000000000000000 --- a/test/freeflow/staggeredni/CMakeLists.txt +++ /dev/null @@ -1,15 +0,0 @@ -add_input_file_links() - -add_dumux_test(test_channel_ni test_channel_ni test_channel_ni.cc ./test_channel_ni) - -add_dumux_test(test_diffusion_ni test_diffusion_ni test_diffusion_ni.cc ./test_diffusion_ni) - -set(CMAKE_BUILD_TYPE Debug) - -#install sources -install(FILES -test_channel_ni.cc -channeltestproblem.hh -test_diffusion_ni.cc -diffusiontestproblem.hh -DESTINATION ${CMAKE_INSTALL_INCLUDEDIR}/dumux/test/freeflow/staggeredni) diff --git a/test/freeflow/staggeredni/channeltestproblem.hh b/test/freeflow/staggeredni/channeltestproblem.hh deleted file mode 100644 index d0faf7462c363a99f90965f9fe37df0bc0be65fc..0000000000000000000000000000000000000000 --- a/test/freeflow/staggeredni/channeltestproblem.hh +++ /dev/null @@ -1,287 +0,0 @@ -// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- -// vi: set et ts=4 sw=4 sts=4: -/***************************************************************************** - * See the file COPYING for full copying permissions. * - * * - * This program is free software: you can redistribute it and/or modify * - * it under the terms of the GNU General Public License as published by * - * the Free Software Foundation, either version 2 of the License, or * - * (at your option) any later version. * - * * - * This program is distributed in the hope that it will be useful, * - * but WITHOUT ANY WARRANTY; without even the implied warranty of * - * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * - * GNU General Public License for more details. * - * * - * You should have received a copy of the GNU General Public License * - * along with this program. If not, see <http://www.gnu.org/licenses/>. * - *****************************************************************************/ -/*! - * \file - * - * \brief Channel flow test for the nonisothermal staggered grid (Navier-)Stokes model - */ -#ifndef DUMUX_CHANNEL_NI_TEST_PROBLEM_HH -#define DUMUX_CHANNEL_NI_TEST_PROBLEM_HH - -#include <dumux/implicit/staggered/properties.hh> -#include <dumux/freeflow/staggeredni/model.hh> -#include <dumux/implicit/problem.hh> -#include <dumux/material/components/simpleh2o.hh> -#include <dumux/material/fluidsystems/liquidphase.hh> - -namespace Dumux -{ -template <class TypeTag> -class ChannelNITestProblem; - -namespace Capabilities -{ - template<class TypeTag> - struct isStationary<ChannelNITestProblem<TypeTag>> - { static const bool value = false; }; -} - -namespace Properties -{ -NEW_TYPE_TAG(ChannelNITestProblem, INHERITS_FROM(StaggeredModel, NavierStokesNI)); - -SET_PROP(ChannelNITestProblem, Fluid) -{ -private: - typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar; -public: - typedef FluidSystems::LiquidPhase<Scalar, Dumux::SimpleH2O<Scalar> > type; -}; - -// Set the grid type -SET_TYPE_PROP(ChannelNITestProblem, Grid, Dune::YaspGrid<2>); - -// Set the problem property -SET_TYPE_PROP(ChannelNITestProblem, Problem, Dumux::ChannelNITestProblem<TypeTag> ); - -SET_BOOL_PROP(ChannelNITestProblem, EnableGlobalFVGeometryCache, true); - -SET_BOOL_PROP(ChannelNITestProblem, EnableGlobalFluxVariablesCache, true); -SET_BOOL_PROP(ChannelNITestProblem, EnableGlobalVolumeVariablesCache, true); - -// Enable gravity -SET_BOOL_PROP(ChannelNITestProblem, ProblemEnableGravity, true); - -//#if ENABLE_NAVIERSTOKES -SET_BOOL_PROP(ChannelNITestProblem, EnableInertiaTerms, true); -//#else -//SET_BOOL_PROP(ChannelNITestProblem, EnableInertiaTerms, false); -//#endif -} - -/*! - * \brief Test problem for the one-phase model: - \todo doc me! - */ -template <class TypeTag> -class ChannelNITestProblem : public NavierStokesProblem<TypeTag> -{ - typedef NavierStokesProblem<TypeTag> ParentType; - - typedef typename GET_PROP_TYPE(TypeTag, GridView) GridView; - typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar; - using FluidSystem = typename GET_PROP_TYPE(TypeTag, FluidSystem); - - // copy some indices for convenience - typedef typename GET_PROP_TYPE(TypeTag, Indices) Indices; - enum { - // Grid and world dimension - dim = GridView::dimension, - dimWorld = GridView::dimensionworld - }; - enum { - massBalanceIdx = Indices::massBalanceIdx, - momentumBalanceIdx = Indices::momentumBalanceIdx, - momentumXBalanceIdx = Indices::momentumXBalanceIdx, - momentumYBalanceIdx = Indices::momentumYBalanceIdx, - energyBalanceIdx = Indices::energyBalanceIdx, - pressureIdx = Indices::pressureIdx, - velocityXIdx = Indices::velocityXIdx, - velocityYIdx = Indices::velocityYIdx, - temperatureIdx = Indices::temperatureIdx - }; - - typedef typename GET_PROP_TYPE(TypeTag, BoundaryTypes) BoundaryTypes; - typedef typename GET_PROP_TYPE(TypeTag, TimeManager) TimeManager; - - typedef typename GridView::template Codim<0>::Entity Element; - typedef typename GridView::Intersection Intersection; - - typedef typename GET_PROP_TYPE(TypeTag, FVElementGeometry) FVElementGeometry; - typedef typename GET_PROP_TYPE(TypeTag, SubControlVolume) SubControlVolume; - - typedef Dune::FieldVector<Scalar, dimWorld> GlobalPosition; - - using CellCenterPrimaryVariables = typename GET_PROP_TYPE(TypeTag, CellCenterPrimaryVariables); - using FacePrimaryVariables = typename GET_PROP_TYPE(TypeTag, FacePrimaryVariables); - - using BoundaryValues = typename GET_PROP_TYPE(TypeTag, BoundaryValues); - using InitialValues = typename GET_PROP_TYPE(TypeTag, BoundaryValues); - using SourceValues = typename GET_PROP_TYPE(TypeTag, BoundaryValues); - -public: - ChannelNITestProblem(TimeManager &timeManager, const GridView &gridView) - : ParentType(timeManager, gridView) - { - name_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, - std::string, - Problem, - Name); - - inletVelocity_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, - Scalar, - Problem, - InletVelocity); - inletTemp_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, - Scalar, - Problem, - InletTemp); - FluidSystem::init(); - } - - /*! - * \name Problem parameters - */ - // \{ - - /*! - * \brief The problem name. - * - * This is used as a prefix for files generated by the simulation. - */ - std::string name() const - { - return name_; - } - - bool shouldWriteRestartFile() const - { - return false; - } - - /*! - * \brief Return the sources within the domain. - * - * \param globalPos The global position - */ - SourceValues sourceAtPos(const GlobalPosition &globalPos) const - { - return SourceValues(0.0); - } - // \} - /*! - * \name Boundary conditions - */ - // \{ - - /*! - * \brief Specifies which kind of boundary condition should be - * used for which equation on a given boundary control volume. - * - * \param globalPos The position of the center of the finite volume - */ - BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const - { - BoundaryTypes values; - // set Dirichlet values for the velocity everywhere - values.setDirichlet(momentumBalanceIdx); - values.setDirichlet(energyBalanceIdx); - - // set a fixed pressure in one cell - if (isOutlet(globalPos)) - { - values.setDirichlet(massBalanceIdx); - values.setOutflow(momentumBalanceIdx); - } - else - values.setOutflow(massBalanceIdx); - - return values; - } - - /*! - * \brief Evaluate the boundary conditions for a dirichlet - * control volume. - * - * \param globalPos The center of the finite volume which ought to be set. - */ - BoundaryValues dirichletAtPos(const GlobalPosition &globalPos) const - { - BoundaryValues values; - values[pressureIdx] = 1.1e+5; - values[temperatureIdx]= 283.15; - - if(isInlet(globalPos)) - { - values[velocityXIdx] = inletVelocity_; - values[velocityYIdx] = 0.0; - values[temperatureIdx] = inletTemp_; - } - else if(isWall(globalPos)) - { - values[velocityXIdx] = 0.0; - values[velocityYIdx] = 0.0; - } - else if(isOutlet(globalPos)) - { - values[velocityXIdx] = 1.0; - values[velocityYIdx]= 0.0; - } - return values; - } - - // \} - - /*! - * \name Volume terms - */ - // \{ - - /*! - * \brief Evaluate the initial value for a control volume. - * - * \param globalPos The global position - */ - InitialValues initialAtPos(const GlobalPosition &globalPos) const - { - InitialValues values; - values[pressureIdx] = 1.1e+5; - values[velocityXIdx] = 0.0; - values[velocityYIdx] = 0.0; - values[temperatureIdx] = 283.15; - - return values; - } - - // \} - -private: - - bool isInlet(const GlobalPosition& globalPos) const - { - return globalPos[0] < eps_; - } - - bool isOutlet(const GlobalPosition& globalPos) const - { - return globalPos[0] > this->bBoxMax()[0] - eps_; - } - - bool isWall(const GlobalPosition& globalPos) const - { - return globalPos[0] > eps_ || globalPos[0] < this->bBoxMax()[0] - eps_; - } - - const Scalar eps_{1e-6}; - Scalar inletVelocity_, inletTemp_; - std::string name_; -}; -} //end namespace - -#endif diff --git a/test/freeflow/staggeredni/diffusiontestproblem.hh b/test/freeflow/staggeredni/diffusiontestproblem.hh deleted file mode 100644 index a395416384332a11612a99cf3126852ff519ac07..0000000000000000000000000000000000000000 --- a/test/freeflow/staggeredni/diffusiontestproblem.hh +++ /dev/null @@ -1,273 +0,0 @@ -// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- -// vi: set et ts=4 sw=4 sts=4: -/***************************************************************************** - * See the file COPYING for full copying permissions. * - * * - * This program is free software: you can redistribute it and/or modify * - * it under the terms of the GNU General Public License as published by * - * the Free Software Foundation, either version 2 of the License, or * - * (at your option) any later version. * - * * - * This program is distributed in the hope that it will be useful, * - * but WITHOUT ANY WARRANTY; without even the implied warranty of * - * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * - * GNU General Public License for more details. * - * * - * You should have received a copy of the GNU General Public License * - * along with this program. If not, see <http://www.gnu.org/licenses/>. * - *****************************************************************************/ -/*! - * \file - * - * \brief Diffusion test for the nonisothermal staggered grid (Navier-)Stokes model - */ -#ifndef DUMUX_DIFFUSION_NI_TEST_PROBLEM_HH -#define DUMUX_DIFFUSION_NI_TEST_PROBLEM_HH - -#include <dumux/implicit/staggered/properties.hh> -#include <dumux/freeflow/staggeredni/model.hh> -#include <dumux/implicit/problem.hh> -#include <dumux/material/fluidsystems/liquidphase.hh> -#include <dumux/material/components/simpleh2o.hh> - -namespace Dumux -{ -template <class TypeTag> -class DiffusionNITestProblem; - -namespace Capabilities -{ - template<class TypeTag> - struct isStationary<DiffusionNITestProblem<TypeTag>> - { static const bool value = false; }; -} - -namespace Properties -{ -NEW_TYPE_TAG(DiffusionNITestProblem, INHERITS_FROM(StaggeredModel, NavierStokesNI)); - -SET_PROP(DiffusionNITestProblem, Fluid) -{ -private: - typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar; -public: - typedef FluidSystems::LiquidPhase<Scalar, Dumux::SimpleH2O<Scalar> > type; -}; - -// Set the grid type -SET_TYPE_PROP(DiffusionNITestProblem, Grid, Dune::YaspGrid<2>); - -// Set the problem property -SET_TYPE_PROP(DiffusionNITestProblem, Problem, Dumux::DiffusionNITestProblem<TypeTag> ); - -SET_BOOL_PROP(DiffusionNITestProblem, EnableGlobalFVGeometryCache, true); - -SET_BOOL_PROP(DiffusionNITestProblem, EnableGlobalFluxVariablesCache, true); -SET_BOOL_PROP(DiffusionNITestProblem, EnableGlobalVolumeVariablesCache, true); - -// Enable gravity -SET_BOOL_PROP(DiffusionNITestProblem, ProblemEnableGravity, true); - -//#if ENABLE_NAVIERSTOKES -SET_BOOL_PROP(DiffusionNITestProblem, EnableInertiaTerms, true); -//#else -//SET_BOOL_PROP(DiffusionNITestProblem, EnableInertiaTerms, false); -//#endif -} - -/*! - * \brief Test problem for the one-phase model: - \todo doc me! - */ -template <class TypeTag> -class DiffusionNITestProblem : public NavierStokesProblem<TypeTag> -{ - typedef NavierStokesProblem<TypeTag> ParentType; - - typedef typename GET_PROP_TYPE(TypeTag, GridView) GridView; - typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar; - using FluidSystem = typename GET_PROP_TYPE(TypeTag, FluidSystem); - - // copy some indices for convenience - typedef typename GET_PROP_TYPE(TypeTag, Indices) Indices; - enum { - // Grid and world dimension - dim = GridView::dimension, - dimWorld = GridView::dimensionworld - }; - enum { - massBalanceIdx = Indices::massBalanceIdx, - momentumBalanceIdx = Indices::momentumBalanceIdx, - momentumXBalanceIdx = Indices::momentumXBalanceIdx, - momentumYBalanceIdx = Indices::momentumYBalanceIdx, - energyBalanceIdx = Indices::energyBalanceIdx, - pressureIdx = Indices::pressureIdx, - velocityXIdx = Indices::velocityXIdx, - velocityYIdx = Indices::velocityYIdx, - temperatureIdx = Indices::temperatureIdx - }; - - typedef typename GET_PROP_TYPE(TypeTag, BoundaryTypes) BoundaryTypes; - typedef typename GET_PROP_TYPE(TypeTag, TimeManager) TimeManager; - - typedef typename GridView::template Codim<0>::Entity Element; - typedef typename GridView::Intersection Intersection; - - typedef typename GET_PROP_TYPE(TypeTag, FVElementGeometry) FVElementGeometry; - typedef typename GET_PROP_TYPE(TypeTag, SubControlVolume) SubControlVolume; - - typedef Dune::FieldVector<Scalar, dimWorld> GlobalPosition; - - using CellCenterPrimaryVariables = typename GET_PROP_TYPE(TypeTag, CellCenterPrimaryVariables); - using FacePrimaryVariables = typename GET_PROP_TYPE(TypeTag, FacePrimaryVariables); - - using BoundaryValues = typename GET_PROP_TYPE(TypeTag, BoundaryValues); - using InitialValues = typename GET_PROP_TYPE(TypeTag, BoundaryValues); - using SourceValues = typename GET_PROP_TYPE(TypeTag, BoundaryValues); - -public: - DiffusionNITestProblem(TimeManager &timeManager, const GridView &gridView) - : ParentType(timeManager, gridView) - { - name_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, - std::string, - Problem, - Name); - } - - /*! - * \name Problem parameters - */ - // \{ - - /*! - * \brief The problem name. - * - * This is used as a prefix for files generated by the simulation. - */ - std::string name() const - { - return name_; - } - - bool shouldWriteRestartFile() const - { - return false; - } - - /*! - * \brief Return the sources within the domain. - * - * \param globalPos The global position - */ - SourceValues sourceAtPos(const GlobalPosition &globalPos) const - { - return SourceValues(0.0); - } - // \} - /*! - * \name Boundary conditions - */ - // \{ - - /*! - * \brief Specifies which kind of boundary condition should be - * used for which equation on a given boundary control volume. - * - * \param globalPos The position of the center of the finite volume - */ - BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const - { - BoundaryTypes values; - - if(isInlet(globalPos)) - { - values.setDirichlet(momentumBalanceIdx); - values.setOutflow(massBalanceIdx); - values.setDirichlet(energyBalanceIdx); - } - else if(isOutlet(globalPos)) - { - values.setOutflow(momentumBalanceIdx); - values.setDirichlet(massBalanceIdx); - values.setOutflow(energyBalanceIdx); - } - - else - { - // set Dirichlet values for the velocity everywhere - values.setDirichlet(momentumBalanceIdx); - values.setOutflow(massBalanceIdx); - values.setOutflow(energyBalanceIdx); - } - - return values; - } - - /*! - * \brief Evaluate the boundary conditions for a dirichlet - * control volume. - * - * \param globalPos The center of the finite volume which ought to be set. - */ - BoundaryValues dirichletAtPos(const GlobalPosition &globalPos) const - { - BoundaryValues values = initialAtPos(globalPos); - - if(isInlet(globalPos)) - { - values[energyBalanceIdx] = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Problem, InletTemp); - } - - return values; - } - - // \} - - /*! - * \name Volume terms - */ - // \{ - - /*! - * \brief Evaluate the initial value for a control volume. - * - * \param globalPos The global position - */ - InitialValues initialAtPos(const GlobalPosition &globalPos) const - { - InitialValues values; - values[pressureIdx] = 1.1e+5; - values[velocityXIdx] = 0.0; - values[velocityYIdx] = 0.0; - values[temperatureIdx] = 283.15; - - return values; - } - - // \} - -private: - - bool isInlet(const GlobalPosition& globalPos) const - { - return globalPos[0] < eps_; - } - - bool isOutlet(const GlobalPosition& globalPos) const - { - return globalPos[0] > this->bBoxMax()[0] - eps_; - } - - bool isWall(const GlobalPosition& globalPos) const - { - return globalPos[0] > eps_ || globalPos[0] < this->bBoxMax()[0] - eps_; - } - - const Scalar eps_{1e-6}; - Scalar inletVelocity_; - std::string name_; -}; -} //end namespace - -#endif diff --git a/test/freeflow/staggeredni/test_channel_ni.cc b/test/freeflow/staggeredni/test_channel_ni.cc deleted file mode 100644 index 345d91423ee7efde68dbe53feb8709e36e6c424e..0000000000000000000000000000000000000000 --- a/test/freeflow/staggeredni/test_channel_ni.cc +++ /dev/null @@ -1,64 +0,0 @@ -// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- -// vi: set et ts=4 sw=4 sts=4: -/***************************************************************************** - * See the file COPYING for full copying permissions. * - * * - * This program is free software: you can redistribute it and/or modify * - * it under the terms of the GNU General Public License as published by * - * the Free Software Foundation, either version 2 of the License, or * - * (at your option) any later version. * - * * - * This program is distributed in the hope that it will be useful, * - * but WITHOUT ANY WARRANTY; without even the implied warranty of * - * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * - * GNU General Public License for more details. * - * * - * You should have received a copy of the GNU General Public License * - * along with this program. If not, see <http://www.gnu.org/licenses/>. * - *****************************************************************************/ -/*! - * \file - * - * \brief Test for the staggered grid non-isothermal (Navier-)Stokes model - */ -#include <config.h> -#include "channeltestproblem.hh" -#include <dumux/common/start.hh> - -/*! - * \brief Provides an interface for customizing error messages associated with - * reading in parameters. - * - * \param progName The name of the program, that was tried to be started. - * \param errorMsg The error message that was issued by the start function. - * Comprises the thing that went wrong and a general help message. - */ -void usage(const char *progName, const std::string &errorMsg) -{ - if (errorMsg.size() > 0) { - std::string errorMessageOut = "\nUsage: "; - errorMessageOut += progName; - errorMessageOut += " [options]\n"; - errorMessageOut += errorMsg; - errorMessageOut += "\n\nThe list of mandatory arguments for this program is:\n" - "\t-TimeManager.TEnd End of the simulation [s] \n" - "\t-TimeManager.DtInitial Initial timestep size [s] \n" - "\t-Grid.File Name of the file containing the grid \n" - "\t definition in DGF format\n" - "\t-SpatialParams.LensLowerLeftX x-coordinate of the lower left corner of the lens [m] \n" - "\t-SpatialParams.LensLowerLeftY y-coordinate of the lower left corner of the lens [m] \n" - "\t-SpatialParams.LensUpperRightX x-coordinate of the upper right corner of the lens [m] \n" - "\t-SpatialParams.LensUpperRightY y-coordinate of the upper right corner of the lens [m] \n" - "\t-SpatialParams.Permeability Permeability of the domain [m^2] \n" - "\t-SpatialParams.PermeabilityLens Permeability of the lens [m^2] \n"; - - std::cout << errorMessageOut - << "\n"; - } -} - -int main(int argc, char** argv) -{ - typedef TTAG(ChannelNITestProblem) ProblemTypeTag; - return Dumux::start<ProblemTypeTag>(argc, argv, usage); -} diff --git a/test/freeflow/staggeredni/test_channel_ni.input b/test/freeflow/staggeredni/test_channel_ni.input deleted file mode 100644 index b8656914852ae32359b905679a29e8e3020f8264..0000000000000000000000000000000000000000 --- a/test/freeflow/staggeredni/test_channel_ni.input +++ /dev/null @@ -1,20 +0,0 @@ -[TimeManager] -DtInitial = 5e-3 # [s] -TEnd = 1e-2 # [s] - -[Grid] -UpperRight = 5 1 -Cells = 50 25 - -[Problem] -Name = test_channel_ni # name passed to the output routines -EnableGravity = false -InletVelocity = 5 #1000 -InletTemp = 293.15 - -[ Newton ] -MaxSteps = 10 -MaxRelativeShift = 1e-8 - -[Vtk] -WriteFaceData = false diff --git a/test/freeflow/staggeredni/test_diffusion_ni.cc b/test/freeflow/staggeredni/test_diffusion_ni.cc deleted file mode 100644 index cf5b516cc5351da358098940a82d3a39bab50929..0000000000000000000000000000000000000000 --- a/test/freeflow/staggeredni/test_diffusion_ni.cc +++ /dev/null @@ -1,64 +0,0 @@ -// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- -// vi: set et ts=4 sw=4 sts=4: -/***************************************************************************** - * See the file COPYING for full copying permissions. * - * * - * This program is free software: you can redistribute it and/or modify * - * it under the terms of the GNU General Public License as published by * - * the Free Software Foundation, either version 2 of the License, or * - * (at your option) any later version. * - * * - * This program is distributed in the hope that it will be useful, * - * but WITHOUT ANY WARRANTY; without even the implied warranty of * - * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * - * GNU General Public License for more details. * - * * - * You should have received a copy of the GNU General Public License * - * along with this program. If not, see <http://www.gnu.org/licenses/>. * - *****************************************************************************/ -/*! - * \file - * - * \brief Test for the staggered grid non-isothermal (Navier-)Stokes model - */ -#include <config.h> -#include "diffusiontestproblem.hh" -#include <dumux/common/start.hh> - -/*! - * \brief Provides an interface for customizing error messages associated with - * reading in parameters. - * - * \param progName The name of the program, that was tried to be started. - * \param errorMsg The error message that was issued by the start function. - * Comprises the thing that went wrong and a general help message. - */ -void usage(const char *progName, const std::string &errorMsg) -{ - if (errorMsg.size() > 0) { - std::string errorMessageOut = "\nUsage: "; - errorMessageOut += progName; - errorMessageOut += " [options]\n"; - errorMessageOut += errorMsg; - errorMessageOut += "\n\nThe list of mandatory arguments for this program is:\n" - "\t-TimeManager.TEnd End of the simulation [s] \n" - "\t-TimeManager.DtInitial Initial timestep size [s] \n" - "\t-Grid.File Name of the file containing the grid \n" - "\t definition in DGF format\n" - "\t-SpatialParams.LensLowerLeftX x-coordinate of the lower left corner of the lens [m] \n" - "\t-SpatialParams.LensLowerLeftY y-coordinate of the lower left corner of the lens [m] \n" - "\t-SpatialParams.LensUpperRightX x-coordinate of the upper right corner of the lens [m] \n" - "\t-SpatialParams.LensUpperRightY y-coordinate of the upper right corner of the lens [m] \n" - "\t-SpatialParams.Permeability Permeability of the domain [m^2] \n" - "\t-SpatialParams.PermeabilityLens Permeability of the lens [m^2] \n"; - - std::cout << errorMessageOut - << "\n"; - } -} - -int main(int argc, char** argv) -{ - typedef TTAG(DiffusionNITestProblem) ProblemTypeTag; - return Dumux::start<ProblemTypeTag>(argc, argv, usage); -} diff --git a/test/freeflow/staggeredni/test_diffusion_ni.input b/test/freeflow/staggeredni/test_diffusion_ni.input deleted file mode 100644 index 46ba61fd961df98e6c64ebe1790b834f075d8da2..0000000000000000000000000000000000000000 --- a/test/freeflow/staggeredni/test_diffusion_ni.input +++ /dev/null @@ -1,19 +0,0 @@ -[TimeManager] -DtInitial = 1e-0 # [s] -TEnd = 1.2e6 # [s] - -[Grid] -UpperRight = 5 5 -Cells = 15 15 - -[Problem] -Name = test_diffusion_ni # name passed to the output routines -InletTemp = 293.15 -EnableGravity = false - -[ Newton ] -MaxSteps = 10 -MaxRelativeShift = 1e-8 - -[Vtk] -WriteFaceData = false