四阶抛物偏微分方程的H1-Galerkin混合元方法及数值模拟

到目前为止, H¹-Galerkin 混合有限元方法研究的问题仅局限于二阶发展方程. 然而对于高阶发展方程, 特别是重要的四阶发展方程问题的研究却没有出现. 本文首次提出四阶发展方程的H¹-Galerkin 混合有限元方法, 为了给出理论分析的需要, 我们考虑四阶抛物型发展方程. 通过引进三个适当的中间辅助变量, 形成四个一阶方程组成的方程组系统, 提出四阶抛物型方程的H¹-Galerkin 混合有限元方法. 得到了一维情形下的半离散和全离散格式的最优收敛阶误差估计和多维情形的半离散格式误差估计, 并采用迭代方法证明了全离散格式的稳定性. 最后, 通过数值例子验证了提出算法的可行性. 在一维情况下我们能够同时得到未知纯量函数、一阶导数、负二阶导数和负三阶导数的最优逼近解, 这一点是以往混合元方法所不能得到的.     

四阶强阻尼波动方程的混合控制体积法

本文利用混合控制体积方法在三角网格剖分下求解四阶强阻尼波动方程.通过使用最低阶Raviart-Thomas混合有限元空间和引入迁移算子把解函数空间映射成试探函数空间,构造了半离散和全离散的混合控制体积格式,得到了最优阶误差估计.     

离散分数阶和分算子与差分算子及互逆性

对连续函数进行离散化,给出离散序列的分数阶和分算子与分数阶差分算子的解析表达式,证明了两算子满足交换律、指数率与互逆性.     

广义神经传播方程的一种修正混合有限元方法的误差分析

利用修正的H~1-Galerkin混合有限元方法研究了广义神经传播方程,论证了其半离散解的存在唯一性,得到了半离散解的最优阶误差估计,该方法的优点是不需验证LBB相容性条件.     

一类推广的离散分数阶Gronwall不等式

离散不等式,特别是离散的Gronwall不等式已被广泛应用于差分方程的研究.近年来,分数阶微分方程引起很多学者的关注.因此,利用一种新的分数阶和分的定义和不等式的方法,讨论一类更一般的离散分数阶Gronwall不等式.     

一个四阶抛物方程整体弱解的正则性

本文研究一个四阶抛物方程的非负大初值混合Dirichlet-Neumann边值问题.使用半离散化解的精细熵估计与插值技巧,得到了正则性更好的整体弱解.     

半线性Sobolev方程的H~1-Galerkin混合有限元方法

利用H~1-Galerkin混合有限元方法研究了一维半线性Sobolev方程,得到了半离散解的最优阶误差估计,优点是不需验证LBB相容性条件.     

四阶强阻尼非线性波动方程的Hermite型矩形混合有限元分析

讨论了四阶强阻尼非线性波动方程的Hermite型混合有限元方法,并证明了半离散格式下解的存在唯一性.基于该元积分恒等式结果,利用插值与Ritz投影之间的误差估计,可得到半离散格式下O(h~3)阶的超逼近性质,再借助于插值后处理技术导出整体超收敛.进而,通过构造一个新的金离散格式,得到了O(h~3+τ~2)的超逼近和超收敛结果.     

非线性抛物型偏积分微分方程的H1-Galerkin 混合有限元方法

收稿给出一类非线性抛物型偏积分微分方程的H1-Galerkin混合有限元方法.给出了一维空间的半离散、全离散格式及最优阶误差估计,并将该方法推广到二维和三维空间.     

离散分数阶神经网络的全局Mittag-Leffler稳定性

研究了一类离散分数阶神经网络的Mittag-Leffler稳定性问题.首先, 基于离散分数阶微积分理论、神经网络理论,提出了一类离散分数阶神经网络.其次,利用不等式技巧和离散Laplace变换,通过构造合适的Lyapunov函数,得到了离散分数阶神经网络全局Mittag-Leffler稳定的充分性判据.最后,通过一个数值仿真算例验证了所提出理论的有效性.     

Sharp Growth Estimates for Modified Poisson Integrals in a Half Space

For continuous boundary data, including data of polynomial growth, modified Poisson integrals are used to write solutions to the half space Dirichlet and Neumann problems in Rⁿ. Pointwise growth estimates for these integrals are given and the estimates are proved sharp in a strong sense. For decaying data, a new type of modified Poisson integral is introduced and used to develop asymptotic expansions for solutions of these half space problems.     

Dirichlet and Neumann problems for Laplace and heat equations in domains with right angles

The Dirichlet and Neumann problems are considered in the n-dimensional cube and in a right angle. The right-hand side is assumed to be bounded, and the boundary conditions are assumed to be zero. The author obtains a priori bounds for solutions in the Zygmund space, which is wider than the Lipschitz space C ^1,1 but narrower than the Hölder space C ^{1, α}, 0 < α < 1. Also, the first and second boundary-value problems are considered for the heat equation with similar conditions. It is shown that the solutions belong to the corresponding Zygmund space.     

