一维气液两相漂移模型的AUSMV算法研究

将AUSMV(advection upstream splitting method V)格式从计算气体动力学问题扩展至一维等温瞬态气液两相管流.阐述了采用AUSMV格式构建气液两相漂移模型数值通量的方法及边界单元的处理方法.采用Runge Kutta方法与经典的保单调MUSCL(monotone upstream centred schemes for conservation laws)方法结合Van Leer限制器,构建具有二阶时间和空间精度的数值计算方法.计算经典Zuber-Findlay激波管问题和复杂漂移关系变质量流动问题并与可靠的参考结果进行了对比.分析表明:AUSMV格式应用于气液两相流动漂移模型时计算效率高、精度高、耗散效应和色散效应小,低流速条件下能够精确地描述间断.     

一维气液两相漂移模型全隐式AUSMV算法研究

气液两相漂移模型显式AUSMV（advection upstream splitting method combined with flux vector splitting method）算法的时间步长受限于CFL（Courant-Friedrichs-Lewy）条件，为了提高计算效率，建立了一种全隐式AUSMV算法求解气液两相漂移模型.采用AUSM格式结合FVS（flux vector splitting）格式构造连续方程和运动方程的对流项数值通量，AUSM格式构造压力项数值通量.离散控制方程是非线性方程组，采用六阶Newton(牛顿)法结合数值Jacobi矩阵求解.计算经典算例Zuber-Findlay激波管问题和复杂漂移关系变质量流动问题，结果分析表明:全隐式AUSMV算法，色散效应小，无数值震荡，计算精度高.在压力波波速高的条件下，可以显著提高计算效率，耗散效应小.     

密相液固两相湍流K-ε-T模型及其在管道两相流中的应用

为了准确预测密相液固两相湍流流动,建立了Ｋ－ε－Ｔ模型,推导了控制方程组。利用该模型对竖直上升管中的密相液固两相流动进行了数值模拟,得到了与实验值吻合较好的结果。     

具弱色散误差的单调差分格式

本文基于[1]、[2]的理论分析,分别给出了计算(一维、二维)非线性双曲方程组的两类具弱色散误差的单调格式。§2针对一维气动力学方程组,通过数值结果。就激波的逼近问题,与已有的格式进行了比较。     

非牛顿流体内流传热研究

在本文中,研究了注入轴对称模腔非牛顿流体非定常流动.本文的第二部份研究了上随体Maxwell流体管内热流动.对于注入模腔流动.其本构方程采用幂律流体模型方程.为了避免在表现粘度中温度关系引起的非线性.引进了一特征粘度的概念.描述本力学过程的基本方程是,本构方程、定常状态的运动方程、非定常能量方程及连续方程.该方程组在空间是二维问题,在数学上是三维问题.采用分裂差分格式求得本方程组的数值解答.分裂法曾成功应用于求解牛顿流体问题.在本文中,首次将分裂法成功地应用解决非牛顿流体流动问题.对于圆管内热流,给出了差分格式,使基本方程组化为一个三对角方程组.其结果,给出了不同时刻的模腔内二维温度分布.     

关于非守恒形式差分格式的能量守恒问题

在建立流体力学方程组的差分格式时,对能量方程有两种不同的选择:一种是采用关于总能量(即内能与动能之和)的守恒形式的方程;另一种是采用关于内能的非守恒形式的方程.对于守恒形式的方程,容易建立能量守恒的差分格式(下面称之为守恒形式的差分格式),而对非守恒形式的方程建立的格式(下面称之为非守恒形式的差分格式),则在     

三级四阶显式辛R-K-N格式的完全构造

s级p阶辛Runge-Kutta-Nystr\"om(R-K-N)方法的一种充要条件是用关于参数的非线性方程组来表示的,辛R-K-N格式的构造问题因而转化为该方程组的求解问题. 在一些特殊的限定条件下, 已有该方程组在s=3,p=4时的两组解,即得到了两个三级四阶显式辛格式. 对于s=3,p=4情形,基于吴方法,利用计算机代数系统Maple及软件包wsolve给出了对应的非线性方程组的全部解, 这样就构造了所有的三级四阶显式辛R-K-N格式, 并证明了三级四阶显式辛R-K-N方法所满足的条件方程有冗余. 数值实验结果显示出新的辛格式在一定的条件下有着较好的误差精度.     

关于解一维抛物型方程组的差分格式

