On a stable calculation of the normal to a surface given approximately

The paper proposes a stable method for constructing a normal to a surface given approximately. The normal is calculated as the gradient of the function in the surface equation. As is known, the problem of calculating the derivative is ill-posed. In the paper, an approach is adopted to solving this problem as to the problem of calculating the values of an unbounded operator. To construct its stable solution, the principle of minimum of the smoothing functional in Morozov’s formulation is used. The normal is obtained in the form of a Fourier series in the expansion in terms of eigenfunctions of the Laplace operator in a rectangle with boundary conditions of the second kind. The functional stabilizer uses the Laplacian, which makes it possible to obtain a normal in the form of a Fourier series that converges uniformly to the exact normal vector as the error in the surface definition tends to zero. The resulting approximate normal vector can be used to solve various problems of mathematical physics using surface integrals, normal derivatives, simple and double layer potentials.

В работе предлагается устойчивый метод построения нормали к поверхности, заданной приближённо. Нормаль вычисляется как градиент функции в уравнении поверхности. Как известно, задача вычисления производной является некорректно поставленной. В работе принят подход к решению этой задачи как к задаче вычисления значений неограниченного оператора. Для построения её устойчивого решения используется принцип минимума сглаживающего функционала в формулировке Морозова. Нормаль получена в виде ряда Фурье в разложении по собственным функциям оператора Лапласа в прямоугольнике с краевыми условиями второго рода. В стабилизаторе функционала используется лапласиан, что позволяет получить нормаль в виде ряда Фурье, равномерно сходящегося к точному вектору нормали при стремлении к нулю погрешности в задании поверхности. Полученный приближенный вектор нормали может использоваться при решении различных задач математической физики, использующих поверхностные интегралы, нормальные производные, потенциалы простого и двойного слоя.