首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The flow of water–kaolinite mixtures exhibits a non‐Newtonian nature that differs from the flow of Newtonian fluid. The varying viscosities and shear history of non‐Newtonian fluid flows necessitate the use of a rheology model in moving particle semi‐implicit (MPS) for the numerical studies. On the other hand, the Lagrangian method has the advantage of handling free surface flows with large deformation and fragmentation. This study proposes a mesh‐free Lagrangian method, namely, the MPS method, together with a simple rheology model to investigate the non‐Newtonian free surface flows. The rheological parameters required in the rheology model are determined based upon experiments. The proposed method is applied to a water–kaolinite mixture collapse problem and is proved to be capable of reproducing the significant flow features observed in laboratory experiments. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

2.
Characteristics of fluid flow through various engineering structures, such as granular filters and asphalt pavements, influence their design life. Numerical simulation of fluid flow is useful for evaluating the hydraulic characteristics of these materials. Among various techniques, the lattice Boltzmann (LB) method is widely accepted due to the ease of implementing boundary conditions and the numerical stability in a wide variety of flow conditions. It has proven to be extremely efficient in the simulation of fluid flow through the complex geometries of granular materials. In this study, two-dimensional and three-dimensional LB models were developed to represent pore-scale monophasic Newtonian incompressible fluid flow in granular materials. Three-dimensional geometries of compacted aggregates and asphalt specimens were generated from X-ray Computed Tomography technique and used as input for the LB model. The accuracy of the models was verified by comparing the results with analytical solutions of simple geometries and hydraulic conductivity measurements on the compacted aggregates and hot mix asphalt specimens. The results of LB simulations were in excellent agreement with those obtained from analytical calculations and laboratory measurements.  相似文献   

3.
A fully coupled transient two‐dimensional model was employed to study fundamentals of flood‐induced surface erosion in a particle bed. The interaction of the liquid and solid phases is the key mechanism related to surface erosion. The solid phase was idealized at a particle scale by using the discrete element method. The fluid phase was modeled at a mesoscale level and solved using the lattice Boltzmann method. The fluid forces applied on the particles were calculated on the basis of the momentum the fluid exchanges with the particle. The proposed approach was used to model both single particles and particle beds subjected to Couette flow conditions. The behavior of both the single particle and the particle bed depended on particle diameter and surface shear fluid velocity. The conducted simulations show that the fluid flow profile penetrates the bed for a small distance. This penetration initiates sheet‐flow and surface erosion as the fluid interacts with particles. The effect of suppressing particle rotation on the fluid‐induced forces on the particle was also examined. Suppressing particle spinning may lead to underestimated erosion rate. Results of fluid and particle velocities were compared against experimental results and appeared to agree with the observed trends.Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

4.
This paper presents a 3D bonded discrete element and lattice Boltzmann method for resolving the fluid‐solid interaction involving complicated fluid‐particle coupling in geomaterials. In the coupled technique, the solid material is treated as an assembly of bonded and/or granular particles. A bond model accounting for strain softening in normal contact is incorporated into the discrete element method to simulate the mechanical behaviour of geomaterials, whilst the fluid flow is solved by the lattice Boltzmann method based on kinetic theory and statistical mechanics. To provide a bridge between theory and application, a 3D algorithm of immersed moving boundary scheme was proposed for resolving fluid‐particle interaction. To demonstrate the applicability and accuracy of this coupled method, a benchmark called quicksand, in which particles become fluidised under the driving of upward fluid flow, is first carried out. The critical hydraulic gradient obtained from the numerical results matches the theoretical value. Then, numerical investigation of the performance of granular filters generated according to the well‐acknowledged design criteria is given. It is found that the proposed 3D technique is promising, and the instantaneous migration of the protected soils can be readily observed. Numerical results prove that the filters which comply with the design criteria can effectively alleviate or eliminate the appearance of particle erosion in dams.  相似文献   

5.
泥石流冲毁桥墩是桥梁在遭受泥石流冲击时的常见破坏形式。为了研究泥石流对桥墩的冲击力大小,通过调整黏土、沙、石子、水的不同含量,配置不同流变特性、不同密度的泥石流,使用所配置的原料在泥石流槽内对两种形状（圆形、方形）的桥墩缩尺模型进行冲击,综合考察了流变特性、流速、桥墩形状以及冲击力的关系。试验表明：试验配置的泥石流原料流变特性差异明显,且可以用简单的选择流变仪测得,用牛顿流体或宾汉体描述。泥石流的流速可用曼宁公式求得,而公式中的糙率系数与泥石流黏度满足幂函数关系。相同工况下,不同形状桥墩所受的冲击力差异明显,方形桥墩阻力系数普遍大于圆形桥墩。使用非牛顿流体雷诺数（Re）可以综合反映流变特性和流速,因此,圆墩的阻力系数可表达为Re的函数,而方墩则没有明显关系。为方便工程应用,可根据黏性泥石流、稀性泥石流对圆墩的阻力系数分别为2.3、0.9,对方墩分别为2.6、1.9进行选用。  相似文献   

