Theoretical analysis of regional groundwater flow: 1. Analytical and numerical solutions to the mathematical model

1966 ◽  
Vol 2 (4) ◽  
pp. 641-656 ◽  
Author(s):  
R. Allan Freeze ◽  
P. A. Witherspoon
2004 ◽  
Vol 36 (4) ◽  
pp. 2021
Author(s):  
A. Παναγόπουλος ◽  
E. Δρακοπούλου ◽  
V. Περλέρος

MODFLOW is a very well verified code of mathematical modeling for simulation of saturated groundwater flow in porous medium. Groundwater flow simulation in discontinuity media (i.e. media characterized by dominance of secondary and tertiary porosity as opposed to primary porosity), such as karstic aquifers, utilizing specialized models is problematic. Due to existing impedes the use of the conventional model MODFLOW was attempted for the simulation of the karstic system of the Viotikos Kifisos river aiming predominantly at assessing the potential, restrictions, particularities and conditions under which such a modelling code could be implemented, especially when relatively restricted volume of raw data is available. Compilation and calibration of the model suggest that MODFLOW may in general be implemented and can provide useful results. As in every mathematical model, knowledge of the assumptions made and the intrinsic restrictions involved is required, along with in-depth evaluation of its capabilities. The mathematical model of Viotikos Kifisos basin forms a valuable tool for management of its water resources and study of regional groundwater evolution.


2014 ◽  
Vol 986-987 ◽  
pp. 1418-1421
Author(s):  
Jun Shan Li

In this paper, we propose a meshless method for solving the mathematical model concerning the leakage problem when the pressure is tested in the gas pipeline. The method of radial basis function (RBF) can be used for solving partial differential equation by writing the solution in the form of linear combination of radius basis functions, that is, when integrating the definite conditions, one can find the combination coefficients and then the numerical solution. The leak problem is a kind of inverse problem that is focused by many engineers or mathematical researchers. The strength of the leak can find easily by the additional conditions and the numerical solutions.


Author(s):  
Quyang Ma ◽  
Guoan Yang ◽  
Mengjun Li

An elbow-shaped surge tank is proposed to suppress the pressure pulsations. The transfer matrix method was developed and the mathematical model was established to predict the distribution of pressure pulsations in the piping system (on which a surge tank was already installed) with an elbow-shaped surge tank. Simulation work of the whole piping system was performed. The results show that the elbow-shaped surge tank has good performance to attenuate the pressure pulsations. The frequency analysis shows that the amplitude for the first pulsation frequency is attenuated to a low level. The impulse response was analyzed to examine the efficiency of suppressing pulsations by using the suppressor. The theoretical analysis showed that there exists the optimal suppression performance when setting the distance between the elbow-shaped surge tank and the existing one. Meanwhile, modifying the ratio of length to diameter with a fixed surge volume could also impact the pressure pulsations. The analysis results can be used as a reference in designing and installing the elbow-shaped surge tank.


2021 ◽  
Author(s):  
Ping-Cheng Hsieh ◽  
Tzu-Ting Huang

Abstract. This study discussed water storage in aquifers of hillslopes under temporally varied rainfall recharge by employing a hillslope-storage equation to simulate groundwater flow. The hillslope width was assumed to vary exponentially to denote the following complex hillslope types: uniform, convergent, and divergent. Both analytical and numerical solutions were acquired for the storage equation with a recharge source. The analytical solution was obtained using an integral transform technique. The numerical solution was obtained using a finite difference method in which the upwind scheme was used for space derivatives and the third-order Runge–Kutta scheme was used for time discretization. The results revealed that hillslope type significantly influences the drains of hillslope storage. Drainage was the fastest for divergent hillslopes and the slowest for convergent hillslopes. The results obtained from analytical solutions require the tuning of a fitting parameter to better describe the groundwater flow. However, a gap existed between the analytical and numerical solutions under the same scenario owing to the different versions of the hillslope-storage equation. The study findings implied that numerical solutions are superior to analytical solutions for the nonlinear hillslope-storage equation, whereas the analytical solutions are better for the linearized hillslope-storage equation. The findings thus can benefit research on and have application in soil and water conservation.


