User Forum of Software BASEMENT

Tamagnone Paolo

User

Registered: 2019-09-25

Posts: 1

Coupling baseplace_2D with basesub

Hello everyone,
is it possible couple baseplace_2D with basesub? I've seen that there is the possibility to introduce in basesub as boundary condition "coupling_baseplane" to use the depth of the hydrodynamic simulation deriving from baseplane_2D.
I'm trying to implement a rainfall-runoff model in basement.
I've already implemented the hydraulic part in baseplane_2D. I'm using precipitation as source.
Now I have problems to implement the basesub part and coupling it with baseplane_2D.
First of all:
- I'm using a triangular mesh in baseplane_2D; it cannot be used as geometry in basesub, isn't it? Does it require a cartesian grid mesh since in the description it's written (under geometry-zone-mesh) "filename of the .2dm gridfile"?;
- I didn't understand how to set the coupling part;

It would be very interesting if there were a chance to simulate infiltration coupling these 2 modules together, in order to send for each time step the depth calculated in baseplane_2D to basesub; then, once calculated the rate of seepage, this amount of water will be subtracted from the 2D domain in the following time step.   

this is my uncompleted code that provide me this error message
"""""""""""""""""""""""""""""""""""""""""""""
// --------------------------------------------------
// BASEMENT log file:
// Version:  2.8 R// bmc file: C:\Users\a\Desktop\dl_sfals\Simulazione/prova sup-subsurface.bmc
// started:  Fri Oct 11 16:19:05 2019
// --------------------------------------------------
/*

WARNING -> Warning at the definition of the boundary serie.
   The simulation time (t = '10800') is longer than time serie (tmax = '7680').
    The values of the time serie are assumed constant after the last time.

EXCEPTION -> soil_default != NULL, in file: ..\src\BASEsub\BASEsub.cpp on line 154,
   Version: 2.8 R5771, compiled at: May 31 2018, 15:29:20

"""""""""""""""""

the code:

PROJECT {
    title = RWHT
}
DOMAIN {
    PARALLEL {
        number_threads = 0
    }
    PHYSICAL_PROPERTIES {
        gravity   = 9.81
        viscosity = 1e-06
        rho_fluid = 1000
    }
    BASEPLANE_2D {
        region_name = dl_sfalsate
        GEOMETRY {
            type = 2dm
            file = dl sfals.2dm
            STRINGDEF {
                name               = outlet
                node_ids           = (78 79 80 87 89 88 68 54 46 34 33 32 31 21 13 7 2 )
                upstream_direction = right
            }
        }
        HYDRAULICS {
            PARAMETER {
                minimum_water_depth            = 0.001
                riemann_solver                 = exact
                riemann_tolerance              = 1.0e-6
                simulation_scheme              = exp
                dynamic_depth_solver           = on
                dynamic_depth_solver_precision = 0.001
            }
            FRICTION {
                type             = manning
                default_friction = 0.06
                wall_friction    = off
                input_type       = index_table
                index            = (01)
                friction         = (0.06)
            }
            INITIAL {
                type = dry
            }
            SOURCE {
                EXTERNAL_SOURCE {
                    name  = extreme_rainfall
                    type  = source_discharge
                    file  = Exteme rainfall.txt
                    index = (01)
                }
            }
            BOUNDARY {
                name        = outlet
                type        = hqrelation
                string_name = outlet
                slope       = 10
            }
        }
        TIMESTEP {
            total_run_time = 10800
            CFL            = 1
            start_time     = 0
        }
        OUTPUT {
            console_time_step = 60
            reference_time    = 01-09-2008-12-00-00
            restart_time_step = 1E32
            SPECIAL_OUTPUT {
                type             = BASEviz
                output_time_step = 60
            }
            SPECIAL_OUTPUT {
                type             = node_centered
                output_time_step = 60
                values           = (depth wse velocity)
                format           = sms
            }
            SPECIAL_OUTPUT {
                type             = stringdef_history
                output_time_step = 60
                stringdefs       = (outlet)
                stringdef_values = (integrate[Q] wse)
            }
            SPECIAL_OUTPUT {
                type             = balance
                balance_values   = (water_volume)
                output_time_step = 60
            }
            SPECIAL_OUTPUT {
                boundary_values  = (Q)
                type             = boundary_history
                output_time_step = 60
            }
        }
    }
    BASESUB {
        region_name = sub-surface
        GEOMETRY {
            ref_coordinates = (331076.711 1505479.27 13.4)
            base_vector_i   = (1 0 0    )
            base_vector_j   = (0 1 0)
            cell_size       = 0.1
            no_cells_i      = 80
            no_cells_j      = 80
            no_cells_k      = 10
            ZONE {
                mesh           = dl sfals.2dm
                mesh_direction = below
                type           = meshes
                name           = sotto-suolo
            }
        }
        TIMESTEP {
            time_step      = 2
            total_run_time = 10800
            start_time     = 0
        }
        PARAMETER {
            type      = theta
            cs_factor = 3.0
            theta_s   = 1.0
            theta_r   = 0.0
        }
        INITIAL {
            type  = hydrostatic
            value = 0
        }
        BOUNDARY {
            type           = coupling_baseplane
            name           = dl_sfalsate
            material_index = (01)
            depth_min      = 0.0001
        }
        MATERIAL {
            SOIL_DEF {
                kf   = 0.000000556
                zone = sotto-suolo
                RETENTION_MODEL {
                    type  = van_genuchten
                    alpha = 0.8
                    n     = 1.09
                    hs    = -0.001
                }
            }
        }
        OUTPUT {
            output_time_step  = 60
            console_time_step = 60
            SPECIAL_OUTPUT {
                type               = cut_plane
                format             = qgis
                output_time_step   = 60
                cut_plane          = yz
                cut_plane_position = 4
                values             = (theta hlb)
            }
        }
    }
    COUPLINGS {
        PARAMETER {
        }
        COUPLING {
            type               = external
            basement_subdomain = dl_sfalsate
            external_subdomain = sub-surface
            SUBSURFACE {
                basement_interface = sub-surface
                external_interface =
            }
        }
    }
}

Last edited by Tamagnone Paolo (2019-10-11 15:20:56)