A multilevel correction method for Stokes eigenvalue problems and its applications

In this paper, a new multilevel correction scheme is proposed to solve Stokes eigenvalue problems by the finite element method. This new scheme contains a series of correction steps, and the accuracy of eigenpair approximation can be improved after each step. In each correction step, we only need to solve a Stokes problem on the corresponding fine finite element space and a Stokes eigenvalue problem on the coarsest finite element space. This correction scheme can improve the efficiency of solving Stokes eigenvalue problems by the finite element method. As applications of this multilevel correction method, a multigrid method and an adaptive finite element technique are introduced for Stokes eigenvalue problems. Some numerical results are given to validate our schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

The continuation of solutions for certain nonlinear differential equations and its applications

We give some sufficient conditions for the continuation of solutions for some nonlinear differential equations. As an application, we obtain a new criterion for the oscillation of solutions of the Liénard equation.     

Parametric analytical preconditioning and its applications to the reduced collocation methods

In this paper, we extend the recently developed reduced collocation method [3] to the nonlinear case, and propose two analytical preconditioning strategies. One is parameter independent and easy to implement, the other one has the traditional affinity with respect to the parameters, which allows an efficient implementation through an offline–online decomposition. Overall, preconditioning improves the quality of the error estimation uniformly on the parameter domain, and speeds up the convergence of the reduced solution to the truth approximation.     

含参变量的三阶方向牛顿法及其收敛性

通过递推关系归纳迭代公式的讨论,研究含多个未知数的非光滑方程组及其收敛性,并以此证明希尔伯特空间上的含参变量的实系数非线性方程组的三阶方向牛顿法的半局部收敛性,给出解的存在性以及先验误差界.     

Unique continuation for the two-dimensional anisotropic elasticity system and its applications to inverse problems

Under some generic assumptions we prove the unique continuation property for the two-dimensional inhomogeneous anisotropic elasticity system. Having established the unique continuation property, we then investigate the inverse problem of reconstructing the inclusion or cavity embedded in a plane elastic body with inhomogeneous anisotropic medium by infinitely many localized boundary measurements.

     

Parametric continuation of the solitary traveling pulse solution in the reaction-diffusion system using the Newton-Krylov method

The matrix-free Newton-Krylov method that uses the GMRES algorithm (an iterative algorithm for solving systems of linear algebraic equations) is used for the parametric continuation of the solitary traveling pulse solution in a three-component reaction-diffusion system. Using the results of integration on a short time interval, we replace the original system of nonlinear algebraic equations by another system that has more convenient (from the viewpoint of the spectral properties of the GMRES algorithm) Jacobi matrix. The proposed parametric continuation proved to be efficient for large-scale problems, and it made it possible to thoroughly examine the dependence of localized solutions on a parameter of the model.     

Parametric Borwein-Preiss variational principle and applications

A parametric version of the Borwein-Preiss smooth variational principle is presented, which states that under suitable assumptions on a given convex function depending on a parameter, the minimum point of a smooth convex perturbation of it depends continuously on the parameter. Some applications are given: existence of a Nash equilibrium and a solution of a variational inequality for a system of partially convex functions, perturbed by arbitrarily small smooth convex perturbations when one of the functions has a non-compact domain; a parametric version of the Kuhn-Tucker theorem which contains a parametric smooth variational principle with constraints; existence of a continuous selection of a subdifferential mapping depending on a parameter.

The tool for proving this parametric smooth variational principle is a useful lemma about continuous -minimizers of quasi-convex functions depending on a parameter, which has independent interest since it allows direct proofs of Ky Fan's minimax inequality, minimax equalities for quasi-convex functions, Sion's minimax theorem, etc.

     

Modified Rayleigh conjecture method and its applications

The Rayleigh conjecture about convergence up to the boundary of the series representing the scattered field in the exterior of an obstacle D$D$ is widely used by engineers in applications. However this conjecture is false for some obstacles. AGR introduced the Modified Rayleigh Conjecture (MRC), which is an exact mathematical result. In this paper we present the theoretical basis for the MRC method for 2D and 3D obstacle scattering problems, for static problems, and for scattering by periodic structures. We also present successful numerical algorithms based on the MRC for various scattering problems. The MRC method is easy to implement for both simple and complex geometries. It is shown to be a viable alternative for other obstacle scattering methods. Various direct and inverse scattering problems require finding global minima of functions of several variables. The Stability Index Method (SIM) combines stochastic and deterministic method to accomplish such a minimization.     

A new continuation method for complementarity problems with uniformP-functions

