The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function (DDF) lattice Boltzmann method (LBM) for the compressible Navier-Stokes equations. Firstly, ...The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function (DDF) lattice Boltzmann method (LBM) for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma:0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.展开更多
The main goal of this paper is to investigate natural convective heat transfer and flow characteristics of non-Newtonian nanofluid streaming between two infinite vertical flat plates in the presence of magnetic field ...The main goal of this paper is to investigate natural convective heat transfer and flow characteristics of non-Newtonian nanofluid streaming between two infinite vertical flat plates in the presence of magnetic field and thermal radiation.Initially,a similarity transformation is used to convert momentum and energy conservation equations in partial differential forms into non-linear ordinary differential equations (ODE) applying meaningful boundary conditions.In order to obtain the non-linear ODEs analytically,Galerkin method (GM) is employed.Subsequently,the ODEs are also solved by a reliable numerical solution.In order to test the accuracy,precision and reliability of the analytical method,results of the analytical analysis are compared with the numerical results.With respect to the comparisons,fairly good compatibilities with insignificant errors are observed.Eventually,the impacts of effective parameters including magnetic and radiation parameters and nanofluid volume fraction on the velocity,skin friction coefficient and Nusselt number distributions are comprehensively described.Based on the results,it is revealed that with increasing the role of magnetic force,velocity profile,skin friction coefficient and thermal performance descend.Radiation parameter has insignificant influence on velocity profile while it obviously has augmentative and decreasing effects on skin friction and Nusselt number,respectively.展开更多
The discrete particle method was used to simulate the distribution of gas holdup in a gas-liquid standard Rushton stirred tank. The gas phase was treated as a large number of bubbles and their trajectories were tracke...The discrete particle method was used to simulate the distribution of gas holdup in a gas-liquid standard Rushton stirred tank. The gas phase was treated as a large number of bubbles and their trajectories were tracked with the results of motion equations. The two-way approach was performed to couple the interphase momentum exchange. The turbulent dispersion of bubbles with a size distribution was modeled using a stochastic tracking model, and the added mass force was involved to account for the effect of bubble acceleration on the surrounding fluid. The predicted gas holdup distribution showed that this method could give reasonable prediction comparable to the reported experimental data when the effect of turbulence was took into account in modification for drag coefficient.展开更多
The supercritical fluid crystallization technique is a novel technology for preparing ultrafine particles. This paper introduced the concept and features of the technique with an emphasis on three kinds of supercritic...The supercritical fluid crystallization technique is a novel technology for preparing ultrafine particles. This paper introduced the concept and features of the technique with an emphasis on three kinds of supercritical fluid crystallization techniques, i.e. rapid expansion of supercritical solutions, supercritical fluid anti-solvent and particles from gas saturated solutions Some questions and the prospect of this technique were also discussed.展开更多
A two dimensional mathematical model was developed to predict the performance characteristics for direct current, linear channel MHD propulsion system in a closed loop environment. The results of analytical and exper...A two dimensional mathematical model was developed to predict the performance characteristics for direct current, linear channel MHD propulsion system in a closed loop environment. The results of analytical and experimental studies of the linear channel MHD propulsor are described. Compared with the data of experiment, the correctness of the computation program is validated.展开更多
Convective pore-fluid flow (CPFF) plays a critical role in generating mineral deposits and oil reservoirs within the deep Earth. Therefore, theoretical understanding and numerical modeling of the thermodynamic process...Convective pore-fluid flow (CPFF) plays a critical role in generating mineral deposits and oil reservoirs within the deep Earth. Therefore, theoretical understanding and numerical modeling of the thermodynamic process that triggers and controls the CPFF are extremely important for the exploration of new mineral deposits and underground oil resources. From the viewpoint of science, the CPFF within the upper crust can be treated as a kind of thermodynamic instability problem of pore-fluid in fluid-saturated porous media. The key issue of dealing with this kind of problem is to assess whether a nonlinear thermodynamic system under consideration is supercritical. To overcome limitations of using theoretical analysis and experimental methods in dealing with the CPFF problems within the upper crust, finite element modeling has been broadly employed for solving this kind of problem over the past two decades. The main purpose of this paper is to overview recent developments and applications of finite element modeling associated with solving the CPFF problems in large length-scale geological systems of complicated geometries and complex material distributions. In particular, two kinds of commonly-used finite element modeling approaches, namely the steady-state and transient-state approaches, and their advantages/disadvantages are thoroughly presented and discussed.展开更多
A numerical method was developed to directly simulate the compressible, particle-laden turbulent jets.The fourth order compact finite difference schemes were used to discretize the space derivatives. The Lagrangian me...A numerical method was developed to directly simulate the compressible, particle-laden turbulent jets.The fourth order compact finite difference schemes were used to discretize the space derivatives. The Lagrangian method was adopted to simulate the particle motion based on one-way coupling. It is found that the turbulent intensity profiles attain self-similar status in the jet downstream regions. At the Stokes number of 1, particles are concentrated largely in the outer boundaries of the large-scale vortex structures with the most uneven distribution and the widest dispersion in the lateral direction. Particles at the much smaller Stokes numbers are distributed evenly in the flow field, and the lateral dispersion is also considerable. Distribution of particles at much larger Stokes numbers is more uniform and the lateral dispersion becomes small. In addition, the inflow conditions have different effects on the particle dispersion. The direct numerical simulation (DNS) results accord with the previous experiments and numerical studies.展开更多
The Xiangshan deposit in Jiangxi province is one of the most important uranium deposits in China. The aim of this study is to achieve a better understanding of mineralization in the Xiangshan deposit through numerical...The Xiangshan deposit in Jiangxi province is one of the most important uranium deposits in China. The aim of this study is to achieve a better understanding of mineralization in the Xiangshan deposit through numerical simulation. In order to find the most favorable locations of mineralization and to help further mineral exploration, a coupling deforma- tion and fluid flow model has been established to describe the mineralization process. In this model, the simulation re- constructs the strata deformations under fields of compressive stress and thrust structure on the hanging wall of the Zou-Shi fault. Compared with practical information, the simulation results are consistent with the No. 51 exploration section of the western Xiangshan. In addition, on the basis of geological information provided by previous investigators, the model simulates the flow process of fluids under compressive stress fields. The result suggests that many tensional areas are formed, which can help the fluid flowing upward from deeper parts. The fluid is easy to concentrate on the breccia fractured zone between two volcanic layers, especially on the intersection parts with faults, resulting in the for- mation of favourable locations of mineralization. In addition, the model is significant in guiding the exploration of ura- nium deposits in the western Xiangshan and provides clues for further exploration of deposits.展开更多
The steady, laminar, incompressible and two dimensional micropolar flow between two porous disks was investigated using optimal homotopy asymptotic method(OHAM) and fourth order Runge–Kutta numerical method. Comparis...The steady, laminar, incompressible and two dimensional micropolar flow between two porous disks was investigated using optimal homotopy asymptotic method(OHAM) and fourth order Runge–Kutta numerical method. Comparison between OHAM and numerical method shows that OHAM is an exact and high efficient method for solving these kinds of problems. The results are presented to study the velocity and rotation profiles for different physical parameters such as Reynolds number, vortex viscosity parameter, spin gradient viscosity and microinertia density parameter. As an important outcome, the magnitude of the microrotation increases with an increase in the values of injection velocity while it decreases by increasing the values of suction velocity.展开更多
This paper summarizes some of the typical riser vortex-induced vibration (VIV) problems in subsea oil and gas developments, and presents the corresponding computational fluid dynamics (CFD) time domain simula- tio...This paper summarizes some of the typical riser vortex-induced vibration (VIV) problems in subsea oil and gas developments, and presents the corresponding computational fluid dynamics (CFD) time domain simula- tion results to address these problems. First, the CFD time domain simulation approach was applied to analyze the wake field behind a stationary cylinder and a vibrating cylinder. Then a vertical riser VIV response under uniform current was studied. The VIV response time histories revealed some valuable clues that could lead to explanation of the higher harmonics. After that, a vertical riser VIV response under shear current was investigated. A 3 000 ft (1 ft=-0.304 8 m) water depth top tensioned riser was sized, and its VIV responses under uniform and shear current were studied. Then this paper continues to discuss one catenary flexible riser VIV response during normal lay. Last, the time domain simulation approach was applied to a partially submerged flexible jumper, to study the jumper VIV behavior, and dynamic motion envelopes. It was demonstrated that the time domain simulation ap- proach is able to disclose details of the flow field, vortex shedding pattern, and riser dynamic behavior, and han- dle different tvoes of risers under different Woe of currents.展开更多
According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in...According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.展开更多
The hydrodynamic forces and flow field of artificial reef models in steady flow were numerically investigated using the RNG κ-ε turbulent model. The numerical simulation results are consistent with results observed ...The hydrodynamic forces and flow field of artificial reef models in steady flow were numerically investigated using the RNG κ-ε turbulent model. The numerical simulation results are consistent with results observed by experimental means. A comparative study indicates that the corresponding errors of forces between calculated values and values observed in the experiment vary in the range of2.3%-11.2% and that the corresponding errors of velocities vary in the range of 1.3%-15.8%. The flow field numerical results show that upstream and vortices exist when the current passes over and through the surface of the reef model. This study suggests that the numerical simulation method can be applied to predict the forces and flow field associated with artificial reefs.展开更多
A three-dimensional, nineteen-velocity(D3Q19) Lattice Boltzmann Method(LBM) model was developed to simulate the fluid flow of a laminar square jet in cross flows based on the single relaxation time algorithm. The code...A three-dimensional, nineteen-velocity(D3Q19) Lattice Boltzmann Method(LBM) model was developed to simulate the fluid flow of a laminar square jet in cross flows based on the single relaxation time algorithm. The code was validated by the mathematic solution of the Poiseuille flow in a square channel, and was further validated with a previous well studied empirical correlation for the central trajectory of a jet in cross flows. The developed LBM model was found to be able to capture the dominant vortex, i.e. the Counter-rotating Vortex Pair(CVP) and the upright wake vortex. Results show that the incoming fluid in the cross flow channel was entrained into the leeside of the jet fluid, which contributes to the blending of the jet. That the spread width of the transverse jet decreases with the velocity ratio. A layer-organized entrainment pattern was found indicating that the incoming fluid at the lower position is firstly entrained into the leeside of the jet, and followed by the incoming fluid at the upper position.展开更多
基金Supported by the Aeronautical Science Foundation of China(20061453020)Foundation for Basic Research of Northwestern Polytechnical University(03)~~
文摘The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function (DDF) lattice Boltzmann method (LBM) for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma:0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.
文摘The main goal of this paper is to investigate natural convective heat transfer and flow characteristics of non-Newtonian nanofluid streaming between two infinite vertical flat plates in the presence of magnetic field and thermal radiation.Initially,a similarity transformation is used to convert momentum and energy conservation equations in partial differential forms into non-linear ordinary differential equations (ODE) applying meaningful boundary conditions.In order to obtain the non-linear ODEs analytically,Galerkin method (GM) is employed.Subsequently,the ODEs are also solved by a reliable numerical solution.In order to test the accuracy,precision and reliability of the analytical method,results of the analytical analysis are compared with the numerical results.With respect to the comparisons,fairly good compatibilities with insignificant errors are observed.Eventually,the impacts of effective parameters including magnetic and radiation parameters and nanofluid volume fraction on the velocity,skin friction coefficient and Nusselt number distributions are comprehensively described.Based on the results,it is revealed that with increasing the role of magnetic force,velocity profile,skin friction coefficient and thermal performance descend.Radiation parameter has insignificant influence on velocity profile while it obviously has augmentative and decreasing effects on skin friction and Nusselt number,respectively.
基金Supported by the National Natural Science Foundation of China (No.20776121) and the Specialized Research Fund for the Doctoral Program of Higher Education (No.20050530001), and the Scientific Research Fund of Hunan Provincial Education Department (No.07C765).
文摘The discrete particle method was used to simulate the distribution of gas holdup in a gas-liquid standard Rushton stirred tank. The gas phase was treated as a large number of bubbles and their trajectories were tracked with the results of motion equations. The two-way approach was performed to couple the interphase momentum exchange. The turbulent dispersion of bubbles with a size distribution was modeled using a stochastic tracking model, and the added mass force was involved to account for the effect of bubble acceleration on the surrounding fluid. The predicted gas holdup distribution showed that this method could give reasonable prediction comparable to the reported experimental data when the effect of turbulence was took into account in modification for drag coefficient.
文摘The supercritical fluid crystallization technique is a novel technology for preparing ultrafine particles. This paper introduced the concept and features of the technique with an emphasis on three kinds of supercritical fluid crystallization techniques, i.e. rapid expansion of supercritical solutions, supercritical fluid anti-solvent and particles from gas saturated solutions Some questions and the prospect of this technique were also discussed.
文摘A two dimensional mathematical model was developed to predict the performance characteristics for direct current, linear channel MHD propulsion system in a closed loop environment. The results of analytical and experimental studies of the linear channel MHD propulsor are described. Compared with the data of experiment, the correctness of the computation program is validated.
基金Project(11272359)supported by the National Natural Science Foundation of China
文摘Convective pore-fluid flow (CPFF) plays a critical role in generating mineral deposits and oil reservoirs within the deep Earth. Therefore, theoretical understanding and numerical modeling of the thermodynamic process that triggers and controls the CPFF are extremely important for the exploration of new mineral deposits and underground oil resources. From the viewpoint of science, the CPFF within the upper crust can be treated as a kind of thermodynamic instability problem of pore-fluid in fluid-saturated porous media. The key issue of dealing with this kind of problem is to assess whether a nonlinear thermodynamic system under consideration is supercritical. To overcome limitations of using theoretical analysis and experimental methods in dealing with the CPFF problems within the upper crust, finite element modeling has been broadly employed for solving this kind of problem over the past two decades. The main purpose of this paper is to overview recent developments and applications of finite element modeling associated with solving the CPFF problems in large length-scale geological systems of complicated geometries and complex material distributions. In particular, two kinds of commonly-used finite element modeling approaches, namely the steady-state and transient-state approaches, and their advantages/disadvantages are thoroughly presented and discussed.
基金Supported by the Natural Science Foundation of Zhejiang Province (No. 502047 and No. M503094)National Basic Research Program of China (No. 2003CB214500).
文摘A numerical method was developed to directly simulate the compressible, particle-laden turbulent jets.The fourth order compact finite difference schemes were used to discretize the space derivatives. The Lagrangian method was adopted to simulate the particle motion based on one-way coupling. It is found that the turbulent intensity profiles attain self-similar status in the jet downstream regions. At the Stokes number of 1, particles are concentrated largely in the outer boundaries of the large-scale vortex structures with the most uneven distribution and the widest dispersion in the lateral direction. Particles at the much smaller Stokes numbers are distributed evenly in the flow field, and the lateral dispersion is also considerable. Distribution of particles at much larger Stokes numbers is more uniform and the lateral dispersion becomes small. In addition, the inflow conditions have different effects on the particle dispersion. The direct numerical simulation (DNS) results accord with the previous experiments and numerical studies.
基金Projects GPMR0547 supported by the State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, and 2002 CB 412601 by the Chinese Ministry of Science and Technology
文摘The Xiangshan deposit in Jiangxi province is one of the most important uranium deposits in China. The aim of this study is to achieve a better understanding of mineralization in the Xiangshan deposit through numerical simulation. In order to find the most favorable locations of mineralization and to help further mineral exploration, a coupling deforma- tion and fluid flow model has been established to describe the mineralization process. In this model, the simulation re- constructs the strata deformations under fields of compressive stress and thrust structure on the hanging wall of the Zou-Shi fault. Compared with practical information, the simulation results are consistent with the No. 51 exploration section of the western Xiangshan. In addition, on the basis of geological information provided by previous investigators, the model simulates the flow process of fluids under compressive stress fields. The result suggests that many tensional areas are formed, which can help the fluid flowing upward from deeper parts. The fluid is easy to concentrate on the breccia fractured zone between two volcanic layers, especially on the intersection parts with faults, resulting in the for- mation of favourable locations of mineralization. In addition, the model is significant in guiding the exploration of ura- nium deposits in the western Xiangshan and provides clues for further exploration of deposits.
文摘The steady, laminar, incompressible and two dimensional micropolar flow between two porous disks was investigated using optimal homotopy asymptotic method(OHAM) and fourth order Runge–Kutta numerical method. Comparison between OHAM and numerical method shows that OHAM is an exact and high efficient method for solving these kinds of problems. The results are presented to study the velocity and rotation profiles for different physical parameters such as Reynolds number, vortex viscosity parameter, spin gradient viscosity and microinertia density parameter. As an important outcome, the magnitude of the microrotation increases with an increase in the values of injection velocity while it decreases by increasing the values of suction velocity.
文摘This paper summarizes some of the typical riser vortex-induced vibration (VIV) problems in subsea oil and gas developments, and presents the corresponding computational fluid dynamics (CFD) time domain simula- tion results to address these problems. First, the CFD time domain simulation approach was applied to analyze the wake field behind a stationary cylinder and a vibrating cylinder. Then a vertical riser VIV response under uniform current was studied. The VIV response time histories revealed some valuable clues that could lead to explanation of the higher harmonics. After that, a vertical riser VIV response under shear current was investigated. A 3 000 ft (1 ft=-0.304 8 m) water depth top tensioned riser was sized, and its VIV responses under uniform and shear current were studied. Then this paper continues to discuss one catenary flexible riser VIV response during normal lay. Last, the time domain simulation approach was applied to a partially submerged flexible jumper, to study the jumper VIV behavior, and dynamic motion envelopes. It was demonstrated that the time domain simulation ap- proach is able to disclose details of the flow field, vortex shedding pattern, and riser dynamic behavior, and han- dle different tvoes of risers under different Woe of currents.
基金Project(51074027)supported by the National Natural Science Foundation of China
文摘According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.
基金Supported by the National High Technology Research and Development Program of China(863 Programs)(No.2006AA100301)Science and Technology Development Program of Shandong Province(No.2005GG3205102)
文摘The hydrodynamic forces and flow field of artificial reef models in steady flow were numerically investigated using the RNG κ-ε turbulent model. The numerical simulation results are consistent with results observed by experimental means. A comparative study indicates that the corresponding errors of forces between calculated values and values observed in the experiment vary in the range of2.3%-11.2% and that the corresponding errors of velocities vary in the range of 1.3%-15.8%. The flow field numerical results show that upstream and vortices exist when the current passes over and through the surface of the reef model. This study suggests that the numerical simulation method can be applied to predict the forces and flow field associated with artificial reefs.
基金Supported by the National Natural Science Foundation of China(51476145,51476146)
文摘A three-dimensional, nineteen-velocity(D3Q19) Lattice Boltzmann Method(LBM) model was developed to simulate the fluid flow of a laminar square jet in cross flows based on the single relaxation time algorithm. The code was validated by the mathematic solution of the Poiseuille flow in a square channel, and was further validated with a previous well studied empirical correlation for the central trajectory of a jet in cross flows. The developed LBM model was found to be able to capture the dominant vortex, i.e. the Counter-rotating Vortex Pair(CVP) and the upright wake vortex. Results show that the incoming fluid in the cross flow channel was entrained into the leeside of the jet fluid, which contributes to the blending of the jet. That the spread width of the transverse jet decreases with the velocity ratio. A layer-organized entrainment pattern was found indicating that the incoming fluid at the lower position is firstly entrained into the leeside of the jet, and followed by the incoming fluid at the upper position.