STABILITY FOR IMPOSING ABSORBING BOUNDARY CONDITIONS IN THE FINITE ELEMENT SIMULATION OF ACOUSTIC WAVE PROPAGATION＊

It is well-known that artificial boundary conditions are crucial for the efficient and accurate computations of wavefields on unbounded domains. In this paper, we investigate stability analysis for the wave equation coupled with the first and the second order absorbing boundary conditions. The computational scheme is also developed. The approach allows the absorbing boundary conditions to be naturally imposed, which makes it easier for us to construct high order schemes for the absorbing boundary conditions. A thirdorder Lagrange finite element method with mass lumping is applied to obtain the spatial discretization of the wave equation. The resulting scheme is stable and is very efficient since no matrix inversion is needed at each time step. Moreover, we have shown both abstract and explicit conditional stability results for the fully-discrete schemes. The results are helpful for designing computational parameters in computations. Numerical computations are illustrated to show the efficiency and accuracy of our method. In particular, essentially no boundary reflection is seen at the artificial boundaries.     

On Second-order Sufficient Conditions in Constrained Nonsmooth Optimization

In this paper, we establish a second-order sufficient condition for constrained optimization problems of a class of so called t-stable functions in terms of the first-order and the second-order Dini type directional derivatives. The result extends the corresponding result of [D. Bednarik and K. Pastor, Math. Program. Ser. A, 113（2008）, 283-298] to constrained optimization problems.     

ON A FUNCTIONAL EQUATION

In this article, the author derives a functional equation η（s）=[（π/4）^s-1/2√2/π Г（1-s）sin（πs/2）]η（1-s） （1） of the analytic function η（s） which is defined by η（s）=1^-s-3^-s-5^-s＋7^-s＋… （2） for complex variable s with Re s 〉 1, and is defined by analytic continuation for other values of s. The author proves （1） by Ramanujan identity （see [1], [3]）. Her method provides a new derivation of the functional equation of Riemann zeta function by using Poisson summation formula.     

Schur Convexity for Two Classes of Symmetric Functions and Their Applications

For x =(x1, x2, ···, xn) ∈ Rn+∪ Rn-, the symmetric functions Fn(x, r) and Gn(x, r) are defined by r1 + xFij n(x, r) = Fn(x1, x2, ···, xn; r) =x1≤iij1i2···ir ≤n j=1and r1- xGij n(x, r) = Gn(x1, x2, ···, xn; r) =,x1≤i1i2···ir ≤n j=1ij respectively, where r = 1, 2, ···, n, and i1, i2, ···, in are positive integers. In this paper,the Schur convexity of Fn(x, r) and Gn(x, r) are discussed. As applications, by a bijective transformation of independent variable for a Schur convex function, the authors obtain Schur convexity for some other symmetric functions, which subsumes the main results in recent literature; and by use of the theory of majorization establish some inequalities. In particular, the authors derive from the results of this paper the Weierstrass inequalities and the Ky Fan's inequality, and give a generalization of Safta's conjecture in the n-dimensional space and others.     

A NOTE ON THE WALSH-FEJER MEANS

The main aim of this paper is to prove that the maximal operator σ# is not bounded from the martingale Hardy space Hp （G） to the martingale Hardy space Hp （G） for 0〈p≤1.     

Characterization of arithmetic functions that preserve the sum-of-squares operation

We characterize all functions f :N → C such thatf(m2+n2)=f(m)2+f(n)2for all m,n∈ N.It turns out that all such functions can be grouped into three families,namelyf≡0,f(n)=±n(subject to some restrictions on when the choice of the sign is possible) and f(n)=±1/2(again subject to some restrictions on when the choice of the sign is possible).     

WEIGHTED A PRIORI ESTIMATES FOR SOLUTION OF （-△）^mu = f WITH HOMOGENEOUS DIRICHLET CONDITIONS

Let u be a weak solution of （-△）mu = f with Dirichlet boundary conditions in a smooth bounded domain Ω C Rn. Then, the main goal of this paper is to prove the following a priori estimate：｜｜u｜｜w2m/ω·p（Ω）≤C｜｜f｜｜L^pω（Ω）,where ω is a weight in the Muckenhoupt class Ap.     

Quadratic Lyapunov Function and Exponential Dichotomy on Time Scales

In this paper, we study the relationship between exponential dichotomy and quadratic Lyapunov function for the linear equation x△=A（t）x on time scales. Moreover, for the nonlinear perturbed equation x△= A（t）x ＋ f（t, x） we give the instability of the zero solution when f is sufficiently small.     

Atomic Decompositions of Triebel-Lizorkin Spaces with Local Weights and Applications

In this paper,the authors characterize the inhomogeneous Triebel-Lizorkin spaces Fp,q s,w(Rn)with local weight w by using the Lusin-area functions for the full ranges of the indices,and then establish their atomic decompositions for s ∈ R,p ∈(0,1] and q ∈ [p,∞).The novelty is that the weight w here satisfies the classical Muckenhoupt condition only on balls with their radii in(0,1].Finite atomic decompositions for smooth functions in Fp,q s,w(Rn)are also obtained,which further implies that a(sub)linear operator that maps smooth atoms of Fp,q s,w(Rn)uniformly into a bounded set of a(quasi-)Banach space is extended to a bounded operator on the whole Fp,q s,w(Rn).As an application,the boundedness of the local Riesz operator on the space Fp,q s,w(Rn)is obtained.     

