We propose a simple fast spectral method for the Boltzmann collision operator with general collision kernels. In contrast to the direct spectral method [L. Pareschi and G. Russo, SIAM J. Numer. Anal., 37 (2000), pp. 1217–1245; I. M. Gamba and S. H. Tharkabhushanam, J. Comput. Phys., 228 (2009), pp. 2012–2036], which requires O(N6) memory to store precomputed weights and has O(N6) numerical complexity, the new method has complexity O(MN4 log N), where N is the number of discretization points in each of the three velocity dimensions and M is the total number of discretization points on the sphere and M ≪ N2. Furthermore, it requires no precomputation for the variable hard sphere model and only O(MN4) memory to store precomputed functions for more general collision kernels. Although a faster spectral method is available [C. Mouhot and L. Pareschi, Math. Comp., 75 (2006), pp. 1833–1852] (with complexity O(MN3 logN)), it works only for hard sphere molecules, thus limiting its use for practical problems. Our new method, on the other hand, can apply to arbitrary collision kernels. A series of numerical tests is performed to illustrate the efficiency and accuracy of the proposed method.

We consider the Schrodinger equation with a logarithmic nonlinearity and a repulsive harmonic potential. Depending on the parameters of the equation, the solution may or may not be dispersive. When dispersion occurs, it does with an exponential rate in time. To control this, we change the unknown function through a generalized lens transform. This approach neutralizes the possible boundary effects, and could be used in the case of the Schrodinger equation without potential. We then employ standard splitting methods on the new equation via a nonuniform grid, after the logarithmic nonlinearity has been regularized. We also discuss the case of a power nonlinearity and give some results concerning the error estimates of the first-order Lie-Trotter splitting method for both cases of nonlinearities. Finally extensive numerical experiments are reported to investigate the dynamics of the equations.

In this paper, we develop and analyze a new multiscale discontinuous Galerkin (DG) method for one-dimensional stationary Schr\"{o}dinger equations with open boundary conditions \bo{which have highly oscillating solutions}. Our method uses a smaller finite element space than the WKB local DG (WKB-LDG) method while achieving the same order of accuracy with no resonance errors. We prove that the DG approximation converges optimally with respect to the mesh size $h$ in $L^2$ norm without the typical constraint \bo{that $h$ has to be smaller than the wave length}. Numerical experiments were carried out to verify the second order optimal convergence rate of the method and to demonstrate its ability to capture oscillating solutions on coarse meshes in the applications to Schr\"{o}dinger equations.

No abstract uploaded!

For the steady-state solution of an integraldifferential equation from a two-dimensional model in transport theory, we shall derive and study a nonsymmetric algebraic Riccati equation B--XF--F+ X+ XB+ X= 0, where F I-s PD, B-(b I+ s P) D-and B+ b I+ s PD+ with a nonnegative matrix P, positive diagonal matrices D, and nonnegative parameters f, b b/(1-f) and s s/(1-f). We prove the existence of the minimal nonnegative solution X under the physically reasonable assumption f+ b+ s P (D++ D-)< 1, and study its numerical computation by fixed-point iteration, Newtons method and doubling. We shall also study several special cases; eg when b = 0 and P is low-ranked, then X= s 2 UV is low-ranked and can be computed using more efficient iterative processes in U and V. Numerical examples will be given to illustrate our theoretical results.