低阶格式在大涡模拟计算中的适用性

采用三维Taylor-Green涡作为研究对象,利用工程中常用的低阶数值格式,研究格式本身的数值误差对大涡模拟计算的影响.结果表明:三种数值格式的数值耗散行为都与亚格子模型行为类似,即在小雷诺数下,流场比较光滑时,耗散很小,当雷诺数增加,流动转捩为湍流,流场梯度增大,耗散显著增大.对于MUSCL格式和二阶有界中心格式,在高雷诺数下,亚格子尺度模型没有明显改善计算结果,但也没有使计算结果恶化.中心格式相比其它两种格式,数值耗散最小,但是在高雷诺数湍流情况下,中心格式的数值耗散仍然主导了能量的耗散,再添加亚格子模型,计算结果反而变得稍差.对于工程中的低阶格式而言,采用中心格式计算大涡模拟是比较好的选择,而且在计算不存在稳定性问题时,采用不添加亚格子模型的隐式大涡模拟效果更好.     

四阶精度显式滤波大涡模拟方法研究

在传统大涡模拟方法的基础上引入显式滤波技术,有效地控制了最小解析尺度区域的各种数值误差.对有限差分格式的守恒特性进行了研究,发展了适用于三维平板通道大涡模拟的四阶精度守恒格式.分别利用二阶和四阶精度的传统与显式滤波大涡模拟程序模拟平板通道内的湍流流动,将计算结果与DNS数据进行了对比,得到了合理结果.     

K(m,n,P)方程多-Compacton相互作用的数值研究

采用有限差分法对非线性色散K(m,n,p)方程的多-Compacton之间的相互作用进行了数值研究.该差分方法为二阶精度且线性意义下绝对稳定的无耗散格式,通过添加人工耗散项有效防止了数值解的爆破现象.首先对单-Compacton的长时间演化行为进行了数值模拟,验证了数值方法的有效性.然后对双-CompaCton和三-Compacton的碰撞过程进行了数值研究,发现多-Compacton碰撞之后基本保持碰撞之前的波形和波速,但在波后产生小振幅的Compacton-Anticompacton对.     

可压涡卷空间演化的迎风紧致差分数值模拟

从数值算法的耗散和色散特征的时空全离散Fourier分析出发,通过直接求解二维非定常可压Navier Stokes方程,将发展的5阶迎风紧致差分格式用于无约束可压平面受迫剪切层中基频涡卷空间演化过程的数值模拟.采用被动守恒标量等方法显示了基频涡卷的饱和、一次对并、二次对并等现象,据此探讨了入口来流亚谐扰动引起的初值效应问题,表明可压大尺度涡结构空间演化形态与受迫扰动方式之间存在关联.     

一种色散介质FDTD通用吸收边界

由单轴各向异性介质所满足的场方程出发,并根据相位匹配原理以得到介质中的色散关系和无反射条件.结合频域到时域的转换关系（即用/t代替jω）和移位算子时域有限差分（SO-FDTD）方法提出了一种适用于各向同性常见色散介质模型,包括德拜模型、洛伦兹模型、德鲁德模型等的通用UPML吸收边界.数值结果表明了该算法的通用性和高效性. 关键词： 时域有限差分 吸收边界 移位算子 色散介质     

针对典型分离流动数值模拟的自适应耗散调节方法

分离流动是一种复杂湍流现象,其中微小尺度结构的发展演化有着重要的影响,但是湍流数值模拟中数值方法的固有耗散会抑制小尺度流动结构的发展.因此,本文结合延迟分离涡模拟方法,基于五阶耗散紧致格式,通过引入流场相关的调节因子,构造自适应耗散调节方法,使其在大涡模拟区域降低耗散影响以提高对小尺度流动结构的辨识能力,在雷诺平均区域恢复正常耗散水平以保持稳定求解.本文首先通过近似色散关系和涡输运算例说明该方法可有效降低数值耗散影响,同时具有更高分辨率的涡保持能力;随后通过各向同性衰减湍流和平板槽道流动展示了该方法可有效提高对微小流动结构的辨识能力,同时在存在大梯度的流场中可保持求解稳定性;最后求解Re=3900亚临界状态下的圆柱绕流,通过对比流场云图、时均速度和雷诺应力分布曲线,说明了该方法通过减小分离区耗散影响,可有效提高对典型分离流动求解的准确性.     

混合通量差分格式的比较

系统研究了几种混合通量差分格式的构造方法和耗散模型,分别对低速平板绕流、二维跨音速喷管流动和高超音速钝头体无粘绕流进行了数值模拟,结合先进的EASM湍流模型对格式的粘性分辨率和激波稳定性进行了细致的比较分析.结果表明混合通量差分格式兼顾了FDS和FVS格式的优点,具有较高的间断分辨率和数值稳定性.     

色散介质电磁特性时域有限差分分析的Newmark方法

