AN IMPROVED NON-TRADITIONAL FINITE ELEMENT FORMULATION FOR SOLVING THE ELLIPTIC INTERFACE PROBLEMS

We propose a non-traditional finite element method with non-body-fitting grids to solve the matrix coefficient elliptic equations with sharp-edged interfaces. All possible situations that the interface cuts the grid are considered. Both Diriehlet and Neumann boundary conditions are discussed. The coefficient matrix data can be given only on the grids, rather than an analytical function. Extensive numerical experiments show that this method is second order accurate in the L∞ norm.     

ALTERNATELY LINEARIZED IMPLICIT ITERATION METHODS FOR SOLVING QUADRATIC MATRIX EQUATIONS

A numerical solution of the quadratic matrix equations associated with a nonsingular M-matrix by using the alternately linearized implicit iteration method is considered. An iteration method for computing a nonsingular M-matrix solution of the quadratic matrix equations is developed, and its corresponding theory is given. Some numerical examples are provided to show the efficiency of the new method.     

AN IMPROVED ERROR ANALYSIS FOR FINITE ELEMENT APPROXIMATION OF BIOLUMINESCENCE TOMOGRAPHY

This paper is concerned with an ill-posed problem which results from the area of molecular imaging and is known as BLT problem. Using Tikhonov regularization technique, a quadratic optimization problem can be formulated. We provide an improved error estimate for the finite element approximation of the regularized optimization problem. Some numerical examples are presented to demonstrate our theoretical results.     

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.     

FIXED-POINT CONTINUATION APPLIED TO COMPRESSED SENSING： IMPLEMENTATION AND NUMERICAL EXPERIMENTS

Fixed-point continuation （FPC） is an approach, based on operator-splitting and continuation, for solving minimization problems with l1-regularization：min ｜｜x｜｜1＋uf（x）.We investigate the application of this algorithm to compressed sensing signal recovery, in which f（x） = 1/2｜｜Ax-b｜｜2M,A∈m×n and m≤n. In particular, we extend the original algorithm to obtain better practical results, derive appropriate choices for M and u under a given measurement model, and present numerical results for a variety of compressed sensing problems. The numerical results show that the performance of our algorithm compares favorably with that of several recently proposed algorithms.     

MULTIGRID METHOD FOR A MODIFIED CURVATURE DRIVEN DIFFUSION MODEL FOR IMAGE INPAINTING

Digital inpainting is a fundamental problem in image processing and many variational models for this problem have appeared recently in the literature. Among them are the very successfully Total Variation （TV） model [11] designed for local inpainting and its improved version for large scale inpainting： the Curvature-Driven Diffusion （CDD） model [10]. For the above two models, their associated Euler Lagrange equations are highly nonlinear partial differential equations. For the TV model there exists a relatively fast and easy to implement fixed point method, so adapting the multigrid method of [24] to here is immediate. For the CDD model however, so far only the well known but usually very slow explicit time marching method has been reported and we explain why the implementation of a fixed point method for the CDD model is not straightforward. Consequently the multigrid method as in [Savage and Chen, Int. J. Comput. Math., 82 （2005）, pp. 1001-1015] will not work here. This fact represents a strong limitation to the range of applications of this model since usually fast solutions are expected. In this paper, we introduce a modification designed to enable a fixed point method to work and to preserve the features of the original CDD model. As a result, a fast and efficient multigrid method is developed for the modified model. Numerical experiments are presented to show the very good performance of the fast algorithm.     

Walsh-Type Wavelet Packet Expansions

We consider a family of basic nonstationary wavelet packets generated using the Haar filters except for a finite number of scales where we allow the use of arbitrary filters. Such a system, which we call a system of Walsh-type wavelet packets, can be considered as a smooth generalization of the Walsh functions. We show that the basic Walsh-type wavelet packets share a number of metric properties with the Walsh system. We prove that the system constitutes a Schauder basis for L^p( ), 1<p<∞, and we construct an explicit function in L¹( ) for which the expansion fails. Then we prove that expansions of L^p( )-functions, 1<p<∞, in the Walsh-type wavelet packets converge pointwise a.e. Finally, we prove that the analogous results are true for periodic Walsh-type wavelet packets in L^p[0,1).     

Cartesian closed categories of FƵ-domains

