7阶WENO-S格式的计算效率研究

WENO-S格式是一类适合于含间断问题数值模拟的加权本质无振荡格式.这类格式的光滑因子满足对单频波为常数,这使得其近似色散关系与线性基底格式一致,并且具有良好的小尺度波动模拟能力．计算效率是数值方法性能指标的一个重要方面.由于WENO-S格式的光滑因子在各子模板上的计算公式除下标不同外形式一致,在计算线性对流方程相邻数值通量时,部分光滑因子完全相同．为此提出一种消除WENO-S格式冗余光滑因子计算的方法.该方法要求一条网格线上用于重构或插值的量可以表示为一个序列.基于此要求分析其对于几种不同物理问题的可行性和使用方法.以7阶WENO-S格式为例介绍了格式性质和去除冗余光滑因子计算的方法.该方法中预先计算和存储一条网格线上的所有光滑因子,在网格点较多的情况下,光滑因子计算次数约为原7阶WENO-S格式的1/4.对一维对流问题、球面波传播问题、二维旋转问题、二维小扰动传播问题及一维和二维无黏流动问题进行了数值模拟.结果表明该格式对多种流动结构具有良好的捕捉能力,并且同时具有良好的计算效率,去除冗余计算后又降低了约20%的计算时间.     

一类非线性奇异积分方程及其数值方法研究

探讨了一类非线性奇异积分方程的数理性质以及在颗粒雷诺数Rep＜1时此类方程解的存在条件，然后详细研究了该方程的数值计算方法并构造称之为P（EC）^k多步法的差分格式，分析了该格式的收敛性和代数精度，得到时域离散步长的约束关系。运用该格式计算了静止流场和均匀振荡流中球形小颗粒的非恒定运动，将计算结果与其解析解及有关实验数据的比较表明，上述数值方法具有良好的计算精度。     

A new weighted upwind finite volume element method based on non‐standard covolume for time‐dependent convection–diffusion problems

In modern numerical simulation of problems in energy resources and environmental science, it is important to develop efficient numerical methods for time‐dependent convection–diffusion problems. On the basis of nonstandard covolume grids, we propose a new kind of high‐order upwind finite volume element method for the problems. We first prove the stability and mass conservation in the discrete forms of the scheme. Optimal second‐order error estimate in L₂‐norm in spatial step is then proved strictly. The scheme is effective for avoiding numerical diffusion and nonphysical oscillations and has second‐order accuracy. Numerical experiments are given to verify the performance of the scheme. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical analysis and construction of limiter of high resolution difference scheme

I.IntroductionAdjustingandgoverningproperlythenumericaldissipationanddispersionarethekeytotheconstructionofhighresolutionandnon-oscillationscheme.SinceVanLcerll]introducednuxlimiterandobtainedhighresolutionandnon-oscillationscheme,choisillgandconstructingsuitablelimiterbecomesanimportantwayofdesigninghighresolutiollalld11on-oscillationscheme,andadjustsnumericaldissipationeffects,dispersioneffectsaswellasgroupvelocityeffects.Roe121'l'l,Chakravarthy'andOSherl4]andetal.providedvariouslimitersan…     

Numerical method for the system of reaction-diffusion equations with a small parameter

This paper deals with the numerical method for the system of reaction-diffusionequations with a small parameter.It is difficult to solve the problems of this kindnumerically because of the boundary layer efect.Besed on singular perturbed theory andGreens’function,we have established the difference scheme that is suited for the solution tothe problems.It e introduce an idea of feasilbe equidistant degree a here.And this provesthat if a≥2.the scheme converges in l,m norm with speed O.h \t)uniformly.     

机器学习在求解一维双曲守恒律方程中的应用

双曲守恒律方程对空气动力学、物理学和海洋学等众多领域问题的计算有着重大意义,本文应用机器学习框架下的BP神经网络对双曲守恒律方程近似求解.首先,采用熵稳定格式及基于自适应移动网格的熵稳定格式所得多个时间层的数值解构造网络输入,采用高分辨率熵稳定格式所得对应的多个时间层的数值解构造网络输出,并对数据集作归一化处理.随后,...     

一类非线性发展方程的交替分段显隐并行数值方法

给出了一类非线性发展方程的交替分段显隐并行数值方法 ,得到了方法的无条件稳定性和并行性兼顾的结果。数值例子说明理论分析的正确性和格式的有效性     

Use of characteristics method with cubic interpolation for unsteady-flow computation

The specified-time-interval (STI) scheme has been used commonly in applying the method of characteristics (MOC) to unsteady open-channel flow problems. However, with the use of STI scheme, the numerical error for the simulation results can always be induced due to the interpolation used to approximate the characteristics trajectory. Hence, in order to remedy the numerical errors caused by the interpolation, one needs to seek some kind of interpolation technique with higher-order accuracy. Instead of the linear interpolation technique, which has been used very commonly and can induce serious numerical diffusion, the Holly--Preissmann two-point, method, which is a cubic interpolation technique with fourth-order of accuracy, is proposed here to integrate with the method of characteristics for the computation of one-dimensional unsteady flow in open channel. The concept of reachback and reachout in space and time directions for the characteristics is also introduced to assure the model stability. The computed results from this new model are compared with those computed by using the Preissmann four-point scheme and the multimode method of characteristics with linear interpolation.     

