KdV方程的Chebyshev-Hermite谱配置法

针对无界区域上Korteweg.-de Vries(KdV)方程构造了时空全离散的ChebyshevHermite谱配置格式,即在空间方向上采用Hermite谱配置方法离散,时间方向上采用Chebyshev谱配置方法离散.提出了一个简单迭代算法,该算法非常适合并行计算.数值结果显示了此算法的有效性.     

路径跟踪方法求解无界非凸区域上的不动点问题

给出了求解一类无界非凸区域上不动点问题的路径跟踪方法.在适当的条件下,给出了不动点存在性的构造性证明,从而得到了路径跟踪方法的全局收敛性结果.研究结果为计算无界非凸区域上不动点问题提供了一种全局收敛性方法.     

多维区域中非线性偏微分方程的修正Laguerre谱与拟谱方法

研究多维区域中非线性偏微分方程的谱与拟谱方法．建立了修正Laguerre正交逼近与插值结果,这些结果对于建立和分析无界区域中的数值方法起着重要的作用．作为结果的一个应用,研究了二维无界区域中的Logistic方程的修正Laguerre谱格式,证明了它的稳定性和收敛性．数值试验结果表明所提出方法具有很高的精度,与理论分析结果完全吻合．     

不可压缩Navier-Stokes方程数值模拟

建立不可压缩Navier-Stokes方程的Crank-Nicolson有限差分方法,数值模拟了在初始正弦波下的二维水槽内流体受到倾斜激励时流场和涡的演变机制.数值结果表明,初始波为正弦波时,流场内出现一个单涡,单涡下沉变成了两个小涡,两个小涡消失后流场内部出现三条规则的流场带,最后这三条流场带演变成一个尖涡,尖涡在周围流体的作用下演变成一个单涡,最后单涡在自由面消失,当耗散系数和Reynolds数增大时,流场和涡演化的周期变小.     

多维区域中非线性偏微分方程的修正Uaguerre谱与拟谱方法

研究多维区域中非线性偏微分方程的谱与拟谱方法.建立了修正Laguerre正交逼近与插值结果,这些结果对于建立和分析无界区域中的数值方法起着重要的作用.作为结果的一个应用,研究了二维无界区域中的Logistic方程的修正Laguerre谱格式,证明了它的稳定性和收敛性.数值试验结果表明所提出方法具有很高的精度,与理论分析结果完全吻合.     

粘弹性二阶流体混合层流场拟序结构的数值研究

本文用拟谱方法对随时间发展的二维粘弹性二阶流体混合层流场进行了直接数据值模拟,给出在高雷诺数和低Deborah数下大涡的卷起、配对和合并等过程,通过与相同雷诺数下牛顿流体的比较,揭示了弱粘弹性对混合层中大涡拟序结构演变的影响.     

圆外区域求解Helmholtz方程

1引言许多科学和工程计算问题都可以归结为无界区域上的偏微分方程边值问题．而求解椭圆方程边值问题的常用技术是有限元方法,可是对于无界区域,在用有限元方法求解时,往往遇到困难．最简单的办法显然是直接略去区域的无界部分求解,但这样做或者导致过低的计算精度,或者要付出很高的计算代价．边界归化,即将求解偏微分方程边值问题转化为边界积分方程,是求解某些无界区域问题的强有力的手段．自70年代以来,有限元和     

螺旋型旋风分离器两相流场的数值模拟

对螺旋型旋风分离器进行了两相流场的三维数值模拟．气体流场通过求解三维N-S方程得到,湍流模型采用了雷诺应力模型．计算结果表明,旋风分离器内部的流场分为两部分:螺旋通道内比较稳定的流场和筒体中心区域的复合涡结构流场．对颗粒运动轨迹的计算表明,颗粒在入口处的初始位置对颗粒分离有比较显著的影响．同时得到了不同入口速度下颗粒的分级效率曲线,并给出了气体流量对旋风分离器性能的影响,结果显示:气体流量的增加会提高分离效率,但同时导致压力损失的急剧增加．     

半无界非线性热传导方程的Laguerre拟谱方法

以Laguerre-Gauss-Radau节点为配置点,利用拟谱方法求数值解,逼近半无界非线性热传导方程非齐次Neumann边界条件的正确解.给出算法格式和相应的数值例子,表明所提算法格式的有效性和高精度.这里所用方法也可用于求解其他非线性问题.     

FENE-P流体混合层中拟序结构的研究

采用FENE-P模型对高雷诺数下二维粘弹性混合过程中拟序结构的演变进行了数值研究.由于对聚合物应力场采用了适当的滤波措施,得到了较AH更宽参数范围内FENE-P流体混合层流场中拟序结构的运动特性,计算结果表明:在流场中加入聚合物,基波和次谐波的发展受到抑制,涡量扩散加强,减慢了配对时两涡核的旋转运动,这种影响随Weissenberg数的增大而减小,但却随参数b的增大而加强.另外,与相同溶液粘度下的牛顿流体相比,配对过程中涡量完全合并的时刻延迟了.     

The eigenvalues of Hermite and rational spectral differentiation matrices

Summary We derive expressions for the eigenvalues of spectral differentiation matrices for unbounded domains. In particular, we consider Galerkin and collocation methods based on Hermite functions as well as rational functions (a Fourier series combined with a cotangent mapping). We show that (i) first derivative matrices have purely imaginary eigenvalues and second derivative matrices have real and negative eigenvalues, (ii) for the Hermite method the eigenvalues are determined by the roots of the Hermite polynomials and for the rational method they are determined by the Laguerre polynomials, and (iii) the Hermite method has attractive stability properties in the sense of small condition numbers and spectral radii.     

