共查询到10条相似文献,搜索用时 31 毫秒
1.
《力学学报》2012,44(4)
采用五阶有限差分WENO格式直接模拟了高初始湍流Mach数的可压缩均匀各向同性湍流,主要分析了湍流的统计特性和压缩性的影响,包括能谱特征、激波串、耗散率、标度律等.研究表明,湍动能主要来自于速度场螺旋分量的贡献;各向同性湍流的小尺度脉动对压缩性更为敏感,并且压缩性的增强加快了湍流大尺度脉动向小尺度脉动的湍动能输运;随着湍流Mach数的升高,胀量(压缩)耗散率所占比率也显著增长.标度律分析表明,强可压缩湍流的横向速度结构函数仍然具有扩展自相似性;当阶数较高(p≥5)时,纵向速度结构函数的扩展自相似性则不再成立.对于压缩性较弱的湍流,与不可压缩湍流一致,横向湍流脉动的间歇性要强于纵向湍流脉动;而对于强可压缩湍流,纵向湍流脉动的间歇性要强于横向湍流脉动. 相似文献
2.
可压缩各向同性衰减湍流直接数值模拟研究 总被引:2,自引:3,他引:2
采用五阶有限差分WENO格式直接模拟了高初始湍流Mach数的可压缩均匀各向同性湍流,主要分析了湍流的统计特性 和压缩性的影响,包括能谱特征、激波串、耗散率、标度律等. 研究表明,湍动能主要来自于速度场螺旋分量的贡献;各向同性湍流的小尺度脉动对压缩性更为敏感,并且压缩性的增强加快了湍流大 尺度脉动向小尺度脉动的湍动能输运;随着湍流Mach数的升高,胀量(压缩)耗散率所占比率也显著增长. 标度律分析表明,强可压缩湍流的横向速度结构函数仍然具有扩展自相似性;当阶数较高(p ≥ 5)时,纵向速度结构函数的扩展自相似性则不再成立. 对于压缩性较弱的湍流,与不可压缩湍流一致,横向湍流脉动的间歇性要强于纵向湍流脉动;而对于强可压缩湍流,纵向湍流脉动的 间歇性要强于横向湍流脉动. 相似文献
3.
Hybrid schemes are very efficient for complex compressible flow simulation. However, for most existing hybrid schemes in literature, empirical problem‐dependent parameters are always needed to detect shock waves and hence greatly decrease the robustness and accuracy of the hybrid scheme. In this paper, based on the nonlinear weights of the weighted essentially non‐oscillatory (WENO) scheme, a novel weighting switch function is proposed. This function approaches 1 with high‐order accuracy in smooth regions and 0 near discontinuities. Then, with the new weighting switch function, a seventh‐order hybrid compact‐reconstruction WENO scheme (HCCS) is developed. The new hybrid scheme uses the same stencil as the fifth‐order WENO scheme, and it has seventh‐order accuracy in smooth regions even at critical points. Numerical tests are presented to demonstrate the accuracy and robustness of both the switch function and HCCS. Comparisons also reveal that HCCS has lower dissipation and less computational cost than the seventh‐order WENO scheme. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
4.
This paper describes a new class of three‐dimensional finite difference schemes for high‐speed turbulent flows in complex geometries based on the high‐order monotonicity‐preserving (MP) method. Simulations conducted for various 1D, 2D, and 3D problems indicate that the new high‐order MP schemes can preserve sharp changes in the flow variables without spurious oscillations and are able to capture the turbulence at the smallest computed scales. Our results also indicate that the MP method has less numerical dissipation and faster grid convergence than the weighted essentially non‐oscillatory method. However, both of these methods are computationally more demanding than the COMP method and are only used for the inviscid fluxes. To reduce the computational cost for reacting flows, the scalar equations are solved by the COMP method, which is shown to yield similar results to those obtained by the MP in supersonic turbulent flows with strong shock waves. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
5.
Jinwei Bai;Meiliang Mao;Yankai Ma;Zhen-Guo Yan;Yaobing Min; 《国际流体数值方法杂志》2024,96(9):1592-1616
Weighted compact nonlinear schemes (WCNSs) are a popular family of high-resolution shock-capturing schemes for simulating compressible flows, of which the nonlinear interpolation procedure is dominant for the performance. In this work, a simplified weighting strategy is introduced for the nonlinear interpolation procedure. Firstly, an equivalent weighting formulation of WCNS is presented by explicitly including the whole-point stencil into the set of candidate stencils. Secondly, motivated by the reorganization of WCNS, the WCNS-CU6 scheme is achieved in a more straightforward way. Thirdly, by introducing a TENO selection procedure in the framework of WCNS-CU6-Simplified, a TCNS6-Simplified scheme is proposed, the resolution of which is comparable with the excellent TENO6 scheme, while the computational cost is much lower. The simplified schemes exhibit more outstanding, at least comparable, fidelity than the original schemes, however, with superior characteristics in terms of efficiency and simplicity. A variety of benchmark test problems are studied to demonstrate the behaviour of the simplified weighting strategy. 相似文献
6.
The purpose of the present paper is to evaluate very‐high‐order upwind schemes for the direct numerical simulation (DNS ) of compressible wall‐turbulence. We study upwind‐biased (UW ) and weighted essentially nonoscillatory (WENO ) schemes of increasingly higher order‐of‐accuracy (J. Comp. Phys. 2000; 160 :405–452), extended up to WENO 17 (AIAA Paper 2009‐1612, 2009). Analysis of the advection–diffusion equation, both as Δx→0 (consistency), and for fixed finite cell‐Reynolds‐number ReΔx (grid‐resolution), indicates that the very‐high‐order upwind schemes have satisfactory resolution in terms of points‐per‐wavelength (PPW ). Computational results for compressible channel flow (Re∈[180, 230]; M?CL ∈[0.35, 1.5]) are examined to assess the influence of the spatial order of accuracy and the computational grid‐resolution on predicted turbulence statistics, by comparison with existing compressible and incompressible DNS databases. Despite the use of baseline O(Δt2) time‐integration and O(Δx2) discretization of the viscous terms, comparative studies of various orders‐of‐accuracy for the convective terms demonstrate that very‐high‐order upwind schemes can reproduce all the DNS details obtained by pseudospectral schemes, on computational grids of only slightly higher density. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
7.
Jeffrey R. Chasnov 《Applied Scientific Research》1993,51(1-2):313-317
We consider the generation of passive scalar fluctuations by decaying isotropic turbulence in the presence of a uniform mean scalar gradient. At high Reynolds numbers, two distinct similarity states may be established depending on the form of the energy spectrum at low wavenumber magnitude (k). In the first similarity state characterized by a low wavenumber magnitude energy spectrum proportional tok
2, the mean-square scalar fluctuation grows liket
4/5, while in the second similarity state characterized by a spectrum proportional tok
4, the mean-square scalar fluctuation grows approximately liket
4/7. These two high Reynolds number asymptotic similarity states have been subsequently confirmed by large-eddy numerical simulations. As a consequence of the decreasing flow Reynolds number as the turbulence decays, these similarity states do not continue indefinitely. At very long times, a final period of decay of the turbulence occurs, and in this final period, the mean-square scalar fluctuation in the first state continues to grow liket
1/2, while that in the second state ultimately decays liket
–1/2. 相似文献
8.
9.
This work describes the implementation and analysis of high‐order accurate schemes applied to high‐speed flows on unstructured grids. The class of essentially non‐oscillatory schemes (ENO), that includes weighted ENO schemes (WENO), is discussed in the paper with regard to the implementation of third‐ and fourth‐order accurate methods. The entire reconstruction process of ENO and WENO schemes is described with emphasis on the stencil selection algorithms. The stencils can be composed by control volumes with any number of edges, e.g. triangles, quadrilaterals and hybrid meshes. In the paper, ENO and WENO schemes are implemented for the solution of the dimensionless, 2‐D Euler equations in a cell centred finite volume context. High‐order flux integration is achieved using Gaussian quadratures. An approximate Riemann solver is used to evaluate the fluxes on the interfaces of the control volumes and a TVD Runge–Kutta scheme provides the time integration of the equations. Such a coupling of all these numerical tools, together with the high‐order interpolation of primitive variables provided by ENO and WENO schemes, leads to the desired order of accuracy expected in the solutions. An adaptive mesh refinement technique provides better resolution in regions with strong flowfield gradients. Results for high‐speed flow simulations are presented with the objective of assessing the implemented capability. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
10.
Qian Jian 《Acta Mechanica Sinica》1994,10(1):12-15
Using nonequilibrium statistical mechanics closure method, it is shown that the skewness factor of the velocity derivative
of isotropic turbulence approaches a constant −0.515 when the Reynolds number is very high, which is in agreement with the
DNS (direct numerical simulation) result of Vincent and Meneguzzi (1991).
The project supported by the National Basic Research Program “Non-linear Science” 相似文献