In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through poro...In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.展开更多
A numerical method is presented that simulates 3D explosive field problems. A code MMIC3D using this method can be used to simulate the propagation and reflected effects of all kinds of rigid boundaries to shock waves...A numerical method is presented that simulates 3D explosive field problems. A code MMIC3D using this method can be used to simulate the propagation and reflected effects of all kinds of rigid boundaries to shock waves produced by an explosive source. These numerical results indicate that the code MMIC3D has the ability in computing cases such as 3D shock waves produced by air explosion, vortex region of the shock wave, the Mach wave, and reflected waves behind rigid boundaries.展开更多
In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numer...In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases： The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver（used for solving the discretized governing equations in the computational domain） upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.展开更多
The developing processes of stress and deformation fields of a protected layer after mining an upper-protective layer with a bow pseudo-incline technique were simulated to locate the protection region. The pressure re...The developing processes of stress and deformation fields of a protected layer after mining an upper-protective layer with a bow pseudo-incline technique were simulated to locate the protection region. The pressure relief of the protected layer was analyzed after mining the upper-protective layer. The pressure relief angle along the strike and incline were located according to the roles of protection of the deformation and stress pressure-relief of the protective layer after mining. This results show that the upper-protective layer with the bow pseudo-incline technique have an upper and downside pressure relief angle of 85 and 68 degrees respectively; the distribution of strike pressure relief angles along the pseudo-incline working face is uneven and their values range from 38.3 to 51 degrees. The pressure relief angle of the inclined middle location was the largest. The distribution of the protection region of the upper-protective layer with the bow pseudo-incline teelmique located by practical tests and numerical simulation is essentially consistent, compared with the results obtained by these methods.展开更多
We employ the parallel computing technology to study numerically the three-dimensional structure of quantized vortices of Bose-Einstein condensates, For anisotropic cases, the bending process of vortices is described ...We employ the parallel computing technology to study numerically the three-dimensional structure of quantized vortices of Bose-Einstein condensates, For anisotropic cases, the bending process of vortices is described in detail by the decrease of Gross-Pitaevskii energy. A completely straight vortex and the steady and symmetrical multiple-vortex configurations are obtained. We analyse the effect of initial conditions and angular velocity on the number and shape of vortices.展开更多
Three-dimension isothermal flows of polymer melt in the kneading blocks of triangularly-arranged and parallelly-arranged intermeshing co-rotating three-screw extruders are simulated using the finite element package PO...Three-dimension isothermal flows of polymer melt in the kneading blocks of triangularly-arranged and parallelly-arranged intermeshing co-rotating three-screw extruders are simulated using the finite element package POLYFLOW. Based on the velocity fields calculated, the particle trajectories in both machines are visualized using particle tracking technique. The numerical results indicate that the flow patterns in three-screw extruders are similar to those in twin-screw extruders. The triangularly-arranged three-screw extruder has the largest pumping capacity and also the highest extrusion stability in terms of flowrate fluctuation with screw rotation. The instantaneous mixing and cumulative residence time distribution （RTD） characteristics are also analyzed and compared with traditional intermeshing co-rotating twin-screw extruders. It is shown that the start section of the cumulative RTD curve for the triangularly-arranged machine has a small shoulder, which is attributed to the faster flow in the central region of this type of extruder.展开更多
For the engineering geology conditions of bad mine roadway roof and floor lithology in extremely weak cemented strata, the best section shape of the roadway is determined from the study of tunnel surrounding rock disp...For the engineering geology conditions of bad mine roadway roof and floor lithology in extremely weak cemented strata, the best section shape of the roadway is determined from the study of tunnel surrounding rock displacement, plastic zone and stress distribution in rectangular, circle arch and arch wall sections, respectively. Based on the mining depth and thickness of the coal seam, roadway support technology solutions with different buried depth and thickness of coal seam are proposed. Support schemes are amended and optimized in time through monitoring data of the deformation of roadway, roof separation, l-beam bracket, bolt and anchor cable force to ensure the long-term stability and security of the roadway surrounding rock and support structure. The monitoring results show that mine roadway support schemes for different buried depth and section can be adapted to the characteristics of ground pressure and deformation of the surrounding rock in different depth well, effectively control the roadway surrounding rock deformation and the floor heave and guarantee the safety of construction and basic stability of surrounding rock and support structure.展开更多
针对混凝土板桩在干粗砂地基中的静压贯入问题,为了明确贯入过程中的挤土效应及贯入阻力的发展规律,利用离散元数值模拟软件EDEM(Engineering Discrete Element Method)对干粗砂中混凝土板桩静压贯入过程进行了三维数值模拟。为此,提出...针对混凝土板桩在干粗砂地基中的静压贯入问题,为了明确贯入过程中的挤土效应及贯入阻力的发展规律,利用离散元数值模拟软件EDEM(Engineering Discrete Element Method)对干粗砂中混凝土板桩静压贯入过程进行了三维数值模拟。为此,提出了适用于类似试验的参数标定方法,分析了贯入过程中颗粒的运动特征。针对板桩结构的特点,对长边和短边方向的挤土效应进行了对比,并以短边为例,分析了不同高度处的挤土效应。最后,根据颗粒与桩结构的接触力,提取了桩端阻力和桩侧摩阻力的值,分析了两者在贯入过程中的变化规律。结果显示:在贯入过程中,颗粒的运动呈现"涡流状",并随贯入过程演化;板桩沿短边方向的表面隆起量约为沿长边方向的隆起量的2倍;随着长宽比的增加,板桩结构在短边方向和长边方向的隆起值相应增加,且短边方向的隆起值与长边方向的隆起值的比值与相应板桩的长宽比近似相等;随着贯入的进行,桩侧摩阻力所占的比重逐渐增加,在完全贯入后,侧摩阻力值约占桩端阻力的1/25;不同长宽比的板桩结构,其桩端阻力和桩侧摩阻力的单位值近似相等,并与准二维计算结果相近。展开更多
基金Projects 50534090 and 50674090 supported by the National Natural Science Foundation of China2005CB221503 by the National Key Basic ResearchDevelopment Program (973 Program)
文摘In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.
文摘A numerical method is presented that simulates 3D explosive field problems. A code MMIC3D using this method can be used to simulate the propagation and reflected effects of all kinds of rigid boundaries to shock waves produced by an explosive source. These numerical results indicate that the code MMIC3D has the ability in computing cases such as 3D shock waves produced by air explosion, vortex region of the shock wave, the Mach wave, and reflected waves behind rigid boundaries.
文摘In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases： The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver（used for solving the discretized governing equations in the computational domain） upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.
基金Projects PLN0610 supported by the Open Fund of State Key Lab of Oil and Gas Reservoir Geology and Exploitation (Southwest Petroleum University)HKLGF200706 by the Opening Project of Henan Key Laboratory of Coal Mine Methane and Fire Prevention+3 种基金50334060, 50474025 and 50774106 by the National Natural Science Foundation of China2005CB221502 by the National Basic Research Program of China50621403 by the Natural Science Innova-tion Group Foundation of ChinaCSTC, 2006BB7147, 2006AA7002 by the Natural Science Foundation of Chongqing
文摘The developing processes of stress and deformation fields of a protected layer after mining an upper-protective layer with a bow pseudo-incline technique were simulated to locate the protection region. The pressure relief of the protected layer was analyzed after mining the upper-protective layer. The pressure relief angle along the strike and incline were located according to the roles of protection of the deformation and stress pressure-relief of the protective layer after mining. This results show that the upper-protective layer with the bow pseudo-incline technique have an upper and downside pressure relief angle of 85 and 68 degrees respectively; the distribution of strike pressure relief angles along the pseudo-incline working face is uneven and their values range from 38.3 to 51 degrees. The pressure relief angle of the inclined middle location was the largest. The distribution of the protection region of the upper-protective layer with the bow pseudo-incline teelmique located by practical tests and numerical simulation is essentially consistent, compared with the results obtained by these methods.
基金Project supported partly by the National Natural Science Foundation of China （Grant Nos 10301034 and 40574069）, The authors thank Professor Du Q very much for his important discussions.
文摘We employ the parallel computing technology to study numerically the three-dimensional structure of quantized vortices of Bose-Einstein condensates, For anisotropic cases, the bending process of vortices is described in detail by the decrease of Gross-Pitaevskii energy. A completely straight vortex and the steady and symmetrical multiple-vortex configurations are obtained. We analyse the effect of initial conditions and angular velocity on the number and shape of vortices.
文摘Three-dimension isothermal flows of polymer melt in the kneading blocks of triangularly-arranged and parallelly-arranged intermeshing co-rotating three-screw extruders are simulated using the finite element package POLYFLOW. Based on the velocity fields calculated, the particle trajectories in both machines are visualized using particle tracking technique. The numerical results indicate that the flow patterns in three-screw extruders are similar to those in twin-screw extruders. The triangularly-arranged three-screw extruder has the largest pumping capacity and also the highest extrusion stability in terms of flowrate fluctuation with screw rotation. The instantaneous mixing and cumulative residence time distribution （RTD） characteristics are also analyzed and compared with traditional intermeshing co-rotating twin-screw extruders. It is shown that the start section of the cumulative RTD curve for the triangularly-arranged machine has a small shoulder, which is attributed to the faster flow in the central region of this type of extruder.
基金Financial support for this work, provided by the Major Program of the National Natural Science Foundation of China (Nos. 51174196 and 51204168)the Program for New Century Excellent Talents in University by Ministry of Education of China (No. NCET-07-0519)
文摘For the engineering geology conditions of bad mine roadway roof and floor lithology in extremely weak cemented strata, the best section shape of the roadway is determined from the study of tunnel surrounding rock displacement, plastic zone and stress distribution in rectangular, circle arch and arch wall sections, respectively. Based on the mining depth and thickness of the coal seam, roadway support technology solutions with different buried depth and thickness of coal seam are proposed. Support schemes are amended and optimized in time through monitoring data of the deformation of roadway, roof separation, l-beam bracket, bolt and anchor cable force to ensure the long-term stability and security of the roadway surrounding rock and support structure. The monitoring results show that mine roadway support schemes for different buried depth and section can be adapted to the characteristics of ground pressure and deformation of the surrounding rock in different depth well, effectively control the roadway surrounding rock deformation and the floor heave and guarantee the safety of construction and basic stability of surrounding rock and support structure.
文摘针对混凝土板桩在干粗砂地基中的静压贯入问题,为了明确贯入过程中的挤土效应及贯入阻力的发展规律,利用离散元数值模拟软件EDEM(Engineering Discrete Element Method)对干粗砂中混凝土板桩静压贯入过程进行了三维数值模拟。为此,提出了适用于类似试验的参数标定方法,分析了贯入过程中颗粒的运动特征。针对板桩结构的特点,对长边和短边方向的挤土效应进行了对比,并以短边为例,分析了不同高度处的挤土效应。最后,根据颗粒与桩结构的接触力,提取了桩端阻力和桩侧摩阻力的值,分析了两者在贯入过程中的变化规律。结果显示:在贯入过程中,颗粒的运动呈现"涡流状",并随贯入过程演化;板桩沿短边方向的表面隆起量约为沿长边方向的隆起量的2倍;随着长宽比的增加,板桩结构在短边方向和长边方向的隆起值相应增加,且短边方向的隆起值与长边方向的隆起值的比值与相应板桩的长宽比近似相等;随着贯入的进行,桩侧摩阻力所占的比重逐渐增加,在完全贯入后,侧摩阻力值约占桩端阻力的1/25;不同长宽比的板桩结构,其桩端阻力和桩侧摩阻力的单位值近似相等,并与准二维计算结果相近。