Spectral methods based on Hermite functions for linear hyperbolic equations

We consider the approximation by spectral and pseudo‐spectral methods of the solution of the Cauchy problem for a scalar linear hyperbolic equation in one space dimension posed in the whole real line. To deal with the fact that the domain of the equation is unbounded, we use Hermite functions as orthogonal basis. Under certain conditions on the coefficients of the equation, we prove the spectral convergence rate of the approximate solutions for regular initial data in a weighted space related to the Hermite functions. Numerical evidence of this convergence is also presented. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

Some progress in spectral methods

In this paper,we review some results on the spectral methods.We frst consider the Jacobi spectral method and the generalized Jacobi spectral method for various problems,including degenerated and singular diferential equations.Then we present the generalized Jacobi quasi-orthogonal approximation and its applications to the spectral element methods for high order problems with mixed inhomogeneous boundary conditions.We also discuss the related spectral methods for non-rectangular domains and the irrational spectral methods for unbounded domains.Next,we consider the Hermite spectral method and the generalized Hermite spectral method with their applications.Finally,we consider the Laguerre spectral method and the generalized Laguerre spectral method for many problems defned on unbounded domains.We also present the generalized Laguerre quasi-orthogonal approximation and its applications to certain problems of non-standard type and exterior problems.     

Zur Theorie der Integraltransformationen

Let A be a selfadjoint uniformly elliptic differential operator. Let the underlying domain be bounded. Eigenvalue problems can be solved then, and an arbitrary square integrable function may be developed in a Fourier series relative to the eigenfunctions. In general elliptic differential operators have a continuous spectrum, if the underlying domain is unbounded. In this case the spectral theorem provides a representation of a given function by an integral transformation. The spectral projector can be calculated, if the outgoing and incoming solutions are known (radiation condition). Thus integral transformations may be derived very easily. Four examples will be given: the Fourier sine transform, the Lebedev transform, the transformation belonging to the Dirichlet problem of the plate equation and, finally, the Fourier transformation.     

A stabilized Hermite spectral method for second‐order differential equations in unbounded domains

A stabilized Hermite spectral method, which uses the Hermite polynomials as trial functions, is presented for the heat equation and the generalized Burgers equation in unbounded domains. In order to overcome instability that may occur in direct Hermite spectral methods, a time‐dependent scaling factor is employed in the Hermite expansions. The stability of the scheme is examined and optimal error estimates are derived. Numerical experiments are given to confirm the theoretical results.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Interaction between Kirchhoff vortices and point vortices in an ideal fluid

We consider the interaction of two vortex patches (elliptic Kirchhoff vortices) which move in an unbounded volume of an ideal incompressible fluid. A moment second-order model is used to describe the interaction. The case of integrability of a Kirchhoff vortex and a point vortex by the variable separation method is qualitatively analyzed. A new case of integrability of two Kirchhoff vortices is found. A reduced form of equations for two Kirchhoff vortices is proposed and used to analyze their regular and chaotic behavior.     

Solution to the boundary value problem for an arbitrary elliptic operator subject to a radiation condition

It is shown that the generalized Fourier transform can be extended to an arbitrary elliptic operator in a cylindrical domain with a Robin boundary condition. In this case, the existence of the Fourier image is a completely correct radiation condition determining a solution to the problem that is a superposition of waves traveling away from the source.     

Bingham Flows in Periodic Domains of Infinite Length

The Bingham fluid model has been successfully used in modeling a large class of non-Newtonian fluids. In this paper, the authors extend to the case of Bingham fluids the results previously obtained by Chipot and Mardare, who studied the asymptotics of the Stokes flow in a cylindrical domain that becomes unbounded in one direction, and prove the convergence of the solution to the Bingham problem in a finite periodic domain, to the solution of the Bingham problem in the infinite periodic domain, as the length of the finite domain goes to infinity. As a consequence of this convergence, the existence of a solution to a Bingham problem in the infinite periodic domain is obtained, and the uniqueness of the velocity field for this problem is also shown. Finally, they show that the error in approximating the velocity field in the infinite domain with the velocity in a periodic domain of length 2? has a polynomial decay in ?, unlike in the Stokes case(see [Chipot,M. and Mardare, S., Asymptotic behaviour of the Stokes problem in cylinders becoming unbounded in one direction, Journal de Math′ematiques Pures et Appliqu′ees, 90(2), 2008,133–159]) where it has an exponential decay. This is in itself an important result for the numerical simulations of non-Newtonian flows in long tubes.     

TWO PROBLEMS OF HERMITE ELLIPTIC EQUATIONS

In this article, the author investigates some Hermite elliptic equations in a modified Sobolev space introduced by X. Ding [2]. First, the author shows the existence of a ground state solution of semilinear Hermite elliptic equation. Second, the author studies the eigenvalue problem of linear Hermite elliptic equation in a bounded or unbounded domain.     

A global existence result for a spatially extended 3D Navier-Stokes problem with non-small initial data

We consider the Navier-Stokes equations in a two- or three-dimensional unbounded cylindrical domain. The existence and uniqueness of solutions is discussed in the space of uniformly local square integrable functions. We show for small initial data and small forcing term that the solutions exist globally in time. This result is extended to a non-small data result in the sense that the high frequency modes of the initial conditions and of the forcing terms are allowed to be large. Moreover, we show the existence of a local attractor for this 3D Navier-Stokes problem in an unbounded domain. In contrast to previous results the spaces used are no Hilbert spaces, and secondly we have a linear operator possessing continuous spectrum without spectral gap. Received March 1999