Numerical Simulation of Gas Injection in Vertical Water Saturated Porous Media

We present numerical simulations of drainage induced by air injection in a vertical water-saturated Hele-Shaw cell filled with glass microbeads. We use the macroscale Subsurface Transport Over Multiple Phases (STOMP) simulator developed by the Pacific Northwest National Laboratory’s Hydrology Group. To trigger fingering, we use random permeability fields consistent to capillary entry pressure fields. We compare the numerical results to our own experimental results shown in a previous study. We analyze the effects of the microheterogeneity degree as well as the macroscopic parameters on the gas saturation results. The main objective of the work is to investigate how microscopic effects could be accounted for by macroscopic variables during drainage.     

Saturated and Unsaturated Water Flow in Inclined Porous Media

This paper considers the two-dimensional saturated and unsaturated flow of water through inclined porous media, namely a waste dump or hill slope. Since the partial differential equation governing this water flow transforms from being parabolic to elliptic as the water flow varies from unsaturated to saturated, an iterative, finite differencing scheme is used to develop a numerical solution. The model can be used to investigate the effects that hill slope angle, depth of soil cover and hilltop width have on water accumulation in the dump and the time required for saturation to occur at different areas in the dump domain. The accuracy and reliability of the computer based solution is tested for two different boundary conditions – (1) no flow on all boundaries (i.e., the internal redistribution of soil moisture to steady state) and (2) a constant rainfall flux on the dump surface. Numerical studies then show the effects of changing the hill slope angle, depth of layer, and dump geometry on the flow characteristics in the dump.     

Humic Acid Transport in Water-Saturated Porous Media

The presence of humic acid (HA) in water poses environmental problems because it enhances the transport of several contaminants. A series of column experiments was conducted toward studying HA transport in different porous media under various pH, ionic strength, and flow rate conditions. The results showed that decreasing pH and increasing ionic strength increased adsorption and therefore delayed the transport of HA in porous media. However, increasing flow rate accelerated the transport of HA in porous media. The effects of pH, ionic strength, and flow rate varied with the solid matrix and were more evident in sands of smaller average particle diameter than in those with larger ones. These factors must be considered when predicting the environmental fate of HA. The results also suggested that HA adsorption was an important process controlling HA transport and should be considered in studies of HA behavior in porous media.     

Numerical Modeling on the Effect of Dissolved Oxygen on Nitrogen Transformation and Transport in Unsaturated Porous System

Nitrogen pollution in groundwater resulting from wastewater application to land is a common problem, and it causes a major threat to groundwater-based drinking water supplies. In this study, a numerical model is developed to study the nitrogen species transport and transformation in unsaturated porous media. Further, a new mass transfer module for dissolved oxygen (DO) is incorporated in the one-dimensional numerical model for nitrogen species transport to describe the fate and transport of nitrogen species, dissolved oxygen, dissolved organic carbon (DOC), and biomass. The spatial and temporal variation of dissolved oxygen is incorporated in the model through the mass transfer from gaseous phase to water phase in an unsaturated porous system. The numerical results of the water flow model and single species and multispecies transport model in an unsaturated zone developed for this purpose have been validated with the available analytical/numerical solution. The developed model is applied in clay loam, silt, and sand soils to analyze the transport behavior of nitrogen species under unsaturated condition. The numerical results suggest that the high rate of oxygen mass transfer from the air phase to the water phase positively increases the dissolved oxygen in the applied wastewater and enhances the nitrification process. Because of this high oxygen mass transfer, the nitrate nitrogen concentration significantly increases in the unsaturated zone and the same is transported to a larger depth at higher simulation period. On the other hand, the low rate of oxygen mass transfer implicitly enhances the denitrification process and finally reduces the nitrate nitrogen concentration in the unsaturated zone. The numerical results also show that the nitrate nitrogen transport is rapid in sandy soil when compared with clay loam and silty soils under high oxygen mass transfer rate. In essence, the high oxygen mass transfer rate significantly increases the nitrate nitrogen in the unsaturated zone, especially at a greater depth at larger time levels and eventually affects the groundwater quality.     

Performance Evaluation of Multiple Groundwater Flow and Nitrate Mass Transport Numerical Models

Environmental Modeling & Assessment - Benchmarking of different numerical models simulating groundwater flow and contaminant mass transport is the aim of the present study, in order to...     

Mathematical Modeling and Numerical Algorithms for Simulation of Oil Pollution

