二维恒定各向同性介质渗透系数反演的遗传算法

给出了利用遗传算法求解二维恒定各项同性介质渗透系数反演的一种新方法,该方法把参数反演问题转化为优化问题通过遗传算法求解.数值模拟结果表明:该方法具有精度高、收敛速度快、编程简单、易于计算机实现等优点,值得在实际工作采用.     

正交各向异性轴对称位势问题的Trefftz有限元分析

Trefftz有限元法（Trefftz finite element method, TFEM）因其独特的优良品质而备受关注.针对正交各向异性轴对称位势问题,提出了一种4节点四边形环状单元．在该单元模型中,首先假设两套独立的位势插值模式:即单元域内场和网线场,然后代入修正变分泛函并利用Gauss散度定理消除区域积分,最后根据驻值原理导得只含边界积分的单元刚度方程．数值算例表明了该单元的准确性、稳定性以及对网格畸变的不敏感性.     

求解二维波动方程正演反演问题的半离散方法

本文用半离散方法将二维波动方程离散为一维耦合波动方程组．给出了离散的收敛性及波动方程组的适定性．利用这种方法可以求解波动方程系数及演问题．     

二维抛物型方程参数反演的遗传算法

本文研究了二维抛物型方程参数反演问题.利用遗传算法求解此反演问题的方法,把参数反演问题转化为优化问题,通过演化计算方法求解.它从多个初始点开始寻优,借助交叉和变异算子来获得参数的全局最优解.且数值模拟结果表明,具有精度高、编程简单、易于计算机实现等特点.     

二维稳态热传导问题的正六边形流形元研究

发展了用于分析二维稳态热传导问题的多边形数值流形方法（numerical manifold method,NMM）.根据热传导问题的控制方程、边界条件以及多边形NMM的温度近似函数,采用修正变分原理导出了多边形NMM求解稳态热传导问题的总体方程,给出了多边形单元上的域积分策略.考虑到NMM中数学覆盖系统可不与物理域边界一致以及规则单元的精度优势,采用Wachspress正六边形数学单元对两个典型热传导问题进行了仿真,计算结果与参考解能较好地吻合,表明多边形NMM可以很好地模拟平面稳态热传导问题.     

多孔介质驱动问题的混合元最小二乘特征有限元方法及其收敛分析

１．引言多孔介质二相驱动问题的数学模型是偶合的非线性偏微分方程组的初边值问题．该问题可转化为压力方程和浓度方程［１－４］．浓度方程一般是对流占优的对流扩散方程,它的对流速度依赖于比浓度方程的扩散系数大得多的Ｆａｒｃｙ速度．因此Ｄａｒｃｙ速度的求解精度直接影响着浓度的求解精度．为了提高速度的求解精度,７０年代Ｐ．Ａ．Ｒａｖｉａｔ和Ｊ．Ｍ．Ｔｈｏｍａｓ提出混合有限元方法［５］．Ｊ．ＤｏｕｇｌａｓＪｒ,Ｔ．Ｆ．Ｒｕｓｓｅｌｌ,Ｒ．Ｅ．Ｅｗｉｎｇ,Ｍ．Ｆ．Ｗｈｅｅｌｅｒ［１］－［４］,［９］,［１２］袁…     

二维非稳态对流扩散边界控制问题的简化算法

针对二维非稳态对流扩散边界控制问题计算量大的问题,提出了基于降阶模型的最优实时控制方法.利用POD(the Proper Orthogonal Decomposition)和奇异值分解以及Galerkin投影方法得到了具有高精度离散形式的状态空间降阶模型.在所得的降阶状态空间模型中,利用离散时间线性二次调节器方法设计出了最优控制器.对流-扩散过程的控制模拟结果说明了所提方法的有效性和准确性.     

二维变系数对流扩散问题的特征配置法

1引言 在多种问题的数值模拟中均涉及抛物型对流扩散方程的数值求解问题.由于配置法 无需计算数值积分，计算简便，收敛阶高等优点，使之在工程技术和计算数学的许多领域 得到广泛的应用，但范围一般局限在一维常系数{1，21和二维常系数问题降，4]，90年代[s] 提出了二维变系数     

一类二维单个守恒律方程的Riemann问题

本文研究一类二维单个守恒律方程的Riemann问题.用广义特征分析方法研究了这类方程,给出了基本波的分类,解决了初值为两片常数的二维Riemann问题,给出了Riemann解.     

ERROR ESTIMATES FOR THE RECURSIVE LINEARIZATION OF INVERSE MEDIUM PROBLEMS

