О предельных решениях разностных уравнений, содержащих оператор Лапласа

Бурное развитие вычислительной техники облегчило применение математического моделирования во многих областях человеческой деятельности. В настоящее время математическим моделированием занимается огромное количество специалистов, использующих уже известные и хорошо обоснованные численные методы. При обосновании ряда методов были доказаны теоремы эквивалентности [1-4] о связи сходимости численного решения к точному с аппроксимацией и устойчивостью, а также теорема [5] об условиях монотонности численного решения. Следует, однако, помнить, что все эти теоремы были доказаны для линейных или линеаризованных уравнений на равномерных сетках. В случае же нелинейных уравнений и при применении неоднородных (адаптивных) сеток погрешность аппроксимации может оказаться несходящейся и предельное решение при сколь угодно большом увеличении числа точек сетки N может отличаться от точного решения. На возникновение несходящейся аппроксимации уравнения теплопроводности разностным уравнением на неравномерной сетке было обращено внимание в [1, 2]. В механике жидкости и газа нарушение равномерности сетки приводит к образованию «энтропийных» следов [6, 7]. В работе рассматривается проблема различия точного решения и предельного при N → ∞ решения нелинейного уравнения теплопроводности и уравнения электростатики в случае, когда соседние ячейки сетки сильно различаются.