Development and evaluation of a numerical model for steady state interface and/or free surface groundwater flow

2012 ◽  
Vol 434-435 ◽  
pp. 110-120 ◽  
Author(s):  
A. Aharmouch ◽  
B. Amaziane
Water ◽  
2018 ◽  
Vol 10 (12) ◽  
pp. 1787 ◽  
Author(s):  
Luka Malenica ◽  
Hrvoje Gotovac ◽  
Grgo Kamber ◽  
Srdjan Simunovic ◽  
Srikanth Allu ◽  
...  

A novel numerical model for groundwater flow in karst aquifers is presented. A discrete-continuum (hybrid) approach, in which a three-dimensional matrix flow is coupled with a one-dimensional conduit flow, was used. The laminar flow in the karst matrix is described by a variably saturated flow equation to account for important hydrodynamic effects in both the saturated and unsaturated zones. Turbulent conduit flow for both free surface and pressurized flow conditions was captured via the noninertia wave equation, whereas the coupling of two flow domains was established through an exchange term proportional to head differences. The novel numerical approach based on Fup basis functions and control-volume formulation enabled us to obtain smooth and locally conservative numerical solutions. Due to its similarity to the isogeometric analysis concept (IGA), we labeled it as control-volume isogeometric analysis (CV-IGA). Since realistic verification of the karst flow models is an extremely difficult task, the particular contribution of this work is the construction of a specially designed 3D physical model ( dimensions: 5.66 × 2.95 × 2.00 m) in order to verify the developed numerical model under controlled laboratory conditions. Heterogeneous porous material was used to simulate the karst matrix, and perforated pipes were used as karst conduits. The model was able to capture many flow characteristics, such as the interaction between the matrix and conduit, rainfall infiltration through the unsaturated zone, direct recharge through sinkholes, and both free surface and pressurized flow in conduits. Two different flow experiments are presented, and comparison with numerical results confirmed the validity of the developed karst flow model under complex laboratory conditions.


2018 ◽  
Vol 10 (2) ◽  
pp. 276-284 ◽  
Author(s):  
Gang Chen ◽  
Shiguang Xu ◽  
Chunxue Liu ◽  
Lei Lu ◽  
Liang Guo

Abstract Mine water inrush is one of the important factors threatening safe production in mines. The accurate understanding of the mine groundwater flow field can effectively reduce the hazards of mine water inrush. Numerical simulation is an important method to study the groundwater flow field. This paper numerically simulates the groundwater seepage field in the GaoSong ore field. In order to ensure the accuracy of the numerical model, the research team completed 3,724 field fissure measurements in the study area. The fracture measurement results were analyzed using the GEOFRAC method and the whole-area fracture network data were generated. On this basis, the rock mass permeability coefficient tensor of the aquifer in the study area was calculated. The tensor calculation results are used in the numerical model of groundwater flow. After calculation, the obtained numerical model can better represent the groundwater seepage field in the study area. In addition, we designed three different numerical models for calculation, mainly to explore the influence of the tensor assignment of permeability coefficient on the calculation results of water yield of the mine. The results showed that irrational fathom tensor assignment would cause a significant deviation in calculation results.


2021 ◽  
Vol 36 (3) ◽  
pp. 165-176
Author(s):  
Kirill Nikitin ◽  
Yuri Vassilevski ◽  
Ruslan Yanbarisov

Abstract This work presents a new approach to modelling of free surface non-Newtonian (viscoplastic or viscoelastic) fluid flows on dynamically adapted octree grids. The numerical model is based on the implicit formulation and the staggered location of governing variables. We verify our model by comparing simulations with experimental and numerical results known from the literature.


2003 ◽  
Vol 41 (1) ◽  
pp. 110-112
Author(s):  
ZhixiaN. Cao ◽  
Rodney Day ◽  
Sarah Liriano

1969 ◽  
Vol 36 (3) ◽  
pp. 505-515 ◽  
Author(s):  
D. C. Gakenheimer ◽  
J. Miklowitz

The propagation of transient waves in a homogeneous, isotropic, linearly elastic half space excited by a traveling normal point load is investigated. The load is suddenly applied and then it moves rectilinearly at a constant speed along the free surface. The displacements are derived for the interior of the half space and for all load speeds. Wave-front expansions are obtained from the exact solution, in addition to results pertaining to the steady-state displacement field. The limit case of zero load speed is considered, yielding new results for Lamb’s point load problem.


2021 ◽  
Author(s):  
Liqun Jiang ◽  
Ronglin Sun ◽  
Xing Liang

<p>Protection and management of groundwater resources demand high-resolution distributions of hydraulic parameters (e.g., hydraulic conductivity (K) and specific storage (Ss)) of aquifers. In the past, these parameters were obtained by traditional analytical solutions (e.g., Theis (1935) or Cooper and Jacob (1946)). However, traditional methods assume the aquifer to be homogeneous and yield the equivalent parameter, which are averages over a large volume and are insufficient for predicting groundwater flow and solute transport process (Butler & Liu, 1993). For obtaining the aquifer heterogeneity, some scholars have used kriging (e.g., Illman et al., 2010) and hydraulic tomography (HT) (e.g., Yeh & Liu, 2000; Zhu & Yeh, 2005) to describe the K distribution.</p><p>In this study, the laboratory heterogeneous aquifer sandbox is used to investigate the effect of different hydraulic parameter estimation methods on predicting groundwater flow and solute transport process. Conventional equivalent homogeneous model, kriging and HT are used to characterize the heterogeneity of sandbox aquifer. A number of the steady-state head data are collected from a series of single-hole pumping tests in the lab sandbox, and are then used to estimate the K fields of the sandbox aquifer by the steady-state inverse modeling in HT survey which was conducted using the SimSLE algorithm (Simultaneous SLE, Xiang et al., 2009), a built-in function of the software package of VSAFT2. The 40 K core samples from the sandbox aquifer are collected by the Darcy experiments, and are then used to obtain K fields through kriging which was conducted using the software package of Surfer 13. The role of prior information on improving HT survey is then discussed. The K estimates by different methods are used to predict the process of steady-state groundwater flow and solute transport, and evaluate the merits and demerits of different methods, investigate the effect of aquifer heterogeneity on groundwater flow and solute transport.</p><p>According to lab sandbox experiments results, we concluded that compared with kriging, HT can get higher precision to characterize the aquifer heterogeneity and predict the process of groundwater flow and solute transport. The 40 K fields from the K core samples are used as priori information of HT survey can promote the accuracy of K estimates. The conventional equivalent homogeneous model cannot accurately predict the process of groundwater flow and solute transport in heterogeneous aquifer. The enhancement of aquifer heterogeneity will lead to the enhancement of the spatial variability of tracer distribution and migration path, and the dominant channel directly determines the migration path and tracer distribution.</p>


Sign in / Sign up

Export Citation Format

Share Document