In this paper,an augmented two-scale finite element method is proposed for a class of linear and nonlinear eigenvalue problems on tensor-product *** a correction step,the augmented two-scale finite element solution is...
详细信息
In this paper,an augmented two-scale finite element method is proposed for a class of linear and nonlinear eigenvalue problems on tensor-product *** a correction step,the augmented two-scale finite element solution is obtained by solving an eigenvalue problem on a low-dimensional augmented *** analysis and numerical experiments show that the augmented two-scale finite element solution achieves the same order of accuracy as the standard finite element solution on a fine grid,but the computational cost required by the former solution is much lower than that demanded by the *** augmented two-scale finite element method also improves the approximation accuracy of eigenfunctions in the L^(2)(Ω)norm compared with the two-scale finite element method.
作者:
Xianmin XuLSEC
Institute of Computational Mathematics and Scientific/Engineering ComputingNCMISAMSSChinese Academy of SciencesBeijing 100190China School of Mathematical Sciences
University of Chinese Academy of SciencesBeijing 100049China
By using the Onsager principle as an approximation tool,we give a novel derivation for the moving finite element method for gradient flow *** show that the discretized problem has the same energy dissipation structure...
详细信息
By using the Onsager principle as an approximation tool,we give a novel derivation for the moving finite element method for gradient flow *** show that the discretized problem has the same energy dissipation structure as the continuous *** enables us to do numerical analysis for the stationary solution of a nonlinear reaction diffusion equation using the approximation theory of free-knot piecewise *** show that under certain conditions the solution obtained by the moving finite element method converges to a local minimizer of the total energy when time goes to *** global minimizer,once it is detected by the discrete scheme,approximates the continuous stationary solution in optimal *** examples for a linear diffusion equation and a nonlinear Allen-Cahn equation are given to verify the analytical results.
The positive definiteness of real quadratic forms with convolution structures plays an important rolein stability analysis for time-stepping schemes for nonlocal operators. In this work, we present a novel analysistoo...
详细信息
The positive definiteness of real quadratic forms with convolution structures plays an important rolein stability analysis for time-stepping schemes for nonlocal operators. In this work, we present a novel analysistool to handle discrete convolution kernels resulting from variable-step approximations for convolution *** precisely, for a class of discrete convolution kernels relevant to variable-step L1-type time discretizations, weshow that the associated quadratic form is positive definite under some easy-to-check algebraic conditions. Ourproof is based on an elementary constructing strategy by using the properties of discrete orthogonal convolutionkernels and discrete complementary convolution kernels. To our knowledge, this is the first general result onsimple algebraic conditions for the positive definiteness of variable-step discrete convolution kernels. Using theunified theory, we obtain the stability for some simple nonuniform time-stepping schemes straightforwardly.
In the present paper,we study the restricted inexact Newton-type method for solving the generalized equation 0∈f(x)+F(x),where X and Y are Banach spaces,f:X→Y is a Frechet differentiable function and F:X■Y is a set...
详细信息
In the present paper,we study the restricted inexact Newton-type method for solving the generalized equation 0∈f(x)+F(x),where X and Y are Banach spaces,f:X→Y is a Frechet differentiable function and F:X■Y is a set-valued mapping with closed *** establish the convergence criteria of the restricted inexact Newton-type method,which guarantees the existence of any sequence generated by this method and show this generated sequence is convergent linearly and quadratically according to the particular assumptions on the Frechet derivative of ***,we obtain semilocal and local convergence results of restricted inexact Newton-type method for solving the above generalized equation when the Frechet derivative of f is continuous and Lipschitz continuous as well as f+F is metrically *** application of this method to variational inequality is *** addition,a numerical experiment is given which illustrates the theoretical result.
In this paper,we develop a new sixth-order WENO scheme by adopting a convex combina-tion of a sixth-order global reconstruction and four low-order local *** the classical WENO schemes,the associated linear weights of ...
详细信息
In this paper,we develop a new sixth-order WENO scheme by adopting a convex combina-tion of a sixth-order global reconstruction and four low-order local *** the classical WENO schemes,the associated linear weights of the new scheme can be any positive numbers with the only requirement that their sum equals ***,a very simple smoothness indicator for the global stencil is *** new scheme can achieve sixth-order accuracy in smooth *** tests in some one-and two-dimensional bench-mark problems show that the new scheme has a little bit higher resolution compared with the recently developed sixth-order WENO-Z6 scheme,and it is more efficient than the classical fifth-order WENO-JS5 scheme and the recently developed sixth-order WENO6-S scheme.
This paper concerns the reconstruction of a scalar coefficient of a second-order elliptic equation in divergence form posed on a bounded domain from internal *** problem finds applications in multi-wave imaging,greedy...
详细信息
This paper concerns the reconstruction of a scalar coefficient of a second-order elliptic equation in divergence form posed on a bounded domain from internal *** problem finds applications in multi-wave imaging,greedy methods to approximate parameter-dependent elliptic problems,and image treatment with partial differential *** first show that the inverse problem for smooth coefficients can be rewritten as a linear transport *** that the coefficient is known near the boundary,we study the well-posedness of associated transport equation as well as its numerical resolution using discontinuous Galerkin *** propose a regularized transport equation that allow us to derive rigorous convergence rates of the numerical method in terms of the order of the polynomial approximation as well as the regularization *** finally provide numerical examples for the inversion assuming a lower regularity of the coefficient,and using synthetic data.
作者:
Li, SiruiXu, JieSchool of Mathematics and Statistics
Guizhou University Guiyang 550025 China LSEC
NCMIS Institute of Computational Mathematics and Scientific/Engineering Computing (ICMSEC) Academy of Mathematics and Systems Science (AMSS) Chinese Academy of Sciences Beijing China
We consider a two-tensor hydrodynamics derived from the molecular model, where high-order tensors are determined by closure approximation through the maximum entropy state or the quasi-entropy. We prove the existence ...
详细信息
This is one of our series works on discrete energy analysis of the variable-step BDF *** this part,we present stability and convergence analysis of the third-order BDF(BDF3)schemes with variable steps for linear diffu...
详细信息
This is one of our series works on discrete energy analysis of the variable-step BDF *** this part,we present stability and convergence analysis of the third-order BDF(BDF3)schemes with variable steps for linear diffusion equations,see,e.g.,[SIAM ***.,58:2294-2314]and[***.,90:1207-1226]for our previous works on the BDF2 *** this aim,we first build up a discrete gradient structure of the variable-step BDF3 formula under the condition that the adjacent step ratios are less than 1.4877,by which we can establish a discrete energy dissipation ***-robust stability and convergence analysis in the L^(2) norm are then *** the mesh robustness means that the solution errors are well controlled by the maximum time-step size but independent of the adjacent time-step *** also present numerical tests to support our theoretical results.
The paper investigates the robustness and parallel scaling properties of a novel physical factorization preconditioner with algebraic multigrid subsolves in the iterative solution of a cell-centered finite volume disc...
详细信息
The paper investigates the robustness and parallel scaling properties of a novel physical factorization preconditioner with algebraic multigrid subsolves in the iterative solution of a cell-centered finite volume discretization of the threedimensional multi-group radiation diffusion *** key idea is to take advantage of a particular kind of block factorization of the resulting system matrix and approximate the left-hand block matrix selectively spurred by parallel processing *** spectral property of the preconditioned matrix is then *** practical strategy is considered sequentially and in ***,numerical results illustrate the numerical robustness,computational efficiency and parallel strong and weak scalabilities over the real-world structured and unstructured coupled problems,showing its competitiveness with many existing block preconditioners.
作者:
Xie, PengchengState Key Laboratory of Scientific and Engineering Computing
Institute of Computational Mathematics and Scientific/Engineering Computing Academy of Mathematics and Systems Science Chinese Academy of Sciences University of Chinese Academy of Sciences ZhongGuanCun East Road No. 55 Beijing China
Optimization methods play a crucial role in various fields and applications. In some optimization problems, the derivative information of the objective function is unavailable. Such black-box optimization problems nee...
详细信息
暂无评论