Dunkl’s theory and best approximation by entire functions of exponential type in L 2-metric with power weight

In this paper, we study the sharp Jackson inequality for the best approximation of f ∈L2,κ(Rd) by a subspace E2κ(σ)(SE2κ(σ)), which is a subspace of entire functions of exponential type(spherical exponential type) at most σ. Here L2,κ(Rd) denotes the space of all d-variate functions f endowed with the L2-norm with the weight v2κ(x) =ξ∈R, which is defined by a positive+|(ξ, x)|κ(ξ)subsystem R+ of a finite root system RRdand a function κ(ξ) : R → R+ invariant under the reflection group G(R) generated by R. In the case G(R) = Zd2, we get some exact results. Moreover,the deviation of best approximation by the subspace E2κ(σ)(SE2κ(σ)) of some class of the smooth functions in the space L2,κ(Rd) is obtained.     

Sto chastic Nonlinear Beam Equations with Levy Jump

In this paper, we study stochastic nonlinear beam equations with Levy jump, and use Lyapunov functions to prove existence of global mild solutions and asymptotic stability of the zero solution.     

Weighted Endpoint Estimates for Multilinear θ-type Calderon-Zygmund Operator

In this paper, weighted endpoint boundedness of multilinear θ-type Calderon-Zygmund operator is obtained.     

SOLVING A CLASS OF INVERSE QP PROBLEMS BY A SMOOTHING NEWTON METHOD

We consider an inverse quadratic programming （IQP） problem in which the parameters in the objective function of a given quadratic programming （QP） problem are adjusted as little as possible so that a known feasible solution becomes the optimal one. This problem can be formulated as a minimization problem with a positive semidefinite cone constraint and its dual （denoted IQD（A, b）） is a semismoothly differentiable （SC^1） convex programming problem with fewer variables than the original one. In this paper a smoothing Newton method is used for getting a Karush-Kuhn-Tucker point of IQD（A, b）. The proposed method needs to solve only one linear system per iteration and achieves quadratic convergence. Numerical experiments are reported to show that the smoothing Newton method is effective for solving this class of inverse quadratic programming problems.     

A Nine-modes Truncation of the Plane Incompressible Navier-Stokes Equations

In this paper a nine-modes truncation of Navier-Stokes equations for a two-dimensional incompressible fluid on a torus is obtained. The stationary solutions, the existence of attractor and the global stability of the equations are firmly proved. What is more, that the force f acts on the mode ks and k7 respectively produces two systems, which lead to a much richer and varied phenomenon. Numerical simulation is given at last, which shows a stochastic behavior approached through an involved sequence of bifurcations.     

A NEW TRUST-REGION ALGORITHM FOR NONLINEAR CONSTRAINED OPTIMIZATION

We propose a new trust region algorithm for nonlinear constrained optimization problems. In each iteration of our algorithm, the trial step is computed by minimizing a quadratic approximation to the augmented Lagrange function in the trust region. The augmented Lagrange function is also used as a merit function to decide whether the trial step should be accepted. Our method extends the traditional trust region approach by combining a filter technique into the rules for accepting trial steps so that a trial step could still be accepted even when it is rejected by the traditional rule based on merit function reduction. An estimate of the Lagrange multiplier is updated at each iteration, and the penalty parameter is updated to force sufficient reduction in the norm of the constraint violations. Active set technique is used to handle the inequality constraints. Numerical results for a set of constrained problems from the CUTEr collection are also reported.     

WEIGHTED L^p BOUNDEDNESS FOR MULTILINEAR FRACTIONAL INTEGRAL ON PRODUCT SPACES

For 0 〈 α 〈 mn and nonnegative integers n ≥ 2, m≥ 1, the multilinear fractional integral is defined bywhere →y= （y1, Y2,…, ym） and 7 denotes the m-tuple （f1, f2,…, fm）. In this note, the one- weighted and two-weighted boundedness on Lp （JRn） space for multilinear fractional integral operator I（am） and the fractional multi-sublinear maximal operator Mα（m） are established re- spectively. The authors also obtain two-weighted weak type estimate for the operator Mα（m）.     

ESTIMATIONS OF THE CONSTANTS IN INVERSE INEQUALITIES FOR FINITE ELEMENT FUNCTIONS

In this paper, we estimate the constants in the inverse inequalities for the finite ele- ment functions. Furthermore, we obtain the least upper bounds of the constants in inverse inequalities for the low-order finite element functions. Such explicit estimates of the con- stants can be used as computable error bounds for the finite element method.     

ON SMOOTH LU DECOMPOSITIONS WITH APPLICATIONS TO SOLUTIONS OF NONLINEAR EIGENVALUE PROBLEMS

We study the smooth LU decomposition of a given analytic functional A-matrix A（A） and its block-analogue. Sufficient conditions for the existence of such matrix decompositions are given, some differentiability about certain elements arising from them are proved, and several explicit expressions for derivatives of the specified elements are provided. By using these smooth LU decompositions, we propose two numerical methods for computing multiple nonlinear eigenvalues of A（A）, and establish their locally quadratic convergence properties. Several numerical examples are provided to show the feasibility and effectiveness of these new methods.     

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.