根据Debye模型、Drude模型和Lorentz模型3种常见色散介质模型频域极化率的特点,利用频域到时域的转换关系jω→?/?t,将极化矢量P与电场强度E的频域关系转换成时域内关于P的二阶微分方程,其对3种色散介质模型皆适用,具有统一的形式.然后采用相比于中心差分具有更高精度的Newmark两步算法(Newmark-β-γ法)求解该方程,进而得到E→P的递推公式,再结合本构关系得到D→E的时域递推式.实现了色散介质电磁场量的时域有限差分迭代计算.数值计算结果表明该方法是适用于3种色散介质模型的通用算法,并且相比于移位算子时域有限差分方法等以中心差分为基础的离散方案具有更高的计算精度.     

一种处理色散介质问题的通用时域有限差分方法

色散介质的介电系数是频率的函数,使本构关系在时域成为卷积关系.这就给用时域有限差分方法计算色散介质中波的散射和传播带来了困难.现有算法往往要针对不同色散介质模型推导相应的递推公式,算法的通用性较差.本文完善和发展了移位算子-时域有限差分方法,使之成为一种处理色散介质电磁问题的通用方法.首先,证明了常见的三种色散介质模型(德拜模型、洛伦兹模型和德鲁模型)的介电系数均可以写成适于移位算子法计算的有理分式函数形式.然后,用/t代替jω,过渡到时域,再引入时域移位算子z_t代替时间微分算子来处理有理分式函数形式的介电系数,给出离散时域本构关系的表示式,进而导出时域有限差分方法当中电位移矢量和电场强度之间的关系.最后,计算了几种色散介质的电磁散射,数值结果表明了本文方法和程序的通用性和正确有效性. 关键词： 时域有限差分方法 色散介质 移位算子     

非结构网格上大涡模拟数值方法研究

随着计算能力的飞速发展,在非结构网格上进行大涡模拟是研究叶轮机械真实几何结构下复杂流动问题的有效途径。本文针对非结构网格上数值格式人工黏性过大的缺点,利用湍流中动能级联特性,发展了适用于大涡模拟的低耗散数值格式,并结合亚格子模型,考察了其在基础湍流算例中的有效性。通过与传统迎风格式的对比,说明本文发展的低耗散格式可以显著提高非结构网格上大涡模拟的精度和可靠性。     

A class of finite difference schemes with low dispersion and controllable dissipation for DNS of compressible turbulence

In this paper, a class of finite difference schemes which achieves low dispersion and controllable dissipation in smooth region and robust shock-capturing capabilities in the vicinity of discontinuities is presented. Firstly, a sufficient condition for semi-discrete finite difference schemes to have independent dispersion and dissipation is derived. This condition enables a novel approach to separately optimize the dissipation and dispersion properties of finite difference schemes and a class of schemes with minimized dispersion and controllable dissipation is thus obtained. Secondly, for the purpose of shock-capturing, one of these schemes is used as the linear part of the WENO scheme with symmetrical stencils to constructed an improved WENO scheme. At last, the improved WENO scheme is blended with its linear counterpart to form a new hybrid scheme for practical applications. The proposed scheme is accurate, flexible and robust. The accuracy and resolution of the proposed scheme are tested by the solutions of several benchmark test cases. The performance of this scheme is further demonstrated by its application in the direct numerical simulation of compressible turbulent channel flow between isothermal walls.     

A high-resolution,hybrid compact-WENO scheme with minimized dispersion and controllable dissipation

In this paper,a high-resolution,hybrid compact-WENO scheme is developed based on the minimized dispersion and controllable dissipation reconstruction technique.Firstly,a suffcient condition for a family of tri-diagonal compact schemes to have independent dispersion and dissipation is derived.Then,a specific 4th order compact scheme with low dispersion and adjustable dissipation is constructed and analyzed.Finally,the optimized compact scheme is blended with the WENO scheme to form the hybrid scheme.Moreover,the approximation dispersion relation approach is employed to optimize the spectral properties of the nonlinear scheme to yield the true wave propagation behavior of the finite difference scheme.Several test cases are carried out to verify the highresolution as well as the robust shock-capturing capabilities of the proposed scheme.     

Analyses of the Dispersion Overshoot and Inverse Dissipation of the High-Order Finite Difference Scheme

Analyses were performed on the dispersion overshoot and inverse dissipation of the high-order finite difference scheme using Fourier and precision analysis. Schemes under discussion included the pointwise- and staggered-grid type, and were presented in weighted form using candidate schemes with third-order accuracy and three-point stencil. All of these were commonly used in the construction of difference schemes. Criteria for the dispersion overshoot were presented and their critical states were discussed. Two kinds of instabilities were studied due to inverse dissipation, especially those that occur at lower wave numbers. Criteria for the occurrence were presented and the relationship of the two instabilities was discussed. Comparisons were made between the analytical results and the dispersion/dissipation relations by Fourier transformation of typical schemes. As an example, an application of the criteria was given for the remedy of inverse dissipation in Weirs and Martin's third-order scheme.     

Generalized finite compact difference scheme for shock/complex flowfield interaction

