Amanzi->Walkabout Mesh 2 July 2018 hextotet 5 instead of connect, Amanzi interpolates velocities from cell centers to the tet nodes. This mesh has the same cell vertices (nodes) as the input Amanzi hex mesh. Each hex cell has been converted to 5 tets. tet5_flts_wtr.exo ExodusII mesh for Amanzi 1 ss_1_bottom.faceset 2 ss_2_south.faceset 3 ss_3_east.faceset 4 ss_4_north.faceset 5 ss_5_west.faceset 6 ss_6_top.faceset tet5_fehm.inp Walkabout sorted mesh with FEHM attributes tet5.inp sorted for walkabout, no attributes Files for FEHM and Walkabout: walkabout.tet5.fehmn walkabout.tet5.graph walkabout.tet5_material.zone walkabout.tet5.stor outside_top_back.zone outside_top.zone ExodusII tet5 mesh number of nodes = 7055100 number of elements = 34650000 element type = tet (hextotet 5) NAME MIN MAX DIFFERENCE LENGTH xic 5.350000000E+05 5.550000000E+05 2.000000000E+04 7055100 yic 4.100000000E+06 4.135000000E+06 3.500000000E+04 7055100 zic -1.000000000E+03 1.468456754E+03 2.468456754E+03 7055100 Note this hextotet5 mesh has negative ccoefs on outside nodes, but all tet volumes and node voronoi volumes are positive. *** Construct and Compress Sparse Matrix:3D *** *** Compress Area Coefficient Values *** AMatbld3d_stor: *****Negative Coefficients ****** AMatbld3d_stor: Total Number of Negative Coefficients 12730985 AMatbld3d_stor: Number of Significant Negative Coefs 12730985 AMatbld3d_stor: Number of 'zero' (< 10e-8 *max) coefs 0 AMatbld3d_stor: npoints = 7055100 ncoefs = 88523102 AMatbld3d_stor: Number of unique coefs = 4489848 AMatbld3d_stor: Maximum num. connections to a node = 19 AMatbld3d_stor: Volume min = 2.8171515E+04 AMatbld3d_stor: Volume max = 2.4928255E+05 AMatbld3d_stor: Total Volume: 1.6217985E+12 AMatbld3d_stor: abs(Aij/xij) min = 0.0000000E+00 AMatbld3d_stor: abs(Aij/xij) max = 4.4355664E+02 AMatbld3d_stor: (Aij/xij) max = 4.0941225E+00 AMatbld3d_stor: (Aij/xij) min = -4.4355664E+02 AMatbld3d_stor Matrix coefficient values stored as scalar area/distance AMatbld3d_stor Matrix compression used for graph and coefficient values ascii STOR file written with name walkabout.tet5.stor tet5_fehm.inp FEHM style with attributes: ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH imt1 1 181 180 7055100 itetclr 1 181 180 34650000 vorvol 2.817125000E+04 2.492853333E+05 2.211140833E+05 7055100 tetvol 3.756202054E+04 8.312708371E+04 4.556506317E+04 34650000 out_zone 0 6 6 7055100 --------------------------------------------------------------------------------------------- May 2017 connected particleDriver tet mesh for walkabout (cell centerss plus bndry nodes) This mesh has nodes at hex cell centers plus hex vertices at boundaries. Particle tracking has issues along the boundaries of this mesh. Walkabout files: /scratch/nts/tam/pahute_mesa/mesh_ER20_uniform/meshA_100_100_25/Amanzi_fit_watertable/walkabout Walkabout work: /scratch/nts/tam/ascem/walkabout/PM_walkaout_work Amanzi hex mesh converted to walkabout tet mesh using particleDriver and connect using points at cell centers plus boundary vertices. Tet elements connected across non-convex top surface. This mesh has many neg ccoefs and neg voronoi volumes AMatbld3d_stor: Matrix compress_eps: 0.1000000E-07 AMatbld3d_stor: Local epsilon: 0.1000000E-14 AMatbld3d_stor: *****Negative Coefficients ****** AMatbld3d_stor: Total Number of Negative Coefficients 281449 AMatbld3d_stor: Number of Significant Negative Coefs 281360 AMatbld3d_stor: npoints = 7178902 ncoefs = 54217304 AMatbld3d_stor: Number of unique coefs = 8566 AMatbld3d_stor: Maximum num. connections to a node = 27 AMatbld3d_stor: Volume min = -3.8502792E+11 AMatbld3d_stor: Volume max = 7.0088824E+11 AMatbld3d_stor: Total Volume: 1.7681825E+12 AMatbld3d_stor: abs(Aij/xij) min = 4.7430981E-06 AMatbld3d_stor: abs(Aij/xij) max = 1.1625565E+08 AMatbld3d_stor: (Aij/xij) max = 5.9501893E+07 AMatbld3d_stor: (Aij/xij) min = -1.1625565E+08 AMatbld3d_stor Matrix coefficient values stored as scalar area/distance AMatbld3d_stor Matrix compression used for graph and coefficient values ascii STOR file written with name walkabout.h5.stor Interpolate materials from Amanzi mesh on to Walkabout nodes because current Amanzi-Walkabout method does not have material properties on center/edge nodes. Use method from /scratch/nts/tam/rainier_mesa/RMSM/slice_2A_20x5m/walkabout_from_hex Work directory /scratch/nts/tam/ascem/walkabout/PM_walkaout_work Source Amanzi hex Mesh: https://meshing.lanl.gov/proj/PM_ER20_MeshA/catalog.html Mesh top directory: /scratch/nts/tam/pahute_mesa/mesh_ER20_uniform/meshA_100_100_25 This is the mesh used for the Test Bed Modeling. Hex Mesh Fit to watertable surface (smooth top): number of dimension: 3 number of nodes: 7055100 number of elements: 6930000 number of element blocks: 95 number of side sets: 6 Elevation -1000 to 1468.45 meters /scratch/nts/tam/pahute_mesa/mesh_ER20_uniform/meshA_100_100_25/Amanzi_fit_watertable/hex_flts_wtrfit.exo material blocks: 1-95 Source points (Amanzi cell center and edge points): /scratch/fwo/zhiming/PM_ASCEM/multiwells/inverse_results_from_wolf/superpara_rerun_newBounds/walk00003.h5.inp Use tetgen to connect source points into a tet mesh. Input points: 7178902 Mesh points: 7178902 Mesh tetrahedra: 47082849 Mesh faces: 94345147 Mesh edges: 54441199 Convex hull faces: 358898 Use LaGriT to interpolate materials and remove extra tets from connect. Connected tet mesh for walkabout files: number of nodes = 7178902 number of elements = 46865773 material zones: 1-181 -rw-r--r-- 1 tamiller sft 2431198794 May 15 15:07 walkabout.h5.fehmn -rw-r--r-- 1 tamiller sft 186227840 May 15 15:39 walkabout.h5_multi_mat.zone -rw-r--r-- 1 tamiller sft 2779682 May 15 15:39 walkabout.h5_outside.zone -rw-r--r-- 1 tamiller sft 16045517 May 15 15:39 walkabout.h5_outside_vor.area -rw-r--r-- 1 tamiller sft 78973907 May 16 11:29 walkabout.h5_material.zone -rw-r--r-- 1 tamiller sft 1470894826 May 15 15:50 walkabout.h5.stor -rw-r--r-- 1 tamiller sft 5951953171 May 15 15:53 walkabout.h5.graph Note: stor file does have neg ccoefs, all located on outside boundary nodes. ************* Summary Table for Cell and Node Materials ***************** EXO_blk Zone Name Cell Count Fraction/Vol Node Count Percent/Total 1 1 ATICU 475547 0.06722 496963 0.06920 2 2 RMICU 73402 0.01026 77973 0.01090 3 5 BMICU 6030 0.00090 7288 0.00100 4 8 SCICU 451 0.00007 524 0.00010 5 10 LCA 550662 0.07981 581309 0.08100 6 11 UCCU 49445 0.00730 52671 0.00730 7 14 SCVCU 249423 0.03539 259295 0.03610 8 15 PBRCM 1442648 0.21352 1484439 0.20680 9 16 BRA 486359 0.07167 500990 0.06980 10 17 BFCU 340052 0.04929 342312 0.04770 11 21 CFCU 80694 0.01167 81869 0.01140 12 22 CFCM 111456 0.01609 111491 0.01550 13 23 IA 50009 0.00721 50651 0.00710 14 24 CHLFA5 53936 0.00780 53942 0.00750 15 25 CHLFA4 27923 0.00405 28054 0.00390 16 26 CHLFA3 18837 0.00277 20160 0.00280 17 27 CHLFA2 7041 0.00104 8481 0.00120 18 28 CHLFA1 3663 0.00052 4248 0.00060 19 29 CHCU 383347 0.05572 392748 0.05470 20 30 CHVTA 4449 0.00066 5610 0.00080 21 32 TSA 20159 0.00286 20285 0.00280 22 33 PLFA 5054 0.00072 5056 0.00070 23 34 LPCU 65473 0.00939 66094 0.00920 24 35 PVTA 304 0.00004 489 0.00010 25 36 TCA 36351 0.00519 38617 0.00540 26 37 MPCU 4996 0.00071 5140 0.00070 27 38 SPA 6214 0.00088 6656 0.00090 28 39 UPCU 28471 0.00409 31054 0.00430 29 40 BA 8507 0.00120 9224 0.00130 30 41 PBPCU 671 0.00010 685 0.00010 31 42 CPA 16857 0.00240 17047 0.00240 32 44 WWA 3 0.00000 6 0.00000 33 45 FCCU 13493 0.00191 15233 0.00210 34 46 RMWTA 1034409 0.14610 1056923 0.14720 35 47 TMLVTA 32985 0.00490 38906 0.00540 36 48 TMWTA 21681 0.00313 24774 0.00350 37 49 THCU 69989 0.00989 73322 0.01020 38 50 THCM 14043 0.00199 15687 0.00220 39 51 THLFA 8837 0.00126 10008 0.00140 40 52 TMUWTA 1373 0.00019 1599 0.00020 41 53 ATCCU 7411 0.00105 7414 0.00100 42 54 ATWTA 361288 0.05083 376373 0.05240 43 55 BWCU 39497 0.00555 39843 0.00560 44 56 BWWTA 11973 0.00168 12116 0.00170 45 59 FCLLFA 219 0.00003 306 0.00000 46 61 FCWTA1 7904 0.00111 7904 0.00110 47 65 FCULFA4 10807 0.00153 10807 0.00150 48 67 FCULFA2 3828 0.00054 4999 0.00070 49 68 FCULFA1 91 0.00001 141 0.00000 50 70 FCCM 175630 0.02471 183785 0.02560 51 74 TCVA 63243 0.00938 71554 0.01000 52 101 20_7 12741 0.00180 13152 0.00180 53 102 20_8 7799 0.00111 8039 0.00110 54 103 ATCSM 37734 0.00532 39636 0.00550 55 105 BMCSM 9920 0.00148 10798 0.00150 56 109 Boxcar 25253 0.00365 26157 0.00360 57 110 Buteo 16644 0.00243 17232 0.00240 58 115 ETCSZ 2852 0.00042 3102 0.00040 59 117 EastGreeley 25869 0.00377 26749 0.00370 60 119 GCCSMNW 8795 0.00131 9393 0.00130 61 120 Handley 15362 0.00230 15957 0.00220 62 123 M1 4928 0.00070 5110 0.00070 63 126 M2 14256 0.00202 14699 0.00200 64 128 M23 1081 0.00015 1406 0.00020 65 130 M25 648 0.00009 1005 0.00010 66 134 M29 4456 0.00062 4626 0.00060 67 135 M3 14203 0.00201 14701 0.00200 68 136 M30 5863 0.00082 6056 0.00080 69 137 M31 5747 0.00080 6853 0.00100 70 140 M4 9798 0.00138 10140 0.00140 71 141 M5N 2863 0.00041 2963 0.00040 72 142 M5S 8450 0.00118 8775 0.00120 73 144 M7 3521 0.00050 3615 0.00050 74 145 M9 5363 0.00077 6297 0.00090 75 147 NTMMSZ 26214 0.00375 27412 0.00380 76 150 Purse 13413 0.00198 14098 0.00200 77 151 RD1 5662 0.00080 5897 0.00080 78 152 RD10 4194 0.00060 4300 0.00060 79 153 RD2 2267 0.00032 3074 0.00040 80 154 RD4 2993 0.00042 3080 0.00040 81 155 RD5 3051 0.00043 3155 0.00040 82 156 RD6 1196 0.00017 1516 0.00020 83 157 RD7 8280 0.00117 8541 0.00120 84 158 RD9 4103 0.00058 4201 0.00060 85 161 RMCSM 33610 0.00476 35835 0.00500 86 163 RibbonCliff 7571 0.00111 7791 0.00110 87 165 SCNSZ 3002 0.00045 3368 0.00050 88 167 SCSZW 3938 0.00059 4332 0.00060 89 172 SPMSZ 7337 0.00106 7713 0.00110 90 173 SWA20CSM 9815 0.00139 10105 0.00140 91 176 WSCSZ 10753 0.00160 11060 0.00150 92 178 WestBoxcar 11642 0.00168 12065 0.00170 93 179 WestEstuary 680 0.00010 1197 0.00020 94 180 WestGreeley 32675 0.00475 34081 0.00470 95 181 WestPurse 6323 0.00090 6555 0.00090