Volume 8, Issue 1

A numerical upscaling approach, NU, for solving multiscale elliptic problems is discussed. The main components of this NU are: i) local solve of aux-iliary problems in grid blocks and formal upscaling of the obtained results to build a coarse scale equation; ii) global solve of the upscaled coarse scale equation; and iii) reconstruction of a fine scale solution by solving local block problems on a dual coarse grid. By its structure NU is similar to other methods for solving multiscale elliptic problems, such as the multiscale finite element method, the multiscale mixed finite element method, the numerical subgrid upscaling method, heterogeneous mul-tiscale method, and the multiscale finite volume method. The difference with those methods is in the way the coarse scale equation is build and solved, and in the way the fine scale solution is reconstructed. Essential components of the presented here NU approach are the formal homogenization in the coarse blocks and the usage of so called multipoint flux approximation method, MPFA. Unlike the usual usage of MPFA as a discretization method for single scale elliptic problems with tensor discontinuous coefficients, we consider its usage as a part of a numerical upscaling approach. An aim of this paper is to compare the performance of NU with the one of MsFEM for ceratin multiscale problems. In particular, it is shown that the resonance effect, which limits the application of the Multiscale FEM, does not appear, or it is significantly relaxed, when the presented here numerical upscaling approach is applied.

Classical solutions of mixed problems for first order partial functional differential equations in several independent variables are approximated by solutions of an Euler-type difference problem. The mesh for the approximate solutions is obtained by the numerical solution of equations of bicharacteristics. The convergence of explicit difference schemes is proved by means of consistency and stability arguments. It is assumed that the given functions satisfy the nonlinear estimates of the Perron type. Differential systems with deviated variables and differential integral systems can be obtained from the general model by specializing the given operators.

In this paper we introduce an adaptive method for the numerical solu-tion of the Pocklington integro-differential equation with exact kernel for the current induced in a smoothly curved thin wire antenna. The hp-adaptive technique is based on the representation of the discrete solution, which is expanded in a piecewise p-hierarchical basis. The key element in the strategy is an element-by-element criterion that controls the h- or p-refinement. Numerical results demonstrate both the simplicity and efficiency of the approach

A numerical upscaling approach, NU, for solving multiscale elliptic problems is discussed. The main components of this NU are: i) local solve of aux-iliary problems in grid blocks and formal upscaling of the obtained results to build a coarse scale equation; ii) global solve of the upscaled coarse scale equation; and iii) reconstruction of a fine scale solution by solving local block problems on a dual coarse grid. By its structure NU is similar to other methods for solving multiscale elliptic problems, such as the multiscale finite element method, the multiscale mixed finite element method, the numerical subgrid upscaling method, heterogeneous mul-tiscale method, and the multiscale finite volume method. The difference with those methods is in the way the coarse scale equation is build and solved, and in the way the fine scale solution is reconstructed. Essential components of the presented here NU approach are the formal homogenization in the coarse blocks and the usage of so called multipoint flux approximation method, MPFA. Unlike the usual usage of MPFA as a discretization method for single scale elliptic problems with tensor discontinuous coefficients, we consider its usage as a part of a numerical upscaling approach. An aim of this paper is to compare the performance of NU with the one of MsFEM for ceratin multiscale problems. In particular, it is shown that the resonance effect, which limits the application of the Multiscale FEM, does not appear, or it is significantly relaxed, when the presented here numerical upscaling approach is applied

In this paper we develop a new method to find a numerical solution for the system of non-linear Volterra integro-differential equations (SNVE). To this end, we present our method based on the matrix form of SNVE. The corresponding unknown coefficients of our method have been determined by using the computational aspects of matrices. Finally the accuracy of the method has been verified by presenting some numerical computations.

The problem of approximate solution of severely ill-posed problems given in the form of linear operator equations of the first kind with approximately known right-hand sides was considered. We have studied a strategy for solving this type of problems, which consists in combinating of Morozov’s discrepancy principle and a finite-dimensional version of the Tikhonov regularization. It is shown that this combination provides an optimal order of accuracy on source sets