共查询到20条相似文献，搜索用时 46 毫秒
1.
《Advances in water resources》2002,25(1):67-84
A discontinuous Galerkin (DG) finite element method is described for the two-dimensional, depth-integrated shallow water equations (SWEs). This method is based on formulating the SWEs as a system of conservation laws, or advection–diffusion equations. A weak formulation is obtained by integrating the equations over a single element, and approximating the unknowns by piecewise, possibly discontinuous, polynomials. Because of its local nature, the DG method easily allows for varying the polynomial order of approximation. It is also “locally conservative”, and incorporates upwinded numerical fluxes for modeling problems with high flow gradients. Numerical results are presented for several test cases, including supercritical flow, river inflow and standard tidal flow in complex domains, and a contaminant transport scenario where we have coupled the shallow water flow equations with a transport equation for a chemical species. 相似文献
2.
Wave breaking and wave runup/rundown have a major influence on nearshore hydrodynamics, morphodynamics and beach evolution. In the case of wave breaking, there is significant mixing of air and water at the wave crest, along with relatively high kinetic energy, so prediction of the free surface is complicated. Most hydrodynamic studies of surf and swash zone are derived from single-phase flow, in which the role of air is ignored. Two-phase flow modeling, consisting of both phases of water and air, may be a good alternative numerical modeling approach for simulating nearshore hydrodynamics and, consequently, sediment transport. A two-phase flow tool can compute more realistically the shape of the free surface, while the effects of air are accounted for. This paper used models based on two-dimensional, two-phase Reynolds-averaged Navier–Stokes equations, the volume-of-fluid surface capturing technique and different turbulence closure models, i.e., k–ε, k–ω and re-normalized group (RNG). Our numerical results were compared with the available experimental data. Comparison of the employed method with a model not utilizing a two-phase flow modeling demonstrates that including the air phase leads to improvement in simulation of wave characteristics, especially in the vicinity of the breaking point. The numerical results revealed that the RNG turbulence model yielded better predictions of nearshore zone hydrodynamics, although the k–ε model also gave satisfactory predictions. The model provides new insights for the wave, turbulence and means flow structure in the surf and swash zones. 相似文献
3.
Clint Dawson Ethan J. KubatkoChristopher Mirabito Craig MichoskiNishant Panda 《Advances in water resources》2011,34(9):1165-1176
Storm surge due to hurricanes and tropical storms can result in significant loss of life, property damage, and long-term damage to coastal ecosystems and landscapes. Computer modeling of storm surge can be used for two primary purposes: forecasting of surge as storms approach land for emergency planning and evacuation of coastal populations, and hindcasting of storms for determining risk, development of mitigation strategies, coastal restoration and sustainability.Storm surge is modeled using the shallow water equations, coupled with wind forcing and in some events, models of wave energy. In this paper, we will describe a depth-averaged (2D) model of circulation in spherical coordinates. Tides, riverine forcing, atmospheric pressure, bottom friction, the Coriolis effect and wind stress are all important for characterizing the inundation due to surge. The problem is inherently multi-scale, both in space and time. To model these problems accurately requires significant investments in acquiring high-fidelity input (bathymetry, bottom friction characteristics, land cover data, river flow rates, levees, raised roads and railways, etc.), accurate discretization of the computational domain using unstructured finite element meshes, and numerical methods capable of capturing highly advective flows, wetting and drying, and multi-scale features of the solution.The discontinuous Galerkin (DG) method appears to allow for many of the features necessary to accurately capture storm surge physics. The DG method was developed for modeling shocks and advection-dominated flows on unstructured finite element meshes. It easily allows for adaptivity in both mesh (h) and polynomial order (p) for capturing multi-scale spatial events. Mass conservative wetting and drying algorithms can be formulated within the DG method.In this paper, we will describe the application of the DG method to hurricane storm surge. We discuss the general formulation, and new features which have been added to the model to better capture surge in complex coastal environments. These features include modifications to the method to handle spherical coordinates and maintain still flows, improvements in the stability post-processing (i.e. slope-limiting), and the modeling of internal barriers for capturing overtopping of levees and other structures. We will focus on applications of the model to recent events in the Gulf of Mexico, including Hurricane Ike. 相似文献
4.
Applying Quantitative Molecular Tools for Virus Transport Studies: Opportunities and Challenges 
下载免费PDF全文
下载免费PDF全文 Bacteriophages have been used in soil column studies for the last several decades as surrogates to study the fate and transport behavior of enteric viruses in groundwater. However, recent studies have shown that the transport behavior of bacteriophages and enteric viruses in porous media can be very different. The next generation of virus transport science must therefore provide more data on mobility of enteric viruses and the relationship between transport behaviors of enteric viruses and bacteriophages. To achieve this new paradigm, labor intensity devoted to enteric virus quantification method must be reduced. Recent studies applied quantitative polymerase chain reaction (qPCR) to column filtration experiments to study the transport behavior of human adenovirus (HAdV) in porous media under a variety of conditions. A similar approach can be used to study the transport of other enteric viruses such as norovirus. Analyzing the column samples with both qPCR and culture assays and applying multiplex qPCR to study cotransport behavior of more than one virus will provide information to under‐explored areas in virus transport science. Both nucleic acid extraction kits and one‐step lysis protocols have been used in these column studies to extract viral nucleic acid for qPCR quantification. The pros and cons of both methods are compared herein and solutions for overcoming problems are suggested. As better understanding of the transport behavior of enteric viruses is clearly needed, we strongly advocate for application of rapid molecular tools in future studies as well as optimization of protocols to overcome their current limitations. 相似文献
5.
《Advances in water resources》2005,28(7):729-744
In this paper, we discuss the local discontinuous Galerkin (LDG) method applied to elliptic flow problems and give details on its implementation, focusing specifically on the case of piecewise linear approximating functions. The LDG method is one a family of discontinuous Galerkin (DG) methods proposed for diffusion models. These DG methods allow for very general hp finite element meshes, and produce locally conservative fluxes which can be used in coupling flow with transport. The drawback to DG methods, when compared to their continuous counterparts, is the number of degrees of freedom required to compute the solution. This motivates a coupled approach, discussed herein, where the solution is allowed to be continuous or discontinuous on a node-by-node basis. This coupled approximation is locally conservative in regions where the numerical solution is discontinuous. Numerical results for fully discontinuous, continuous and coupled discontinuous/continuous solutions are given, where we compare solution accuracy, matrix condition numbers and mass balance errors for the various approaches. 相似文献
6.
We present an extended law of mass-action (xLMA) method for multiphase equilibrium calculations and apply it in the context of reactive transport modeling. This extended LMA formulation differs from its conventional counterpart in that (i) it is directly derived from the Gibbs energy minimization (GEM) problem (i.e., the fundamental problem that describes the state of equilibrium of a chemical system under constant temperature and pressure); and (ii) it extends the conventional mass-action equations with Lagrange multipliers from the Gibbs energy minimization problem, which can be interpreted as stability indices of the chemical species. Accounting for these multipliers enables the method to determine all stable phases without presuming their types (e.g., aqueous, gaseous) or their presence in the equilibrium state. Therefore, the here proposed xLMA method inherits traits of Gibbs energy minimization algorithms that allow it to naturally detect the phases present in equilibrium, which can be single-component phases (e.g., pure solids or liquids) or non-ideal multi-component phases (e.g., aqueous, melts, gaseous, solid solutions, adsorption, or ion exchange). Moreover, our xLMA method requires no technique that tentatively adds or removes reactions based on phase stability indices (e.g., saturation indices for minerals), since the extended mass-action equations are valid even when their corresponding reactions involve unstable species. We successfully apply the proposed method to a reactive transport modeling problem in which we use PHREEQC and GEMS as alternative backends for the calculation of thermodynamic properties such as equilibrium constants of reactions, standard chemical potentials of species, and activity coefficients. Our tests show that our algorithm is efficient and robust for demanding applications, such as reactive transport modeling, where it converges within 1–3 iterations in most cases. The proposed xLMA method is implemented in Reaktoro, a unified open-source framework for modeling chemically reactive systems. 相似文献
7.
间断有限元（Discontinuous Galerkin：DG）方法具有低数值频散、网格剖分灵活、能模拟地震波在复杂介质中传播等优点.因此,本文将一种新的DG方法推广到双相和黏弹性等复杂介质的地震波场模拟,发展了求解Biot弹性波方程和D'Alembert介质波动方程的DG方法.首先通过引入辅助变量将Biot双相介质弹性波方程和D'Alembert介质波动方程转化为关于时间-空间的一阶偏微分方程组,然后对该方程组进行DG空间离散,得到半离散化的常微分方程组.最后,对此常微分方程组,应用加权的Runge-Kutta格式进行时间推进计算.数值结果表明,DG方法可以有效地求解Biot双相介质弹性波方程和D'Alembert介质波动方程,并能很好地压制因离散求解波动方程而产生的数值频散,获得清晰的各种地震波震相. 相似文献
8.
Moment equation methods are popular and powerful tools for modeling transport processes in randomly heterogeneous porous media, but the application of these methods to advection-dispersion equations often leads to erroneous oscillations. Perturbative methods, required to close systems of moment equations, become inaccurate for large perturbations; however, little quantitative theory exists for determining when this occurs for advection-dispersion equations. We consider three different methods (asymptotic approximation, Eulerian truncation, and iterative solution) for closing and solving advection-dispersion moment equations describing transport in stratified porous media with random permeability. We obtain approximate analytical expressions for time above which the asymptotic approximation to the mean diverges, in particular quantifying the impact that dispersion has on delaying—but not eliminating—divergence. We demonstrate that Eulerian truncation and iterative solution methods do not eliminate divergent behavior either. Our divergence criteria provide a priori estimates that signal a warning to the practitioner of stochastic advection-dispersion equations to carefully consider whether to apply perturbative approaches. 相似文献
9.
1 INTRODUCTION In recent years, due to the increase in population and industrial developments, mankind has faced manyproblems associated with rivers, coastal waters and reservoirs. Some of these problems are flood control,water supply, power generation, and irrigation. In addition, making new hydraulic structures changesnatural conditions. Prediction of these changes is necessary for designing such constructions. For solutionof these problems usually an assessment of flow pattern, sedim… 相似文献
10.
We present advances in compositional modeling of two-phase multi-component flow through highly complex porous media. Higher-order methods are used to approximate both mass transport and the velocity and pressure fields. We employ the Mixed Hybrid Finite Element (MHFE) method to simultaneously solve, to the same order, the pressure equation and Darcy's law for the velocity. The species balance equation is approximated by the discontinuous Galerkin (DG) approach, combined with a slope limiter. In this work we present an improved DG scheme where phase splitting is analyzed at all element vertices in the two-phase regions, rather than only as element averages. This approximation is higher-order than the commonly employed finite volume method and earlier DG approximations. The method reduces numerical dispersion, allowing for an accurate capture of shock fronts and lower dependence on mesh quality and orientation. Further new features are the extension to unstructured grids and support for arbitrary permeability tensors (allowing for both scalar heterogeneity, and shear anisotropy). The most important advancement in this work is the self-consistent modeling of two-phase multi-component Fickian diffusion. We present several numerical examples to illustrate the powerful features of our combined MHFE–dg method with respect to lower-order calculations, ranging from simple two component fluids to more challenging real problems regarding CO2 injection into a vertical domain saturated with a multi-component petroleum fluid. 相似文献
11.
《Advances in water resources》1998,22(1):33-57
In this article we consider the transport of an adsorbing solute in a two-region model of a chemically and mechanically heterogeneous porous medium when the condition of large-scale mechanical equilibrium is valid. Under these circumstances, a one-equation model can be used to predict the large-scale averaged velocity, but a two-equation model may be required to predict the regional velocities that are needed to accurately describe the solute transport process. If the condition of large-scale mass equilibrium is valid, the solute transport process can be represented in terms of a one-equation model and the analysis is simplified greatly. The constraints associated with the condition of large-scale mass equilibrium are developed, and when these constraints are satisfied the mass transport process can be described in terms of the large-scale average velocity, an average adsorption isotherm, and a single large-scale dispersion tensor. When the condition of large-scale mass equilibrium is not valid, two equations are required to describe the mass transfer process, and these two equations contain two adsorption isotherms, two dispersion tensors, and an exchange coefficient. The extension of the analysis to multi-region models is straight forward but tedious. 相似文献
12.
Upscaling pore-scale processes into macroscopic quantities such as hydrodynamic dispersion is still not a straightforward matter for porous media with complex pore space geometries. Recently it has become possible to obtain very realistic 3D geometries for the pore system of real rocks using either numerical reconstruction or micro-CT measurements. In this work, we present a finite element–finite volume simulation method for modeling single-phase fluid flow and solute transport in experimentally obtained 3D pore geometries. Algebraic multigrid techniques and parallelization allow us to solve the Stokes and advection–diffusion equations on large meshes with several millions of elements. We apply this method in a proof-of-concept study of a digitized Fontainebleau sandstone sample. We use the calculated velocity to simulate pore-scale solute transport and diffusion. From this, we are able to calculate the a priori emergent macroscopic hydrodynamic dispersion coefficient of the porous medium for a given molecular diffusion Dm of the solute species. By performing this calculation at a range of flow rates, we can correctly predict all of the observed flow regimes from diffusion dominated to convection dominated. 相似文献
13.
Recent work with stochastic inverse modeling techniques has led to the development of efficient algorithms for the construction of transmissivity (T) fields conditioned to measurements of T and head. Small numbers of calibration targets and correlation between model parameters in these inverse solutions can lead to a relatively large region in parameter space that will produce a near optimal calibration of the T field to measured heads. Most applications of these inverse techniques have not considered the effects of non-unique calibration on subsequent predictions made with the T fields. Use of these T fields in predictive contaminant transport modeling must take into account the non-uniqueness of the T field calibration. A recently developed ‘predictive estimation’ technique is presented and employed to create T fields that are conditioned to observed heads and measured T values while maximizing the conservatism of the associated predicted advective travel time. Predictive estimation employs confidence and prediction intervals calculated simultaneously on the flow and transport models, respectively. In an example problem, the distribution of advective transport results created with the predictive estimation technique is compared to the distribution of results created under traditional T field optimization where model non-uniqueness is not considered. The predictive estimation technique produces results with significantly shorter travel times relative to traditional techniques while maintaining near optimal calibration. Additionally, predictive estimation produces more accurate estimates of the fastest travel times. 相似文献
14.
《国际泥沙研究》2023,38(5):698-710
Every dam or barrage construction affects the watercourse and the retention of sediment that previously was carried by the river, which can lead to siltation of the reservoir and obstruction of water intakes over time, reducing their capacities. However, the information available regarding the effect of sediment and drawdown parameters, sediment management at reservoirs, as well as different equational approaches, is scarce. The current research aims to evaluate the effect of parameters associated with numerical modeling of sediment management in reservoirs considering scenarios with different drawdowns, transport equations, sediment size distributions, and thickness of the initial sediment layer. The case study of the Aimorés Hydropower Plant (HPP) is used, applying the Delft3D-FLOW model for two-dimensional modeling. All parameters influenced the volume of mobilized sediment, among which the initial layer thickness was the parameter that resulted in the greatest changes in simulated results. In general, the results show that the uncertainties in the input parameters outweigh the uncertainties between the techniques, which found large variations in results when evaluating the use of different transport equations. These results indicate the importance of proper estimation of model parameters for predicting effects with accuracy and the need for such studies before planning and management operations are evaluated to avoid environmental harm and energy waste. 相似文献
15.
To serve as a tool in the long term evaluation of the risk of accumulation of microbial contaminants (bacteria and viruses) entering soil and groundwater, a mathematical model is developed to predict the spatial and temporal distribution of pollutant concentration. The governing equation for bacterial transport is coupled with a transport equation for the bacterial nutrient present in the seeping wastewater. The deposition and declogging mechanisms are incorporated into the model as a rate process for bacteria and as an equilibrium partitioning for viruses. While the decay is assumed to be a first order reaction and the growth of bacteria is assumed to follow the Monod equation, the model equations exhibit nonlinearity and coupling. A simplified set of equations is solved analytically to test the numerical results. Coupled numerical solutions in one and two dimensions are obtained by the Galerkin method at spatial and temporal locations of interest. Cases studied included a soil column and a horizontal two-dimensional field coupled with the one dimensional solution. For these examples, the bacteria are removed almost totally within the top 7 cm of soil with minimal risk of clogging. 相似文献
16.
《Advances in water resources》1996,19(1):29-47
Transport processes in heterogeneous porous media are often treated in terms of one-equation models. Such treatment assumes that the velocity, pressure, temperature, and concentration can be represented in terms of a single large-scale averaged quantity in regions having significantly different mechanical, thermal, and chemical properties. In this paper we explore the process of single-phase flow in a two-region model of heterogeneous porous media. The region-averaged equations are developed for the case of a slightly compressible flow which is an accurate representation for a certain class of liquid-phase flows. The analysis leads to a pair of transport equations for the region averaged pressures that are coupled through a classic exchange term, in addition to being coupled by a diffusive cross effect. The domain of validity of the theory has been identified in terms of a series of length and timescale constraints.In Part II the theory is tested, in the absence of adjustable parameters, by comparison with numerical experiments for transient, slightly compressible flow in both stratified and nodular models of heterogeneous porous media. Good agreement between theory and experiment is obtained for nodular and stratified systems, and effective transport coefficients for a wide range of conditions are presented on the basis of solutions of the three closure problems that appear in the theory. Part III of this paper deals with the principle of large-scale mechanical equilibrium and the region-averaged form of Darcy's law. This form is necessary for the development and solution of the region-averaged solute transport equations that are presented in Part IV. Finally, in Part V we present results for the dispersion tensors and the exchange coefficient associated with the two-region model of solute transport with adsorption. 相似文献
17.
This paper presents a coupling of an ensemble Kalman filter (EnKF) with a discontinuous Galerkin-based, two-dimensional circulation model (DG ADCIRC-2DDI) to improve the state estimation of tidal hydrodynamics including water surface elevations and depth-integrated velocities. The methodology in this paper using EnKF perturbs the modeled hydrodynamics and bottom friction parameterization in the model while assimilating data with inherent error, and demonstrates a capability to apply EnKF within DG ADCIRC-2DDI for data assimilation. Parallel code development presents a unique aspect of the approach taken and is briefly described in the paper, followed by an application to a real estuarine system, the lower St. Johns River in north Florida, for the state estimation of tidal hydrodynamics. To test the value of gauge observations for improving state estimation, a tide modeling case study is performed for the lower St. Johns River successively using one of the four available tide gauging stations in model-data comparison. The results are improved simulations of water surface elevations and depth-integrated velocities using DG ADCIRC-2DDI with EnKF, both locally where data are available and non-locally where data are not available. The methodology, in general, is extensible to other modeling and data applications, for example, the use of remote sensing data, and specifically, can be readily applied as is to study other tidal systems. 相似文献
18.
A conceptual mathematical model was developed to describe the simultaneous transport (cotransport) of viruses and colloids in three-dimensional, water saturated, homogeneous porous media with uniform flow. The model accounts for the migration of individual virus and colloid particles as well as viruses attached onto colloids. Viruses can be suspended in the aqueous phase, attached onto suspended colloids and the solid matrix, and attached onto colloids previously attached on the solid matrix. Colloids can be suspended in the aqueous phase or attached on the solid matrix. Viruses in all four phases (suspended in the aqueous phase, attached onto suspended colloid particles, attached on the solid matrix, and attached onto colloids previously attached on the solid matrix) may undergo inactivation with different inactivation coefficients. The governing coupled partial differential equations were solved numerically using finite difference methods, which were implemented explicitly or implicitly so that both stability and speed factors were satisfied. Furthermore, the experimental data collected by Syngouna and Chrysikopoulos [1] were satisfactorily fitted by the newly developed cotransport model. 相似文献
19.
Transport problems occurring in porous media and including convection, diffusion and chemical reactions, can be well represented by systems of Partial Differential Equations. In this paper, a numerical procedure is proposed for the fast and robust solution of flow and transport problems in 2D heterogeneous saturated media. The governing equations are spatially discretized with unstructured triangular meshes that must satisfy the Delaunay condition. The solution of the flow problem is split from the solution of the transport problem and it is obtained with an approach similar to the Mixed Hybrid Finite Elements method, that always guarantees the M-property of the resulting linear system. The transport problem is solved applying a prediction/correction procedure. The prediction step analytically solves the convective/reactive components in the context of a MAST Finite Volume scheme. The correction step computes the anisotropic diffusive components in the context of a recently proposed Finite Elements scheme. Massa balance is locally and globally satisfied in all the solution steps. Convergence order and computational costs are investigated and model results are compared with literature ones. 相似文献