共查询到20条相似文献，搜索用时 78 毫秒
1.
The axisymmetric formulation of the governing equations for geomechanics in the framework of smoothed particle hydrodynamics (SPH) is presented in this study. Two forms of SPH discretization for the motion equations, which are labeled as form I and form II, are proposed, and the methods to compute the hoop stress and strain terms including hoop strain rate and the acceleration introduced by the hoop stress are compared. To avoid possible singularity problem near the axis of symmetry, a perfectly smooth contact along with ghost particles are applied to prevent the real particles from overly approaching the axis of symmetry to remove this potential singularity. In addition, the Mohr–Coulomb constitutive model is implemented into the SPH formulation in describing soil behavior. Four numerical tests are carried out to validate and compare the accuracy and stability of the proposed algorithms, and their results are compared with analytical solutions and results from FEM analysis. The performance in these comparisons suggests that SPH II with hoop terms computed through direct hoop method is more stable than the others, and the adoption of contact for the symmetric axis is efficient in eliminating the singularity problem. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
2.
3.
We consider two-component (typically, water and hydrogen) compressible liquid–gas porous media flows including mass exchange between phases possibly leading to gas-phase (dis)appearance, as motivated by hydrogen production in underground repositories of radioactive waste. Following recent work by Bourgeat, Jurak, and Smaï, we formulate the governing equations in terms of liquid pressure and dissolved hydrogen density as main unknowns, leading mathematically to a nonlinear elliptic–parabolic system of partial differential equations, in which the equations degenerate when the gas phase disappears. We develop a discontinuous Galerkin method for space discretization, combined with a backward Euler scheme for time discretization and an incomplete Newton method for linearization. Numerical examples deal with gas-phase (dis)appearance, ill-prepared initial conditions, and heterogeneous problem with different rock types. 相似文献
4.
In this paper, we present a semi-implicit method for the incompressible three-phase flow equations in two dimensions. In particular, a high-order discontinuous Galerkin spatial discretization is coupled with a backward Euler discretization in time. We consider a pressure-saturation formulation, decouple the pressure and saturation equations, and solve them sequentially while still keeping each equation implicit in its respective unknown. We present several numerical examples on both homogeneous and heterogeneous media, with varying permeability and porosity. Our results demonstrate the robustness of the scheme. In particular, no slope limiters are required and a relatively large time step may be taken. 相似文献
5.
A Lagrangian particle‐based method, smooth particle hydrodynamics (SPH), is used in this paper to model the flow of self‐compacting concretes (SCC) with or without short steel fibres. An incompressible SPH method is presented to simulate the flow of such non‐Newtonian fluids whose behaviour is described by a Bingham‐type model, in which the kink in the shear stress vs shear strain rate diagram is first appropriately smoothed out. The viscosity of the SCC is predicted from the measured viscosity of the paste using micromechanical models in which the second phase aggregates are treated as rigid spheres and the short steel fibres as slender rigid bodies. The basic equations solved in the SPH are the incompressible mass conservation and Navier–Stokes equations. The solution procedure uses prediction–correction fractional steps with the temporal velocity field integrated forward in time without enforcing incompressibility in the prediction step. The resulting temporal velocity field is then implicitly projected on to a divergence‐free space to satisfy incompressibility through a pressure Poisson equation derived from an approximate pressure projection. The results of the numerical simulation are benchmarked against actual slump tests carried out in the laboratory. The numerical results are in excellent agreement with test results, thus demonstrating the capability of SPH and a proper rheological model to predict SCC flow and mould‐filling behaviour. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
6.
土体的大变形流滑导致了许多地质灾害的发生,对人们的生命财产安全构成了极大的威胁,因此越来越多的研究开始关注土体的大变形流滑特性。其中,光滑粒子流体动力学（SPH）方法是常用的模拟方法之一,但SPH方法的粒子特性导致其计算时间过长,影响了在工程地质领域的进一步应用。对此,本研究基于SPH方法的基本原理、非牛顿流体理论和等效黏度概念,提出了适用于土体大变形流滑分析的三维SPH仿真模型。结合OpenMP并行计算原理,实现了SPH算法的并行优化。在此基础上,对土体流滑模型试验进行了二维和三维分析,得到了滑动距离、滑动冲击力和冲击力峰值等动力学参数,分析了计算维数和边界条件对流滑特性的影响机制。通过不同线程数下计算时间的对比,获得了计算效率随线程数的变化规律。结果证明了本文的OpenMP并行优化具有较高的计算效率,显著降低了三维SPH模拟的计算耗时,对工程地质数值方法的效率提升具有重要的借鉴意义。 相似文献
7.
Coupled formulation and algorithms for the simulation of non‐planar three‐dimensional hydraulic fractures using the generalized finite element method 
下载免费PDF全文
下载免费PDF全文 This paper presents a coupled hydro‐mechanical formulation for the simulation of non‐planar three‐dimensional hydraulic fractures. Deformation in the rock is modeled using linear elasticity, and the lubrication theory is adopted for the fluid flow in the fracture. The governing equations of the fluid flow and elasticity and the subsequent discretization are fully coupled. A Generalized/eXtended Finite Element Method (G/XFEM) is adopted for the discretization of the coupled system of equations. A Newton–Raphson method is used to solve the resulting system of nonlinear equations. A discretization strategy for the fluid flow problem on non‐planar three‐dimensional surfaces and a computationally efficient strategy for handling time integration combined with mesh adaptivity are also presented. Several three‐dimensional numerical verification examples are solved. The examples illustrate the generality and accuracy of the proposed coupled formulation and discretization strategies. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
8.
考虑海冰热力因素对其厚度、密集度的影响,在光滑质点流体动力学(SPH)基础上发展了一个海冰热力-动力数值模式。该模式既解决了传统欧拉有限差分法和质点网格法存在的数值扩散问题,同时弥补了光滑质点动力学海冰动力模式未考虑热力因素的不足,具有精确模拟冰缘线运动、计算精度高等优点。首先介绍了光滑质点流体动力学的基本原理,并对海冰生消的热力因素进行了分析,将影响冰厚和密集度的热力因素引入到光滑质点流体动力学的海冰动力模式中,得到该热力-动力模式的控制方程。应用该数值模式对渤海海冰进行了48 h数值模拟,得到了海冰厚度和速度矢量的分布规律;对JZ20-2海域的海冰厚度、冰内温度场分布以及热力因素的变化特性进行了讨论。数值模拟结果表明,该数值模式能够很好地适用于渤海海冰数值模拟,是一种有效的海冰数值模拟方法。 相似文献
9.
A new SPH-based approach to simulation of granular flows using viscous damping and stress regularisation 总被引：2，自引：1，他引：1
Cuong T. Nguyen Chi T. Nguyen Ha H. Bui Giang D. Nguyen Ryoichi Fukagawa 《Landslides》2017,14(1):69-81
The smoothed particle hydrodynamics (SPH) method was recently extended to simulate granular materials by the authors and demonstrated to be a powerful continuum numerical method to deal with the post-flow behaviour of granular materials. However, most existing SPH simulations of granular flows suffer from significant stress oscillation during the post-failure process, despite the use of an artificial viscosity to damp out stress fluctuation. In this paper, a new SPH approach combining viscous damping with stress/strain regularisation is proposed for simulations of granular flows. It is shown that the proposed SPH algorithm can improve the overall accuracy of the SPH performance by accurately predicting the smooth stress distribution during the post-failure process. It can also effectively remove the stress oscillation issue in the standard SPH model without having to use the standard SPH artificial viscosity that requires unphysical parameters. The predictions by the proposed SPH approach show very good agreement with experimental and numerical results reported in the literature. This suggests that the proposed method could be considered as a promising continuum alternative for simulations of granular flows. 相似文献
10.
Run-out analysis of flow-like landslides triggered by the Ms 8.0 2008 Wenchuan earthquake using smoothed particle hydrodynamics 总被引：8，自引：5，他引：3
Flow-like landslides have caused significant damage and casualties worldwide. However, studying such phenomena with traditional simulation methods is made difficult by their complex fluidization characteristics. In this paper, we use smoothed-particle hydrodynamics (SPH) for the run-out analysis of flow-like landslides. Compared with conventional methods, the proposed SPH modeling technique is the combination of a Bingham flow model and Navier?CStokes equations in the framework of computational fluid dynamics. At first, two benchmark problems of dam break and granular flow are simulated and verified to evaluate the accuracy of the SPH model. Then, run-out analyses are performed for flow-like landslides triggered by the Ms 8.0 Wenchuan earthquake that occurred on 12 May 2008 in Sichuan Province, China. Run-out analyses of the Tangjiashan, Wangjiayan, and Donghekou landslides are conducted by the application of SPH models to real flow-like landslides. All simulations show good agreement with characteristics of flow-like landslides observed in the field. We have found that numerical modeling can capture the fundamental dynamic behavior of these flow-like landslides and produce preliminary results for hazard assessment and site selection for reconstruction in earthquake-prone areas. 相似文献
11.
12.
Simulation of large deformation and post‐failure of geomaterial in the framework of smoothed particle hydrodynamics (SPH) are presented in this study. The Drucker–Prager model with associated and non‐associated plastic flow rules is implemented into the SPH code to describe elastic–plastic soil behavior. In contrast to previous work on SPH for solids, where the hydrostatic pressure is often estimated from density by an equation of state, this study proposes to calculate the hydrostatic pressure of soil directly from constitutive models. Results obtained in this paper show that the original SPH method, which has been successfully applied to a vast range of problems, is unable to directly solve elastic–plastic flows of soil because of the so‐called SPH tensile instability. This numerical instability may result in unrealistic fracture and particles clustering in SPH simulation. For non‐cohesive soil, the instability is not serious and can be completely removed by using a tension cracking treatment from soil constitutive model and thereby give realistic soil behavior. However, the serious tensile instability that is found in SPH application for cohesive soil requires a special treatment to overcome this problem. In this paper, an artificial stress method is applied to remove the SPH numerical instability in cohesive soil. A number of numerical tests are carried out to check the capability of SPH in the current application. Numerical results are then compared with experimental and finite element method solutions. The good agreement obtained from these comparisons suggests that SPH can be extended to general geotechnical problems. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
13.
The method of smoothed particle hydrodynamics (SPH) has recently been applied to computational geomechanics and has been shown to be a powerful alternative to the standard numerical method, that is, the finite element method, for handling large deformation and post‐failure of geomaterials. However, very few studies apply the SPH method to model saturated or submerged soil problems. Our recent studies of this matter revealed that significant errors may be made if the gradient of the pore‐water pressure is handled using the standard SPH formulation. To overcome this problem and to enhance the SPH applications to computational geomechanics, this article proposes a general SPH formulation, which can be applied straightforwardly to dry and saturated soils. For simplicity, the current work assumes hydrostatic pore‐water pressure. It is shown that the proposed formulation can remove the numerical error mentioned earlier. Moreover, this formulation automatically satisfies the dynamic boundary conditions at a submerged ground surface, thereby saving computational cost. Discussions on the applications of the standard and new SPH formulations are also given through some numerical tests. Furthermore, techniques to obtain the correct SPH solution are also proposed and discussed throughout. As an application of the proposed method, the effect of the dilatancy angle on the failure mechanism of a two‐sided embankment subjected to a high groundwater table is presented and compared with that of other solutions. Finally, the proposed formulation can be considered a basic formulation for further developments of SPH for saturated soils. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
14.
15.
In this paper, we present a computational framework for the simulation of coupled flow and reservoir geomechanics. The physical
model is restricted to Biot’s theory of single-phase flow and linear poroelasticity, but is sufficiently general to be extended
to multiphase flow problems and inelastic behavior. The distinctive technical aspects of our approach are: (1) the space discretization
of the equations. The unknown variables are the pressure, the fluid velocity, and the rock displacements. We recognize that
these variables are of very different nature, and need to be discretized differently. We propose a mixed finite element space
discretization, which is stable, convergent, locally mass conservative, and employs a single computational grid. To ensure
stability and robustness, we perform an implicit time integration of the fluid flow equations. (2) The strategies for the
solution of the coupled system. We compare different solution strategies, including the fully coupled approach, the usual
(conditionally stable) iteratively coupled approach, and a less common unconditionally stable sequential scheme. We show that
the latter scheme corresponds to a modified block Jacobi method, which also enjoys improved convergence properties. This computational
model has been implemented in an object-oriented reservoir simulator, whose modular design allows for further extensions and
enhancements. We show several representative numerical simulations that illustrate the effectiveness of the approach. 相似文献
16.
Liquefaction can result in the damage or collapse of structures during an earthquake and can therefore be a great threat to life and property. Many site investigations of liquefaction disasters are needed to study the large-scale deformation and flow mechanisms of liquefied soils that can be used for performance assessments and infrastructure improvement. To overcome the disadvantages of traditional flow analysis methods for liquefied soils, a soil–water-coupled smoothed particle hydrodynamics (SPH) modeling method was developed to analyze flow in liquefied soils. In the proposed SPH method, water and soil were simulated as different layers, while permeability, porosity, and interaction forces could be combined to model water-saturated porous media. A simple shear test was simulated using the SPH method with an elastic model to verify its application to solid phase materials. Subsequently, the applicability of the proposed SPH modeling method to the simulation of interaction forces between water and soil was verified by a falling-head permeability test. The coupled SPH method produced good simulations for both the simple shear and falling-head permeability tests. Using a fit-for-purpose experimental apparatus, a physical flow model test of liquefied sand has been designed and conducted. To complement the physical test, a numerical simulation has been undertaken based on the soil–water-coupled SPH method. The numerical results correspond well with the physical model test results in observed configurations and velocity vectors. An embankment failure in northern Sweden was selected so that the application of the soil–water-coupled SPH method could be extended to an actual example of liquefaction. The coupled SPH method simulated the embankment failure with the site investigation well. They have also estimated horizontal displacements and velocities, which can be used to greatly improve the seismic safety of structures. 相似文献
17.
Smoothed particle hydrodynamics (SPH) is a meshfree, Lagrangian particle method which has advantages in handling solids with extremely large deformation. Like any other numerical methods, cares must be taken to ensure its desirable accuracy and stability through considering several correction techniques in calculation. The selection of values for parameters in those correction approaches is a key step in SPH simulation, which is always difficult for new beginners to deal well with effectively. This paper examines the common inconsistency and instability problems in SPH method and studies its computational efficiency when applied to hydrodynamics problems with material strength like soil column collapse. We analyzed in detail how the correction techniques mitigate these inconsistency and instability problems. Also, the numerical testing results associate with different values for the parameters used in the correction techniques are provided for better understanding the influence of these parameters and for finding out the desirable values. It is found that (1) the SPH method is easily subjected to an inconsistency problem in the boundary area due to the boundary deficiency, and it can be treated well by adopting “virtual particles” contributing to the particle summations. (2) The numerical oscillation in SPH simulation can be mitigated effectively by artificial viscosity with the suggested parameter values. (3) The tension cracking treatment, artificial viscosity and artificial stress work well in removing the tensile instability problem in SPH method. In addition, the nearest neighboring particle searching (NNPS) algorithm, spacing ratio, smoothing length and time step influence the efficiency and accuracy of SPH method significantly. It is shown that SPH method with suggested parameters values can produce a very good result compared with the experimental result. 相似文献
18.
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. 相似文献
19.