The complementarity problem with a nonlinear continuous mappingf from the nonnegative orthantR ₊ ⁿ ofR ⁿ intoR ⁿ can be written as the system of equationsF(x, y) = 0 and(x, y) R ₊ ²ⁿ , whereF denotes the mapping from the nonnegative orthantR ₊ ²ⁿ ofR ²ⁿ intoR ₊ ⁿ × Rⁿ defined byF(x, y) = (x ₁y₁,,x_ny_n, f₁(x) – y₁,, f_n(x) – y_n) for every(x, y) R ₊ ²ⁿ . Under the assumption thatf is a uniformP-function, this paper establishes that the mappingF is a homeomorphism ofR ₊ ²ⁿ ontoR ₊ ⁿ × Rⁿ. This result provides a theoretical basis for a new continuation method of tracing the solution curve of the one parameter family of systems of equationsF(x, y) = tF(x ⁰, y⁰) and(x, y) R ₊ ²ⁿ from an arbitrary initial point(x ⁰, y⁰) R ₊ ²ⁿ witht = 1 until the parametert attains 0. This approach is an extension of the one used in the polynomially bounded algorithm recently given by Kojima, Mizuno and Yoshise for solving linear complementarity problems with positive semi-definite matrices.     

An improved generalized F-expansion method and its applications

An improved generalized F-expansion method is proposed to seek exact solutions of nonlinear evolution equations. With the aid of symbolic computation, we choose the generalized nonlinear Schrodinger equation with a source to illustrate the validity and advantages of the proposed method. Many new and more general non-travelling wave solutions are obtained, including single and combined non-degenerate Jacobi elliptic function solutions, solitons-like solutions, trigonometric function solutions and rational solutions. The method is straightforward and concise, and it can also be applied to other nonlinear evolution equations in physics.     

Extended Jacobin elliptic function method and its applications

An extended Jacobin elliptic function method is presented for constructing exact travelling wave solutions of nonlinear partial differential equations (PDEs) in a unified way. The main idea of this method is to take full advantage of the elliptic equation that Jacobin elliptic functions satisfy and use its solutions to replace Jacobin elliptic functions in Jacobin elliptic function method. It is interesting that many other methods are special cases of our method. Some illustrative equations are investigated by this means.     

Unique continuation for first-order systems with integrable coefficients and applications to elasticity and plasticity

Let $\Omega ?{\mathbb{R}}^N$ be a Lipschitz domain and Γ be a relatively open and non-empty subset of its boundary ?Ω. We show that the solution to the linear first-order system:(1)$\mathrm{?}\zeta =G\zeta ,\zeta {|}_{\Gamma }=0,$ vanishes if $G\in {\mathsf{L}}^1(\Omega ;{\mathbb{R}}^{(N\times N)\times N})$ and $\zeta \in {\mathsf{W}}^{1,1}(\Omega ;{\mathbb{R}}^N)$. In particular, square-integrable solutions ζ of (1) with $G\in {\mathsf{L}}^1\cap {\mathsf{L}}^2(\Omega ;{\mathbb{R}}^{(N\times N)\times N})$ vanish. As a consequence, we prove that:$???:{\mathsf{C}}_{°}^{\infty }(\Omega ,\Gamma ;{\mathbb{R}}^3)\to [0,\infty ),u?6\mathrm{sym}\left(\mathrm{?}u{P}^{?1}\right)6_{{\mathsf{L}}^2(\Omega )}$ is a norm if $P\in {\mathsf{L}}^{\infty }(\Omega ;{\mathbb{R}}^{3\times 3})$ with $\mathrm{Curl}P\in {\mathsf{L}}^p(\Omega ;{\mathbb{R}}^{3\times 3})$, $\mathrm{Curl}{P}^{?1}\in {\mathsf{L}}^q(\Omega ;{\mathbb{R}}^{3\times 3})$ for some $p,q>1$ with $1/p+1/q=1$ as well as $\det P?c^{+}>0$. We also give a new and different proof for the so-called ‘infinitesimal rigid displacement lemma’ in curvilinear coordinates: Let $\Phi \in {\mathsf{H}}^1(\Omega ;{\mathbb{R}}^3)$, $\Omega ?{\mathbb{R}}^3$, satisfy $\mathrm{sym}(\mathrm{?}{\Phi }^{?}\mathrm{?}\Psi )=0$ for some $\Psi \in {\mathsf{W}}^{1,\infty }(\Omega ;{\mathbb{R}}^3)\cap {\mathsf{H}}^2(\Omega ;{\mathbb{R}}^3)$ with $\det \mathrm{?}\Psi ?c^{+}>0$. Then there exists a constant translation vector $a\in {\mathbb{R}}^3$ and a constant skew-symmetric matrix $A\in \mathfrak{so}(3)$, such that $\Phi =A\Psi +a$.     

A CONTINUATION HOMOTOPY METHOD FOR THE INVERSE PROBLEM OF OPERATOR IDENTIFICATION AND ITS APPLICATION

A new widly convergent method for solving the problem of operator kientification is illustrated.Numerical simulations are carried out to test the feasibllity and to study the general characteristics of the technique without the real measurement data.This technique is a direct application of the continuation homotopy method for solving nonlinear systems of equations.It is found that this method does give excellent results in solving the inverse problem of the elliptic differential equations.     

Variational method in wiener-hopf equations and its applications

Further discussions are made on equivalence of variational inequalities with Wiener—Hopf equations especially with Toeplitz equations. equivalent propositions are given. Starting from the results, variational method for solving Toeplitz equation is considered and several applications of Toeplitz equation on space crack and contact problems are also discussed.