关于高阶精度WCNS格式的无粘通量分裂方法

高阶精度加权紧致非线性格式(WCNS)越来越广泛地应用于复杂流动数值模拟.WCNS可以与多种无粘通量分裂方法结合起来使用.但是,常见的通量分裂方法都是基于低阶格式发展起来的,目前还不清楚哪些通量分裂方法最适合WCNS,也不知道这些方法与高阶格式结合时将会产生什么效果.表面热流计算是高超声速流动数值模拟的难点之一,为了在热流计算时选择合适的通量,研究了多种通量分裂方法的耗散大小.每种通量都可以表示成中心部分与耗散部分之和.这些通量的中心部分相同且非常简单,但是耗散部分较为复杂,且不同的通量分裂方法可导致不同的耗散表达式.通过对通量耗散进行分析可以发现耗散大小与网格界面两侧的物理量跳跃近似线性正相关.数值计算表明高阶格式得到的网格界面左右两侧的物理量跳跃通常远比低阶格式小,因而带来的通量耗散小.通过3个典型算例考察了通量耗散对热流计算的影响,其中包括高超激波/边界层干扰算例.基于对van Leer通量、Steger-Warming通量、KFVS通量、Roe通量、AUSM类通量和HLL类通量的考察,给出了通量选择建议.     

低Mach数翼型绕流数值模拟的低耗散预处理格式

低Mach数流动中,基于可压缩流动的数值模拟算法存在严重的刚性问题,预处理方法可以有效地解决这一问题,但其计算结果不稳定．基于原有的预处理Roe格式,引入可调节参数,得到一种新的低耗散格式．该格式可以减弱边界层以及极低速区域的过度耗散,使得整个流场计算稳定．低Mach数、低Reynolds数定常圆柱绕流和低Mach数、高Reynolds数翼型(NACA0012和S809)绕流3个验证算例表明,带可调节参数的低耗散预处理方法正确可靠,是低速流动数值模拟的有效方法．     

High-Order Bicompact Schemes for Shock-Capturing Computations of Detonation Waves

Computational Mathematics and Mathematical Physics - An implicit scheme with splitting with respect to physical processes is proposed for a stiff system of two-dimensional Euler gas dynamics...     

脉冲波作用下的高超音速非定常流动模拟及影响研究

采用高阶精度有限差分方法模拟了快声波脉冲扰动作用下的高超音速非定常流场,分析了脉冲波与高超音速流场的相互干扰,并应用Fourier频谱分析研究扰动波在边界层的发展．结果表明:来流脉冲扰动波与激波及边界层强烈相互作用,弓形激波明显向内弯曲,激波后扰动波被显著放大;来流扰动波与弓形激波干扰形成的边界层外的扰动波和近壁面内形成的边界层扰动波存在明显分界．钝锥头部参数扰动幅值要远大于其他位置参数扰动幅值．在边界层内的发展阶段,一些扰动模态持续增长,一些扰动模态被过滤掉,不再增长,甚至衰减,而也有一些扰动模态先衰减再增长．总的来说,在钝锥头部低频扰动模态为主导模态,随着扰动从流场上游向下游发展,总扰动模态中的低频模态成份和高频模态成份所占的比例开始转变,高频模态成分显著地增大．     

High-Order Numerical Methods for Electromagnetic Induction

The present paper establishes an axisymmetric benchmark model of a conducting loop, which implies an electromagnetic induction. Therein, the fully coupled MAXWELL equations are demonstrated in a monolithic solution strategy. This dynamic problem is solved with a high order finite element discretization using GALERKIN's method in space as well as in time. Furthermore, high order RUNGE-KUTTA time integration methods are analyzed. Studies regarding an h error estimation and the order of convergence are examined. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Compact High-Order Finite Differences for Interface Problems in One Dimension

