三维定常问题的高阶紧致差分格式向非定常问题的推广

将已经建立的求解三维定常对流扩散方程的高阶紧致差分格式直接推广到三维非定常对流扩散方程的数值求解,时间导数项利用二阶向后欧拉差分公式,所得到的高阶隐式紧致差分格式时间为二阶精度,空间为四阶精度,并且是无条件稳定的.数值实验结果验证了本文方法的精确性和稳健性.     

求解可压Navier-Stokes方程的对称超紧致差分格式及其并行算法

考查了超紧致差分方法,并将其精度同传统差分格式和紧致差分格式做了比较,结果显示超紧致方法具有良好求解效率.用分块流水线方法设计了超紧致差分格式的并行算法,进行数值实验及并行性能分析.     

非定常对流扩散方程的高精度多重网格方法

由已有的求解定常对流扩散方程的高阶紧致差分格式出发,直接推导出了数值求解非定常对流扩散方程的一种高阶隐式紧致差分格式,其时间为二阶精度,空间为四阶精度,并且是无条件稳定的。为了加快传统迭代法在求解隐格式时在每一个时间步上的迭代收敛速度,采用了多重网格加速技术。数值实验结果验证了本文方法的高阶精度、高效性及高稳定性。     

求解对流扩散方程的紧致修正方法

提出了求解对流扩散方程的紧致修正方法,该方法是在低阶离散格式的源项中,引入紧致修正项,从而构造高阶紧致修正格式,并进行求解.采用紧致修正方法对典型的对流扩散方程进行计算.结果表明,紧致修正方法虽然与二阶经典差分方法建立在相同的结点数上,但紧致修正方法的精度与紧致方法的精度相同,均具有四阶精度.所以紧致修正方法可以在少网...     

气相爆轰高阶中心差分-WENO组合格式自适应网格方法

研究一种高阶中心差分-WENO组合格式,并采用自适应网格方法进行二维和三维气相爆轰波的数值模拟.采用ZND爆轰模型的控制方程为包含化学反应源项的Euler方程组.组合格式在大梯度区采用WENO格式捕捉间断,在光滑区采用高阶中心差分格式提高计算效率.采用一种基于流场结构特征的自适应网格.计算结果,表明这种方法同时具有高精度、高分辨率和高效率的特点.     

一种混合优化差分格式

优化差分格式一般用于计算气动声学和小尺度的湍流数值模拟,这类格式为了获取更好的短波分辨率通常牺牲了部分收敛精度.文章尝试结合最高阶精度格式与优化格式的特点,构造混合优化格式,提高优化格式的收敛精度以及谱分辨率.混合优化格式由模板上的最高阶精度格式与优化格式加权组合得到,权系数由当前模板上的值确定,这使得该格式为非线性格式.对于单色波问题,通过优化权的设计可大幅度减小相位误差.但是加权混合过程使得计算时间有所增加.数值计算证明了该格式的特点.      

耗散紧致格式的频谱特性研究与应用

采用Fourier分析方法,给出线性耗散紧致格式和基于m级Runge-Kutta时间积分方法的全离散格式的频谱特性,并应用五阶耗散紧致格式模拟典型高频波传播和超声速平面Couette流动的特征值问题及其稳定性边值问题,展示耗散紧致差分格式良好的频谱特性.     

基于紧致差分格式的高效时域有限差分算法

探讨一种基于紧致差分格式的高效时域有限差分算法(high-order compact-FDTD),该方法不仅提高计算精度,而且网格结点少、内存使用率和CPU时间大为降低.利用紧致格式FDTD方法实现无耗波导系统及光子晶体光纤中电磁波传播的数值模拟.通过计算实例验证算法的高效性.     

三维对流扩散方程非等距网格上的四阶紧致格式及其多重网格方法

提出了数值求解三维变系数对流扩散方程非等距网格上的四阶精度19点紧致差分格式,为了提高求解效率,采用多重网格方法求解高精度格式所形成的大型代数方程组。数值实验结果表明本文方法对于不同的网格雷诺数问题,在精确性、稳定性和减少计算工作量方面均明显优于7点中心差分格式。     

边界层流动中湍斑的直接数值模拟