This paper deals with the mathematical modeling and algorithms for the problem of oil pollution. For solving this task, we derive the adjoint problem for the advection–diffusion equation describing the propagation of oil slick after an accident, which we call the main problem. We prove a fundamental equality between the solutions of the main and the adjoint problems. Based on this equality, we propose a novel method for the identification of the pollution source location and the accident time of oil emission. This approach is illustrated on an example for an accident in the offshore of the central part of the Vietnamese coast. Numerical simulations demonstrate the effectiveness of the proposed method. Besides, the method is verified for 1D model of substance propagation.     

Modelling the Competitive Sorption Process of Multiple Solutes During their Transport in Porous Media

A general mathematical model to solve the advection–dispersion transport equation for multiple solutes was developed, where the dual porosity mobile–immobile mass transfer, the two-site non-equilibrium model and first-order transformation reactions were included. The two-site model was expressed with an equilibrium sorption term and a kinetic term. One of three kinetic models could be used: the non-linear, the bilinear and the pore diffusion model. The traditional Freundlich or Langmuir isotherms were employed to simulate no-interaction between the solutes, but with the extended Freundlich or extended Langmuir isotherms, a competitive sorption could be simulated. The transport equation was solved with the Moving Concentration Slope method. The mathematical model was tested and further simplified by using real data from soil column experiments, with 1,2-cis-dichloroethene and trichloroethene as model contaminants and silica gel and real soil samples as porous medium.     

Response Surface Method Coupled with First-Order Reliability Method Based Methodology for Groundwater Flow and Contaminant Transport Model for the Uranium Tailings Pond Site

Groundwater flow and contaminant transport modelling is carried out to predict the concentrations of radionuclides in groundwater beneath the tailings pond. This is used as a tool for radiological impact assessment studies. The geo-hydrological parameters involved in the modelling exhibit inherent uncertainty associated with their values. This uncertainty may get propagated to the model output, i.e. the dose to the public. The propagation of the uncertainty in the input parameters to the output is modelled using suitable methods of probabilistic analysis. Response surface method coupled with first-order reliability method is used to develop a methodology for estimating the probability that the dose rate value through drinking water pathway at a location around the tailings pond exceeds the WHO guidelines for drinking water, termed as the probability of exceedance of acceptable levels. This method also gives the estimate of sensitivity of the probability of exceedance to the different input parameters. It is observed that the probability of failure decreases as the distance from tailings pond centre is increased and beyond a distance of around 0.5 km from tailings pond centre, the probability reduces to zero. The work also brings out the importance of quantifying the uncertainty in case of actual field problems where there is wide variation in the values of the various parameters within the domain under study.     

复杂岩溶条件下锰矿尾矿库地下水溶质运移特征数值模拟研究

受构造、地层等复杂地质情况综合影响的岩溶地区地下水,其溶质运移及流场模拟向来是水文地质、环境地质研究中的难点。以贵州省松桃县某锰矿尾矿库为例,在充分分析研究区周边水文地质、环境地质条件的基础上,建立地质模型,概化边界条件,利用有限元软件FEFLOW进行流场拟合。在此基础上,模拟特征污染物锰及氨氮的弥散情况。模拟结果表明,在不做任何防渗处理的情况下,特征污染物将在F201断层、地层产状等边界条件的影响下汇入地下暗河,最终流入地表水系及地下水系,流向东部的松江河。因此,必须对锰矿尾矿库采取防渗措施,从而降低其对地下水、地表水的污染风险。     

Characterization of Infiltration Rates from Landfills: Supporting Groundwater Modeling Efforts

The purpose of this paper is to review the literature to characterize infiltration rates from landfill liners to support groundwater modeling efforts. The focus of this investigation was on collecting studies that describe the performance of liners 'as installed' or 'as operated'. This document reviews the state of the science and practice on the infiltration rate through compacted clay liner (CCL) for 149 sites and geosynthetic clay liner (GCL) for 1 site. In addition, it reviews the leakage rate through geomembrane (GM) liners and composite liners for 259 sites. For compacted clay liners (CCL), there was limited information on infiltration rates (i.e., only 9 sites reported infiltration rates.), thus, it was difficult to develop a national distribution. The field hydraulic conductivities for natural clay liners range from 1 x 10(-9) cm s(-1) to 1 x 10(-4) cm s(-1), with an average of 6.5 x 10(-8) cm s(-1). There was limited information on geosynthetic clay liner. For composite lined and geomembrane systems, the leak detection system flow rates were utilized. The average monthly flow rate for composite liners ranged from 0-32 lphd for geomembrane and GCL systems to 0 to 1410 lphd for geomembrane and CCL systems. The increased infiltration for the geomembrane and CCL system may be attributed to consolidation water from the clay.     

