动力学的流矢量分裂法模拟浅水波方程

Based on the analogy to gas dynamics, the kinetic flux vector splitting (KFVS) method is used to stimulate the shallow water wave equations. The flux vectors of the equations are split on the basis of the local equilibrium Maxwell-Boltzmann distribution. One dimensional examples including a dam breaking wave and flows over a ridge are calculated. The solutions exhibit second-order accuracy with no spurious oscillation.     

解浅水方程离散化线性方程组的PEk方法

定义于球面的浅水方程能够很好地描述浅齐次的不可压缩非黏滞流体层的性状,它在全球大气模型、海洋数字模型和天气预报的数值计算中都有广泛的应用,浅水方程的一般形式如下：     

求解浅水波方程的熵相容格式

提出了一种求解浅水波方程组的熵相容格式．在熵稳定通量中添加特征速度差分绝对值的项来抵消解在跨过激波时所产生的熵增,从而实现熵相容．新的数值差分格式具有形式简单、计算效率高、无需添加任何的人工数值粘性的特点．数值算例充分说明了其显著的优点．利用新格式成功地模拟了不同类型溃坝问题的激波、稀疏波传播及溃坝两侧旋涡的形成,是求解浅水波方程组较为理想的方法.     

基于Gauss全局径向基函数的近岸浅水变形波高数值计算新方法

应用Gauss全局径向基函数来模拟波浪浅水变形波高变化方程中的未知函数,经实例分析探讨得到了一种可用于求解该方程数值解的新方法,并将其计算结果与常用数值分析方法得到的数值解相互对比印证,证明了基于Gauss全局径向基函数法计算结果的正确性.经验证,Gauss径向基函数法的平均计算误差相比其他方法均要小,表明该方法拥有更高的计算精度.同时,根据Gauss全局径向基函数的逼近结果,得出了浅水变形波高变化微分方程数值解的拟合函数,在实际工程中,可以利用该拟合函数来代替原方程的解析解,研究成果可为求解近岸浅水区域波浪运动提供一种新思路.     

声波方程高阶差分方法数值模拟结果研究

首先,给出了声波方程时空域高阶差分格式.证明了差分格式在时间上是二阶精度,在空间上最高可以达到2M精度.其次,对该高阶差分格式优缺点评价,新的差分格式是一种关于一个方向上所用到的临近网格节点数目一半的可变长度的差分算子,为不同精度需求提供了灵活可变的差分算子公式.另外一种意义上说,新的差分格式是依赖于频率的,同一种差分格式对于不同频率的震源,其频散的程度也不一样,因此只有采用了合适的差分形式才能够达到最好的效果.最后,给出了一系列数值模拟实例,验证了文中所得到的结论.     

Burgers方程的混合元分析及其数值模拟

1．引言混合有限元法在高阶偏微分方程和含有两个战者两个以上）的未知国数的偏微分方程的数值解的研究中起着重要的作用．但是,到目前为止,混合有限元法主要是用于2n阶或一阶偏微分方程（组）,如二阶椭圆型方程、平面弹性力学方程、双调和方程、Stokes和Navier－stokes方程、抛物型方程以及电磁场方程修见＞到以及当中的参考文献）．然而,R前混合有限元法还没有被用于对非线性的Burgers方程作数值研究．而过去对Burgers方程的数值研究主要采用标准有限元法、差分方法和谱方法修见【IO－12］以及当中的参考文献）．本文的目的是用混…     

基于N-S方程的裂缝酸刻蚀数值模拟

酸化压裂是碳酸盐岩储层增产改造的主要措施.通过酸压数值模拟预测碳酸盐岩裂缝的刻蚀形态对酸化压裂施工有指导意义.酸压数值模拟多以立方定律为基础,然而随着油气勘探技术进步,油气开采难度增加,酸压数值模拟需要更清晰的流场分布和更高的精度.针对此问题,以N-S方程为基础,推导裂缝酸刻蚀数学模型,随后采用有限差分求解速度场,迎风...     

浅水波方程的一种基于特征方向的Galerkin方法

本文给出求解二维浅水波方程组的一种Crank-Nicolson型基于特征方向的Galerkin有限元方法，证明了该方法的一个误差估计结果，并给出了该方法的一个算例．     

浅水波方程的初边值问题

本文用Kato关于拟线性演化方程的初值问题的存在性定理证明了浅水波方程在半无界直线上初值问题局部解的存在性。用解的先估计证明了整体解的存在性或解的Blow-up性质，并给出了解关于x的渐近估计。     

基于拟Shannon小波浅水长波近似方程组的数值解

本文研究了浅水长波近似方程组初边值问题的数值解.利用小波多尺度分析和区间拟Shannon小波,对浅水长波近似方程组空间导数实施空间离散,用时间步长自适应精细积分法对其变换所的非线性常微分方程组进行求解,得到了浅水长波近似方程组的数值解,并将此方法计算的结果与其解析解进行比较和验证.     