Caapck曾经研究过解多维抛物型方程组的经济格式.用他的方法解一维问题时,是将抛物型方程组的系数矩阵写成一个下三角形矩阵和一个上三角形矩阵之和,然后采用分数步长法求解.如果未知函数的个数为M,则对于每一个时间步长,需要用2M次追赶法.格式的收敛速度为Ο(τ~(1/2) h~2),这里τ和h分别为时间和空间步长.本文提出一种解一维抛物型方程组的绝对稳定的差分格式.对于每一个时间步长,求解差分方程组只要用M次追赶法,它的收敛速度为Ο(τ h~2)。     

非线性耦合Schrdinger-KdV方程组的一个局部能量守恒格式

本文利用平均值离散梯度给出了一个构造哈密尔顿偏微分方程的局部能量守恒格式的系统方法.并用非线性耦合Schrdinger-KdV方程组加以说明.证明了格式满足离散的局部能量守恒律,在周期边界条件下,格式也保持离散整体能量及系统的其它两个不变量.最后数值实验验证了理论结果的正确性.     

三维两相渗流驱动问题迎风区域分裂显隐差分法

对三维两相渗流驱动问题提出了两种迎风区域分裂显隐差分格式．压力方程采用了七点差分格式,为了能达到实际并行计算的要求,对饱和度方程采用了迎风区域分裂差分法,内边界处和各子区域分别对应显隐格式．得到了离散l2模收敛性分析,最后给出数值试验,支撑了理论分析结果．     

Some modifications of the quasilinearization method with higher-order convergence for solving nonlinear BVPs

In this paper, modifications of the quasilinearization method with higher-order convergence for solving nonlinear differential equations are constructed. A general technique for systematically obtaining iteration schemes of order m (?>?2) for finding solutions of highly nonlinear differential equations is developed. The proposed iterative schemes have convergence rates of cubic, quartic and quintic orders. These schemes were further applied to bifurcation problems and to obtain critical parameter values for the existence and uniqueness of solutions. The accuracy and validity of the new schemes is tested by finding accurate solutions of the one-dimensional Bratu and Frank-Kamenetzkii equations.     

Additive schemes for certain operator-differential equations

Unconditionally stable finite difference schemes for the time approximation of first-order operator-differential systems with self-adjoint operators are constructed. Such systems arise in many applied problems, for example, in connection with nonstationary problems for the system of Stokes (Navier-Stokes) equations. Stability conditions in the corresponding Hilbert spaces for two-level weighted operator-difference schemes are obtained. Additive (splitting) schemes are proposed that involve the solution of simple problems at each time step. The results are used to construct splitting schemes with respect to spatial variables for nonstationary Navier-Stokes equations for incompressible fluid. The capabilities of additive schemes are illustrated using a two-dimensional model problem as an example.     

SM-stability of operator-difference schemes

The spectral mimetic (SM) properties of operator-difference schemes for solving the Cauchy problem for first-order evolutionary equations concern the time evolution of individual harmonics of the solution. Keeping track of the spectral characteristics makes it possible to select more appropriate approximations with respect to time. Among two-level implicit schemes of improved accuracy based on Padé approximations, SM-stability holds for schemes based on polynomial approximations if the operator in an evolutionary equation is self-adjoint and for symmetric schemes if the operator is skew-symmetric. In this paper, additive schemes (also called splitting schemes) are constructed for evolutionary equations with general operators. These schemes are based on the extraction of the self-adjoint and skew-symmetric components of the corresponding operator.     

Flux-splitting schemes for parabolic problems

Splitting with respect to space variables can be used in solving boundary value problems for second-order parabolic equations. Classical alternating direction methods and locally one-dimensional schemes could be examples of this approach. For problems with rapidly varying coefficients, a convenient tool is the use of fluxes (directional derivatives) as independent variables. The original equation is written as a system in which not only the desired solution but also directional derivatives (fluxes) are unknowns. In this paper, locally one-dimensional additional schemes (splitting schemes) for second-order parabolic equations are examined. By writing the original equation in flux variables, certain two-level locally one-dimensional schemes are derived. The unconditional stability of locally one-dimensional flux schemes of the first and second approximation order with respect to time is proved.     

Two difference schemes for the numerical solution of Maxwell’s equations as applied to extremely and super low frequency signal propagation in the Earth-ionosphere waveguide

Two explicit two-time-level difference schemes for the numerical solution of Maxwell’s equations are proposed to simulate propagation of small-amplitude extremely and super low frequency electromagnetic signals (200 Hz and lower) in the Earth-ionosphere waveguide with allowance for the tensor conductivity of the ionosphere. Both schemes rely on a new approach to time approximation, specifically, on Maxwell’s equations represented in integral form with respect to time. The spatial derivatives in both schemes are approximated to fourth-order accuracy. The first scheme uses field equations and is second-order accurate in time. The second scheme uses potential equations and is fourth-order accurate in time. Comparative test computations show that the schemes have a number of important advantages over those based on finite-difference approximations of time derivatives. Additionally, the potential scheme is shown to possess better properties than the field scheme.     