A subset system Z assigns to each partially ordered set P a certain collection Z（P） of subsets. In this paper, a new kind of subset systems called directable subset systems is introduced. For a directable subset system Z, the concepts of FZ-way-below relation and FZ-domain are introduced. The well-known Scott topology is naturally generalized to the Z-level and the resulting topology is called FZ-Scott topology, and the continuous functions with respect to this topology are characterized by preserving the suprema of directed Z-sets. Then, we mainly consider a generalization of the cartesian closedness of the categories DCPO of directed complete posets, BF of bifinite domains and FS of FS-domains to the Z-level. Corresponding to them, it is proved that, for a suitable subset system Z, the categories FZCPO of Z-complete posets, FSFZ of finitely separated FZ-domains and BFFZ of bifinite FZ-domains are all cartesian closed. Some examples of these categories are given.     

A NEW PRECONDITIONING STRATEGY FOR SOLVING A CLASS OF TIME-DEPENDENT PDE-CONSTRAINED OPTIMIZATION PROBLEMS

In this paper, by exploiting the special block and sparse structure of the coefficient matrix, we present a new preconditioning strategy for solving large sparse linear systems arising in the time-dependent distributed control problem involving the heat equation with two different functions. First a natural order-reduction is performed, and then the reduced- order linear system of equations is solved by the preconditioned MINRES algorithm with a new preconditioning techniques. The spectral properties of the preconditioned matrix are analyzed. Numerical results demonstrate that the preconditioning strategy for solving the large sparse systems discretized from the time-dependent problems is more effective for a wide range of mesh sizes and the value of the regularization parameter.     

Bifurcation in a Class of Planar Piecewise Smooth Systems with 3-parameters

This paper is concerned with the bifurcation properties on the line of discontinuity of planar piecewise smooth systems. The existence of equilibria and periodic solutions with sliding motion in a class of planar piecewise smooth systems with 3-parameters is investigated in this paper using the theory of differential inclu-sion and tools of Poincar′e maps.     

SPECIAL DYNAMIC BEHAVIORS OF A TEMPORAL CHAOTIC SYSTEM

When dynamic behaviors of temporal chaotic system are analyzed, we find that a temporal chaotic system has not only generic dynamic behaviors of chaotic reflection, but also has phenomena influencing two chaotic attractors by original values. Along with the system parameters changing to certain value, the system will appear a break in chaotic region, and jump to another orbit of attractors. When it is opposite that the system parameters change direction, the temporal chaotic system appears complicated chaotic behaviors.     

Almost Everywhere Convergence of Sequences of Cesàro and Riesz Means of Integrable Functions with Respect to the Multidimensional Walsh System

