Further developments in transient_fsi_debug.jl

Project.toml

@@ -4,20 +4,25 @@ authors = ["Santiago Badia <[email protected]>", "Francesc Verdugo <fver
 version = "0.10.0"
 
 [deps]
+ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
 Gridap = "56d4f2e9-7ea1-5844-9cf6-b9c51ca7ce8e"
 GridapODEs = "55e38337-5b6e-4c7c-9cfc-e00dd49102e6"
 IJulia = "7073ff75-c697-5162-941a-fcdaad2a7d2a"
+IncompleteLU = "40713840-3770-5561-ab4c-a76e7d0d7895"
+IterativeSolvers = "42fd0dbc-a981-5370-80f2-aaf504508153"
 LineSearches = "d3d80556-e9d4-5f37-9878-2ab0fcc64255"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 Literate = "98b081ad-f1c9-55d3-8b20-4c87d4299306"
 NLsolve = "2774e3e8-f4cf-5e23-947b-6d7e65073b56"
 Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
+Preconditioners = "af69fa37-3177-5a40-98ee-561f696e4fcd"
 Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
+Revise = "295af30f-e4ad-537b-8983-00126c2a3abe"
+TimerOutputs = "a759f4b9-e2f1-59dc-863e-4aeb61b1ea8f"
 WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192"
 
 [compat]
-Gridap = "0.11"
 julia = "1.0"
 
 [extras]

src/transient_fsi_debug.jl

@@ -128,7 +128,7 @@ const ρ_f = 1.0e3
 const μ_f = ρ_f * Um * ⌀ / Re
 const β_m = 4 # β_m = λ_m / μ_m
 const E_m = 1.0
-const ν_m = 0.2
+const ν_m = 0.25
 
 # ### Constitutive laws
 @law F(∇u) = ∇u + one(∇u)
@@ -144,7 +144,7 @@ const ν_m = 0.2
 #@law σ_m2(α,∇ut,Finv) = 16.0*α*μ_s*tr(∇ut⋅Finv)*one(Finv) + α*μ_s*(∇ut⋅Finv + (Finv')⋅(∇ut')) 
 @law conv(c,∇v) = (∇v') ⋅ c
 @law Sm(∇u) = 2*μ_m*E(∇u) + λ_m*tr(E(∇u))*one(E(∇u))
-@law α(ve,J) = 5.0e-4/(J)
+@law α(ve,J) = 1.0e-5/(J)
 @law λ_m(α) = α * (E_m*ν_m) / ((1.0+ν_m)*(1.0-2.0*ν_m))
 @law μ_m(α) = α * E_m / (2.0*(1.0+ν_m))
 #@law λ_m(J) = (E_m*ν_m) / ((1.0+ν_m)*(1.0-2.0*ν_m))
@@ -159,7 +159,7 @@ const ν_m = 0.2
 @law dFinvT(∇u,∇du) = (dFinv(∇u,∇du)')
 @law dconv(dc,∇dv,c,∇v) = conv(c,∇dv) + conv(dc,∇v)
 @law dSm(∇u,∇du) = 2*μ_m*dE(∇u,∇du) + λ_m*tr(dE(∇u,∇du))*one(E(∇u))
-@law dα(ve,J,dJ) = - 5.0e-4 * 1.0 / (J*J) * dJ
+@law dα(ve,J,dJ) = - 1.0e-5 * 1.0 / (J*J) * dJ
 @law dλ_m(dα) = dα * (E_m*ν_m) / ((1.0+ν_m)*(1.0-2.0*ν_m))
 @law dμ_m(dα) = dα * E_m / (2.0*(1.0+ν_m))
 @law dσ_m(dα,ε) = dλ_m(dα)*tr(ε)*one(ε) + 2.0*dμ_m(dα)*ε