A High-Order Method for Stiff Boundary Value Problems with Turning Points

This paper describes some high-order collocation-like methods for the numerical solution of stiff boundary-value problems with turning points. The presentation concentrates on the implementation of these methods in conjunction with the implementation of the a priori mesh construction algorithm introduced by Kreiss, Nichols and Brown [SIAM J. Numer. Anal., 23 (1986), pp. 325–368] for such problems. Numerical examples are given showing the high accuracy which can be obtained in solving the boundary value problem for singularly perturbed ordinary differential equations with turning points

Fast iterative solution of reaction-diffusion control problems arising from chemical processes

PDE-constrained optimization problems, and the development of preconditioned iterative methods for the efficient solution of the arising matrix system, is a field of numerical analysis that has recently been attracting much attention. In this paper, we analyze and develop preconditioners for matrix systems that arise from the optimal control of reaction-diffusion equations, which themselves result from chemical processes. Important aspects in our solvers are saddle point theory, mass matrix representation and effective Schur complement approximation, as well as the outer (Newton) iteration to take account of the nonlinearity of the underlying PDEs

Symmetry breaking and phase coexistence in a driven diffusive two-channel system

We consider classical hard-core particles moving on two parallel chains in the same direction. An interaction between the channels is included via the hopping rates. For a ring, the stationary state has a product form. For the case of coupling to two reservoirs, it is investigated analytically and numerically. In addition to the known one-channel phases, two new regions are found, in particular the one, where the total density is fixed, but the filling of the individual chains changes back and forth, with a preference for strongly different densities. The corresponding probability distribution is determined and shown to have an universal form. The phase diagram and general aspects of the problem are discussed.Comment: 12 pages, 10 figures, to appear in Phys.Rev.