diff --git a/doc/airfoil.png b/doc/airfoil.png index 9091cb9..3404eea 100644 Binary files a/doc/airfoil.png and b/doc/airfoil.png differ diff --git a/doc/thin_fracture.png b/doc/thin_fracture.png new file mode 100644 index 0000000..7df2570 Binary files /dev/null and b/doc/thin_fracture.png differ diff --git a/input/06_airfoil_from_csv.para b/input/06_airfoil_from_csv.para index 1c5518e..c1c0038 100644 --- a/input/06_airfoil_from_csv.para +++ b/input/06_airfoil_from_csv.para @@ -3,10 +3,10 @@ #=========================================================== #----------------------------------------------------------- -# MESH DEFINITION +# OUTER MESH DEFINITION #----------------------------------------------------------- Define mesh: - Element size: 0.004 + Element size: 0.003 Element color: 1 Meshing algorithm: Triangulation Number of quad refinements: 0 @@ -22,21 +22,21 @@ Define mesh: Define interior boundary from CSV file: ../auxiliary/test_data/Airfoil.csv - Add quad layers: 0.69132, 0.09754, 0.69132, 0.09754, 10, 0.0005, 1.10 + Add quad layers: 0.69132, 0.09754, 0.69132, 0.09754, 10, 0.0004, 1.10 #--------------------------------------------------------- # LOCAL REFINEMENT #--------------------------------------------------------- Define fixed vertices: - 0.69095,0.09625, 0.0020,0.02 - 0.85985,0.07582, 0.0008,0.005 + 0.69095, 0.09625, 0.0020, 0.02 + 0.85985, 0.07582, 0.0008, 0.005 End fixed vertices End mesh #----------------------------------------------------------- -# Mesh #3 +# INNER MESH DEFINITION #----------------------------------------------------------- Define mesh: Element size: 0.003 diff --git a/input/08_thin_fracture.para b/input/08_thin_fracture.para index f6dcd27..fbf4972 100644 --- a/input/08_thin_fracture.para +++ b/input/08_thin_fracture.para @@ -14,8 +14,8 @@ Define mesh: Number of smoothing iterations: 3 # Output formats: COUT, TXT, VTU - Output file format: TXT - Output file prefix: ./frac2d-quad + Output file format: VTU + Output file prefix: ./thin_fracture #--------------------------------------------------------- # BOUNDARIES diff --git a/src/examples/06_airfoil_from_csv.cpp b/src/examples/06_airfoil_from_csv.cpp index 0bb51fe..88bbf5e 100644 --- a/src/examples/06_airfoil_from_csv.cpp +++ b/src/examples/06_airfoil_from_csv.cpp @@ -32,12 +32,18 @@ void airfoil_from_csv() /*------------------------------------------------------------------ | Define the size function ------------------------------------------------------------------*/ - UserSizeFunction f_outer = [](const Vec2d& p) { return 0.003; }; + UserSizeFunction f_outer = [](const Vec2d& p) { return 0.004; }; Domain outer_domain { f_outer }; outer_domain.add_exterior_boundary().set_shape_circle(1, {0.77, 0.09}, 0.14, 60); outer_domain.add_interior_boundary().set_shape_from_csv(csv_file); + /*------------------------------------------------------------------ + | We use some additional fixed vertices to locally refine the mesh + ------------------------------------------------------------------*/ + outer_domain.add_fixed_vertex(0.69095, 0.09625, 0.0020, 0.02); + outer_domain.add_fixed_vertex(0.85985, 0.07582, 0.0008, 0.005); + /*------------------------------------------------------------------ | Initialize the outer mesh ------------------------------------------------------------------*/ @@ -53,7 +59,7 @@ void airfoil_from_csv() ------------------------------------------------------------------*/ generator.quad_layer_generation(outer_mesh) .n_layers(10) - .first_height(0.0003) + .first_height(0.0004) .growth_rate(1.10) .starting_position({0.69132,0.09754}) .ending_position({0.69132,0.09754})