Minor

gridapapps · Nov 8, 2024 · ee08203 · ee08203
1 parent f7b6ec0
commit ee08203
Show file tree

Hide file tree

Showing 3 changed files with 11 additions and 5 deletions.
diff --git a/src/Applications/channel.jl b/src/Applications/channel.jl
@@ -70,7 +70,8 @@ function _channel(;
   np_per_level = nothing,
   rt_scaling = false,
   formulation = :mhd,
-  simplexify = false
+  simplexify = false,
+  fluid_disc = ifelse(!simplexify,:Qk_dPkm1,:SV)
 )
   @assert inlet ∈ [:parabolic,:shercliff,:constant]
   @assert initial_value ∈ [:zero,:inlet,:solve]
@@ -151,7 +152,8 @@ function _channel(;
   params[:fespaces] = Dict(
     :order => order,
     :order_j => order_j,
-    :rt_scaling => rt_scaling ? 1.0/get_mesh_size(model) : nothing
+    :rt_scaling => rt_scaling ? 1.0/get_mesh_size(model) : nothing,
+    :fluid_disc => fluid_disc,
   )
 
   # Fluid parameters

diff --git a/src/Applications/expansion.jl b/src/Applications/expansion.jl
@@ -63,6 +63,7 @@ function _expansion(;
   solid_coupling = :none,
   niter = nothing,
   convection = :newton,
+  fluid_disc = :Qk_dPkm1,
   rt_scaling = false,
   savelines = false,
   petsc_options = "",
@@ -123,7 +124,8 @@ function _expansion(;
   params[:fespaces] = Dict{Symbol,Any}(
     :order_u => order,
     :order_j => order_j,
-    :rt_scaling => rt_scaling ? 1.0/get_mesh_size(model) : nothing
+    :rt_scaling => rt_scaling ? 1.0/get_mesh_size(model) : nothing,
+    :fluid_disc => fluid_disc,
   )
 
   params[:fluid] = Dict{Symbol,Any}(

diff --git a/src/Applications/hunt.jl b/src/Applications/hunt.jl
@@ -74,7 +74,8 @@ function _hunt(;
   kmap_x = 1,
   kmap_y = 1,
   ranks_per_level = nothing,
-  simplexify = false
+  simplexify = false,
+  fluid_disc = ifelse(!simplexify,:Qk_dPkm1,:SV)
 )
   @assert formulation ∈ [:cfd,:mhd]
   @assert initial_value ∈ [:zero,:solve]
@@ -160,7 +161,8 @@ function _hunt(;
   params[:fespaces] = Dict{Symbol,Any}(
     :order_u => order,
     :order_j => order_j,
-    :rt_scaling => rt_scaling ? 1.0/get_mesh_size(model) : nothing
+    :rt_scaling => rt_scaling ? 1.0/get_mesh_size(model) : nothing,
+    :fluid_disc => fluid_disc
   )
 
   # Boundary conditions