Branch data Line data Source code
1 : : // ***************************************************************************** 2 : : /*! 3 : : \file src/PDE/Transport/Problem/SlotCyl.hpp 4 : : \copyright 2012-2015 J. Bakosi, 5 : : 2016-2018 Los Alamos National Security, LLC., 6 : : 2019-2021 Triad National Security, LLC. 7 : : All rights reserved. See the LICENSE file for details. 8 : : \brief Problem configuration for transport equations 9 : : \details This file declares a Problem policy class for the transport 10 : : equations, defined in PDE/Transport/CGTransport.h implementing 11 : : node-centered continuous Galerkin (CG) and PDE/Transport/DGTransport.h 12 : : implementing cell-centered discontinuous Galerkin (DG) discretizations. 13 : : See PDE/Transport/Problem.h for general requirements on Problem policy 14 : : classes for cg::Transport and dg::Transport. 15 : : */ 16 : : // ***************************************************************************** 17 : : #ifndef TransportProblemSlotCyl_h 18 : : #define TransportProblemSlotCyl_h 19 : : 20 : : #include <vector> 21 : : #include <array> 22 : : 23 : : #include "Inciter/InputDeck/InputDeck.hpp" 24 : : #include "Inciter/Options/Problem.hpp" 25 : : #include "EoS/EOS.hpp" 26 : : 27 : : namespace inciter { 28 : : 29 : : //! Transport PDE problem: rotation of Zalesak's slotted cylinder 30 : : //! \see Zalesak, S. (1978) Fully Multidimensional Flux-Corrected Transport 31 : : //! Algorithms for Fluids. Journal of Computational Physics 31, 335-362. 32 : : //! \see Leveque R.J. (1996) High-Resolution Conservative Algorithms for 33 : : //! Advection in Incompressible Flow. SIAM Journal on Numerical Analyis 33, 34 : : //! 627-665. 35 : : class TransportProblemSlotCyl { 36 : : private: 37 : : using ncomp_t = tk::ncomp_t; 38 : : using eq = tag::transport; 39 : : 40 : : public: 41 : : //! Initialize numerical solution 42 : : static std::vector< tk::real > 43 : : initialize( ncomp_t ncomp, 44 : : const std::vector< EOS >& mat_blk, tk::real x, tk::real y, 45 : : tk::real, tk::real t ); 46 : : 47 : : //! Evaluate analytical solution at (x,y,z,t) for all components 48 : : static std::vector< tk::real > 49 : : analyticSolution( ncomp_t ncomp, 50 : : const std::vector< EOS >& mat_blk, tk::real x, 51 : : tk::real y, tk::real z, tk::real t ) 52 : 403420 : { return initialize( ncomp, mat_blk, x, y, z, t ); } 53 : : 54 : : //! Do error checking on PDE parameters 55 : : void errchk( ncomp_t ) const {} 56 : : 57 : : //! Assign prescribed velocity at a point 58 : : static std::vector< std::array< tk::real, 3 > > 59 : : prescribedVelocity( ncomp_t ncomp, 60 : : tk::real x, 61 : : tk::real y, 62 : : tk::real, 63 : : tk::real ); 64 : : 65 : : //! Return problem type 66 : : static ctr::ProblemType type() noexcept 67 : : { return ctr::ProblemType::SLOT_CYL; } 68 : : }; 69 : : 70 : : } // inciter:: 71 : : 72 : : #endif // TransportProblemSlotCyl_h