Least-squares spectral method for solving advective population balance problems

The combined CFD-PBE (population balance models) are computationally intensive requiring efficient numerical methods for solving practical problems. In this paper, a high order method is presented based on the least-squares method (LSM) for the solution of a spatial-dependent population balance equation which includes advective processes. Numerical experiments are performed in order to study the behavior of the proposed method for one-dimensional cases using model problems with analytical solutions.     

On distributional solutions of local and nonlocal problems of porous medium type

We present a theory of well-posedness and a priori estimates for bounded distributional (or very weak) solutions of

(0.1)${?}_{t}u?{\mathfrak{L}}^{\sigma ,\mu }[\phi (u)]=g(x,t)\text{in}{\mathbb{R}}^N\times (0,T),$

where φ is merely continuous and nondecreasing, and ${\mathfrak{L}}^{\sigma ,\mu }$ is the generator of a general symmetric Lévy process. This means that ${\mathfrak{L}}^{\sigma ,\mu }$ can have both local and nonlocal parts like, e.g., ${\mathfrak{L}}^{\sigma ,\mu }=\mathrm{\Delta }?{(?\mathrm{\Delta })}^{\frac{1}{2}}$. New uniqueness results for bounded distributional solutions to this problem and the corresponding elliptic equation are presented and proven. A key role is played by a new Liouville type result for ${\mathfrak{L}}^{\sigma ,\mu }$. Existence and a priori estimates are deduced from a numerical approximation, and energy-type estimates are also obtained.     

The discontinuous Galerkin method for fractional degenerate convection-diffusion equations

We propose and study discontinuous Galerkin methods for strongly degenerate convection-diffusion equations perturbed by a fractional diffusion (Lévy) operator. We prove various stability estimates along with convergence results toward properly defined (entropy) solutions of linear and nonlinear equations. Finally, the qualitative behavior of solutions of such equations are illustrated through numerical experiments.     

On Error Bounds for Approximation Schemes for Non-Convex Degenerate Elliptic Equations

In this paper we provide estimates of the rates of convergence of monotone approximation schemes for non-convex equations in one space-dimension. The equations under consideration are the degenerate elliptic Isaacs equations with x-depending coefficients, and the results applies in particular to certain finite difference methods and control schemes based on the dynamic programming principle. Recently, Krylov, Barles, and Jakobsen obtained similar estimates for convex Hamilton–Jacobi–Bellman equations in arbitrary space-dimensions. Our results are only valid in one space-dimension, but they are the first results of this type for non-convex second-order equations.     

An efficient step size control for continuation methods

In solving a nonlinear equation by the use of a continuation method one of the crucial problems is the choice of the step sizes. We present a model for the total computational cost of a standard numerical continuation process and solve the problem of optimal step size control for this model. Using the theoretical results as a basis, we develop an adaptive step size algorithm for Newton's method. This procedure is computationally inexpensive and it gives quite satisfactory results compared to some other numerical experiments found in the literature.