set initial values of k_s, c_s to be zero

Author: keiyamamo

turtleFSI/problems/RobinBC_validation.py

@@ -53,7 +53,7 @@ def set_problem_parameters(default_variables, **namespace):
         lambda_s=lambda_s_val,              # Solid 1rst Lamé coef. [Pa]
         robin_bc = True,                    # Robin BC
         k_s = 1.0E5,                        # elastic response necesary for RobinBC
-        c_s = 1.0E2,                            # viscoelastic response necesary for RobinBC 
+        c_s = 1.0E2,                        # viscoelastic response necesary for RobinBC 
         dx_f_id=1,                          # ID of marker in the fluid domain
         dx_s_id=1002,                       # ID of marker in the solid domain
         extrapolation="laplace",            # laplace, elastic, biharmonic, no-extrapolation

turtleFSI/problems/Stenosis_2D.py

@@ -43,6 +43,8 @@ def set_problem_parameters(default_variables, **namespace):
         mu_s=mu_s_val,                       # Solid shear modulus or 2nd Lame Coef. [Pa]
         nu_s=nu_s_val,                       # Solid Poisson ratio [-]
         lambda_s=lambda_s_val,               # Solid 1rst Lamé coef. [Pa]
+        k_s = 1.0E8,                         # elastic response necesary for RobinBC
+        c_s = 1.0E2,                         # viscoelastic response necesary for RobinBC
         u_max= 0.75,                         # max inlet flow velocity value [m/s]
         p_val= 5000,                         # inner pressure for initialisation [Pa]
         vel_t_ramp= 0.2,                     # time for velocity ramp 
@@ -53,7 +55,7 @@ def set_problem_parameters(default_variables, **namespace):
         recompute=30,                        # Number of iterations before recompute Jacobian. 
         recompute_tstep=10,                  # Number of time steps before recompute Jacobian. 
         save_step=1,                         # Save frequency of files for visualisation
-        folder="stenosis_test",              # Folder where the results will be stored
+        folder="stenosis_2d",              # Folder where the results will be stored
         checkpoint_step=50,                  # checkpoint frequency
         kill_time=100000,                    # in seconds, after this time start dumping checkpoints every timestep
         save_deg=1                           # Default could be 1. 1 saves the nodal values only while 2 takes full advantage of the mide side nodes available in the P2 solution. P2 for nice visualisations
@@ -66,7 +68,7 @@ def set_problem_parameters(default_variables, **namespace):
 def get_mesh_domain_and_boundaries(folder, **namespace):
     # Import mesh file
     mesh = Mesh()
-    with XDMFFile("mesh/2D_stenosis/mesh_correct_id.xdmf") as infile:
+    with XDMFFile("mesh/Stenosis_2D/mesh_correct_id.xdmf") as infile:
         infile.read(mesh)
     # Rescale the mesh coordinated from [mm] to [m]
     x = mesh.coordinates()

turtleFSI/problems/__init__.py

@@ -57,7 +57,7 @@ def info_red(s, check=True):
     # Domain settings
     dx_f_id=1,       # Domain id of the fluid domain
     dx_s_id=2,       # Domain id of the solid domain
-    ds_s_id=None,   # Domain id of the solid boundary necessary for 
+    ds_s_id=None,   # Domain id of the solid boundary necessary for RobinBC
 
     # Meterial settings
     fluid_properties=[],
@@ -77,15 +77,15 @@ def info_red(s, check=True):
     mu_s=5.0E4,      # Shear modulus or 2nd Lame Coef. for the solid
     nu_s=0.45,       # Poisson ratio
     lambda_s=4.5E5,  # 1st Lame Coef. for the solid
-    k_s = 1.0E8,     # elastic response necesary for RobinBC
-    c_s = 1.0E2,     # viscoelastic response necesary for RobinBC 
+    k_s = 0.0,       # elastic response necesary for RobinBC
+    c_s = 0.0,       # viscoelastic response necesary for RobinBC 
     gravity=None,    # Gravitational force on the solid
     Um=0.8,          # Maximum velocity at inlet
 
     # Variational formulations
     fluid="fluid",                             # ["fluid", "no_fluid"] Turn off fluid and only solve the solid problem
     solid="solid",                             # ["solid", "no_solid"] Turn off solid and only solve the fluid problem
-    robin_bc = False,                         # [True, False] Use Robin boundary conditions for solid
+    robin_bc = False,                          # [True, False] Use Robin boundary conditions for solid
     extrapolation="laplace",                   # laplace, elastic, biharmonic, no_extrapolation
     extrapolation_sub_type="constant",         # small_constant, volume, constant, constrained_disp, constrained_disp_vel
     bc_ids=[],                                 # List of ids for weak form of biharmonic mesh lifting operator with 'constrained_disp_vel'