3d_ploton

<< Click to Display Table of Contents >>

Navigation:  Sample Problems > Usage > Plotting >

3d_ploton

Previous pageReturn to chapter overviewNext page

{ 3D_PLOTON.PDE

 

 This problem shows some of the possible 'ON' qualifiers for 3D plots.

 

}

 

title '3D Test -- Plot Qualifiers'

 

coordinates

   cartesian3

 

Variables

   u

 

definitions

   k = 0.1

   heat = 4

 

equations

   U:     div(K*grad(u)) + heat   = 0

 

extrusion

  surface z = 0

  surface z = 0.8-0.3*(x^2+y^2)

  surface z = 1.0-0.3*(x^2+y^2)  

 

boundaries

  region 1 'outer'

      layer 2 k = 1

      start(-1,-1)

          value(u) = 0

      line to (1,-1) to (1,1) to (-1,1) to close

  region 2 'plug'

      layer 2 k = 1

      start 'dot' (0.5,0.5) arc(center=0,0) angle=360

 

plots

    grid(x,y,z) on region 1 as "Only Region 1, both layers"

    grid(x,y,z) on region 'plug' on layer 2 as "Region 2 Layer 2"

    grid(x,y,z) on region 'plug' on layers 1,2 paintregions as "Region 2, both layers"

    grid(y,z) on x=0 on 'plug'   as "Cut plane on region 2"

    contour(u) on x=0.51   on layer 2 as "Solution on X-cut in layer 2"  

    contour(u) on z=0.51   on region 2 as "Solution on Z-cut in region 2"

    contour(u) on surface 2 on region 2         as "Solution on paraboloidal layer interface"

    vector(grad(u)) on surface 2 on 'outer' as "Flux on layer interface in region 1"

 

end