1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 // ***************************************************************************** /*! \file src/PDE/MultiMat/Problem/SodShocktube.hpp \copyright 2012-2015 J. Bakosi, 2016-2018 Los Alamos National Security, LLC., 2019-2021 Triad National Security, LLC. All rights reserved. See the LICENSE file for details. \brief Problem configuration for Sod's shock-tube \details This file defines a policy class for the multi-material compressible flow equations, defined in PDE/MultiMat/MultiMat.hpp. See PDE/MultiMat/Problem.hpp for general requirements on Problem policy classes for MultiMat. */ // ***************************************************************************** #ifndef MultiMatProblemSodShocktube_h #define MultiMatProblemSodShocktube_h #include #include "Types.hpp" #include "Fields.hpp" #include "FunctionPrototypes.hpp" #include "SystemComponents.hpp" #include "Inciter/Options/Problem.hpp" namespace inciter { //! MultiMat system of PDEs problem: Sod shock-tube //! \see G. A. Sod. A Survey of Several Finite Difference Methods for Systems of //! Nonlinear Hyperbolic Conservation Laws. J. Comput. Phys., 27:1–31, 1978. class MultiMatProblemSodShocktube { protected: using ncomp_t = tk::ctr::ncomp_t; using eq = tag::multimat; public: //! Initialize numerical solution static tk::InitializeFn::result_type initialize( ncomp_t system, ncomp_t ncomp, tk::real x, tk::real, tk::real, tk::real ); //! Evaluate analytical solution at (x,y,z,t) for all components static std::vector< tk::real > analyticSolution( ncomp_t system, ncomp_t ncomp, tk::real x, tk::real y, tk::real z, tk::real t ) { return initialize( system, ncomp, x, y, z, t ); } //! Compute and return source term for this problem static tk::MultiMatSrcFn::result_type src( ncomp_t, ncomp_t ncomp, tk::real, tk::real, tk::real, tk::real ) { std::vector< tk::real > s( ncomp, 0.0 ); }<--- Variable 's' is assigned a value that is never used. //! Return names of integral variables to be output to diagnostics file static std::vector< std::string > names( ncomp_t ); //! Return problem type static ctr::ProblemType type() noexcept { return ctr::ProblemType::SOD_SHOCKTUBE; } }; } // inciter:: #endif // MultiMatProblemSodShocktube_h