In order to visualize the 3-D field of explosion and describe the complex physical phenomena of explosion, the 3-D data resulting from numerical simulation by 3-D multi-material in cell (MMIC), and the application of ...In order to visualize the 3-D field of explosion and describe the complex physical phenomena of explosion, the 3-D data resulting from numerical simulation by 3-D multi-material in cell (MMIC), and the application of volume visualization is explored, based on the characteristics of explosion and shock. Based on this, a visualization system for 3-D explosion--ViSC3D is designed. Approaches for the visualization of 3-D field of explosion are presented. The algorithm and the functions of ViSC3D are also presented. ViSC3D is thus a useful tool to observe and analyze either the full picture or the details of a 3-D field of explosion, that are difficult to observe and analyze directly. With ViSC3D, the field of explosion between the hill slopes is visualized. The cutaway views and 2-D slices are also given. The full picture and partial details of 3-D field of explosion can be observed clearly. Furthermore, ViSC3D can be used to visualize other similar 3-D data fields.展开更多
Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately ...Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately simulate three-dimensional （3-D） full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity-stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity Vx, Vy, and Vz for the same node in 3-D staggered-grid finite difference models by calculating the average value of Vy, and Vz of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.展开更多
The liquid flow on a single-pass sieve distillation tray is simulated with a three-dimensional computational fluid dynamics (CFD) program with the K-ε turbulence model. In the model, a source term SMi is formulatedin...The liquid flow on a single-pass sieve distillation tray is simulated with a three-dimensional computational fluid dynamics (CFD) program with the K-ε turbulence model. In the model, a source term SMi is formulatedin the Navier-Stokes equations to represent the interfacial momentum transfer and another term Sc is added to themass transfer equation as the source of interfacial mass transfer. The simulation provides the detailed informationof the three-dimensional distribution of liquid velocity on the tray, the circulation area and the concentration profilealong the height of liquid layer.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
Room and pillar sizes are key factors for safe mining and ore recovery in open-stope mining. To investigate the influence of room and pillar configurations on stope stability in highly fractured and weakened areas, an...Room and pillar sizes are key factors for safe mining and ore recovery in open-stope mining. To investigate the influence of room and pillar configurations on stope stability in highly fractured and weakened areas, an orthogonal design with two factors, three levels and nine runs was proposed, followed by three-dimensional numerical simulation using ANSYS and FLAC3~. Results show that surface settlement after excavation is concentrically ringed, and increases with the decrease of pillar width and distances to stope gobs. In the meantime, the ore-control fault at the ore-rock boundary and the fractured argillaceous dolomite with intercalated slate at the hanging wall deteriorate the roof settlement. Additionally, stope stability is challenged due to pillar rheological yield and stress concentration, and both are induced by redistribution of stress and plastic zones after mining. Following an objective function and a constraint function, room and pillar configuration with widths of 14 m and 16 m, respectively, is presented as the optimization for improving the ore recovery rate while maintaining a safe working environment.展开更多
Backfill is increasingly used in underground mines to reduce the surface impact from the wastes produced by the mining operations. But the main objectives of backfilling are to improve ground stability and reduce ore ...Backfill is increasingly used in underground mines to reduce the surface impact from the wastes produced by the mining operations. But the main objectives of backfilling are to improve ground stability and reduce ore dilution. To this end, the backfill in a stope must possess a minimum strength to remain self-standing during mining of an adjacent stope. This required strength is often estimated using a solution proposed by Mitchell and co-workers, which was based on a limit equilibrium analysis of a wedge exposed by the open face. In this paper, three dimensional numerical simulations have been performed to assess the behavior of the wedge model. A new limit equilibrium solution is proposed, based on the backfill displacements obtained from the simulations. Comparisons are made between the proposed solution and experimental and numerical modeling results. Compared with the previous solution, a better agreement is obtained between the new solution and experimental results for the required cohesion and factor of safety. For large scale(field) conditions, the results also show that the required strength obtained from the proposed solution corresponds quite well to the simulated backfill response.展开更多
文摘In order to visualize the 3-D field of explosion and describe the complex physical phenomena of explosion, the 3-D data resulting from numerical simulation by 3-D multi-material in cell (MMIC), and the application of volume visualization is explored, based on the characteristics of explosion and shock. Based on this, a visualization system for 3-D explosion--ViSC3D is designed. Approaches for the visualization of 3-D field of explosion are presented. The algorithm and the functions of ViSC3D are also presented. ViSC3D is thus a useful tool to observe and analyze either the full picture or the details of a 3-D field of explosion, that are difficult to observe and analyze directly. With ViSC3D, the field of explosion between the hill slopes is visualized. The cutaway views and 2-D slices are also given. The full picture and partial details of 3-D field of explosion can be observed clearly. Furthermore, ViSC3D can be used to visualize other similar 3-D data fields.
基金supported by National Natural Science Foundation of China(Nos.41204077,41372290,41572244,51034003,51174210,and 51304126)natural science foundation of Shandong Province(Nos.ZR2011EEZ002 and ZR2013EEQ019)State Key Research Development Program of China(No.2016YFC0600708-3)
文摘Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately simulate three-dimensional （3-D） full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity-stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity Vx, Vy, and Vz for the same node in 3-D staggered-grid finite difference models by calculating the average value of Vy, and Vz of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.
基金Supported by Visiting Scholar Foundation of Key Lab in Tianjin University.
文摘The liquid flow on a single-pass sieve distillation tray is simulated with a three-dimensional computational fluid dynamics (CFD) program with the K-ε turbulence model. In the model, a source term SMi is formulatedin the Navier-Stokes equations to represent the interfacial momentum transfer and another term Sc is added to themass transfer equation as the source of interfacial mass transfer. The simulation provides the detailed informationof the three-dimensional distribution of liquid velocity on the tray, the circulation area and the concentration profilealong the height of liquid layer.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Projects(50934002,51074013,51104100)supported by the National Natural Science Foundation of ChinaProject(IRT0950)supported by the Program for Changjiang Scholars and Innovative Research Team in University of China
文摘Room and pillar sizes are key factors for safe mining and ore recovery in open-stope mining. To investigate the influence of room and pillar configurations on stope stability in highly fractured and weakened areas, an orthogonal design with two factors, three levels and nine runs was proposed, followed by three-dimensional numerical simulation using ANSYS and FLAC3~. Results show that surface settlement after excavation is concentrically ringed, and increases with the decrease of pillar width and distances to stope gobs. In the meantime, the ore-control fault at the ore-rock boundary and the fractured argillaceous dolomite with intercalated slate at the hanging wall deteriorate the roof settlement. Additionally, stope stability is challenged due to pillar rheological yield and stress concentration, and both are induced by redistribution of stress and plastic zones after mining. Following an objective function and a constraint function, room and pillar configuration with widths of 14 m and 16 m, respectively, is presented as the optimization for improving the ore recovery rate while maintaining a safe working environment.
基金financial support of the Natural Sciences and Engineering Research Council (NSERC) of Canada and the partners of Research Institute on Mines and the Environment (RIME UQAT-Polytechnique http://rime-irme.ca)
文摘Backfill is increasingly used in underground mines to reduce the surface impact from the wastes produced by the mining operations. But the main objectives of backfilling are to improve ground stability and reduce ore dilution. To this end, the backfill in a stope must possess a minimum strength to remain self-standing during mining of an adjacent stope. This required strength is often estimated using a solution proposed by Mitchell and co-workers, which was based on a limit equilibrium analysis of a wedge exposed by the open face. In this paper, three dimensional numerical simulations have been performed to assess the behavior of the wedge model. A new limit equilibrium solution is proposed, based on the backfill displacements obtained from the simulations. Comparisons are made between the proposed solution and experimental and numerical modeling results. Compared with the previous solution, a better agreement is obtained between the new solution and experimental results for the required cohesion and factor of safety. For large scale(field) conditions, the results also show that the required strength obtained from the proposed solution corresponds quite well to the simulated backfill response.
