The internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence mod...The internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence model and pressure Poisson equation were discretized by upwind difference scheme.A new full implicit difference scheme of 5-point was constructed by using finite volume method and finite difference method.A large sparse matrix with five diagonals was formed and was stored by three arrays of one dimension in a compressed mode.General iterative methods do not work wel1 with large sparse matrix.With algebraic multigrid method(AMG),linear algebraic system of equations was solved and the precision was set at 10-6.The computation results were compared with the experimental results.The results show that the computation results have a good agreement with the experiment data.The precision of computational results and numerical simulation efficiency are greatly improved.展开更多
The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete elem...The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
Mine safety have top-five disasters,which including the water,gas,fire,dust and geological dynamic disaster.The coal mine water disaster is one of the important factors which restricted the development of China’s coa...Mine safety have top-five disasters,which including the water,gas,fire,dust and geological dynamic disaster.The coal mine water disaster is one of the important factors which restricted the development of China’s coal production.It is showed by statistics that 60%of mine accidents are affected by groundwater,which not only result in the production losses,casualties and a variety of展开更多
The main purpose of this research is the second-order modeling of flow and turbulent heat flux in nonpremixed methane-air combustion.A turbulent stream of non-premixed combustion in a stoichiometric condition,is numer...The main purpose of this research is the second-order modeling of flow and turbulent heat flux in nonpremixed methane-air combustion.A turbulent stream of non-premixed combustion in a stoichiometric condition,is numerically analyzed through the Reynolds averaged Navier-Stokes(RANS) equations.For modeling radiation and combustion,the discrete ordinates(DO) and eddy dissipation concept model have been applied.The Reynolds stress transport model(RSM) also was used for turbulence modeling.For THF in the energy equation,the GGDH model and high order algebraic model of HOGGDH with simple eddy diffusivity model have been applied.Comparing the numerical results of the SED model(with the turbulent Prandtl 0.85) and the second-order heat flux models with available experimental data follows that applying the second-order models significantly led to the modification of predicting temperature distribution and species mass fraction distribution in the combustion chamber.Calculation of turbulent Prandtl number in the combustion chamber shows that the assumption of Pr_(t) of 0.85 is far from reality and Pr_(t) in different areas varies from 0.4 to 1.2.展开更多
This work provides a method to predict the three-dimensional equivalent elastic properties of the filament-wound composites based on the multi-scale homogenization principle.In the meso-scale,a representative volume e...This work provides a method to predict the three-dimensional equivalent elastic properties of the filament-wound composites based on the multi-scale homogenization principle.In the meso-scale,a representative volume element(RVE)is defined and the bridging model is adopted to establish a theoretical predictive model for its three-dimensional equivalent elastic constants.The results obtained through this method for the previous experimental model are compared with the ones gained respectively by experiments and classical laminate theory to verify the reliability of this model.In addition,the effects of some winding parameters,such as winding angle,on the equivalent elastic behavior of the filament-wound composites are analyzed.The rules gained can provide a theoretical reference for the optimum design of filament-wound composites.展开更多
In order to calculate the air temperature of the near surface layer in urban environment,the surface layer air was divided into several sections in the vertical direction,and some energy balance equations were develop...In order to calculate the air temperature of the near surface layer in urban environment,the surface layer air was divided into several sections in the vertical direction,and some energy balance equations were developed for each air layer,in which the heat exchange due to vertical turbulence and horizontal air flow was taken into account.Then,the vertical temperature distribution of the surface layer air was obtained through the coupled calculation using the energy balance equations of underlying surfaces and building walls.Moreover,the measured air temperatures in a small area(with a horizontal scale of less than 500 m) and a large area(with a horizontal scale of more than 1 000 m) in Guangzhou in summer were used to validate the proposed model.The calculated results accord well with the measured ones,with a maximum relative error of 4.18%.It is thus concluded that the proposed model is a high-accuracy method to theoretically analyze the urban heat island and the thermal environment.展开更多
Calculation grid and turbulence model for numerical simulating pressure fluctuations in a high-speed train tunnel are studied through the comparison analysis of numerical simulation and moving model test.Compared the ...Calculation grid and turbulence model for numerical simulating pressure fluctuations in a high-speed train tunnel are studied through the comparison analysis of numerical simulation and moving model test.Compared the waveforms and peak-peak values of pressure fluctuations between numerical simulation and moving model test,the structured grid and the SST k-ωturbulence model are selected for numerical simulating the process of high-speed train passing through the tunnel.The largest value of pressure wave amplitudes of numerical simulation and moving model test meet each other.And the locations of the largest value of the initial compression and expansion wave amplitude of numerical simulation are in agreement with that of moving model test.The calculated pressure at the measurement point fully conforms to the propagation law of compression and expansion waves in the tunnel.展开更多
A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-indu...A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.展开更多
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.展开更多
The influence of rheological parameters on vortex dynamics and the extent of drag reduction (DR) were deciphered via extensively analyzing the hi-fidelity direct numerical simulation results of the turbulent channel f...The influence of rheological parameters on vortex dynamics and the extent of drag reduction (DR) were deciphered via extensively analyzing the hi-fidelity direct numerical simulation results of the turbulent channel flow with polymer solutions. It has been observed that in all drag reduction regimes from the onset of DR to maximum drag reduction (MDR) limit, the Deborah number is defined as the product of an effective Weissenberg number, and the root mean square streamwise vorticity fluctuation remains O(1) in the near wall region. The ratio of the average lifetime of axial vortices to the vortex rotating duration decreases with increasing DR, and MDR is achieved when these time scales become nearly equal. Based on these observations a simple framework is proposed adequately to describe the influence of polymer additives on the extent of DR from onset to MDR as well as the universality of the MDR in flow systems with polymer additives.展开更多
In order to calculate the air temperature of the near surface layer in urban environment,the surface layer air was divided into several layers in the vertical direction,and some energy balance equations were developed...In order to calculate the air temperature of the near surface layer in urban environment,the surface layer air was divided into several layers in the vertical direction,and some energy balance equations were developed for each air layer,in which the heat exchange due to vertical turbulence and horizontal air flow was taken into account.Then,the vertical temperature distribution of the surface layer air was obtained through the coupled calculation using the energy balance equations of underlying surfaces and building walls.Moreover,the measured air temperatures in a small area(with a horizontal scale of less than 500 m)and a large area(with a horizontal scale of more than 1 000 m)in Guangzhou in summer were used to validate the proposed model.The calculated results agree well with the measured ones,with a maximum relative error of 4.18%.It is thus concluded that the proposed model is a high-accuracy method to theoretically analyze the urban heat island and the thermal environment.展开更多
弧形闸门是水利工程中重要的挡水、泄水建筑物,虽然弧形闸门被设计成具有足够的刚度来承受设计水压,但泄水过程中水压的脉动可能会使闸门产生较大的振动,进而造成闸门的破坏。结合湍流模型和流体体积方法(volume of fluid,VOF)建立弧形...弧形闸门是水利工程中重要的挡水、泄水建筑物,虽然弧形闸门被设计成具有足够的刚度来承受设计水压,但泄水过程中水压的脉动可能会使闸门产生较大的振动,进而造成闸门的破坏。结合湍流模型和流体体积方法(volume of fluid,VOF)建立弧形闸门泄流三维湍流流场数值模型,数值求解采用稳态计算和瞬态计算两个连续的运算步骤进行,以便更好地确定入口处的流速。分别采用k-ε湍流模型和k-ω湍流模型对闸门周围流场和作用在弧形面板上的流体压力进行计算,结果表明:k-ε湍流模型与壁函数相结合的方法不能捕捉到稳定泄流阶段闸门面板上的压力脉动行为,而k-ω湍流模型结合壁面积分不仅能够得到闸门周围的流场变化,而且能准确计算闸门面板上的脉动压力。基于k-ω湍流模型计算结果,分析下游水位变化对闸门周围流场和脉动压力的影响,闸门面板上压力的脉动主要是由于闸门前的漩涡引起的,压力脉动的优势频率取决于闸孔出流形式,与上下游水位差无关,自由出流时的压力脉动优势频率比淹没出流时的大。展开更多
基金Projects(59375211,10771178,10676031) supported by the National Natural Science Foundation of ChinaProject(07A068) supported by the Key Project of Hunan Education CommissionProject(2005CB321702) supported by the National Key Basic Research Program of China
文摘The internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence model and pressure Poisson equation were discretized by upwind difference scheme.A new full implicit difference scheme of 5-point was constructed by using finite volume method and finite difference method.A large sparse matrix with five diagonals was formed and was stored by three arrays of one dimension in a compressed mode.General iterative methods do not work wel1 with large sparse matrix.With algebraic multigrid method(AMG),linear algebraic system of equations was solved and the precision was set at 10-6.The computation results were compared with the experimental results.The results show that the computation results have a good agreement with the experiment data.The precision of computational results and numerical simulation efficiency are greatly improved.
基金Project(51378006) supported by National Natural Science Foundation of ChinaProject(141076) supported by Huoyingdong Foundation of the Ministry of Education of China+1 种基金Project(2242015R30027) supported by Excellent Young Teacher Program of Southeast University,ChinaProject(BK20140109) supported by the Natural Science Foundation of Jiangsu Province,China
文摘The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
文摘Mine safety have top-five disasters,which including the water,gas,fire,dust and geological dynamic disaster.The coal mine water disaster is one of the important factors which restricted the development of China’s coal production.It is showed by statistics that 60%of mine accidents are affected by groundwater,which not only result in the production losses,casualties and a variety of
文摘The main purpose of this research is the second-order modeling of flow and turbulent heat flux in nonpremixed methane-air combustion.A turbulent stream of non-premixed combustion in a stoichiometric condition,is numerically analyzed through the Reynolds averaged Navier-Stokes(RANS) equations.For modeling radiation and combustion,the discrete ordinates(DO) and eddy dissipation concept model have been applied.The Reynolds stress transport model(RSM) also was used for turbulence modeling.For THF in the energy equation,the GGDH model and high order algebraic model of HOGGDH with simple eddy diffusivity model have been applied.Comparing the numerical results of the SED model(with the turbulent Prandtl 0.85) and the second-order heat flux models with available experimental data follows that applying the second-order models significantly led to the modification of predicting temperature distribution and species mass fraction distribution in the combustion chamber.Calculation of turbulent Prandtl number in the combustion chamber shows that the assumption of Pr_(t) of 0.85 is far from reality and Pr_(t) in different areas varies from 0.4 to 1.2.
文摘This work provides a method to predict the three-dimensional equivalent elastic properties of the filament-wound composites based on the multi-scale homogenization principle.In the meso-scale,a representative volume element(RVE)is defined and the bridging model is adopted to establish a theoretical predictive model for its three-dimensional equivalent elastic constants.The results obtained through this method for the previous experimental model are compared with the ones gained respectively by experiments and classical laminate theory to verify the reliability of this model.In addition,the effects of some winding parameters,such as winding angle,on the equivalent elastic behavior of the filament-wound composites are analyzed.The rules gained can provide a theoretical reference for the optimum design of filament-wound composites.
基金Project(50808083) supported by the National Natural Science Foundation of China
文摘In order to calculate the air temperature of the near surface layer in urban environment,the surface layer air was divided into several sections in the vertical direction,and some energy balance equations were developed for each air layer,in which the heat exchange due to vertical turbulence and horizontal air flow was taken into account.Then,the vertical temperature distribution of the surface layer air was obtained through the coupled calculation using the energy balance equations of underlying surfaces and building walls.Moreover,the measured air temperatures in a small area(with a horizontal scale of less than 500 m) and a large area(with a horizontal scale of more than 1 000 m) in Guangzhou in summer were used to validate the proposed model.The calculated results accord well with the measured ones,with a maximum relative error of 4.18%.It is thus concluded that the proposed model is a high-accuracy method to theoretically analyze the urban heat island and the thermal environment.
文摘Calculation grid and turbulence model for numerical simulating pressure fluctuations in a high-speed train tunnel are studied through the comparison analysis of numerical simulation and moving model test.Compared the waveforms and peak-peak values of pressure fluctuations between numerical simulation and moving model test,the structured grid and the SST k-ωturbulence model are selected for numerical simulating the process of high-speed train passing through the tunnel.The largest value of pressure wave amplitudes of numerical simulation and moving model test meet each other.And the locations of the largest value of the initial compression and expansion wave amplitude of numerical simulation are in agreement with that of moving model test.The calculated pressure at the measurement point fully conforms to the propagation law of compression and expansion waves in the tunnel.
文摘A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.
基金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.
基金Project (10672069) supported by the National Natural Science Foundation of China
文摘The influence of rheological parameters on vortex dynamics and the extent of drag reduction (DR) were deciphered via extensively analyzing the hi-fidelity direct numerical simulation results of the turbulent channel flow with polymer solutions. It has been observed that in all drag reduction regimes from the onset of DR to maximum drag reduction (MDR) limit, the Deborah number is defined as the product of an effective Weissenberg number, and the root mean square streamwise vorticity fluctuation remains O(1) in the near wall region. The ratio of the average lifetime of axial vortices to the vortex rotating duration decreases with increasing DR, and MDR is achieved when these time scales become nearly equal. Based on these observations a simple framework is proposed adequately to describe the influence of polymer additives on the extent of DR from onset to MDR as well as the universality of the MDR in flow systems with polymer additives.
基金Supported by National Natural Science Foundation of China(50538040,50720165805,50808083)the 111 project(111-2-13)State Key Laboratory of Subtropical Building(2008ZB14))
文摘In order to calculate the air temperature of the near surface layer in urban environment,the surface layer air was divided into several layers in the vertical direction,and some energy balance equations were developed for each air layer,in which the heat exchange due to vertical turbulence and horizontal air flow was taken into account.Then,the vertical temperature distribution of the surface layer air was obtained through the coupled calculation using the energy balance equations of underlying surfaces and building walls.Moreover,the measured air temperatures in a small area(with a horizontal scale of less than 500 m)and a large area(with a horizontal scale of more than 1 000 m)in Guangzhou in summer were used to validate the proposed model.The calculated results agree well with the measured ones,with a maximum relative error of 4.18%.It is thus concluded that the proposed model is a high-accuracy method to theoretically analyze the urban heat island and the thermal environment.
文摘弧形闸门是水利工程中重要的挡水、泄水建筑物,虽然弧形闸门被设计成具有足够的刚度来承受设计水压,但泄水过程中水压的脉动可能会使闸门产生较大的振动,进而造成闸门的破坏。结合湍流模型和流体体积方法(volume of fluid,VOF)建立弧形闸门泄流三维湍流流场数值模型,数值求解采用稳态计算和瞬态计算两个连续的运算步骤进行,以便更好地确定入口处的流速。分别采用k-ε湍流模型和k-ω湍流模型对闸门周围流场和作用在弧形面板上的流体压力进行计算,结果表明:k-ε湍流模型与壁函数相结合的方法不能捕捉到稳定泄流阶段闸门面板上的压力脉动行为,而k-ω湍流模型结合壁面积分不仅能够得到闸门周围的流场变化,而且能准确计算闸门面板上的脉动压力。基于k-ω湍流模型计算结果,分析下游水位变化对闸门周围流场和脉动压力的影响,闸门面板上压力的脉动主要是由于闸门前的漩涡引起的,压力脉动的优势频率取决于闸孔出流形式,与上下游水位差无关,自由出流时的压力脉动优势频率比淹没出流时的大。
