{
This problem demonstrates automatic mesh densification due to curvature and
proximity to small features.
}
title '3D Layer Pinch-out Test'
coordinates
cartesian3
select
paintregions
variables
Tp
definitions
long = 1
wide = 1
K = 1
Q = 0
narrow = 0.2
z1 = 0
z2 = 0.1+0.3*exp(-(x^2+y^2)/narrow^2)
z3 = 0.5
z4= 1
initial values
Tp = 0.
equations
div[k*grad(Tp)] + Q = 0
extrusion z = z1,z2,z3,z4
boundaries
surface 1 value (Tp)=0
surface 4 value (Tp)=1
Region 1
layer 1 k=10
layer 3 k=5
start (-wide,-wide)
line to (wide,-wide) to (wide,wide) to (-wide,wide) to close
monitors
grid (x,z) on y=0
contour (Tp) on z=0.38 painted
plots
grid(x,y,z) on layer 1
grid (x,z) on y=0
grid(x,y) on surface 2
contour (Tp) on y=0 as "ZX Temp"
contour (Tp) on z=0.38 painted
end