6.
Flow slides in municipal solid waste (MSW) dumps have caused serious damage to structures and casualties all over the world. Therefore, much attention should be paid to this type of disaster to elucidate the flow mechanisms and fluidization characteristics of MSW, which are essential for the assessment and prevention of flowlike hazards. To bypass the deficiencies of the traditional analysis methods that use the mesh method and are based on a framework of solid mechanics, the moving particle semi-implicit (MPS) method, which is a purely Lagrangian meshless method and proposed for incompressible flow, is introduced to study flow slides in MSW landfills. Considering the no-physical pressure fluctuation that affects the simulation accuracy in the original MPS, the original MPS is revised in three ways: the kernel function, the source term of the Poisson equation and the search for free surface particles. Two benchmark problems, the dam break problem and the static pressure problem, are computed to illustrate the improvement of the pressure stability of the modified MPS. The Bingham constitutive model combined with the Mohr–Coulomb failure criterion is adopted to depict the dynamic features of MSW flow slides, and the equivalent viscosity is employed to bridge the gap between Bingham fluid models and Newtonian fluid models. This method, ultimately, is applied to simulate real flow slides in the Umraniye–Hekimbashi landfill and the Payatas waste dump. The numerical results show good consistency with the field data, indicating that the modified MPS method is capable of capturing the essential dynamic behavior and reproducing the entire process of complicated flow slides in MSW dumps.  相似文献   

7.
In this article, we present a numerical method to deal with fluid–solid interactions and simulate particle–fluid systems as encountered in soils. This method is based on a coupling between two methods, now widely used in mechanics of granular media and fluid dynamics respectively: the discrete element (DE) method and the lattice Boltzmann (LB) method. The DE method is employed to model interactions between particles, whereas the LB method is used to describe an interstitial Newtonian fluid flow. The coupling presented here is a full one in the sense that particle motions act on fluid flow and reciprocally. This article presents in details each of the two methods and the principle of the coupling scheme. Determination of hydrodynamic forces and torques is also detailed, and the treatment of boundaries is explained. The coupled method is finally illustrated on a simple example of piping erosion, which puts in evidence that the combined LB–DE scheme constitutes a promising tool to study coupled problems in geomechanics. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

8.
用格子波耳兹曼方法模拟双重孔隙介质中的流体迁移   总被引：3，自引：0，他引：3  
作者在本文中介绍了基于格子波耳兹曼模型的双重孔隙介质中流体运移的数值模拟计算方法。我们从格子波耳兹曼碰撞模型出发，利用格子波耳兹曼方程、Chapman-Enskog展开，以及多尺度技术，得到了描述双重孔隙介质中流体迁移的二维扩散方程。利用格子气自动机方法计算该扩散方程，实现了对双重孔隙介质中流体运移过程的数值模拟仿真。数值实验表明，我们所使用的方法正确、有效。  相似文献   

9.
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.  相似文献   

10.
This paper describes the flow behaviour of certain non-Newtonian fluids through a porous medium. A generalized Bingham rheological model of power-law in the presence of a yield stress has been considered. Several problems of fluid mechanics, which appear currently in oil reservoir engineering, have been investigated and the rheological behaviour effect has been emphasized. The short time solutions have been formulated in terms of a moving boundary problem. The approximate solutions in a closed formwere obtained by means of the integral method. Several dimensionless groups have been found to be relevant in evaluating the rheological effect on the steady and unsteady and unsteady flow behaviour. The deviation from Newtonian flow behaviour has been illustrated using several numerical examples.  相似文献   

11.
The coupled discrete element method and lattice Boltzmann method (DEMLBM) has increasingly drawn attention of researchers in geomechanics due to its mesoscopic nature since 2000. Immersed boundary method (IBM) and immersed moving boundary (IMB) are two popular schemes for coupling fluid particle in DEMLBM. This work aims at coupling DEM and LBM using the latest IBM algorithm and investigating its accuracy, computational efficiency, and applicability. Two benchmark tests, interstitial fluid flow in an ideal packing and single particle sedimentation in viscous fluid, are carried out to demonstrate the accuracy of IBM through semi-empirical Ergun equation, finite element method (FEM), and IMB. Then, simulations of particle migration with relatively large velocity in Poiseuille flow are utilized to address limitations of IBM in DEMLBM modeling. In addition, advantages and deficiencies of IBM are discussed and compared with IMB. It is found that the accuracy of IBM can be only guaranteed when sufficient boundary points are used and it is not suitable for geomechanical problems involving large fluid or particle velocity.  相似文献   