Mathematics ◽  
2020 ◽  
Vol 8 (11) ◽  
pp. 1961
Author(s):  
Kanognudge Wuttanachamsri

Cilia on the surface of ciliated cells in the respiratory system are organelles that beat forward and backward to generate metachronal waves to propel mucus out of lungs. The layer that contains the cilia, coating the interior epithelial surface of the bronchi and bronchiolesis, is called the periciliary layer (PCL). With fluid nourishment, cilia can move efficiently. The fluid in this region is named the PCL fluid and is considered to be an incompressible, viscous, Newtonian fluid. We propose there to be a free boundary at the tips of cilia underlining a gas phase while the cilia are moving forward. The Brinkman equation on a macroscopic scale, in which bundles of cilia are considered rather than individuals, with the Stefan condition was used in the PCL to determine the velocity of the PCL fluid and the height/shape of the free boundary. Regarding the numerical methods, the boundary immobilization technique was applied to immobilize the moving boundaries using coordinate transformation (working with a fixed domain). A finite element method was employed to discretize the mathematical model and a finite difference approach was applied to the Stefan problem to determine the free interface. In this study, an effective stroke is assumed to start when the cilia make a 140∘ angle to the horizontal plane and the velocitiesof cilia increase until the cilia are perpendicular to the horizontal plane. Then, the velocities of the cilia decrease until the cilia make a 40∘ angle with the horizontal plane. From the numerical results, we can see that although the velocities of the cilia increase and then decrease, the free interface at the tips of the cilia continues increasing for the full forward phase. The numerical results are verified and compared with an exact solution and experimental data from the literature. Regarding the fixed boundary, the numerical results converge to the exact solution. Regarding the free interface, the numerical solutions were compared with the average height of the PCL in non-cystic fibrosis (CF) human tissues and were in excellent agreement. This research also proposes possible values of parameters in the mathematical model in order to determine the free interface. Applications of these fluid flows include animal hair, fibers and filter pads, and rice fields.


2015 ◽  
Vol 2015 ◽  
pp. 1-8 ◽  
Author(s):  
Yong Song ◽  
Yibin Li ◽  
Caihong Li ◽  
Xin Ma

This paper presents a mathematical model of multirobot cooperative hunting behavior. Multiple robots try to search for and surround a prey. When a robot detects a prey it forms a following team. When another “searching” robot detects the same prey, the robots form a new following team. Until four robots have detected the same prey, the prey disappears from the simulation and the robots return to searching for other prey. If a following team fails to be joined by another robot within a certain time limit the team is disbanded and the robots return to searching state. The mathematical model is formulated by a set of rate equations. The evolution of robot collective hunting behaviors represents the transition between different states of robots. The complex collective hunting behavior emerges through local interaction. The paper presents numerical solutions to normalized versions of the model equations and provides both a steady state and a collaboration ratio analysis. The value of the delay time is shown through mathematical modeling to be a strong factor in the performance of the system as well as the relative numbers of the searching robots and the prey.


2013 ◽  
Vol 390 ◽  
pp. 242-245 ◽  
Author(s):  
Alexander V. Chekanin

The article deals with the actual problem of improving the accuracy of determining the dynamic characteristics of beam structures. To solve such problems the displacement method is used. Defining matrices are calculated with the Godunovs scheme. Numerical solutions in this case can be obtained practically with any accuracy within accepted hypotheses of the mathematical model of the calculated object. This suggests that the resulting solutions are standard. The examples of determining the natural frequencies of vibrations of beam structures that demonstrate an extremely high accuracy of the proposed algorithm are given.


2010 ◽  
Vol 51 ◽  
Author(s):  
Donatas Švitra ◽  
Giedrius Žemaitis

This article describes the authors’ work the existing ecosystem in the Curonian Lagoon. Using the mathematical model of the ichtyocenose (1)–(2) is simulated the dynamics of the ichtyocenose in the Curonian Lagoon. It is done by using Runge–Kut IV method from this simulation program ModelMaker. The model numerical solutions of F1 - F8 are compared with the experimental data for the monitoring of fish. The dynamics is projected to the year 2016.


Sign in / Sign up

Export Citation Format

Share Document