A high-order accurate hybrid scheme using a central flux scheme and a WENO scheme for compressible flowfield analysis

A high-order accurate hybrid central-WENO scheme is proposed. The fifth order WENO scheme [G.S. Jiang, C.W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys. 126 (1996) 202–228] is divided into two parts, a central flux part and a numerical dissipation part, and is coupled with a central flux scheme. Two sub-schemes, the WENO scheme and the central flux scheme, are hybridized by means of a weighting function that indicates the local smoothness of the flowfields. The derived hybrid central-WENO scheme is written as a combination of the central flux scheme and the numerical dissipation of the fifth order WENO scheme, which is controlled adaptively by a weighting function. The structure of the proposed hybrid central-WENO scheme is similar to that of the YSD-type filter scheme [H.C. Yee, N.D. Sandham, M.J. Djomehri, Low-dissipative high-order shock-capturing methods using characteristic-based filters, J. Comput. Phys. 150 (1999) 199–238]. Therefore, the proposed hybrid scheme has also certain merits that the YSD-type filter scheme has. The accuracy and efficiency of the developed hybrid central-WENO scheme are investigated through numerical experiments on inviscid and viscous problems. Numerical results show that the proposed hybrid central-WENO scheme can resolve flow features extremely well.     

随机扰动对拟小波方法求解对流扩散方程的影响

引进拟小波方法数值求解对流扩散方程，研究结果表明，计算带宽W有一个极值，当计算带宽W取该极值时，该方程的拟小波解的精度最高，且好于迎风格式。当边界发生随机不等幅扰动时，对于积分时间较长的情况，拟小波格式的效果要稍逊于迎风格式；当边界发生随机等幅扰动时，若计算带宽W取大于等于20的整数时，方程拟小波解的精度与迎风格式相同；当参数受到随机扰动时，W取10时的拟小波解的均方根误差要小于迎风格式；在初值发生随机扰动且计算带宽W取10时，方程的拟小波解的精度最高，好于迎风格式。     

A Fifth-Order Low-Dissipative Conservative Upwind Compact Scheme Using Centered Stencil

In this paper, a conservative fifth-order upwind compact scheme using centered stencil is introduced. This scheme uses asymmetric coefficients to achieve the upwind property since the stencil is symmetric. Theoretical analysis shows that the proposed scheme is low-dissipative and has a relatively large stability range. To maintain the convergence rate of the whole spatial discretization, a proper non-periodic boundary scheme is also proposed. A detailed analysis shows that the spatial discretization implemented with the boundary scheme proposed by Pirozzoli [J. Comput. Phys., 178 (2001), pp. 81-117] is approximately fourth-order. Furthermore, a hybrid methodology, coupling the compact scheme with WENO scheme, is adopted for problems with discontinuities. Numerical results demonstrate the effectiveness of the proposed scheme.     

柱几何Rayleigh-Taylor不稳定性的数值模拟

 给出了柱几何中流体力学方程组及其在数值模拟中采用的计算方法。对二维柱几何Rayleigh-Taylor不稳定性进行数值模拟，在线性阶段与线性理论符合得很好；不稳定性增长进入非线性区域的阈值依赖于界面的位置，并且明显不同于平面情况。     

金属锡Rayleigh-Taylor不稳定性对模型参数敏感性的数值分析

利用自研的爆轰与冲击动力学欧拉计算程序和Steinberg-Guinan(SG)本构模型,数值模拟分析了样品初始参数(初始振幅、初始波长、样品初始厚度)和SG本构模型初始参数对爆轰驱动锡Rayleigh-Taylor(RT)不稳定性增长的影响。结果表明金属锡样品的初始参数对其RT不稳定性增长有很大的影响。RT不稳定性增长随着初始振幅的减小而减小,且存在一个截止初始振幅;存在一个最不稳定的模态(波长),当初始波长大于该波长时,RT不稳定性增长随着初始波长的减小而增大,反之,RT不稳定性增长随着初始波长的减小而减小;样品厚度的增大可以抑制RT不稳定性增长,而且存在一个样品截止厚度。金属锡的RT不稳定性增长对其SG本构模型应变硬化系数和应变硬化指数的变化不敏感,而对压力硬化系数和热软化系数比较敏感。从采用扰动增长法预估材料强度的角度来说,修正压力硬化系数以获得锡合理的材料强度是合理的途径。     

Aggregation Behaviors of a Two-Species System with Lose-Lose Interactions