Higher order finite difference schemes for the magnetic induction equations

We describe high order accurate and stable finite difference schemes for the initial-boundary value problem associated with the magnetic induction equations. These equations model the evolution of a magnetic field due to a given velocity field. The finite difference schemes are based on Summation by Parts (SBP) operators for spatial derivatives and a Simultaneous Approximation Term (SAT) technique for imposing boundary conditions. We present various numerical experiments that demonstrate both the stability as well as high order of accuracy of the schemes.      

Mono-implicit Runge-Kutta schemes for the parallel solution of initial value ODEs

Among the numerical techniques commonly considered for the efficient solution of stiff initial value ordinary differential equations are the implicit Runge-Kutta (IRK) schemes. The calculation of the stages of the IRK method involves the solution of a nonlinear system of equations usually employing some variant of Newton's method. Since the costs of the linear algebra associated with the implementation of Newton's method generally dominate the overall cost of the computation, many subclasses of IRK schemes, such as diagonally implicit Runge-Kutta schemes, singly implicit Runge-Kutta schemes, and mono-implicit (MIRK) schemes, have been developed to attempt to reduce these costs. In this paper we are concerned with the design of MIRK schemes that are inherently parallel in that smaller systems of equations are apportioned to concurrent processors. This work builds on that of an earlier investigation in which a special subclass of the MIRK formulas were implemented in parallel. While suitable parallelism was achieved, the formulas were limited to some extent because they all had only stage order 1. This is of some concern since in the application of a Runge-Kutta method to a system of stiff ODEs the phenomenon of order reduction can arise; the IRK method can behave as if its order were only its stage order (or its stage order plus one), regardless of its classical order. The formulas derived in the current paper represent an improvement over the previous investigation in that the full class of MIRK formulas is considered and therefore it is possible to derive efficient, parallel formulas of orders 2, 3, and 4, having stage orders 2 or 3.     

Alternating triangular schemes for convection–diffusion problems

Explicit–implicit approximations are used to approximate nonstationary convection–diffusion equations in time. In unconditionally stable two-level schemes, diffusion is taken from the upper time level, while convection, from the lower layer. In the case of three time levels, the resulting explicit–implicit schemes are second-order accurate in time. Explicit alternating triangular (asymmetric) schemes are used for parabolic problems with a self-adjoint elliptic operator. These schemes are unconditionally stable, but conditionally convergent. Three-level modifications of alternating triangular schemes with better approximating properties were proposed earlier. In this work, two- and three-level alternating triangular schemes for solving boundary value problems for nonstationary convection–diffusion equations are constructed. Numerical results are presented for a two-dimensional test problem on triangular meshes, such as Delaunay triangulations and Voronoi diagrams.     

WENO schemes with Lax–Wendroff type time discretizations for Hamilton–Jacobi equations

In this paper, a class of weighted essentially non-oscillatory (WENO) schemes with a Lax–Wendroff time discretization procedure, termed WENO-LW schemes, for solving Hamilton–Jacobi equations is presented. This is an alternative method for time discretization to the popular total variation diminishing (TVD) Runge–Kutta time discretizations. We explore the possibility in avoiding the nonlinear weights for part of the procedure, hence reducing the cost but still maintaining non-oscillatory properties for problems with strong discontinuous derivative. As a result, comparing with the original WENO with Runge–Kutta time discretizations schemes (WENO-RK) of Jiang and Peng [G. Jiang, D. Peng, Weighted ENO schemes for Hamilton–Jacobi equations, SIAM J. Sci. Comput. 21 (2000) 2126–2143] for Hamilton–Jacobi equations, the major advantages of WENO-LW schemes are more cost effective for certain problems and their compactness in the reconstruction. Extensive numerical experiments are performed to illustrate the capability of the method.     

RBF‐ENO/WENO schemes with Lax–Wendroff type time discretizations for Hamilton–Jacobi equations

In this research, a class of radial basis functions (RBFs) ENO/WENO schemes with a Lax–Wendroff time discretization procedure, named as RENO/RWENO‐LW, for solving Hamilton–Jacobi (H–J) equations is designed. Particularly the multi‐quadratic RBFs are used. These schemes enhance the local accuracy and convergence by locally optimizing the shape parameters. Comparing with the original WENO with Lax–Wendroff time discretization schemes of Qiu for HJ equations, the new schemes provide more accurate reconstructions and sharper solution profiles near strong discontinuous derivative. Also, the RENO/RWENO‐LW schemes are easy to implement in the existing original ENO/WENO code. Extensive numerical experiments are considered to verify the capability of the new schemes.