In this mini-review we summarize the progress of Lattice Boltzmann （LB） modeling and simulating compressible flows in our group in recent years. Main contents include （i） Single-Relaxation-Time （SRT） LB model su...In this mini-review we summarize the progress of Lattice Boltzmann （LB） modeling and simulating compressible flows in our group in recent years. Main contents include （i） Single-Relaxation-Time （SRT） LB model supplemented by additional viscosity, （ii） Multiple-Relaxation-Time （MRT） LB model, and （iii） LB study on hydrodynamic instabilities. The former two belong to improvements of physical modeling and the third belongs to simulation or application. The SRT-LB model sup- plemented by additional viscosity keeps the original framework of Lattice Bhatnagar-Gross Krook （LBGK）. So, it is easier and more convenient for previous SRT-LB users. The MRT-LB is a com- pletely new framework for physical modeling. It significantly extends the range of LB applications. The cost is longer computational time. The developed SRT-LB and MRT-LB are complementary from the sides of convenience and applicability.展开更多
The use of X-ray microtomographic （XMT） methods in analysing particulate systems has expanded rapidly in recent years with the availability of affordable desk-top apparatus. This review presents a summary of the maj...The use of X-ray microtomographic （XMT） methods in analysing particulate systems has expanded rapidly in recent years with the availability of affordable desk-top apparatus. This review presents a summary of the major applications in which computer simulations are explicitly coupled with XMT in the area of granular and porous materials. We envisage two main ways of establishing the coupling between both techniques, based on the transference or exchange of information by using physical or geometrical parameters （i.e. a parametric link through fitting to a process model） or through the direct use of3D XMT digital images （i.e. comparing image pixels and features directly）. Examples of coupled applications are shown for the study of transport properties of rocks, particle packing, mechanical loading and sintering. Often, the link between XMT and computer simulations is based on visual comparisons and we conclude that the use of quantitative parameters such as the number of interparticle contacts, force networks or granule shape to link both techniques is still underrepresented in the literature. Strategies to provide a more robust and quantitative approach to optimise the information obtained from such tomography analyses are proposed. 2010 Chinese Society of Particuology and Institute of Process Engineering, Chinese Academy of Sciences. Published by Elsevier B.V. All rights reserved.展开更多
This article deals with the solute transport in a single fracture with the combination of the Lattice Boltzmann Method （LBM） and Modified Moment Propagation （MMP） method, and this mixed method is proved to have se...This article deals with the solute transport in a single fracture with the combination of the Lattice Boltzmann Method （LBM） and Modified Moment Propagation （MMP） method, and this mixed method is proved to have several advantages over the LBM and Moment Propagation （MP） mixed method which leads to negative concentrations under some conditions in computation. The disadvantage of LBM/MP has been overcome to a certain extent. Also, this work presents an LBM solution of modeling single fractures with uniformly or randomly distributed grains, which can provide a new path of applying the LBM in solute transport simulation in fractures.展开更多
A supercomputer with 1.0 Petaflops peak performance in single precision, designed and established by Institute of Process Engineering, Chinese Academy of Sciences, is introduced in this brief communication. A designin...A supercomputer with 1.0 Petaflops peak performance in single precision, designed and established by Institute of Process Engineering, Chinese Academy of Sciences, is introduced in this brief communication. A designing philosophy utilizing the similarity between hardware, software and the problems to be solved is embodied, based on the multi-scale method and discrete simulation approaches developed at Institute of Process Engineering （IPE） and implemented in a graphic processing unit （GPU）-based hybrid computing mode. The preliminary applications of this machine in areas of multi-phase flow, molecular dynamics and so on are reported, demonstrating the supercomputer as a paradigm of green computation in new architecture.展开更多
How to solve the hypersonic aerothermodynamics around large-scale uncontrolled spacecraft during falling disintegrated process from outer space to earth,is the key to resolve the problems of the uncontrolled Tiangong-...How to solve the hypersonic aerothermodynamics around large-scale uncontrolled spacecraft during falling disintegrated process from outer space to earth,is the key to resolve the problems of the uncontrolled Tiangong-No.1 spacecraft reentry crash.To study aerodynamics of spacecraft reentry covering various flow regimes,a Gas-Kinetic Unified Algorithm(GKUA)has been presented by computable modeling of the collision integral of the Boltzmann equation over tens of years.On this basis,the rotational and vibrational energy modes are considered as the independent variables of the gas molecular velocity distribution function,a kind of Boltzmann model equation involving in internal energy excitation is presented by decomposing the collision term of the Boltzmann equation into elastic and inelastic collision terms.Then,the gas-kinetic numerical scheme is constructed to capture the time evolution of the discretized velocity distribution functions by developing the discrete velocity ordinate method and numerical quadrature technique.The unified algorithm of the Boltzmann model equation involving thermodynamics non-equilibrium effect is presented for the whole range of flow regimes.The gas-kinetic massive parallel computing strategy is developed to solve the hypersonic aerothermodynamics with the processor cores 500~45,000 at least 80%parallel efficiency.To validate the accuracy of the GKUA,the hypersonic flows are simulated including the reentry Tiangong-1 spacecraft shape with the wide range of Knudsen numbers of 220~0.00005 by the comparison of the related results from the DSMC and N-S coupled methods,and the low-density tunnel experiment etc.For uncontrolling spacecraft falling problem,the finite-element algorithm for dynamic thermalforce coupling response is presented,and the unified simulation of the thermal structural response and the hypersonic flow field is tested on the Tiangong-1 shape under reentry aerodynamic environment.Then,the forecasting analysis platform of end-of-life largescale spacecraft flying track i展开更多
The influences of fractal pore structure in coal reservoir on coalbed methane(CBM) migration were analyzed in detail by coupling theoretical models and numerical methods.Different types of fractals were generated base...The influences of fractal pore structure in coal reservoir on coalbed methane(CBM) migration were analyzed in detail by coupling theoretical models and numerical methods.Different types of fractals were generated based on the construction thought of the standard Menger Sponge to model the 3D nonlinear coal pore structures.Then a correlation model between the permeability of fractal porous medium and its pore-size-distribution characteristics was derived using the parallel and serial modes and verified by Lattice Boltzmann Method(LBM).Based on the coupled method,porosity(ф),fractal dimension of pore structure(Db),pore size range(rmin,rmax) and other parameters were systematically analyzed for their influences on the permeability(ф) of fractal porous medium.The results indicate that:① the channels connected by pores with the maximum size(rmax) dominate the permeability,approximating in the quadratic law;② the greater the ratio of r max and r min is,the higher is;③ the relationship between D b and follows a negative power law model,and breaks into two segments at the position where Db ≌2.5.Based on the results above,a predicting model of fractal porous medium permeability was proposed,formulated as k=cfrnmax,where C and n(approximately equal to 2) are constants and f is an expression only containing parameters of fractal pore structure.In addition,the equivalence of the new proposed model for porous medium and the Kozeny-Carman model k=Crn was verified at Db =2.0.展开更多
The studies of bionics reveal that some aquatic animals and winged insects have developed an unsmoothed surface possessing good characteristics of drag reduction. In this paper, four types of bionic surfaces, placoid-...The studies of bionics reveal that some aquatic animals and winged insects have developed an unsmoothed surface possessing good characteristics of drag reduction. In this paper, four types of bionic surfaces, placoid-shaped, V-shaped, rib- let-shaped, and ridge-shaped grooved surfaces, are employed as the microchannel surfaces for the purpose of reducing pressure loss. Lattice Boltzmann Method （LBM）, a new numerical approach on mescoscopic level, is used to conduct the numerical investigationsr. The results show that the micro-grooved surfaces possess the drag reduction performance. The existence of the vortices formed within the grooves not only decrease the shear force between fluid and wall but also minimize the contact area between fluid and walls, which can lead to a reduction of pressure loss. The drag reduction coefficient （η） for these four types of micro-structures could be generalized as follows： ηridge_shaped 〉 ηv-shaped 〉 ηplacoid-shaped 〉 ηriblet-shaped. Besides, the geometrical optimizations for the ridge-shaped grooves, which have the highest drag reduction performance, are performed as well. The results suggest that, for the purpose of drag reduction, the ridge-shaped grooves with smaller width to height ratio are recommended for the lower Reynolds number flow, while the ridge-shaped grooves with larger width to height ratio are be more suitable for the larger Reynolds number flow.展开更多
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.展开更多
A 2-D Lattice Boltzmann Method（LBM） coupled with a Sub-Grid Stress（SGS） model is proposed and validated by flows around a non-submerged spur dike in a channel.And then the LBM-SGS model is further applied to flows...A 2-D Lattice Boltzmann Method（LBM） coupled with a Sub-Grid Stress（SGS） model is proposed and validated by flows around a non-submerged spur dike in a channel.And then the LBM-SGS model is further applied to flows in a pumping-station forebay.Shallow water equations are numerically solved by the LBM and the turbulence can be taken into account and modeled efficiently by the Large Eddy Simulation（LES） model.The bounce-back scheme of the non-equilibrium part of the distribution function is used at the inlet boundary,the normal gradient of the distribution function is set as zero at the outlet boundary and the bounce-back scheme is applied to the solid wall to ensure non-slip boundary conditions.Firstly,the model successfully predicts the flow characteristics around a spur dike,such as circulating flow,velocity and water depth distributions.The results are verified by the experimental data and compared to the results obtained by conventional Smagoringsky Model（SM） of LES.Finally,the LBM-SGS model is used to further predict the flow characteristics in a forebay,such as secondary flow and water level.The comparisons show that the model scheme has the capacity to simulate complex flows in shallow water with reasonable accuracy and reliability.展开更多
Fuel cells and flow batteries are promising technologies to address climate change and air pollution problems. An understanding of the complex multiscale and multiphysics transport phenomena occurring in these electro...Fuel cells and flow batteries are promising technologies to address climate change and air pollution problems. An understanding of the complex multiscale and multiphysics transport phenomena occurring in these electrochemical systems requires powerful numerical tools. Over the past decades, the lattice Boltzmann (LB) method has attracted broad interest in the computational fluid dynamics and the numerical heat transfer communities, primarily due to its kinetic nature making it appropriate for modeling complex multiphase transport phenomena. More importantly, the LB method fits well with parallel computing due to its locality feature, which is required for large-scale engineering applications. In this article, we review the LB method for gas-liquid two-phase flows, coupled fluid flow and mass transport in porous media, and particulate flows. Examples of applications are provided in fuel cells and flow batteries. Further developments of the LB method are also outlined.展开更多
A single-phase free surface tracking model based on the Lattice Boltzmann Method （LBM）, which has capability of simulating liquid-gas system with the assumption that the gas phase has only negligible influence on th...A single-phase free surface tracking model based on the Lattice Boltzmann Method （LBM）, which has capability of simulating liquid-gas system with the assumption that the gas phase has only negligible influence on the liquid phase, is utilized to simulate the flow of a drop impacting on a liquid film. Three typical outcomes in the flows, i.e., deposition, crown and splash, which have been observed in the previous experiments, are obtained in the present three dimensional numerical simulations. The numerical results are consistent with the experimental and analytical results available.展开更多
This paper investigates the convergence proof of the Direct Simulation Monte Carlo(DSMC) method and the Gas-Kinetic Unified Algorithm in simulating the Boltzmann equation.It can be shown that the particle velocity dis...This paper investigates the convergence proof of the Direct Simulation Monte Carlo(DSMC) method and the Gas-Kinetic Unified Algorithm in simulating the Boltzmann equation.It can be shown that the particle velocity distribution function obtained by the DSMC method converges to a modified form of the Boltzmann equation,which is the equation of the gas-kinetic unified algorithm to directly solve the molecular velocity distribution function.Their convergence is derived through mathematical treatment.The collision frequency is presented using various molecular models and the local equilibrium distribution function is obtained by Enskog expansion using the converged equation of the DSMC method.These two expressions agree with those used in the unified algorithm.Numerical validation of the converging consistency between these two approaches is illustrated by simulating the pressure driven Poiseuille flow in the slip transition flow regime and the two-dimensional and three-dimensional flows around a circular cylinder and spherical-cone reentry body covering the whole flow regimes from low speed micro-channel flow to high speed non-equilibrium aerothermodynamics.展开更多
Two-phase flow in two digital cores is simulated by the color-gradient lattice Boltzmann method.This model can be applied totwo-phase flow with high-density ratio(on order of 1000).The first digital core is an artific...Two-phase flow in two digital cores is simulated by the color-gradient lattice Boltzmann method.This model can be applied totwo-phase flow with high-density ratio(on order of 1000).The first digital core is an artificial sandstone core,and itsthree-dimensional gray model is obtained by Micro-CT scanning.The gray scale images are segmented into discrete phases(solid particles and pore space) by the Otsu algorithm.The second one is a digital core of shale,which is reconstructed usingMarkov Chain Monte Carlo method with segmented SEM scanning image as input.The wettability of solid wall and relativepermeability of a cylindrical tube are simulated to verify the model.In the simulations of liquid and gas two phase flow in digital cores,density ratios of 100,200,500 and 1000 between liquid and gas are chosen.Based on the gas distribution in the digital core at different times,it is found that the fingering phenomenon is more salient at high density ratio.With the density ratioincreasing,the displacement efficiency decreases.Besides,due to numerous small pores in the shale,the displacement efficiency is over 20% less than that in the artificial sandstone and the difference is even about 30% when density ratio is greaterthan 500.As the density ratio increases,the gas saturation decreases in big pores,and even reaches zero in some small pores orbig pores with small throats.Residual liquid mainly distributes in the small pores and the edge of big pores due to the wettability of liquid.Liquid recovery can be enhanced effectively by decreasing its viscosity.展开更多
The dynamics of two-phase flows with a constant driving force inside a micro-channel is studied by using the Lattice Boltzmann Method （LBM） and the Shan-Chen multiphase model in this article. Flow regimes under diff...The dynamics of two-phase flows with a constant driving force inside a micro-channel is studied by using the Lattice Boltzmann Method （LBM） and the Shan-Chen multiphase model in this article. Flow regimes under different wall wettabilities and over smooth and grooved geometric surfaces are investigated. It is found that flow behaviors are strongly affected by the wall wettability and topography. Our results show that the LBM is efficient and accurate, and has very good application prospect in the study of drag reduction of microscopic seepage of reservoir.展开更多
The apparent slip between solid wall and liquid is studied by using the Lattice Boltzmann Method （LBM） and the Shan-Chen multiphase model in this paper. With a no-slip bounce-back scheme applied to the interface, fl...The apparent slip between solid wall and liquid is studied by using the Lattice Boltzmann Method （LBM） and the Shan-Chen multiphase model in this paper. With a no-slip bounce-back scheme applied to the interface, flow regimes under different wall wettabilities are investigated. Because of the wall wettability, liquid apparent slip is observed. Slip lengths for different wall wettabilities are found to collapse nearly onto a single curve as a function of the static contact angle, and thereby a relationship between apparent slip length and contact angle is suggested. Our results also show that the wall wettability leads to the formation of a low-density layer between solid wall and liquid, which produced apparent slip in the micro-scale.展开更多
Concentration gradient induced Rayleigh convection can influence effectively interracial mass transfer processes, but the convection phenomena are known as mesoscopic and complex. In order to investigate this phenomen...Concentration gradient induced Rayleigh convection can influence effectively interracial mass transfer processes, but the convection phenomena are known as mesoscopic and complex. In order to investigate this phenomenon, a two-equation Lattice Boltzmann Method （LBM） is proposed to simulate the velocity and the concentra-tion distributions of Rayleigh convection generated in the CO2 absorptlon into ethanol liquid.The simulated results on velocity distributions are experimentally verified by PIV （particle image velocimetry technique） measurements. In order to simplify the analysis, the convection in the simulation as well as in the experiment, the Rayleigh convection was manipulated into a single down flow pattern, The simulated results show that the concentration contours agree qualitatively with the schlieren images in the literature. The experimental and simulated results show that theRayleigh convection under investigation is dominated by the flow in the downward direction and impels exchange of the liquid between the interfacial vicinity and the liquid bulk promoting the renewal of interfacial liquid, and hence enhances mass transfer. The comparison between the simulated and experimental results demonstrated that the proposed LBM is a promising alternative for simulating mass transfer induced Rayleigh convection.展开更多
The two-dimensional Rayleigh-Taylor instability problem is simulated with a multiple-relaxation-time discrete Boltzmann model with a gravity term. Viscosity, heat conductivity, and Prandtl number effects are probed fr...The two-dimensional Rayleigh-Taylor instability problem is simulated with a multiple-relaxation-time discrete Boltzmann model with a gravity term. Viscosity, heat conductivity, and Prandtl number effects are probed from macroscopic and nonequilibrium viewpoints. In the macro sense, both viscosity and heat conduction show a significant inhibitory effect in the reacceleration stage, which is mainly achieved by inhibiting the development of the Kelvin-Helmholtz instability. Before this, the Prandtl number effect is not sensitive. Viscosity, heat conductivity, and Prandtl number effects on nonequilibrium manifestations and the degree of correlation between the nonuniformity and the nonequilibrium strength in the complex flow are systematically investigated.展开更多
Map is one of the communication means created by human being.Cartographers have been making efforts on the comparison of maps to natural languages so as to establish a"cartographic language"or"map langu...Map is one of the communication means created by human being.Cartographers have been making efforts on the comparison of maps to natural languages so as to establish a"cartographic language"or"map language".One of such efforts is to adopt the Shannon’s Information Theory originated in digital communication into cartography so as to establish an entropy-based cartographic communication theory.However,success has been very limited although research work had started as early as the mid-1960 s.It is then found that the bottleneck problem was the lack of appropriate measures for the spatial(configurational)information of(graphic and image)maps,as the classic Shannon entropy is only capable of characterizing statistical information but fails to capture the configurational information of(graphic and image)maps.Fortunately,after over 40-year development,some bottleneck problems have been solved.More precisely,generalized Shannon entropies for metric and thematic information of(graphic)maps have been developed and the first feasible solution for computing the Boltzmann entropy of image maps has been invented,which is capable of measuring the spatial information of not only numerical images but also categorical maps.With such progress,it is now feasible to build the"Information Theory of Cartography".In this paper,a framework for such a theory is proposed and some key issues are identified.For these issues,some have already been tackled while others still need efforts.As a result,a research agenda is set for future action.After all these issues are tackled,the theory will become matured so as to become a theoretic basis of cartography.It is expected that the Information Theory of Cartography will play an increasingly important role in the discipline of cartography because more and more researchers have advocated that information is more fundamental than matter and energy.展开更多
The lattice Boltzmann method （LBM） is used to examine free convection of nanofluids. The space between the cold outer square and heated inner circular cylinders is filled with water including various kinds of nanopa...The lattice Boltzmann method （LBM） is used to examine free convection of nanofluids. The space between the cold outer square and heated inner circular cylinders is filled with water including various kinds of nanoparticles： TiO2, Ag, Cu, and A1203. The Brinkman and Maxwell-Garnetts models are used to simulate the viscosity and the effective thermal conductivity of nanofluids, respectively. Results from the performed numerical analysis show good agreement with those obtained from other numerical meth- ods. A variety of the Rayleigh number, the nanoparticle volume fraction, and the aspect ratio are examined. According to the results, choosing copper as the nanoparticle leads to obtaining the highest enhancement for this problem. The results also indicate that the maximum value of enhancement occurs at λ =2.5 when Ra = 106 while at A = 1.5 for other Rayleigh numbers.展开更多
文摘In this mini-review we summarize the progress of Lattice Boltzmann （LB） modeling and simulating compressible flows in our group in recent years. Main contents include （i） Single-Relaxation-Time （SRT） LB model supplemented by additional viscosity, （ii） Multiple-Relaxation-Time （MRT） LB model, and （iii） LB study on hydrodynamic instabilities. The former two belong to improvements of physical modeling and the third belongs to simulation or application. The SRT-LB model sup- plemented by additional viscosity keeps the original framework of Lattice Bhatnagar-Gross Krook （LBGK）. So, it is easier and more convenient for previous SRT-LB users. The MRT-LB is a com- pletely new framework for physical modeling. It significantly extends the range of LB applications. The cost is longer computational time. The developed SRT-LB and MRT-LB are complementary from the sides of convenience and applicability.
基金support of this work (Grant EP/D031257/1)and the contribution of illustrations from Prof. U. Tuzun (University of Surrey), Prof. C. Lin(University of Utah) and Dr C. Selomulya (Monash University)
文摘The use of X-ray microtomographic （XMT） methods in analysing particulate systems has expanded rapidly in recent years with the availability of affordable desk-top apparatus. This review presents a summary of the major applications in which computer simulations are explicitly coupled with XMT in the area of granular and porous materials. We envisage two main ways of establishing the coupling between both techniques, based on the transference or exchange of information by using physical or geometrical parameters （i.e. a parametric link through fitting to a process model） or through the direct use of3D XMT digital images （i.e. comparing image pixels and features directly）. Examples of coupled applications are shown for the study of transport properties of rocks, particle packing, mechanical loading and sintering. Often, the link between XMT and computer simulations is based on visual comparisons and we conclude that the use of quantitative parameters such as the number of interparticle contacts, force networks or granule shape to link both techniques is still underrepresented in the literature. Strategies to provide a more robust and quantitative approach to optimise the information obtained from such tomography analyses are proposed. 2010 Chinese Society of Particuology and Institute of Process Engineering, Chinese Academy of Sciences. Published by Elsevier B.V. All rights reserved.
基金the National Natural Science Foundation of China (Grant No. 50579012)
文摘This article deals with the solute transport in a single fracture with the combination of the Lattice Boltzmann Method （LBM） and Modified Moment Propagation （MMP） method, and this mixed method is proved to have several advantages over the LBM and Moment Propagation （MP） mixed method which leads to negative concentrations under some conditions in computation. The disadvantage of LBM/MP has been overcome to a certain extent. Also, this work presents an LBM solution of modeling single fractures with uniformly or randomly distributed grains, which can provide a new path of applying the LBM in solute transport simulation in fractures.
文摘A supercomputer with 1.0 Petaflops peak performance in single precision, designed and established by Institute of Process Engineering, Chinese Academy of Sciences, is introduced in this brief communication. A designing philosophy utilizing the similarity between hardware, software and the problems to be solved is embodied, based on the multi-scale method and discrete simulation approaches developed at Institute of Process Engineering （IPE） and implemented in a graphic processing unit （GPU）-based hybrid computing mode. The preliminary applications of this machine in areas of multi-phase flow, molecular dynamics and so on are reported, demonstrating the supercomputer as a paradigm of green computation in new architecture.
基金The National Key Basic Research and Development Program(2014CB744100)and the National Natural Science Foundation of China(91530319 and 11325212)support the present researches in the design of the study and collection,analysis,and interpretation of data and in writing the manuscript.
文摘How to solve the hypersonic aerothermodynamics around large-scale uncontrolled spacecraft during falling disintegrated process from outer space to earth,is the key to resolve the problems of the uncontrolled Tiangong-No.1 spacecraft reentry crash.To study aerodynamics of spacecraft reentry covering various flow regimes,a Gas-Kinetic Unified Algorithm(GKUA)has been presented by computable modeling of the collision integral of the Boltzmann equation over tens of years.On this basis,the rotational and vibrational energy modes are considered as the independent variables of the gas molecular velocity distribution function,a kind of Boltzmann model equation involving in internal energy excitation is presented by decomposing the collision term of the Boltzmann equation into elastic and inelastic collision terms.Then,the gas-kinetic numerical scheme is constructed to capture the time evolution of the discretized velocity distribution functions by developing the discrete velocity ordinate method and numerical quadrature technique.The unified algorithm of the Boltzmann model equation involving thermodynamics non-equilibrium effect is presented for the whole range of flow regimes.The gas-kinetic massive parallel computing strategy is developed to solve the hypersonic aerothermodynamics with the processor cores 500~45,000 at least 80%parallel efficiency.To validate the accuracy of the GKUA,the hypersonic flows are simulated including the reentry Tiangong-1 spacecraft shape with the wide range of Knudsen numbers of 220~0.00005 by the comparison of the related results from the DSMC and N-S coupled methods,and the low-density tunnel experiment etc.For uncontrolling spacecraft falling problem,the finite-element algorithm for dynamic thermalforce coupling response is presented,and the unified simulation of the thermal structural response and the hypersonic flow field is tested on the Tiangong-1 shape under reentry aerodynamic environment.Then,the forecasting analysis platform of end-of-life largescale spacecraft flying track i
基金supported by National Natural Science Foundation of China(Grant Nos.41102093&41072153)CBM Union Foundation of Shanxi Province (Grant No.2012012002)Doctoral Scientific Foundation of Henan Polytechnic University(Grant No.648706)
文摘The influences of fractal pore structure in coal reservoir on coalbed methane(CBM) migration were analyzed in detail by coupling theoretical models and numerical methods.Different types of fractals were generated based on the construction thought of the standard Menger Sponge to model the 3D nonlinear coal pore structures.Then a correlation model between the permeability of fractal porous medium and its pore-size-distribution characteristics was derived using the parallel and serial modes and verified by Lattice Boltzmann Method(LBM).Based on the coupled method,porosity(ф),fractal dimension of pore structure(Db),pore size range(rmin,rmax) and other parameters were systematically analyzed for their influences on the permeability(ф) of fractal porous medium.The results indicate that:① the channels connected by pores with the maximum size(rmax) dominate the permeability,approximating in the quadratic law;② the greater the ratio of r max and r min is,the higher is;③ the relationship between D b and follows a negative power law model,and breaks into two segments at the position where Db ≌2.5.Based on the results above,a predicting model of fractal porous medium permeability was proposed,formulated as k=cfrnmax,where C and n(approximately equal to 2) are constants and f is an expression only containing parameters of fractal pore structure.In addition,the equivalence of the new proposed model for porous medium and the Kozeny-Carman model k=Crn was verified at Db =2.0.
文摘The studies of bionics reveal that some aquatic animals and winged insects have developed an unsmoothed surface possessing good characteristics of drag reduction. In this paper, four types of bionic surfaces, placoid-shaped, V-shaped, rib- let-shaped, and ridge-shaped grooved surfaces, are employed as the microchannel surfaces for the purpose of reducing pressure loss. Lattice Boltzmann Method （LBM）, a new numerical approach on mescoscopic level, is used to conduct the numerical investigationsr. The results show that the micro-grooved surfaces possess the drag reduction performance. The existence of the vortices formed within the grooves not only decrease the shear force between fluid and wall but also minimize the contact area between fluid and walls, which can lead to a reduction of pressure loss. The drag reduction coefficient （η） for these four types of micro-structures could be generalized as follows： ηridge_shaped 〉 ηv-shaped 〉 ηplacoid-shaped 〉 ηriblet-shaped. Besides, the geometrical optimizations for the ridge-shaped grooves, which have the highest drag reduction performance, are performed as well. The results suggest that, for the purpose of drag reduction, the ridge-shaped grooves with smaller width to height ratio are recommended for the lower Reynolds number flow, while the ridge-shaped grooves with larger width to height ratio are be more suitable for the larger Reynolds number flow.
基金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. 50779069)the Beijing Natural Science Foundation (Grant No. 3083022)+3 种基金the Open Fund of Key Laboratory of Yellow River Sediment Research of the Ministry of Water Resources (Grant No. 200903)the National Science and Technology Supporting Programs in the 11th Five-Year Plan (Grant No. 2006BAB06B02)the Chinese Universities Scientific Fund (Grant Nos. 2009-1-90, 2009-2-12)the Program for Changjiang Scholars and Innovative Research Team in University (Grant No. PCSIRT0657)
文摘A 2-D Lattice Boltzmann Method（LBM） coupled with a Sub-Grid Stress（SGS） model is proposed and validated by flows around a non-submerged spur dike in a channel.And then the LBM-SGS model is further applied to flows in a pumping-station forebay.Shallow water equations are numerically solved by the LBM and the turbulence can be taken into account and modeled efficiently by the Large Eddy Simulation（LES） model.The bounce-back scheme of the non-equilibrium part of the distribution function is used at the inlet boundary,the normal gradient of the distribution function is set as zero at the outlet boundary and the bounce-back scheme is applied to the solid wall to ensure non-slip boundary conditions.Firstly,the model successfully predicts the flow characteristics around a spur dike,such as circulating flow,velocity and water depth distributions.The results are verified by the experimental data and compared to the results obtained by conventional Smagoringsky Model（SM） of LES.Finally,the LBM-SGS model is used to further predict the flow characteristics in a forebay,such as secondary flow and water level.The comparisons show that the model scheme has the capacity to simulate complex flows in shallow water with reasonable accuracy and reliability.
基金supported by the Research Grants Council of the Hong Kong Special Administrative Region, China (Grant 623313)
文摘Fuel cells and flow batteries are promising technologies to address climate change and air pollution problems. An understanding of the complex multiscale and multiphysics transport phenomena occurring in these electrochemical systems requires powerful numerical tools. Over the past decades, the lattice Boltzmann (LB) method has attracted broad interest in the computational fluid dynamics and the numerical heat transfer communities, primarily due to its kinetic nature making it appropriate for modeling complex multiphase transport phenomena. More importantly, the LB method fits well with parallel computing due to its locality feature, which is required for large-scale engineering applications. In this article, we review the LB method for gas-liquid two-phase flows, coupled fluid flow and mass transport in porous media, and particulate flows. Examples of applications are provided in fuel cells and flow batteries. Further developments of the LB method are also outlined.
基金the National Natural Science Foundation of China (Grant No. 10472062)the Special Scientific Foundation for Selection and Cultivation of Excellent Young Scholars in Shanghai and the Shanghai Leading Academic Discipline Project (Grant No.Y0103)
文摘A single-phase free surface tracking model based on the Lattice Boltzmann Method （LBM）, which has capability of simulating liquid-gas system with the assumption that the gas phase has only negligible influence on the liquid phase, is utilized to simulate the flow of a drop impacting on a liquid film. Three typical outcomes in the flows, i.e., deposition, crown and splash, which have been observed in the previous experiments, are obtained in the present three dimensional numerical simulations. The numerical results are consistent with the experimental and analytical results available.
基金supported by the National Natural Science Foundation of China (Grant Nos. 91016027 and 91130018)
文摘This paper investigates the convergence proof of the Direct Simulation Monte Carlo(DSMC) method and the Gas-Kinetic Unified Algorithm in simulating the Boltzmann equation.It can be shown that the particle velocity distribution function obtained by the DSMC method converges to a modified form of the Boltzmann equation,which is the equation of the gas-kinetic unified algorithm to directly solve the molecular velocity distribution function.Their convergence is derived through mathematical treatment.The collision frequency is presented using various molecular models and the local equilibrium distribution function is obtained by Enskog expansion using the converged equation of the DSMC method.These two expressions agree with those used in the unified algorithm.Numerical validation of the converging consistency between these two approaches is illustrated by simulating the pressure driven Poiseuille flow in the slip transition flow regime and the two-dimensional and three-dimensional flows around a circular cylinder and spherical-cone reentry body covering the whole flow regimes from low speed micro-channel flow to high speed non-equilibrium aerothermodynamics.
基金supported by the National Natural Science Foundation of China(Grant No.51234007,51404291)Program for Changjiang Scholars and Innovative Research Team in University(Grant No.IRT1294)Introducing Talents of Discipline to Universities(Grant No.B08028)
文摘Two-phase flow in two digital cores is simulated by the color-gradient lattice Boltzmann method.This model can be applied totwo-phase flow with high-density ratio(on order of 1000).The first digital core is an artificial sandstone core,and itsthree-dimensional gray model is obtained by Micro-CT scanning.The gray scale images are segmented into discrete phases(solid particles and pore space) by the Otsu algorithm.The second one is a digital core of shale,which is reconstructed usingMarkov Chain Monte Carlo method with segmented SEM scanning image as input.The wettability of solid wall and relativepermeability of a cylindrical tube are simulated to verify the model.In the simulations of liquid and gas two phase flow in digital cores,density ratios of 100,200,500 and 1000 between liquid and gas are chosen.Based on the gas distribution in the digital core at different times,it is found that the fingering phenomenon is more salient at high density ratio.With the density ratioincreasing,the displacement efficiency decreases.Besides,due to numerous small pores in the shale,the displacement efficiency is over 20% less than that in the artificial sandstone and the difference is even about 30% when density ratio is greaterthan 500.As the density ratio increases,the gas saturation decreases in big pores,and even reaches zero in some small pores orbig pores with small throats.Residual liquid mainly distributes in the small pores and the edge of big pores due to the wettability of liquid.Liquid recovery can be enhanced effectively by decreasing its viscosity.
基金supported by the National Natural Science Foundation of China(Grant No.50874071)the National High Technology Research and Development Program of China(863 Program,Grant No.2008AA06Z201)+3 种基金the Program of Science and Technology Commission of Shanghai Municipality(Grant No.071605102)the project of Shanghai Education Commission Research Innovation(Grand No.08ZZ45)the Program for Changjiang Scholars and Innovative Research Team in Universities(Grant No.IRT0844)the China Postdoctoral Science Foundation(Grant No.20090450687)
文摘The dynamics of two-phase flows with a constant driving force inside a micro-channel is studied by using the Lattice Boltzmann Method （LBM） and the Shan-Chen multiphase model in this article. Flow regimes under different wall wettabilities and over smooth and grooved geometric surfaces are investigated. It is found that flow behaviors are strongly affected by the wall wettability and topography. Our results show that the LBM is efficient and accurate, and has very good application prospect in the study of drag reduction of microscopic seepage of reservoir.
基金the National Natural Science Foundation of China (Grant No. 50874071)the National High Technology Research and Development of China (863 Program,Grant No. 2008AA06Z201)+1 种基金the Key Program of Science and Technology Commission of Shanghai Municipality (Grant No.071605102)the Leading Talent Funding of Shanghai
文摘The apparent slip between solid wall and liquid is studied by using the Lattice Boltzmann Method （LBM） and the Shan-Chen multiphase model in this paper. With a no-slip bounce-back scheme applied to the interface, flow regimes under different wall wettabilities are investigated. Because of the wall wettability, liquid apparent slip is observed. Slip lengths for different wall wettabilities are found to collapse nearly onto a single curve as a function of the static contact angle, and thereby a relationship between apparent slip length and contact angle is suggested. Our results also show that the wall wettability leads to the formation of a low-density layer between solid wall and liquid, which produced apparent slip in the micro-scale.
基金Supported by the National Natural Science Foundation of China （20736005）.
文摘Concentration gradient induced Rayleigh convection can influence effectively interracial mass transfer processes, but the convection phenomena are known as mesoscopic and complex. In order to investigate this phenomenon, a two-equation Lattice Boltzmann Method （LBM） is proposed to simulate the velocity and the concentra-tion distributions of Rayleigh convection generated in the CO2 absorptlon into ethanol liquid.The simulated results on velocity distributions are experimentally verified by PIV （particle image velocimetry technique） measurements. In order to simplify the analysis, the convection in the simulation as well as in the experiment, the Rayleigh convection was manipulated into a single down flow pattern, The simulated results show that the concentration contours agree qualitatively with the schlieren images in the literature. The experimental and simulated results show that theRayleigh convection under investigation is dominated by the flow in the downward direction and impels exchange of the liquid between the interfacial vicinity and the liquid bulk promoting the renewal of interfacial liquid, and hence enhances mass transfer. The comparison between the simulated and experimental results demonstrated that the proposed LBM is a promising alternative for simulating mass transfer induced Rayleigh convection.
文摘The two-dimensional Rayleigh-Taylor instability problem is simulated with a multiple-relaxation-time discrete Boltzmann model with a gravity term. Viscosity, heat conductivity, and Prandtl number effects are probed from macroscopic and nonequilibrium viewpoints. In the macro sense, both viscosity and heat conduction show a significant inhibitory effect in the reacceleration stage, which is mainly achieved by inhibiting the development of the Kelvin-Helmholtz instability. Before this, the Prandtl number effect is not sensitive. Viscosity, heat conductivity, and Prandtl number effects on nonequilibrium manifestations and the degree of correlation between the nonuniformity and the nonequilibrium strength in the complex flow are systematically investigated.
基金National Natural Science Foundation of China(Nos.41930104,41971330)Hong Kong Research Grants Council General Research Fund(No.152219/18E)。
文摘Map is one of the communication means created by human being.Cartographers have been making efforts on the comparison of maps to natural languages so as to establish a"cartographic language"or"map language".One of such efforts is to adopt the Shannon’s Information Theory originated in digital communication into cartography so as to establish an entropy-based cartographic communication theory.However,success has been very limited although research work had started as early as the mid-1960 s.It is then found that the bottleneck problem was the lack of appropriate measures for the spatial(configurational)information of(graphic and image)maps,as the classic Shannon entropy is only capable of characterizing statistical information but fails to capture the configurational information of(graphic and image)maps.Fortunately,after over 40-year development,some bottleneck problems have been solved.More precisely,generalized Shannon entropies for metric and thematic information of(graphic)maps have been developed and the first feasible solution for computing the Boltzmann entropy of image maps has been invented,which is capable of measuring the spatial information of not only numerical images but also categorical maps.With such progress,it is now feasible to build the"Information Theory of Cartography".In this paper,a framework for such a theory is proposed and some key issues are identified.For these issues,some have already been tackled while others still need efforts.As a result,a research agenda is set for future action.After all these issues are tackled,the theory will become matured so as to become a theoretic basis of cartography.It is expected that the Information Theory of Cartography will play an increasingly important role in the discipline of cartography because more and more researchers have advocated that information is more fundamental than matter and energy.
文摘The lattice Boltzmann method （LBM） is used to examine free convection of nanofluids. The space between the cold outer square and heated inner circular cylinders is filled with water including various kinds of nanoparticles： TiO2, Ag, Cu, and A1203. The Brinkman and Maxwell-Garnetts models are used to simulate the viscosity and the effective thermal conductivity of nanofluids, respectively. Results from the performed numerical analysis show good agreement with those obtained from other numerical meth- ods. A variety of the Rayleigh number, the nanoparticle volume fraction, and the aspect ratio are examined. According to the results, choosing copper as the nanoparticle leads to obtaining the highest enhancement for this problem. The results also indicate that the maximum value of enhancement occurs at λ =2.5 when Ra = 106 while at A = 1.5 for other Rayleigh numbers.
