二次锥规划的预估-校正光滑方法

基于光滑Fischer-Burmeister函数,给出一个求解二次锥规划的预估-校正光滑牛顿法.该算法构造一个等价于最优性条件的非线性方程组,再用牛顿法求解此方程组的扰动.在适当的假设下,证明算法是全局收敛且是局部二阶收敛的.数值试验表明算法的有效性.     

旋转流体中圆球沿旋转轴运动产生的扰动

本文用线化理论分析了整体旋转的理想流体中有一个圆球沿旋转轴作匀速运动时流体的扰动,基于旋轴对称流动的假设导出了决定运动稳定性的扰动压力方程和扰动流函数方程.用简正模法分析了扰动流函数方程,得出了非平凡中性扰动的波数与波速必须满足的约束条件,并求出了扰动的精确表达式.文中得出结论,中性扰动共有三种可能的形式.     

求解广义非线性互补问题的光滑化拟牛顿法

利用光滑对称扰动Fischer-Burmeister函数将广义非线性互补问题转化为非线性方程组,提出新的光滑化拟牛顿法求解该方程组.然后证明该算法是全局收敛的,且在一定条件下证明该算法具有局部超线性(二次)收敛性.最后用数值实验验证了该算法的有效性.     

具有常外力项的零压欧拉方程组初值涉及狄拉克函数的黎曼问题

研究具有常外力项的一维零压欧拉方程组初值涉及狄拉克函数的黎曼问题.首先,研究对应的扰动初值问题;其次,基于一维非线性守恒律方程组弱解的稳定性理论,通过分析对应的扰动初值问题解的极限,最终得到了6类显示解.特别地,对于某些初值,捕捉到了在密度变量和内能变量上同时涉及狄拉克函数的狄拉克接触间断,这是一类重要的非线性现象.此外,这些结果清晰地刻画了常外力项对解结构的影响.     

外区域上可压缩Euler方程组的有限时间破裂

本文研究高维(n=2, 3)外区域中可压缩Euler方程组初边值问题解的长时间行为.假定初始密度和速度场在常状态附近有一个紧支集小扰动,本文证明带不可渗透边界条件的初边值问题没有整体解,并进一步给出解关于初值小扰动参数的生命跨度上界估计;引入一种比较简单的试探函数方法,并结合能量恒等式证明主要结论.     

求解大型非对称线性方程组的不完全广义最小向后扰动法

本文给出了求解大型非对称线性方程组的广义最小向后扰动法(GMBACK)的截断版本——不完全广义最小向后扰动法(IGMBACK).该方法基于Krylov向量的不完全正交化,从而在Krylov子空间上求出一个近似的或者拟最小向后扰动解.本文对新算法IGMBACK做了一些理论研究,包括算法的有限终止、解的存在性和唯一性等方面的研究;且给出了IGMBACK的执行.数值实验表明:IGMBACK通常比GMBACK和广义最小残量法(GMRES)更有效;且IGMBACK和GMBACK经常比GMRES收敛得更好.特殊地,如果系数矩阵是敏感矩阵,且方程组右侧的向量平行于系数矩阵的最小奇异值对应的左奇异向量时,重新开始的GMRES不一定收敛,而IGMBACK和GMBACK一般收敛,且比GMRES收敛得更好.     

圆锥规划问题的光滑牛顿方法

给出求解圆锥规划问题的一种新光滑牛顿方法.基于圆锥互补函数的一个新光滑函数,将圆锥规划问题转化成一个非线性方程组,然后用光滑牛顿方法求解该方程组.该算法可从任意初始点开始,且不要求中间迭代点是内点.运用欧几里得代数理论,证明算法具有全局收敛性和局部超线性收敛速度.数值算例表明算法的有效性.     

求解病态线性方程组的模拟退火算法

病态方程组的条件数较大,当输入数据有微小扰动或计算过程中的舍入误差都可能引起输出数据的很大扰动,使得解严重失真,因此求解此类方程组是相当困难的.本文尝试使用模拟退火算法来求解病态线性方程组,得到了较好的结果,并与传统的求解方法作了简单的比较.     

一类带有空间时间白噪音随机弹性方程的全离散差分格式

