数值天气模式参数选择的全局优化方法

１引言数值天气预报模式中关于参数的选择直接影响到天气预报的准确率,在建立一个数值天气预报系统时,为了得到好的预报效果,必须对模式参数进行优化．在这方面已有许多文献［１］－［７］作过有益的探讨,提供了许多有效的方法,在文献［２］中,给出了一种参数反演的方法．并应用广义线性反演,获得较稳定的计算格式．然而,此方法在每一次迭代时,至少需要解ｎ＋１个正问题（其中ｎ为参数的个数）．又在文献［６］中。引进了四维同化的共轭梯度法,适宜于求解高维问题．然而,共轭梯度法只能求得局部最优解,对初始参数的选取很敏感,…     

由于波响应反演声阻抗的特征带方法

1．引言子波激发下的反演问题通常是不适定的，如何构造稳定、高效的算法是反问题研究中的重要课题．当前的波动方程反演方法主要有两类：特征线方法和最优化方法［1］．特征线方法是数值求解波动方程反问题的一种重要而有效的方法，它的基本思想是沿着波动方程上、下行波的特征传播方向逐层推进，并按照因果律求解．关于这方面的早期工作可参看［7］．在［2］中证明了脉冲激发下一维波动方程系数反问题的适定性，为这一方法提供了理论基础．随后，［4］讨论了特征线方法的差分计算的收敛性，［5，6］提供了成功的数值计算实例．近来人们逐…     

延迟动力系统线性θ-方法的散逸性

１．引言 科学与工程中的许多问题具有散逸性,即系统具有一有界吸引集,从任意初始条件出发的解经过有限时间后进入该吸引集并随后保持在里面．如 ２维的 Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程、Ｌｏｒｅｎｚ方程等许多重要系统都是散逸的．散逸性研究一直是动力系统研究中的重要课题（参见Ｔｅｍａｍ［７］）．当数值求解这些系统时,自然希望数值方法能保持系统的该重要特性．１９９４年, Ｈｕｍｐｈｒｉｅｓ和 Ｓｔｕａｒｔ［６］首次研究了 Ｒｕｎｇｅ－Ｋｕｔｔａ方法对有限维系统的散逸性．１９９７年Ｈｉｌｌ［２］研究了其无穷维散逸性…     

一个求解多维守恒律方程组的二阶显式有限元格式

１．引言 近年来,在非结构网格上求解双曲型守恒律的数值方法引起了较为广泛的关注,出现了有限体积方法［１］,间断 Ｇａｌｅｒｋｉｎ方法 ［２］,流线扩散方法［３］,以及 ＮＮＤ格式 ［４］等．我们在［６,７］中提出了一种求解双曲型守恒律方程式的有限元方法,它是在一个求解对流扩散问题的有限元方法 ［５］的基础上发展起来的．它是一个显式有限元方法,因此计算量很小．在这个方法中,我们将任意维的问题归结为在单元棱边上的一维计算,引入了积分因子,因此在单元内部可以容纳边界层．这样,它特别适合于对流占优问题以及双曲…     

浅水流动与污染物扩散的高分辨率计算模型

将组合型TVD格式应用于守恒型的浅水方程和污染物扩散方程，建立了二者耦合求解的高分辨率有限体积计算模型。给出了溃坝水流、明渠突扩流和污染物输运计算的典型算例，并与实验数据或其它数值结果进行了比较，证实了该模型的有效性，表明它不但能处理有激波的非恒定流问题，也能较好地计算具有任意边界的一般的浅水流动和污染物扩散问题，为浅水流动和水环境模拟提供了精度高、稳定性好、普适性强的数值方法。     

求解二维浅水方程的一种高分辨率有限体积法

本文将一种van Albada型可微的限制器函数引入到二维浅水方程的求解中，发展了一种求解二维浅水方程的有限体积法．数值实验结果表明，该方法不仅计算精度高，而且较其它求解二维浅水方程的高精度有限体积法，在数值解的收敛性能方面大有改善．     

W空间中有界线性算子的最佳逼近及其应用