用Navier-Stokes方程直接数值模拟平板边界层流动中湍斑的形成和演化过程.发展了模拟湍斑的高精度、高分辨率的高效计算方法,包括推出四阶时间分裂法以提高精度;提出三维耦合差分方法,用于关于压力的泊松方程和关于速度的亥姆霍兹方程的空间离散,建立其四阶三维耦合中心差分格式;并采用四阶紧致迎风差分格式,避免了一般四阶中心差分格式不适用于边界邻域的困难和提高了分辨率;精心地处理各种边界条件,以保持精度和稳定.该方法适用于包含边界邻域的整个区域内的湍斑模拟.通过模拟平板边界层流动中湍斑的复杂演化过程,显示了湍斑的基本特征.     

Efficient kinetic schemes for steady and unsteady flow simulations on unstructured meshes

This paper presents efficient second-order kinetic schemes on unstructured meshes for both compressible unsteady and incompressible steady flows. For compressible unsteady flows, a time-dependent gas distribution function with a discontinuous particle velocity space at a cell interface is constructed and used for the evaluations of both numerical fluxes and conservative flow variables. As a result, a compact scheme on the unstructured meshes is developed. For incompressible steady flows, a continuous second-order gas-kinetic BGK type scheme is presented, for which the time-dependent gas distribution function with a continuous particle velocity is used on unstructured meshes. The efficiency of the schemes lies in the fact that the slopes of the flow variables inside each cell can be constructed using values of the flow variables within that cell only without involving neighboring cells. Therefore, even with the stencil of a first-order scheme, a high resolution method is constructed. Numerical examples are presented which are compared with the benchmark solutions and the experimental measurements.     

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.     

Hybrid weighted essentially non-oscillatory schemes with different indicators

A key idea in finite difference weighted essentially non-oscillatory (WENO) schemes is a combination of lower order fluxes to obtain a higher order approximation. The choice of the weight to each candidate stencil, which is a nonlinear function of the grid values, is crucial to the success of WENO schemes. For the system case, WENO schemes are based on local characteristic decompositions and flux splitting to avoid spurious oscillation. But the cost of computation of nonlinear weights and local characteristic decompositions is very high. In this paper, we investigate hybrid schemes of WENO schemes with high order up-wind linear schemes using different discontinuity indicators and explore the possibility in avoiding the local characteristic decompositions and the nonlinear weights for part of the procedure, hence reducing the cost but still maintaining non-oscillatory properties for problems with strong shocks. The idea is to identify discontinuity by an discontinuity indicator, then reconstruct numerical flux by WENO approximation in discontinuous regions and up-wind linear approximation in smooth regions. These indicators are mainly based on the troubled-cell indicators for discontinuous Galerkin (DG) method which are listed in the paper by Qiu and Shu (J. Qiu, C.-W. Shu, A comparison of troubled-cell indicators for Runge–Kutta discontinuous Galerkin methods using weighted essentially non-oscillatory limiters, SIAM Journal of Scientific Computing 27 (2005) 995–1013). The emphasis of the paper is on comparison of the performance of hybrid scheme using different indicators, with an objective of obtaining efficient and reliable indicators to obtain better performance of hybrid scheme to save computational cost. Detail numerical studies in one- and two-dimensional cases are performed, addressing the issues of efficiency (less CPU time and more accurate numerical solution), non-oscillatory property.     

Nonlinear evolution of turbulent spots in the near-wall shear flow

The initial model of turbulent spots with the wall disturbance using the pulse form was proposed. A group of three-dimensional coupling compact difference schemes with high accuracy and high resolution were developed, and implemented to simulate the formation and development of turbulent spots in the near-wall shear flow by means of direct numerical simulation of the Navier-Stokes equations. Growing and decaying modes were used to analyze nonlinear evolution characteristics of the turbulent spots.     

非结构网格下非定常流场双时间步长的加权ENO-强紧致杂交高分辨率格式

针对三维非定常、可压缩流场的Navier-Stokes方程组,本文提出一种新的双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间离散上不低于三阶。在对流项与粘性项的处理上,本格式分别采用了加权ENO-强紧致格式与紧致四阶精度格式的思想。几个典型算例的实践表明:计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下具有高效率与高分辨率的特征。     

一维Gross-Pitaevskii方程的高阶紧致分裂步多辛格式

首先把一维Gross-Pitaevskli方程改写成多辛Hamiltonian系统的形式,把形式通过分裂变成2个子哈密尔顿系统.然后,对这些子系统用辛或者多辛算法进行离散.通过对子系统数值算法的不同组合方式,得到不同精度的具有多辛算法特征数值格式.这些格式不仅具有多辛格式、分裂步方法和高阶紧致格式的特征,而且是质量守恒的.数值实验验证了新格式的数值行为.     