This paper is devoted to the mathematical analysis of a general recursive linearization algorithm for solving inverse medium problems with multi-frequency measurements. Under some reasonable assumptions, it is shown that the algorithm is convergent with error estimates. The work is motivated by our effort to analyze recent significant numerical results for solving inverse medium problems. Based on the uncertainty principle, the recursive linearization allows the nonlinear inverse problems to be reduced to a set of linear problems and be solved recursively in a proper order according to the measurements. As an application, the convergence of the recursive linearization algorithm [Chen, Inverse Problems 13（1997）, pp.253-282] is established for solving the acoustic inverse scattering problem.     

A TRUST REGION METHOD FOR SOLVING DISTRIBUTED PARAMETER IDENTIFICATION PROBLEMS

This paper is concerned with the ill-posed problems of identifying a parameter in an elliptic equation which appears in many applications in science and industry. Its solution is obtained by applying trust region method to a nonlinear least squares error problem.Trust region method has long been a popular method for well-posed problems. This paper indicates that it is also suitable for ill-posed problems, Numerical experiment is given to compare the trust region method with the Tikhonov regularization method. It seems that the trust region method is more promising.     

QUANTUM COMPLEXITY OF THE INTEGRATION PROBLEM FOR ANISOTROPIC CLASSES

We obtain the optimal order of high-dimensional integration complexity in the quantumcomputation model in anisotropic Sobolev classes W∞^r ([0, 1]^d) and Hǒlder Nikolskii classes H∞^r([0, 1]^d). It is proved that for these classes of functions there is a speed-up of quantum algorithms over deterministic classical algorithms due to factor n^-1 and over randomized classical methods due to factor n^-1/2. Moreover, we give an estimation for optimal query complexity in the class H∞^∧ (D) whose smoothness index is the boundary of some complete set in Z ^d.     

两类各向异性非协调元的某些超收敛性质分析

在各向异性网格下,讨论了两类非协调矩形元对二阶椭圆边值问题的某些超逼近性和超收敛性,并证明了在单元中心点这种超收敛性仅为一种点态现象.数值结果验证了我们理论分析的正确性.     

AN ITERATIVE ALGORITHM FOR DETERMINING TIME－DEPENDENT COEFFICIENT OF TWO DIMENSIONAL LINEAR WAVE EQUATION

ＡＮＩＴＥＲＡＴＩＶＥＡＬＧＯＲＩＴＨＭＦＯＲＤＥＴＥＲＭＩＮＩＮＧＴＩＭＥ－ＤＥＰＥＮＤＥＮＴＣＯＥＦＦＩＣＩＥＮＴＯＦＴＷＯＤＩＭＥＮＳＩＯＮＡＬＬＩＮＥＡＲＷＡＶＥＥＱＵＡＴＩＯＮ（苏超伟）西北工业大学，邮编：７１００７２ＳｕＣｈａｏｗｅｉ（Ｎ...     

三维种群生态动力系统的参数识别

研究具有反应一扩散现象的三维种群生态动力系统的参数识别问题，依该系统正问题解的性质，建立了参数识别的数学模型；论证了系统正问题解关于待识别参数的连续依赖性与参数识别问题最优解的存在性。     

ACM元各向异性分析的Newton方法

本文将一维Lagrange插值多项式的Newton表达式推广到二维非标准的Hermite插值,给出著名板元-ACM元插值多项式的Newton表达式,由此给出ACM元对四阶和二阶椭圆问题的各向异性插值误差估计,为复杂单元的各向异性分析开辟了新的途径.     

AN ANISOTROPIC NONCONFORMING FINITE ELEMENT METHOD FOR APPROXIMATING A CLASS OF NONLINEAR SOBOLEV EQUATIONS

An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approximate schemes, which are the same as the traditional finite element methods. In addition, the global superconvergence is derived through the postprocessing technique. Numerical experiments are included to illustrate the feasibility of the proposed method.     

VORTEX DYNAMICS OF THE ANISOTROPIC GINZBURG-LANDAU EQUATION

In this article, using coordinate transformation and Gronwall inequality, we study the vortex motion law of the anisotropic Ginzburg-Landau equation in a smooth bounded domain Ω  R2, that is, ■tuε = 2Σ/j,k=1 (ajkj■xjkuε)xj + b(x)(1-ε|ε|2)uε/2u, x ∈Ω, and conclude that each vortex bj(t) (j=1, 2,···, N) satisfies dbdjt(t)= -(a1k(bj(t)b)■(xk))a(a (bj(t)), a2k(bj (t))/xk a(bj (t)) a(bj (t)) , where a(x) =(a11a22-a122(1/2)). We prove that all the vortices are pinned together to the critical points of a(x). Furthermore, we prove that these critical points can not be the maximum points.