Transient nsi

Project.toml

@@ -5,6 +5,7 @@ version = "0.10.0"
 
 [deps]
 Gridap = "56d4f2e9-7ea1-5844-9cf6-b9c51ca7ce8e"
+GridapODEs = "55e38337-5b6e-4c7c-9cfc-e00dd49102e6"
 IJulia = "7073ff75-c697-5162-941a-fcdaad2a7d2a"
 LineSearches = "d3d80556-e9d4-5f37-9878-2ab0fcc64255"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
@@ -13,6 +14,7 @@ NLsolve = "2774e3e8-f4cf-5e23-947b-6d7e65073b56"
 Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
+WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192"
 
 [compat]
 Gridap = "0.10"

deps/build.jl

@@ -15,7 +15,8 @@ files = [
   "Incompressible Navier-Stokes"=>"inc_navier_stokes.jl",
   "Stokes equation" => "stokes.jl",
   "Isotropic damage model" => "isotropic_damage.jl",
-  "Fluid-Structure Interaction"=>"fsi_tutorial.jl"]
+  "Fluid-Structure Interaction"=>"fsi_tutorial.jl",
+  "Laminar flow around a cylinder"=>"transient_inc_navier_stokes.jl"]
 
 Sys.rm(notebooks_dir;recursive=true,force=true)
 for (i,(title,filename)) in enumerate(files)

models/cylinder_NSI.geo

@@ -0,0 +1,23 @@
+// Gmsh project created on Fri Mar 20 14:01:03 2020
+SetFactory("OpenCASCADE");
+
+// Define Geometry
+Rectangle(1) = {0, 0, 0, 2.2, 0.41, 0};
+Disk(2) = {0.2, 0.2, 0, 0.05, 0.05};
+BooleanDifference{ Surface{1}; Delete; }{ Surface{2}; Delete; }
+Coherence;
+
+// Define mesh sizes
+Transfinite Curve {7} = 40 Using Progression 1;
+Transfinite Curve {5} = 40 Using Progression 1;
+Transfinite Curve {8} = 20 Using Progression 1;
+Transfinite Curve {6, 9} = 80 Using Progression 1;
+
+// Define Physical groups//+
+Physical Surface("fluid", 1) = {1};
+Physical Curve("inlet", 2) = {7};
+Physical Curve("outlet", 3) = {8};
+Physical Curve("noslip", 4) = {6, 9};
+Physical Curve("cylinder", 5) = {5};
+Physical Point("noslip", 4) = {6, 8, 9, 7};
+Physical Point("cylinder", 5) = {5};

models/cylinder_NSI.jl

@@ -0,0 +1,10 @@
+using Gridap
+using Gridap.Io
+using GridapGmsh
+
+model = GmshDiscreteModel("cylinder_NSI.msh")
+
+writevtk(model,"cilinder_NSI")
+
+fn = "cylinder_NSI.json"
+to_json_file(model,fn)