3d_post_processing

<< Click to Display Table of Contents >>

Navigation:  Sample Problems > Usage > Import_Export >

3d_post_processing

Previous pageReturn to chapter overviewNext page

{ 3D_POST_PROCESSING.PDE

 

 This example demonstrates the use of the TRANSFERMESH facility to import

 both data and mesh structure from 3D_MESH_EXPORT.PDE and perform

 post-processing without gridding or solving any equations.

 

 This is easily accomplished in a step-wise process:

 1) make a copy of the script that generated the exported data

 2) remove the VARIABLES and EQUATIONS sections

 3) remove any boundary conditions stated in the BOUNDARIES section

 4) add the TRANSFERMESH statement in the DEFINITIONS section

 5) add any new plots that you desire

 

 Note that the domain structure must exactly match that of the exporting problem.

 

3D_MESH_EXPORT.PDE must be run before running this problem.

 

}  

 

title 'Using TRANSFERMESH for post-processing'  

 

coordinates cartesian3  

 

definitions  

   k = 1  

   an = pi/4         { this is the angular size of the repeated segment }  

   crot = cos(an)     { the sine and cosine needed in the transformation }  

   srot = sin(an)  

   H = 0  

   xc = 1.5  

   yc = 0.2  

   rc = 0.1  

  transfermesh("3d_mesh_export_output/mesh3u.xfr",U)

 

extrusion z=0,0.4,0.6,1  

 

boundaries  

  Region 1  

 

      start(1,0) line to (2,0)  

      arc(center=0,0) to (2*crot,2*srot)  

      line to (crot,srot)  

      arc(center= 0,0) to close  

 

  Limited Region 2  

      layer 2 H=1  

      start(xc-rc,0) line to (xc+rc,0) to (xc+rc,rc) to (xc-rc,rc) to close  

 

  Limited Region 3  

      layer 2 H=-1  

      start((xc-rc)*crot,(xc-rc)*srot)  

      line to ((xc+rc)*crot,(xc+rc)*srot)  

            to ((xc+rc)*crot+rc*srot,(xc+rc)*srot-rc*crot)  

            to ((xc-rc)*crot+rc*srot,(xc-rc)*srot-rc*crot) to close  

 

plots  

  grid(x,y,z)  

  grid(x,y) on z=0.5

  contour(u) on z=0.5 zoom(1.3,0,0.4,0.4)  

  contour(u) on z=0.5 zoom(1.4,0,0.2,0.2) paint

 

end