THE GLOBAL ARTIFICIAL BOUNDARY CONDITIONS FOR NUMERICAL SIMULATIONS OF THE 3D FLOW AROUND A SUBMERGED BODY

We consider the numerical approximations of the three-dimensional steady potential flow around a body moving in a liquid of finite constant depth at constant speed and distance below a free surface in a channel. One vertical side is introduced as the up-stream artificial boundary and two vertical sides are introduced as the downstream arti-ficial boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the original problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. After solving the variational problem by the finite element method, we obtain the numerical approximation of the original problem. The numerical examples show that the artificial boundary conditions given in this paper are very effective.     

粘性浅水方程柯西问题解的Lp-估计

本文通过对相应线性方程格林函数的估计，而获得粘性浅水方程柯西问题解的Lp渐近估计。     

THE APPROXIMATIONS OF THE EXACT BOUNDARY CONDITION AT AN ARTIFICIAL BOUNDARY FOR LINEARIZED INCOMPRESSIBLE VISCOUS FLOWS

1.IntroductionManyproblemsarisinginfluidmechanicsaregiveninanunboundeddomain,suchasfluidflowaroundobstacles.Whencomputingthenumericalsolutionsoftheseproblems,oneoftenintroducesartificialboundariesandsetsupaxtificialboundaryconditionsonthem.Thentheoriginal…     

带非线性边界条件的反应扩散方程的数值方法

１引言近年来关于非线性抛物型方程数值解法的研究取得了许多好的结果,其中以Ｃ．Ｖ．Ｐａｏ为主的研究者们利用上、下解方法对带线性边界条件的半线性抛物型方程的有限差分系统进行了广泛的研究,提出了一系列有效的迭代算法（见［１］、［２］、［３］、［４］）．但对带非线性边界条件的半线性抛物型方程初边值问题,作者至今尚未见到有研究者将上、下解方法用在相应的差分系统上,求得数值解．其主要原因是由于边界上函数的非线性,解在边界网格点上的值未知且无法用内部网格点上的值直接表示,相应的差分系统表示形式受到影响,边界网…     

ARTIFICIAL BOUNDARY METHOD FOR THE THREE-DIMENSIONAL EXTERIOR PROBLEM OF ELASTICITY

The exact boundary condition on a spherical artificial boundary is derived for thethree-dimensional exterior problem of linear elasticity in this paper. After this bound-ary condition is imposed on the artificial boundary, a reduced problem only defined in abounded domain is obtained. A series of approximate problems with increasing accuracycan be derived if one truncates the series term in the variational formulation, which isequivalent to the reduced problem. An error estimate is presented to show how the errordepends on the finite element discretization and the accuracy of the approximate problem.In the end, a numerical example is given to demonstrate the performance of the proposedmethod.     

A NUMERICAL STUDY FOR THE PERFORMANCE OF THE WENO SCHEMES BASED ON DIFFERENT NUMERICAL FLUXES FOR THE SHALLOW WATER EQUATIONS

In this paper we investigate the performance of the weighted essential non-oscillatory （WENO） methods based on different numerical fluxes, with the objective of obtaining better performance for the shallow water equations by choosing suitable numerical fluxes. We consider six numerical fluxes, i.e., Lax-Friedrichs, local Lax-Friedrichs, Engquist-Osher, Harten-Lax-van Leer, HLLC and the first-order centered fluxes, with the WENO finite volume method and TVD Runge-Kutta time discretization for the shallow water equations. The detailed numerical study is performed for both one-dimensional and two-dimensional shallow water equations by addressing the property, and resolution of discontinuities. issues of CPU cost, accuracy, non-oscillatory     

ON DISCRETE PROJECTION AND NUMERICAL BOUNDARY CONDITIONS FOR THE NUMERICAL SOLUTION OF THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

1. IntroductionLet us consider the unsteady incompressible Navier--Stokes equations (INSE)on a two--dimensional rectangular region fl with boundary 0fl. Here w = (u, v)" is tl1e velocityvector, p is the pressure, and f a known vector function of x) y, and…     

A POSTERIORI ERROR ESTIMATES OF A NON-CONFORMING FINITE ELEMENT METHOD FOR PROBLEMS WITH ARTIFICIAL BOUNDARY CONDITIONS

An a posteriori error estimator is obtained for a nonconforming finite element approximation of a linear elliptic problem, which is derived from a corresponding unbounded domain problem by applying a nonlocal approximate artificial boundary condition. Our method can be easily extended to obtain a class of a posteriori error estimators for various conforming and nonconforming finite element approximations of problems with different artificial boundary conditions. The reliability and efficiency of our a posteriori error estimator are rigorously proved and are verified by numerical examples.