In this paper, two multiscale time integrators (MTIs), motivated from two types of multiscale decomposition by either frequency or frequency and amplitude, are proposed and analyzed for solving highly oscillatory second order differential equations with a dimensionless parameter $0 < \varepsilon≤ 1.$ In fact, the solution to this equation propagates waves with wavelength at $O(\varepsilon^2)$ when $0<\varepsilon≪1,$ which brings significantly numerical burdens in practical computation. We rigorously establish two independent error bounds for the two MTIs at $O(\tau^2/\varepsilon^2)$ and $O(\varepsilon^2)$ for $\varepsilon ∈ (0,1]$ with $\tau > 0$ as step size, which imply that the two MTIs converge uniformly with linear convergence rate at $O(\tau)$ for $ε ∈ (0,1]$ and optimally with quadratic convergence rate at $O(\tau^2)$ in the regimes when either $ε=O(1)$ or $0<ε≤\tau.$ Thus the meshing strategy requirement (or $ε$-scalability) of the two MTIs is $\tau =O(1)$ for $0<ε≪1,$ which is significantly improved from $\tau =O(ε^3)$ and $\tau =O(ε^2)$ requested by finite difference methods and exponential wave integrators to the equation, respectively. Extensive numerical tests and comparisons with those classical numerical integrators are reported, which gear towards better understanding on the convergence and resolution properties of the two MTIs. In addition, numerical results support the two error bounds very well.

The linear finite element approximation of a general linear diffusion problem with arbitrary anisotropic meshes is considered. The conditioning of the resultant stiffness matrix and the Jacobi preconditioned stiffness matrix is investigated using a density function approach proposed by Fried in 1973. It is shown that the approach can be made mathematically rigorous for general domains and used to develop bounds on the smallest eigenvalue and the condition number that are sharper than existing estimates in one and two dimensions and comparable in three and higher dimensions. The new results reveal that the mesh concentration near the boundary has less influence on the condition number than the mesh concentration in the interior of the domain. This is especially true for the Jacobi preconditioned system where the former has little or almost no influence on the condition number. Numerical examples are presented.

In this paper, we propose a finite difference/collocation method for two-dimensional time fractional diffusion equation with generalized fractional operator. The main purpose of this paper is to design a high order numerical scheme for the new generalized time fractional diffusion equation. First, a finite difference approximation formula is derived for the generalized time fractional derivative, which is verified with order $2-\alpha$ $(0<\alpha<1)$. Then, collocation method is introduced for the two-dimensional space approximation. Unconditional stability of the scheme is proved. To make the method more efficient, the alternating direction implicit method is introduced to reduce the computational cost. At last, numerical experiments are carried out to verify the effectiveness of the scheme.

A simple and efficient spectral method for solving the second, third order and fourth order elliptic equations with variable coefficients and nonlinear differential equations is presented. It is different from spectral-collocation method which leads to dense, ill-conditioned matrices. The spectral method in this paper solves for the coefficients of the solution in a Chebyshev series, leads to discrete systems with special structured matrices which can be factorized and solved efficiently. We also extend the method to boundary value problems in two space dimensions and solve 2-D separable equation with variable coefficients. As an application, we solve Cahn-Hilliard equation iteratively via first-order implicit time discretization scheme. Ample numerical results indicate that the proposed method is extremely accurate and efficient.

In this paper, we are concerned with the properties of positive solutions of the following nonlinear integral systems on the Heisenberg group $\mathbb{H}^n$, \begin{equation} \left\{\begin{array}{ll} u(x)=\int_{\mathbb{H}^n}\frac{v^{q}(y)w^{r}(y)}{|x^{-1}y|^\alpha|y|^\beta}\,dy,\\ v(x)=\int_{\mathbb{H}^n}\frac{u^{p}(y)w^{r}(y)}{|x^{-1}y|^\alpha|y|^\beta}\,dy,\\ w(x)=\int_{\mathbb{H}^n}\frac{u^{p}(y)v^{q}(y)}{|x^{-1}y|^\alpha|y|^\beta}\,dy,\\ \end{array}\right.\end{equation}
for $x\in \mathbb{H}^n$, where $0<\alpha<Q=2n+2$, $n\geq3$, $\beta\geq0$, $\alpha+\beta<Q$, and $p,q,r > 1$ satisfying $\frac{1}{p+1} $+ $\frac{1}{q+1} + \frac{1}{r+1} = \frac{Q+α+β}{Q}.$ We show that positive solution triples $(u,v,w)\in L^{p+1}(\mathbb{H}^n)\times L^{q+1}(\mathbb{H}^n)\times L^{r+1}(\mathbb{H}^n)$ are bounded and they converge to zero when $|x|→∞.$