3 #import step "/Users/likask/MyBuild/mofem-bitbucket/mofem/meshes/HSNLF-0213.stp" heal
5 create sphere radius 2000
6 #move Volume 1 location volume 2 include_merged
8 create sphere radius 500
9 #subtract volume 2 from volume 1
10 webcut volume 1 tool volume 2
12 #block 1 name 'PotentialFlow'
14 #sideset 1 name 'NormalVelocity'
16 #sideset 2 name 'BodySurface'
19 #chop volume 2 with volume 1
22 ##split curve 7 fraction 0.5
23 ##split curve 9 fraction 0.5
24 ##split surface 9 across location vertex 9 location vertex 10
26 ##move Vertex 6 location vertex 5 include_merged
27 ##move Vertex 8 location vertex 7 include_merged
30 #block 1 name 'PotentialFlow'
32 ##sideset 2 surface 4 5 7 8 9 10
34 #sideset 2 name 'NormalVelocity'
37 #sideset 3 name 'BodySurface'
39 #volume all scheme Tetmesh
40 #volume all size auto factor 8
41 #surface 9 sizing function type skeleton scale 10 time_accuracy_level 1
42 #surface 9 sizing function type skeleton
46 ##refine surface 9 numsplit 1 bias 1.0 depth 1 smooth
49 #nodeset 1 name 'ZeroPressure'
51 #set duplicate block elements on
53 #block 2 element type tetra10
55 #save as "/Users/likask/MyBuild/mofem-bitbucket/mofem/meshes/airfoil.cub" overwrite