随机弹性方程在结构工程中有许多应用.本文研究一类由空间时间白噪音扰动的随机弹性方程的全离散有限差分格式.通过引入新的函数,将随机弹性方程表示成一阶方程组的形式,然后对噪音项进行分片常数逼近,构造了带有空间时间白噪音随机弹性方程的全离散差分格式.基于对Gronwall不等式和Burkholder不等式的应用,证明了格式的L~p收敛性并得到了收敛阶.在数值实验中结合Monte-Carlo方法,所得实验结果与理论分析是一致的.     

奇异摄动对流扩散方程组的参数一致的二阶有限差分格式

针对一类奇异摄动对流扩散方程组,在■网格上构造了经典的迎风有限差分格式,并利用闸函数方法证明了数值方法为一阶收敛.在此基础上,设计了一个Richardson外推格式,并严格证明了外推方法的精度为二阶一致收敛.数值实验验证了本文的理论结果.     

TRAVELING WAVES CONNECTING EQUILIBRIUM AND PERIODIC ORBIT FOR DELAYED LATTICE DIFFERENTIAL EQUATION

A class of lattice with time delay and nonlocal response is considered.By transforming the lattice delay differential system into an integral equations in a Banach space,we reduces a singular perturbation problem to a regular perturbation problem.Traveling wave solution therefore is obtained by applying Liapunov-Schmidt method and the implicit function theorem.     

2.5维介质Born近似波速反演唯一性

考虑脉冲源引起的２．５维弱不均匀介质波速反演问题，利用线性化方法得到了波速的二维小扰动满足的积分方程，这是一个积分几何的问题，进而由Ｆｏｕｒｉｅｒ变换和脉冲函数的性质将此二维积分方程化为单变量的积分方程，最后用压缩映象理论证明了积分方程解的唯一性。本文给出了二给波速反演的一种新算法。同时，唯一性结果证明了已有的迭代算法的合理性。     

A perturbation-boundary integral equation method for transient dynamics of nonlinear elastic materials

This study presents a formulation for solving the transient dynamics of nonlinear elastic materials. By using a perturbation expansion to linearize the basic equations and applying the Laplace transform to the subsequent perturbation equations, the boundary value problem of the transformed equations is further reduced to various boundary integral equations. After discretization of the integral equations, these are solved numerically, completing the solution in the Laplace transform space. Performing a numerical inversion of the Laplace transform yields the solution of the problem in the time domain.     

Linearization Ill-Posedness for 2.5-D Ware Equation Inversion Model

Abstract For the weakly inhomogeneous acoustic medium in Ω={(x,y,z):z>0},we consider the inverse problemof determining the density function ρ(x,y).The inversion input for our inverse problem is the wave field givenon a line.We get an integral equation for the 2-D density perturbation from the linearization.By virtue of theintegral transform,we prove the uniqueness and the instability of the solution to the integral equation.Thedegree of ill-posedness for this problem is also given.     

Detection of Wave Front Set Perturbations via Correlation: Foundation for Wave-equation Tomography

We discuss the mathematical aspects of wave field measurements used in traveltime inversion from seismograms. The primary information about the medium is assumed to be carried by the wave front set and its perturbation with respect to a hypothetical background medium is to be estimated. By a convincing heuristics a detection procedure for this perturbation was proposed based on optimization of wave field correlations. We investigate its theoretical foundation in simple mathematical case studies using the distribution theoretic definition of oscillatory integrals.     

On the linearized theory op flow around bodies. Method of force sources

The method of force sources is proposed for solving linear problems related to the interaction between rigid bodies, and fluids, or gases. Method is based on the introduction of perturbation force sources into equation of motion of fluid media. Boundary conditions at the rigid body surface make it possible to reduce the problem of hydrodynamic reactions to an integral equation defining the function of force sources. Method is illustrated by the solution of three simple problems in the field of acoustics, and of viscous, and compressible media flow around bodies.

In the linearized theory of flow around rigid bodies, as well as in acoustics, an important part of the sound wave generation analysis concerns the determination of hydrodynamic reactions of the medium on moving, pulsating, or oscillating bodies. Such reactions make themselves felt as constant, or variable mechanical forces, such as drag and lift, or in the case of sound wave emitters, as the wave resistance. Various methods had been proposed for the computation of such forces, as for example, in the monographs [1 to 6].

