mipwelcona: mipwelcona: calculate the development of the concentrations...

In KeesVanImmerzeel/mipwelcona: Calculate The Development Of Concentrations Of Solvents In Pumped Groundwater

Description

The development of the concentrations of solvents in pumped groundwater is calculated based on:

• streamline calculations (backwards from well screens) in a MIPWA groundwater model;

• Initial concentrations in specified layers;

• Levels of concentration layers;

• Parameters of different processes.

Details

The mipwelcona package is a further development of the program WELCONA.

The following processes can be specified:

• Dispersion;

• Retardation;

• Decay.

The following functions are exported:

• mw_read_streamlines

• mw_read_well_filters

• mw_conc_init

• mw_conc_streamlines

• mw_create_sl_fltr_table

• mw_conc_filters

• mw_conc_wells

Sample data sets are made available by calling the appropiate functions:

• mw_example_concentrations

• mw_example_conc_layer_levels

• mw_example_base_streamline_conc_table

• mw_example_conc_streamlines

• mw_example_sl_fltr_table

Example work flow:

• Step 1: read streamlines.

         fname <- system.file("extdata","streamlines.iff",package="mipwelcona")
         strm_lns <- mw_read_streamlines(fname)

• Step 2: Read well filters.

         fname <- system.file("extdata","well_filters.ipf",package="mipwelcona")
         well_fltrs <- mw_read_well_filters(fname)

• Step 3: Create concentration layer levels (8 example rasters layers in this case).

         conc_l_lev <- mw_example_conc_layer_levels()

• Step 4: Create Initial concentrations of different layers in the subsoil (9 example rasters layers in this case).

         conc_l <- mw_example_concentrations()

• Step 5: Initialize base streamline concentration table.

         conc_strm_lns <- mw_conc_init(strm_lns, conc_l_lev, conc_l)

• Step 6: Calculate concentrations on streamlines at specified times.

         conc_streamlines <- mw_conc_streamlines(conc_strm_lns, times=c(1*365,5*365,10*365,25*365), processes=c("dispersion","decay", "retardation"), alpha=0.3, rho=3, labda=0.0001, retard=1)

• Step 7: Create table with indices linking the streamlines to well filters.

         sl_fltr_table <- mw_create_sl_fltr_table(strm_lns, well_fltrs, maxdist=100)

• Step 8: Calculate concentrations (mixed) of selected filters.

         # Overrule previously created table first for this example. Normally, you wouln't do this!
         sl_fltr_table <- mw_example_sl_fltr_table()
         conc <- mw_conc_filters( fltr_nrs=c(1,2), well_fltrs, sl_fltr_table, conc_streamlines )

KeesVanImmerzeel/mipwelcona documentation built on Nov. 11, 2020, 10:35 p.m.