Integra(ng an open source dynamic river model in hydrology modeling frameworks

Transcription

1 Integra(ng an open source dynamic river model in hydrology modeling frameworks Simula(on of Guadalupe and San Antonio River basin during a flood event with 1.3 x 10 5 computa(onal nodes at 100 m resolu(on. On a desktop computer this runs 300x faster than real (me. 1/4 next

2 Integra(ng an open source dynamic river model in hydrology modeling frameworks 2D Catchment- based runoff models must communicate with 1D river network models, but we lack consistent simple standards for model communica(ons. Integrated runoff/river models are limi(ng, because we cannot readily adopt new advances. prev. 2/4 next We have created the Simula(on Program for River Networks (SPRNT) as an open source dynamic river model (Saint- Venant equa(ons) that can communicate with a catchment model either through a simple 250 text- based hours Netlist or through an 274 Applica(on hours Programming Interface (API).

3 Integra(ng an open source dynamic river model in hydrology modeling frameworks SPRNT is unique in adap(ng methods developed in microprocessor electronic circuit design, known as Very Large Scale Integra(on (VLSI), to provide more rapid simula(on of river networks. We are just beginning to explore the myriad of possibili(es. prev. 3/4 next

4 Integra(ng an open source dynamic river model in hydrology modeling frameworks By clicking on the boxes below you can delve a li_le deeper into the world of SPRNT. prev. 4/4 Simula(on Program for River Networks (SPRNT) Electric circuits and river networks Netlist communica(on API communica(on Movie of model Publica(ons

5 Integra(ng an open source dynamic river model in hydrology Simula;on modeling Program frameworks for River Networks (SPRINT) Fully dynamic Saint- Venant equa(ons Large capacity, high performance, Full nonlinear opera(on with Newton- Raphson solu(on Standard Preismann implicit discre(za(on Jacobian bypass for solu(on accelera(on Opera(on using a Opera(on with Wri_en in C/C++ Netlist API Open Source available on for simple text- based input for integra(on with other models GitHub modular design (search on SPRNT)

6 Integra(ng an open source dynamic river model in hydrology modeling frameworks Netlist Similar to ideas used for semiconductor circuit Advantages of text- based format: Easy to read (for human and machine) Univ. of Texas at Aus.n sodware Easy to modify and create from different types of sodware Provides a simple standard for iden(fying Format is easy to translate from/to other models Disadvantages of text- based format: Data exchange through files, which can be large Forcing terms (boundary condi(ons) are defined in Netlist, so updates require rewri(ng file nodes and elements

7 Integra(ng an open source dynamic river model in hydrology Applica;on modeling frameworks Programming Interface A set of C++ func(ons that other programs can use to... Construct river network nodes and elements Call SPRNT to compute steady or unsteady solu(ons Obtain model results Advantages of API Works as run- stop- run- stop interface with hydrologic model SPRNT can be hosted on a computer accessed through web, i.e. does not require (ght coupling to hydrologic model. Disadvantages of API Requires C++ customiza(on to interface with another model Greater development effort to use

8 Integra(ng an open source dynamic river model in hydrology Open modeling source frameworks version control with Git Git is a sodware version control and management system Wri_en by Linus Torvalds (THE Linux guy) Allows easy collabora(ve code development git checkout, git add, git branch control modifica(on git clone, git push and git pull for mul(ple user copies Git is easy with GitHub clone John s git push, pull cloned repository on Jane s computer

9 Integra(ng an open source dynamic river model in hydrology modeling Git frameworks on the GitHub GitHub is a hos(ng service for Git repositories Free for open source code SPRNT is available on GitHub (just search GitHub for SPRNT) Clone SPRNT to get compile- ready code that can be readily updated with future bugfixes Fork SPRNT if you want to build something new using the code.

10 Integra(ng an open source dynamic river model in hydrology modeling SPRNT Modular frameworks Design command- line run_sprnt Netlist Parser File I/O API Model Construc;on Run Controls Steady Solve Unsteady Solve Nonlinear Solver Linear Solver

11 Integra(ng an open source dynamic river model in hydrology modeling Nodes frameworks and elements The con(nuous river must be broken into nodes and elements, connected at junc(ons. Q1 Q2 Communica(ng with a hydrological model requires transla(on of the 1D network of node- element- junc(on space to the 2D space of catchments. We use our Netlist format, but the community really needs to develop a standard! Q3

12 Integra(ng an open source dynamic river model in hydrology Electric modeling circuits frameworks and river networks Electric circuits are more complex, but the dynamic equa(ons with backwards and forwards propaga(ng characteris(cs are more difficult to solve. 3D connec(ons Connec;on network of a microprocessor electric circuit. h_p://visicad.eecs.umich.edu/bk/fgr Connec;on network of the Amazon River h_p://daac.ornl.gov/lba/guides/cd06_camre.html

13 Frank Liu Integra(ng an open source dynamic 3D circuit structure on a silicon chip river model IBM Research Aus.n R. Hodges in hydrology modeling frameworks Electric circuits and river networks Univ. oben f Texas at Aus.n Electric circuits are more complex, but the dynamic equa(ons with backwards and forwards propaga(ng characteris(cs are more difficult to solve. b hours 3D connec(ons 274 hours network of a microprocessor Connec(on network of the Amazon River Connec(on h_p://en.wikipedia.org/wiki/standard_cell electric circuit. h_p://daac.ornl.gov/lba/guides/cd06_camre.html This material is based upon work supported in part by the U.S. Na(onal Science Founda(on under Grant No. CCF

14 ## model of Possum Creek; done by John Doe, P.E. ## based on survey of YZ consul(ng, June 1999 def node id=1 sr= n=0.04 zr= hr=0.0 def xy x=0.0 y=6.0 x=2.5 y=1.0 x=3.5 y=1.0 x=6.0 y=6 def node id=2 sr= n=0.04 zr= hr=0.0 def trapezoidal bo_omwidth=1 slope= def segment up=1 down=2 length=39.1 def segment up=2 down=3 length=45.2 def qsource loca(on=1 (meunit=minute def (meseries (meunit=minute t=0 v=1.0 t=10 v= def op(ons metric=1 ## more comments can be added Netlist format Any lines star(ng with # are comments def pairs define nodes Complex cross- sec(on geometry defined by xy pairs Simple trapezoidal cross- sec(ons can be used segments: connec(ng nodes qsources: inflow terms op(ons: control Complete Netlist

15 ## automa(cally generated by a transla(on script on Jul :40:09 CDT ## ## ** one segment creek, test lateral inflow def op(ons metric=0 def op(ons (mestep=60 (mestepunit=second def op(ons prtq=1 prta=1 prtdepth=1 prtsurfelevz=1 def node id=1 sr= n=0.04 zr= hr=0.0 def trapezoidal bo_omwidth=0.1 slope=1.5 def node id=2 sr= n=0.04 zr= hr=0.0 def trapezoidal bo_omwidth=0.1 slope=1.50 def node id=3 sr= n=0.04 zr= hr=0.0 def trapezoidal bo_omwidth=0.1 slope=1.5 def node id=4 sr= n=0.04 zr= hr=0.0 def trapezoidal bo_omwidth=0.1 slope=1.5 Netlist format def node id=5 sr= n=0.04 zr= hr=0.0 def trapezoidal bo_omwidth=0.1 slope=1.5 def node id=6 sr= n=0.04 zr= hr=0.0 def trapezoidal bo_omwidth=0.1 slope=1.5 def segment up=1 down=2 length=40 def segment up=2 down=3 length=50 def segment up=3 down=4 length=60 def segment up=4 down=5 length=70 def segment up=5 down=6 length=80 def qsource loca(on=1 (meunit=minute def (meseries (meunit=minute t=0 v=1.0 t=10 v=1.0 t=1400 v=1.0 t=2500 v=1.0 t=5000 v=1.0 Con(nue

16 def lateralsource loca(on=2 (meunit=minute def (meseries (meunit=minute t=0 v=1.0 t=100 v=1.0 t=200 v=2.0 t=2000 v=1.0 Netlist format (con;nued) def op(ons stop(me=600 stop(meunit=minute pr(nterval=5 pr(ntervalunit=minute ## ** Done def lateralsource loca(on=3 (meunit=minute def (meseries (meunit=minute t=0 v=1.2 t=2000 v=1.2 def boundarycondi(on loca(on=6 type=area def (meseries (meunit=minute t=0 v=8 t=2800 v=8

17 #include <stdio.h> #include "base.h" #include "waveforms.h" #include "op(ons.h" #include "subcatch.h int main() { Subcatchment *BASIN;. header files for SPRNT a subcatchment object holds SPRNT data API model Used in place of Netlist to automate connec(ons to other models BASIN- >MakeNode(1, , 0.04, 0, 0, 0, 0, , 0, TRAP, 0.1,1.5); BASIN- >MakeNode(2, , 0.04, 0, 0, 0, 0, , 0, TRAP, 0.1,1.5); define nodes BASIN- >MakeStvEqua(on(4, 5, 80); BASIN- >MakeStvEqua(on(5, 6, 80); BASIN- >MakeQrSrcEqua(on(0, idx1 ); BASIN- >MakeQrSrcEqua(on(4, idx2); define segments define forcing terms int rc = BASIN- >TopologyCheck(); int rc2 = BASIN- >SteadySolve(120, 120, 5*OPT.Tol()); checking all nodes are connected steady solve : ini(al condi(ons Con(nue

18 // con(nued from previous page. start unsteady solve // the (me interval is 0 to 150 BASIN- >UnsteadySolve(150, 1, 25,.);... WV- >GetQatNode(6, &rep_val1[0]); WV- >GetZatNode(6, &rep_val2[0]); query results... // now we move forward to the (me point 420 sec (from 150 sec) BASIN- >UnsteadySolve(420, 1, 25,...); // assume we update forcing terms from 420 forward BASIN- >UpdateSource(idx1, 16, ); BASIN- >UpdateSource(idx2, 26, ); update forcing terms con(nue unsteady solve API model Used in place of Netlist to automate connec(ons to other models // simulate from 420 second to 680 sec BASIN- >UnsteadySolve(680, 1, 25, ); con(nue unsteady // repeat if (BASIN) delete BASIN; clean up return 0; }

19 Integra(ng an open source dynamic river model in hydrology modeling frameworks

20 Integra(ng an open source dynamic river model in hydrology modeling frameworks Hodges and Liu (2014) Rivers and Electric Networks: Crossing Disciplines in Modeling and Simula(on, Founda.ons and Trs in Electronic Design Automa.on 8:1:1-116 Liu and Hodges (2014) Applying microprocessor analysis methods to river network modeling, Environmental Modelling & SoBware, 52: Hodges (2013) Challenges in con(nental river dynamics, Environmental Modelling & SoBware, 50: Downloads at h_p:// Frank Liu IBM Research Aus.n Ben R. Hodges Univ. of Texas at Aus.n

21 Integra(ng an open source dynamic river model in hydrology modeling frameworks End of presenta(on