Here, a different approach to the problem of determination of surface forces exerted by liquids and gases on the rigid body is proposed. By resorting to the formalism of the generalized functions it is possible to introduce into the equations of motion of fluid media a perturbation source in the form volume density of forces exercised by the body on the gas. The distribution of surface tension entering into the expression of this force is selected in such a manner as to satisfy boundary conditions at the body surface. It becomes possible with the use of this device to reduce the problem of determination of forces acting on the body surface to the solution of certain Integral equations. The proposed method is in all respects completely analogous to the well-known method of sources and sinks [1 to 1]. Both methods reduce the problem of interaction between body and gas to the solution of Integral equations. The method of sources and sinks, however, leads to an integral equation which describes the distribution of fictitious sources and sinks in the volume of the body having the density of the medium, while the method of force sources yields an integral equation which directly defines the distribution of mechanical forces over the surface of the body (*).

We may note that the method of force sources had to a certain extent been already used in papers [6 and 7] for the determination of sound radiation by means of point-force sources.     

Traveling waves connecting equilibrium and periodic orbit for reaction-diffusion equations with time delay and nonlocal response

A class of reaction-diffusion equations with time delay and nonlocal response is considered. Assuming that the corresponding reaction equations have heteroclinic orbits connecting an equilibrium point and a periodic solution, we show the existence of traveling wave solutions of large wave speed joining an equilibrium point and a periodic solution for reaction-diffusion equations. Our approach is based on a transformation of the differential equations to integral equations in a Banach space and the rigorous analysis of the property for a corresponding linear operator. Our approach eventually reduces a singular perturbation problem to a regular perturbation problem. The existence of traveling wave solution therefore is obtained by the application of Liapunov-Schmidt method and the Implicit Function Theorem.     

一类对偶积分方程组正则化为Cauchy奇异积分方程组解法

本文基于Mellin变换法求解复杂更一般形式的对偶积分方程组．通过积分变换，由实数域化成复数域上的方程组，引入未知函数的积分变换，移动积分路径，应用Cauchy积分定理，实现退耦正则化为Cauchy奇异积分方程组，由此给出一般性解，并严格证明了对偶积分方程组退耦正则化为Cauchy奇异积分方程组与原对偶积分方程组等价性，以及对偶积分方程组解的存在性和唯一性．给出的解法和理论解，作为求解复杂对偶积分方程组一种有效解法，可供求解复杂的数学、物理、力学中的混合边值问题应用．     

Anisotropic transient thermoelasticity analysis in a two-dimensional decagonal quasicrystal using meshless local Petrov–Galerkin (MLPG) method

The meshless local Petrov–Galerkin (MLPG) method is employed for anisotropic transient thermoelasticity analysis of 2D decagonal quasicrystals (QCs) subjected to transient thermal and mechanical shock loadings. The wave type model and the elasto-hydrodynamic model are applied to derive the phonon and phason governing equations, respectively. The temperature affects only the phonon field. To find the temperature distributions on the assumed 2D domain, the anisotropic heat conduction problem is solved using the MLPG method. Also, the MLPG method is successfully employed to obtain the transient behaviors of both phonon and phason displacements by solving the governing equations in local integral equations (LIEs) forms. Making use a unit step function as the test functions in the local weak-form of governing equations, we derived the local integral equations (LIEs) considered on small subdomains identical with support domains of test functions around each node. The radial basis functions are used for approximation of the spatial variation of field variables. The Laplace-transform technique is utilized to discretize the time variations.     

The limiting behaviour of solutions of the acoustic transmission problems

The boundary integral equations for the transmission problems for the scalar Helmholtz equation have the property that the dimension of the null spaces changes as the transmission parameter tends to zero in the case where the wave number is an interior eigenvalue. Therefore, the investigation of the continuous dependence of solutions of the transmission problems leads to a certain singular perturbation problem. In order to investigate this problem, we generalize a perturbation theorem for parameter dependent linear operator equations of the second kind in Banach spaces given by Kress [6, 7]. In our study we also introduce a new integral equation formulation for the transmission problem which is better situated for our purposes than the classical approach.