<< Click to Display Table of Contents >> 3d_planespec |
{ 3D_PLANESPEC.PDE
This problem demonstrates the use of the PLANE generating function in
3D domain specifications.
We construct a hexahedron using two PLANE statements. One from explicit
points and another from named points in the definitions section.
}
title 'PLANE surface generation'
coordinates
cartesian3
variables
Tp
select
ngrid=5 { reduce mesh size for example }
definitions
long = 1
wide = 1
K = 1
Q = 10*exp(-x^2-y^2-z^2)
{ define three points in the plane surface }
bll = point(-1,-1,0)
blr = point(1,-1,0.2)
bul = point(-1,1,0.3)
initial values
Tp = 0.
equations
Tp: div(k*grad(Tp)) + Q = 0
extrusion
{ bottom surface using named points }
surface 'bottom' z = PLANE(bll,blr,bul)
{ top surface using explicit points }
surface 'top' z = PLANE((-1,-1,1), (1,-1,1.2), (1,1,2))
boundaries
surface 1 value(Tp)=0
surface 2 value(Tp)=0
Region 1
start(-wide,-wide)
value(Tp) = 0
line to (wide,-wide)
to (wide,wide)
to (-wide,wide)
to close
monitors
grid(x,z) on y=0
plots
grid(x,y,z) viewpoint(-7,-9,10)
grid(x,z) on y=0
contour(Tp) on y=0 as "ZX Temp"
contour(Tp) on x=0 as "YZ Temp"
end