def main(case_dir='coreDartStability', runRhoCentral = False, parallel = True): #Create a new case file, raise an error if the directory already exists case = create_new_case(case_dir) write_control_dict(case) #write the base block mesh make_block_mesh(case) parts = load_multiple_geometries(GeometryFormat.STL,path_list,part_list,case) snap = SnappyHexMesh(parts,4,case) ###############un-comment out the code snippet below if you want to apply extra refinement######### snap.refinementSurfaceMin =7 snap.maxGlobalCells=100000000 snap.refinementSurfaceMax =8 snap.distanceLevels = [7,6,4,2] snap.distanceRefinements = [0.010,0.020,0.06,0.2] snap.snap=True snap.snapTolerance = 8 snap.edgeRefinementLevel = 8 #############un-comment out the code snippet below if you want to apply extra refinement######### snap.locationToKeep = [0.0012,0.124,0.19] #odd numbers to ensure not on face snap.addLayers=False #we need to write fvSchemes and fvSolution to be able to use paraForm and run snappy? write_fv_schemes(case) write_fv_solution(case) write_thermophysical_properties(case) write_turbulence_properties(case) write_initial_conditions(case) write_decompose_settings(case, processors, "simple") snap.generate_mesh()
def main(case_dir='snappy'):
#Create a new case file, raise an error if the directory already exists
case = create_new_case(case_dir)
write_control_dict(case)
#write the base block mesh
make_block_mesh(case)
rocket = Geometry(GeometryFormat.STL,'STLS/example.stl','example',case)
rocket.scale(0.5);
rocket.translate([0.5,2,2])
snap = SnappyHexMesh([rocket],4,case)
snap.snap=True
snap.snapTolerance = 8;
snap.locationToKeep = [0.0012,0.124,0.19] #odd numbers to ensure not on face
snap.addLayers=False
#we need to write fvSchemes and fvSolution to be able to use paraForm and run snappy?
write_fv_schemes(case)
write_fv_solution(case)
snap.generate_mesh()