共查询到20条相似文献，搜索用时 46 毫秒
1.
Iterative solvers preconditioned with algebraic multigrid have been devised as an optimal technology to speed up the response of large sparse linear systems. In this work, this technique was implemented in the framework of the dual delineation approach. This involves a single groundwater flow linear solution and a pure advective transport solution with different right-hand sides. The new solver was compared with other preconditioned iterative methods, the MODFLOW's GMG solver, and direct sparse solvers. Test problems include two- and three-dimensional benchmarks spanning homogeneous and highly heterogeneous and anisotropic formations. For the groundwater flow problems, using the algebraic multigrid preconditioning speeds up the numerical solution by one to two orders of magnitude. The algebraic multigrid preconditioner efficiency was preserved for the three dimensional heterogeneous and anisotropic problem unlike for the MODFLOW's GMG solver. Contrarily, a sparse direct solver was the most efficient for the pure advective transport processes such as the forward travel time simulations. Hence, the best sparse solver for the more general advection-dispersion transport equation is likely to be Péclet number dependent. When equipped with the best solvers, processing multimillion grid blocks by the dual delineation approach is a matter of seconds. This paves the way for its routine application to large geological models. The paper gives practical hints on the strategies and conditions under which algebraic multigrid preconditioning would remain competitive for the class of nonlinear and/or transient problems. 相似文献
2.
The χMD matrix solver package is incorporated into USGS groundwater modeling software, such as MODFLOW-NWT, MODFLOW-USG, and MT3D. The solver is used to solve matrices assembled through numerical discretization of the groundwater flow equation, and solute transport equations. χMD has demonstrated its higher robustness, faster execution speed, and more efficient memory usage compared to the existing solvers for many types of groundwater flow problems. χMD uses preconditioned iterative Krylov-subspace methods and consists of preconditioning and acceleration modules. Because the solver package uses a variety of preconditioning features including level-based incomplete lower-upper (ILU) factorization method with a drop tolerance scheme, users must choose optimal preconditioning parameters to improve execution speed and robustness. In order to examine how the preconditioning parameters, ILU factorization level, and drop tolerance values affect the overall performance of the matrix solver, we evaluated five different groundwater model applications using MODFLOW-USG that include different numerical complexities. For those five cases, the number of discretization nodes varied from 10,000 cells to 730,300 cells. From the analysis, we found that the preconditioning parameters greatly affect execution times and memory usage of the preconditioning and acceleration procedures. In addition, a combination of the ILU level between five to seven and the drop tolerance value between 10−2 and 10−3 usually resulted in shorter overall execution time. Our study suggests that the users can elicit higher performance and robustness of the χMD matrix solver using this combination of the parameters and enhance computational efficiency of solving groundwater and solute transport problems. 相似文献
3.
The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning. 相似文献
4.
三维反演解释是电磁法勘探发展的重要趋势,而如何提高三维反演的可靠性、稳定性和计算效率是算法开发者们目前的研究重点.本文实现了一种频率域可控源电磁(CSEM)三维反演算法.其中正演基于拟态有限体积法离散化,利用直接矩阵分解技术来求解大型线性系统方程,不仅准确、稳定,而且特别有利于含有大量发射场源位置的CSEM勘探情况;对目标函数的最优化采用高斯牛顿法(GN),具有近似二次的收敛性;使用预条件共轭梯度法(PCG)求解每次GN迭代所得到的法方程,避免了显式求解和存储灵敏度矩阵,减小了计算量.以上这些方法的结合应用,使得本文的三维反演算法准确、稳定且高效.通过陆地和海洋CSEM勘探场景中的典型理论模型的反演测试,验证了本文算法的有效性. 相似文献
5.
Comparison of an algebraic multigrid algorithm to two iterative solvers used for modeling ground water flow and transport 总被引：1，自引：0，他引：1
Numerical solution of large-scale ground water flow and transport problems is often constrained by the convergence behavior of the iterative solvers used to solve the resulting systems of equations. We demonstrate the ability of an algebraic multigrid algorithm (AMG) to efficiently solve the large, sparse systems of equations that result from computational models of ground water flow and transport in large and complex domains. Unlike geometric multigrid methods, this algorithm is applicable to problems in complex flow geometries, such as those encountered in pore-scale modeling of two-phase flow and transport. We integrated AMG into MODFLOW 2000 to compare two- and three-dimensional flow simulations using AMG to simulations using PCG2, a preconditioned conjugate gradient solver that uses the modified incomplete Cholesky preconditioner and is included with MODFLOW 2000. CPU times required for convergence with AMG were up to 140 times faster than those for PCG2. The cost of this increased speed was up to a nine-fold increase in required random access memory (RAM) for the three-dimensional problems and up to a four-fold increase in required RAM for the two-dimensional problems. We also compared two-dimensional numerical simulations of steady-state transport using AMG and the generalized minimum residual method with an incomplete LU-decomposition preconditioner. For these transport simulations, AMG yielded increased speeds of up to 17 times with only a 20% increase in required RAM. The ability of AMG to solve flow and transport problems in large, complex flow systems and its ready availability make it an ideal solver for use in both field-scale and pore-scale modeling. 相似文献
6.
电磁法的三维数值模拟是一个对数值算法和计算机硬件要求都非常高的问题.对常用的微分类方法如有限单元法和有限差分法而言,求解最后所得的大型线性方程组是至关重要的一步,直接影响到正演算法的实用性.如何高效、稳定且准确地解线性方程长期以来一直是被探讨的问题.本文实现了基于线性系统直接求解技术的频率域可控源电磁(CSEM)三维正演.使用交错网格有限体积法(FV)来离散化关于二次电场的Helmholtz方程;使用直接解法取代传统的迭代解法来求解离散线性系统,即对系统矩阵进行完全LU分解,具体通过调用大规模并行矩阵直接求解器(MUMPS)来实现.基于理论模型做了一系列数值实验,首先证明了直接解法的高精度和稳定性,并考察了其内存需求、计算时间和并行可伸缩性等主要计算性能,最后检验了所开发的算法快速模拟多场源CSEM问题的能力以及对常规海洋和陆地CSEM模拟的有效性. 相似文献
7.
《Advances in water resources》2001,24(6):621-630
One of the more advanced approaches for simulating groundwater flow in fractured porous media is the discrete-fracture approach. This approach is limited by the large computational overheads associated with traditional modeling methods. In this work, we apply the Lanczos reduction method to the modeling of groundwater flow in fractured porous media using the discrete-fracture approach. The Lanczos reduction method reduces a finite element equation system to a much smaller tridiagonal system of first-order differential equations. The reduced system can be solved by a standard tridiagonal algorithm with little computational effort. Because solving the reduced system is more efficient compared to solving the original system, the simulation of groundwater flow in discretely fractured media using the reduction method is very efficient. The proposed method is especially suitable for the problem of large-scale and long-term simulation. In this paper, we develop an iterative version of Lanczos algorithm, in which the preconditioned conjugate gradient solver based on ORTHOMIN acceleration is employed within the Lanczos reduction process. Additional efficiency for the Lanczos method is achieved by applying an eigenvalue shift technique. The “shift” method can improve the Lanczos system convergence, by requiring fewer modes to achieve the same level of accuracy over the unshifted case. The developed model is verified by comparison with dual-porosity approach. The efficiency and accuracy of the method are demonstrated on a field-scale problem and compared to the performance of classic time marching method using an iterative solver on the original system. In spite of the advances, more theoretical work needs to be carried out to determine the optimal value of the shift before computations are actually carried out. 相似文献
8.
Three-dimensional finite-element modelling of magnetotelluric data with a divergence correction 总被引：4，自引：0，他引：4
A finite-element method for computing the electric field in a 3-D conductivity model of the Earth for plane wave sources, thus enabling magnetotelluric responses to be calculated, is presented. The method incorporates in the iterative solution of the electric-field system of equations the divergence correction technique introduced for finite-difference solutions by Smith (1996). The correction technique accelerates the development of the discontinuity of the normal component of the approximate electric field across conductivity discontinuities. The convergence rate of the iterative solution is improved significantly, especially for low frequencies. The correction technique involves computing the divergence of the current density for the approximate electric field, computing the static potential whose source is this divergence of the current density, and ‘correcting’ the approximate electric field by subtracting from it the gradient of the potential. This is repeated at regular intervals during the iterative solution of the electric-field system of equations. For the method presented here, the Earth model is discretised using a rectilinear mesh comprising uniform cells. Edge-element basis functions are used to approximate the electric field and nodal basis functions are used to approximate the correction potential. The Galerkin method is used to derive the systems of equations for the approximate electric field and correction potential from the respective differential equations. A bi-conjugate gradient solver was found to be adequate for the system of equations for the correction potential; a generalised minimum residual solver was found to be better for the electric-field system of equations. The method is illustrated using the COMMEMI 3D-1A and 3D-2A models. 相似文献
9.
《Advances in water resources》2001,24(6):595-605
Adaptive time stepping with embedded error control is applied to the mixed form of Richards equation. It is the first mathematically based adaptive scheme applied to this form of Richards equation. The key to the method is the approximation of the local truncation error of the scheme in terms of the pressure head, although, to enforce mass conservation, the principal time approximation is based on the moisture content. The time stepping scheme is closely related to an implicit Thomas–Gladwell approximation and is unconditionally stable and second-order accurate. Numerical trials demonstrate that the new algorithm fully automates stepsize selection and robustly constrains temporal discretisation errors given a user tolerance. The adaptive mechanism is shown to improve the performance of the non-linear solver, providing accurate initial solution estimates for the iterative process. Furthermore, the stepsize variation patterns reflect the adequacy of the spatial discretisation, here accomplished by linear finite elements. When sufficiently dense spatial grids are used, the time step varies smoothly, while excessively coarse grids induce stepsize oscillations. 相似文献
10.
For wave propagation simulation in piecewise heterogeneous media, Gaussian-elimination-based full-waveform solutions to the generalized Lippmann–Schwinger integral equation (GLSIE) are highly accurate, but involved with extremely time-consuming computations because of the very large size of the resulting boundary–volume integral equation matrix to be inverted. Several flexible approximations to the GLSIE are scaled in an iterative way to adapt numerical solutions to the smoothness of heterogeneous media in terms of incident wavelengths, with a great saving of computing time and memory. Among various typical iterative schemes to the GLSIE matrix, the generalized minimal residual method (GMRES) is an efficient approach to reduce the computational intensity to some degree. The most efficient approximation can be obtained using a Born series, as an alternative iterative solution, to both the boundary-scattering and volume-scattering waves, leading to the Born-series approximation (BSA) scheme and the improved Born-series approximation (IBSA) scheme. These iteration schemes are validated by dimensionless frequency responses to a heterogeneous semicircular alluvial valley, and then applied to a heterogeneous multilayered model by calculating synthetic seismograms to evaluate approximation accuracies. Numerical experiments, compared with the full-waveform numerical solution, indicate that the convergence rates of these methods decrease gradually with increasing velocity perturbations. The comparison also shows that the BSA scheme has a faster convergence than the GMRES method for velocity perturbations less than 10 percent, but converges slowly and even hardly achieves convergence for velocity perturbations greater than 15 percent. The IBSA scheme gives a superior performance over the other methods, with the least iterations to achieve the necessary convergence. 相似文献
11.
巨大的计算量是制约全波形反演（FWI）生产实用化的难题之一.为此,本文提出了一种高效的波场迭代解法,将其应用于频率域常密度声波方程FWI,并给出了详细的反演流程.通过建立用于波场迭代的目标函数,推导相应梯度、步长公式,新方法将反演中波场正传和残差波场反传过程转化为无约束优化问题,从理论上分析了新方法的计算效率显著高于常规FWI.在数值试验中,本文方法通过几次迭代便能获得高精度的正传、残差反传波场,收敛速度明显高于未经预处理的GMRES方法.进一步引入高效编码策略,新方法的计算时间约为常规编码FWI的1/8,与理论分析结果吻合（波场迭代次数为8,模型未知量个数约为7万）,且波场迭代次数为6时,反演效果已与常规编码FWI相近. 相似文献
12.
基于网格精度、形成系统方程的方式、边界条件以及预条件线性算子等,文中对大地电磁（MT）有限差分数值解作了对比.对不同网格剖分方式下的三个均匀半空间模型的一维MT响应对比显示,在降低首层厚度的同时保持层间厚度变化在合理范围可以同时提高主场和辅助场的精度.在利用正常中心网格法（主场和辅助场都定义在单元顶面的中心）计算二维（2-D）TM模式响应时,应该从Maxwell一次差分方程开始组建二次差分方程,这样可以更充分考虑模型电阻率的变化.在对边界值如何影响数值解的测试表明,仅仅提高一维（1-D）边界值的精度对提高2-D MT有限差分数值解的精度是有限的.线性算子对提高MT解的效率十分重要,简单的对比进一步表明合适的预条件再配合好的线性算子（如文中求解2-D MT时所采用的DILU-BICGSTAB方法）不仅可以加速收敛,而且可以降低迭代次数. 相似文献
13.
Qianqian Liu Gulimire Hanati Sulitan Danierhan Yin Zhang Zhiping Zhang 《Ground water》2019,57(6):969-979
As a component of arid ecosystems, groundwater plays an important role in plant growth; therefore, it is essential to use deterministic models to reconstruct the process of groundwater level change. Typically, the linearized solution of the one-dimensional (1-D) Boussinesq equation yields acceptable performance in simulating transient conditions over short recharge periods in ephemeral stream systems, but the ability of this solution to simulate multiyear changes in groundwater levels is limited. In this study, an improved groundwater hydraulics (GH-D2) model is built based on the groundwater hydraulics (GH) solution of the 1-D Boussinesq equation to simulate multiyear changes in the groundwater level in ephemeral stream systems. The model is validated in the lower reaches of the Tarim River to simulate groundwater level fluctuations within the scope of influence of the river (300, 500, 750, 1050 m) over a 16-year period (2000 to 2015). To evaluate the performance of the models, the bias, mean absolute error, root mean squared error, Nash-Sutcliffe efficiency (NSE), and coefficient of determination (R2) are calculated. The results show that the improved GH-D2 model, which considers ephemeral streamflow, unsteady flow theory and the delayed response effect of groundwater level changes, performs well in simulating multiyear changes in the groundwater level in the ephemeral stream system. The observed and simulated values of the groundwater level at different river distances are consistent, and the model provides a new basis for multiyear simulations of groundwater level fluctuations in ephemeral stream systems. 相似文献
14.
本研究针对大地电磁测深法有限元数值模拟中,迭代法求解线性方程组效率较低的问题,利用亥姆霍兹分解原理,将电场矢量双旋度方程的预条件问题转化为基于矢量位的泊松问题和基于标量位的拉普拉斯问题,并在四面体非结构化棱边元离散的情况下,借助节点元辅助网格离散上述预条件问题,进一步利用代数多重网格方法（AMG）实施求解,最终实现预条件算法.利用经典的COMMEMI理论模型进行试算并与前人的积分方程解进行对比,验证了本文数值模拟程序与预条件方法的正确性和可靠性.此外,利用不同自由度规模的实验模型对这一预条件算法的效率进行了测试.结果表明,这一算法可以有效地提升大地电磁测深法棱边有限元数值模拟迭代法的收敛性,计算效率较通用的不完全LU分解预条件算法明显更高;在较大自由度网格（>1000万）数值模拟计算中,其算法效率及内存占用相对直接解法有较大优势,也使小型工作站上利用较大自由度的有限元网格进行大地电磁测深数值模拟计算成为可能. 相似文献
15.
地震直达波走时层析成像可归结为求解一个大型的、稀疏的、常常是病态的线性方程组.求解方程组常用的迭代法,需要一个比较合理的初始猜测解,也即是初始速度模型.初始模型关系到反演的效率甚至成像的正确性.本文在前人研究基础上提出一种生成模型网格节点初始速度方法,假定震源到检波点路径为直线,记录每条射线穿过的单元和统计每个网格单元穿过的射线数目、自动拾取网格节点所在单元的数目等.实例中,由原始数学模型的正演旅行时资料生成节点初始速度模型,效果可以.最后,分别采用均匀模型和本文方法生成的初始模型进行迭代反演,通过比较,证实该自动生成节点初始模型的可行性和可靠性,并对存在的问题进行讨论和解释. 相似文献
16.
Increasing numerical efficiency of iterative solution for total least-squares in datum transformations 总被引：1，自引：0，他引：1
Cartesian coordinate transformation between two erroneous coordinate systems is considered within the Errors-In-Variables (EIV) model. The adjustment of this model is usually called the total Least-Squares (LS). There are many iterative algorithms given in geodetic literature for this adjustment. They give equivalent results for the same example and for the same user-defined convergence error tolerance. However, their convergence speed and stability are affected adversely if the coefficient matrix of the normal equations in the iterative solution is ill-conditioned. The well-known numerical techniques, such as regularization, shifting-scaling of the variables in the model, etc., for fixing this problem are not applied easily to the complicated equations of these algorithms. The EIV model for coordinate transformations can be considered as the nonlinear Gauss-Helmert (GH) model. The (weighted) standard LS adjustment of the iteratively linearized GH model yields the (weighted) total LS solution. It is uncomplicated to use the above-mentioned numerical techniques in this LS adjustment procedure. In this contribution, it is shown how properly diminished coordinate systems can be used in the iterative solution of this adjustment. Although its equations are mainly studied herein for 3D similarity transformation with differential rotations, they can be derived for other kinds of coordinate transformations as shown in the study. The convergence properties of the algorithms established based on the LS adjustment of the GH model are studied considering numerical examples. These examples show that using the diminished coordinates for both systems increases the numerical efficiency of the iterative solution for total LS in geodetic datum transformation: the corresponding algorithm working with the diminished coordinates converges much faster with an error of at least 10-5 times smaller than the one working with the original coordinates. 相似文献
17.
To assess the post-earthquake seismic safety of buildings, it is crucial to predict seismic response, and it is necessary to set the appropriate physical parameters of the response analysis model. Numerous methods have been proposed to identify physical parameters. However, most of them are limited to linear systems, and previous researches on nonlinear systems have difficulties in practical applications. In this paper, a nonlinear response analysis model is identified for a full-scale ten-story reinforced concrete building with the degrading tri-linear stiffness model by the modal iterative error correction (MIEC) method, and the accuracy of this technique is discussed by comparing with the shaking table test. 相似文献
18.
An iterative solution method is presented and illustrated to analyse the dynamic response of bridge–vehicle systems. The method consists in dividing the whole system into 2 subsystems at the interface of the bridge and vehicles; these 2 subsystems are solved separately; their compatibility at the interface is achieved by an iterative procedure with under-relaxation or with Aitken acceleration. The characteristics of this method are explained on a simplified system with 2 degrees of freedom (DOF). The numerical results for a simple example demonstrate the high performances of the proposed method: good convergence rate and high accuracy. Finally, the method is applied to a practical example: the linear dynamic response of the Yangtze-River Bridge at Wuhan under a moving train with 2 locomotives and 4 freight cars. The efficiency is attained because neither formation nor factorisation of the coefficient matrices for the equations of the system are needed at every time step in linear analysis. The Aitken acceleration technique is more efficient in systems with multi-degrees of freedom than the relaxation technique. The proposed method will be even more efficient in non-linear dynamic response because, in this case, the iterations are necessary whether the system is solved as a whole or not. 相似文献
19.
As a result of rock dissolution processes, karst aquifers exhibit highly conductive features such as caves and conduits. Within these structures, groundwater flow can become turbulent and therefore be described by nonlinear gradient functions. Some numerical groundwater flow models explicitly account for pipe hydraulics by coupling the continuum model with a pipe network that represents the conduit system. In contrast, the Conduit Flow Process Mode 2 (CFPM2) for MODFLOW-2005 approximates turbulent flow by reducing the hydraulic conductivity within the existing linear head gradient of the MODFLOW continuum model. This approach reduces the practical as well as numerical efforts for simulating turbulence. The original formulation was for large pore aquifers where the onset of turbulence is at low Reynolds numbers (1 to 100) and not for conduits or pipes. In addition, the existing code requires multiple time steps for convergence due to iterative adjustment of the hydraulic conductivity. Modifications to the existing CFPM2 were made by implementing a generalized power function with a user-defined exponent. This allows for matching turbulence in porous media or pipes and eliminates the time steps required for iterative adjustment of hydraulic conductivity. The modified CFPM2 successfully replicated simple benchmark test problems. 相似文献