A novel second-order adaptive grid method for singularly perturbed convection–diffusion equations

In this paper, a singularly perturbed convection–diffusion equation is studied. At first, the original problem is transformed into a parameterized singularly perturbed Volterra integro-differential equation by using an integral transform. Then, a second-order finite difference method on an arbitrary mesh is given. The stability and local truncation error estimates of the discrete schemes are analyzed. Based on the mesh equidistribution principle and local truncation error estimation, an adaptive grid algorithm is given. In addition, in order to calculate the parameters of the transformation equation, a nonlinear unconstrained optimization problem is constructed. Numerical experiments are given to illustrate the effectiveness of our presented adaptive grid algorithm.

On the algorithmic solution of optimization problems subject to probabilistic/robust (probust) constraints

We present an adaptive grid refinement algorithm to solve probabilistic optimization problems with infinitely many random constraints. Using a bilevel approach, we iteratively aggregate inequalities that provide most information not in a geometric but in a probabilistic sense. This conceptual idea, for which a convergence proof is provided, is then adapted to an implementable algorithm. The efficiency of our approach when compared to naive methods based on uniform grid refinement is illustrated for a numerical test example as well as for a water reservoir problem with joint probabilistic filling level constraints.

Development of algorithms for constructing two-dimensional optimal boundary-adaptive grids and their software implementation

Introduction. It is noted that the use of adaptive grids in calculations makes it possible to improve the accuracy and efficiency of computational algorithms without increasing the number of nodes. This approach is especially efficient when calculating nonstationary problems. The objective of this study is the development, construction and software implementation of methods for constructing computational two-dimensional optimal boundary-adaptive grids for complex configuration regions while maintaining the specified features of the shape and boundary of the region. The application of such methods contributes to improving the accuracy, efficiency, and cost-effectiveness of computational algorithms.Materials and Methods. The problem of automatic construction of an optimal boundary-adaptive grid in a simply connected region of arbitrary geometry, topologically equivalent to a rectangle, is considered. A solution is obtained for the minimum set of input information: the boundary of the region in the physical plane and the number of points on it are given. The creation of an algorithm and a mesh generation program is based on a model of particle dynamics. This provides determining the trajectories of individual particles and studying the dynamics of their pair interaction in the system under consideration. The interior and border nodes of the grid are separated through using the mask tool, and this makes it possible to determine the speed of movement of nodes, taking into account the specifics of the problem being solved.Results. The developed methods for constructing an optimal boundary-adaptive grid of a complex geometry region provides solving the problem on automatic grid construction in two-dimensional regions of any configuration. To evaluate the results of the algorithm research, a test problem was solved, and the solution stages were visualized. The computational domain of the test problem and the operation of the function for calculating the speed of movement of interior nodes are shown in the form of figures. Visualization confirms the advantage of this meshing method, which separates the border and interior nodes.Discussion and Conclusions. The theoretical and numerical studies results are important both for the investigation of the grids qualitative properties and for the computational grid methods that provide solving numerical modeling problems efficiently and with high accuracy.

Travelling Waves for Adaptive Grid Discretizations of Reaction Diffusion Systems I: Well-Posedness

In this paper we consider a spatial discretization scheme with an adaptive grid for the Nagumo PDE. In particular, we consider a commonly used time dependent moving mesh method that aims to equidistribute the arclength of the solution under consideration. We assume that the discrete analogue of this equidistribution is strictly enforced, which allows us to reduce the effective dynamics to a scalar non-local problem with infinite range interactions. We show that this reduced problem is well-posed and obtain useful estimates on the resulting nonlinearities. In the sequel papers (Hupkes and Van Vleck in Travelling waves for adaptive grid discretizations of reaction diffusion systems II: linear theory; Travelling waves for adaptive grid discretizations of reaction diffusion systems III: nonlinear theory) we use these estimates to show that travelling waves persist under these adaptive spatial discretizations.

An adaptive moving grid finite difference method

The grid generation is very crucial for the accuracy of the numerical solution of PDEs, especially for problems with very rapid variations or sharp layers, such as shock waves, wing leading and trailing edges, regions of separation, and boundary layers. The adaptive grid generation is an iterative approach to accommodate these complex structures. In this paper, we introduce a deformation based adaptive grid generation method, in which a differentiable and invertible transformation from computational domain to physical domain is constructed such that the cell volume (Jacobian determinant) of the new grid is equal to a prescribed monitor function. A vector field is obtained by solving the div-curl system and can be used to move the grids to the desired locations. By computing the inverse of Jacobian, any deformed grids can also be transformed back to the uniform grid. Several numerical results in two dimensions are presented. Some applications in image registration are discussed.