diff --git a/src/parameters.jl b/src/parameters.jl index 05a61a9..035e160 100644 --- a/src/parameters.jl +++ b/src/parameters.jl @@ -374,6 +374,11 @@ const FLUID_DISCRETIZATIONS = (; :p_conformity => :H1, :p_order => (k) -> k-1, ), + :RT => (; # Raviart-Thomas + :p_space => :Q, + :p_conformity => :L2, + :p_order => (k) -> k-1, + ), ), :TET => (; # Tetrahedra :Pk_Pkm1 => (; # Taylor Hood @@ -395,6 +400,21 @@ const FLUID_DISCRETIZATIONS = (; ) ) +function get_fluid_reffe_pair(disc::Symbol,poly::Polytope,feparams) + if poly == HEX + + if disc == :Qk_dPkm1 + + else + + end + elseif poly == TET + + else + @notimplemented "Only HEX and TET supported for now." + end +end + """ const CURRENT_DISCRETIZATIONS diff --git a/test/seq/channel_tests.jl b/test/seq/channel_tests.jl index 6a06d30..af7445a 100644 --- a/test/seq/channel_tests.jl +++ b/test/seq/channel_tests.jl @@ -15,7 +15,7 @@ solver = Dict( channel( nc=(4,2,2), order = 3, - order_j = 4, + order_j = 3, L = 2.0, w = 1.0, ΞΆ = 50.0,