We propose an aggregation evolution model of two-species （A- and B-species） aggregates to study the prevalent aggregation phenomena in social and economic systems. In this model, A- and B-species aggregates perform self-exchange-driven growths with the exchange rate kernels K（k, l） = Kkl and L（k, l） = Lkl, respectively, and the two species aggregates perform self-birth processes with the rate kernels J1（k） = J1 k and J2（ k ） = J2k, and meanwhile the interaction between the aggregates of different species A and B causes a lose-lose scheme with the rate kernel H（k,l） = Hkl. Based on the mean-field theory, we investigated the evolution behaviors of the two species aggregates to study the competitions among above three aggregate evolution schemes on the distinct initial monomer concentrations A0 and B0 of the two species. The results show that the evolution behaviors of A- and B-species are crucially dominated by the competition between the two self-birth processes, and the initial monomer concentrations Ao and Bo play important roles, while the lose-lose scheme play important roles in some special cases.     

积分梯度法的两种格式

在二维柱坐标系下Lagrange流体力学的计算中,积分梯度法是动量方程的一种有效离散方法.积分梯度法中,IGT(Integral Gradient Total)格式不能保持柱几何下一维球对称性;IGA(Integral Gradient Average)格式可以保持一维球对称性,但当相邻网格质量相差比较大时,会得到远远脱离真实物理现象的加速度.深入研究IGA和IGT格式发现,当相邻网格边界压力取为质量加权时,即使相邻网格质量相差较大,对于一维平面和一维柱问题,IGT与IGA等价;在二维情形下,可以缩小IGT和IGA之间的差异.理论证明,IGA格式不能保持系统的动量守恒,IGT格式能保持系统的动量守恒性.数值模拟结果进一步显示了这两个格式的优缺点.     

一维强场模型研究中的非齐线性正则方程的辛算法

就一维强场模型,采用对称差商代替空间变量的2阶偏导数,将含有SchrÖdinger方程的初边值问题离散成"非齐线性正则方程",它的齐方程的通解和非齐方程特解都由"辛变换生成",分别采用辛格式计算.采用这种辛算法和R-K法计算了一个数值例子,并与精确解作了比较.结果表明,经长时间计算后,辛算法保持解的固有特征,而R-K法则面目全非.     

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

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

价键优选法及其在纳米结构预测与物性研究中的应用

本文论述价键优选法作为一种新颖的理论方法在材料结构预测与物性研究中的应用, 特别是在低维数纳米结构如团簇与纳米线研究中所展示的优势. 价键优选法以原子几何构型和电子云(主要是由费米能级附近的分子轨道组成, 即广义前线轨道)空间分布来合理决定纳米结构的稳定构型的选取. 本文以硅团簇为例说明价键优选法的特点, 以及锂、钠、铍、镁等金属团簇为例说明价键优选法在结构预测与材料 物性随尺寸演化规律研究中的应用, 以锂离子在MoS纳米线中的吸附为例说明价键优选法在储能材料离子传导研究中的应用, 最后总结价键优选法的进一步发展方向. 关键词： 价键优选法 团簇 纳米结构 离子传导     

Large-eddy simulation of turbulent channel flows with conservative IDO scheme

The resolution of a numerical scheme in both physical and Fourier spaces is one of the most important requirements to calculate turbulent flows. A conservative form of the interpolated differential operator (IDO-CF) scheme is a multi-moment Eulerian scheme in which point values and integrated average values are separately defined in one cell. Since the IDO-CF scheme using high-order interpolation functions is constructed with compact stencils, the boundary conditions are able to be treated as easy as the 2nd-order finite difference method (FDM). It is unique that the first-order spatial derivative of the point value is derived from the interpolation function with 4th-order accuracy and the volume averaged value is based on the exact finite volume formulation, so that the IDO-CF scheme has higher spectral resolution than conventional FDMs with 4th-order accuracy. The computational cost to calculate the first-order spatial derivative with non-uniform grid spacing is one-third of the 4th-order FDM. For a large-eddy simulation (LES), we use the coherent structure model (CSM) in which the model coefficient is locally obtained from a turbulent structure extracted from a second invariant of the velocity gradient tensor, and the model coefficient correctly satisfies asymptotic behaviors to walls.     

自忆模式中差分格式的稳定性研究

基于大气运动是一种不可逆过程的观点，引进了忆及过去时次资料的记忆函数，导出热传导的自忆性方程，研究了方程分别取Richardson和DuFort-Frankel格式，回溯阶p取1时的稳定性.探讨了多时刻模式中数值积分有时发散的问题，揭示了由过去时次资料动态求取记忆函数，改变了原定设计的差分格式，且它是一个时间平滑因子的本质. 关键词： 回溯阶差分格式 多时次 稳定性 记忆性     

Large deformation measurement scheme for 3D digital image correlation method

Difficulties often arise for digital image correlation (DIC) technique when serious de-correlation occurs between the reference image and the deformed image due to large deformation. An updating reference image scheme could be employed to deal with large deformation situation, however that will introduce accumulated errors. A large deformation measurement scheme, combining improved coarse search method and updating reference image scheme, is proposed in this paper. For a series of deformation images, the correlation calculation begins with a seed point and spreads out. An improved coarse search method is developed to calculate the initial correlation parameters for the seed point, which guarantees that the correlation calculation can be carried out successfully even in large deformation situation. Only for extremely large deformation, the reference image is updated. Using this method, not only extremely large deformation can be measured successfully but also the accumulated error could be controlled. A polymer material tensile test and a foam compression test are used to verify the proposed scheme. Experimental results show that up to 450% tensile deformation and 83% compression deformation can be measured successfully.     

