Numerical Modeling of Porous Flow in Fractured Rock and Its Applications in Geothermal Energy Extraction
|
|
- Osborn Payne
- 5 years ago
- Views:
Transcription
1 Journal of Earth Science, Vol. 6, No., p , February 05 ISSN X Printed in China DOI: 0.007/s Numerical Modeling of Porous Flow in Fractured Rock and Its Applications in Geothermal Energy Extraction Yucang Wang, Shimin Wang*, Sheng Xue, Deepak Adhikary. Earth Science and Resource Engineering, CSIRO, Kenmore Qld 4069, Brisbane, Australia. Key Laboratory of Computational Geodynamics, University of Chinese Academy of Sciences, Beijing 00049, China ABSTRACT: Understanding the characteristics of hydraulic fracture, porous flow and heat transfer in fractured rock is critical for geothermal power generation applications, and numerical simulation can provide a powerful approach for systematically and thoroughly investigating these problems. In this paper, we present a fully coupled solid-fluid code using discrete element method (DEM) and lattice Boltzmann method (LBM). The DEM with bonded particles is used to model the deformation and fracture in solid, while the LBM is used to model the fluid flow. The two methods are two-way coupled, i.e., the solid part provides a moving boundary condition and transfers momentum to fluid, while the fluid exerts a dragging force to the solid. Two widely used open source codes, the ESyS_Particle and the OpenLB, are integrated into one code and paralleled with Message Passing Interface (MPI) library. Some preliminary D simulations, including particles moving in a fluid and hydraulic fracturing induced by injection of fluid into a borehole, are carried out to validate the integrated code. The preliminary results indicate that the new code is capable of reproducing the basic features of hydraulic fracture and thus offers a promising tool for multiscale simulation of porous flow and heat transfer in fractured rock. KEY WORDS: discrete element method, lattice Boltzmann method, hydraulic fracturing, geothermal energy extraction, multiscale modelling. 0 INTRODUCTION Geothermal energy is clean, renewable, and steadily supplied by the Earth. The use of geothermal energy can help in both meeting the increasing energy needs over the world and protecting the environment of our planet. Among the various uses of geothermal energy, geothermal power generation is of particular importance. As the capacity factor of geothermal power plants can be up to 95% (much higher than wind power, solar power, and the direct usage of geothermal energy), geothermal power can serve as the base load of a county s power grid, and is also capable of cogeneration of electricity and heat (Pang et al., 0). In China, the potential of geothermal power is enormous, but the operating geothermal power capacity has not increased during the past two decades. Geothermal energy extraction process typically involves the injection of high pressure fluid to exchange heat from rocks at depths. In order to effectively and economically extract geothermal energy from the Earth, natural and artificial fractures are commonly used as the working fluid passage. Generally hydraulic fracturing is used to break the rock and increase permeability for the fluid flow. Although hydraulic *Corresponding author: smwang@ucas.ac.cn China University of Geosciences and Springer-Verlag Berlin Heidelberg 05 Manuscript received February 7, 04. Manuscript accepted May 9, 04. fracturing is a mature technique in the oil industry, there remain some issues which are not fully understood. The challenges in this process include how to stimulate and sustain the flow of fluid through the geothermal field and how to generate an efficient hydraulic subsurface heat exchanger system. In-depth research is needed to offer satisfactory prediction for types of fracturing as well as flow rates with respect to different operating conditions, such as the fluid viscosity, propellants used, pump rates, etc.. Therefore, understanding the characteristics of hydraulic fracture, porous flow and heat transfer in fractured rock is critical for geothermal power generation applications, including conventional and enhanced geothermal systems. Due to limitations of existing technologies and high costs, the mechanically-hydrodynamically-thermally coupled behaviors of fractured rocks cannot be studied solely relying on the in-situ geophysical measurements and laboratory tests, while numerical simulation provides an alternative approach for systematically and thoroughly investigating the characteristics of porous flow and heat transfer in fractured rocks. In this paper, we mainly focus on a microscopic model based on combining the discrete element method (DEM) and lattice Boltzmann method (LBM) to provide two-way coupled simulations for solid deformation and frature as well as thermal fluid flow. The DEM has emerged as a powerful numerical tool to model granular flow. Some DEM models allow particles to be bonded and the bonds to break, explicitly modeling microscopic fracturing events. Using DEM, many problems which are highly dynamic with large deformations and a large Wang, Y. C., Wang, S. M., Xue, S., et al., 05. Numerical Modeling of Porous Flow in Fractured Rock and Its Applications in Geothermal Energy Extraction. Journal of Earth Science, 6(): 0 7. doi:0.007/s
2 Numerical Modeling of Porous Flow in Fractured Rock and Its Applications in Geothermal Energy Extraction number of frequently changing contacts can be modeled easily. On the other hand, the LBM is a powerful tool to model fluid flow. It is based on the kinetic gas theory, which simulates fluid flows by tracking the evolution of fluid particle distribution. Some advantages of LBM over the classical Navier-Stokes approach include its ease in implementation, paralleling, and handling boundary conditions of complicated geometries. In this study, an integrated code coupling DEM and LBM is developed by combining two well developed open source codes: the ESyS_Particle and the OpenLB. The coupled approach is tested by some preliminary D simulations, including particles moving in a fluid and hydraulic fracturing induced by injection of fluid into a borehole. DISCRETE ELEMENT METHOD The discrete element method (Cundall and Strack, 979) is based on the concept that the modeled material can be represented as a collection of discrete solid particles interacting with one another at their contacts. The precise nature of the interaction depends on the scale of interest and the details of the simulation. At each time step, the calculations performed in DEM alternate between integrating equations of motion for each particle and applying the force-displacement law at each contact, through which the contact forces are updated based on the relative motions between two particles and their relevant contact stiffness. One kind of DEM allows particles to be bonded so that tensile forces can be transmitted. Fracturing is represented explicitly as broken bonds, which form and coalesce into macroscopic fractures. The bonded DEM model is often used to model wave propagation and rock fractures. The ESyS_Particle is an open source DEM software developed at the University of Queensland, Australia. Written in C++ and Python, the ESyS_Particle is designed for execution on parallel supercomputers based on spatial domain decomposition via the Message Passing Interface (MPI) (Abe et al., 004). The ESyS_Particle has been extended by the first author to include single particle rotation and a full set of interactions between particles (Wang, 009; Wang and Alonso-Marroquin, 009; Wang and Mora, 009; Wang et al., 006). The other major features of the ESyS_Particle include the explicit representation of particle orientations using unit quaternion, complete interactions (six kinds of independent relative movements are transmitted between two 3D interacting particles, see Fig. ) and a new method of decomposing the relative rotations between two rigid bodies such that the resulting torques and forces can be uniquely determined. Theoretical studies have been done to investigate the relationship between particle-scale stiffness and macroscopic elastic parameters (Wang and Mora, 008). The ESyS_Particle has been utilised in the study of rock fracture and earthquake dynamics (Mora and Place, 994, 993). The details of the ESyS_Particle code can be found in open literature (Wang et al., 0, 006; Wang, 009; Wang and Alonso-Marroquin, 009; Wang and Mora, 009, 008). Here two examples are given to show the capability of the code in modeling rock fracturing. Figure shows simulations of brittle crack extension in rock-like materials under impact by a hard ball. The brittle rock is modeled by bonded particles, and bonds between particles can break if a criterion is reached, explicitly modeling microscopic fracturing events. When the ball hits the rock, cracks induced by the impact are generated from the hit area and progressively propagated away. In Fig. b, the colors represent horizontal displacement of particles. A crushing area near the impact can be seen where the material is highly fragmented. Several groups of meridian cracks radiated from the impact area (color changes indicate displacement disconinuity) in vetical planes are observed. Figure 3 presents f s X Figure. Six kinds of interactions between bonded particles. f r is normal force, f s and f s are shear forces, τ t is twisting torque, and τ b and τ b are bending torques. (a) (b) τ b Figure. Simulation of brittle crack extension of rock-like materials under impact by a hard ball. (a) D; (b) 3D. Z f r τ t τ b f s Y
3 Yucang Wang, Shimin Wang, Sheng Xue and Deepak Adhikary (a) (b) (c) Figure 3. Progressive fracture development under bi-axial compression. Dog-eared shaped breakouts are formed. the results for a simulation of typical progressive borehole breakout development in the case of σ h /σ H =0., where σ H (in vertical direction) and σ h are, respectively, the maximum and minimum principle stresses. Two groups of conjugate cracks, dominated by shear fractures, start and intersect at both sides of the borehole surface, and then release the crushed grains (Fig. 3a). With the removal of some particles and continuous loading, stresses concentrate at the corner of the previous fault intersection, inducing several groups of shear fracture bands, and producing another major layer of fractures. Finally at both sides of the wall, a group of shear conjugate fractures develops progressively, forming a dog eared shape (Figs. 3b and 3c). Such dog eared fracture patterns have been widely observed in laboratory tests and in-situ surveys. LATTICE BOLTZMANN METHOD In contrast to the conventional computational fluid dynamics (CFD) techniques that solve macroscopic Navier- Stokes equations, LBM is built on a mesoscopic scale in which fluid is described by a group of discrete particles that propagate along a regular lattice and collide with each other. LBM solves the fluid particle distribution function f. The completely discretized equation, with the time step t and space step x, is given by BGK model (Chen and Doolen, 998) eq f e Δt, t Δt) () where τ denotes the lattice relaxation time, e α is the discrete lattice velocity in direction α, x i is a point in the discretized physical space, and f α eq is the equilibrium distribution function (see Eq. 4). Equation is usually solved in the following two steps (collision step and streaming step) ~ eq f, t Δt) () ~ f e Δt, t Δt), t Δt) (3) where f ~ represents the post-collision state. LBM uses the Eulerian approach to describe the motion of fluid particles; therefore the fixed Euler grids are adopted. A square and nine-velocity lattice, the DQ9 model (Fig. 4), is used in this study. The equilibrium distribution function is of the form (DQ9) f w ( ) c c c e u e u u u (4) eq 4 where c= x/ t is the lattice speed, ρ is the lattice fluid density, u is the macroscopic velocity, and w α is the weighting factor given by w 4 / 9, 0, /9,,3,5,7, /36,,4,6,8. The macroscopic quantities such as density, momentum, fluid pressure and viscosity can be obtained by 8 0 where c s c / 3 f, u 8 0 f s s (5) e, p c, v ( / ) c Δt is the speed of sound in this model. 3 TWO-WAY COUPLING OF DEM AND LBM In this study, two open source codes, the OpenLB and the ESyS_Particle, are used in coupling of LBM and DEM. At each time step, the fluid and solid components run in a separate and independent fashion, while the two codes are integrated into one combined system to facilitate a fully coupled solution. In order to implement such coupling, the following issues need to be considered: moving boundary conditions for a curved shape for LBM which handles the momentum transfer from solid particle to fluid and the reflection of fluid at the solid boundaries, as well as force and momentum transfer from fluid to solid particles. The most widely used boundary condition in LBM is no-slip bounce-back condition in which an incoming fluid particle from a fluid node is reflected back to the node it comes from, and velocity or pressure boundary condition (Zou and He, 997) in which velocity or density along the boundary node is fixed. These boundary conditions are only applicable in the case of non-moving objects and straight boundaries. For moving objects inside the fluid it is necessary to adapt the moving boundary conditions. Curved boundary treatments provide a means of improving the computational accuracy of the stair-shaped approximation used in LBM simulations. Due to the arbitrary position and curved surfaces of the solid particles, (6)
4 Numerical Modeling of Porous Flow in Fractured Rock and Its Applications in Geothermal Energy Extraction 3 the particle surface can intersect the link between two nodes of LBM at an arbitrary distance with a ratio of x f xw 0, (7) x x f b where x f and x b denote the lattice node on the fluid side of the boundary and that on the solid side. When δ=0, x w =x f, the boundary wall locates at the point x f ; when δ=, x w =x b, the boundary wall locates at the point x b. The reflected distribution function at node x f can be calculated using an interpolation scheme (Yu et al., 003) f ( xf, t Δ t) ( ) ( xf, t Δ t) ( xb, tδ t) f ( xf, t Δ t) 6 w we uw / c where w α is the same weight factor as in Eqs. 4 and 5, ρ w is the fluid density at node x f, u w is the velocity of the solid particle, and x f =x f e α t (Fig. 5). The last term in Eq. 8 represents the momentum transferred from solid particle to fluid. The fluid force acting on the solid particle surface can be obtained using the formula 9 F ( b, ) ( f, Δ ) (Δ ) /Δ xb (8) f e f x t f x t t x t (9) where the first summation is taken over all solid boundary nodes at x b adjacent to the fluid nodes and the second summation is taken over all possible lattice directions pointing from x b towards all possible neighboring fluid nodes around the solid node x b. This force is added to the solid particle force in DEM code. Hence, a combination of Eq. 8 and Eq. 9 enables fully two-way coupled solid-fluid simulations. To implementing this boundary condition, a single interpolation is used for both δ<0.5 and δ 0.5, which is found to be very stable. 4 PRELIMINARY RESULTS To assess the performance of the coupled DEM-LBM approach, two simple preliminary numerical results are presented here. At this stage, validations of large scale simulations and detailed comparisons with physical experiments have not yet 6 3 e 3 e 7 e 6 e Figure 4. A widely used -D 9-velocity lattice model (DQ9). The cell in the center is stationary, representing fluid particles at rest. The fluid particles are only allowed to move along the eight moving directions to neighboring cells or to rest in this cell. e e 5 e 8 e 5 x f e ~ x f Boundary wall e ~ e x w x b Figure 5. The moving curved wall boundary condition. The curved line represents the boundary wall between a solid particle and the fluid. x f and x f are the lattice nodes on the fluid side of the boundary; and x b is the node on the solid side; x w denotes the intersection of the wall with lattice link.
5 4 Yucang Wang, Shimin Wang, Sheng Xue and Deepak Adhikary been carried out since the coupled code is just integrated. These small scaled models, however, can be used as a qualitative display of the capability and potential of the coupled approach. Figure 6 shows the modeling results of solid particle motion in a fluid. Three particles bonded as a single rigid body with unbroken bonds. These particles are driven by the fluid flow in a long pipe. A boundary condition of a given parabolic velocity profile is specified for the inlet (the left boundary), while a constant pressure condition is applied at the outlet (the right boundary) to facilitate a Poiseuille-type flow. The motion and rotation of the rigid body, as well as vortex behind it, are clearly observed. In this simulation, Reynolds number Re is 00 and the solid particle size is 5 times greater than the grid size of the fluid. This example verifies the algorithm of the two-way coupling of solid and fluid, including moving boundary condition and dragging forces. Figures 7 to 9 show a D simulation of a hydraulic fracturing process. In this test, rock is modeled by 7 37 DEM particles with a hole in the middle. Particles with sizes ranging from 0. to unit are bonded by brittle elastic bonds. When the stress in a bond is large enough, the bond breaks, explicitly modeling a micro-fracturing event. The release of stress on the broken bond will transfer stress to the neighboring bonds, possibly causing them to break in a later time. In this way, a progressive failure process is modeled by extensions of micro-cracks. The fluid flow is simulated using 0 0 LBM grids. The size of the fluid grid is half of the minimum solid particle radius. Four motionless walls provide constant displacement boundary conditions, and the friction between walls and particles is set to zero to minimize influence of the boundary. Constant pressure boundary conditions are set for the LBM fluid. Fluid pressure in the center of the hole is increased slowly, modeling the injection of water. A pressure increase is realized by adding a source term on the right-hand side of Eq., resulting in f e Δt, t Δt ) f, t ) f, t ) f eq, t ) Δf, t ) (0) where Δfα(xi,t) in each direction is proportional to fα(xi,t) and Figure 6. Three particles bonded as a rigid body and move in the fluid. Particles are driven by the drag forces of fluid flow. Colors represent fluid scalar velocity, with blue for low and red for high. Figure 7. Snapshots of hydraulic fracture. Fluid pressure in the centre of the hole is increased slowly, modeling the injection of water. Colors represent solid particle displacement (blue for low and red for high) in horizontal direction. Fracture initiation and propagation due to hydraulic loading is clearly seen.
6 Numerical Modeling of Porous Flow in Fractured Rock and Its Applications in Geothermal Energy Extraction 5 Figure 8. Fracturing event distribution. Tensile fractures are dominant in the process of the crack initiation and extension. Figure 9. Fluid pressure (blue for low and red for high) during hydraulic fracturing. Only those LBM grids covered by fluid are visualized, representing the process of cracking and fluid flowing in the hydraulically generated fracturing tunnels. the increase of fluid pressure in one time step is determined by Δ p cs 8 0 Δ f. Figure 7 shows several snapshots of fracture initiation and propagation due to hydraulic loading. Different colors in this figure represent the different magnitudes of solid particle displacement in horizontal direction. The cracks are observed to start from the surface of the borehole, caused by the increasing fluid pressure on the surface of the borehole, and propagate inside the rock with the continuous breakage of bonds and infiltration of pressured fluid. Figure 8 plots evolution of the micro-fracturing of bonds. The sizes of the event represent fracturing energy released. It is clearly seen that tensile fractures are dominant in the beginning of the crack initiation. The fluid pressure is visualized in Fig. 9. Since there are many LBM grids ( 0 0), these plots can rather be interpreted as the presence and extension of the fluid. In the initial state, fluid only occupies the borehole area. With the occurrence of fractures, fluid is driven by pressure into the cracks, and then exerts further pressure on the crack surface, which results in more fractures generated. Although the current simulation is still rudimental and of small scale, and the influence of input parameters, such as particle-scale stiffness, fracturing parameters, confining pressure and injection rate, on the fracturing behavior is not investigated, the simulation includes the major mechanisms of hydraulic fracturing process: (i) mechanical deformation and fracturing induced by the fluid pressure; (ii) flow of fluid within the fracture; and (iii) fracture propagation. Furthermore, there are no specific assumptions about how fluid flows inside cracks and how a solid breaks, i.e., they are treated in a straightforward way in the new coupled model. 5 DISCUSSION Although the simulation results presented above remain un-calibrated against laboratory and in-situ testing, they do reproduce the basic features of hydraulic fracturing and demonstrate the promising applications of the new modeling approach. The coupled code requires significant improvements to realize its full potentials. Firstly, the current coupled code is still D and needs to be extended to 3D, which should be straightforward as both the ESyS_Particle and the OpenLB are 3D codes and the coupled code is structured to have provisions for 3D coupling. Secondly, thermal effect is required to achieve a full hydro- and mechanical-coupling capability. Currently, the ESyS_Particle has included thermal effects, including heat transfer, thermal expansion and friction generated heat, and the OpenLB also has this capability. Once these improvements have been incorporated, the fully developed model is capable of efficiently simulating the dynamics and evolution of real deep geothermal reservoir systems. The existing numerical models for extracting geothermal energy are almost all based on the continuum approach (Hunsweck et al., 03; Secchi and Schrefler, 0; Zhang et al., 009). Although these models have been shown to be valuable in geothermal energy exploration, as reported in a number of case studies in the literature (Zhou and Hou, 03; Gong et al., 0; Blocher et al., 00; Bataille et al., 006), they are incapable of capturing the full physics of the modeled fractured rocks. As fractures in rock normally involve various geometric
7 6 Yucang Wang, Shimin Wang, Sheng Xue and Deepak Adhikary scales, the continuum models may be adequate to represent the large scale fractures, but definitely not for complicated micro-fracture systems. A physically sound and computational affordable numerical approach for modeling the fractured rocks is in need but not available yet. We belive multiscale numerical modeling is an ideal approach to investigate the mechanism and characteristics of porous flow and heat transfer in fractured rock as applied to the real world problems of geothermal power generation. In the mesoscale simulation, the DEM can be adopted to model the deformation in solid as well as the formation, enlargement, connection and closure behaviors of microfractures, while the LBM can be adopted to model the fluid flow in pores and fractures; the combination of these two methods provides a fully coupled simulation of solid deformation, fluid flow as well as heat conduction and convection. Based on the analysis and integration of the mesoscale modeling results, quantitative relations of parameters like permeability and thermal dispersion in fractured rock as functions of pressure, deviatoric stress and temperature can be determined, and then the macroscale constitutive relations of the porous flow and heat transfer in fractured rock can be derived, which form the basis for constructing continuum models to simulate the macroscale hydro-thermal behaviors of fractured rocks. Such multiscale modeling appoach may be employed to simulate the porous flow and heat transfer processes between injection and production wells in geothermal power plants, to study the effects of injection pressure, injection rate and natural convection on the geothermal reservoir temperature and its serving life, and to investigate the optimal distribution of the injection and production wells as well as the technological strategy for wisely and economically extracting deep geothermal resources. 6 CONCLUSIONS This paper presents a two-way coupled DEM-LBM scheme based on two open source codes, the ESyS_Particle for DEM and the OpenLB for LBM. The ESyS_Particle solves solid particle motion, interaction between particles, and fracturing of bonds, while the OpenLB is responsible for fluid motion on the fixed Euler grids. The two-way coupling is achieved through drag forces applied to DEM particles by fluid, moving boundaries, and momentum transfer from solid particles to fluid grids. The preliminary results of D simulations for particles moving in a fluid and hydraulic fracturing induced by injection of fluid into a borehole indicate that the new code is capable of solving coupled fluid-solid interaction problems and reproducing the basic features of hydraulic fracturing. As a result, the numerical approach proposed in this study offers a promising tool for multiscale simulation of porous flow and heat transfer in fractured rock. ACKNOWLEDGMENTS This study was financially supported by the CSIRO Earth Science and Resource Engineering Geothermal Energy Capability Development Fund, the Huainan Coal Mining Group in China, and the National Natural Science Foundation of China (No ). It was also supported by the NCI National Facility at the ANU and ivec at the UWA through the use of advanced computing resources. Two open source codes, the ESyS_Particle and the OpenLB, are used in this study. We thank the developers of these codes. REFERENCES CITED Abe, S., Place, D., Mora, P., 004. A Parallel Implementation of the Lattice Solid Model for the Simulation of Rock Mechanics and Earthquake Dynamics. Pure Appl. Geophys., 6: Bataille, A., Genthon, P., Rabinowicz, M., et al., 006. Modeling the Coupling between Free and Forced Convection in a Vertical Permeable Slot: Implications for the Heat Production of an Enhanced Geothermal System. Geothermics, 35: Blocher, M. G., Zimmermann, G., Moeck, I., et al., 00. 3D Numerical Modeling of Hydrothermal Processes during the Lifetime of a Deep Geothermal Reservoir. Geofluids, 0: Chen, S., Doolen, G., 998. Lattice Boltzmann Method for Fluid Flows. Ann. Rev. Fluid Mech., 30: Cundall, P. A., Strack, O. D. L., 979. Discrete Numerical Model for Granular Assemblies. Geotechnique, 9: Gong, B., Liang, H., Xin, S., et al., 0. Effect of Water Injection on Reservoir Temperature during Power Generation in Oil Fields. Proceedings of 36th Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford. 3 January February, 0 Hunsweck, M. J., Shen, Y., Lew, A. J., 03. A Finite Element Approach to the Simulation of Hydraulic Fractures with Lag. Int. J. Numer. Anal. Meth. Geomech., 37: Mora, P., Place, D., 993. A Lattice Solid Model for the Nonlinear Dynamics of Earthquakes. Int. J. Mod. Phys., C4: Mora, P., Place, D., 994. Simulation of the Frictional Stick-Slip Instability. Pure Appl. Geophys., 43: 6 87 Pang, Z. H., Hu, S. B., Wang, J. Y., 0. A Roadmap to Geothermal Energy Development in China. Science & Technology Review, 30(3): 8 4 (in Chinese with English Abstract) Secchi, S., Schrefler, B. A., 0. A Method for 3-D Hydraulic Fracturing Simulation. Int. J. Fract., 78: Wang, Y. C., 009. A New Algorithm to Model the Dynamics of 3-D Bonded Rigid Bodies with Rotations. Acta Geotechnica, 4: 7 7 Wang, Y. C., Abe, S., Latham, S., et al., 006. Implementation of Particle-Scale Rotation in the 3D Lattice Solid Model. Pure Appl. Geophys., 63: Wang, Y. C., Alonso-Marroquin, F., 009. A New Discrete Element Model: Particle Rotation and Parameter Calibration. Granular Matter, : Wang, Y. C., Mora, P., 008. Elastic Properties of Regular Lattices. J. Mech. Phys. Solids, 56: Wang, Y. C., Mora, P., 009. ESyS_Particle: A New 3-D Discrete Element Model with Single Particle Rotation. In: Xing, H. L., ed., Advances in Geocomputing. Springer. 83 8
8 Numerical Modeling of Porous Flow in Fractured Rock and Its Applications in Geothermal Energy Extraction 7 Wang, Y. C., Xue, S., Xie, J., 0. Discrete Element Method and Its Applications in Earthquake and Rock Fracture Modeling. In: Li, Y. G., ed., Imaging, Modeling and Assimilation in Seismology. China High Education Press, Beijing; De Gruyter, Boston Yu, D., Mei, R., Luo, L., et al., 003. Viscous Flow Computations with the Method of Lattice Boltzmann Equation. Prog. Aerospace Sci., 39: Zhang, X., Jeffrey, R. G., Thiercelin, M., 009. Mechanics of Fluid-Driven Fracture Growth in Naturally Fractured Reservoirs with Simple Network Geometries. J. Geophys. Res., 4: B406. doi:0.09/009jb Zhou, L., Hou, M. Z., 03. A New Numerical 3D-Model for Simulation of Hydraulic Fracturing in Consideration of Hydro-Mechanical Coupling Effects. International Journal of Rock Mechanics & Mining Sciences, 60: Zou, Q., He, X., 997. On Pressure and Velocity Boundary Conditions for the Lattice Boltzmann BGK Model. Phys. Fluids, 9:
Research of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method
Research Journal of Applied Sciences, Engineering and Technology 6(14): 50-55, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: November 08, 01 Accepted: December 8,
More informationLattice Boltzmann Method for Fluid Simulations
1 / 16 Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas Simon Fraser University April 7, 2011 2 / 16 Ludwig Boltzmann and His Kinetic Theory of Gases The Boltzmann Transport
More informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationDiscrete Element Modeling of Thermo-Hydro-Mechanical Coupling in Enhanced Geothermal Reservoirs
PROCEEDINGS, Thirty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California Discrete Element Modeling of Thermo-Hydro-Mechanical Coupling in Enhanced Geothermal Reservoirs
More informationOutline. Advances in STAR-CCM+ DEM models for simulating deformation, breakage, and flow of solids
Advances in STAR-CCM+ DEM models for simulating deformation, breakage, and flow of solids Oleh Baran Outline Overview of DEM in STAR-CCM+ Recent DEM capabilities Parallel Bonds in STAR-CCM+ Constant Rate
More informationExternal and Internal Incompressible Viscous Flows Computation using Taylor Series Expansion and Least Square based Lattice Boltzmann Method
Available online at http://ijim.srbiau.ac.ir/ Int. J. Industrial Mathematics (ISSN 2008-5621) Vol. 10, No. 2, 2018 Article ID IJIM-00726, 8 pages Research Article External and Internal Incompressible Viscous
More informationModeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures
Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures ENRU LIU, SERAFEIM VLASTOS, and XIANG-YANG LI, Edinburgh
More informationModeling of Propagation of Interacting Cracks Under Hydraulic Pressure Gradient. Hai Huang, Earl Mattson and Rob Podgorney
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Modeling of Propagation of Interacting Cracks Under Hydraulic Pressure Gradient Hai Huang, Earl Mattson and Rob Podgorney
More informationSimulation of floating bodies with lattice Boltzmann
Simulation of floating bodies with lattice Boltzmann by Simon Bogner, 17.11.2011, Lehrstuhl für Systemsimulation, Friedrich-Alexander Universität Erlangen 1 Simulation of floating bodies with lattice Boltzmann
More informationTable of Contents. Preface... xiii
Preface... xiii PART I. ELEMENTS IN FLUID MECHANICS... 1 Chapter 1. Local Equations of Fluid Mechanics... 3 1.1. Forces, stress tensor, and pressure... 4 1.2. Navier Stokes equations in Cartesian coordinates...
More informationThe Influence of Contact Friction on the Breakage Behavior of Brittle Granular Materials using DEM
The Influence of Contact Friction on the Breakage Behavior of Brittle Granular Materials using DEM *Yi-Ming Liu 1) and Hua-Bei Liu 2) 1), 2) School of Civil Engineering and Mechanics, Huazhong University
More informationSimulation of 2D non-isothermal flows in slits using lattice Boltzmann method
Simulation of 2D non-isothermal flows in slits using lattice Boltzmann method João Chambel Leitão Department of Mechanical Engineering, Instituto Superior Técnico, Lisboa, Portugal Abstract: The present
More informationBrittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm
Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd
More informationLATTICE BOLTZMANN SIMULATION OF FLUID FLOW IN A LID DRIVEN CAVITY
LATTICE BOLTZMANN SIMULATION OF FLUID FLOW IN A LID DRIVEN CAVITY M. Y. Gokhale, Ignatius Fernandes Maharashtra Institute of Technology, Pune 4 38, India University of Pune, India Email : mukundyg@yahoo.co.in,
More informationDiscrete Element Modelling of a Reinforced Concrete Structure
Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr
More informationLattice Bhatnagar Gross Krook model for the Lorenz attractor
Physica D 154 (2001) 43 50 Lattice Bhatnagar Gross Krook model for the Lorenz attractor Guangwu Yan a,b,,liyuan a a LSEC, Institute of Computational Mathematics, Academy of Mathematics and System Sciences,
More informationPREDICTION OF INTRINSIC PERMEABILITIES WITH LATTICE BOLTZMANN METHOD
PREDICTION OF INTRINSIC PERMEABILITIES WITH LATTICE BOLTZMANN METHOD Luís Orlando Emerich dos Santos emerich@lmpt.ufsc.br Carlos Enrique Pico Ortiz capico@lmpt.ufsc.br Henrique Cesar de Gaspari henrique@lmpt.ufsc.br
More informationA Momentum Exchange-based Immersed Boundary-Lattice. Boltzmann Method for Fluid Structure Interaction
APCOM & ISCM -4 th December, 03, Singapore A Momentum Exchange-based Immersed Boundary-Lattice Boltzmann Method for Fluid Structure Interaction Jianfei Yang,,3, Zhengdao Wang,,3, and *Yuehong Qian,,3,4
More informationSimulation of the cutting action of a single PDC cutter using DEM
Petroleum and Mineral Resources 143 Simulation of the cutting action of a single PDC cutter using DEM B. Joodi, M. Sarmadivaleh, V. Rasouli & A. Nabipour Department of Petroleum Engineering, Curtin University,
More informationAPPRAISAL OF FLOW SIMULATION BY THE LATTICE BOLTZMANN METHOD
APPRAISAL OF FLOW SIMULATION BY THE LATTICE BOLTZMANN METHOD Guillermo Izquierdo Bouldstridge Imperial College of London Department of Aeronautics Master s Thesis Supervisor: Dr. Joaquim Peiró September
More informationDISCRETE ELEMENT SIMULATIONS OF WATER FLOW THROUGH GRANULAR SOILS
15th ASCE Engineering Mechanics Conference June 2-5, 2002, Columbia University, New York, NY EM 2002 DISCRETE ELEMENT SIMULATIONS OF WATER FLOW THROUGH GRANULAR SOILS Usama El Shamy 1, Student Member,
More informationThe Mine Geostress Testing Methods and Design
Open Journal of Geology, 2014, 4, 622-626 Published Online December 2014 in SciRes. http://www.scirp.org/journal/ojg http://dx.doi.org/10.4236/ojg.2014.412046 The Mine Geostress Testing Methods and Design
More informationPredictor-Corrector Finite-Difference Lattice Boltzmann Schemes
Applied Mathematical Sciences, Vol. 9, 2015, no. 84, 4191-4199 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.5138 Predictor-Corrector Finite-Difference Lattice Boltzmann Schemes G. V.
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
More information3D HM-DEM model for Hydro-Fracturing
3D HM-DEM model for Hydro-Fracturing E. Papachristos, F.V. Donzé & B. Chareyre Laboratoire Sols, Solides, Structures, Grenoble, France efthymios.papachristos@3sr-grenoble.fr,, frederic.donze@3srgrenoble.fr,
More informationSimulation of Rarefied Gas Flow in Slip and Transitional Regimes by the Lattice Boltzmann Method
www.cfdl.issres.net Vol. 2 (2) June 2010 Simulation of Rarefied Gas Flow in Slip and Transitional Regimes by the Lattice Boltzmann Method N. Azwadi C. Sidik 1C, N.C. Horng 1, M.A. Mussa 2, and S. Abdullah
More informationThe Use of Lattice Boltzmann Numerical Scheme for Contaminant Removal from a Heated Cavity in Horizontal Channel
www.cfdl.issres.net Vol. 6 (3) September 2014 The Use of Lattice Boltzmann Numerical Scheme for Contaminant Removal from a Heated Cavity in Horizontal Channel Nor Azwadi Che Sidik C and Leila Jahanshaloo
More informationRole of lithological layering on spatial variation of natural and induced fractures in hydraulic fracture stimulation
Role of lithological layering on spatial variation of natural and induced fractures in hydraulic fracture stimulation Vincent Roche *, Department of Physics, University of Alberta, Edmonton roche@ualberta.ca
More informationCritical Borehole Orientations Rock Mechanics Aspects
Critical Borehole Orientations Rock Mechanics Aspects By R. BRAUN* Abstract This article discusses rock mechanics aspects of the relationship between borehole stability and borehole orientation. Two kinds
More informationEffect Of The In-Situ Stress Field On Casing Failure *
Effect Of The In-Situ Stress Field On Casing Failure * Tang Bo Southwest Petroleum Institute, People's Republic of China Lian Zhanghua Southwest Petroleum Institute, People's Republic of China Abstract
More informationCopyright 2014 The Authors. Deposited on: 25 June 2014
Sachau, T., and Koehn, D. (2014) A new mixed-mode fracture criterion for large-scale lattice models. Geoscientific Model Development, 7 (1). pp. 243-247. ISSN 1991-959X Copyright 2014 The Authors http://eprints.gla.ac.uk/94706
More informationEFFECT OF A CLUSTER ON GAS-SOLID DRAG FROM LATTICE BOLTZMANN SIMULATIONS
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 EFFECT OF A CLUSTER ON GAS-SOLID DRAG FROM LATTICE BOLTZMANN SIMULATIONS Milinkumar
More informationTHE EFFECT OF THERMOELASTIC STRESS CHANGE IN THE NEAR WELLBORE REGION ON HYDRAULIC FRACTURE GROWTH
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, 30 Jan 2011-1 Feb 2012 THE EFFECT OF THERMOELASTIC STRESS CHANGE IN THE NEAR WELLBORE
More informationOn pressure and velocity boundary conditions for the lattice Boltzmann BGK model
On pressure and velocity boundary conditions for the lattice Boltzmann BGK model Qisu Zou Theoretical Division, Los Alamos National Lab, Los Alamos, New Mexico 7545 and Department of Mathematics, Kansas
More informationLattice Boltzmann Modeling of Wave Propagation and Reflection in the Presence of Walls and Blocks
Lattice Boltzmann Modeling of Wave Propagation and Reflection in the Presence of Walls and Blocs Maysam Saidi, Hassan Basirat Tabrizi, Member, IAENG, and Reza Sepahi Samian Abstract Lattice Boltzmann method
More informationSTUDY ON FAILURE PROCESS OF TAILING DAMS BASED ON PARTICLE FLOW THEORIES
ISSN 176-459 Int j simul model 14 (015) 4, 658-668 Original scientific paper STUDY ON FAILURE PROCESS OF TAILING DAMS BASED ON PARTICLE FLOW THEORIES Yuan, L.-W. * ; Li, S.-M. *# ; Peng, B. ** & Chen,
More information68 Guo Wei-Bin et al Vol. 12 presented, and are thoroughly compared with other numerical data with respect to the Strouhal number, lift and drag coeff
Vol 12 No 1, January 2003 cfl 2003 Chin. Phys. Soc. 1009-1963/2003/12(01)/0067-08 Chinese Physics and IOP Publishing Ltd Lattice-BGK simulation of a two-dimensional channel flow around a square cylinder
More informationSimulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang,
Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, yfwang09@stanford.edu 1. Problem setting In this project, we present a benchmark simulation for segmented flows, which contain
More informationLattice Boltzmann Method for Moving Boundaries
Lattice Boltzmann Method for Moving Boundaries Hans Groot March 18, 2009 Outline 1 Introduction 2 Moving Boundary Conditions 3 Cylinder in Transient Couette Flow 4 Collision-Advection Process for Moving
More informationLattice Boltzmann Method
3 Lattice Boltzmann Method 3.1 Introduction The lattice Boltzmann method is a discrete computational method based upon the lattice gas automata - a simplified, fictitious molecular model. It consists of
More informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
More informationMain Means of Rock Stress Measurement
Real Stress Distributions through Sedimentary Strata and Implications for Reservoir Development and Potential Gas and Coal Development Strategies Ian Gray Sigra Pty Ltd 93 Colebard St West, Acacia Ridge,
More informationarxiv: v1 [physics.flu-dyn] 10 Aug 2015
Slip velocity of lattice Boltzmann simulation using bounce-back boundary scheme Jianping Meng, Xiao-Jun Gu, and David R Emerson Scientific Computing Department, STFC Daresbury laboratory, arxiv:1508.009v1
More informationAssessment of the Accuracy of the Multiple-Relaxation-Time Lattice Boltzmann Method for the Simulation of Circulating Flows
Mathematical Modelling and Applications 2017; 25): 47-51 http://www.sciencepublishinggroup.com/j/mma doi: 10.11648/j.mma.20170205.11 ISSN: 2575-1786 Print); ISSN: 2575-1794 Online) Assessment of the Accuracy
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
More informationLattice Boltzmann model for the Elder problem
1549 Lattice Boltzmann model for the Elder problem D.T. Thorne a and M.C. Sukop a a Department of Earth Sciences, Florida International University, PC 344, University Park, 11200 SW 8th Street, Miami,
More information2. Governing Equations. 1. Introduction
Multiphysics Between Deep Geothermal Water Cycle, Surface Heat Exchanger Cycle and Geothermal Power Plant Cycle Li Wah Wong *,1, Guido Blöcher 1, Oliver Kastner 1, Günter Zimmermann 1 1 International Centre
More informationAnalysis of Fracture Propagation under Thermal Stress in Geothermal Reservoirs
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Analysis of Fracture Propagation under Thermal Stress in Geothermal Reservoirs Ahmad Ghassemi, Sergej Tarasovs Mailing
More informationFinite element modelling of fault stress triggering due to hydraulic fracturing
Finite element modelling of fault stress triggering due to hydraulic fracturing Arsalan, Sattari and David, Eaton University of Calgary, Geoscience Department Summary In this study we aim to model fault
More informationAnalysis and boundary condition of the lattice Boltzmann BGK model with two velocity components
Analysis and boundary condition of the lattice Boltzmann BGK model with two velocity components Xiaoyi He and Qisu Zou arxiv:comp-gas/9507002v1 22 Jul 1995 Abstract In this paper, we study the two dimensional
More informationA TIME-DEPENDENT DAMAGE LAW IN DEFORMABLE SOLID: A HOMOGENIZATION APPROACH
9th HSTAM International Congress on Mechanics Limassol, Cyprus, - July, A TIME-DEPENDENT DAMAGE LAW IN DEFORMABLE SOLID: A HOMOGENIZATION APPROACH Cristian Dascalu, Bertrand François, Laboratoire Sols
More informationStudy on lattice Boltzmann method/ large eddy simulation and its application at high Reynolds number flow
Research Article Study on lattice Boltzmann method/ large eddy simulation and its application at high Reynolds number flow Advances in Mechanical Engineering 1 8 Ó The Author(s) 2015 DOI: 10.1177/1687814015573829
More informationHybrid Finite-Discrete Element Modelling of Dynamic Fracture of Rock and Resultant Fragment Arching by Rock Blast
Hybrid Finite-Discrete Element Modelling of Dynamic Fracture of Rock and Resultant Fragment Arching by Rock Blast H.M. An a, b and H.Y. Liu a * a School of Engineering and ICT, the University of Tasmania,
More informationIntegrating Lab and Numerical Experiments to Investigate Fractured Rock
Integrating Lab and Numerical Experiments to Investigate Fractured Rock Bradford H. Hager Director, Earth Resources Laboratory and Cecil and Ida Green Professor Department of Earth, Atmospheric and Planetary
More informationROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT
ROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT PETER C. LICHTNER 1, QINJUN KANG 1 1 Los Alamos National Laboratory,
More informationNUMERICAL ANALYSIS OF THE THREE-MATERIAL DOWNHOLE FLOW FIELD IN HYDROTHERMAL JET DRILLING
2017 WJTA-IMCA Conference and Expo October 25-27, 2017 New Orleans, Louisiana Paper NUMERICAL ANALYSIS OF THE THREE-MATERIAL DOWNHOLE FLOW FIELD IN HYDROTHERMAL JET DRILLING Xianzhi Song, Zehao Lyu, Haizhu
More informationElectromagnetic excitation of particle suspensions in hydraulic fractures using a coupled lattice Boltzmanndiscrete
Electromagnetic excitation of particle suspensions in hydraulic fractures using a coupled lattice Boltzmanndiscrete element model The MIT Faculty has made this article openly available. Please share how
More informationEXTENDED FREE SURFACE FLOW MODEL BASED ON THE LATTICE BOLTZMANN APPROACH
METALLURGY AND FOUNDRY ENGINEERING Vol. 36, 2010, No. 2 Micha³ Szucki*, Józef S. Suchy***, Pawe³ ak*, Janusz Lelito**, Beata Gracz* EXTENDED FREE SURFACE FLOW MODEL BASED ON THE LATTICE BOLTZMANN APPROACH
More informationAdvantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations
The Hilltop Review Volume 7 Issue 1 Winter 2014 Article 10 December 2014 Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations Tai-Hsien Wu Western Michigan University
More informationA micromechanical approach to describe internal erosion effects in soils
A micromechanical approach to describe internal erosion effects in soils Luc Scholtès, Pierre-Yves Hicher, Luc Sibille To cite this version: Luc Scholtès, Pierre-Yves Hicher, Luc Sibille. A micromechanical
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Surface and body forces Tensors, Mohr circles. Theoretical strength of materials Defects Stress concentrations Griffith failure
More informationParticle flow simulation of sand under biaxial test
5th International Conference on Civil Engineering and Transportation (ICCET 2015) Particle flow simulation of sand under biaxial test Xiao-li Dong1,2, a *,Wei-hua Zhang1,a 1 Beijing City University, China
More informationAn Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation
An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation Nachiket Patil, Deepankar Pal and Brent E. Stucker Industrial Engineering, University
More informationNUMERICAL SIMULATION OF THREE DIMENSIONAL GAS-PARTICLE FLOW IN A SPIRAL CYCLONE
Applied Mathematics and Mechanics (English Edition), 2006, 27(2):247 253 c Editorial Committee of Appl. Math. Mech., ISSN 0253-4827 NUMERICAL SIMULATION OF THREE DIMENSIONAL GAS-PARTICLE FLOW IN A SPIRAL
More informationEnabling Technologies
Enabling Technologies Mechanical Modelling 1 Key Parameter Mineral System Exploration is reflected in scale-dependent translation A. Gradient in hydraulic potential B. Permeability C. Solubility sensitivity
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationComputational Simulation of the Hydraulic Fracturing Process
Computational Simulation of the Hydraulic Fracturing Process Sergio Andrés Galindo Torrres Jose Daniel Muñoz Castaño Simulation of Physical Systems Group Universidad Nacional de Colombia sagalindot@unal.edu.co
More informationModeling of Humidification in Comsol Multiphysics 4.4
Modeling of Humidification in Comsol Multiphysics 4.4 Indrajit Wadgaonkar *1 and Suresh Arikapudi 1 1 Tata Motors Ltd. Pimpri, Pune, India, 411018. *Corresponding author: Indrajit Wadgaonkar, Tata Motors
More informationThe Lift Force on a Spherical Particle in Rectangular Pipe Flow. Houhui Yi
CHINESE JOURNAL OF PHYSICS VOL. 52, NO. 1-I February 2014 The Lift Force on a Spherical Particle in Rectangular Pipe Flow Houhui Yi Institute of Theoretical Physics, Department of Opto-Electronic Engineering,
More informationApplications of parabolized stability equation for predicting transition position in boundary layers
Appl. Math. Mech. -Engl. Ed., 33(6), 679 686 (2012) DOI 10.1007/s10483-012-1579-7 c Shanghai University and Springer-Verlag Berlin Heidelberg 2012 Applied Mathematics and Mechanics (English Edition) Applications
More informationNatural frequency analysis of fluid-conveying pipes in the ADINA system
Journal of Physics: Conference Series OPEN ACCESS Natural frequency analysis of fluid-conveying pipes in the ADINA system To cite this article: L Wang et al 2013 J. Phys.: Conf. Ser. 448 012014 View the
More informationNUMERICAL MODELING OF THE GAS-PARTICLE FLUID FLOW AND HEAT TRANSFER IN THE SLIP REGIME
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00394 NUMERICAL MODELING OF THE GAS-PARTICLE FLUID FLOW AND HEAT TRANSFER IN THE SLIP
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationSimulation of Lid-driven Cavity Flow by Parallel Implementation of Lattice Boltzmann Method on GPUs
Simulation of Lid-driven Cavity Flow by Parallel Implementation of Lattice Boltzmann Method on GPUs S. Berat Çelik 1, Cüneyt Sert 2, Barbaros ÇETN 3 1,2 METU, Mechanical Engineering, Ankara, TURKEY 3 METU-NCC,
More informationCalculating Method for the Axial Force of Washover String During Extracting Casing in Directional Well
Advances in Petroleum Exploration and Development Vol. 9, No., 05, pp. 86-9 DOI:0.3968/6634 ISSN 95-54X [Print] ISSN 95-5438 [Online] www.cscanada.net www.cscanada.org Calculating Method for the Axial
More informationLattice Boltzmann Simulation of One Particle Migrating in a Pulsating Flow in Microvessel
Commun. Theor. Phys. 56 (2011) 756 760 Vol. 56, No. 4, October 15, 2011 Lattice Boltzmann Simulation of One Particle Migrating in a Pulsating Flow in Microvessel QIU Bing ( ), 1, TAN Hui-Li ( Û), 2 and
More informationSlow viscous flow in a microchannel with similar and different superhydrophobic walls
Journal of Physics: Conference Series PAPER OPEN ACCESS Slow viscous flow in a microchannel with similar and different superhydrophobic walls To cite this article: A I Ageev and A N Osiptsov 2018 J. Phys.:
More informationDEM simulation of fracture process of inherently anisotropic rock under Brazilian test condition
Title DEM simulation of fracture process of inherently anisotropic rock under Brazilian test condition Author(s) Kwok, CY; Duan, K Citation The 49th US Rock Mechanics / Geomechanics Symposium, San Francisco,
More information3D dynamic crack propagation analysis with PDS-FEM
JAMSTEC-R IFREE Special Issue, November 2009 3D dynamic crack propagation analysis with PDS-FEM 1*, Hide Sakaguchi 1, Kenji Oguni 2, Muneo Hori 2 discretizations. Keywords +81-45-778-5972 lalith@jamstec.go.jp
More information3D simulations of an injection test done into an unsaturated porous and fractured limestone
3D simulations of an injection test done into an unsaturated porous and fractured limestone A. Thoraval *, Y. Guglielmi, F. Cappa INERIS, Ecole des Mines de Nancy, FRANCE *Corresponding author: Ecole des
More informationPLASTICITY FOR CRUSHABLE GRANULAR MATERIALS VIA DEM
Plasticity for crushable granular materials via DEM XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII E. Oñate, D.R.J. Owen, D. Peric and M. Chiumenti
More informationSTRESS DROP AS A RESULT OF SPLITTING, BRITTLE AND TRANSITIONAL FAULTING OF ROCK SAMPLES IN UNIAXIAL AND TRIAXIAL COMPRESSION TESTS
Studia Geotechnica et Mechanica, Vol. 37, No. 1, 2015 DOI: 10.1515/sgem-2015-0003 STRESS DROP AS A RESULT OF SPLITTING, BRITTLE AND TRANSITIONAL FAULTING OF ROCK SAMPLES IN UNIAXIAL AND TRIAXIAL COMPRESSION
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lectures & 3, 9/31 Aug 017 www.geosc.psu.edu/courses/geosc508 Discussion of Handin, JGR, 1969 and Chapter 1 Scholz, 00. Stress analysis and Mohr Circles Coulomb Failure
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationNumerical Simulation of Viscous Fingering Phenomenon in Immiscible Displacement of Two Fluids in Porous Media Using Lattice Boltzmann Method
nd Micro and Nano Flows Conference West London, UK, 1- September 009 Numerical Simulation of Viscous Fingering Phenomenon in Immiscible Displacement of Two Fluids in Porous Media Using Lattice Boltzmann
More informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationMechanics of Granular Matter
Mechanics of Granular Matter Mechanics of Granular Matter Qicheng Sun & Guangqian Wang Tsinghua University, Beijing, China Qicheng Sun & Guangqian Wang Tsinghua University, Beijing, China Published by
More informationSome New Developments in Rock Mechanics Research and Application
Some New Developments in Rock Mechanics Research and Application Jian Zhao Wuhan University, 30 October 014 Developments in Rock Mechanics Research and Application Research Micromechanics and critical
More informationSEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS
The Seventh Nobeyama Workshop on CFD: To the Memory of Prof. Kuwahara Sanjo Kaikan, the University of Tokyo, Japan, September. 23-24, 2009 SEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS Jaw-Yen
More informationFluid-soil multiphase flow simulation by an SPH-DEM coupled method
Fluid-soil multiphase flow simulation by an SPH-DEM coupled method *Kensuke Harasaki 1) and Mitsuteru Asai 2) 1), 2) Department of Civil and Structural Engineering, Kyushu University, 744 Motooka, Nishi-ku,
More informationDrag Force Simulations of Particle Agglomerates with the Lattice-Boltzmann Method
Drag Force Simulations of Particle Agglomerates with the Lattice-Boltzmann Method Christian Feichtinger, Nils Thuerey, Ulrich Ruede Christian Binder, Hans-Joachim Schmid, Wolfgang Peukert Friedrich-Alexander-Universität
More informationMaterial is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f
Material is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f Material undergoes plastic deformation when stress exceeds yield stress σ 0 Permanent strain results from
More informationParallel Simulations of Self-propelled Microorganisms
Parallel Simulations of Self-propelled Microorganisms K. Pickl a,b M. Hofmann c T. Preclik a H. Köstler a A.-S. Smith b,d U. Rüde a,b ParCo 2013, Munich a Lehrstuhl für Informatik 10 (Systemsimulation),
More informationModelling Progressive Failure with MPM
Modelling Progressive Failure with MPM A. Yerro, E. Alonso & N. Pinyol Department of Geotechnical Engineering and Geosciences, UPC, Barcelona, Spain ABSTRACT: In this work, the progressive failure phenomenon
More informationConstitutive model development for granular and porous materials and modelling of particulate processes
Constitutive model development for granular and porous materials and modelling of particulate processes Csaba Sinka ics4@le.ac.uk Department of Engineering Mechanics of Materials Group Geophysics modelling
More informationVORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS
THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1395-1399 1395 VORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS by Li-Zhong HUANG a* and De-Ming NIE b a State Key Laboratory
More informationNumerical Investigation on Shear Driven Cavity Flow by the Constrained Interpolated Profile Lattice Boltzmann Method
Numerical Investigation on Shear Driven Cavity Flow by the Constrained Interpolated Profile Lattice Boltzmann Method *C. S. NOR AZWADI, M. H. Al-MOLA and S. AGUS Department of Thermo-fluid Universiti Teknologi
More informationA multiscale framework for lubrication analysis of bearings with textured surface
A multiscale framework for lubrication analysis of bearings with textured surface *Leiming Gao 1), Gregory de Boer 2) and Rob Hewson 3) 1), 3) Aeronautics Department, Imperial College London, London, SW7
More information