def main(case_dir='joinSTL', runRhoCentral = False): #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=20000000 #snap.refinementSurfaceMax =8 #snap.distanceLevels = [7,6,4,2] #snap.distanceRefinements = [0.01,0.02,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) snap.generate_mesh() getTrueMesh(case) if runRhoCentral: case.run_tool('rhoCentralFoam')
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()