Superconvergence of recovered gradients of discrete time/piecewise linear Galerkin approximations for linear and nonlinear parabolic problems

Superconvergent error estimates in l₂(H¹) and l∞(H¹) norms are derived for recovered gradients of finite difference in time/piecewise linear Galerkin approximations in space for linear and quasinonlinear parabolic problems in two space dimensions. The analysis extends previous results for elliptic problems to the parabolic context, and covers problems in regions with nonsmooth boundaries under certain assumptions on the regularity of the solutions. © 1994 John Wiley & Sons, Inc.     

A numerical comparison of Chebyshev methods for solving fourth order semilinear initial boundary value problems

In solving semilinear initial boundary value problems with prescribed non-periodic boundary conditions using implicit-explicit and implicit time stepping schemes, both the function and derivatives of the function may need to be computed accurately at each time step. To determine the best Chebyshev collocation method to do this, the accuracy of the real space Chebyshev differentiation, spectral space preconditioned Chebyshev tau, real space Chebyshev integration and spectral space Chebyshev integration methods are compared in the L² and W^2,2 norms when solving linear fourth order boundary value problems; and in the L^∞([0,T];L²) and L^∞([0,T];W^2,2) norms when solving initial boundary value problems. We find that the best Chebyshev method to use for high resolution computations of solutions to initial boundary value problems is the spectral space Chebyshev integration method which uses sparse matrix operations and has a computational cost comparable to Fourier spectral discretization.     

Left-Definite Regular Hamiltonian Systems

Linear Hamiltonian systems allow us to generalize, as well as consider, self-adjoint problems of any even order. Such left-definite problems are interesting, not only because of the generalization, but also because of the new intricacies they expose, some of which have made it possible to go beyond fourth order scale problems. We explore the left definite Sobolev settings for such problems, which are in general subspaces determined by boundary conditions. We show that the Hamiltonian operator remains self-adjoint, and inherits the same resolvent and spectral resolution from its original L² space when set in the left-definite Sobolev space.     

Méthode de compacité et de décomposition applications: Minimisation,convergence des martingales,lemma de Fatou multivoque

Summary New results of decomposition for bounded sequences in L _E ¹ and in the space of integrably bounded multifunctions with non empty convex weakly compact values in a Banach space E and its applications to problems of Minimization, convergence of martingales, Multivalued Fatou lemma are presented.     

Domain decomposition methods for hyperbolic problems

In this paper a method is developed for solving hyperbolic initial boundary value problems in one space dimension using domain decomposition, which can be extended to problems in several space dimensions. We minimize a functional which is the sum of squares of the L ² norms of the residuals and a term which is the sum of the squares of the L ² norms of the jumps in the function across interdomain boundaries. To make the problem well posed the interdomain boundaries are made to move back and forth at alternate time steps with sufficiently high speed. We construct parallel preconditioners and obtain error estimates for the method. The Schwarz waveform relaxation method is often employed to solve hyperbolic problems using domain decomposition but this technique faces difficulties if the system becomes characteristic at the inter-element boundaries. By making the inter-element boundaries move faster than the fastest wave speed associated with the hyperbolic system we are able to overcome this problem.     

On polynomial solvability of some problems of a vector subset choice in a Euclidean space of fixed dimension

The problems under study are connected with the choice of a vector subset from a given finite set of vectors in the Euclidean space ℝ ^k . The sum norm and averaged square of the sumnorm are considered as the target functions (to be maximized). The optimal combinatorial algorithms with time complexity O(k ² n ^2k ) are developed for these problems. Thus, the polynomial solvability of these problems is proved for k fixed.     

Optimal control problems on manifolds: a dynamic programming approach

Dynamic programming identifies the value function of continuous time optimal control with a solution to the Hamilton-Jacobi equation, appropriately defined. This relationship in turn leads to sufficient conditions of global optimality, which have been widely used to confirm the optimality of putative minimisers. In continuous time optimal control, the dynamic programming methodology has been used for problems with state space a vector space. However there are many problems of interest in which it is necessary to regard the state space as a manifold. This paper extends dynamic programming to cover problems in which the state space is a general finite-dimensional C^∞ manifold. It shows that, also in a manifold setting, we can characterise the value function of a free time optimal control problem as a unique lower semicontinuous, lower bounded, generalised solution of the Hamilton-Jacobi equation. The application of these results is illustrated by the investigation of minimum time controllers for a rigid pendulum.     

To the Sobolev embedding theorem for the limiting exponent

We establish embeddings of the Sobolev space W _p ^s and the space B _pq ^s (with the limiting exponent) in certain spaces of locally integrable functions of zero smoothness. This refines the embedding of the Sobolev space in the Lorentz and Lorentz-Zygmund spaces. Similar problems are considered for the case of irregular domains and for the potential space.