﻿ Sample Problems > Usage > 3D_domains > 3d_planespec

# 3d_planespec

Navigation:  Sample Problems > Usage > 3D_domains >

# 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