１．引言 对线性算子的有限秩算子逼近是最经典的问题．并且它的应用极广．如数值积分公式、函数的逼近、数值原函数、方程的数值解法等．１９８６年,在文［１］中,首次给出了在再生核空间中函数的最佳逼近算子（恒等算子的有限秩算子逼近）．之后;在文［２］中给出了数值原函数．又在文［３］、［５］、［６］等中利用有限秩算子逼近（并非是最佳逼近）给出了一些方程的数值解法．但这些讨论都是在一元函数空间上只对特殊算子进行的．１９９７年,虽然在文［４］中给出了完备的二元再生核空间及二元函数的最佳逼近插值算子．但是对多元…     

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

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

θ-方法的非线性渐近稳定性

１引言 数值求解延迟微分方程时,方法的稳定性具有无容置疑的重要性．自１９７５年Ｂａｒｗｅｌｌ引入Ｐ－稳定性与ＧＰ－稳定性概念以来,该领域研究已获许多重要成果（如［７］［８］）．它们大多是基于下面标量线性模型方程：其中λ,μ为复数且满足延迟量τ（＞０）为常数,函数θ（ｔ）连续． 我们首先回忆Ｂａｒｗｅｌｌ［１］的定义． 定义１．１一个数值方法称为是Ｐ－稳定的,如果对任意正整数ｒ用该方法按步长ｈ＝τ／ｒ求解（１．１）时在节点ｔｎ＝ｎｈ的数值解ｙｎ满足 定义１．２一个数值方法称为是ＧＰ－稳定的,如果用该方…     

数值求解延时微分方程的步长准则

１．引言 用一个数值方法求解下列延时微分方程：其中, ｆ： Ｒ × Ｃｄ × Ｃｄ → Ｃｄ为给定函数, Ｕ（ｔ）当上＞ ０时为未知函数,τ＞ ０为常数延时量,ф（ｔ）∈Ｃｄ为已知向量值函数．为了检验一个数值方法的数值稳定性,常用如下试验方程：来观察方法的数值稳定性,这里ａ,ｂ∈Ｃ（Ｃ为复数集）为已知常数,ф（ｔ）为给定的连续函数（ｔ≤０）． 定义 １［2］．延时微分方程（简记为ＤＤＥＳ）（３）被称为是渐近稳定的,如果（３）的每一个解Ｕ（ｔ）满足 方程(３)的特征方程为： 定义 ２［２］．一数值方法求解ＤＤＥＳ称为…     

The Boundary Value Problem for the Nonlinear Shallow Water Equations

We propose an approximate analytical solution of the boundary value problem (BVP) for the nonlinear shallow waters equations. Our work, based on the Carrier and Greenspan [ 1 ] hodograph transformation, focuses on the propagation of nonlinear nonbreaking waves over a uniformly plane beach. Available results are briefly discussed with specific emphasis on the comparison between the Initial Value Problem and the BVP; the latter more completely representing the physical phenomenon of wave propagation on a beach. The solution of the BVP is achieved through a perturbation approach solely using the assumption of small waves incoming at the seaward boundary of the domain. The most significant results, i.e., the shoreline position estimation, the actual wave height and velocity at the seaward boundary, the reflected wave height and velocity at the seaward boundary are given for three specific input waves and compared with available solutions.     

Local Artificial Boundary Conditions for the Incompressible Viscous Flow in a Slip Channel

1.IntroductionManyboundaxyvaJueproblemsofpartialdiffereotialequationsinvo1vingunboundeddomainoccurinmanyareasofapplications,e-g.lfluidflowaroundobstacles,couplingofstructureswithfoundationandsoon.Forgettingthenumericalsolutionsoftheproblemsonunboundeddomian,anaturalapproachistocutoffanunboundedpartofthedomainbyintroducinganartificialboundaryandsetupanaPpropriatear-tificialboundaryconditiononthearti%ialboundaryThentheoriginalproblemisapproximatedbyaproblemonbou.d.dfdomain.Inthelastteny6aJrs,b…     

Analysis of artificial boundary conditions for exterior boundary value problems in three dimensions

Summary. Simple boundary conditions on an artificial boundary are discussed, then an exact boundary condition on the artificial boundary is obtained. Approximation to this boundary condition with high accuracy is given, and the error estimates are obtained. A numerical example is presented, and the numerical results are compared with the exact solution. Received January 27, 1997 / Revised version received May 14, 1999 / Published online February 17, 2000     

