POLYMIN - University of Waterloo

POLYMIN 2005
Option for running POLYMIN with HYDRUS2D:
If running POLYMIN from a HYDRUS flow simulation, Steps 1 & 2 above are replaced
with the following steps as shown below:
1: Run HYDRUS2D to get these files: meshtria.txt, th.txt, v.txt
which contain, respectively, the grid, moisture content and velocities
The files th.txt and v.txt may contain data for several time levels. The number of time
levels should be kept less than or equal to 13 as POLYMIN will only read a maximum of
13, and will only use the data from the last time level.
2: Edit the file process_hydrus-polymin.in to set default host and secondary material
properties, and to set the nodes for the breakthrough points:
host material: dgrav,porgrav,fsgrav (diameter, porosity, sulphide fraction)
secondary: dsable,porsable,fssable
breakthrough nodes: ipoint(i), i=1,10
The host material is the default material, initially assigned to all elements. The secondary
material properties are assigned to elements that lie within quadrilaterals defined in the file:
geofile.dat – which is created in the next step.
3: Create the file: "geofile.dat" which contains the coordinates of at least one quadrilateral
which defines the locations of the secondary material
Note: you can create "geofile.dat" manually, or by using Tecplot.
If using Tecplot:
(a) follow steps 1-2, then run the processing code to get a Tecplot file of the grid.
(b) use the Tecplot file to plot the mesh, then define the quadrilaterals using the
mouse and export them to a geometry file, and rename this file to geofile.dat
(c) edit the geofile.dat file to remove unecessary text
Here is a sample format for the file geofile.dat:
gs112dH
; title
GEOMETRY
0.00, 20.00 ;origin of quadrilateral 1 – vertex 1
0.000 -0.5 ;vertex 2 (dx,dz)
6. -0.5 ;vertex 3
5.000 0.00 ;vertex 4
0.0 0.0 ;vertex 1 (dx=0, dz=0)
GEOMETRY ;origin of quadrilateral 2 – vertex 1
0.00, 10.01
0.000 -0.5
31. -0.5
30. 0.00
0.0 0.0 ;vertex 1 (dx=0, dz=0)
Up to 100 quadrilaterals can be defined. Note that the coordinates for the quadrilateral vertices
are relative to each origin, i.e. the coordinates are given as dx, dz relative to vertex.
38