The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the n...The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment （SOM） combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment （ASOM） combustion model is well validated by DNS.展开更多
On the basis of the full velocity difference （FVD） model, an improved multiple car-following （MCF） model is proposed by taking into account multiple information inputs from preceding vehicles. The linear stability...On the basis of the full velocity difference （FVD） model, an improved multiple car-following （MCF） model is proposed by taking into account multiple information inputs from preceding vehicles. The linear stability condition of the model is obtained by using the linear stability theory. Through nonlinear analysis, a modified Korteweg-de Vries equation is constructed and solved. The traffic jam can thus be described by the klnk-antikink soliton solution for the mKdV equation. The improvement of this new model over the previous ones lies in the fact that it not only theoretically retains many strong points of the previous ones, but also performs more realistically than others in the dynamical evolution of congestion. Furthermore, numerical simulation of traffic dynamics shows that the proposed model can avoid the disadvantage of negative velocity that occurs at small sensitivity coefficients λ in the FVD model by adjusting the information on the multiple leading vehicles. No collision occurs and no unrealistic deceleration appears in the improved model.展开更多
A mathematical model for coupled multiphase fluid flow and sedi- mentation deformation is developed based on fluid-solid interaction mechanism.A finite difference-finite element numerical approach is presented.The res...A mathematical model for coupled multiphase fluid flow and sedi- mentation deformation is developed based on fluid-solid interaction mechanism.A finite difference-finite element numerical approach is presented.The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances,and the coupled model has practical significance for oilfield development.展开更多
Based on Xue＇s lattice model, an extended lattice model is proposed by considering the relative current information about next-nearest-neighbour sites ahead. The linear stability condition of the presented model is o...Based on Xue＇s lattice model, an extended lattice model is proposed by considering the relative current information about next-nearest-neighbour sites ahead. The linear stability condition of the presented model is obtained by employing the linear stability theory. The density wave is investigated analytically with the perturbation method. The results show that the occurrence of traffic jamming transitions can be described by the kink-antikink solution of the modified Korteweg-de Vries （mKdV） equation. The simulation results are in good agreement with the analytical results, showing that the stability of traffic flow can be enhanced when the relative current of next-nearest-neighbour sites ahead is considered.展开更多
It is assumed in this paper that for a high Reynolds number nearly homogeneouswind flow, the Reynolds stresses are uniquely related to the mean velocity gradientsand the two independent turbulent scaling parameters k ...It is assumed in this paper that for a high Reynolds number nearly homogeneouswind flow, the Reynolds stresses are uniquely related to the mean velocity gradientsand the two independent turbulent scaling parameters k and E. By applying dimensionalanalysis and owing to the Cayley-Hamilton theorem for tensors, a new turbulenceenclosure model so-called the axtended k-ε model has been developed. The coefficientsof the model expression were detemined by the wind tunnel experimental data ofhomogeneous shear turbulent flow. The model was compared with the standard k-εmodel in in composition and the prediction of the Reynold's normal Stresses. Using thenew model the numerical simulation of wind flow around a square cross-section tallbuilding was performed. The results show that the extended k-ε model improves theprediction of wind velocities around the building the building and wind pressures on the buildingenvelope.展开更多
Direct numerical simulation of decaying homogeneous isotropic turbulence （DHIT） of a polymer solution is performed. In order to understand the polymer effect on turbulence or additive-turbulence interaction, we dire...Direct numerical simulation of decaying homogeneous isotropic turbulence （DHIT） of a polymer solution is performed. In order to understand the polymer effect on turbulence or additive-turbulence interaction, we directly investigate the influence of polymers on velocity gradient tensor including vorticity and strain. By visualizing vortex tubes and sheets, we observe a remarkable inhibition of vortex structures in an intermediate-scale field and a small-scale field but not for a large scale field in DHIT with polymers. The geometric study indicates a strong relevance among the vorticity vector, rate-of-strain tensor, and polymer conformation tensor. Joint probability density functions show that the polymer effect can increase ＂strain generation resistance＂ and ＂vorticity generation resistance＂, i.e., inhibit the generation of vortex sheets and tubes, ultimately leading to turbulence inhibition effects.展开更多
The 3-dimensional fluid flow in a water model of the continuous casting tundish is simulated with the k-ε two-equation turbulence model. The methods to decide the size of wall-adjacent grids and the effects of resid...The 3-dimensional fluid flow in a water model of the continuous casting tundish is simulated with the k-ε two-equation turbulence model. The methods to decide the size of wall-adjacent grids and the effects of residuals and mesh size on the simulation accuracy are discussed. The current investigation concludes that the following condition should be satisfied to get accurate enough simulation results: (1) If the dimension of the domain is in the order of cubic meters, the average size of cells in the mesh system should be at least smaller than 30 mm; (2) The normalized non-scaled residual should be reached at least smaller than 10-4.展开更多
Numerical simulation on the flow,heat transfer and cracking reactions in commercial fluid catalyticcracking(FCC)riser reactors were carried out employing the developed turbulent gas-solid two-phase flow-reac-tion mode...Numerical simulation on the flow,heat transfer and cracking reactions in commercial fluid catalyticcracking(FCC)riser reactors were carried out employing the developed turbulent gas-solid two-phase flow-reac-tion model for FCC riser reactors given in Part Ⅰ of the present paper.Detailed information about the turbulentflow fields in the riser reactor obtained revealed the basic characteristics of the gas-solid two-phase turbulentflows when heat transfer and catalytic cracking reactions were co-existing in the riser.Results showed that thedistributions of the flow,the turbulence kinetic energy and the catalyst particle concentration are not uniform inthe axial,radial and tangential directions.The most complicated part of the riser reactor is the feed injectingzone.The complicated configuration of the turbulent gas-solid two-phase flows would exert a great influence onthe results of interphase heat transfer and cracking reactions.展开更多
A nonlinear flow reservoir mathematical model was established based on the flow characteristic of low-permeability reservoir.The well-grid equations were deduced and the dimensionless permeability coefficient was intr...A nonlinear flow reservoir mathematical model was established based on the flow characteristic of low-permeability reservoir.The well-grid equations were deduced and the dimensionless permeability coefficient was introduced to describe the permeability variation of nonlinear flow.The nonlinear flow numerical simulation program was compiled based on black-oil model.A quarter of five-spot well unit was simulated to study the effect of nonlinear flow on the exploitation of low-permeability reservoir.The comprehensive comparison and analysis of the simulation results of Darcy flow,quasi-linear flow and nonlinear flow were provided.The dimensionless permeability coefficient distribution was gained to describe the nonlinear flow degree.The result shows that compared with the results of Darcy flow,when considering nonlinear flow,the oil production is low,and production decline is rapid.The fluid flow in reservoir consumes more driving energy,which reduces the water flooding efficiency.Darcy flow model overstates the reservoir flow capability,and quasi-linear flow model overstates the reservoir flow resistance.The flow ability of the formation near the well and artificial fracture is strong while the flow ability of the formation far away from the main streamline is weak.The nonlinear flow area is much larger than that of quasi-linear flow during the fluid flow in low-permeability reservoir.The water propelling speed of nonlinear flow is greatly slower than that of Darcy flow in the vertical direction of artificial fracture,and the nonlinear flow should be taken into account in the well pattern arrangement of low-permeability reservoir.展开更多
Simulation of flow in fractured aquifers is a complex issue. The problem of mathematical modeling of highly heterogeneous porous media, typical of natural systems, couples with the needs making proper simplifying assu...Simulation of flow in fractured aquifers is a complex issue. The problem of mathematical modeling of highly heterogeneous porous media, typical of natural systems, couples with the needs making proper simplifying assumptions and approximations. In this kind of groundwater systems, studying contamination spreading and analyzing risk are challenging tasks. The main difficulty stems from determining both the travel times and the maximum distances covered by pollutants. In this context, the risk of contamination in the deep carbonate aquifer of the Augusta coastal area is presented. We used a geostatistical approach and numerical codes (MODFLOW-2000, MT3DMS) to reconstruct the complex geological framework of the study area where several contamination scenarios of hypothetical point source in a risk assessment framework were simulated. Results on the contaminant spreading are discussed and the effect of the geological structures characterizing the zone under study, namely horst and graben, are described.展开更多
A three-dimensional numerical model is presented for studying the convection-condensation of mixture with vapor in a tube with edgefold-twisted-tape inserts under transition flow.According to the diffusion layer theor...A three-dimensional numerical model is presented for studying the convection-condensation of mixture with vapor in a tube with edgefold-twisted-tape inserts under transition flow.According to the diffusion layer theory and laminar species transport,a condensation model with user defined function is proposed and compared with heat and mass transfer analogy and experimental test.With the condensation model,the influences of gap width and op-erating parameters on thermal-hydrodynamics performance are simulated.As the gap width increases,convection and condensation heat transfer increase initially and then decrease,while convection heat transfer increases sharply and then decreases slightly.Increasing vapor fraction has a significant effect on condensation heat transfer but it has little effect on convective heat transfer.With the increase of inner wall temperature both convection and condensa-tion heat transfer all decrease and the ratio of condensation to total heat decrease dramatically.Increases inlet tem-perature mainly affects convection heat transfer.展开更多
The conventional Shear Stress Transport(SST)k–ωturbulence model often exhibits substantial inaccu-racies when applied to the prediction of flow behavior in complex regions within axial flow control valves.To enhance...The conventional Shear Stress Transport(SST)k–ωturbulence model often exhibits substantial inaccu-racies when applied to the prediction of flow behavior in complex regions within axial flow control valves.To enhance its predictive fidelity for internal flow fields,this study introduces a novel calibration framework that integrates an artificial neural network(ANN)surrogate model with a particle swarm optimization(PSO)algorithm.In particular,an optimal Latin hypercube sampling strategy was employed to generate representative sample points across the empirical parameter space.For each sample,numerical simulations using ANSYS Fluent were conducted to evaluate the flow characteristics,with empirical turbulence model parameters as inputs and flow rate as the target output.These data were used to construct the high-fidelity ANN surrogate model.The PSO algorithm was then applied to this surrogate to identify the optimal set of empirical parameters tailored specifically to axial flow control valve configurations.A revealed by the presented results,the calibrated SST k–ωmodel significantly improves prediction accuracy:deviations from large eddy simulation(LES)benchmarks at small valve openings were reduced from 7.6%to under 3%.Furthermore,the refined model maintains the computational efficiency characteristic of Reynolds-averaged Navier-Stokes(RANS)simulations while substantially enhancing the accuracy of both pressure and velocity field predictions.Overall,the proposed methodology effectively reconciles the trade-off between computational cost and predictive accuracy,offering a robust and scalable approach for turbulence model calibration in complex internal flow scenarios.展开更多
Wet steam condensing flow in low-pressure steam turbine leads to efficiencylosses and blade erosions. In order to investigate this problem by numerical approach, anEulerian/Eulerian model has been developed, in which ...Wet steam condensing flow in low-pressure steam turbine leads to efficiencylosses and blade erosions. In order to investigate this problem by numerical approach, anEulerian/Eulerian model has been developed, in which the wet steam is regarded as mixture comprisingtwo coupled systems: the vapor phase and the liquid phase. These two systems are both described byconservation equations. High resolution TVD scheme is employed to capture condensing phenomena inwet steam flow. This model has been validated by numerical simulations of condensing flows in 1D and2D nozzles. Compared with experimental data, a good agreement is observed. This Eulerian/Eulerianmodel can be extended to 3D calculation of condensing flow.展开更多
ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged u...ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged unsteady flow field.Meanwhile,drag and lift coefficients of the four different low-Reynolds number turbulence models were analyzed.The simulated results of YANG-SHIH model are close to the large eddy simulation results and experimental results,and they are significantly better than those of ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMR models.The modification of the generation of turbulence kinetic energy is the key factor to a successful simulation for YANG-SHIH model,while the correction of the turbulence near the wall has minor influence on the simulation results.For ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMA models satisfactory simulation results cannot be obtained due to lack of the modification of the generation of turbulence kinetic energy.With the joint force of wall function and the turbulence models with the adoption of corrected swirl stream,flow around a square cylinder can be fully simulated with less grids by the near-wall.展开更多
Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster managem...Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.展开更多
In order to evaluate CCFL （countercurrent flow limitation） characteristics in a PWR （pressurized water reactor） hot leg under reflux condensation, numerical simulations have been conducted using a 2F （two-fluid）...In order to evaluate CCFL （countercurrent flow limitation） characteristics in a PWR （pressurized water reactor） hot leg under reflux condensation, numerical simulations have been conducted using a 2F （two-fluid） model and a VOF （volume of fluid） method implemented in the CFD （computational fluid dynamics） software, FLUENT6.3.26. The 2F model gave good agreement with CCFL data in low pressure conditions but did not give good results for high pressure steam-water conditions. In the previous study, the computational grid and schemes were improved in the VOF method to improve calculations in circular tubes, and the calculated CCFL characteristics agreed well with the UPTF （Upper Plenum Test Facility） data at 1.5 MPa. In this study, therefore, using the 2F model and the computational grid previously improved for the VOF calculations, numerical simulations were conducted for steam-water flows at 1.5 MPa under PWR full-scale conditions. In the range of medium gas volumetric fluxes, the calculated CCFL characteristics agreed well with the values calculated by the VOF method and the UPTF data at 1.5 MPa. This indicated that the reference set of the interfacial drag correlations employed in this study could be applied not only to low pressures but also to high pressures.展开更多
Dense gas-solid flows are very common in actual production and industrial fields,so it is significant to understand their hydrodynamic characteristics and heat and mass transfer behaviors.This article provides a brief...Dense gas-solid flows are very common in actual production and industrial fields,so it is significant to understand their hydrodynamic characteristics and heat and mass transfer behaviors.This article provides a brief review of multi-scale numerical simulation of flow,heat and mass transfer behaviors in dense gas-solid flows.It describes multiscale models(direct numerical simulation,discrete particle model,and two-fluid model)and the results of related research.Finally,it discusses possible future developments in research on the flow,heat and mass transfer characteristics of dense gas-solid two-phase flows.展开更多
Slurry jets in a static uniform environment were simulated with a two-phase mixture model in which flow-particle interactions were considered. A standard k-e turbulence model was chosen to close the governing equation...Slurry jets in a static uniform environment were simulated with a two-phase mixture model in which flow-particle interactions were considered. A standard k-e turbulence model was chosen to close the governing equations. The computational results were in agreement with previous laboratory measurements. The characteristics of the two-phase flow field and the influences of hydraulic and geometric parameters on the distribution of the slurry jets were analyzed on the basis of the computational results. The calculated results reveal that if the initial velocity of the slurry jet is high, the jet spreads less in the radial direction. When the slurry jet is less influenced by the ambient fluid （when the Stokes number St is relatively large）, the turbulent kinetic energy k and turbulent dissipation rate e, which are relatively concentrated around the jet axis, decrease more rapidly after the slurry jet passes through the nozzle. For different values of St, the radial distributions of streamwise velocity and particle volume fraction are both self-similar and fit a Gaussian profile after the slurry jet fully develops. The decay rate of the particle velocity is lower than that of water velocity along the jet axis, and the axial distributions of the centerline particle streamwise velocity are self-similar along the jet axis. The pattern of particle dispersion depends on the Stokes number St. When St = 0.39, the panicle dispersion along the radial direction is considerable, and the relative velocity is very low due to the low dynamic response time. When St = 3.08, the dispersion of particles along the radial direction is very little, and most of the particles have high relative velocities along the streamwise direction.展开更多
An accurate critical heat flux（CHF） prediction method is the key factor for realizing the steady-state operation of a water-cooled divertor that works under one-sided high heating flux conditions.An improved CHF pre...An accurate critical heat flux（CHF） prediction method is the key factor for realizing the steady-state operation of a water-cooled divertor that works under one-sided high heating flux conditions.An improved CHF prediction method based on Euler＇s homogeneous model for flow boiling combined with realizable k-ε model for single-phase flow is adopted in this paper in which time relaxation coefficients are corrected by the Hertz-Knudsen formula in order to improve the calculation accuracy of vapor-liquid conversion efficiency under high heating flux conditions.Moreover,local large differences of liquid physical properties due to the extreme nonuniform heating flux on cooling wall along the circumference direction are revised by formula IAPWSIF97.Therefore,this method can improve the calculation accuracy of heat and mass transfer between liquid phase and vapor phase in a CHF prediction simulation of water-cooled divertors under the one-sided high heating condition.An experimental example is simulated based on the improved and the uncorrected methods.The simulation results,such as temperature,void fraction and heat transfer coefficient,are analyzed to achieve the CHF prediction.The results show that the maximum error of CHF based on the improved method is 23.7%,while that of CHF based on uncorrected method is up to 188%,as compared with the experiment results of Ref.[12].Finally,this method is verified by comparison with the experimental data obtained by International Thermonuclear Experimental Reactor（ITER）,with a maximum error of 6% only.This method provides an efficient tool for the CHF prediction of water-cooled divertors.展开更多
基金Supported by the National Natural Science Foundation of China （50606026, 50736006）.
文摘The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment （SOM） combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment （ASOM） combustion model is well validated by DNS.
基金Project supported by the National High Tech Research and Development Program of China (Grant No 511-0910-1031)the National "10th Five-year" Science and Technique Important Program of China (Grant No 2002BA404A07)
文摘On the basis of the full velocity difference （FVD） model, an improved multiple car-following （MCF） model is proposed by taking into account multiple information inputs from preceding vehicles. The linear stability condition of the model is obtained by using the linear stability theory. Through nonlinear analysis, a modified Korteweg-de Vries equation is constructed and solved. The traffic jam can thus be described by the klnk-antikink soliton solution for the mKdV equation. The improvement of this new model over the previous ones lies in the fact that it not only theoretically retains many strong points of the previous ones, but also performs more realistically than others in the dynamical evolution of congestion. Furthermore, numerical simulation of traffic dynamics shows that the proposed model can avoid the disadvantage of negative velocity that occurs at small sensitivity coefficients λ in the FVD model by adjusting the information on the multiple leading vehicles. No collision occurs and no unrealistic deceleration appears in the improved model.
文摘A mathematical model for coupled multiphase fluid flow and sedi- mentation deformation is developed based on fluid-solid interaction mechanism.A finite difference-finite element numerical approach is presented.The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances,and the coupled model has practical significance for oilfield development.
基金Project supported by the National High Technology Research and Development Program of China(Grant No.511-0910-1031)
文摘Based on Xue＇s lattice model, an extended lattice model is proposed by considering the relative current information about next-nearest-neighbour sites ahead. The linear stability condition of the presented model is obtained by employing the linear stability theory. The density wave is investigated analytically with the perturbation method. The results show that the occurrence of traffic jamming transitions can be described by the kink-antikink solution of the modified Korteweg-de Vries （mKdV） equation. The simulation results are in good agreement with the analytical results, showing that the stability of traffic flow can be enhanced when the relative current of next-nearest-neighbour sites ahead is considered.
文摘It is assumed in this paper that for a high Reynolds number nearly homogeneouswind flow, the Reynolds stresses are uniquely related to the mean velocity gradientsand the two independent turbulent scaling parameters k and E. By applying dimensionalanalysis and owing to the Cayley-Hamilton theorem for tensors, a new turbulenceenclosure model so-called the axtended k-ε model has been developed. The coefficientsof the model expression were detemined by the wind tunnel experimental data ofhomogeneous shear turbulent flow. The model was compared with the standard k-εmodel in in composition and the prediction of the Reynold's normal Stresses. Using thenew model the numerical simulation of wind flow around a square cross-section tallbuilding was performed. The results show that the extended k-ε model improves theprediction of wind velocities around the building the building and wind pressures on the buildingenvelope.
基金supported by the National Natural Science Foundation of China (Grant No. 10872060)the Fundamental Research Funds for the Central Universities (Grant No. HIT.BRET2.2010008)
文摘Direct numerical simulation of decaying homogeneous isotropic turbulence （DHIT） of a polymer solution is performed. In order to understand the polymer effect on turbulence or additive-turbulence interaction, we directly investigate the influence of polymers on velocity gradient tensor including vorticity and strain. By visualizing vortex tubes and sheets, we observe a remarkable inhibition of vortex structures in an intermediate-scale field and a small-scale field but not for a large scale field in DHIT with polymers. The geometric study indicates a strong relevance among the vorticity vector, rate-of-strain tensor, and polymer conformation tensor. Joint probability density functions show that the polymer effect can increase ＂strain generation resistance＂ and ＂vorticity generation resistance＂, i.e., inhibit the generation of vortex sheets and tubes, ultimately leading to turbulence inhibition effects.
文摘The 3-dimensional fluid flow in a water model of the continuous casting tundish is simulated with the k-ε two-equation turbulence model. The methods to decide the size of wall-adjacent grids and the effects of residuals and mesh size on the simulation accuracy are discussed. The current investigation concludes that the following condition should be satisfied to get accurate enough simulation results: (1) If the dimension of the domain is in the order of cubic meters, the average size of cells in the mesh system should be at least smaller than 30 mm; (2) The normalized non-scaled residual should be reached at least smaller than 10-4.
文摘Numerical simulation on the flow,heat transfer and cracking reactions in commercial fluid catalyticcracking(FCC)riser reactors were carried out employing the developed turbulent gas-solid two-phase flow-reac-tion model for FCC riser reactors given in Part Ⅰ of the present paper.Detailed information about the turbulentflow fields in the riser reactor obtained revealed the basic characteristics of the gas-solid two-phase turbulentflows when heat transfer and catalytic cracking reactions were co-existing in the riser.Results showed that thedistributions of the flow,the turbulence kinetic energy and the catalyst particle concentration are not uniform inthe axial,radial and tangential directions.The most complicated part of the riser reactor is the feed injectingzone.The complicated configuration of the turbulent gas-solid two-phase flows would exert a great influence onthe results of interphase heat transfer and cracking reactions.
基金Project(10672187) supported by the National Natural Science Foundation of ChinaProject(2008ZX05000-013-02) supported by the National Science and Technology Major Program of China
文摘A nonlinear flow reservoir mathematical model was established based on the flow characteristic of low-permeability reservoir.The well-grid equations were deduced and the dimensionless permeability coefficient was introduced to describe the permeability variation of nonlinear flow.The nonlinear flow numerical simulation program was compiled based on black-oil model.A quarter of five-spot well unit was simulated to study the effect of nonlinear flow on the exploitation of low-permeability reservoir.The comprehensive comparison and analysis of the simulation results of Darcy flow,quasi-linear flow and nonlinear flow were provided.The dimensionless permeability coefficient distribution was gained to describe the nonlinear flow degree.The result shows that compared with the results of Darcy flow,when considering nonlinear flow,the oil production is low,and production decline is rapid.The fluid flow in reservoir consumes more driving energy,which reduces the water flooding efficiency.Darcy flow model overstates the reservoir flow capability,and quasi-linear flow model overstates the reservoir flow resistance.The flow ability of the formation near the well and artificial fracture is strong while the flow ability of the formation far away from the main streamline is weak.The nonlinear flow area is much larger than that of quasi-linear flow during the fluid flow in low-permeability reservoir.The water propelling speed of nonlinear flow is greatly slower than that of Darcy flow in the vertical direction of artificial fracture,and the nonlinear flow should be taken into account in the well pattern arrangement of low-permeability reservoir.
文摘Simulation of flow in fractured aquifers is a complex issue. The problem of mathematical modeling of highly heterogeneous porous media, typical of natural systems, couples with the needs making proper simplifying assumptions and approximations. In this kind of groundwater systems, studying contamination spreading and analyzing risk are challenging tasks. The main difficulty stems from determining both the travel times and the maximum distances covered by pollutants. In this context, the risk of contamination in the deep carbonate aquifer of the Augusta coastal area is presented. We used a geostatistical approach and numerical codes (MODFLOW-2000, MT3DMS) to reconstruct the complex geological framework of the study area where several contamination scenarios of hypothetical point source in a risk assessment framework were simulated. Results on the contaminant spreading are discussed and the effect of the geological structures characterizing the zone under study, namely horst and graben, are described.
基金Supported by the Technology Development Program of Jinan City (201102039,201202087)the Technology Development Program of Shandong Province (2011GNC11401)
文摘A three-dimensional numerical model is presented for studying the convection-condensation of mixture with vapor in a tube with edgefold-twisted-tape inserts under transition flow.According to the diffusion layer theory and laminar species transport,a condensation model with user defined function is proposed and compared with heat and mass transfer analogy and experimental test.With the condensation model,the influences of gap width and op-erating parameters on thermal-hydrodynamics performance are simulated.As the gap width increases,convection and condensation heat transfer increase initially and then decrease,while convection heat transfer increases sharply and then decreases slightly.Increasing vapor fraction has a significant effect on condensation heat transfer but it has little effect on convective heat transfer.With the increase of inner wall temperature both convection and condensa-tion heat transfer all decrease and the ratio of condensation to total heat decrease dramatically.Increases inlet tem-perature mainly affects convection heat transfer.
基金funded by Gansu Provincial Department of Education(Industrial Support Plan Project:2025CYZC-048).
文摘The conventional Shear Stress Transport(SST)k–ωturbulence model often exhibits substantial inaccu-racies when applied to the prediction of flow behavior in complex regions within axial flow control valves.To enhance its predictive fidelity for internal flow fields,this study introduces a novel calibration framework that integrates an artificial neural network(ANN)surrogate model with a particle swarm optimization(PSO)algorithm.In particular,an optimal Latin hypercube sampling strategy was employed to generate representative sample points across the empirical parameter space.For each sample,numerical simulations using ANSYS Fluent were conducted to evaluate the flow characteristics,with empirical turbulence model parameters as inputs and flow rate as the target output.These data were used to construct the high-fidelity ANN surrogate model.The PSO algorithm was then applied to this surrogate to identify the optimal set of empirical parameters tailored specifically to axial flow control valve configurations.A revealed by the presented results,the calibrated SST k–ωmodel significantly improves prediction accuracy:deviations from large eddy simulation(LES)benchmarks at small valve openings were reduced from 7.6%to under 3%.Furthermore,the refined model maintains the computational efficiency characteristic of Reynolds-averaged Navier-Stokes(RANS)simulations while substantially enhancing the accuracy of both pressure and velocity field predictions.Overall,the proposed methodology effectively reconciles the trade-off between computational cost and predictive accuracy,offering a robust and scalable approach for turbulence model calibration in complex internal flow scenarios.
基金This project is supported by National Natural Science Foundation of China (No.50176035)Doctorate Foundation of Xi'an Jiaotong University,China (No.DFXJU2000-18)
文摘Wet steam condensing flow in low-pressure steam turbine leads to efficiencylosses and blade erosions. In order to investigate this problem by numerical approach, anEulerian/Eulerian model has been developed, in which the wet steam is regarded as mixture comprisingtwo coupled systems: the vapor phase and the liquid phase. These two systems are both described byconservation equations. High resolution TVD scheme is employed to capture condensing phenomena inwet steam flow. This model has been validated by numerical simulations of condensing flows in 1D and2D nozzles. Compared with experimental data, a good agreement is observed. This Eulerian/Eulerianmodel can be extended to 3D calculation of condensing flow.
基金Project(2006BAJ04B04)supported by the National Science and Technology Pillar Program in the Eleventh Five-year Plan PeriodProject(2006AA05Z229)supported by the National High Technology Research and Development Program of China+1 种基金Project supportedby the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education MinistryProject(06wk3023)supported by Hunan Science and Technology Office
文摘ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged unsteady flow field.Meanwhile,drag and lift coefficients of the four different low-Reynolds number turbulence models were analyzed.The simulated results of YANG-SHIH model are close to the large eddy simulation results and experimental results,and they are significantly better than those of ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMR models.The modification of the generation of turbulence kinetic energy is the key factor to a successful simulation for YANG-SHIH model,while the correction of the turbulence near the wall has minor influence on the simulation results.For ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMA models satisfactory simulation results cannot be obtained due to lack of the modification of the generation of turbulence kinetic energy.With the joint force of wall function and the turbulence models with the adoption of corrected swirl stream,flow around a square cylinder can be fully simulated with less grids by the near-wall.
基金financially supported by Department of Space,India(Grant No.ISRO/RES/4/663/18-19)。
文摘Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.
文摘In order to evaluate CCFL （countercurrent flow limitation） characteristics in a PWR （pressurized water reactor） hot leg under reflux condensation, numerical simulations have been conducted using a 2F （two-fluid） model and a VOF （volume of fluid） method implemented in the CFD （computational fluid dynamics） software, FLUENT6.3.26. The 2F model gave good agreement with CCFL data in low pressure conditions but did not give good results for high pressure steam-water conditions. In the previous study, the computational grid and schemes were improved in the VOF method to improve calculations in circular tubes, and the calculated CCFL characteristics agreed well with the UPTF （Upper Plenum Test Facility） data at 1.5 MPa. In this study, therefore, using the 2F model and the computational grid previously improved for the VOF calculations, numerical simulations were conducted for steam-water flows at 1.5 MPa under PWR full-scale conditions. In the range of medium gas volumetric fluxes, the calculated CCFL characteristics agreed well with the values calculated by the VOF method and the UPTF data at 1.5 MPa. This indicated that the reference set of the interfacial drag correlations employed in this study could be applied not only to low pressures but also to high pressures.
基金financially supported by the National Natural Science Foundation of China(U20A20304)。
文摘Dense gas-solid flows are very common in actual production and industrial fields,so it is significant to understand their hydrodynamic characteristics and heat and mass transfer behaviors.This article provides a brief review of multi-scale numerical simulation of flow,heat and mass transfer behaviors in dense gas-solid flows.It describes multiscale models(direct numerical simulation,discrete particle model,and two-fluid model)and the results of related research.Finally,it discusses possible future developments in research on the flow,heat and mass transfer characteristics of dense gas-solid two-phase flows.
基金supported by the National Natural Science Foundation of China (Grant No. 11172218)the Fundamental Research Funds for the Central Universities (Grant No. 2012206020209)
文摘Slurry jets in a static uniform environment were simulated with a two-phase mixture model in which flow-particle interactions were considered. A standard k-e turbulence model was chosen to close the governing equations. The computational results were in agreement with previous laboratory measurements. The characteristics of the two-phase flow field and the influences of hydraulic and geometric parameters on the distribution of the slurry jets were analyzed on the basis of the computational results. The calculated results reveal that if the initial velocity of the slurry jet is high, the jet spreads less in the radial direction. When the slurry jet is less influenced by the ambient fluid （when the Stokes number St is relatively large）, the turbulent kinetic energy k and turbulent dissipation rate e, which are relatively concentrated around the jet axis, decrease more rapidly after the slurry jet passes through the nozzle. For different values of St, the radial distributions of streamwise velocity and particle volume fraction are both self-similar and fit a Gaussian profile after the slurry jet fully develops. The decay rate of the particle velocity is lower than that of water velocity along the jet axis, and the axial distributions of the centerline particle streamwise velocity are self-similar along the jet axis. The pattern of particle dispersion depends on the Stokes number St. When St = 0.39, the panicle dispersion along the radial direction is considerable, and the relative velocity is very low due to the low dynamic response time. When St = 3.08, the dispersion of particles along the radial direction is very little, and most of the particles have high relative velocities along the streamwise direction.
基金supported by the National Magnetic Confinement Fusion Science Program of China(No.2010GB104005)National Natural Science Foundation of China(No.51406085)
文摘An accurate critical heat flux（CHF） prediction method is the key factor for realizing the steady-state operation of a water-cooled divertor that works under one-sided high heating flux conditions.An improved CHF prediction method based on Euler＇s homogeneous model for flow boiling combined with realizable k-ε model for single-phase flow is adopted in this paper in which time relaxation coefficients are corrected by the Hertz-Knudsen formula in order to improve the calculation accuracy of vapor-liquid conversion efficiency under high heating flux conditions.Moreover,local large differences of liquid physical properties due to the extreme nonuniform heating flux on cooling wall along the circumference direction are revised by formula IAPWSIF97.Therefore,this method can improve the calculation accuracy of heat and mass transfer between liquid phase and vapor phase in a CHF prediction simulation of water-cooled divertors under the one-sided high heating condition.An experimental example is simulated based on the improved and the uncorrected methods.The simulation results,such as temperature,void fraction and heat transfer coefficient,are analyzed to achieve the CHF prediction.The results show that the maximum error of CHF based on the improved method is 23.7%,while that of CHF based on uncorrected method is up to 188%,as compared with the experiment results of Ref.[12].Finally,this method is verified by comparison with the experimental data obtained by International Thermonuclear Experimental Reactor（ITER）,with a maximum error of 6% only.This method provides an efficient tool for the CHF prediction of water-cooled divertors.