NUMERICAL SOLUTIONS OF PARABOLIC PROBLEMS ON UNBOUNDED 3-D SPATIAL DOMAIN

In this paper,the numerical solutions of heat equation on 3-D unbounded spatial do-main are considered. n artificial boundary Γ is introduced to finite the computationaldomain.On the artificial boundary Γ,the exact boundary condition and a series of approx-imating boundary conditions are derived,which are called artificial boundary conditions.By the exact or approximating boundary condition on the artificial boundary,the originalproblem is reduced to an initial-boundary value problem on the bounded computationaldomain,which is equivalent or approximating to the original problem.The finite differencemethod and finite element method are used to solve the reduced problems on the finitecomputational domain.The numerical results demonstrate that the method given in thispaper is effective and feasible.     

THE ARTIFICIAL BOUNDARY CONDITIONS FOR NUMERICAL SIMULATIONS OF THE COMPLEX AMPLITUDE IN A COUPLED BAY-RIVER SYSTEM

We consider the numerical approximations of the complex amplitude in a coupled bayriver system in this work. One half-circumference is introduced as the artificial boundary in the open sea, and one segment is introduced as the artificial boundary in the river if the river is semi-infinite. On the artificial boundary a sequence of high-order artificial boundary conditions are proposed. Then the original problem is solved in a finite computational domain, which is equivalent to a variational problem. The numerical approximations for the original problem are obtained by solving the variational probiem with the finite element method. The numerical examples show that the artificial boundary conditions given in this work are very effective.     

Error bounds for the finite element approximation of an incompressible material in an unbounded domain

Summary. In this paper we design high-order local artificial boundary conditions and present error bounds for the finite element approximation of an incompressible elastic material in an unbounded domain. The finite element approximation is formulated in a bounded computational domain using a nonlocal approximate artificial boundary condition or a local one. In fact there are a family of nonlocal approximate artificial boundary conditions with increasing accuracy (and computational cost) and a family of local ones for a given artificial boundary. Our error bounds indicate how the errors of the finite element approximations depend on the mesh size, the terms used in the approximate artificial boundary condition and the location of the artificial boundary. Numerical examples of an incompressible elastic material outside a circle in the plane is presented. Numerical results demonstrate the performance of our error bounds. Received August 31, 1998 / Revised version received November 6, 2001 / Published online March 8, 2002     

EXACT NONREFLECTING BOUNDARY CONDITIONS FOR AN ACOUSTIC PROBLEM IN THREE DIMENSIONS

In this paper,nonreflecting artificial boundary,conditions are considered for an acoustic problem in three dimensions.With the technique of Fourier dexompostition under the orthogonal basis of spherical harmonics,three kinds of equivalemt exact artificial boundary conditions are obtained on a spherical artificail boundary.A numerical test is presented to show the performance of the method.     

Splitting as an approach to constructing local exact artificial boundary conditions

We employ the technique of splitting for constructing artificial boundary conditions (ABCs) for the linear advection–diffusion–reaction equation when the computational domain is an nD open set with a piecewise smooth artificial boundary. The splitting is performed both by the physical processes and by coordinates. The former permits to construct ABCs for each of the processes separately, which provides local exact boundary conditions; the latter leads to ABCs much less exigent to the shape of artificial boundary in comparison with many others. We also prove that the corresponding boundary value problems are well-posed, and present results of the numerical experiments that confirm the theoretical study.     

Error estimates for the finite element approximation of linear elastic equations in an unbounded domain

In this paper we present error estimates for the finite element approximation of linear elastic equations in an unbounded domain. The finite element approximation is formulated on a bounded computational domain using a nonlocal approximate artificial boundary condition or a local one. In fact there are a family of nonlocal approximate boundary conditions with increasing accuracy (and computational cost) and a family of local ones for a given artificial boundary. Our error estimates show how the errors of the finite element approximations depend on the mesh size, the terms used in the approximate artificial boundary condition, and the location of the artificial boundary. A numerical example for Navier equations outside a circle in the plane is presented. Numerical results demonstrate the performance of our error estimates.

     

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.