Optimized prefactored compact schemes

The numerical simulation of aeroacoustic phenomena requires high-order accurate numerical schemes with low dispersion and dissipation errors. In this paper we describe a strategy for developing high-order accurate prefactored compact schemes, requiring very small stencil support. These schemes require fewer boundary stencils and offer simpler boundary condition implementation than existing compact schemes. The prefactorization strategy splits the central implicit schemes into forward and backward biased operators. Using Fourier analysis, we show it is possible to select the coefficients of the biased operators such that their dispersion characteristics match those of the original central compact scheme and their numerical wavenumbers have equal and opposite imaginary components. This ensures that when the forward and backward stencils are added, the original central compact scheme is recovered. To extend the resolution characteristic of the schemes, an optimization strategy is employed in which formal order of accuracy is sacrificed in preference to enhanced resolution characteristics across the range of wavenumbers realizable on a given mesh. The resulting optimized schemes yield improved dispersion characteristics compared to the standard sixth- and eighth-order compact schemes making them more suitable for high-resolution numerical simulations in gas dynamics and computational aeroacoustics. The efficiency, accuracy and convergence characteristics of the new optimized prefactored compact schemes are demonstrated by their application to several test problems.     

Finite volume treatment of dispersion-relation-preserving and optimized prefactored compact schemes for wave propagation

In developing suitable numerical techniques for computational aero-acoustics, the dispersion-relation-preserving (DRP) scheme by Tam and co-workers and the optimized prefactored compact (OPC) scheme by Ashcroft and Zhang have shown desirable properties of reducing both dissipative and dispersive errors. These schemes, originally based on the finite difference, attempt to optimize the coefficients for better resolution of short waves with respect to the computational grid while maintaining pre-determined formal orders of accuracy. In the present study, finite volume formulations of both schemes are presented to better handle the nonlinearity and complex geometry encountered in many engineering applications. Linear and nonlinear wave equations, with and without viscous dissipation, have been adopted as the test problems. Highlighting the principal characteristics of the schemes and utilizing linear and nonlinear wave equations with different wavelengths as the test cases, the performance of these approaches is documented. For the linear wave equation, there is no major difference between the DRP and OPC schemes. For the nonlinear wave equations, the finite volume version of both DRP and OPC schemes offers substantially better solutions in regions of high gradient or discontinuity.     

Nonlinear spectral-like schemes for hybrid schemes

In spectral-like resolution-WENO hybrid schemes,if the switch function takes more grid points as discontinuity points,the WENO scheme is often turned on,and the numerical solutions may be too dissipative.Conversely,if the switch function takes less grid points as discontinuity points,the hybrid schemes usually are found to produce oscillatory solutions or just to be unstable.Even if the switch function takes less grid points as discontinuity points,the final hybrid scheme is inclined to be more stable,provided the spectral-like resolution scheme in the hybrid scheme has moderate shock-capturing capability.Following this idea,we propose nonlinear spectral-like schemes named weighted group velocity control(WGVC)schemes.These schemes show not only high-resolution for short waves but also moderate shock capturing capability.Then a new class of hybrid schemes is designed in which the WGVC scheme is used in smooth regions and the WENO scheme is used to capture discontinuities.These hybrid schemes show good resolution for small-scales structures and fine shock-capturing capabilities while the switch function takes less grid points as discontinuity points.The seven-order WGVC-WENO scheme has also been applied successfully to the direct numerical simulation of oblique shock wave-turbulent boundary layer interaction.     

三角翼涡破裂的高精度数值模拟

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟65°后掠角尖前缘三角翼的大攻角跨声速绕流流场,考察低耗散、高分辨率的WCNS-E-5格式对于三角翼涡破裂模拟的适用性,及激波旋涡干扰对涡破裂点位置的影响,重点研究三角翼大攻角旋涡破裂点的突然前移.通过求解任意坐标系下的非定常雷诺平均N-S方程,采用WCNS-E-5和SST两方程湍流模型,与试验结果和文献计算结果对比,表明既有高阶精度又能光滑捕捉激波的WCNS格式在模拟三角翼旋涡破裂方面具有一定优势,其数值结果与试验结果吻合较好,三角翼大攻角旋涡破裂点的突然前移是由于跨声速流场的激波旋涡干扰.