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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
为探究双组元空间发动机内自燃推进剂喷雾、燃烧和传热特性,首先对现有MMH/N_(2)O_(4)(甲基肼/四氧化二氮)详细燃烧反应机理进行适当简化,提出适用于三维模拟的MMH/N_(2)O_(4)化学反应动力学模型,并在不同工况下对简化机理的准确性进行...为探究双组元空间发动机内自燃推进剂喷雾、燃烧和传热特性,首先对现有MMH/N_(2)O_(4)(甲基肼/四氧化二氮)详细燃烧反应机理进行适当简化,提出适用于三维模拟的MMH/N_(2)O_(4)化学反应动力学模型,并在不同工况下对简化机理的准确性进行了验证;随后,采用流体体积(volume of fluid,VOF)模型对空间发动机内MMH/N_(2)O_(4)射流撞击雾化过程进行了非稳态模拟,并重点分析了直流互击作用下扇形液膜的形成和演变规律;基于以上化学反应机理及液雾分布模拟结果,在欧拉-拉格朗日体系下构建了离散液滴初始分布,并结合部分搅拌反应器湍流燃烧模型,开展了空间发动机内MMH/N_(2)O_(4)喷雾燃烧及流固耦合传热过程的模拟研究.结果表明,空间发动机内推进剂湍流喷射雾化燃烧过程对壁面冷却液膜的形成和发展具有重要影响,冷却液膜主要出现在燃烧室直线段,在高温燃气与固体域之间建立了一个明显的温度缓冲层,从而实现对发动机壁温的保护.展开更多
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.展开更多
弧形闸门是水利工程中重要的挡水、泄水建筑物,虽然弧形闸门被设计成具有足够的刚度来承受设计水压,但泄水过程中水压的脉动可能会使闸门产生较大的振动,进而造成闸门的破坏。结合湍流模型和流体体积方法(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.
文摘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.
基金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.
基金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.
文摘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.
文摘为探究双组元空间发动机内自燃推进剂喷雾、燃烧和传热特性,首先对现有MMH/N_(2)O_(4)(甲基肼/四氧化二氮)详细燃烧反应机理进行适当简化,提出适用于三维模拟的MMH/N_(2)O_(4)化学反应动力学模型,并在不同工况下对简化机理的准确性进行了验证;随后,采用流体体积(volume of fluid,VOF)模型对空间发动机内MMH/N_(2)O_(4)射流撞击雾化过程进行了非稳态模拟,并重点分析了直流互击作用下扇形液膜的形成和演变规律;基于以上化学反应机理及液雾分布模拟结果,在欧拉-拉格朗日体系下构建了离散液滴初始分布,并结合部分搅拌反应器湍流燃烧模型,开展了空间发动机内MMH/N_(2)O_(4)喷雾燃烧及流固耦合传热过程的模拟研究.结果表明,空间发动机内推进剂湍流喷射雾化燃烧过程对壁面冷却液膜的形成和发展具有重要影响,冷却液膜主要出现在燃烧室直线段,在高温燃气与固体域之间建立了一个明显的温度缓冲层,从而实现对发动机壁温的保护.
基金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.
文摘弧形闸门是水利工程中重要的挡水、泄水建筑物,虽然弧形闸门被设计成具有足够的刚度来承受设计水压,但泄水过程中水压的脉动可能会使闸门产生较大的振动,进而造成闸门的破坏。结合湍流模型和流体体积方法(volume of fluid,VOF)建立弧形闸门泄流三维湍流流场数值模型,数值求解采用稳态计算和瞬态计算两个连续的运算步骤进行,以便更好地确定入口处的流速。分别采用k-ε湍流模型和k-ω湍流模型对闸门周围流场和作用在弧形面板上的流体压力进行计算,结果表明:k-ε湍流模型与壁函数相结合的方法不能捕捉到稳定泄流阶段闸门面板上的压力脉动行为,而k-ω湍流模型结合壁面积分不仅能够得到闸门周围的流场变化,而且能准确计算闸门面板上的脉动压力。基于k-ω湍流模型计算结果,分析下游水位变化对闸门周围流场和脉动压力的影响,闸门面板上压力的脉动主要是由于闸门前的漩涡引起的,压力脉动的优势频率取决于闸孔出流形式,与上下游水位差无关,自由出流时的压力脉动优势频率比淹没出流时的大。
