We present a novel energy-based numerical analysis of semilinear diffusion-reaction boundary value problems. Based on a suitable variational setting, the proposed computational scheme can be seen as an energy minimisation approach. More specifically, this procedure aims to generate a sequence of numerical approximations, which results from the iterative solution of related (stabilised) linearised discrete problems, and tends to a local minimum of the underlying energy functional. Simultaneously, the finite-dimensional approximation spaces are adaptively refined; this is implemented in terms of a new mesh refinement strategy in the context of finite element discretisations, which again relies on the energy structure of the problem under consideration, and does not involve any a posteriori error indicators. In combination, the resulting adaptive algorithm consists of an iterative linearisation procedure on a sequence of hierarchically refined discrete spaces, which we prove to converge towards a solution of the continuous problem in an appropriate sense. Numerical experiments demonstrate the robustness and reliability of our approach for a series of examples.
翻译:我们对半线性扩散-反应边界值问题进行了新型的基于能源的数字分析。根据适当的变异设置,拟议的计算方法可被视为一种能源最小化方法。更具体地说,这一程序旨在产生一系列数字近似,这是相关(稳定)线性离散问题的迭接解决方案的结果,而且趋向于当地最低基本能量功能。与此同时,有限维度近距离空间正在适应性地改进;这是在有限元素分解背景下采用的新网状精细化战略实施的,这一战略再次依赖所审议问题的能源结构,并不涉及任何事后误差指标。综合起来,由此产生的适应算法包括按分级精细化离散空间序列的迭接线化程序,我们证明这种程序在适当意义上可以汇合到持续问题的解决方案。数字实验表明,我们的方法对于一系列例子来说是稳健和可靠的。