Volume 7, Issue 3

A new approach has been developed for solving the tide dynamics problem, based on the splitting methods and the optimal control theory. We first apply the classical splitting method for solving the problem. Then we use the optimal control theory approaches to the system of equations obtained after the splitting. An optimal control problem has been formulated for realizing a splitting method step. We prove that the optimal control problem is well-posed and we propose an iterative process of the minimization problem. The results of the numerical experiments are presented.

A priory estimates of the stability in the sense of the initial data of the difference scheme approximating weakly compressible liquid equations in the Riemann invariants have been obtained. These estimates have been proved without any assumptions about the properties of the solution of the differential problem and depend only on the behavior of the initial conditions. As distinct from linear problems, the obtained estimates of stability in the general case exist only for a finite instant of time t≤t_0. In particular, this is confirmed by the fact, that nonfulfilment of these stability conditions lead to the appearance of supersonic flows or domains with large gradients. The questions of uniqueness and convergence of the difference solution are considered also. The results of the computating experiment confirming the theoretical conclusions are given.

This paper presents methods to compute integrals of the Jacobi polynomials by the representation in terms of the Bernstein — B´ezier basis. We do this because the integration of the Bernstein — B´ezier form simply corresponds to applying the de Casteljau algorithm in an easy way. Formulas for the definite integral of the weighted Bernstein polynomials are also presented. Bases transformations are used. In this paper, the methods of integration enable us to gain from the properties of the Jacobi and Bernstein bases.

In the proposed method, the variation of displacement in each time step is assumed to be a fourth order polynomial in time and its five unknown coefficients are calculated based on: two initial conditions from the previous time step; satisfying the equation of motion at both ends of the time step; and the zero weighted residual within the time step. This method is non-dissipative and its dispersion is considerably less than in other popular methods. The stability of the method shows that the critical time step is more than twice of that for the linear acceleration method and its convergence is of fourth order.

Finite-order weights have been introduced in recent years to describe the often occurring situation that multivariate integrands can be approximated by a sum of functions each depending only on a small subset of the variables. The aim of this paper is to demonstrate the danger of relying on this structure when designing lattice integration rules, if the true integrand has components lying outside the assumed finiteorder function space. It does this by proving, for weights of order two, the existence of 3-dimensional lattice integration rules for which the worst case error is of order O(N¯½), where N is the number of points, yet for which there exists a smooth 3- dimensional integrand for which the integration rule does not converge.

The approximation properties of a fully discrete projection method for Symm’s integral equation with a infinite smooth boundary have been investigated. For the method, error bounds have been found in the metric of Sobolev’s spaces. The method turns out to be more accurate compared to the fully discrete collocation method known before.

A new fourth order accurate centered finite difference scheme for the solution of hyperbolic conservation laws is presented. A technique of making the fourth order scheme TVD is presented. The resulting scheme can avoid spurious oscillations and preserve fourth order accuracy in smooth parts. We discuss the extension of the TVD scheme to the nonlinear scalar hyperbolic conservation laws. For nonlinear systems, the TVD constraint is applied by solving shallow water equations. Then, we propose to use this fourth order flux as a building block in spatially fifth order weighted essentially non-oscillatory (WENO) schemes. The numerical solution is advanced in time by the third order TVD Runge — Kutta method. The performance of the scheme is assessed by solving test problems. The numerical results are presented and compared to the exact solutions and other methods.