Enforcing the non-negativity constraint and maximum principles for diffusion with decay on general computational grids

In this paper, we consider anisotropic diffusion with decay, and the diffusivity coefficient to be a second-order symmetric and positive definite tensor. It is well-known that this particular equation is a second-order elliptic equation, and satisfies a maximum principle under certain regularity assumptions. However, the finite element implementation of the classical Galerkin formulation for both anisotropic and isotropic diffusion with decay does not respect the maximum principle. We first show that the numerical accuracy of the classical Galerkin formulation deteriorates dramatically with increase in the decay coefficient for isotropic medium and violates the discrete maximum principle. However, in the case of isotropic medium, the extent of violation decreases with mesh refinement. We then show that, in the case of anisotropic medium, the classical Galerkin formulation for anisotropic diffusion with decay violates the discrete maximum principle even at lower values of decay coefficient and does not vanish with mesh refinement. We then present a methodology for enforcing maximum principles under the classical Galerkin formulation for anisotropic diffusion with decay on general computational grids using optimization techniques. Representative numerical results (which take into account anisotropy and heterogeneity) are presented to illustrate the performance of the proposed formulation

Discrete Maximum Principle for Nonsmooth Optimal Control Problems with Delays

We consider optimal control problems for discrete-time systems with delays. The main goal is to derive necessary optimality conditions of the discrete maximum principle type in the case of nonsmooth minimizing functions. We obtain two independent forms of the discrete maximum principle with transversality conditions described in terms of subdifferentials and superdifferentials, respectively. The superdifferential form is new even for non-delayed systems and may be essentially stronger than a more conventional subdifferential form in some situations

Optimal production scheduling and inventory control by the discrete maximum principle

Call number: LD2668 .R4 1967 K4

On nonlinear artificial viscosity, discrete maximum principle and hyperbolic conservation laws

A finite element method for Burgers' equation is studied. The method is analyzed using techniques from stabilized finite element methods and convergence to entropy solutions is proven under certain hypotheses on the artificial viscosity. In particular we assume that a discrete maximum principle holds. We then construct a nonlinear artificial viscosity that satisfies the assumptions required for convergence and that can be tuned to minimize artificial viscosity away from local extrema. The theoretical results are exemplified on a numerical exampl