-
Notifications
You must be signed in to change notification settings - Fork 3
/
ThreeD_Donut.jl
33 lines (28 loc) · 968 Bytes
/
ThreeD_Donut.jl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
using WaterLily
using StaticArrays
function donut(p=6;Re=1e3,mem=Array,U=1)
# Define simulation size, geometry dimensions, viscosity
n = 2^p
center,R,r = SA[n/2,n/2,n/2], n/4, n/16
ν = U*R/Re
# Apply signed distance function for a torus
norm2(x) = √sum(abs2,x)
body = AutoBody() do xyz,t
x,y,z = xyz - center
norm2(SA[x,norm2(SA[y,z])-R])-r
end
# Initialize simulation and return center for flow viz
Simulation((2n,n,n),(U,0,0),R;ν,body,mem),center
end
# using CUDA
sim,center = donut()#mem=CUDA.CuArray);
dat = sim.flow.σ[inside(sim.flow.σ)] |> Array;
function ω_θ!(dat,sim,center=center)
dt, a = sim.L/sim.U, sim.flow.σ
@inside a[I] = WaterLily.ω_θ(I,(1,0,0),center,sim.flow.u)*dt
copyto!(dat,a[inside(a)])
end
include("../src/ThreeD_plots.jl")
@time makie_video!(sim,dat,ω_θ!,name="donut.mp4",duration=10,step=0.25) do obs
contour(obs, levels=[-5,5], colormap=:balance)
end