@@ -9,14 +9,15 @@ grid = "cascade_3D_periodic_2p5mm.unv"
99mesh_file = joinpath (grids_dir, grid)
1010mesh = UNV3D_mesh (mesh_file, scale= 0.001 )
1111
12- backend = CUDABackend (); workgroup = 32
13- # backend = CPU(); workgroup = 1024; activate_multithread(backend)
12+ # backend = CUDABackend(); workgroup = 32
13+ backend = CPU (); workgroup = 1024 ; activate_multithread (backend)
1414
1515hardware = Hardware (backend= backend, workgroup= workgroup)
1616mesh_dev = adapt (backend, mesh)
1717
1818periodic1 = construct_periodic (mesh, backend, :top , :bottom )
19- periodic2 = construct_periodic (mesh, backend, :side1 , :side2 )
19+ # periodic2 = construct_periodic(mesh, backend, :side1, :side2)
20+ periodic2 = Symmetry .([:side1 , :side2 ])
2021
2122velocity = [0.25 , 0.0 , 0.0 ]
2223nu = 1e-3
@@ -59,7 +60,7 @@ BCs= assign(
5960 )
6061)
6162
62- divergence = LUST # Upwind Linear LUST
63+ divergence = Linear # Upwind Linear LUST
6364schemes = (
6465 # # transient schemes
6566 # U = Schemes(time=Euler, divergence=divergence, gradient=Gauss),
0 commit comments