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