12.
Three porous media flow problems, in which the fluid mechanical interactions are critical, are studied in a mesoscopic–microscopic coupling system. In this system, fluid flow in the pore space is explicitly modeled at mesoscopic level by the lattice Boltzmann method, the geometrical representation and the mechanical behavior of the solid skeleton are modeled at microscopic level by the particulate distinct element method (DEM), and the interfacial interaction between the fluid and the solids is resolved by an immersed boundary scheme. In the first benchmark problem, the well‐known and frequently utilized Ergun equation is validated in periodic particle and periodic pore models. In the second problem, the upward seepage problem is simulated over three stages: The settlement of the column of sphere under gravity loading is measured to illustrate the accuracy of the DEM scheme; the system is solved to hydrostatic state with pore space filled with fluid, showing that the buoyancy effect is captured correctly in the mesoscopic–microscopic coupling system; then, the flow with constant rate is supplied at the bottom of the column; the swelling of the ground surface and pore pressure development from the numerical simulation are compared with the predictions of the macroscopic consolidation theory. In the third problem, the fluid‐flow‐induced collapse of a sand arch inside a perforation cavity is tested to illustrate a more practical application of the developed system. Through comparing simulation results with analytical solutions, empirical law and physical laboratory observations, it is demonstrated that the developed lattice Boltzmann–distinct element coupling system is a powerful fundamental research tool for investigating hydromechanical physics in porous media flow. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

13.
Interaction between solid particles and fluid is of fundamental interest to scientists and engineers in many different applications—cardiopulmonary flows, aircraft and automobile aerodynamics, and wind loading on buildings to name a few. In geomechanics, particle shape significantly affects both particle-particle and particle-fluid interaction. Herein, we present a generalized method for modeling the interaction of arbitrarily shaped polyhedral particles and particle assemblages with fluid using a coupled discrete element method (DEM) and lattice Boltzmann method (LBM) formulation. The coupling between DEM and LBM is achieved through a new algorithm based on a volume-fraction approach to consider three-dimensional convex polyhedral particles moving through fluid. The algorithm establishes the interaction using linear programming and simplex integration and is validated against experimental data. This approach to modeling the interaction between complex polyhedral particles and fluid is shown to be accurate for directly simulating hydrodynamic forces on the particles.  相似文献   

14.
To investigate the movement mechanism of debris flow, a two‐dimensional, two‐phase, depth‐integrated model is introduced. The model uses Mohr‐Coulomb plasticity for the solid rheology, and the fluid stress is modeled as a Newtonian fluid. The interaction between solid and liquid phases, which plays a major role in debris flow movement, is assumed to consist of drag and buoyancy forces. The applicability of drag force formulas is discussed. Considering the complex interaction between debris flow and the bed surface, a combined friction boundary condition is imposed on the bottom, and this is also discussed. To solve the complex model equations, a numerical method with second‐order accuracy based on the finite volume method is proposed. Several numerical experiments are performed to verify the feasibilities of model and numerical schemes. Numerical results demonstrate that different solid volume fractions substantially affect debris flow movement.  相似文献   

15.
Two‐phase fluid distributions in fractured porous media were studied using a single‐component multiphase (SCMP) lattice Boltzmann method (LBM), which was selected among three commonly used numerical approaches through a comparison against the available results of micro X‐ray computed tomography. The influence of the initial configuration and the periodic boundary conditions in the SCMP LBM for the fluid distribution analysis were investigated as well. It was revealed that regular porous media are sensitive to the initial distribution, whereas irregular porous media are insensitive. Moreover, to eliminate the influence of boundaries, the model's buffer size of an SCMP LBM simulation was suggested to be taken as approximately 12.5 times the average particle size. Then, the two‐phase fluid distribution of a porous medium was numerically studied using the SCMP LBM. Both detailed distribution patterns and macroscopic morphology parameters were reasonably well captured. Finally, the two‐phase fluid distributions in a fractured porous media were investigated. The influence of the degree of saturation, fracture length, and fracture width on the fluid distributions and migration was explored. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

16.
路基冻胀问题严重影响寒区高速铁路的安全服役,而成冰相变过程是解释冻胀机制的关键。基于介观尺度的格子Boltzmann方法,将修正的孔隙水冻结温度算法与焓法固液相变格子Boltzmann模型相结合,模拟了悬浮液滴冻结和冻土孔隙水成冰两个过程,分别揭示了液态水在自由状态和孔隙束缚状态下冰水相变的细观机制。计算结果表明：土体孔隙中冰晶由中心向外生长的过程与悬浮在空气中的液滴冻结过程截然不同,并且孔隙水越接近颗粒表面,其冻结温度越低。相同粒径颗粒按照不同排列方式得到的冻结特征曲线（soil freezing characteristic curves,简称SFCC）具有明显差异;不同粒径的SFCC随着颗粒增大残余水含量逐渐变少,形态更加陡峭。通过与文献试验结果对比,验证了格子Boltzmann方法的有效性,表明该方法能够为研究多孔介质水气迁移与相变过程提供介观尺度的新手段。  相似文献   

