Search

Current state-of-the-art methods for solving discrete optimization problems are usually restricted to convex settings. In this paper, we propose a general approach based on cutting planes for solving nonlinear, possibly nonconvex, binary optimization problems. We provide a rigorous convergence analysis that quantifies the number of iterations required under different conditions. This is different to most other work in discrete optimization where only finite convergence is proved. Moreover, using tools from variational analysis, we provide necessary and sufficient dual optimality conditions.

This paper proposes a new gradient-based optimization approach for designing optimal feedback kernels for parabolic distributed parameter systems with boundary control. Unlike traditional kernel optimization methods for parabolic systems, our new method does not require solving non-standard Riccati-type or Klein-Gorden-type partial differential equations (PDEs). Instead, the feedback kernel is parameterized as a second-order polynomial whose coefficients are decision variables to be tuned via gradient-based dynamic optimization, where the gradient of the system cost functional (which penalizes both kernel and output magnitude) is computed by solving a so-called costate PDE instandard form. Special constraints are imposed on the kernel coefficients to ensure that, under mild conditions, the optimized kernel yields closed-loop stability. Numerical simulations demonstrate the effectiveness of the proposed approach., Comment: 26 pages, 17 figures

Most current prevalent iterative methods can be classified into the so-called extended Krylov subspace methods, a class of iterative methods which do not fall into this category are also proposed in this paper. Comparing with traditional Krylov subspace methods which always depend on the matrix-vector multiplication with a fixed matrix, the newly introduced methods(the so-called (progressively) accumulated projection methods, or AP (PAP) for short) use a projection matrix which varies in every iteration to form a subspace from which an approximate solution is sought. More importantly an accelerative approach(called APAP) is introduced to improve the convergence of PAP method. Numerical experiments demonstrate some surprisingly improved convergence behavior. Comparison between benchmark extended Krylov subspace methods(Block Jacobi and GMRES) are made and one can also see remarkable advantage of APAP in some examples. APAP is also used to solve systems with extremely ill-conditioned coefficient matrix (the Hilbert matrix) and numerical experiments shows that it can bring very satisfactory results even when the size of system is up to a few thousands.

This paper introduces DD calculus and describes the basic calculus concepts of derivative and integral in a direct and non-traditional way, without limit definition: Derivative is computed from the point-slope equation of a tangent line and integral is defined as the height increment of a curve. This direct approach to calculus has three distinct features: (i) it defines derivative and (definite) integral without using limits, (ii) it defines derivative and antiderivative simultaneously via a derivative-antiderivative (DA) pair, and (iii) it posits the fundamental theorem of calculus as a natural corollary of the definitions of derivative and integral. The first D in DD calculus attributes to Descartes for his method of tangents and the second D to DA-pair. The DD calculus, or simply direct calculus, makes many traditional notations and procedures unnecessary, a plus when introducing calculus to the non-mathematics majors. It has few intermediate procedures, which can help dispel the mystery of calculus as perceived by the general public. The materials in this paper are intended for use in a two-hour introductory lecture on calculus., Comment: 19 pages, 11 figures

In this paper we will prove that the vorticity belongs to L1(0; T ; L2(R3)) for the Cauchy problem of 3D incompressible Navier-Stokes equation, then the existence of a global smooth solution is obtained. Our approach is to construct a set of auxiliary problems to approximate the original one of vorticity equation., Comment: 34 pages. arXiv admin note: substantial text overlap with arXiv:1308.2297

In this paper we will prove that the vorticity belongs to L1(0; T ; L2(\Omega)) for 3D incompressible Navier-Stokes equation with periodic initial-boundary value conditions, then the existence of a global smooth solution is obtained. Our approach is to construct a set of auxiliary problems to approximate the original one of vorticity equation., Comment: modification

In this paper, we present a parallel scheme to solve the population balance equations based on the method of characteristics and the finite element discretization. The application of the method of characteristics transform the higher dimensional population balance equation into a series of lower dimensional convection-diffusion-reaction equations which can be solved in a parallel way.Some numerical results are presented to show the accuracy and efficiency., Comment: 10 pages, 0 figure

The aim of the paper is to introduce a new systematic method that can produce lower bounds for eigenvalues. The main idea is to use nonconforming finite element methods. The general conclusion herein is that if local approximation properties of nonconforming finite element spaces $V_h$ are better than global continuity properties of $V_h$, corresponding methods will produce lower bounds for eigenvalues. More precisely, under three conditions on continuity and approximation properties of nonconforming finite element spaces we first show abstract error estimates of approximate eigenvalues and eigenfunctions. Subsequently, we propose one more condition and prove that it is sufficient to guarantee nonconforming finite element methods to produce lower bounds for eigenvalues of symmetric elliptic operators. As one application, we show that this condition hold for most nonconforming elements in literature. As another important application, this condition provides a guidance to modify known nonconforming elements in literature and to propose new nonconforming elements. In fact, we enrich locally the Crouzeix-Raviart element such that the new element satisfies the condition; we propose a new nonconforming element for second order elliptic operators and prove that it will yield lower bounds for eigenvalues. Finally, we prove the saturation condition for most nonconforming elements., Comment: 24 pages

This article is devoted to computing the lower and upper bounds of the Laplace eigenvalue problem. By using the special nonconforming finite elements, i.e., enriched Crouzeix-Raviart element and extension $Q_1^{\rm rot}$, we get the lower bound of the eigenvalue. Additionally, we also use conforming finite elements to do the postprocessing to get the upper bound of the eigenvalue. The postprocessing method need only to solve the corresponding source problems and a small eigenvalue problem if higher order postprocessing method is implemented. Thus, we can obtain the lower and upper bounds of the eigenvalues simultaneously by solving eigenvalue problem only once. Some numerical results are also presented to validate our theoretical analysis., Comment: 19 pages, 4 figures

In this paper, a new type of multi-level correction scheme is proposed for solving eigenvalue problems by finite element method. With this new scheme, the accuracy of eigenpair approximations can be improved after each correction step which only needs to solve a source problem on finer finite element space and an eigenvalue problem on the coarsest finite element space. This correction scheme can improve the efficiency of solving eigenvalue problems by finite element method., Comment: 16 pages, 5 figures

Assume that $V_h$ is a space of piecewise polynomials of degree less than $r\geq 1$ on a family of quasi-uniform triangulation of size $h$. Then the following well-known upper bound holds for a sufficiently smooth function $u$ and $p\in [1, \infty]$ $$ \inf_{v_h\in V_h}\|u-v_h\|_{j,p,\Omega,h} \le C h^{r-j} |u|_{r,p,\Omega},\quad 0\le j\le r. $$ In this paper, we prove that, roughly speaking, if $u\not\in V_h$, the above estimate is sharp. Namely, $$ \inf_{v_h\in V_h}\|u-v_h\|_{j,p,\Omega,h} \ge c h^{r-j},\quad 0\le j\le r, \ \ 1\leq p\leq \infty, $$ for some $c>0$. The above result is further extended to various situations including more general Sobolev space norms, general shape regular grids and many different types of finite element spaces. As an application, the sharpness of finite element approximation of elliptic problems and the corresponding eigenvalue problems is established., Comment: 15