共查询到20条相似文献，搜索用时 468 毫秒
1.
We present a high-order method for miscible displacement simulation in porous media. The method is based on discontinuous Galerkin discretization with weighted average stabilization technique and flux reconstruction post processing. The mathematical model is decoupled and solved sequentially. We apply domain decomposition and algebraic multigrid preconditioner for the linear system resulting from the high-order discretization. The accuracy and robustness of the method are demonstrated in the convergence study with analytical solutions and heterogeneous porous media, respectively. We also investigate the effect of grid orientation and anisotropic permeability using high-order discontinuous Galerkin method in contrast with cell-centered finite volume method. The study of the parallel implementation shows the scalability and efficiency of the method on parallel architecture. We also verify the simulation result on highly heterogeneous permeability field from the SPE10 model. 相似文献
2.
Joachim Berdal Haga Harald Osnes Hans Petter Langtangen 《Computational Geosciences》2012,16(3):723-734
Large-scale simulations of coupled flow in deformable porous media require iterative methods for solving the systems of linear
algebraic equations. Construction of efficient iterative methods is particularly challenging in problems with large jumps
in material properties, which is often the case in realistic geological applications, such as basin evolution at regional
scales. The success of iterative methods for such problems depends strongly on finding effective preconditioners with good
parallel scaling properties, which is the topic of the present paper. We present a parallel preconditioner for Biot’s equations
of coupled elasticity and fluid flow in porous media. The preconditioner is based on an approximation of the exact inverse
of the two-by-two block system arising from a finite element discretisation. The approximation relies on a highly scalable
approximation of the global Schur complement of the coefficient matrix, combined with generally available state-of-the-art
multilevel preconditioners for the individual blocks. This preconditioner is shown to be robust on problems with highly heterogeneous
material parameters. We investigate the weak and strong parallel scaling of this preconditioner on up to 512 processors and
demonstrate its ability on a realistic basin-scale problem in poroelasticity with over eight million tetrahedral elements. 相似文献
3.
Graphic Processing Unit (GPU), as a computing device, has upgraded from single-subject graphical processors to multi-core processors with tremendous computational horsepower. This paper proposes to accelerate the DDA using parallel Jacobi Preconditioned Conjugate Gradient (JPCG) technique on GPUs. Based on the results of two numerical examples, the calculation accuracies of the DDA with serial and parallel solvers are validated, and we found that the DDA with parallel solvers exhibits a much higher execution efficiency. The movement process of Daguangbao landslide triggered by the Wenchuan earthquake is replicated and the modeled deposit pattern coincides well with the actual topography after earthquake. 相似文献
4.
Multiscale methods can in many cases be viewed as special types of domain decomposition preconditioners. The localisation
approximations introduced within the multiscale framework are dependent upon both the heterogeneity of the reservoir and the
structure of the computational grid. While previous works on multiscale control volume methods have focused on heterogeneous
elliptic problems on regular Cartesian grids, we have tested the multiscale control volume formulations on two-dimensional
elliptic problems involving heterogeneous media and irregular grid structures. Our study shows that the tangential flow approximation
commonly used within multiscale methods is not suited for problems involving rough grids. We present a more robust mass conservative
domain decomposition preconditioner for simulating flow in heterogeneous porous media on general grids. 相似文献
5.
Soil–structure interaction problems are commonly encountered in engineering practice, and the resulting linear systems of equations are difficult to solve due to the significant material stiffness contrast. In this study, a novel partitioned block preconditioner in conjunction with the Krylov subspace iterative method symmetric quasiminimal residual is proposed to solve such linear equations. The performance of these investigated preconditioners is evaluated and compared on both the CPU architecture and the hybrid CPU–graphics processing units (GPU) computing environment. On the hybrid CPU–GPU computing platform, the capability of GPU in parallel implementation and high-intensity floating point operations is exploited to accelerate the iterative solutions, and particular attention is paid to the matrix–vector multiplications involved in the iterative process. Based on a pile-group foundation example and a tunneling example, numerical results show that the partitioned block preconditioners investigated are very efficient for the soil–structure interaction problems. However, their comparative performances may apparently depend on the computer architecture. When the CPU computer architecture is used, the novel partitioned block symmetric successive over-relaxation preconditioner appears to be the most efficient, but when the hybrid CPU–GPU computer architecture is adopted, it is shown that the inexact block diagonal preconditioners embedded with simple diagonal approximation to the soil block outperform the others. 相似文献
6.
Smoothed particle hydrodynamics (SPH) is a versatile technique which can be applied to single and multiphase flow through porous media. The versatility of SPH is offset by its computational expense which limits the practicability of SPH for large problems involving low Reynolds number flow. A parallel pore-scale numerical model based on SPH is described for modeling flow phenomena in porous media. Aspects of SPH which complicate parallelization are emphasized. The speed of the method is demonstrated to be proportional to the number of processors for test cases where load balance was achieved. The parallel algorithm permits the application of SPH to more complicated porous media problems than previously considered. For such problems, best performance is achieved when several soil grains are simulated by each processor. Finally, future applications of the method and possible extensions are discussed. 相似文献
7.
Ettore Vidotto Rainer Helmig Martin Schneider Barbara Wohlmuth 《Computational Geosciences》2018,22(6):1487-1502
In this paper, we present a fast streamline-based numerical method for the two-phase flow equations in high-rate flooding scenarios for incompressible fluids in heterogeneous and anisotropic porous media. A fractional flow formulation is adopted and a discontinuous Galerkin method (DG) is employed to solve the pressure equation. Capillary effects can be neglected in high-rate flooding scenarios. This allows us to present an improved streamline approach in combination with the one-dimensional front tracking method to solve the transport equation. To handle the high computational costs of the DG approximation, domain decomposition is applied combined with an algebraic multigrid preconditioner to solve the linear system. Special care at the interior interfaces is required and the streamline tracer has to include a dynamic communication strategy. The method is validated in various two- and three-dimensional tests, where comparisons of the solutions in terms of approximation of flow front propagation with standard fully implicit finite-volume methods are provided. 相似文献
8.
We describe a new approach for simulation of multiphase flows through heterogeneous porous media, such as oil reservoirs.
The method, which is based on the wavelet transformation of the spatial distribution of the single-phase permeabilities, incorporates
in the upscaled computational grid all the relevant data on the permeability, porosity, and other important properties of
a porous medium at all the length scales. The upscaling method generates a nonuniform computational grid which preserves the resolved structure
of the geological model in the near-well zones as well as in the high-permeability sectors and upscales the rest of the geological
model. As such, the method is a multiscale one that preserves all the important information across all the relevant length
scales. Using a robust front-detection method which eliminates the numerical dispersion by a high-order total variation diminishing
method (suitable for the type of nonuniform upscaled grid that we generate), we obtain highly accurate results with a greatly
reduced computational cost. The speed-up in the computations is up to over three orders of magnitude, depending on the degree
of heterogeneity of the model. To demonstrate the accuracy and efficiency of our methods, five distinct models (including
one with fractures) of heterogeneous porous media are considered, and two-phase flows in the models are studied, with and
without the capillary pressure. 相似文献
9.
The problem of finite element simulation of incompressible fluid flow in porous medium is considered. The porous medium is characterized by the X‐ray microtomography technique in three dimensions. The finite calculus‐based stabilization technique is reviewed to implement the equal order finite element interpolation functions for both velocity and pressure. A noble preconditioner, the nodal block diagonal preconditioner, is considered whose performance is thoroughly investigated. Combining this preconditioner with a standard iterative solver during the computational homogenization procedure, it is possible to carry out the large‐scale fluid flow simulation for estimating permeability of the porous medium with reasonable accuracy and reliability. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
10.
The material point method (MPM), which is a combination of the finite element and meshfree methods, suffers from significant computational workload due to the fine mesh that is required in spite of its advantages in simulating large deformations. This paper presents a parallel computing strategy for the MPM on the graphics processing unit (GPU) to boost the method’s computational efficiency. The interaction between a structural element and soil is investigated to validate the applicability of the parallelisation strategy. Two techniques are developed to parallelise the interpolation from soil particles to nodes to avoid a data race; the technique that is based on workload parallelisation across threads over the nodes has a higher computational efficiency. Benchmark problems of surface footing penetration and a submarine landslide are analysed to quantify the speedup of GPU parallel computing over sequential simulations on the central processing unit. The maximum speedup with the GPU used is ∼30 for single-precision calculations and decreases to ∼20 for double-precision calculations. 相似文献
11.
In this paper, a multiscale homogenization approach is developed for fully coupled saturated porous media to represent the idealized sugar cube model, which is generally employed in fractured porous media on the basis of dual porosity models. In this manner, an extended version of the Hill-Mandel theory that incorporates the microdynamic effects into the multiscale analysis is presented, and the concept of the deformable dual porosity model is demonstrated. Numerical simulations are performed employing the multiscale analysis and dual porosity model, and the results are compared with the direct numerical simulation through 2 numerical examples. Finally, a combined multiscale-dual porosity technique is introduced by employing a bridge between these 2 techniques as an alternative approach that reduces the computational cost of numerical simulation in modeling of heterogeneous deformable porous media. 相似文献
12.
Erik O. D. Sevre Monica D. Christiansen Matt Broten Shuo M. Wang David A. Yuen 《Visual Geosciences》2008,13(1):125-132
The Sony PlayStation 3 (PS3) offers the computational power of a parallel processor at low cost, which makes it a great starter
unit for development in parallel programming. To explore the capabilities of the unit, we took a simple ray tracing program
and extended it to render triangulated height field data across the PS3’s 6 synergistic processing units (SPUs). We also implemented
the heat averaging equation as a precursor to CFD analysis on the PS3. In our studies, we found the Cell engine in the PS3
to be a powerful machine, however great care must be taken while developing because its unique platform calls for many levels
of optimization to ensure efficiency.
For insight on early work done on modern stream computing and GPU processing look up the Merrimac project and Pat Hanrahan. 相似文献
13.
JHomogenizer: a computational tool for upscaling permeability for flow in heterogeneous porous media
This paper presents an object-oriented programming approach for the design of numerical homogenization programs, called JHomogenizer. It currently includes five functional modules to compute effective permeability and simple codes for computing solutions for flow in porous media. Examples with graphical output are shown to illustrate some functionalities of the program. A series of numerical examples demonstrates the effectiveness of the methodology for two-phase flow in heterogeneous reservoirs. The software is freely available, and the open architecture of the program facilitates further development and can adapt to suit specific needs easily and quickly. 相似文献
14.
Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method 
下载免费PDF全文
下载免费PDF全文 Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
15.
F. O. Alpak F. Gray N. Saxena J. Dietderich R. Hofmann S. Berg 《Computational Geosciences》2018,22(3):815-832
Digital rock physics (DRP) is a rapidly evolving technology targeting fast turnaround times for repeatable core analysis and multi-physics simulation of rock properties. We develop and validate a rapid and scalable distributed-parallel single-phase pore-scale flow simulator for permeability estimation on real 3D pore-scale micro-CT images using a novel variant of the lattice Boltzmann method (LBM). The LBM code implementation is designed to take maximum advantage of distributed computing on multiple general-purpose graphics processing units (GPGPUs). We describe and extensively test the distributed parallel implementation of an innovative LBM algorithm for simulating flow in pore-scale media based on the multiple-relaxation-time (MRT) model that utilizes a precise treatment of body force. While the individual components of the resulting simulator can be separately found in various references, our novel contributions are (1) the integration of all of the mathematical and high-performance computing components together with a highly optimized code implementation and (2) the delivery of quantitative results with the simulator in terms of robustness, accuracy, and computational efficiency for a variety of flow geometries including various types of real rock images. We report on extensive validations of the simulator in terms of accuracy and provide near-ideal distributed parallel scalability results on large pore-scale image volumes that were largely computationally inaccessible prior to our implementation. We validate the accuracy of the MRT-LBM simulator on model geometries with analytical solutions. Permeability estimation results are then provided on large 3D binary microstructures including a sphere pack and rocks from various sandstone and carbonate formations. We quantify the scalability behavior of the distributed parallel implementation of MRT-LBM as a function of model type/size and the number of utilized GPGPUs for a panoply of permeability estimation problems. 相似文献
16.
深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体(分辨率3.7μm)为研究对象,在空间节点数多达64 000 000的情况下,基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数,并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现,结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD(即10^-13m^2)量级,远低于规则堆叠的圆球结构(>20 000 mD);土壤的迂曲度为1.40~1.60,明显高于规则堆叠的圆球结构。研究结果表明,渗透率大的样品具有较小的迂曲度,这与结构的空间异质性有较强的关系;土壤的渗透率和迂曲度呈现各向异性;在水力梯度一定的前提下,渗透率较大的样品,纵向弥散系数也较大;同时,结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟,可用于土壤多孔介质在孔隙尺度下的水力学特性研究。 相似文献
17.
深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体（分辨率3.7 μm）为研究对象，在空间节点数多达64 000 000的情况下，基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数，并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现，结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD（即10-13m2）量级，远低于规则堆叠的圆球结构（>20 000 mD）；土壤的迂曲度为1.40~1.60，明显高于规则堆叠的圆球结构。研究结果表明，渗透率大的样品具有较小的迂曲度，这与结构的空间异质性有较强的关系；土壤的渗透率和迂曲度呈现各向异性；在水力梯度一定的前提下，渗透率较大的样品，纵向弥散系数也较大；同时，结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟，可用于土壤多孔介质在孔隙尺度下的水力学特性研究。 相似文献
18.
Alexandr A. Nikitin Alexandr S. Serdyukov Anton A. Duchkov 《Computational Geosciences》2018,22(3):775-787
Numerical solution of the eikonal equation is frequently used to compute first-arrival travel times for a given velocity model in seismic applications. Computations for large three-dimensional models become expensive requiring the use of efficient parallel solvers. We present new parallel implementations of the fast sweeping and locking sweeping methods optimized for shared memory systems such as multicore CPUs; we call them block fast sweeping method (BFSM) and block locking sweeping method (BLSM). Proposed methods are based on the domain decomposition approach with a special attention paid to high efficiency of the cache utilization and task execution synchronization. Performance tests on realistic models show high parallel efficiency of 85–95% on modern multicore CPUs and require the same number of iterations to converge as do the serial sweeping methods. We also highlight the importance of properly selecting the stopping criterion in the iterative sweeping methods aiming for a balance between computational time and accuracy of the result required by an application. In particular, we show that in seismic applications one can reach reasonable accuracy of computed travel times while dramatically reducing the number of iterations compared to the case of using the full convergence stopping criterion. 相似文献
19.
Numerical simulation of gas migration driven by compressible two-phase partially miscible flow in porous media is of major importance for safety assessment of deep geological repositories for long-lived high-level nuclear waste. We present modeling of compositional liquid and gas flow for numerical simulations of hydrogen migration in deep geological radioactive waste repository based on persistent primary variables. Two-phase flow is considered, with incompressible liquid and compressible gas, which includes capillary effects, gas dissolution, and diffusivity. After discussing briefly the existing approaches to deal with phase appearance and disappearance problem, including a persistent set of variables already considered in a previous paper (Bourgeat et al., Comput Geosci 13(1):29–42, 2009), we focus on a new variant of the primary variables: dissolved hydrogen mass concentration and liquid pressure. This choice leads to a unique and consistent formulation in liquid saturated and unsaturated regions, which is well adapted to heterogeneous media. We use this new set of variable for numerical simulations and show computational evidences of its adequacy to simulate gas phase appearance and disappearance in different but typical situations for gas migration in an underground radioactive waste repository. 相似文献
20.
A three-dimensional, three-phase numerical model is presented for simulating the movement of immiscible fluids, including
nonaqueous-phase liquids (NAPLs), through porous media. The model is designed to simulate soil flume experiments and for practical
application to a wide variety of contamination scenarios involving light or dense NAPLs in heterogeneous subsurface systems.
The model is derived for the three-phase flow of water, NAPL, and air in porous media. The basic governing equations are based
upon the mass conservation of the constitutents within the phases. The descretization chosen to transform the governing equations
into the approximating equations, although logically regular, is very general. The approximating equations are a set of simultaneous
coupled nonlinear equations which are solved by the Newton-Raphson method. The linear system solutions needed for the Newton-Raphson
method are obtained using a matrix of preconditioner/accelerator iterative methods.
Because of the special way the governing equations are implemented, the model is capable of simulating many of the phenomena
considered necessary for the sucessful simulation of field problems including entry pressure phenomena, entrapment, and preferential
flow paths. The model is verified by comparing it with several exact analytic test solutions and three soil flume experiments
involving the introduction and movement of light nonaqueous-phase liquid (LNAPL) or dense nonaqueous-phase liquid (DNAPL)
in heterogeneous sand containing a watertable.
This revised version was published online in August 2006 with corrections to the Cover Date. 相似文献