17.
The success in previous papers of this series in obtaining closed form solutions to the sedimentation of individual particles and to the suspensions of particles in fluids motivated the application of the same procedure to the turbulent flow of Newtonian fluids in conduits.The purpose of this paper is to present explicit equations for the Fanning friction factor and of the average velocity of a Newtonian fluid in a cylindrical tube, for any value of the Reynolds number and dimensionless roughness.Previously published equations for the friction coefficient are either implicit or explicit in the friction coefficient but dependent on the Reynolds number. Since the Reynolds number depends on the pipe diameter and fluid velocity, in both cases iterative method must be used to calculate the flowrate.By defining a dimensionless pipe diameter and a dimensionless average fluid velocity, and using conduit-flow material properties to define characteristic size and velocity parameters, explicit equations were obtained for the friction coefficient and for the average flow velocity valid for different roughness and flow regimes.  相似文献   

18.
The lattice Boltzmann method is a popular tool for pore-scale simulation of flow. This is likely due to the ease of including complex geometries such as porous media and representing multiphase and multifluid flows. Many advancements, including multiple relaxation times, increased isotropy, and others have improved the accuracy and physical fidelity of the method. Additionally, the lattice Bolzmann method is computationally very efficient, thanks to the explicit nature of the algorithm and relatively large amount of local work. The combination of many algorithmic options and efficiency means that a software framework enabling the usage and comparison of these advancements on computers from laptops to large clusters has much to offer. In this paper, we introduce Taxila LBM, an open-source software framework for lattice Boltzmann simulations. We discuss the design of the framework and lay out the features available, including both methods in the literature and a few new enhancements which generalize methods to complex geometries. We discuss the trade-off of accuracy and performance in various methods, noting how the Taxila LBM makes it easy to perform these comparisons for real problems. And finally, we demonstrate a few common applications in pore-scale simulation, including the characterization of permeability of a Berea sandstone and analysis of multifluid flow in heterogenous micromodels.  相似文献   

19.
This paper presents a two‐dimensional coupled bonded particle and lattice Boltzmann method (BPLBM) developed to simulate the fluid–solid interactions in geomechanics. In this new technique, the bonded particle model is employed to describe the inter‐particle movement and forces, and the bond between a pair of contacting particles is assumed to be broken when the tensile force or tangential force reaches a certain critical value. As a result the fracture process can be delineated based on the present model for the solid phase comprising particles, such as rocks and cohesive soils. In the meantime, the fluid phase is modelled by using the LBM, and the immersed moving boundary scheme is utilized to characterize the fluid–solid interactions. Based on the novel technique case studies have been conducted, which show that the coupled BPLBM enjoys substantially improved accuracy and enlarged range of applicability in characterizing the mechanics responses of the fluid–solid systems. Indeed such a new technique is promising for a wide range of application in soil erosion in Geotechnical Engineering, sand production phenomenon in Petroleum Engineering, fracture flow in Mining Engineering and fracture process in a variety of engineering disciplines. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

20.
This paper analyses the plane strain problem of a fracture, driven by injection of an incompressible viscous Newtonian fluid, which propagates parallel to the free surface of an elastic half‐plane. The problem is governed by a hyper‐singular integral equation, which relates crack opening to net pressure according to elasticity, and by the lubrication equations which describe the laminar fluid flow inside the fracture. The challenge in solving this problem results from the changing nature of the elasticity operator with growth of the fracture, and from the existence of a lag zone of a priori unknown length between the crack tip and the fluid front. Scaling of the governing equations indicates that the evolution problem depends in general on two numbers, one which can be interpreted as a dimensionless toughness and the other as a dimensionless confining stress. The numerical method adopted to solve this non‐linear evolution problem combines the displacement discontinuity method and a finite difference scheme on a fixed grid, together with a technique to track both crack and fluid fronts. It is shown that the solution evolves in time between two asymptotic similarity solutions. The small time asymptotic solution corresponding to the solution of a hydraulic fracture in an infinite medium under zero confining stress, and the large time to a solution where the aperture of the fracture is similar to the transverse deflection of a beam clamped at both ends and subjected to a uniformly distributed load. It is shown that the size of the lag decreases (to eventually vanish) with increasing toughness and compressive confining stress. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号