Field emission electric propulsion(FEEP) thrusters possess excellent characteristics, such as high specific impulse, low power requirements, compact size and precise pointing capabilities,making them ideal propulsion ...Field emission electric propulsion(FEEP) thrusters possess excellent characteristics, such as high specific impulse, low power requirements, compact size and precise pointing capabilities,making them ideal propulsion devices for micro-nano satellites. However, the detection of certain aspects, such as the evolution process of the liquid cone and the physical quantities at the cone apex, proves challenging due to the minute size of the needle tip and the vacuum environment in which they operate. Consequently, this paper introduces a computational fluid dynamics(CFD) model to gain insight into the formation process of the liquid cone on the tip apex of indium FEEP. The CFD model is based on electrohydrodynamic(EHD) equations and the volume of fluid(VOF) method. The entire cone formation process can be divided into three stages, and the time-dependent characteristics of the physical quantities at the cone apex are investigated. The influences of film thickness, apex radius size and applied voltage are compared.The results indicate a gradual increase in the values of electrostatic stress and surface tension stress at the cone apex over an initial period, followed by a rapid escalation within a short duration.Apex configurations featuring a small radius, thick film and high voltage exhibit a propensity for liquid cone formation, and the cone growth time decreases as the film thickness increases.Moreover, some unstable behavior is observed during the cone formation process.展开更多
Dry storage containers must be secured and reliable during long-term storage,and the effect of decay heat released from the internal spent fuel on the cask has become an important research topic.In this paper,a 3D com...Dry storage containers must be secured and reliable during long-term storage,and the effect of decay heat released from the internal spent fuel on the cask has become an important research topic.In this paper,a 3D computational fluid dynamics model is presented,and the accuracy of the calculation is verified,with computational errors of less than 6.2%.The thermal stress of the dry storage cask was estimated by coupling it with a transient temperature field.The total power remained constant and adjusting the power ratio of the inner and outer zones had a small effect on the stress results,with a maximum equivalent stress of approximately 5.2 kPa,which occurred at the lower edge of the shell.In the case of tilt,the temperature gradient varied in a wavy distribution,and the wave crest moved from right to left.Altering the tilt angle affects the air distribution in the annular gap,leading to the shell temperature being transformed,with a maximum equivalent stress of 202 MPa at the bottom of the shell.However,the equivalent stress in both cases was less than the yield stress(205 MPa).展开更多
Computational Fluid Dynamics (CFD) has become an alternative method to experiments for understanding the fluid dynamics of multiphase flow. A two-fluid model, which contains additional terms in both the gas- and sol...Computational Fluid Dynamics (CFD) has become an alternative method to experiments for understanding the fluid dynamics of multiphase flow. A two-fluid model, which contains additional terms in both the gas- and solid-phase momentum equations, is used to investigate the fluidization quality in a fluidized bed. A case study for quartz sand with a density of 2,660 kg/m^3 and a diameter of 500 μm, whose physical property is similar to a new kind of catalyst for producing clean fuels through the residue fluid catalytic cracking process, is simulated in a two-dimensional fluidized bed with 0.57 m width and 1.00 m height. Transient bubbling and collapsing characteristics are numerically investigated in the platform of CFX 4.4 by integrating user-defined Fortran subroutines. The results show that the fluidization and collapse process is in fair agreement with the classical theory of Geldart B classification, but the collapse time is affected by bubbles at the interface between the dense phase and freeboard.展开更多
Turbidity currents represent a major agent for sediment transport in lakes, seas and oceans. In particu-lar, they formulate the most significant clastic accumulations in the deep sea, which become many of the world's...Turbidity currents represent a major agent for sediment transport in lakes, seas and oceans. In particu-lar, they formulate the most significant clastic accumulations in the deep sea, which become many of the world's most important hydrocarbon reservoirs. Several boreholes in the Qiongdongnan Basin, the north-western South China Sea, have recently revealed turbidity current deposits as significant hydrocarbon res-ervoirs. However, there are some arguments for the potential provenances. To solve this problem, it is es-sential to delineate their sedimentary processes as well as to evaluate their qualities as reservoir. Numerical simulations have been developed rapidly over the last several years, offering insights into turbidity current behaviors, as geologically significant turbidity currents are difficult to directly investigate due to their large scale and often destructive nature. Combined with the interpretation of the turbidity system based on high-resolution 3D seismic data, the paleotophography is acquired via a back-stripping seismic profile integrated with a borehole, i.e., Well A, in the western Qiongdongnan Basin; then a numerical model is built on the basis of this back-stripped profile. After defining the various turbidity current initial boundary conditions, includ-ing grain size, velocity and sediment concentration, the structures and behaviors of turbidity currents are investigated via numerical simulation software ANSYS FLUENT. Finally, the simulated turbidity deposits are compared with the interpreted sedimentary bodies based on 3D seismic data and the potential provenances of the revealed turbidites by Well A are discussed in details. The simulation results indicate that a sedimen-tary body develops far away from its source with an average grain size of 0.1 mm, i.e., sand-size sediment. Taking into account the location and orientation of the simulated seismic line, the consistence between normal forward simulation results and the revealed cores in Well A indicates that the turbidites should have been transported from Vietnam instead of Hainan Island. This interpretation has also been verified by the planar maps of sedimentary systems based on integration of boreholes and seismic data. The identification of the turbidity provenance will benefit the evaluation of extensively distributed submarine fans for hydro-carbon exploration in the deepwater areas.展开更多
Computational fluid dynamics(CFD) simulations are adopted to investigate rectangular microchannel flows with various periodic micro-structured wall by introducing velocity slip boundary condition at low Reynolds num...Computational fluid dynamics(CFD) simulations are adopted to investigate rectangular microchannel flows with various periodic micro-structured wall by introducing velocity slip boundary condition at low Reynolds number. The purpose of the current study is to numerically find out the effects of periodic micro-structured wall on the flow resistance in rectangular microchannel with the different spacings between microridges ranging from 15 to 60 pm. The simulative results indicate that pressure drop with different spacing between microridges increases linearly with flow velocity and decreases monotonically with slip velocity; Pressure drop reduction also increases with the spacing between microridges at the same condition of slip velocity and flow velocity. The results of numerical simulation are compared with theoretical predictions and experimental results in the literatures. It is found that there is qualitative agreement between them.展开更多
A simple hydrodynamic model based on two-fluid theory, taking into account the effect of discrete particles on both the gas- and solid-phase momentum equations, was used to numerically investigate the pressure fluctua...A simple hydrodynamic model based on two-fluid theory, taking into account the effect of discrete particles on both the gas- and solid-phase momentum equations, was used to numerically investigate the pressure fluctuation characteristics in a gas-solid fluidized bed with the aid of CFX 4.4, a commercial CFD software package, by adding user-defined Fortran subroutines. Numerical simulations together with typical experimental measurements show that pressure fluctuations originate above the distributor when a gas pulse is injected into the fluidized bed. The pressure above the bubble gradually increases due to the presence of a rising bubble. When the bubble passes through the bed surface, the pressure near the bed surface gradually decreases to a lower value. Moreover, the pressure signals in the bubbling fluidized beds show obviously periodic characteristics. The major frequency of pressure fluctuations at the same vertical position is affected slightly by the operating gas velocity, and the amplitude of pressure fluctuations is related to both the operating gas velocity and the vertical height. In this study, the influence of the operating gas velocity on the pressure wave propagation velocity can be ignored, and only two peak frequencies in the power spectrum of the pressure fluctuations are observed which are associated with the bubble formation above the distributor and its eruption at the bed surface.展开更多
The Issue of mixing efficiency in agitated tanks has drawn serious concern in many industrial processes. The turbulence model is very critical to predicting mixing process in agitated tanks. On the basis of computatio...The Issue of mixing efficiency in agitated tanks has drawn serious concern in many industrial processes. The turbulence model is very critical to predicting mixing process in agitated tanks. On the basis of computational fluid dynamics(CFD) software package Fluent 6.2, the mixing characteristics in a tank agitated by dual six-blade-Rushton-turbines(6-DT) are predicted using the detached eddy simulation(DES) method. A sliding mesh(SM) approach is adopted to solve the rotation of the impeller. The simulated flow patterns and liquid velocities in the agitated tank are verified by experimental data in the literature. The simulation results indicate that the DES method can obtain more flow details than Reynolds-averaged Navier-Stokes(RANS) model. Local and global mixing time in the agitated tank is predicted by solving a tracer concentration scalar transport equation. The simulated results show that feeding points have great influence on mixing process and mixing time. Mixing efficiency is the highest for the feeding point at location of midway of the two impellers. Two methods are used to determine global mixing time and get close result. Dimensionless global mixing time remains unchanged with increasing of impeller speed. Parallel, merging and diverging flow pattern form in the agitated tank, respectively, by changing the impeller spacing and clearance of lower impeller from the bottom of the tank. The global mixing time is the shortest for the merging flow, followed by diverging flow, and the longest for parallel flow. The research presents helpful references for design, optimization and scale-up of agitated tanks with multi-impeller.展开更多
In the radiant section of cracking furnace,the thermal cracking process is highly coupled with turbulent flow,heat transfer and mass transfer.In this paper,a three-dimensional simulation of propane pyrolysis reactor t...In the radiant section of cracking furnace,the thermal cracking process is highly coupled with turbulent flow,heat transfer and mass transfer.In this paper,a three-dimensional simulation of propane pyrolysis reactor tube is performed based on a detailed kinetic radical cracking scheme,combined with a comprehensive rigorous computational fluid dynamics(CFD)model.The eddy-dissipation-concept(EDC)model is introduced to deal with turbulence-chemistry interaction of cracking gas,especially for the multi-step radical kinetics.Considering the high aspect ratio and severe gradient phenomenon,numerical strategies such as grid resolution and refinement,stepping method and relaxation technique at different levels are employed to accelerate convergence.Large scale of radial nonuniformity in the vicinity of the tube wall is investigated.Spatial distributions of each radical reaction rate are first studied,and made it possible to identify the dominant elementary reactions.Additionally,a series of operating conditions including the feedstock feed rate,wall temperature profile and heat flux profile towards the reactor tubes are investigated.The obtained results can be used as scientific guide for further technical retrofit and operation optimization aiming at high conversion and selectivity of pyrolysis process.展开更多
Thermal comfort and indoor air quality as well as the energy efficiency have been recognized as essential parts of sustainable building assessment. This work aims to analyze the energy conservation of the heat recover...Thermal comfort and indoor air quality as well as the energy efficiency have been recognized as essential parts of sustainable building assessment. This work aims to analyze the energy conservation of the heat recovery ventilator and to investigate the effect of the air supply arrangement. Three types of mixing ventilation are chosen for the analysis of coupling ANSYS/FLUENT (a computational fluid dynamics (CFD) program) with TRNSYS (a building energy simulation (BES) software). The adoption of mutual complementary boundary conditions for CFD and BES provides more accurate and complete information of indoor air distribution and thermal performance in buildings. A typical office-space situated in a middle storey is chosen for the analysis. The office-space is equipped with air-conditioners on the ceiling. A heat recovery ventilation system directly supplies flesh air to the office space. Its thermal performance and indoor air distribution predicted by the coupled method are compared under three types of ventilation system. When the supply and return openings for ventilation are arranged on the ceiling, there is no critical difference between the predictions of the coupled method and BES on the energy consumption of HVAC because PID control is adopted for the supply air temperature of the occupied zone. On the other hand, approximately 21% discrepancy for the heat recovery estimation in the maximum between the simulated results of coupled method and BES-only can be obviously found in the floor air supply ventilation case. The discrepancy emphasizes the necessity of coupling CFD with BES when vertical air temperature gradient exists. Our future target is to estimate the optimum design of heat recovery ventilation system to control CO2 concentration by adjusting flow rate of flesh air.展开更多
基金supported by National Natural Science Foundation of China(No.52075334)。
文摘Field emission electric propulsion(FEEP) thrusters possess excellent characteristics, such as high specific impulse, low power requirements, compact size and precise pointing capabilities,making them ideal propulsion devices for micro-nano satellites. However, the detection of certain aspects, such as the evolution process of the liquid cone and the physical quantities at the cone apex, proves challenging due to the minute size of the needle tip and the vacuum environment in which they operate. Consequently, this paper introduces a computational fluid dynamics(CFD) model to gain insight into the formation process of the liquid cone on the tip apex of indium FEEP. The CFD model is based on electrohydrodynamic(EHD) equations and the volume of fluid(VOF) method. The entire cone formation process can be divided into three stages, and the time-dependent characteristics of the physical quantities at the cone apex are investigated. The influences of film thickness, apex radius size and applied voltage are compared.The results indicate a gradual increase in the values of electrostatic stress and surface tension stress at the cone apex over an initial period, followed by a rapid escalation within a short duration.Apex configurations featuring a small radius, thick film and high voltage exhibit a propensity for liquid cone formation, and the cone growth time decreases as the film thickness increases.Moreover, some unstable behavior is observed during the cone formation process.
基金the High-Performance Computing Center of Nanjing Tech University for supporting the computational resources
文摘Dry storage containers must be secured and reliable during long-term storage,and the effect of decay heat released from the internal spent fuel on the cask has become an important research topic.In this paper,a 3D computational fluid dynamics model is presented,and the accuracy of the calculation is verified,with computational errors of less than 6.2%.The thermal stress of the dry storage cask was estimated by coupling it with a transient temperature field.The total power remained constant and adjusting the power ratio of the inner and outer zones had a small effect on the stress results,with a maximum equivalent stress of approximately 5.2 kPa,which occurred at the lower edge of the shell.In the case of tilt,the temperature gradient varied in a wavy distribution,and the wave crest moved from right to left.Altering the tilt angle affects the air distribution in the annular gap,leading to the shell temperature being transformed,with a maximum equivalent stress of 202 MPa at the bottom of the shell.However,the equivalent stress in both cases was less than the yield stress(205 MPa).
基金support from the Major State Basic Research Development Program of China (973 Program,2005CB221205)National Natural Science Foundation of China (No.20490200 and 20576076)
文摘Computational Fluid Dynamics (CFD) has become an alternative method to experiments for understanding the fluid dynamics of multiphase flow. A two-fluid model, which contains additional terms in both the gas- and solid-phase momentum equations, is used to investigate the fluidization quality in a fluidized bed. A case study for quartz sand with a density of 2,660 kg/m^3 and a diameter of 500 μm, whose physical property is similar to a new kind of catalyst for producing clean fuels through the residue fluid catalytic cracking process, is simulated in a two-dimensional fluidized bed with 0.57 m width and 1.00 m height. Transient bubbling and collapsing characteristics are numerically investigated in the platform of CFX 4.4 by integrating user-defined Fortran subroutines. The results show that the fluidization and collapse process is in fair agreement with the classical theory of Geldart B classification, but the collapse time is affected by bubbles at the interface between the dense phase and freeboard.
基金The National Science and Technology Major Project of China under contract No.2011ZX05025-002-02the National Natural Science Foundation of China under contract Nos 41476032,91028009 and 40806019
文摘Turbidity currents represent a major agent for sediment transport in lakes, seas and oceans. In particu-lar, they formulate the most significant clastic accumulations in the deep sea, which become many of the world's most important hydrocarbon reservoirs. Several boreholes in the Qiongdongnan Basin, the north-western South China Sea, have recently revealed turbidity current deposits as significant hydrocarbon res-ervoirs. However, there are some arguments for the potential provenances. To solve this problem, it is es-sential to delineate their sedimentary processes as well as to evaluate their qualities as reservoir. Numerical simulations have been developed rapidly over the last several years, offering insights into turbidity current behaviors, as geologically significant turbidity currents are difficult to directly investigate due to their large scale and often destructive nature. Combined with the interpretation of the turbidity system based on high-resolution 3D seismic data, the paleotophography is acquired via a back-stripping seismic profile integrated with a borehole, i.e., Well A, in the western Qiongdongnan Basin; then a numerical model is built on the basis of this back-stripped profile. After defining the various turbidity current initial boundary conditions, includ-ing grain size, velocity and sediment concentration, the structures and behaviors of turbidity currents are investigated via numerical simulation software ANSYS FLUENT. Finally, the simulated turbidity deposits are compared with the interpreted sedimentary bodies based on 3D seismic data and the potential provenances of the revealed turbidites by Well A are discussed in details. The simulation results indicate that a sedimen-tary body develops far away from its source with an average grain size of 0.1 mm, i.e., sand-size sediment. Taking into account the location and orientation of the simulated seismic line, the consistence between normal forward simulation results and the revealed cores in Well A indicates that the turbidites should have been transported from Vietnam instead of Hainan Island. This interpretation has also been verified by the planar maps of sedimentary systems based on integration of boreholes and seismic data. The identification of the turbidity provenance will benefit the evaluation of extensively distributed submarine fans for hydro-carbon exploration in the deepwater areas.
基金National Natural Science Foundation of China (No.50435030)
文摘Computational fluid dynamics(CFD) simulations are adopted to investigate rectangular microchannel flows with various periodic micro-structured wall by introducing velocity slip boundary condition at low Reynolds number. The purpose of the current study is to numerically find out the effects of periodic micro-structured wall on the flow resistance in rectangular microchannel with the different spacings between microridges ranging from 15 to 60 pm. The simulative results indicate that pressure drop with different spacing between microridges increases linearly with flow velocity and decreases monotonically with slip velocity; Pressure drop reduction also increases with the spacing between microridges at the same condition of slip velocity and flow velocity. The results of numerical simulation are compared with theoretical predictions and experimental results in the literatures. It is found that there is qualitative agreement between them.
基金support from National Basic Research Program of China(No.2009CB219801)National Natural Science Foundation of China(No.20976191)+1 种基金International Cooperative Program of Guizhou Province([2009]700110)Program for New Century Excellent Talents in University(NCET-09-0342)
文摘A simple hydrodynamic model based on two-fluid theory, taking into account the effect of discrete particles on both the gas- and solid-phase momentum equations, was used to numerically investigate the pressure fluctuation characteristics in a gas-solid fluidized bed with the aid of CFX 4.4, a commercial CFD software package, by adding user-defined Fortran subroutines. Numerical simulations together with typical experimental measurements show that pressure fluctuations originate above the distributor when a gas pulse is injected into the fluidized bed. The pressure above the bubble gradually increases due to the presence of a rising bubble. When the bubble passes through the bed surface, the pressure near the bed surface gradually decreases to a lower value. Moreover, the pressure signals in the bubbling fluidized beds show obviously periodic characteristics. The major frequency of pressure fluctuations at the same vertical position is affected slightly by the operating gas velocity, and the amplitude of pressure fluctuations is related to both the operating gas velocity and the vertical height. In this study, the influence of the operating gas velocity on the pressure wave propagation velocity can be ignored, and only two peak frequencies in the power spectrum of the pressure fluctuations are observed which are associated with the bubble formation above the distributor and its eruption at the bed surface.
基金Supported by Key Scientific Research Project of Sichuan Provincial Education Department(Grant No.15ZA0107)Doctor Foundation of Southwest University of Science and Technology(Grant No.11zx7162)
文摘The Issue of mixing efficiency in agitated tanks has drawn serious concern in many industrial processes. The turbulence model is very critical to predicting mixing process in agitated tanks. On the basis of computational fluid dynamics(CFD) software package Fluent 6.2, the mixing characteristics in a tank agitated by dual six-blade-Rushton-turbines(6-DT) are predicted using the detached eddy simulation(DES) method. A sliding mesh(SM) approach is adopted to solve the rotation of the impeller. The simulated flow patterns and liquid velocities in the agitated tank are verified by experimental data in the literature. The simulation results indicate that the DES method can obtain more flow details than Reynolds-averaged Navier-Stokes(RANS) model. Local and global mixing time in the agitated tank is predicted by solving a tracer concentration scalar transport equation. The simulated results show that feeding points have great influence on mixing process and mixing time. Mixing efficiency is the highest for the feeding point at location of midway of the two impellers. Two methods are used to determine global mixing time and get close result. Dimensionless global mixing time remains unchanged with increasing of impeller speed. Parallel, merging and diverging flow pattern form in the agitated tank, respectively, by changing the impeller spacing and clearance of lower impeller from the bottom of the tank. The global mixing time is the shortest for the merging flow, followed by diverging flow, and the longest for parallel flow. The research presents helpful references for design, optimization and scale-up of agitated tanks with multi-impeller.
基金Supported by the National Science&Technology Supporting Plan(2012BAF05B00)the National Basic Research Program(2012CB720500)
文摘In the radiant section of cracking furnace,the thermal cracking process is highly coupled with turbulent flow,heat transfer and mass transfer.In this paper,a three-dimensional simulation of propane pyrolysis reactor tube is performed based on a detailed kinetic radical cracking scheme,combined with a comprehensive rigorous computational fluid dynamics(CFD)model.The eddy-dissipation-concept(EDC)model is introduced to deal with turbulence-chemistry interaction of cracking gas,especially for the multi-step radical kinetics.Considering the high aspect ratio and severe gradient phenomenon,numerical strategies such as grid resolution and refinement,stepping method and relaxation technique at different levels are employed to accelerate convergence.Large scale of radial nonuniformity in the vicinity of the tube wall is investigated.Spatial distributions of each radical reaction rate are first studied,and made it possible to identify the dominant elementary reactions.Additionally,a series of operating conditions including the feedstock feed rate,wall temperature profile and heat flux profile towards the reactor tubes are investigated.The obtained results can be used as scientific guide for further technical retrofit and operation optimization aiming at high conversion and selectivity of pyrolysis process.
基金Project supported by Grant-in-Aid for Scientific Research (JSPS KAKENHI for Young Scientists (S), 21676005)
文摘Thermal comfort and indoor air quality as well as the energy efficiency have been recognized as essential parts of sustainable building assessment. This work aims to analyze the energy conservation of the heat recovery ventilator and to investigate the effect of the air supply arrangement. Three types of mixing ventilation are chosen for the analysis of coupling ANSYS/FLUENT (a computational fluid dynamics (CFD) program) with TRNSYS (a building energy simulation (BES) software). The adoption of mutual complementary boundary conditions for CFD and BES provides more accurate and complete information of indoor air distribution and thermal performance in buildings. A typical office-space situated in a middle storey is chosen for the analysis. The office-space is equipped with air-conditioners on the ceiling. A heat recovery ventilation system directly supplies flesh air to the office space. Its thermal performance and indoor air distribution predicted by the coupled method are compared under three types of ventilation system. When the supply and return openings for ventilation are arranged on the ceiling, there is no critical difference between the predictions of the coupled method and BES on the energy consumption of HVAC because PID control is adopted for the supply air temperature of the occupied zone. On the other hand, approximately 21% discrepancy for the heat recovery estimation in the maximum between the simulated results of coupled method and BES-only can be obviously found in the floor air supply ventilation case. The discrepancy emphasizes the necessity of coupling CFD with BES when vertical air temperature gradient exists. Our future target is to estimate the optimum design of heat recovery ventilation system to control CO2 concentration by adjusting flow rate of flesh air.