POLYMIN - University of Waterloo

POLYMIN 2005
5. STEP BY STEP INSTRUCTIONS
POLYMIN is a mass transport model which assumes a steady state velocity field. The velocity
field can be obtained either from the FLONET model for saturated systems (Molson & Frind,
2004), or from the HYDRUS2D model (Simunek et al., 1999). In either case, the output files
must converted to POLYMIN input format by running a “processing” code (process_flonetpoly.exe,
or process_hydrus-poly.exe).
Steps 1 & 2 below assume the flow field will be derived from the FLONET model. If running
directly from a HYDRUS flow simulation, steps 1 & 2 can be skipped and replaced with steps 1b
& 2b which follow.
1. Simulate the flow system (refer to the FLONET user guide)
• edit the fnpcg.data file
• make sure kp=1 to get required output files
• run fnpcg.exe
• transfer these file to your transport directory: flow.vxyst, flow.incid, flow.nodes
2. Convert flow system file formats for transport:
• edit process_flonet-poly.in to adjust gradius, fs2, and source nodes for oxidation
• run process_flonet-poly.exe
• you should now have these files for input to POLYMIN:
meshtria.txt, bdy.txt, v.txt, nodeprop.txt, th.txt
(mesh, boundary nodes, velocities, nodal properties & water content)
3. Prepare transport input files
• edit polymin.in (must contain the name of your input file; e.g. polymin.dat)
• edit polymin.dat (primary input file), check path to data base
• edit bkgrndpoly.min (contains background chemistry data)
• check database files:
alk.dbm, analy.dbm, compcp.dbm, error.dbm, thermcp.dbm, type6cp.dbm
4. Run polymin5.exe
5. Interpret results:
• Check if job ran completely: check polymino.gen
• Aqueous component output: polymin.ur1
(Tecplot-compatible file of x,z,O2, ph, aqueous components, etc)
• Solid species (minerals) output: polymin.ur2 (x,z, solid species)
• Edit polyminbrkaq_nn.out and polyminbrks_nn.out files – replace “999” in the
third line with the number of data lines (= number of time levels)
• Use Tecplot to visualize the above files
37