Merge pull request #61 from MultiSimOLab/dissipation

Author: jmartfrut

src/PhysicalModels/ThermoElectroMechanicalModels.jl

@@ -153,3 +153,13 @@ struct ThermoElectroMech_Bonet{T<:Thermo,E<:Electro,M<:Mechano} <: ThermoElectro
     return (Ψ, ∂Ψ_∂F, ∂Ψ_∂E, ∂Ψ_∂δθ, ∂2Ψ_∂2F, ∂2Ψ_∂2E, ∂2Ψ_∂2δθ, ∂ΨEF, ∂ΨFδθ, ∂ΨEδθ, η)
   end
 end
+
+function Dissipation(obj::ThermoElectroMech_Bonet, Δt)
+  @unpack Cv,θr, α, κ, γv, γd = obj.thermo
+  Dvis = Dissipation(obj.mechano, Δt)
+  gd(δθ) = 1/(γd+1) * (((δθ+θr)/θr)^(γd+1) -1)
+  ∂gd(δθ) = (δθ+θr)^γd / θr^(γd+1)
+  D(F, E, δθ, A...) = (1 + gd(δθ)) * Dvis(F, A...)
+  ∂D∂θ(F, E, δθ, A...) = ∂gd(δθ) * Dvis(F, A...)
+  return(D, ∂D∂θ)
+end

src/PhysicalModels/ViscousModels.jl

@@ -497,28 +497,3 @@ function ViscousDissipation(obj::ViscousIncompressible, Δt::Float64,
   Dvis = -Se ⊙ (invCCe ⊙ ∂Se)
   Dvis
 end
-
-
-function DissipationDerivative(obj::ViscousIncompressible, Δt::Float64,
-                              Se_::Function, ∂Se∂Ce_::Function,
-                              F::TensorValue, Fn::TensorValue, A::VectorValue)
-  Uvn = TensorValue{3,3}(A[1:9]...)
-  λαn = A[10]
-  #------------------------------------------
-  # Get kinematics
-  #------------------------------------------
-  invUvn  = inv(Uvn)
-  C = Cauchy(F)
-  Cn = Cauchy(Fn)
-  Ceᵗʳ = ElasticCauchy(C, invUvn)
-  Cen  = ElasticCauchy(Cn, invUvn)
-  #------------------------------------------
-  # Return mapping algorithm
-  #------------------------------------------
-  Ce, λα = return_mapping_algorithm!(obj, Δt, Se_, ∂Se∂Ce_, C, Ceᵗʳ, Cen, λαn)
-  #------------------------------------------
-  # Dissipation derivative
-  #------------------------------------------
-  """Derivative with respect to the temperature."""
-  return 0.0
-end