430 likes | 543 Views
Using WPHAST to Model the Central Oklahoma Aquifer—Flow. David Parkhurst, Ken Kipp, and Scott Charlton. Arsenic in the Central Oklahoma Aquifer. Arsenic mostly in confined part of aquifer Arsenic associated with high pH Flow: Unconfined Confined Unconfined. The Wizard.
E N D
Using WPHAST to Model the Central Oklahoma Aquifer—Flow David Parkhurst, Ken Kipp, and Scott Charlton
Arsenic in the Central Oklahoma Aquifer • Arsenic mostly in confined part of aquifer • Arsenic associated with high pH • Flow: • Unconfined • Confined • Unconfined
The Wizard • Procures the minimum amount of information for a simulation • All data can be changed subsequently • EXCEPT location of coordinate origin on map and angle, when using a map bitmap • Heterogeneity can be defined subsequently
years Select Units / year
Importing a Bitmap • Import map c:\student\ok_wphast\ok.bmp • World file (ARC/INFO generated) ok.bpw
X 90000, 16 nodes Y 48000, 9 nodes Z 400, 5 nodes Move grid as shown river at north and south boundaries (east side) east boundary at contact Define Grid Mouse Controls • Inside grid affects grid • Outside grid affects all • Middle mouse—pan • Right mouse—zoom in/out • Cntl left mouse—2D rotate
Leave Default for Chemistry Initial Conditions Chemistry ignored for Flow Only
Define Time Step and Simulation Period Steady-flow calculation has automatic time stepping Default minimum = time step/1000; maximum = time step*1000
WPHAST Main Screen Save in c:\student\okflow Save often File|Save As… Tree • Current definitions • Expand/hide branches • Check to toggle visibility • Click to select • Double click to modify
WPHAST Main Screen • View|View From Direction • Cntl-d, page up/down • Home • Edit|Undo • Cntl-z • Edit|Redo • Cntl-y Mouse Controls • Left mouse—3D rotate • Middle mouse—pan • Right mouse—zoom in/out • Cntl left mouse—2D rotate Tools|Geometry Scale the Z-axis (40x)
Boundary Conditions • Defined by “cell” or “node” • We will define 6 boundary conditions • 1 Flux b.c. for recharge • 2 Rivers • 2 Leaky b.c. • 1 Specified head b.c. • Last definition for a cell takes precedence
Flux Boundary Condition Create new zone • 2D view Home key • Recharge (flux) on Eastern 2/3 of area encloses X=36000 to 90000 • Depth and location of zone can be adjusted later
Flux Boundary Condition Choose Flux_BC
Flux Boundary Condition Modify zone • Z definition encloses a single plane of nodes at 400 m
Define the Little River • View • View|Hide all • View|Set projection to|Parallel • Tools|New river or icon • Draw river ~(0,46000) to ~(0,90000) • Type “enter” after last point
Leaky Boundaries • North • X to edge of confining unit 0-35000 • Y 48000-48000 • Z 0-400 • South • X to edge of confining unit 0-42000 • Y 0-0 • Z 0-400
Specified Head: Lake Stanley Draper • Zone must enclose a node • Head 348 • Z 300-400 Click Apply
Run Run dialogue Run simulation Save
Model Viewer • 3D Viewer • Views heads, velocities, and user-specified chemical data • Slices and dices • Animations • Same mouse control as WPhast • Reads HDF file made by PHAST (.h5)
Model Viewer • File|New • 3 screens to get started
View Head Distribution with Model Viewer • Show|Solid • Show|Color Bar • Tools|Geometry • Set Z to 40 (40x vertical exaggeration)
Mouse Controls • Left mouse— 3D rotate • Middle mouse—pan • Right mouse—zoom in/out • Cntl left mouse—2D rotate
Model Viewer—“Show” menu Select items to be present in the visualization • Visualization • Solid/ • Isosurface • None • Vectors • Model features • Grid lines • Time • Color bar
Model Viewer—“Tools” menu Provides menus that change the look of the features selected by Show. • Data • Color bar • Geometry • Model features • Crop • Animation
Summary • PHAST as a flow model • WPhast and GoPhast are 3D interfaces to PHAST • 3D viewing by Model Viewer Optional exercises follow Next—Transport and Reaction
Optional Exercises • Add a well • Refine grid • Add another river
Adding a Well • XY location • Diameter/Radius • Flow allocation • Open intervals • Pumping/Injection schedule • Positive is injection • Negative is pumping