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()