This paper is concerned with the construction and analysis ofcompact finite difference approximations to the model linearsource problem –(pu')' + qu = f where the functions p,q, and f can have jump discontinuities at a finite number ofpoints. Explicit formulae that give O(h²) O(h³) and O(h⁴) accuracyare derived, and a procedure for computing three-point schemesof any prescribed order of accuracy is presented. A rigoroustruncation and discretization error analysis is offered. Numericalresults are also given.     

三维对流扩散方程的三种高精度分裂格式

在算子分裂法思想的基础上,将两种高精度的离散格式推广应用于三维对流扩散方程,同时对经典ADI格式的对流项做了改进,改进后的格式的对流项对空间具有4阶精度,而经典ADI格式对空间只有2阶精度,由此可见,提高了该格式的实用性．最后对两种典型的浓度场进行了数值模拟,将3种格式的计算结果与解析解以及其它传统差分格式的计算结果进行了对比,得出当Peclet数不大于5时,3种格式均获得了令人满意的数值结果,说明推广的这三种方法具有很高的准确性和可靠性．     

高精度隐式参差光滑格式的构造

参照Lax-Wendroff格式的构造方法,就双曲型方程、抛物型方程和双曲-抛物型方程,构造了一种新的IRS(implicit residual smoothing)格式。该IRS格式有二阶或三阶时间精度且大大地拓宽了解的稳定区域和CFL数。这种新的IRS格式有中心加权型和迎风偏向型两种,并用von-Neumann分析方法分析了格式的稳定范围。讨论了在透平机械中广泛应用的Dawes方法的局限性,发现该方法对稳态问题得出的解与时间步长的选取有关,对粘性问题求解时,时间步长受严格限制。最后,结合TVD(total variation diminishing)格式和四阶Runge-Kutta技术,用IRS格式和Dawes方法对二维反射激波场进行了数值模拟,数值结果支持本文的分析结论。     

求解时间分布阶扩散方程的两个高阶有限差分格式

基于复化Simpson公式和复化两点Gauss-Legendre公式，构造了两个求解时间分布阶扩散方程的高阶有限差分格式.不同于以往文献中提出的时间一阶或二阶格式，这两种格式在时间方向都具有三阶精度，而在分布阶和空间方向可达到四阶精度.数值结果表明，两种算法都是稳定且收敛的，从而是有效的.两种格式的收敛速率也通过数值实验进行了验证，并且通过和文献中的算法对比可以得出其更为高效,     

带激波的一维非定常两相流动的数值模拟

本文将处理带激波的单相气体非定常流动问题的两种高分辨数值方法(随机取样法和二阶GRP有限差分法)推广应用于气固悬浮体(亦称含灰气体)两相情况,计算了含灰气体激波管中两相激波特性、波后流场结构及气固两相流动参数随时间的变化.数值结果表明:这两种方法均能给出带有尖锐间断阵面的两相激波松弛结构.二阶GRP方法在计算精度和机时耗用等方面优于随机取样法.     

Uniform Convergence of Variational Methods for Plane Subsonic Flow

Uniform convergence is proved for certain direct variationalmethods, assuming the existence of a solution with certain properties.In particular one important property used is that of "subsonicity".This is then imposed on the approximating functions in orderto establish the uniform convergence. This overcomes a crucialassumption made in previous work for Rayleigh–Ritz approximations.further, previous results for bounded subdomains are extentedto infinity. The analysis applies mroe generally to plane uniformlyelliptic problems when a priori bounds on the first derivativesof the solution are known to exist.     

湍流度对翼型绕流影响的数值模拟及与实验的对比

研究了NACA0012翼型在低Reynolds数不同湍流度下的空气动力学特性．既运用了有限谱法和QUICK格式相结合的数值计算方法,又通过实验手段,研究了湍流度对流场的影响,并得到了计算与实验相符的结果．计算结果说明了有限谱法的高精度及与其他格式结合后应用的灵活性．对翼型受力和流场的结果分析显示:低湍流度下无明显的失速特征,高湍流度下失速特征明显;湍流度还对边界层后的剪切层有很大影响,高湍流度的受力特征和流场结果与高Reynolds数下的情形相似．     

球Couette流的数值模拟

该文利用谱方法对同心旋转球间轴对称Couette流进行数值模拟.给出Navier Stokes方程的流函数涡度形式,利用Stokes流把边界条件齐次化, 选取Stokes算子的特征函数做为逼近子空间的基函数,对同心旋转球间轴对称Couette流进行谱逼近     

Some Stability Inequalities for Compact Finite Difference Schemes

For finite difference schemes of compact form on nonuniform grids approximating m-th order two-point boundary value problems stability inequalities are proved which use a norm analogous to the Spijker-norm in the case of multistep methods. The results are applied to a number of finite difference schemes for which they establish a higher order of convergence than naively expected.     

The Alternating-Direction Schemes and Numerical Analysis for the Three-dimensional Seawater Intrusion Simulation

Both numerical simulation and theoretical analysis of seawater intrusion in coastal regions are of great theoretical importance in environmental sciences. The mathematical model can be described as a coupled system of three dimensional nonlinear partial differential equations with initial-boundary value problems. In this paper, according to the actual conditions of molecular and three-dimensional characteristic of the problem, we construct the characteristic finite element alternating-direction schemes which can be divided into three continuous one-dimensional problems. By making use of tensor product algorithm, and priori estimation theory and techniques, the optimal order estimates in H1 norm are derived for the error in the approximate solution.     

Sensitivities Calculation for Unsteady Flow Problems

In this work we present an efficient approach for the calculation of the sensitivities in the case of unsteady problems like the compressible high-speed flows. The sensitivities are of a known response with respect to the initial conditions. If the number of parameters is large, the adjoint method is preferred since the cost of computing the sensitivities is independent of the number of design parameters. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Importance of Reasonable Numerical Schemes in Lake Circulation Simulation

Three difference methods ‐ two traditional central and upstream difference schemes (CDS and UDS) and a highresolution total variation diminishing (TVD) method ‐ are applied to discretize the one‐dimensional advection problem as well as the advection terms in the three‐dimensional advection‐diffusion lake circulation problem. It is demonstrated that the choice of reasonable numerical schemes is fairly importance. The traditional numerical treatments with respect of the advection terms in the lake dynamical equations are fairly arguable, although they are still most frequently employed in lake circulation models, and the TVD technique is most preferable and reasonable.     

New Sixth-Order Compact Schemes for Poisson/Helmholtz Equations

Some new sixth-order compact finite difference schemes for Poisson/Helmholtz equations on rectangular domains in both two- and three-dimensions are developed and analyzed. Different from a few sixth-order compact finite difference schemes in the literature, the finite difference and weight coefficients of the new methods have analytic simple expressions. One of the new ideas is to use a weighted combination of the source term at staggered grid points which is important for grid points near the boundary and avoids partial derivatives of the source term. Furthermore, the new compact schemes are exact for 2D and 3D Poisson equations if the solution is a polynomial less than or equal to 6. The coefficient matrices of the new schemes are $M$-matrices for Helmholtz equations with wave number $K≤0,$ which guarantee the discrete maximum principle and lead to the convergence of the new sixth-order compact schemes. Numerical examples in both 2D and 3D are presented to verify the effectiveness of the proposed schemes.     

A High-Order Finite Difference Scheme for 3D Unsteady Convection Diffusion Reaction Equations北大核心CSCD

针对三维非稳态对流扩散反应方程,构造了一种高精度紧致有限差分格式,对空间的离散采用四阶紧致差分方法,对时间的离散采用Taylor级数展开和余项修正技术,所提格式在时间上的精度为二阶、在空间上的精度为四阶。利用Fourier稳定性分析法证明了该格式是无条件稳定的。最后给出数值算例验证了理论结果。