The aim of this paper is to prove the a.e.convergence of sequences of the Cesa`ro and Riesz means of the Walsh–Fourier series of d variable integrable functions.That is,let a=(a1,...,ad):N→Nd(d∈P)such that aj(n+1)≥δsupk≤n aj(k)(j=1,...,d,n∈N)for someδ0 and a1(+∞)=···=ad(+∞)=+∞.Then,for each integrable function f∈L1(Id),we have the a.e.relation for the Ces`aro means limn→∞σαa(n)f=f and for the Riesz means limn→∞σα,γa(n)f=f for any 0αj≤1≤γj(j=1,...,d).A straightforward consequence of our result is the so-called cone restricted a.e.convergence of the multidimensional Ces`aro and Riesz means of integrable functions,which was proved earlier by Weisz.     

Bifurcations of limit cycles from a quintic Hamiltonian system with a heteroclinic cycle

In this paper,we consider Li′enard systems of the form dx/dt=y,dy/dt=x+bx3-x5+ε(α+βx2+γx4)y,where b∈R,0|ε|1,(α,β,γ)∈D∈R3 and D is bounded.We prove that for |b|1(b0) the least upper bound of the number of isolated zeros of the related Abelian integrals I(h)=∮Γh(α+βx2+γx4)ydx is 2(counting the multiplicity) and this upper bound is a sharp one.     

SUPERCONVERGENCE ANALYSIS OF FINITE ELEMENT METHODS FOR OPTIMAL CONTROL PROBLEMS OF THE STATIONARY BENARD TYPE

In this paper, we consider the finite element approximation of the distributed optimal control problems of the stationary Benard type under the pointwise control constraint. The states and the co-states are approximated by polynomial functions of lowest-order mixed finite element space or piecewise linear functions and the control is approximated by piecewise constant functions. We give the superconvergence analysis for the control; it is proved that the approximation has a second-order rate of convergence. We further give the superconvergence analysis for the states and the co-states. Then we derive error estimates in L^∞-norm and optimal error estimates in L^2-norm.     

GLOBAL ASYMPTOTICS TOWARD THE REREFACTION WAVES FOR A PARABOLIC-ELLIPTIC SYSTEM RELATED TO THE CAMASSA-HOLM SHALLOW WATER EQUATION

This article is concerned with the global existence and large time behavior of solutions to the Cauchy problem for a parabolic-elliptic system related to the Camassa-Holm shallow water equation {ut＋（u^2/2）x＋px=εuxx, t〉0,x∈R, -αPxx＋P=f（u）＋α/2ux^2-1/2u^2, t〉0,x∈R, （E） with the initial data u（0,x）=u0（x）→u±, as x→±∞ （I） Here, u_ 〈 u＋ are two constants and f（u） is a sufficiently smooth function satisfying f＂ （u） 〉 0 for all u under consideration. Main aim of this article is to study the relation between solutions to the above Cauchy problem and those to the Riemann problem of the following nonlinear conservation law It is well known that if u_ 〈 u＋, the above Riemann problem admits a unique global entropy solution u^R（x/t） u^R（x/t）={u_,（f′）^-1（x/t）,u＋, x≤f′（u_）t, f′（u_）t≤x≤f′（u＋）t, x≥f′（u＋）t. Let U（t, x） be the smooth approximation of the rarefaction wave profile constructed similar to that of [21, 22, 23], we show that if u0（x） - U（0,x） ∈ H^1（R） and u_ 〈 u＋, the above Cauchy problem （E） and （I） admits a unique global classical solution u（t, x） which tends to the rarefaction wave u^R（x/t） as → ＋∞ in the maximum norm. The proof is given by an elementary energy method.     

T＊-extension of Lie Sup ertriple Systems

In this article,we study the Lie supertriple system(LSTS)T over a field K admitting a nondegenerate invariant supersymmetric bilinear form(call such a T metrisable).We give the definition of T_ω~*-extension of an LSTS T,prove a necessary and sufficient condition for a metrised LSTS(T,φ)to be isometric to a T~*-extension of some LSTS,and determine when two T~*-extensions of an LSTS are"same",i.e.,they are equivalent or isometrically equivalent.     

EXISTENCE OF PERIODIC SOLUTIONS TO A SYSTEM WITH FUNCTIONAL RESPONSE ON TIME SCALES

This paper investigates the existence of periodic solutions of a three-species food-chain diffusive system with Beddington-DeAngelis functional responses and time delays in a two-patch environment on time scales. By using a continuation theorem based on coincidence degree theory, we obtain sufficient criteria for the existence of periodic solutions for the system. Moreover, when the time scale T is chosen as R or Z, the existence of the periodic solutions of the corresponding continuous and discrete models follows. Therefore, the method is unified to provide the existence of the desired solutions for continuous differential equations and discrete difference equations.     

Blow-up vs.Global Finiteness for an Evolution p-Laplace System with Nonlinear Boundary Conditions

In this paper, the authors consider the positive solutions of the system of the evolution p-Laplacian equations
with nonlinear boundary conditions
and the initial data （u0, v0）, where Ω is a bounded domain in Rn with smooth boundary δΩ, p 〉 2, h（·,·） and s（·,· ） are positive C1 functions, nondecreasing in each variable. The authors find conditions on the functions f, g, h, s that prove the global existence or finite time blow-up of positive solutions for every （u0, v0）.     

ON THE FINITE ELEMENT APPROXIMATION OF SYSTEMS OF REACTION-DIFFUSION EQUATIONS BY p/hp METHODS

We consider the approximation of systems of reaction-diffusion equations, with the finite element method. The highest derivative in each equation is multiplied by a parameter ε∈ （0, 1], and as ε → 0 the solution of the system will contain boundary layers. We extend the analysis of the corresponding scalar problem from [Melenk, IMA J. Numer. Anal. 17（1997）, pp. 577-601], to construct a finite element scheme which includes elements of size O（εp） near the boundary, where p is the degree of the approximating polynomials. We show that, under the assumption of analytic input data, the method yields exponential rates of convergence, independently of ε, when the error is measured in the energy norm associated with the problem. Numerical computations supporting the theory are also presented, which also show that the method yields robust exponential convergence rates when the error in the maximum norm is used.     

Analytic regularity for an operator with Treves curves

We study a partial differential operator with analytic coefficients, which is of the form “sum of squares”. is hypoelliptic on any open subset of , yet possesses the following properties: (1) is not analytic hypoelliptic on any open subset of that contains 0. (2) If u is any distribution defined near with the property that is analytic near 0, then u must be analytic near 0. (3) The point 0 lies on the projection of an infinite number of Treves curves (bicharacteristics).These results are consistent with the Treves conjectures. However, it follows that the analog of Treves conjecture, in the sense of germs, is false.As far as we know, is the first example of a “sum of squares” operator which is not analytic hypoelliptic in the usual sense, yet is analytic hypoelliptic in the sense of germs.