We propose an inexact affine scaling Levenberg-Marquardt method for solving bound-constrained semismooth equations under the local error bound assumption which is much weaker than the standard nonsingularity condition...We propose an inexact affine scaling Levenberg-Marquardt method for solving bound-constrained semismooth equations under the local error bound assumption which is much weaker than the standard nonsingularity condition. The affine scaling Levenberg-Marquardt equation is based on a minimization of the squared Euclidean norm of linearized model adding a quadratic affine scaling matrix to find a solution which belongs to the bounded constraints on variable. The global convergence and the superlinear convergence rate are proved.Numerical results show that the new algorithm is efficient.展开更多
An h-adaptive meshless method is proposed in this paper. The error estimation is based on local fit technology, usually confined to Voronoi Cells. The error is achieved by comparison of the computational results with ...An h-adaptive meshless method is proposed in this paper. The error estimation is based on local fit technology, usually confined to Voronoi Cells. The error is achieved by comparison of the computational results with smoothed ones, which are projected with Taylor series. Voronoi Cells are introduced not only for integration of potential energy but also for guidance of refinement. New nodes are placed within those cells with high estimated error. At the end of the paper, two numerical examples with severe stress gradient are analyzed. Through adaptive analysis accurate results are obtained at critical subdomains, which validates the efficiency of the method.展开更多
For a class of non-smooth Fredholm integral equations we prove thatRichardson extrapolation method can be used to increase accuracy for finite ele-ment methods.On the basis of local refinement mesh,a superconvergence ...For a class of non-smooth Fredholm integral equations we prove thatRichardson extrapolation method can be used to increase accuracy for finite ele-ment methods.On the basis of local refinement mesh,a superconvergence result isshown.Our theory also covers a class of non-smooth Wiener-Hopf equations andan application includes the calculation in certain linear elastic fracture problems.展开更多
Allen and Liu (1995) introduced a new method for a time-dependent convection dominated diffusion problem, which combines the modified method of characteristics and method of streamline diffusion. But they ignored the ...Allen and Liu (1995) introduced a new method for a time-dependent convection dominated diffusion problem, which combines the modified method of characteristics and method of streamline diffusion. But they ignored the fact that the accuracy of time discretization decays at half an order when the characteristic line goes out of the domain. In present paper, the author shows that, as a remedy, a simple lumped scheme yields a full accuracy approximation. Forthermore, some local error bounds independent of the small viscosity axe derived for this scheme outside the boundary layers.展开更多
In this paper, a projected gradient trust region algorithm for solving nonlinear equality systems with convex constraints is considered. The global convergence results are developed in a very general setting of comput...In this paper, a projected gradient trust region algorithm for solving nonlinear equality systems with convex constraints is considered. The global convergence results are developed in a very general setting of computing trial directions by this method combining with the line search technique. Close to the solution set this method is locally Q-superlinearly convergent under an error bound assumption which is much weaker than the standard nonsingularity condition.展开更多
In this paper,we apply local discontinuous Galerkin methods to the pattern formation dynamical model in polymerizing action flocks.Optimal error estimates for the density and filament polarization in different norms a...In this paper,we apply local discontinuous Galerkin methods to the pattern formation dynamical model in polymerizing action flocks.Optimal error estimates for the density and filament polarization in different norms are established.We use a semi-implicit spectral deferred correction time method for time discretization,which allows a relative large time step and avoids computation of a Jacobian matrix.Numerical experiments are presented to verify the theoretical analysis and to show the capability for simulations of action wave formation.展开更多
针对电力系统低频振荡非线性时变的特点,提出了一种基于改进局部均值分解(local mean decomposition,LMD)的电力系统低频振荡信号分析方法。利用改进的局部均值分解,电力系统中的单一多模态测量信号可以分解为一组乘积函数(product func...针对电力系统低频振荡非线性时变的特点,提出了一种基于改进局部均值分解(local mean decomposition,LMD)的电力系统低频振荡信号分析方法。利用改进的局部均值分解,电力系统中的单一多模态测量信号可以分解为一组乘积函数(product function,PF)分量的和。每个PF分量可以表示为一个调幅(amplitude modulated,AM)信号和一个调频(frequency modulated,FM)信号的乘积。其中,AM信号可以近似当作相应振荡模态的瞬时幅值,并由此计算阻尼信息;FM信号可以通过直接正交和插值相结合的综合方法,计算PF的瞬时频率。数值仿真和实际测量信号的计算结果证明了所提方法的有效性和可行性。展开更多
In this paper we continue our effort in Liu-Shu (2004) and Liu-Shu (2007) for developing local discontinuous Galerkin (LDG) finite element methods to discretize moment models in semiconductor device simulations. We co...In this paper we continue our effort in Liu-Shu (2004) and Liu-Shu (2007) for developing local discontinuous Galerkin (LDG) finite element methods to discretize moment models in semiconductor device simulations. We consider drift-diffusion (DD) and high-field (HF) models of one-dimensional devices, which involve not only first derivative convection terms but also second derivative diffusion terms, as well as a coupled Poisson potential equation. Error estimates are obtained for both models with smooth solutions. The main technical difficulties in the analysis include the treatment of the inter-element jump terms which arise from the discontinuous nature of the numerical method, the nonlinearity, and the coupling of the models. A simulation is also performed to validate the analysis.展开更多
In this paper, we study the superconvergence of the error for the local discontinuous Galerkin （LDG） finite element method for one-dimensional linear parabolic equations when the alternating flux is used. We prove t...In this paper, we study the superconvergence of the error for the local discontinuous Galerkin （LDG） finite element method for one-dimensional linear parabolic equations when the alternating flux is used. We prove that if we apply piecewise k-th degree polynomials, the error between the LDG solution and the exact solution is （k ＋ 2）-th order superconvergent at the Radau points with suitable initial discretization. Moreover, we also prove the LDG solution is （k ＋ 2）-th order superconvergent for the error to a particular projection of the exact solution. Even though we only consider periodic boundary condition, this boundary condition is not essential, since we do not use Fourier analysis. Our analysis is valid for arbitrary regular meshes and for P^k polynomials with arbitrary k ≥ 1. We perform numerical experiments to demonstrate that the superconvergence rates proved in this paper are sharp.展开更多
基金Supported by National Natural Science Foundation of China(No.11571074)Scientific Research Fund of Hunan Provincial Education Department(No.18A351,17C0393)Natural Science Foundation of Hunan Province(No.2019JJ50105)
文摘We propose an inexact affine scaling Levenberg-Marquardt method for solving bound-constrained semismooth equations under the local error bound assumption which is much weaker than the standard nonsingularity condition. The affine scaling Levenberg-Marquardt equation is based on a minimization of the squared Euclidean norm of linearized model adding a quadratic affine scaling matrix to find a solution which belongs to the bounded constraints on variable. The global convergence and the superlinear convergence rate are proved.Numerical results show that the new algorithm is efficient.
基金Project supported by the National Natural Science Foundation of China (No. 50175060).
文摘An h-adaptive meshless method is proposed in this paper. The error estimation is based on local fit technology, usually confined to Voronoi Cells. The error is achieved by comparison of the computational results with smoothed ones, which are projected with Taylor series. Voronoi Cells are introduced not only for integration of potential energy but also for guidance of refinement. New nodes are placed within those cells with high estimated error. At the end of the paper, two numerical examples with severe stress gradient are analyzed. Through adaptive analysis accurate results are obtained at critical subdomains, which validates the efficiency of the method.
文摘For a class of non-smooth Fredholm integral equations we prove thatRichardson extrapolation method can be used to increase accuracy for finite ele-ment methods.On the basis of local refinement mesh,a superconvergence result isshown.Our theory also covers a class of non-smooth Wiener-Hopf equations andan application includes the calculation in certain linear elastic fracture problems.
文摘Allen and Liu (1995) introduced a new method for a time-dependent convection dominated diffusion problem, which combines the modified method of characteristics and method of streamline diffusion. But they ignored the fact that the accuracy of time discretization decays at half an order when the characteristic line goes out of the domain. In present paper, the author shows that, as a remedy, a simple lumped scheme yields a full accuracy approximation. Forthermore, some local error bounds independent of the small viscosity axe derived for this scheme outside the boundary layers.
基金Supported by the National Natural Science Foundation of China (10871130)the Research Fund for the Doctoral Program of Higher Education of China (20093127110005)the Scientific Computing Key Laboratory of Shanghai Universities
文摘In this paper, a projected gradient trust region algorithm for solving nonlinear equality systems with convex constraints is considered. The global convergence results are developed in a very general setting of computing trial directions by this method combining with the line search technique. Close to the solution set this method is locally Q-superlinearly convergent under an error bound assumption which is much weaker than the standard nonsingularity condition.
基金supported by National Natural Science Foundation of China(Grant Nos.11801569 and 11571367)Natural Science Foundation of Shandong Province(CN)(Grant Nos.ZR2018BA011 and ZR2019MA015)+1 种基金the Fundamental Research Funds for the Central Universities(Grant Nos.18CX02021A and 18CX05003A)National Science Foundation of USA(Grant No.DMS-1818467).
文摘In this paper,we apply local discontinuous Galerkin methods to the pattern formation dynamical model in polymerizing action flocks.Optimal error estimates for the density and filament polarization in different norms are established.We use a semi-implicit spectral deferred correction time method for time discretization,which allows a relative large time step and avoids computation of a Jacobian matrix.Numerical experiments are presented to verify the theoretical analysis and to show the capability for simulations of action wave formation.
基金supported in part by the Program for Natural Science Foundation of Shandong Province, China (Grant No. ZR2009AM015)National Natural Science Foundation of China (Grant Nos. 10501031, 10911140142)+2 种基金supported in part by Army Research Office of USA (Grant No. W911NF-08-1-0520)National Science Foundation of USA (Grant No.DMS-0809086)supported in part by the Abdus Salam International Center for Theoretical Physics during the first author’s visit to the Mathematics Section
文摘In this paper we continue our effort in Liu-Shu (2004) and Liu-Shu (2007) for developing local discontinuous Galerkin (LDG) finite element methods to discretize moment models in semiconductor device simulations. We consider drift-diffusion (DD) and high-field (HF) models of one-dimensional devices, which involve not only first derivative convection terms but also second derivative diffusion terms, as well as a coupled Poisson potential equation. Error estimates are obtained for both models with smooth solutions. The main technical difficulties in the analysis include the treatment of the inter-element jump terms which arise from the discontinuous nature of the numerical method, the nonlinearity, and the coupling of the models. A simulation is also performed to validate the analysis.
文摘In this paper, we study the superconvergence of the error for the local discontinuous Galerkin （LDG） finite element method for one-dimensional linear parabolic equations when the alternating flux is used. We prove that if we apply piecewise k-th degree polynomials, the error between the LDG solution and the exact solution is （k ＋ 2）-th order superconvergent at the Radau points with suitable initial discretization. Moreover, we also prove the LDG solution is （k ＋ 2）-th order superconvergent for the error to a particular projection of the exact solution. Even though we only consider periodic boundary condition, this boundary condition is not essential, since we do not use Fourier analysis. Our analysis is valid for arbitrary regular meshes and for P^k polynomials with arbitrary k ≥ 1. We perform numerical experiments to demonstrate that the superconvergence rates proved in this paper are sharp.
