Discretization of the solutions to Poisson’s equation

Sharp estimates (in the power scale) are obtained for the discretization error in the solutions to Poisson’s equation whose right-hand side belongs to a Korobov class. Compared to the well-known Korobov estimate, the order is almost doubled and has an ultimate value in the power scale.     

Finite-Dimensional Approximations for Scattering Theory Operators in the Wave-Packet Representation

We consider some types of packet discretization for continuous spectra in quantum scattering problems. As we previously showed, this discretization leads to a convenient finite-dimensional (i.e., matrix) approximation for integral operators in the scattering theory and allows reducing the solution of singular integral equations connected with the scattering theory to some suitable purely algebraic equations on an analytic basis. All singularities are explicitly singled out. Our primary emphasis is on realizing the method practically.     

Implementation of some higher-order convection schemes on non-uniform grids

A generalized formulation is applied to implement the quadratic upstream interpolation (QUICK) scheme, the second-order upwind (SOU) scheme and the second-order hybrid scheme (SHYBRID) on non-uniform grids. The implementation method is simple. The accuracy and efficiency of these higher-order schemes on non-uniform grids are assessed. Three well-known bench mark convection-diffusion problems and a fluid flow problem are revisited using non-uniform grids. These are: (1) transport of a scalar tracer by a uniform velocity field; (2) heat transport in a recirculating flow; (3) two-dimensional non-linear Burgers equations; and (4) a two-dimensional incompressible Navier-Stokes flow which is similar to the classical lid-driven cavity flow. The known exact solutions of the last three problems make it possible to thoroughly evaluate accuracies of various uniform and non-uniform grids. Higher accuracy is obtained for fewer grid points on non-uniform grids. The order of accuracy of the examined schemes is maintained for some tested problems if the distribution of non-uniform grid points is properly chosen.     

带跳的非线性随机微分方程的Lyapunov指数的估计

本文研究了带跳的非线性随机微分方程Lyapunov指数的估计,在适当的条件下,确定其Lyapunov指数q的值．对于给定的步长h,考虑此微分系统的Euler离散化模型,给出了的理论误差估计.     

Polynomial expansion of basis sets: Alternatives to fully optimized exponents

Alternatives based on polynomial expansions of gaussian basis set exponents are introduced and evaluated. The formulas presented here outperform methods based upon the even-tempered formula or combinations of it. They closely match the performance of other methods based upon larger polynomial expansions of the logarithm of the exponents using the same or one less parameter per orbital angular symmetry.     

Optimal monotone hysteretic Markov policies in anM/M/1 queueing model with switching costs and finite time horizon

In anM/M/1 queueing model, a decision maker can choosem pairs of arrival- and service rates. He can change his action at any time epoch, a switch of action costs an amount depending on the size of the switch. Besides that there are continuously incurring costs. Over a finite time horizon, there exists an optimal monotone hysteretic Markov policy. This is shown essentially by the technique of time discretization.The work producing this article was done during a half year stay at the University of Leiden, The Netherlands, with Prof. Arie Hordijk. A technical report (a more detailled version of this article) was written there [6]. The opportunity for this stay was given by the University of Bonn, Germany, where the author, at that time, worked as scientific assistant of Prof. M. Schäl.     

Orbital-Free Density Functional Theory for Molecular Structure Calculations

We give here an overview of the orbital-free density functional theory that is used for modeling atoms and molecules.We review typical approximations to the kinetic energy,exchange-correlation corrections to the kinetic and Hartree energies, and constructions of the pseudopotentials.We discuss numerical discretizations for the orbital-free methods and include several numerical results for illustrations.     

Wave-packet continuum discretization for quantum scattering

A general approach to a solution of few- and many-body scattering problems based on a continuum-discretization procedure is described in detail. The complete discretization of continuous spectrum is realized using stationary wave packets which are the normalized states constructed from exact non-normalized continuum states. Projecting the wave functions and all scattering operators like t$t$-matrix, resolvent, etc. on such a wave-packet basis results in a formulation of quantum scattering problem entirely in terms of discrete elements and linear equations with regular matrices. It is demonstrated that there is a close relation between the above stationary wave packets and pseudostates which are employed often to approximate the scattering states with a finite L₂$L_2$ basis. Such a fully discrete treatment of complicated few- and many-body scattering problems leads to significant simplification of their practical solution. Also we get finite-dimensional approximations for complicated operators like effective interactions between composite particles constructed via the Feshbach-type projection formalism. As illustrations to this general approach we consider several important particular problems including multichannel scattering and scattering in the three-nucleon system within the Faddeev framework.     

Multilevel preconditioning and low‐rank tensor iteration for space–time simultaneous discretizations of parabolic PDEs

This paper addresses the solution of parabolic evolution equations simultaneously in space and time as may be of interest in, for example, optimal control problems constrained by such equations. As a model problem, we consider the heat equation posed on the unit cube in Euclidean space of moderately high dimension. An a priori stable minimal residual Petrov–Galerkin variational formulation of the heat equation in space–time results in a generalized least squares problem. This formulation admits a unique, quasi‐optimal solution in the natural space–time Hilbert space and serves as a basis for the development of space–time compressive solution algorithms. The solution of the heat equation is obtained by applying the conjugate gradient method to the normal equations of the generalized least squares problem. Starting from stable subspace splittings in space and in time, multilevel space–time preconditioners for the normal equations are derived. In order to reduce the complexity of the full space–time problem, all computations are performed in a compressed or sparse format called the hierarchical Tucker format, supposing that the input data are available in this format. In order to maintain sparsity, compression of the iterates within the hierarchical Tucker format is performed in each conjugate gradient iteration. Its application to vectors in the hierarchical Tucker format is detailed. Finally, numerical results in up to five spatial dimensions based on the recently developed htucker toolbox for MATLAB are presented. Copyright © 2014 John Wiley & Sons, Ltd.