精细时程积分法及其数值衍生格式应用评估

旋翼气动弹性耦合动力学方程本质上是一组刚性比较大的非线性偏微分方程。在有限元结构离散后，可改写为非齐次微分方程组，其中非齐次项是桨叶运动量（位移与速度）和气动载荷的函数。针对这类方程，本文尝试引入精细积分法及其衍生格式，借助数值方法计算Duhamel积分项。从积分精度与数值稳定性方面比较研究具有代表性的精细库塔法和高精度直接积分法。结合隐式积分算法，评估精细积分法应用于旋翼动力学方程的可行性。算例表明，精细积分法对矩形直桨叶动力学方程具有足够的求解精度。     

Extension of an explicit finite volume method to large time steps (CFL>1): application to shallow water flows

In this work, the explicit first order upwind scheme is presented under a formalism that enables the extension of the methodology to large time steps. The number of cells in the stencil of the numerical scheme is related to the allowable size of the CFL number for numerical stability. It is shown how to increase both at the same time. The basic idea is proposed for a 1D scalar equation and extended to 1D and 2D non‐linear systems with source terms. The importance of the kind of grid used is highlighted and the method is outlined for irregular grids. The good quality of the results is illustrated by means of several examples including shallow water flow test cases. The bed slope source terms are involved in the method through an upwind discretization. Copyright © 2005 John Wiley & Sons, Ltd.     

概率密度演化方程TVD格式的自适应时间步长技术及其初值条件改进

随机性普遍存在于实际工程问题中,而复杂结构的非线性随机响应分析是其中的一个难点,近年发展的概率密度演化方法为此类问题的求解提供了新的途径.由于实际问题的复杂性,概率密度演化方程通常采用数值方法求解,因此提高计算效率和求解精度对实际应用具有重要意义.本文基于变网格技术,推导了概率密度演化方程在非均匀时间步长上的总变差减小(total variation diminishing,TVD)差分格式,算例结果表明通过自适应插值可将迭代次数减少为原来的43.4%,当随机过程样本持续时间增大时均值估计的平均误差基本不变,而标准差估计的平均误差不断增大,但增大幅度不断减小;计算耗时随样本持续时间的增大也呈增大趋势,而由于使用了时间步长自适应插值算法导致有些情况下长持时样本的计算耗时反而比短持时样本的计算耗时短;在传统的脉冲函数型初值条件基础上,提出了一种高阶导数更稳定的余弦函数型初值条件形式.结果表明,脉冲函数型的初值条件是余弦函数型初值条件的一个特例,当参数取值适当时,余弦函数型初值条件的数值求解结果具有更高的精度.本文的工作进一步完善了概率密度演化方程的求解方法,为其在实际工程中的应用提供了基础.      

带时滞的二元神经元网络的OPNCL控制和复杂网络的时间延迟反馈控制

针对一种带时滞的二元神经元网络和一种复杂网络的混沌控制问题,利用开环加非线性闭环(Open Plus Nonlinear Closed Loop,OPNCL)方法和时间延迟反馈控制(Time Delay FeedbackControl,TDFC)方法,分别设计了该混沌网络和混沌系统的控制器。从理论上证明了第一种控制器可使该网络系统的解稳定地传递到选定的目标,并通过数值模拟实验进一步验证了两种方法的有效性。结果表明:该方法避免了开环加非线性闭环控制的一些限制因素;对于任何目标,所控制混沌系统的传递域(Basins of Entrainment)是全局的,避免了有关确定传递域范围的繁琐计算。     

非线性双曲型守恒律的两步二阶TVD差分格式

通过Mac Cormack格式和Warming-Beam的结合,构造了一种非常简单的两步二阶TVD差分格式,该差分格式更适合于使用分量形式差分计算而无须对欧拉方程组进行特征解耦。通过对流体力学方程组的大量数值试验,并与二阶ENO格式进行了比较,充分显示了该格式高精度、高分辨并且极其简单的优良特性。     

双开口压力谐振管的数值计算

本文研究两端开口的压力谐振管的数值计算方法。文中采用格杜诺夫差分格式,对文献[1]中的喷射边界条件作了改进,提出了双开口边界条件的处理方法。计算结果与实验基本一致。计算方法可用于多种形式的准周期性一维非定常流动,如双开口热分离器变压管的设计和特性分析。     

Hybrid finite analytic solutions of shallow water circulation

IntroductionThepredictionoftidalandstormsurgesisessentialfornavigation ,fishingandprotectionofthecoast.Furtherthedeterminationofhydrodynamiccirculationbecomesnecessaryforthemodelingofdispersionofpollutants,sedimenttransportandsalinitydistribution .Recently ,somekindsofnumericalmodelingtechniqueshavebeenextensivelydevelopedinthelasttwodecades.Inordertoimprovethewaterqualityinabays/lagoons ,manymethodshavebeenproposedbyresearchers.Forexample ,Komatsuetal.(2 0 0 0 ) [1] proposedanewmethodtocontr…     

