Based on water inrush accident of 1841 working face of Desheng Coal Mine in Wu'an, Hebei province, China, an evaluation model of hydrodynamic characteristics of the project is set up and simulated using Matlab. It...Based on water inrush accident of 1841 working face of Desheng Coal Mine in Wu'an, Hebei province, China, an evaluation model of hydrodynamic characteristics of the project is set up and simulated using Matlab. It is assumed that the pipe flow would transform into seepage flow when the aggregates are plugged into the water inrush channel and the seepage flow would disappear along with grouting process. The simulation results show that the flow velocity will increase with an increase in height of aggregates accumulation body during the aggregates filling process; the maximum seepage velocity occurs on the top of plugging zone; and the water flow decreases with increasing plugging height of water inrush channel. Finally, the field construction results show that the water inrush channel can be plugged effectively by the compacted body prepared with aggregate and cement slurry.展开更多
The evolution of two spanwise-aligned low-speed streaks in a wall turbulent flow, triggered by the instability of the subharmonic varicose (SV) mode, is studied by a direct numerical simulation (DNS) method in a s...The evolution of two spanwise-aligned low-speed streaks in a wall turbulent flow, triggered by the instability of the subharmonic varicose (SV) mode, is studied by a direct numerical simulation (DNS) method in a small spatial-periodic channel. The results show that the SV low-speed streaks are self-sustained at the early stage, and then transform into subharmonic sinuous (SS) low-speed streaks. Initially, the streamwise vortex sheets are formed by shearing, and then evolve into zigzag vortex sheets due to the mutual induction. As the intensification of the SV low-speed streaks becomes prominent, the tilted streamwise vortex tubes and the V-like streamwise vortex tubes can be formed simultaneously by increasing +~. When the SV low-speed streaks break down, new zigzag streamwise vortices will be generated, thus giving birth to the next sustaining cycle of the SV low-speed streaks. When the second breakdown happens, new secondary V-like streamwise vortices instead of zigzag streamwise vortices will be generated. Because of the sweep motion of the fluid induced by the secondary V-like streamwise vortices, each decayed low-speed streak can be divided into two parts, and each part combines with the part of another streak, finally leading to the formation of SS low-speed streaks.展开更多
A coupled numerical method for the direct numerical simulation of particle-fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by t...A coupled numerical method for the direct numerical simulation of particle-fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by the time-driven hard-sphere model, while the hydrodynamic equations governing fluid flow are solved by the lattice Boltzmann method (LBM), Particle-fluid coupling is realized by an immersed boundary method (IBM), which considers the effect of boundary on surrounding fluid as a restoring force added to the governing equations of the fluid. The proposed scheme is validated in the classical flow-around-cylinder simulations, and preliminary application of this scheme to fluidization is reported, demonstrating it to be a promising computational strategy for better understanding complex behavior in particle-fluid systems.展开更多
Numerical simulation based on computational fluid dynamics (CFD) is a useful approach for quantitatively investigating the underlying thermal-mechanical conditions during FSW, such as temperature field and material ...Numerical simulation based on computational fluid dynamics (CFD) is a useful approach for quantitatively investigating the underlying thermal-mechanical conditions during FSW, such as temperature field and material deformation field. One of the critical issues in CFD simulation of FSW is the use of the frictional boundary condition, which represents the friction between the welding tool and the workpiece in the numerical models. In this study, three-dimensional numerical simulation is conducted to analyze the heat transfer and plastic deformation behaviors during the FSW of AA2024. For comparison purposes, both the boundary velocity (BV) models and the boundary shear stress (BSS) models are employed in order to assess their performances in predicting the temperature and material deformation in FSW. It is interesting to note that different boundary conditions yield similar predictions on temperature, but quite different predictions on material deformation. The numerical predictions are compared with the experimental results. The predicted deformation zone geometry by the BSS model is consistent with the experimental results while there is large difference between the predictions by the BV models and the experimental measurements. The fact that the BSS model yields more reasonable predictions on the deformation zone geometry is attributed to its capacity to automatically adjust the contact state at the tool/workpiece interface. Based on the favorable predictions on both the temperature field and the material deformation field, the BSS model is suggested to have a better performance in numerical simulation of FSW than the BV model.展开更多
The lattice thermal conductivity of boron nitride nanoribbon(BNNR) is calculated by using equilibrium molecular dynamics(EMD) simulation method. The Green–Kubo relation derived from linear response theory is used...The lattice thermal conductivity of boron nitride nanoribbon(BNNR) is calculated by using equilibrium molecular dynamics(EMD) simulation method. The Green–Kubo relation derived from linear response theory is used to acquire the thermal conductivity from heat current auto-correlation function(HCACF). HCACF of the selected BNNR system shows a tendency of a very fast decay and then be followed by a very slow decay process,finally,approaching zero approximately within 3 ps. The convergence of lattice thermal conductivity demonstrates that the thermal conductivity of BNNR can be simulated by EMD simulation using several thousands of atoms with periodic boundary conditions. The results show that BNNR exhibit lower thermal conductivity than that of boron nitride(BN) monolayer,which indicates that phonons boundary scatting significantly suppresses the phonons transport in BNNR. Vacancies in BNNR greatly affect the lattice thermal conductivity,in detail,only 1% concentration of vacancies in BNNR induce a 60% reduction of the lattice thermal conductivity at room temperature.展开更多
The lattice Boltzmann method (LBM) has gained increasing popularity in the last two decades as an alternative numerical approach for solving fluid flow problems. One of the most active research areas in the LBM is i...The lattice Boltzmann method (LBM) has gained increasing popularity in the last two decades as an alternative numerical approach for solving fluid flow problems. One of the most active research areas in the LBM is its application in particle-fluid systems, where the advantage of the LBM in efficiency and parallel scalability has made it superior to many other direct numerical simulation (DNS) techniques. This article intends to provide a brief review of the application of the LBM in particle-fluid systems. The numerical techniques in the LBM pertaining to simulations of particles are discussed, with emphasis on the advanced treatment for boundary conditions on the particle-fluid interface. Other numerical issues, such as the effect of the internal fluid, are also briefly described. Additionally, recent efforts in using the LBM to obtain closures for particle-fluid drag force are also reviewed.展开更多
基金Financial support for this work, provided by the National Natural Science Foundation of China (Nos. 41072031, 40172119)the Natural Science Foundation of Hebei Province of China(No. D2012402008)
文摘Based on water inrush accident of 1841 working face of Desheng Coal Mine in Wu'an, Hebei province, China, an evaluation model of hydrodynamic characteristics of the project is set up and simulated using Matlab. It is assumed that the pipe flow would transform into seepage flow when the aggregates are plugged into the water inrush channel and the seepage flow would disappear along with grouting process. The simulation results show that the flow velocity will increase with an increase in height of aggregates accumulation body during the aggregates filling process; the maximum seepage velocity occurs on the top of plugging zone; and the water flow decreases with increasing plugging height of water inrush channel. Finally, the field construction results show that the water inrush channel can be plugged effectively by the compacted body prepared with aggregate and cement slurry.
基金supported by the National Natural Science Foundation of China(Nos.11372140 and11202102)the Innovation Project for College Graduates of Jiangsu Province(No.CXZZ13-0189)
文摘The evolution of two spanwise-aligned low-speed streaks in a wall turbulent flow, triggered by the instability of the subharmonic varicose (SV) mode, is studied by a direct numerical simulation (DNS) method in a small spatial-periodic channel. The results show that the SV low-speed streaks are self-sustained at the early stage, and then transform into subharmonic sinuous (SS) low-speed streaks. Initially, the streamwise vortex sheets are formed by shearing, and then evolve into zigzag vortex sheets due to the mutual induction. As the intensification of the SV low-speed streaks becomes prominent, the tilted streamwise vortex tubes and the V-like streamwise vortex tubes can be formed simultaneously by increasing +~. When the SV low-speed streaks break down, new zigzag streamwise vortices will be generated, thus giving birth to the next sustaining cycle of the SV low-speed streaks. When the second breakdown happens, new secondary V-like streamwise vortices instead of zigzag streamwise vortices will be generated. Because of the sweep motion of the fluid induced by the secondary V-like streamwise vortices, each decayed low-speed streak can be divided into two parts, and each part combines with the part of another streak, finally leading to the formation of SS low-speed streaks.
基金sponsored by Ministry of Finance under the grant ZDYZ2008-2National Key Science and Technology Project under the grant 2008ZX05014-003-006HZthe Chinese Academy of Sciences under the grant KGCX2-YW-124
文摘A coupled numerical method for the direct numerical simulation of particle-fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by the time-driven hard-sphere model, while the hydrodynamic equations governing fluid flow are solved by the lattice Boltzmann method (LBM), Particle-fluid coupling is realized by an immersed boundary method (IBM), which considers the effect of boundary on surrounding fluid as a restoring force added to the governing equations of the fluid. The proposed scheme is validated in the classical flow-around-cylinder simulations, and preliminary application of this scheme to fluidization is reported, demonstrating it to be a promising computational strategy for better understanding complex behavior in particle-fluid systems.
基金supported by the National Natural Science Foundation of China(Grant No.51375259 and Grant No.51705280)the Ministry of Science and Technology of China(Grant No.2012ZX04012-011)+1 种基金Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund(the second phase,Grant No.U1501501)the Tsinghua National Laboratory for Information Science and Technology
文摘Numerical simulation based on computational fluid dynamics (CFD) is a useful approach for quantitatively investigating the underlying thermal-mechanical conditions during FSW, such as temperature field and material deformation field. One of the critical issues in CFD simulation of FSW is the use of the frictional boundary condition, which represents the friction between the welding tool and the workpiece in the numerical models. In this study, three-dimensional numerical simulation is conducted to analyze the heat transfer and plastic deformation behaviors during the FSW of AA2024. For comparison purposes, both the boundary velocity (BV) models and the boundary shear stress (BSS) models are employed in order to assess their performances in predicting the temperature and material deformation in FSW. It is interesting to note that different boundary conditions yield similar predictions on temperature, but quite different predictions on material deformation. The numerical predictions are compared with the experimental results. The predicted deformation zone geometry by the BSS model is consistent with the experimental results while there is large difference between the predictions by the BV models and the experimental measurements. The fact that the BSS model yields more reasonable predictions on the deformation zone geometry is attributed to its capacity to automatically adjust the contact state at the tool/workpiece interface. Based on the favorable predictions on both the temperature field and the material deformation field, the BSS model is suggested to have a better performance in numerical simulation of FSW than the BV model.
基金Supported by the Natural Science Foundation of Hubei Province(2014CFB610)the Excellent Young Innovation Team Project of Hubei Province(T201429)
文摘The lattice thermal conductivity of boron nitride nanoribbon(BNNR) is calculated by using equilibrium molecular dynamics(EMD) simulation method. The Green–Kubo relation derived from linear response theory is used to acquire the thermal conductivity from heat current auto-correlation function(HCACF). HCACF of the selected BNNR system shows a tendency of a very fast decay and then be followed by a very slow decay process,finally,approaching zero approximately within 3 ps. The convergence of lattice thermal conductivity demonstrates that the thermal conductivity of BNNR can be simulated by EMD simulation using several thousands of atoms with periodic boundary conditions. The results show that BNNR exhibit lower thermal conductivity than that of boron nitride(BN) monolayer,which indicates that phonons boundary scatting significantly suppresses the phonons transport in BNNR. Vacancies in BNNR greatly affect the lattice thermal conductivity,in detail,only 1% concentration of vacancies in BNNR induce a 60% reduction of the lattice thermal conductivity at room temperature.
文摘The lattice Boltzmann method (LBM) has gained increasing popularity in the last two decades as an alternative numerical approach for solving fluid flow problems. One of the most active research areas in the LBM is its application in particle-fluid systems, where the advantage of the LBM in efficiency and parallel scalability has made it superior to many other direct numerical simulation (DNS) techniques. This article intends to provide a brief review of the application of the LBM in particle-fluid systems. The numerical techniques in the LBM pertaining to simulations of particles are discussed, with emphasis on the advanced treatment for boundary conditions on the particle-fluid interface. Other numerical issues, such as the effect of the internal fluid, are also briefly described. Additionally, recent efforts in using the LBM to obtain closures for particle-fluid drag force are also reviewed.