An iterative discontinuous Galerkin(DG)method is proposed to solve the nonlinear Poisson Boltzmann(PB)*** first identify a function space inwhich the solution of the nonlinear PB equation is iteratively approximated t...
详细信息
An iterative discontinuous Galerkin(DG)method is proposed to solve the nonlinear Poisson Boltzmann(PB)*** first identify a function space inwhich the solution of the nonlinear PB equation is iteratively approximated through a series of linear PB equations,while an appropriate initial guess and a suitable iterative parameter are selected so that the solutions of linear PB equations are monotone within the identified solution *** the spatial discretization we apply the direct discontinuous Galerkin method to those linear PB *** precisely,we use one initial guess when the Debye parameter l=O(1),and a special initial guess for l≪1 to ensure *** iterative parameter is carefully chosen to guarantee the existence,uniqueness,and convergence of the *** particular,iteration steps can be reduced for a variable iterative *** one and two-dimensional numerical results are carried out to demonstrate both accuracy and capacity of the iterative DG method for both cases of l=O(1)and l≪***(m+1)th order of accuracy for L2 and mth order of accuracy for H1 for Pm elements are numerically obtained.
We investigate the well-posedness of a perfectly matched layer model developed by Cohen and Monk [Comput Methods Appl Mech Eng 169 (1999), 197-217]. A new time-domain finite element method is proposed to solve the mod...
详细信息
The initial value problems of nonlinear ordinary differential equations which contain stiff and nonstiff terms often arise from many applications. In order to reduce the computation cost, implicit–explicit (IMEX) met...
详细信息
The initial value problems of nonlinear ordinary differential equations which contain stiff and nonstiff terms often arise from many applications. In order to reduce the computation cost, implicit–explicit (IMEX) methods are often applied to these problems, i.e. the stiff and non-stiff terms are discretized by using implicit and explicit methods, respectively. In this paper, we mainly consider the nonlinear stiff initial-value problems satisfying the one-sided Lipschitz condition and a class of singularly perturbed initial-value problems, and present two classes of the IMEX multistep methods by combining implicit one-leg methods with explicit linear multistep methods and explicit one-leg methods, respectively. The order conditions and the convergence results of these methods are obtained. Some efficient methods are constructed. Some numerical examples are given to verify the validity of the obtained theoretical results.
The Gross–Pitaevskii equation is the model equation of the single-particle wave function in a Bose–Einstein condensation. A computation difficulty of the Gross–Pitaevskii equation comes from the semiclassical probl...
详细信息
The Gross–Pitaevskii equation is the model equation of the single-particle wave function in a Bose–Einstein condensation. A computation difficulty of the Gross–Pitaevskii equation comes from the semiclassical problem in supercritical case. In this paper, we apply a diffeomorphism to transform the original one-dimensional Gross–Pitaevskii equation into a modified equation. The adaptive grids are constructed through the interpolating wavelet method. Then, we use the time-splitting finite difference method with the wavelet-adaptive grids to solve the modified Gross–Pitaevskii equation, where the approximation to the second-order derivative is given by the Lagrange interpolation method. At last, the numerical results are given. It is shown that the obtained time-splitting finite difference method with the wavelet-adaptive grids is very efficient for solving the one-dimensional semiclassical Gross–Pitaevskii equation in supercritical case and it is suitable to deal with the local high oscillation of the solution.
Particulate flows are commonly encountered in both engineering and environmental applications. The discrete element method (DEM) has attracted plentiful attentions since it can predict the whole motion of the particul...
详细信息
This paper identifies some scaling relationships between solar activity and geomagnetic activity. We examine the scaling properties of hourly data for two geomagnetic indices (ap and AE), two solar indices (solar X-ra...
详细信息
A kind of equal-task multiple traveling salesman problem (ET-mTSP) was proposed based on the mTSP and its corresponding mathematical model was constructed;Then, a series of discrete operations for firefly algorithm (F...
详细信息
We propose general spectral and pseudo-spectral Jacobi–Galerkin methods for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide rig...
We propose general spectral and pseudo-spectral Jacobi–Galerkin methods for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide rigorous error analysis for spectral and pseudo-spectral Jacobi–Galerkin methods, which show that the errors of the approximate solution decay exponentially in $$L^\infty $$ norm and weighted $$L^2$$ -norm. The numerical examples are given to illustrate the theoretical results.
In this paper, applying majorization inequalities, new upper and lower bounds for summation of eigenvalues (including the trace) of the solution of the continuous algebraic Riccati equation are presented. Correspondin...
详细信息
暂无评论