提出一种适用于非均匀离散颗粒体系的接触检测算法,时间复杂度为O(n),对内存的占用较小,对于大规模、分布密集、大粒径比的颗粒体系具有较高的效率。在NBS (no binary search)接触检测算法的基础上进行改进,提高了对于非均匀颗粒体系的...提出一种适用于非均匀离散颗粒体系的接触检测算法,时间复杂度为O(n),对内存的占用较小,对于大规模、分布密集、大粒径比的颗粒体系具有较高的效率。在NBS (no binary search)接触检测算法的基础上进行改进,提高了对于非均匀颗粒体系的适用性。采用对颗粒进行分组,对检测过程分步的策略提高了邻居检索精度;通过设置基于包围盒的初判环节减少了接触判断次数,解决了颗粒尺寸非均匀造成的效率下降问题。颗粒分组方式简单,易于调整,对于各种不同级配的颗粒体系具有广泛的适用性。通过算例验证了算法的性能,在非均匀颗粒体系中相对于NBS算法具有明显的优势。展开更多
This paper uses various mean-field approaches and the Monte Carlo simulation to calculate asymmetric simple exclusion processes with particles of arbitrary size in the successive defects system. In this system, the ho...This paper uses various mean-field approaches and the Monte Carlo simulation to calculate asymmetric simple exclusion processes with particles of arbitrary size in the successive defects system. In this system, the hopping probability p （p 〈 1） and the size d of particles are not constant, Through theoretical calculation and computer simulation, it obtains the exact theoretical results and finds that the theoretical results are in agreement with the computer simulation. These results are helpful in analysing the effect of traffic with different hopping probabilities p and sizes d of particle.展开更多
Compared with vertical and horizontal wells, the solution and computation of transient pressure responses of slanted wells are more complex. Vertical and horizontal wells are both simplified cases of slanted wells at ...Compared with vertical and horizontal wells, the solution and computation of transient pressure responses of slanted wells are more complex. Vertical and horizontal wells are both simplified cases of slanted wells at particular inclination, so the model for slanted wells is more general and more complex than other models for vertical and horizontal wells. Many authors have studied unsteady-state flow of fluids in slanted wells and various solutions have been proposed. However, until now, few of the published results pertain to the computational efficiency. Whether in the time domain or in the Laplace domain, the computation of integration of complex functions is necessary in obtaining pressure responses of slanted wells, while the computation of the integration is complex and time-consuming. To obtain a perfect type curve the computation time is unacceptable even with an aid of high-speed computers. The purpose of this paper is to present an efficient algorithm to compute transient pressure distributions caused by slanted wells in reservoirs. Based on rigorous derivation, the transient pressure solution for slanted wells of any inclination angle is presented. Assuming an infinite-conductivity wellbore, the location of the equivalent-pressure point is determined. More importantly, according to the characteristics of the integrand in a transient pressure solution for slanted wells, the whole integral interval is partitioned into several small integral intervals, and then the method of variable substitution and the variable step-size piecewise numerical integration are employed. The amount of computation is significantly reduced and the computational efficiency is greatly improved. The algorithm proposed in this paper thoroughly solved the difficulty in the efficient and high-speed computation of transient pressure distribution of slanted wells with any inclination angle.展开更多
The microstructures of a composite determine its macroscopic properties. In this study, microstructures with particles of arbitrary shapes and sizes are constructed by using several developed fractal geometry algorith...The microstructures of a composite determine its macroscopic properties. In this study, microstructures with particles of arbitrary shapes and sizes are constructed by using several developed fractal geometry algorithms implemented in MATLAB. A two-dimensional (2D) quadrilateral fractal geometry algorithm is developed based on the modified Sierpinski carpet algorithm. Square-, rectangle-, circle-, and ellipse-based microstructure constructions are special cases of the 2D quadrilateral fractal geometry algorithm. Moreover, a three-dimensional (3D) random hexahedron geometry algorithm is developed according to the Menger sponge algorithm. Cube-and sphere-based mi-crostructure constructions are special cases of the 3D hexahedron fractal geometry algo-rithm. The polydispersities of fractal shapes and random fractal sub-units demonstrate significant enhancements compared to those obtained by the original algorithms. In ad-dition, the 2D and 3D algorithms mentioned in this article can be combined according to the actual microstructures. The verification results also demonstrate the practicability of these algorithms. The developed algorithms open up new avenues for the constructions of microstructures, which can be embedded into commercial finite element method soft-wares.展开更多
文摘提出一种适用于非均匀离散颗粒体系的接触检测算法,时间复杂度为O(n),对内存的占用较小,对于大规模、分布密集、大粒径比的颗粒体系具有较高的效率。在NBS (no binary search)接触检测算法的基础上进行改进,提高了对于非均匀颗粒体系的适用性。采用对颗粒进行分组,对检测过程分步的策略提高了邻居检索精度;通过设置基于包围盒的初判环节减少了接触判断次数,解决了颗粒尺寸非均匀造成的效率下降问题。颗粒分组方式简单,易于调整,对于各种不同级配的颗粒体系具有广泛的适用性。通过算例验证了算法的性能,在非均匀颗粒体系中相对于NBS算法具有明显的优势。
基金Project supported by the National Key Fundamental Research and Development Project of China (Grant No 2005CB724206)
文摘This paper uses various mean-field approaches and the Monte Carlo simulation to calculate asymmetric simple exclusion processes with particles of arbitrary size in the successive defects system. In this system, the hopping probability p （p 〈 1） and the size d of particles are not constant, Through theoretical calculation and computer simulation, it obtains the exact theoretical results and finds that the theoretical results are in agreement with the computer simulation. These results are helpful in analysing the effect of traffic with different hopping probabilities p and sizes d of particle.
基金financial support from the special fund of China’s central government for the development of local colleges and universities―the project of national first-level discipline in Oil and Gas Engineering, the National Science Fund for Distinguished Young Scholars of China (Grant No. 51125019)the National Program on Key fundamental Research Project (973 Program, Grant No. 2011CB201005)
文摘Compared with vertical and horizontal wells, the solution and computation of transient pressure responses of slanted wells are more complex. Vertical and horizontal wells are both simplified cases of slanted wells at particular inclination, so the model for slanted wells is more general and more complex than other models for vertical and horizontal wells. Many authors have studied unsteady-state flow of fluids in slanted wells and various solutions have been proposed. However, until now, few of the published results pertain to the computational efficiency. Whether in the time domain or in the Laplace domain, the computation of integration of complex functions is necessary in obtaining pressure responses of slanted wells, while the computation of the integration is complex and time-consuming. To obtain a perfect type curve the computation time is unacceptable even with an aid of high-speed computers. The purpose of this paper is to present an efficient algorithm to compute transient pressure distributions caused by slanted wells in reservoirs. Based on rigorous derivation, the transient pressure solution for slanted wells of any inclination angle is presented. Assuming an infinite-conductivity wellbore, the location of the equivalent-pressure point is determined. More importantly, according to the characteristics of the integrand in a transient pressure solution for slanted wells, the whole integral interval is partitioned into several small integral intervals, and then the method of variable substitution and the variable step-size piecewise numerical integration are employed. The amount of computation is significantly reduced and the computational efficiency is greatly improved. The algorithm proposed in this paper thoroughly solved the difficulty in the efficient and high-speed computation of transient pressure distribution of slanted wells with any inclination angle.
基金Project supported by the National Natural Science Foundation of China(Nos.11972218 and11472165)
文摘The microstructures of a composite determine its macroscopic properties. In this study, microstructures with particles of arbitrary shapes and sizes are constructed by using several developed fractal geometry algorithms implemented in MATLAB. A two-dimensional (2D) quadrilateral fractal geometry algorithm is developed based on the modified Sierpinski carpet algorithm. Square-, rectangle-, circle-, and ellipse-based microstructure constructions are special cases of the 2D quadrilateral fractal geometry algorithm. Moreover, a three-dimensional (3D) random hexahedron geometry algorithm is developed according to the Menger sponge algorithm. Cube-and sphere-based mi-crostructure constructions are special cases of the 3D hexahedron fractal geometry algo-rithm. The polydispersities of fractal shapes and random fractal sub-units demonstrate significant enhancements compared to those obtained by the original algorithms. In ad-dition, the 2D and 3D algorithms mentioned in this article can be combined according to the actual microstructures. The verification results also demonstrate the practicability of these algorithms. The developed algorithms open up new avenues for the constructions of microstructures, which can be embedded into commercial finite element method soft-wares.
