Revert "tau_m vms computation options"

This reverts commit b34f9f6.

Author: carlodev

src/Commons/Equations/StabilizationOperations.jl

@@ -29,22 +29,21 @@ function momentum_stabilization(uu, stab_coeff::TensorStabilization,simcase::Sim
     @unpack G,GG,gg = stab_coeff
     @unpack sprob = simcase
     @unpack Ci = sprob.coeff_method
-    @unpack r = sprob.coeff_method
 
     @sunpack ν, dt = simcase
 
     function τm(uun, G, GG)
-        τ₁ = Ci[1] * (2 / dt)^r #Here, you can increse the 2 if CFL high
-        τ₃ = Ci[2] * (ν^r * (GG)^r/2 )
+        τ₁ = Ci[1] * (2 / dt)^2 
+        τ₃ = Ci[2] * (ν^2 * GG )
 
         uu_new = VectorValue(val_u.(uun)...)
 
         if iszero(norm(uu_new))
-            return (τ₁ .+ τ₃) .^ (-1 / r)
+            return (τ₁ .+ τ₃) .^ (-1 / 2)
         end
 
         τ₂ = uu_new ⋅ G ⋅ uu_new
-        return (τ₁ .+ τ₂ .+ τ₃) .^ (-1 / r)
+        return (τ₁ .+ τ₂ .+ τ₃) .^ (-1 / 2)
     end
 
 

src/Commons/MatrixCreation.jl

@@ -99,7 +99,7 @@ end
 
 function initialize_matrices(u_adv, params, simcase)
   @info "allocation matrix and vectors"
-  @time matrices = allocate_all_matrices_vectors(u_adv, params, simcase)
+  matrices = allocate_all_matrices_vectors(u_adv, params, simcase)
   @info "matrix and vectors allocated"
 
   @info "updating matrix and vectors"

src/Commons/ParametersDef/StabilizationStruct.jl

@@ -27,12 +27,7 @@ end
 
 @with_kw struct TensorFormulation <: StabilizationFormulation
     r::Int64 = 2 
-    Ci::Vector{Real} = [4,36]
-    τm_comp::Real = -1 #Standard VMS; τm_comp = 1 SUPG standard
-    @assert length(Ci) == 2 "Ci length must be 2"
-    @assert (Ci[1]>=0 && Ci[2]>=0) "Ci values must be non-negative"
-    @assert τm_comp==1 ||  τm_comp==-1 " τm values +1 or -1 in TensorFormulation"
-
+    Ci::Vector{Int64} = [4,36]
 end