A class of generalized high order finite compact difference schemes is proposed for shock/vortex, shock/boundary layer interaction problems. The finite compact difference scheme takes the region between two shocks as a compact stencil. The high order WENO fluxes on shock stencils are used as the internal boundary fluxes for the compact scheme. A lemma based on the property of smoothness estimators on a 5-points stencil is given to detect the shock position. There is no free parameter introduced to switch the compact scheme and the WENO scheme. Some numerical experiments are given and they demonstrate that the present scheme has low dissipation due to the compact central differencing scheme used in the smooth regions.     

分辨率优化的混合WENO格式

为提高有限差分格式的分辨率,利用傅里叶分析对WENO格式进行色散及耗散优化,并给出优化的线性权重.用优化后的WENO格式与保单调格式（MP）进行加权混合,得到新的加权混合WENO格式（H-WENO）.通过一维激波管问题、Shu-Osher问题及二维双Mach反射问题及R-T不稳定性问题对格式进行数值测试.结果显示,新格式具有强健的激波捕捉能力和对小尺度波结构的高分辨率,与原WENO格式相比改进明显.     

Construction of explicit and implicit dynamic finite difference schemes and application to the large-eddy simulation of the Taylor–Green vortex

A general class of explicit and implicit dynamic finite difference schemes for large-eddy simulation is constructed, by combining Taylor series expansions on two different grid resolutions. After calibration for Re→∞, the dynamic finite difference schemes allow to minimize the dispersion errors during the calculation through the real-time adaption of a dynamic coefficient. In case of DNS resolution, these dynamic schemes reduce to Taylor-based finite difference schemes with formal asymptotic order of accuracy, whereas for LES resolution, the schemes adapt to Dispersion-Relation Preserving schemes. Both the explicit and implicit dynamic finite difference schemes are tested for the large-eddy simulation of the Taylor–Green vortex flow and numerical errors are investigated as well as their interaction with the dynamic Smagorinsky model and the multiscale Smagorinsky model. Very good results are obtained.     

On the role of helicity in complex fluid flows

Results of direct numerical simulations of the Taylor-Green vortex are analysed by conditional sampling. In regions of small energy dissipation, there are tendencies for (1) velocity, u, and vorticity, ω, to be aligned and (2) vorticity and curl of vorticity, ? × ω, to be nearly orthogonal. The fields of dissipation, enstrophy, turbulence production, and vortex stretching exhibit a striking similarity in their spatial structure.     

Scale separation for implicit large eddy simulation

With implicit large eddy simulation (ILES) the truncation error of the discretization scheme acts as subgrid-scale (SGS) model for the computation of turbulent flows. Although ILES is comparably simple, numerically robust and easy to implement, a considerable challenge is the design of numerical discretization schemes resulting in a physically consistent SGS model. In this work, we consider the implicit SGS modeling capacity of the adaptive central-upwind weighted-essentially-non-oscillatory scheme (WENO-CU6) [X.Y. Hu, Q. Wang, N.A. Adams, An adaptive central-upwind weighted essentially non-oscillatory scheme, J. Comput. Phys. 229 (2010) 8952–8965] by incorporating a physically-motivated scale-separation formulation. Scale separation is accomplished by a simple modification of the WENO weights. The resulting modified scheme maintains the shock-capturing capabilities of the original WENO-CU6 scheme while it is also able to reproduce the Kolmogorov range of the kinetic-energy spectrum for turbulence at the limit of infinite Reynolds number independently of grid resolution. For isentropic compressible turbulence the pseudo-sound regime of the dilatational kinetic-energy spectrum and the non-Gaussian probability-density function of the longitudinal velocity derivative are reproduced.     

Nonuniform time-step Runge–Kutta discontinuous Galerkin method for Computational Aeroacoustics

With many superior features, Runge–Kutta discontinuous Galerkin method (RKDG), which adopts Discontinuous Galerkin method (DG) for space discretization and Runge–Kutta method (RK) for time integration, has been an attractive alternative to the finite difference based high-order Computational Aeroacoustics (CAA) approaches. However, when it comes to complex physical problems, especially the ones involving irregular geometries, the time step size of an explicit RK scheme is limited by the smallest grid size in the computational domain, demanding a high computational cost for obtaining time accurate numerical solutions in CAA. For computational efficiency, high-order RK method with nonuniform time step sizes on nonuniform meshes is developed in this paper. In order to ensure correct communication of solutions on the interfaces of grids with different time step sizes, the values at intermediate-stages of the Runge–Kutta time integration on the elements neighboring such interfaces are coupled with minimal dissipation and dispersion errors. Based upon the general form of an explicit p-stage RK scheme, a linear coupling procedure is proposed, with details on the coefficient matrices and execution steps at common time-levels and intermediate time-levels. Applications of the coupling procedures to Runge–Kutta schemes frequently used in simulation of fluid flow and acoustics are given, including the third-order TVD scheme, and low-storage low dissipation and low dispersion (LDDRK) schemes. In addition, an analysis on the stability of coupling procedures on a nonuniform grid is carried out. For validation, numerical experiments on one-dimensional and two-dimensional problems are presented to illustrate the stability and accuracy of proposed nonuniform time-step RKDG scheme, as well as the computational benefits it brings. Application to a one-dimensional nonlinear problem is also investigated.