某化工园区地下水污染迁移及控制数值模拟研究

某化工园区潜水中氟化物超标,为控制受污染地下水扩散,沿下游河堤修建防渗帷幕。通过分析研究区水文地质条件,应用FEFLOW 7.0建立地下水数值模型,预测防渗帷幕修建前后受污染地下水迁移情况,并评估防渗帷幕污染控制效果。结果表明：化工园区受污染地下水在以粉砂、粉细砂等岩性为主的潜水含水层中迁移速度快,第200天时污染物对源强附近及下游大片范围地下水造成污染;在污染源下游沿堤修建防渗帷幕短期虽可有效阻断污染物持续迁移,但会导致污染物在帷幕前端聚集,一段时间后地下水污染羽会绕过防渗帷幕发生渗流污染,须配合其他治理措施综合控制地下水污染。     

A Compartmental Model for Stable Time-Dependent Modeling of Surface Water and Groundwater

Water resource system planners make decisions that guide water management policy. The fundamental tools for assessing management and infrastructure strategies are combined hydro-economic models of river basins (RBHE models). These models have improved the economic efficiency of water use in situations of competition for scarce water resources. In RBHE models, a groundwater model is coupled with surface water models of the various water resources. Today, the groundwater models used in an RBHE model can be of two types: cell models or numerical models. Cell models are easy to use, but they are too simple to realistically describe the geology and hydrology of the area under investigation. Numerical models, in contrast, are closer to the physical behavior of the aquifer. However, the vast quantity of data to be analyzed makes them impractical for many management scenarios. Moreover, the calibrations of these high-resolution models are generally difficult and sensitive to the variation of parameters, especially when boundary conditions are dynamic. This is the case when dynamic river data or dynamic surface lake data are present. In this work, a compartmental cell model is built on the hydrogeology of the aquifer. In this approach, the hydrogeology of the aquifer and the dynamic boundary conditions are treated with separate models. A general mathematical formulation is presented where the calibration stage, the validation stage, and the prediction stage are formulated as a series of sub-model calibrations and solved using a general least squares routine. With this approach, it becomes possible to treat both the water level and the pumping rate in each cell as variables to be predicted. In most of the models, the pumping rates are known and the goals of the computation are to estimate the groundwater level. However, when for political or technical reasons access to some of the wells is difficult, the pumping rates are only partially known. Then, both groundwater levels and pumping rates are variables to be predicted by the groundwater model. A computer program was developed using MATLAB, with a Visual Basic graphical user interface using COM technology to access the advanced mathematical libraries. The approach was implemented with a real case study of the Yarkon–Taninim aquifer in Israel. The results indicate that the method is more stable than the classical approach.     

Selection and Calibration of Numerical Modeling in Flow and Water Quality

Numerical modeling is frequently used in coastal engineering research and application. One possible issue associated with using this method however, is that initial outcomes might differ from expectation. Modeling manipulation addresses this issue by changing the initial parameters which in turn affects the final output results.The advancement of modeling techniques in recent years has seen a gradual but steady accruement of knowledge regarding the relationship between input parameters and possible outputs, which has resulted in improved accuracy and efficiency methods. The advancements within this field have predominantly occurred on an individual basis and as such lack standardisation. An expert numerical modeler may use this knowledge subconsciously, yet may not know how to convey it to the model users. This paper has been written to introduce a systematic intelligent coding schema designed to organise current coastal engineering modeling manipulation knowledge into a standardised format. This system is relevant to the development of appropriate strategies for improved accuracy and efficiency as well as model modification to simulate specific real phenomena in prototype application cases.     

Selection and Calibration of Numerical Modeling in Flow and Water Quality

Numerical modeling is frequently used in coastal engineering research and application. One possible issue associated with using this method however, is that initial outcomes might differ from expectation. Modeling manipulation addresses this issue by changing the initial parameters which in turn affects the final output results. The advancement of modeling techniques in recent years has seen a gradual but steady accruement of knowledge regarding the relationship between input parameters and possible outputs, which has resulted in improved accuracy and efficiency methods. The advancements within this field have predominantly occurred on an individual basis and as such lack standardisation. An expert numerical modeler may use this knowledge subconsciously, yet may not know how to convey it to the model users. This paper has been written to introduce a systematic intelligent coding schema designed to organise current coastal engineering modeling manipulation knowledge into a standardised format. This system is relevant to the development of appropriate strategies for improved accuracy and efficiency as well as model modification to simulate specific real phenomena in prototype application cases.     

Modeling Solute Transport in Volcanic Ash Soils with Cation Exchange and Anion Retardation