A cell-centered Lagrangian scheme with the preservation of symmetry and conservation properties for compressible fluid flows in two-dimensional cylindrical geometry

We develop a new cell-centered control volume Lagrangian scheme for solving Euler equations of compressible gas dynamics in cylindrical coordinates. The scheme is designed to be able to preserve one-dimensional spherical symmetry in a two-dimensional cylindrical geometry when computed on an equal-angle-zoned initial grid. Unlike many previous area-weighted schemes that possess the spherical symmetry property, our scheme is discretized on the true volume and it can preserve the conservation property for all the conserved variables including density, momentum and total energy. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the performance of the scheme in terms of symmetry, accuracy and non-oscillatory properties.     

A class of asymptotic-preserving schemes for kinetic equations and related problems with stiff sources

In this paper, we propose a general time-discrete framework to design asymptotic-preserving schemes for initial value problem of the Boltzmann kinetic and related equations. Numerically solving these equations are challenging due to the nonlinear stiff collision (source) terms induced by small mean free or relaxation time. We propose to penalize the nonlinear collision term by a BGK-type relaxation term, which can be solved explicitly even if discretized implicitly in time. Moreover, the BGK-type relaxation operator helps to drive the density distribution toward the local Maxwellian, thus naturally imposes an asymptotic-preserving scheme in the Euler limit. The scheme so designed does not need any nonlinear iterative solver or the use of Wild Sum. It is uniformly stable in terms of the (possibly small) Knudsen number, and can capture the macroscopic fluid dynamic (Euler) limit even if the small scale determined by the Knudsen number is not numerically resolved. It is also consistent to the compressible Navier–Stokes equations if the viscosity and heat conductivity are numerically resolved. The method is applicable to many other related problems, such as hyperbolic systems with stiff relaxation, and high order parabolic equations.     

A NEW DIFFERENCE SCHEME WITH MULTI-TIME LEVELS

In view of making the best use of information coming from past observational data, a new difference scheme with multi-time levels (p>3) is suggested. Some mathematical characteristics of the scheme, which is called the retrospective scheme, are discussed. The numerical results of some examples show that the calculation accuracy of linear and nonlinear advection equations computed with the retrospective scheme is higher than that obtained via the leapfrog scheme. The scheme can be applied to many fields, such as meteorology, engineering physics, astronautics, environment and economy etc, where systematic observations are made normally.     

2阶精度时间后差隐式格式TVD充分条件

非定常流动问题计算中常用到含三个时间层的二阶精度时间后差隐式格式,并且希望构成TVD格式,然而理论上的问题多年一直没有解决。本文找到了解决办法,构造了这种类型的隐式TVD格式,证明了其为TVD的充分条件。理论结果为计算所验证,并表明通常未采取本文对时间差分处理方法的格式尚不具备TVD性质。     

Robust HLLC Riemann solver with weighted average flux scheme for strong shock

Many researchers have reported failures of the approximate Riemann solvers in the presence of strong shock. This is believed to be due to perturbation transfer in the transverse direction of shock waves. We propose a simple and clear method to prevent such problems for the Harten–Lax–van Leer contact (HLLC) scheme. By defining a sensing function in the transverse direction of strong shock, the HLLC flux is switched to the Harten–Lax–van Leer (HLL) flux in that direction locally, and the magnitude of the additional dissipation is automatically determined using the HLL scheme. We combine the HLLC and HLL schemes in a single framework using a switching function. High-order accuracy is achieved using a weighted average flux (WAF) scheme, and a method for v-shear treatment is presented. The modified HLLC scheme is named HLLC–HLL. It is tested against a steady normal shock instability problem and Quirk’s test problems, and spurious solutions in the strong shock regions are successfully controlled.     

求解Euler方程的高精度的空间--时间守恒格式

空间-时间守恒(STC)格式是近年来发展出的一种计算格式,在现有的STC格式构造过程中,流动变量在解元中的分布都用其一阶Taylor展开式来表示.STC格式的精度与所采用的Taylor展开式的阶数有关.该文采用流动变量的二阶Taylor展开式来表示其在解元上的分布、构造出了求解一维Euler方程的STC格式.用该格式对几个问题进行了计算,将计算结果与精确解进行了比较,比较表明该格式有较高的精度.     

非正交网格上的九点格式在热传导问题中的改进

对文献[1]在二维Lagrange流体力学网格上构造的扩散方程九点差分格式做了进一步的讨论和改进,给出了一般形式边界条件的计算格式。数值试验的结果表明,这些改进提高了原格式数值结果的精度。