共查询到20条相似文献，搜索用时 703 毫秒
1.
We develop a new approach for solving the nonlinear Richards’ equation arising in variably saturated flow modeling. The growing complexity of geometric models for simulation of subsurface flows leads to the necessity of using unstructured meshes and advanced discretization methods. Typically, a numerical solution is obtained by first discretizing PDEs and then solving the resulting system of nonlinear discrete equations with a Newton-Raphson-type method. Efficiency and robustness of the existing solvers rely on many factors, including an empiric quality control of intermediate iterates, complexity of the employed discretization method and a customized preconditioner. We propose and analyze a new preconditioning strategy that is based on a stable discretization of the continuum Jacobian. We will show with numerical experiments for challenging problems in subsurface hydrology that this new preconditioner improves convergence of the existing Jacobian-free solvers 3-20 times. We also show that the Picard method with this preconditioner becomes a more efficient nonlinear solver than a few widely used Jacobian-free solvers. 相似文献
2.
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. 相似文献
3.
《Advances in water resources》1998,21(3):237-249
The transport and fate of reactive chemicals in groundwater is governed by equations which are often difficult to solve due to the nonlinear relationship between the solute concentrations for the liquid and solid phases. The nonlinearity may cause mass balance errors during the numerical simulation in addition to numerical errors for linear transport system. We have generalized the modified Picard iteration algorithm of Celia et al.5 for unsaturated flow to solve the nonlinear transport equation. Written in a ‘mixed-form’ formulation, the total solute concentration is expanded in a Taylor series with respect to the solution concentration to linearize the transport equation, which is then solved with a conventional finite element method. Numerical results of this mixed-form algorithm are compared with those obtained with the concentration-based scheme using conventional Picard iteration. In general, the new solver resulted in negligible mass balance errors (< ∥10−8∥%) and required less computational time than the conventional iteration scheme for the test examples, including transport involving highly nonlinear adsorption under steady-state as well as transient flow conditions. In contrast, mass balance errors resulting from the conventional Picard iteration method were higher than 10% for some highly nonlinear problems. Application of the modified Picard iteration scheme to solve the nonlinear transport equation may greatly reduce the mass balance errors and increase computational efficiency. 相似文献
4.
N. O. Weiss 《地球物理与天体物理流体动力学》2013,107(1-4):229-247
Abstract Cowling investigated the effect of an imposed magnetic field on convection in order to explain the origin of sunspots. After summarizing the classical linear theory of Boussinesq magnetoconvection, this review proceeds to more recent nonlinear results. Weakly nonlinear theory is used to establish the relevant bifurcation structure, which involves steady, oscillatory and chaotic solutions. Behaviour found in numerical experiments can then be related to these analytical results. Thereafter, attention is focused on the astrophysically relevant problem of fully compressible magnetoconvection. Steady two-dimensional nonlinear solutions show two important effects: stratification introduces an asymmetry between rising and falling fluid, while compressibility leads to evacuated magnetic flux sheets. Time-dependent behaviour includes transitions between standing waves and travelling waves, as well as changes in horizontal scale, leading to the development of more complicated spatial structures. Work on three-dimensional models, which is now in progress, will lead to a better understanding of the structure of a sunspot. 相似文献
5.
We examine the implementation of a wave-breaking mechanism into a nonlinear potential flow solver. The success of the mechanism will be studied by implementing it into the numerical model HOS-NWT, which is a computationally efficient, open source code that solves for the free surface in a numerical wave tank using the high-order spectral (HOS) method. Once the breaking mechanism is validated, it can be implemented into other nonlinear potential flow models. To solve for wave-breaking, first a wave-breaking onset parameter is identified, and then a method for computing wave-breaking associated energy loss is determined. Wave-breaking onset is calculated using a breaking criteria introduced by Barthelemy et al. (J Fluid Mech https://arxiv.org/pdf/1508.06002.pdf, submitted) and validated with the experiments of Saket et al. (J Fluid Mech 811:642–658, 2017). Wave-breaking energy dissipation is calculated by adding a viscous diffusion term computed using an eddy viscosity parameter introduced by Tian et al. (Phys Fluids 20(6): 066,604, 2008, Phys Fluids 24(3), 2012), which is estimated based on the pre-breaking wave geometry. A set of two-dimensional experiments is conducted to validate the implemented wave breaking mechanism at a large scale. Breaking waves are generated by using traditional methods of evolution of focused waves and modulational instability, as well as irregular breaking waves with a range of primary frequencies, providing a wide range of breaking conditions to validate the solver. Furthermore, adjustments are made to the method of application and coefficient of the viscous diffusion term with negligible difference, supporting the robustness of the eddy viscosity parameter. The model is able to accurately predict surface elevation and corresponding frequency/amplitude spectrum, as well as energy dissipation when compared with the experimental measurements. This suggests the model is capable of calculating wave-breaking onset and energy dissipation successfully for a wide range of breaking conditions. The model is also able to successfully calculate the transfer of energy between frequencies due to wave focusing and wave breaking. This study is limited to unidirectional waves but provides a valuable basis for future application of the wave-breaking model to a multidirectional wave field. By including parameters for removing energy due to wave-breaking into a nonlinear potential flow solver, the risk of developing numerical instabilities due to an overturning wave is decreased, thereby increasing the application range of the model, including calculating more extreme sea states. A computationally efficient and accurate model for the generation of a nonlinear random wave field is useful for predicting the dynamic response of offshore vessels and marine renewable energy devices, predicting loads on marine structures, and in the study of open ocean wave generation and propagation in a realistic environment. 相似文献
6.
本文利用MHD二维不可压模式,研究了地球磁层顶边界区剪切流引起的Kelvin-Helmholtz（K-H）不稳定性问题,得到了一个新的非线性微分方程组.理论和数值分析表明:该问题的非线性演化对初值非常敏感,而且在雷诺数和磁雷诺数给定的条件下,Alfven马赫数（MA）对K-H不稳定性的非线性演化起决定性作用.这组方程蕴含几个吸引子,如不动点,极限环和奇异吸引子等,这体现了磁层顶非线性系统的复杂性.文中还发现背景磁场在磁层顶K-H不稳定性的非线性演化过程中起很重要的作用. 相似文献
7.
An Experimenting Field Approach for the Numerical Solution of Multiphase Flow in Porous Media 
下载免费PDF全文
下载免费PDF全文 In this work, we apply the experimenting pressure field technique to the problem of the flow of two or more immiscible phases in porous media. In this technique, a set of predefined pressure fields are introduced to the governing partial differential equations. This implies that the velocity vector field and the divergence at each cell of the solution mesh can be determined. However, since none of these fields is the true pressure field entailed by the boundary conditions and/or the source terms, the divergence at each cell will not be the correct one. Rather the residue which is the difference between the true divergence and the calculated one is obtained. These fields are designed such that these residuals are used to construct the matrix of coefficients of the pressure equation and the right‐hand side. The experimenting pressure fields are generated in the solver routine and are fed to the different routines, which may be called physics routines, which return to the solver the elements of the matrix of coefficients. Therefore, this methodology separates the solver routines from the physics routines and therefore results in simpler, easy to construct, maintain, and update algorithms. 相似文献
8.
A Parallel PCG Solver for MODFLOW 总被引：2，自引：0，他引：2
In order to simulate large-scale ground water flow problems more efficiently with MODFLOW, the OpenMP programming paradigm was used to parallelize the preconditioned conjugate-gradient (PCG) solver with in this study. Incremental parallelization, the significant advantage supported by OpenMP on a shared-memory computer, made the solver transit to a parallel program smoothly one block of code at a time. The parallel PCG solver, suitable for both MODFLOW-2000 and MODFLOW-2005, is verified using an 8-processor computer. Both the impact of compilers and different model domain sizes were considered in the numerical experiments. Based on the timing results, execution times using the parallel PCG solver are typically about 1.40 to 5.31 times faster than those using the serial one. In addition, the simulation results are the exact same as the original PCG solver, because the majority of serial codes were not changed. It is worth noting that this parallelizing approach reduces cost in terms of software maintenance because only a single source PCG solver code needs to be maintained in the MODFLOW source tree. 相似文献
9.
全文以地震道反演的数值实验为例,根据混沌理论讨论了非线性地震反演的以下特性:1.对于带误差的地震数据和极其平滑的初始模型,逐次线性化的迭代过程产生的输出序列最终走向无序,这是非线性反演系统本身的特性决定的;2.迭代（非线性反演）系统是由Poin-care映的或系统方程描述的;3.可以根据地震反问题的特点,用多种不同的方法定义相应的Lyapunov指数,它们的数值和组合是非线性反演系统状态的有效指示;4.取决于系统参数的选取,发现在二维相空间有奇异吸引子的对应物.总之,从混沌理论的角度来研究纯数学的非线性地震反问题,可以揭示某些过去鲜为人知的内在规律性.本文为全文的上篇,主要讨论前两个问题,即理论和地震道反演数值试验结果. 相似文献
10.
Bellie Sivakumar 《地球表面变化过程与地形》2007,32(7):969-979
Whether or not river flow exhibits nonlinear determinism remains an unresolved question. While studies on the use of nonlinear deterministic methods for modeling and prediction of river flow series are on the rise and the outcomes are encouraging, suspicions and criticisms of such studies continue to exist as well. An important reason for this situation is that the correlation dimension method, used as a nonlinear determinism identification tool in most of those studies, may possess certain limitations when applied to real river flow series, which are always finite and often short and also contaminated with noise (e.g. measurement error). In view of this, the present study addresses the issue of nonlinear determinism in river flow series using prediction as a possible indicator. This is done by (1) reviewing studies that have employed nonlinear deterministic methods (coupling phase‐space reconstruction and local approximation techniques) for river flow predictions and (2) identifying nonlinear determinism (or linear stochasticity) based on the level of prediction accuracy in general, and on the prediction accuracy against the phase‐space reconstruction parameters in particular (termed as the ‘inverse approach’). The results not only provide possible indications to the presence of nonlinear determinism in the river flow series studied, but also support, both qualitatively and quantitatively, the low correlation dimensions reported for such. Therefore, nonlinear deterministic methods are a viable complement to linear stochastic ones for studying river flow dynamics, if sufficient caution is exercised in their applications and in interpreting the outcomes. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
11.
Tabea Broecker Vahid Sobhi Gollo Aryeh Fox Jörg Lewandowski Gunnar Nützmann Shai Arnon Reinhard Hinkelmann 《Ground water》2021,59(4):488-502
Transport processes that lead to exchange of mass between surface water and groundwater play a significant role for the ecological functioning of aquatic systems, for hydrological processes and for biogeochemical transformations. In this study, we present a novel integral modeling approach for flow and transport at the sediment–water interface. The model allows us to simultaneously simulate turbulent surface and subsurface flow and transport with the same conceptual approach. For this purpose, a conservative transport equation was implemented to an existing approach that uses an extended version of the Navier–Stokes equations. Based on previous flume studies which investigated the spreading of a dye tracer under neutral, losing and gaining flow conditions the new solver is validated. Tracer distributions of the experiments are in close agreement with the simulations. The simulated flow paths are significantly affected by in- and outflowing groundwater flow. The highest velocities within the sediment are found for losing condition, which leads to shorter residence times compared to neutral and gaining conditions. The largest extent of the hyporheic exchange flow is observed under neutral condition. The new solver can be used for further examinations of cases that are not suitable for the conventional coupled models, for example, if Reynolds numbers are larger than 10. Moreover, results gained with the integral solver provide high-resolution information on pressure and velocity distributions at the rippled streambed, which can be used to improve flow predictions. This includes the extent of hyporheic exchange under varying ambient groundwater flow conditions. 相似文献
12.
《Advances in water resources》2001,24(6):667-676
Gradient-based nonlinear programming (NLP) methods can solve problems with smooth nonlinear objectives and constraints. However, in large and highly nonlinear models, these algorithms can fail to find feasible solutions, or converge to local solutions which are not global. Evolutionary search procedures in general, and genetic algorithms (GAs) specifically, are less susceptible to the presence of local solutions. However, they often exhibit slow convergence, especially when there are many variables, and have problems finding feasible solutions in constrained problems with “narrow” feasible regions. In this paper, we describe strategies for solving large nonlinear water resources models management, which combine GAs with linear programming. The key idea is to identify a set of complicating variables in the model which, when fixed, render the problem linear in the remaining variables. The complicating variables are then varied by a GA. This GA&LP approach is applied to two nonlinear models: a reservoir operation model with nonlinear hydropower generation equations and nonlinear reservoir topologic equations, and a long-term dynamic river basin planning model with a large number of nonlinear relationships. For smaller instances of the reservoir model, the CONOPT2 nonlinear solver is more accurate and faster, but for larger instances, the GA&LP approach finds solutions with significantly better objective values. The multiperiod river basin model is much too large to be solved in its entirety. The complicating variables are chosen here so that, when they are fixed, each period's model is linear, and these models can be solved sequentially. This approach allows sufficient model detail to be retained so that long-term sustainability issues can be explored. 相似文献
13.
Multiphase flow in porous media is described by coupled nonlinear mass conservation laws. For immiscible Darcy flow of multiple fluid phases, whereby capillary effects are negligible, the transport equations in the presence of viscous and buoyancy forces are highly nonlinear and hyperbolic. Numerical simulation of multiphase flow processes in heterogeneous formations requires the development of discretization and solution schemes that are able to handle the complex nonlinear dynamics, especially of the saturation evolution, in a reliable and computationally efficient manner. In reservoir simulation practice, single-point upwinding of the flux across an interface between two control volumes (cells) is performed for each fluid phase, whereby the upstream direction is based on the gradient of the phase-potential (pressure plus gravity head). This upwinding scheme, which we refer to as Phase-Potential Upwinding (PPU), is combined with implicit (backward-Euler) time discretization to obtain a Fully Implicit Method (FIM). Even though FIM suffers from numerical dispersion effects, it is widely used in practice. This is because of its unconditional stability and because it yields conservative, monotone numerical solutions. However, FIM is not unconditionally convergent. The convergence difficulties are particularly pronounced when the different immiscible fluid phases switch between co-current and counter-current states as a function of time, or (Newton) iteration. Whether the multiphase flow across an interface (between two control-volumes) is co-current, or counter-current, depends on the local balance between the viscous and buoyancy forces, and how the balance evolves in time. The sensitivity of PPU to small changes in the (local) pressure distribution exacerbates the problem. The common strategy to deal with these difficulties is to cut the timestep and try again. Here, we propose a Hybrid-Upwinding (HU) scheme for the phase fluxes, then HU is combined with implicit time discretization to yield a fully implicit method. In the HU scheme, the phase flux is divided into two parts based on the driving force. The viscous-driven and buoyancy-driven phase fluxes are upwinded differently. Specifically, the viscous flux, which is always co-current, is upwinded based on the direction of the total-velocity. The buoyancy-driven flux across an interface is always counter-current and is upwinded such that the heavier fluid goes downward and the lighter fluid goes upward. We analyze the properties of the Implicit Hybrid Upwinding (IHU) scheme. It is shown that IHU is locally conservative and produces monotone, physically-consistent numerical solutions. The IHU solutions show numerical diffusion levels that are slightly higher than those for standard FIM (i.e., implicit PPU). The primary advantage of the IHU scheme is that the numerical overall-flux of a fluid phase remains continuous and differentiable as the flow regime changes between co-current and counter-current conditions. This is in contrast to the standard phase-potential upwinding scheme, in which the overall fractional-flow (flux) function is non-differentiable across the boundary between co-current and counter-current flows. 相似文献
14.
A two-dimensional (2D) numerical model has been developed to solve shallow water equations for simulation of dam-break flows. The spatial derivatives are discretized using a well-balanced explicit central upwind conservative scheme. The scheme is Riemann solver free and guarantees the positivity of the flow depth over complex topography if the Courant number is kept less than 0.25. The time integration is performed by Euler’s scheme. The model is verified against analytical results for water surface elevation and discharge for three benchmark test cases. A good agreement between analytical solutions and computed results is observed. The property of well-balancing in still water over an uneven bottom is also confirmed. The model is then validated by simulating a laboratory experiment in which a dam break flow propagates over a triangular obstacle. The model performance was found to be satisfactory. A dam break laboratory experimental test case on a frictionless horizontal bottom is also simulated for 2D validation of the model, and good agreement between simulation and the experimental data is observed. The suitability of the proposed model for real life applications is demonstrated by simulating the Malpasset dam-break event, which occurred in 1959 in France. The computed arrival time of the flood wave front and the maximum flow depths at various observation points matched well with the measurements on a 1/400 scale physical model. The overall performance indicates that this model can be applied for simulation of dam-break waves in real life cases. 相似文献
15.
A numerical method of the Godunov type is presented for solving either Two-Phase or Two-Layer forms of Debris Flow Models (DFMs) describing shallow-water flow and sediment dynamics. DFMs explicitly link sediment concentrations to the momentum balance, and thus can be applied to cases involving high sediment concentrations, as in debris flows, in addition to low concentration test cases typically found in surface waters. In this paper, Two-Phase and Two-Layer DFMs are presented in a common mathematical framework to illuminate key similarities and differences and lay a foundation for a general purpose DFM solver. The proposed solver termed LHLLC is shown to achieve good accuracy over a wide range of test cases. Importantly, numerical diffusion of sediment profiles is minimized, particularly on steep slopes, the scheme is shown to preserve stationary solutions involving wet/dry interfaces, and the scheme accounts for gravity-driven slumping (avalanching) which cannot be resolved by classical DFMs. 相似文献
16.
A generalized transformation approach for simulating steady-state variably-saturated subsurface flow
A numerical method is proposed to accurately and efficiently compute a direct steady-state solution of the nonlinear Richards equation. In the proposed method, the Kirchhoff integral transformation and a complementary transformation are applied to the governing equation in order to separate the nonlinear hyperbolic characteristic from the linear parabolic part. The separation allows the transformed governing equation to be applied to partially- to fully-saturated systems with arbitrary constitutive relations between primary (pressure head) and secondary variables (relative permeability). The transformed governing equation is then discretized with control volume finite difference/finite element approximations, followed by inverse transformation. The approach is compared to analytical and other numerical approaches for variably-saturated flow in 1-D and 3-D domains. The results clearly demonstrate that the approach is not only more computationally efficient but also more accurate than traditional numerical solutions. The approach is also applied to an example flow problem involving a regional-scale variably-saturated heterogeneous system, where the vadose zone is up to 1 km thick. The performance, stability, and effectiveness of the transform approach is exemplified for this complex heterogeneous example, which is typical of many problems encountered in the field. It is shown that computational performance can be enhanced by several orders of magnitude with the described integral transformation approach. 相似文献
17.
A correct understanding of the hydrodynamics and morphodynamics of tidal basins is of fundamental importance for the fate of the Venice Lagoon, Italy. If on one hand, the development of sophisticated numerical models is called for in order to reproduce the complexity of the mechanisms governing the morphodynamic evolution of many natural environments, including lagoons, on the other hand, a clear knowledge of the reliability and limits of the results provided by these models is crucial in order to establish the condition under which they can be safely applied. To this aim, researchers involved in numerical modeling in the framework of the recent Corila research programmes, agreed to perform an accurate comparison of results provided by three different numerical models, applying them to the test case offered by the experimental investigations performed under controlled conditions by Tambroni et al. (2005a). Here, we consider the following numerical models: (i) a 2D finite element hydrodynamic model coupled with a 2D finite volume morphodynamic model (5 and 3); (ii) a 2D finite element morphodynamic model (Ferrarin et al., 2008); (iii) a 2D depth-averaged model for the inlet region, coupled with a 1D model for the channel (Tambroni et al., 2005b). A first set of simulations concerns the fixed bed case and shows that all the models provide similar results: in particular, they are able to predict the observed free surface oscillations satisfactorily, while comparison with the measured velocity field is less satisfactory. Moreover, as far as the flow field at the inlet is concerned, the models describe accurately the potential flow into the channel during the flood phase, while they are not able to adequately reproduce the occurrence of the fine structure of the shear layers shed by the inlet edges during the ebb phase. This limit is related to the shallow water character of the models. As for the morphodynamics, the long term equilibrium configurations of the bottom of the channel and of the near inlet region show qualitative agreement with the experimental observations, although in this case the differences between the results provided by the distinct numerical approaches are more marked. 相似文献
18.
电磁法的三维数值模拟是一个对数值算法和计算机硬件要求都非常高的问题.对常用的微分类方法如有限单元法和有限差分法而言,求解最后所得的大型线性方程组是至关重要的一步,直接影响到正演算法的实用性.如何高效、稳定且准确地解线性方程长期以来一直是被探讨的问题.本文实现了基于线性系统直接求解技术的频率域可控源电磁(CSEM)三维正演.使用交错网格有限体积法(FV)来离散化关于二次电场的Helmholtz方程;使用直接解法取代传统的迭代解法来求解离散线性系统,即对系统矩阵进行完全LU分解,具体通过调用大规模并行矩阵直接求解器(MUMPS)来实现.基于理论模型做了一系列数值实验,首先证明了直接解法的高精度和稳定性,并考察了其内存需求、计算时间和并行可伸缩性等主要计算性能,最后检验了所开发的算法快速模拟多场源CSEM问题的能力以及对常规海洋和陆地CSEM模拟的有效性. 相似文献
19.
20.
首先以Lorenz混沌方程产生的非线性时间序列为例,讨论了在不同时间序列长度下各种延迟时间算法对噪声的适用性.研究发现,采用C_C算法计算延迟时间的鲁棒性强.在此基础上,给出了垂直上升管中气水两相流电导波动信号混沌表征结果,发现在较低水相表观速度时,随着气相表观速度增加,泡状流及混状流动力学特性变得愈加复杂,而段塞流动力学特性受液相表观速度影响较大;在较高水相表观速度时,随着气相表观速度增加,当流型从泡状流向段塞流转变时,气液两相流动力学特性变得相对简单.但是,由于受液相湍流作用影响,段塞流的动力学特性表现出了涨落现象,呈现不稳定性,当流型从段塞流向混状流转变时,气液两相流动力学特性则变得愈加复杂.研究结果表明：基于电导波动信号的混沌分析可以较好地表征气液两相流流型变化,是理解流型转变机理及其动力学演变特性的有用工具. 相似文献