Groundwater under agricultural lands is often contaminated by nitrate. In southern Japan, aquifers are covered by volcanic ash soils that can leach nitrate to the groundwater. In this study, column experiments using two volcanic ash soils (Kuroboku and Akahoya) were carried out. A mixed solution of KNO₃ and K₂SO₄ was used in the leaching experiments. Based on the experiments, a reactive transport model was developed using the Constrained Interpolation Profile method for ion transport and used to calculate chemical equilibrium for the cation exchange reactions. Anion adsorption was modeled by retardation in the numerical model. The developed simulation model results were compared to results obtained by the reactive transport model PHAST. The developed model was shown to quite well reproduce general characteristics of the experimental results. Also, the developed model results agreed well with results from PHAST. Slight discrepancy between observed and calculated breakthrough curves was probably caused by ignoring the kinetic reaction in the model calculations.     

Numerical Study of Hydrodynamic and Solute Transport with Discontinuous Flows in Coastal Water

This paper aimed to develop a depth-averaged explicit model for flow and pollutant transport in coastal waters based on the shallow water equations and the mass advection-diffusion equation. The proposed model was discretized using the finite volume method (FVM) with triangular cells. Then, it applied Roe’s approximate Riemann solution to compute the water momentum flux on the grid interfaces. This model enabled the higher accuracy in capturing the dry-wet moving fronts (discontinuous problems for flow and solute). The high-resolution scheme was evaluated to solve the advection and diffusion terms for mass transport. The model was verified by comparing the predictions of analytical solutions, laboratory tests, and other simulations for Gironde estuary with good computational accuracy. The developed model was also used to calculate the circulation and the motion of chemical oxygen demand (COD) pollutants from the sewage outfalls in the Zhuanghe coastal water with dry and wet moving boundaries. The research results showed that the residual current directions of spring and neap tides were basically the same in the Zhuanghe coastal water. However, the tide residual current of spring tide was slightly greater than that of the neap tide. In addition, there were tide residual currents from the northeast to the southwest in nearshore water and from the southwest to the northeast outside the banks, respectively. The tidal flows in the alongshore direction were strong, resulting in highly spread concentration distributions. In particular, the COD concentration reached some parts of the southern water. It could be seen that the excessive pollutant discharge from the sewage outfalls located at Zhuanghe district would cause serious pollution in aquaculture water near Shicheng island.     

A Numerical Model for the Transport of Chromium and Fine Sediments

A mathematical model which represents the transport processes of heavy metals and fine sediments in a fluvial stream was developed. The model consists of a three-equation system: the first one for total chromium concentration in the water column, C _Tw, the second one for total suspended sediment concentration, S _w, and the third one for chromium concentration in bed sediments, r. The third equation represents the chromium exchange between the water column and bed sediments by two processes: diffusion of soluble chromium and erosion/deposition of chromium sorbed to sediments. The basic assumption of the model is the instantaneous equilibrium. The main parameters are the partition coefficients in the water column and bed sediments, the depth of the active bed sediment layer, and the mass transfer coefficient between the water column and sediment pore water. The numerical model approximates the equations of advection–dispersion for chromium in water and suspended sediments by using a Eulerian third-order scheme. Numerical vs. analytical solutions were considered satisfactory for different initial, boundary, and sedimentological conditions. In order to estimate the impact of a chromium side discharge, the model was implemented for the Salado River in a reach of 65.6-km long (Santa Fe, Argentina). The results showed the effect of chromium discharge on almost the whole reach, then the vulnerability of the water quality in the Salado River when the flow was low was evidenced. When comparing the computed and measured results, the former showed a reasonable representation of the presence of chromium in water and bed sediments.     

Mathematical Modeling of Column and Field Dense Nonaqueous Phase Liquid Tracer Tests

Mathematical models for the simulation of dense nonaqueous phase liquid tracer tests (DTTs) in laboratory columns and in the field are developed and examined. The DTT technique is a means of estimating the quantity of dense nonaqueous phase liquid (DNAPL) in a domain of interest in an aquifer. The two-dimensional field DTT model uses the Method of Principal Directions and an asymmetrical upwind algorithm for describing advective transport. Both models include diffusion transport of tracer into and from low-permeability porous structures such as clay lenses, as well as the mass transport kinetics of partitioning tracer to and from the DNAPL droplets. The dependence of the effluent tracer concentration curves on the parameters of the models is explored, and conclusions are drawn regarding the applicability of, and several possible problems with, the DTT technique. Model results indicate that the DTT performs well at locating distributed droplets of DNAPL, but is unlikely to be useful in the assessment of pooled DNAPL.