{ 3D_TOGGLE.PDE
This problem shows the use of curved extrusion surfaces and VOID layers.
}
title '3D CYLINDRICAL VOID LAYER TEST'
coordinates
cartesian3
select
errlim = 0.005
Variables
u
definitions
k = 0.1
h = 1
L = 1
Ro = 1
Ri = Ro/2
Yc = sqrt(Ro^2-Ri^2)
Z4 = L
Z3 = sqrt(max(Ri^2-x^2,0))
Z2 = -Z3
Z1 = -L
equations
div(K*grad(u)) + h = 0
extrusion z=Z1,Z2,Z3,Z4
boundaries
Region 1
start(Ro,0)
value(u)=0 { Force U=0 on perimeter }
arc(center=0,0) angle=360 to close
limited region 2 { the hole exists only in layer 2 }
layer 2 void
start(Ri,Yc) arc(center=0,0) to (-Ri,Yc)
line to (-Ri,-Yc)
arc(center=0,0) to (Ri,-Yc)
line to close
monitors
grid(x,y,z)
elevation(u) from (-Ro,0,0) to (Ro,0,0)
contour(u) on z=0
contour(u) on y=0
plots
grid(x,y,z)
elevation(u) from (-Ro,0,0) to (Ro,0,0)
contour(u) on z=0
contour(u) on y=0
end