The high order compact d if ference method is developed for solving the perturbation equations based on Navi er Stokes equations, and is used in studying complex evolution processes from w all negative pulse to the ...The high order compact d if ference method is developed for solving the perturbation equations based on Navi er Stokes equations, and is used in studying complex evolution processes from w all negative pulse to the turbulent coherent structure in the channel flow. Th is method contains three dimensional coupling difference scheme with high accur acy and high resolution, and the high order time splitting methods. Compared with the general spectral method, the method can be used to research turbule nt coherent structure under more general boundary conditions and in flow domains . In this paper, the generation and evolution of the turbulent coherent structur es ind uced by wall pulse in the channel flow are simulated, and the basic characterist ics and rules of the turbulent coherent structure are shown. Computational r esults indicate that a wall negative pulse is more convenient than the resonant three wave model.展开更多
Making use of Weierstrass＇s theorem and Chebyshev＇s theorem and referring to the equations of state of the scaled-particle theory and the Pereus-Yevick integration equation, we demonstrate that there exists a sequen...Making use of Weierstrass＇s theorem and Chebyshev＇s theorem and referring to the equations of state of the scaled-particle theory and the Pereus-Yevick integration equation, we demonstrate that there exists a sequence of polynomials such that the equation of state is given by the limit of the sequence of polynomials. The polynomials of the best approximation from the third order up to the eighth order are obtained so that the Carnahan-Starling equation can be improved successively. The resulting equations of state are in good agreement with the simulation results on the stable fluid branch and on the metastable fluid branch.展开更多
An attempt was made to extend mild combustion to forward flow furnace, such as the refinery and petrochemical tube furnace. Three dimensional numerical simulation was carried out to study the performance of this furna...An attempt was made to extend mild combustion to forward flow furnace, such as the refinery and petrochemical tube furnace. Three dimensional numerical simulation was carried out to study the performance of this furnace. The Eddy Dissipation Concept(EDC) model coupled with the reaction mechanism DRM-19 was used. The prediction showed a good agreement with the measurement. The effect of air nozzle circle(D), air nozzle diameter(d), air nozzle number(N), and air preheating temperature(Tair) on the flow, temperature and species fields, and the CO and NO emissions was investigated. The results indicate that there are four zones in the furnace, viz.: a central jet zone, an ignition zone, a combustion reaction zone, and a flue gas zone, according to the distribution profiles of H_2 CO and OH. The central jet entrains more flue gas in the furnace upstream with an increasing D while the effect of D is negligible in the downstream. The air jet momentum increases with a decreasing d or an increasing Tair, and entrains more flue gas. The effect of N is mainly identified near the burner exit. More heat is absorbed in the radiant section and less heat is discharged to the atmosphere with a decreasing d and an increasing N as evidenced by the flue gas temperature. The CO and NO emissions are less than 50 μL/L and 10 μL/L, respectively, in most of conditions.展开更多
The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain ...The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain decomposition method for two dimensional time dependent convection diffusion equations with variable coefficients. By combining predictor-corrector technique,modified upwind differences with explicitimplicit coupling,the method under consideration provides intrinsic parallelism while maintaining good stability and accuracy.Moreover,for multi-dimensional problems, the method can be readily implemented on a multi-processor system and does not have the limitation on the choice of subdomains required by some other similar predictor-corrector or stabilized schemes.These properties of the method are demonstrated in this work through both rigorous mathematical analysis and numerical experiments.展开更多
Efficiency in solving the Saint-Venant equations for watershed rainfall-runoff routing is important in flood hydrology. This paper presents a high-efficiency numerical solution of one-dimensional dynamic wave equation...Efficiency in solving the Saint-Venant equations for watershed rainfall-runoff routing is important in flood hydrology. This paper presents a high-efficiency numerical solution of one-dimensional dynamic wave equations(HEDWE) for watershed rainfall-runoff routing, in which the full momentum equation is written as a quadratic equation with only one unknown variable Q, water depth is derived from the continuity equation using the two-step predictorcorrector method, and the discrete scheme is the explicit upwind scheme. The results of numerical tests showed the HEDWE approach has several major advantages. 1) It is a stable numerical method, even for an initially dry area. 2) Its computational efficiency is higher than 4.76E+05 times/s. 3) It can be used for overland flow, river flow, and combinations thereof. The primary disadvantages of the HEDWE approach are its unsuitability for rapidly varying flow, such as dam-break floods.展开更多
The flow over broad-crested weirs was simulated by computational fluid dynamic model. The water surface profile over broad crested weir was measured in a laboratory model and validated using two and three dimensional ...The flow over broad-crested weirs was simulated by computational fluid dynamic model. The water surface profile over broad crested weir was measured in a laboratory model and validated using two and three dimensional Fluent programs. The Reynolds Averaged Navier-Stokes equations coupled with the turbulent standard （k-ε） model and volume of fluid method were applied to estimate the water surface profile. The results of numerical model were compared with experimental results to evaluate the ability of model in describing the behaviour of water surface profile over the weir. The results indicated that the 3D required more time in comparison with 2D results and the flow over weir changed from subcritical flow at the upstream （U/S） face of weir to critical flow over the crest and to supercritical flow at downstream （D/S）. A reasonable agreement was noticed between numerical results and experimental observations with mean error less than 2 %.展开更多
文摘The high order compact d if ference method is developed for solving the perturbation equations based on Navi er Stokes equations, and is used in studying complex evolution processes from w all negative pulse to the turbulent coherent structure in the channel flow. Th is method contains three dimensional coupling difference scheme with high accur acy and high resolution, and the high order time splitting methods. Compared with the general spectral method, the method can be used to research turbule nt coherent structure under more general boundary conditions and in flow domains . In this paper, the generation and evolution of the turbulent coherent structur es ind uced by wall pulse in the channel flow are simulated, and the basic characterist ics and rules of the turbulent coherent structure are shown. Computational r esults indicate that a wall negative pulse is more convenient than the resonant three wave model.
文摘Making use of Weierstrass＇s theorem and Chebyshev＇s theorem and referring to the equations of state of the scaled-particle theory and the Pereus-Yevick integration equation, we demonstrate that there exists a sequence of polynomials such that the equation of state is given by the limit of the sequence of polynomials. The polynomials of the best approximation from the third order up to the eighth order are obtained so that the Carnahan-Starling equation can be improved successively. The resulting equations of state are in good agreement with the simulation results on the stable fluid branch and on the metastable fluid branch.
基金supported by the technology development fund of China Petroleum & Chemical Corporation (Sinopec 312016 and 314054)
文摘An attempt was made to extend mild combustion to forward flow furnace, such as the refinery and petrochemical tube furnace. Three dimensional numerical simulation was carried out to study the performance of this furnace. The Eddy Dissipation Concept(EDC) model coupled with the reaction mechanism DRM-19 was used. The prediction showed a good agreement with the measurement. The effect of air nozzle circle(D), air nozzle diameter(d), air nozzle number(N), and air preheating temperature(Tair) on the flow, temperature and species fields, and the CO and NO emissions was investigated. The results indicate that there are four zones in the furnace, viz.: a central jet zone, an ignition zone, a combustion reaction zone, and a flue gas zone, according to the distribution profiles of H_2 CO and OH. The central jet entrains more flue gas in the furnace upstream with an increasing D while the effect of D is negligible in the downstream. The air jet momentum increases with a decreasing d or an increasing Tair, and entrains more flue gas. The effect of N is mainly identified near the burner exit. More heat is absorbed in the radiant section and less heat is discharged to the atmosphere with a decreasing d and an increasing N as evidenced by the flue gas temperature. The CO and NO emissions are less than 50 μL/L and 10 μL/L, respectively, in most of conditions.
基金the National Natural Science Foundation of China(No.10571017)supported in part by the National Natural Science Foundation of China(No.60533020)supported in part by NSF DMS 0712744
文摘The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain decomposition method for two dimensional time dependent convection diffusion equations with variable coefficients. By combining predictor-corrector technique,modified upwind differences with explicitimplicit coupling,the method under consideration provides intrinsic parallelism while maintaining good stability and accuracy.Moreover,for multi-dimensional problems, the method can be readily implemented on a multi-processor system and does not have the limitation on the choice of subdomains required by some other similar predictor-corrector or stabilized schemes.These properties of the method are demonstrated in this work through both rigorous mathematical analysis and numerical experiments.
基金funded by the National Natural Science Foundation of China (Grant No. 41501046)the Innovation Program of Guangdong Province, China (Grant No. 2016-14)
文摘Efficiency in solving the Saint-Venant equations for watershed rainfall-runoff routing is important in flood hydrology. This paper presents a high-efficiency numerical solution of one-dimensional dynamic wave equations(HEDWE) for watershed rainfall-runoff routing, in which the full momentum equation is written as a quadratic equation with only one unknown variable Q, water depth is derived from the continuity equation using the two-step predictorcorrector method, and the discrete scheme is the explicit upwind scheme. The results of numerical tests showed the HEDWE approach has several major advantages. 1) It is a stable numerical method, even for an initially dry area. 2) Its computational efficiency is higher than 4.76E+05 times/s. 3) It can be used for overland flow, river flow, and combinations thereof. The primary disadvantages of the HEDWE approach are its unsuitability for rapidly varying flow, such as dam-break floods.
文摘The flow over broad-crested weirs was simulated by computational fluid dynamic model. The water surface profile over broad crested weir was measured in a laboratory model and validated using two and three dimensional Fluent programs. The Reynolds Averaged Navier-Stokes equations coupled with the turbulent standard （k-ε） model and volume of fluid method were applied to estimate the water surface profile. The results of numerical model were compared with experimental results to evaluate the ability of model in describing the behaviour of water surface profile over the weir. The results indicated that the 3D required more time in comparison with 2D results and the flow over weir changed from subcritical flow at the upstream （U/S） face of weir to critical flow over the crest and to supercritical flow at downstream （D/S）. A reasonable agreement was noticed between numerical results and experimental observations with mean error less than 2 %.