Exponential dichotomies of nonlinear discrete systems and its application to numerical analysis and computation

In this pepar we consider the upwind difference scheme of a kind of boundary value problems for nonlinear, second order, ordinary differential equations. Singular perturbation method is applied to construct the asymptotic approximation of the solution to the upwind difference equation. Using the theory of exponential dichotomies we show that the solution of an order-reduced equation is a good approximation of the solution to the upwind difference equation except near boundaries. We construct correctors which yield asymptotic approximations by adding them to the solution of the order-reduced equation. Finally, some numerical examples are illustrated.     

Tensor Universal Serendipity Elements and unsteady Taylor-Galerkin finite element method

A new kind of Universal Serendipity Element (USE)— the Tensor Universal Serendipity Element (TUSE) is constructed by using both tensor force finite elements and the basic idea of USE. The formulation of shape functions and their derivatives for TUSE is presented. TUSE can be used to study steady and unsteady transonic flow fields when combined with Taylor-Galerkin Finite Element Methods, the NND scheme in FDM, and four-stage Runge-Kutta methods. As numerical examples the transonic flow in cascades and one kind of complex unsteady transonic axisymmetric flow in engineering are studied. It is shown that the algorithm presented in this paper is efficient and robust. The project supported by the National Natural Science Foundation of China     

Simulation and application of three-dimensional migration-accumulation of oil resources

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｉｌｆｏｒｍａｔｉｏｎｉｎｓｅｄｉｍｅｎｔｂａｓｉｎｓ,ｉｔｓｄｉｓｐｌａｃｅｍｅｎｔ,ｔｒａｎｓｐｏｒｔａｎｄａｃｃｕｍｕｌａｔｉｏｎ,ａｎｄｔｈｅｆｉｎａｌｆｏｒｍａｔｉｏｎｏｆｏｉｌｄｅｐｏｓｉｔｓｈａｖｅｂｅｅｎｔｈｅｋｅｙｐｒｏｂｌｅｍｓｉｎｔｈｅｅｘｐｌｏｒａｔｉｏｎｏｆｏｉｌ＿ｇａｓｒｅｓｏｕｒｃｅｓ．Ｈｏｗｈａｓｏｉｌｂｅｅｎａｃｃｕｍｕｌａｔｅｄｉｎｔｈｅｐｒｅｓｅｎｔｌｏｏｐａｃｃｏｒｄｉｎｇｔｏｔｈｅｍｅｃｈａｎｉｃｓｏｆｉｍｍｉｓｃｉｂｌｅ…     

求解可压缩流的二维通量分裂格式

传统的一维通量分裂格式在计算界面数值通量时,只考虑网格界面法向的波系。采用传统的TV格式分别求解对流通量和压力通量。通过求解考虑了横向波系影响的角点数值通量来构造一种真正二维的TV通量分裂格式。在计算一维数值算例时,该格式与传统的TV格式具有相同的数值通量计算公式,因此其保留了传统的TV格式精确捕捉接触间断和膨胀激波的优点。在计算二维算例时,该格式比传统的TV格式具有更高的分辨率;在计算二维强激波问题时,消除了传统TV格式的非物理现象,表现出更好的鲁棒性;此外,该格式大大提高了稳定性CFL数,从而具有更高的计算效率。因此,本文方法是一种精确、高效并且具有强鲁棒性的数值方法,在可压缩流的数值模拟中具有广阔的应用前景。     

NUMERICAL SIMULATION OF OIL MIGRATION-ACCUMULATION OF MULTILAYER AND ITS APPLICATION

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｉｌｆｏｒｍａｔｉｏｎｉｎｓｅｄｉｍｅｎｔｂａｓｉｎｓ ,ｉｔｓｄｉｓｐｌａｃｅｍｅｎｔ,ｔｒａｎｓｐｏｒｔａｎｄａｃｃｕｍｕｌａｔｉｏｎ ,ａｎｄｔｈｅｆｉｎａｌｆｏｒｍａｔｉｏｎｏｆｏｉｌｄｅｐｏｓｉｔｓｈａｖｅｂｅｅｎｏｎｅｏｆｔｈｅｋｅｙｐｒｏｂｌｅｍｓｉｎｔｈｅｅｘｐｌｏｒａｔｉｏｎｏｆｏｉｌ＿ｇａｓｒｅｓｏｕｒｃｅｓ.Ｈｏｗｈａｓｏｉｌｂｅｅｎａｃｃｕｍｕｌａｔｅｄｉｎｔｈｅｐｒｅｓｅｎｔｌｏｏｐａｃｃｏｒｄｉｎｇｔｏｔｈｅｍｅｃｈａｎｉｃｓｏｆｉｍｍｉｓｃ…