共查询到20条相似文献，搜索用时 15 毫秒
1.
In the evaluation of potential risk from ingestion of groundwater near an impacted site, numerical simulation of fate and transport processes of chemicals of concern is often required. If there is potential concern about multiple chemicals, numerical simulation of each chemical separately is often needed. In this paper, a semi-analytical solution is presented based on a numerical solution of the transport of a conservative and nonreactive tracer. When multiple chemicals undergoing sorption and first-order degradation need to be modeled, we can avoid performing individual numerical simulations for each chemical by applying the semi-analytical solution. Numerical test runs were conducted to verify the semi-analytical solution; simulation results reveal that the concentrations derived from the semi-analytical solution are identical to those derived from the individual numerical fate and transport model simulations. The semi-analytical solution requires steady-state flow conditions, no continuing contaminant source, and similar initial source concentration distributions. 相似文献
2.
Ground water management models require the repeated solution of a simulation model to identify an optimal solution to the management problem. Limited precision in simulation model calculations can cause optimization algorithms to produce erroneous solutions. Experiments are conducted on a transient field application with a streamflow depletion control management formulation solved with a response matrix approach. The experiment consists of solving the management model with different levels of simulation model solution precision and comparing the differences in optimal solutions obtained. The precision of simulation model solutions is controlled by choice of solver and convergence parameter and is monitored by observing reported budget discrepancy. The difference in management model solutions results from errors in computation of response coefficients. Error in the largest response coefficients is found to have the most significant impact on the optimal solution. Methods for diagnosing the adequacy of precision when simulation models are used in a management model framework are proposed. 相似文献
3.
介绍了钢筋混凝土框架结构计算机仿真和并行计算的研究现状。数值仿真主要采用有限元、离散元等数学物理模型,而可视化技术及图形仿真基于图形学和图像处理技术,是驾驭计算过程及理解大体积数据的唯一有效途径。有限元并行计算有SBS、EBE两种策略,非线性求解有直接与迭代解法,动力时程分析有显式、隐式和精细时程积分法。最后介绍了并行计算在钢筋混凝土结构分析中的应用。 相似文献
4.
This paper compares the performance of analytical and numerical approaches for modeling DNAPL dissolution with biodecay. A solution derived from a 1-D advective transport formulation (“Parker” model) is shown to agree very closely with high resolution numerical solutions. A simple lumped source mass balance solution in which with decay is assumed proportional to DNAPL mass (“Falta1” model) over- or underpredicts aqueous phase biodecay depending on the magnitude of the exponential factor governing the relationship between dissolution rate and DNAPL mass. A modification of the Falta model that assumes decay proportional to the source exit concentration is capable of accurately simulating source behavior with strong aqueous phase biodecay if model parameters are appropriately selected or calibrated (“Falta2” model). However, parameters in the lumped models exhibit complex interdependencies that cannot be quantified without consideration of transport processes within the source zone. Combining the Falta2 solution with relationships derived from the Parker model was found to resolve these limitations and track the numerical model results. A method is presented to generalize the analytical solutions to enable simulation of partial mass removal with changes in source parameters over time due to various remedial actions. The algorithm is verified by comparison with numerical simulation results. An example application is presented that demonstrates the interactions of partial mass removal, enhanced biodecay, enhanced mass transfer and source zone flow reduction applied at various time periods on contaminant flux reduction. Increasing errors that arise in numerical solutions with coarse discretization and high decay rates are shown to be controlled by using an adjusted decay coefficient derived from the Parker analytical solution. 相似文献
5.
An Explicit Finite Difference Model for Unconfined Aquifers 总被引：2，自引：0，他引：2
Most of the current simulation models for unconfined aquifers are based on the assumption that the free surface variation is small so that it can be combined with permeability to reduce the nonlinear Boussinesq equation to a linear partial differential equation (PDE). One of the most obvious reasons for using the linearization assumption is for the ease of numerical solution. This work presents a simpler alternative which permits an easy direct solution of the Boussinesq equation. A forward in time, central in space (FTCS) explicit finite difference method is used in the simulation model. The model was first validated by comparing its results with known analytical solution. It was then applied to an actual situation in which the short-term responses (from pumping) of an unconfined aquifer were simulated. The study shows that the stability of the model can be easily controlled, and because of the simple algorithm used, the code can be expeditiously developed and economically run on smaller machines. Due to the uncertainties in the calibration results, it is recommended here that more data be collected to improve the calibration before the model is used as a real-time simulation tool. 相似文献
6.
针对处理起伏地表(或含地下不规则波阻抗界面)条件下发展起来的地震波场数值模拟算法的模拟结果与解析解(大多数情形下无法得到)无法进行对比,且其有效性和正确性难以验证的情况,本文提出了一种可以相互验证波场数值模拟结果与射线追踪数值模拟结果的正确性和有效性的佐证方法,验证了参考射线追踪法.其中,波场数值模拟中采用曲线网格DRP/opt MacCormack有限差分法,射线追踪模拟则采用分区多步三角网格最短路径算法.通过系统对比上述两种方法得到的波场快照、单炮地震记录,以及合成理论地震图的结果显示,本方法相互作证了两种方法所得结果的正确性和有效性.双层和三层起伏层状模型的对比分析结果表明,这种方法不但可以加深理解地震波在复杂介质中的传播规律,同时射线法的引入为清晰识别和标定地震波场数值模拟中各种不同震相提供了一种便捷的途径. 相似文献
7.
The sliding behaviour of a rigid mass supported on a randomly vibrating foundation through a non-symmetric Coulomb-friction contact is studied both analytically and by numerical simulation. The analysis is based on a stationary solution of the associated Fokker-Planck equation, and makes use of equivalent linearization and of a suitable decomposition of the non-zero mean non-stationary sliding process. It is shown that the analytical results yield several exact asymptotic expansions for both small and large values of time. An extensive Monte Carlo type numerical simulation study produces non-stationary response statistics which are in very good accord with the analytical results. Furthermore, it is found that Gumbel's Extreme Value Distribution reproduces with remarkable accuracy the observed cumulative frequency of maximum slip displacement. The results of this paper may find application in seismic design of embankment dams, earth retaining walls and base ‘isolation’ systems. 相似文献
8.
Groundwater temperature is an important water quality parameter that affects species distributions in subsurface and surface environments. To investigate the response of subsurface temperature to atmospheric climate change, an analytical solution is derived for a one‐dimensional, transient conduction–advection equation and verified with numerical methods using the finite element code SUTRA. The solution can be directly applied to forward model the impact of future climate change on subsurface temperature profiles or inversely applied to produce a surface temperature history from measured borehole profiles. The initial conditions are represented using superimposed linear and exponential functions, and the boundary condition is expressed as an exponential function. This solution expands on a classic solution in which the initial and boundary conditions were restricted to linear functions. The exponential functions allow more flexibility in matching climate model projections (boundary conditions) and measured temperature–depth profiles (initial conditions). For example, measured borehole temperature data from the Sendai Plain and Tokyo, Japan, were used to demonstrate the improved accuracy of the exponential function for replicating temperature–depth profiles. Also, the improved accuracy of the exponential boundary condition was demonstrated using air temperature anomaly data from the Intergovernmental Panel on Climate Change. These air temperature anomalies were then used to forward model the effect of surficial thermal perturbations in subsurface environments with significant groundwater flow. The simulation results indicate that recharge can accelerate shallow subsurface warming, whereas upward groundwater discharge can enhance deeper subsurface warming. Additionally, the simulation results demonstrate that future groundwater temperatures obtained from the proposed analytical solution can deviate significantly from those produced with the classic solution. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
9.
10.
Haibing Shao Florian Centler Cecilia De Biase Martin Thullner Olaf Kolditz 《Advances in water resources》2009
In this comment we present a re-analysis of the analytical solution presented by Cirpka and Valocchi for steady-state concentrations of dissolved bioreactive compounds and bacterial biomass in porous media. We discuss the validity range of the analytical solution. In particular, the criterion used to determine the sustainability of biomass is revisited. This re-analysis shows that the ω criterion used by Cirpka and Valocchi is only a necessary but not a sufficient criterion to determine the bioreactive zones. As a consequence, the analytical solution does not provide the exact distribution of compounds throughout the domain, but can serve as upper or lower boundaries for species concentrations at a given location. These conclusions are supported by the simulation results obtained from an established reactive transport model. 相似文献
11.
A combined simulation–genetic algorithm (GA) optimization model is developed to determine optimal reservoir operational rule curves of the Nam Oon Reservoir and Irrigation Project in Thailand. The GA and simulation models operate in parallel over time with interactions through their solution procedure. A GA is selected as an optimization model, instead of traditional techniques, owing to its powerful and robust performance and simplicity in combining with a simulation technique. A GA is different from conventional optimization techniques in the way that it uses objective function information and does not require its derivatives, whereas in real‐world optimization problems the search space may include discontinuities and may often include a number of sub‐optimum peaks. This may cause difficulties for calculus‐based and enumerative schemes, but not in a GA. The simulation model is run to determine the net system benefit associated with state and control variables. The combined simulation–GA model is applied to determine the optimal upper and lower rule curves on a monthly basis for the Nam Oon Reservoir, Thailand. The objective function is maximum net system benefit subject to given constraints for three scenarios of cultivated areas. The monthly release is calculated by the simulation model in accordance with the given release policy, which depends on water demand. The optimal upper and lower rule curves are compared with the results of the HEC‐3 model (Reservoir System Analysis for Conservation model) calculated by the Royal Irrigation Department, Thailand, and those obtained using the standard operating policy. It was found that the optimal rule curves yield the maximum benefit and minimum damages caused by floods and water shortages. The combined simulation–GA model shows an excellent performance in terms of its optimization results and efficient computation. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
12.
随着地震动数据数量的增长和质量的提高,将基于数据驱动的机器学习方法应用到地震动模拟中有重要意义。以2021年5月21日云南漾濞MS6.4地震为例,利用主成分析方法从前震及余震地震动记录中提取特征母波时程,将地震动三要素作为模拟误差约束,在求解母波的线性组合系数时使用多目标优化算法寻优,最终找到帕累托最优解作为模拟目标台站记录时的组合系数,得到模拟地震动时程。结果表明:主成分析法在对实际地震动记录进行特征提取后,得到的特征母波时程可以在一定程度上保留原始数据的主要信息;考虑幅值、频谱和持时这三要素的角度去控制模拟误差,可以使得模拟的地震动时程更加接近真实记录。提出的基于特征提取的地震动模拟方法可以为基于小震数据合成大震地震动提供参考。 相似文献
13.
Derived herein is the integral representation solution of a Rayleigh-damped Bernoulli–Euler beam subjected to multi-support motion, which is free from calculation of a quasi-static solution, and in which the modal participation factor for support motion is formulated as a boundary modal reaction, thus making efficient calculation feasible. Three analytical methods, including (1) the quasi-static decomposition method, (2) the integral representation with the Cesàro sum technique, and (3) the integral representation in conjunction with Stokes' transformation, are presented. Two additional numerical methods of (4) the large mass FEM simulation technique and (5) large stiffness FEM simulation technique are easily incorporated into a commercial program to solve the problem. It is found that the results obtained by using these five methods are in good agreement, and that both the Cesàro sum and Stokes' transformation regularization techniques can extract the finite part of the divergent series of the integral representation. In comparison with the Mindlin method and Cesàro sum technique, Stokes' transformation is the best way because it is not only free of calculation of the quasi-static solution, but also because it can obtain the convergence rate as rapidly as the mode acceleration method can. 相似文献
14.
Baha Y. Mirghani Kumar G. Mahinthakumar Michael E. Tryby Ranji S. Ranjithan Emily M. Zechman 《Advances in water resources》2009
Groundwater characterization involves the resolution of unknown system characteristics from observation data, and is often classified as an inverse problem. Inverse problems are difficult to solve due to natural ill-posedness and computational intractability. Here we adopt the use of a simulation–optimization approach that couples a numerical pollutant-transport simulation model with evolutionary search algorithms for solution of the inverse problem. In this approach, the numerical transport model is solved iteratively during the evolutionary search. This process can be computationally intensive since several hundreds to thousands of forward model evaluations are typically required for solution. Given the potential computational intractability of such a simulation–optimization approach, parallel computation is employed to ease and enable the solution of such problems. In this paper, several variations of a groundwater source identification problem is examined in terms of solution quality and computational performance. The computational experiments were performed on the TeraGrid cluster available at the National Center for Supercomputing Applications. The results demonstrate the performance of the parallel simulation–optimization approach in terms of solution quality and computational performance. 相似文献
15.
16.
在地震动数值模拟方法中,谱元法和有限元法是应用较广泛的两种方法。基于经典的Lamb问题模型,首先推导给出地表竖向位移的解析解答。然后分别利用常用的四阶谱元法和线性有限元法,模拟了地表脉冲力源作用下模型的位移响应。考虑有意义的最短波长内的采样点个数及单元高宽比的变化,对比了两种方法的模拟精度;结果表明:对于谱元法,观测点与波源之间需至少包含两个网格,在此条件下,最短波长内包含一个网格（最短波长内5个采样点）时,数值解与解析解的误差小于1%,已达很高的精度;对于有限元法,最短波长内需包含10个网格时才能达到这一精度。此外,在满足网格尺寸要求的前提下,单元水平向与垂直向尺寸的比值在1∶1到5∶1的范围内时,谱元法和有限元法的模拟精度均变化不大。因此,单位波长内采样点个数相同时,谱元法的模拟精度比有限元法高的多,同时,在一定范围内两种方法的模拟结果对于宽高比的变化不敏感。 相似文献
17.
Advantages of the analytic element method for the solution of groundwater management problems 总被引：1，自引：0，他引：1
In the simulation‐optimization approach, a coupled optimization and groundwater flow/transport model is used to solve groundwater management problems. The efficiency of the numerical method, which is used to simulate the groundwater flow, is one the major reason to obtain the best solution for a management problem. This study was carried out to examine the advantages of the analytic element method (AEM) in the simulation‐optimization approach, for the solution of groundwater management problems. For this study, the AEM and finite difference method (FDM) based flow models were developed and coupled with the particle swarm optimization (PSO)‐based optimization model. Furthermore, the AEM‐PSO and FDM‐PSO models developed were applied in hypothetical as well as real field conditions to address groundwater management problems and the results were compared. For the real field situation, the models developed were applied to the Dore River basin in France to minimize the installation and operational cost of new pumping wells taking the location and discharge of the pumping wells as decision variables. The constraints of the problem were identified with the help of stakeholders and water authority officials. The AEM flow model was developed to facilitate the management model particularly when at each iteration, the optimization model calls for a simulation model to calculate the values of groundwater heads. The results show that, at some points, the AEM‐PSO model is efficient in identifying the optimal location of wells and consequently results in optimal costs, sometimes difficult when using the FDM. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
18.
The interaction between the fundamental mode surface Rayleigh waves and the buried heterogeneities with various sizes and different velocity contrasts was studied on base numerical simulation. The field of surface oscillations in the proximity of the scattering heterogeneities was computed as a function of frequency. The synthetic seismograms were used for numerical simulation of the microseismic sounding technology proposed earlier, implying that the solution of the inverse problem for the structure of the medium containing inclusions can be derived from the information contained in the ambient microseismic field. It is assumed that the depth of the layer to be reconstructed is linked with the frequency of the microseisms by a simple relation with the help of a numerical coefficient equal to 0.4–0.5. The combined results of the simulation of a direct problem together with the simple inverse problem solution show that the microseismic sounding technique ensures adequate estimation of the medium structure. Previously, the technology was based on the experimental data only and was phenomenological in character. Some relations between the velocity parameters of the original model heterogeneities and their reconstructed images were also studied. 相似文献
19.
20.
断层演化过程中的周期和非周期地震滑移及非地震滑移一维断层模型分析 总被引：2，自引：0，他引：2
本文首先根据Dieterich和Ruina提出的含速率和状态的摩擦定律(Dieterich-Ruina定律), 基于一维弹簧-滑块模型推导了地震复发周期的解析表达式, 然后将该近似解与数值模拟结果以及Barbot等的相关研究进行了对比分析． 此外, 本文还利用数值模拟与理论分析研究了断层周期和非周期演化的力学成因机制以及非地震滑移形成的另类力学机制, 并讨论了一维弹簧-滑块模型的优点及其局限性． 结果表明: ① 震后滑移和自加速/成核阶段的持续时间在整个演化过程中不能被忽略; ② 在修正后的复发周期模型中, 复发周期的长短除了与断层特征尺度、 作用于断层面上的有效正应力和远场加载速率相关外, 还受Dieterich-Ruina定律中摩擦参数的取值以及临界滑移距离的影响; ③ 当给定各个物理参数和几何参数时, 目前所得到的解析近似解可以很好地估计地震的复发周期, 其相对误差可小于5%; ④ 在断层演化过程中, 施加剪切应力加载会产生非周期的地震滑移, 而在自加速/成核阶段后期或震后滑移阶段早期, 施加较大的剪切应力加载, 则会出现非地震滑移． 相似文献