Volume 37, Issue 3

The 4D variational data assimilation technique is presented for modelling the sea dynamics problems, developed at the Marchuk Institute of Numerical Mathematics of the Russian Academy of Sciences (INM RAS). The approach is based on the splitting method for the mathematical model of sea dynamics and the minimization of cost functionals related to the observation data by solving an optimality system that involves the adjoint equations and observation and background error covariances. Efficient algorithms for solving the variational data assimilation problems are presented based on iterative processes with a special choice of iterative parameters. The technique is illustrated for the Black Sea dynamics model with variational data assimilation to restore the sea surface heat fluxes.

An iterative method with the number of iterations independent of the coefficient jumps is proposed for the boundary value problem for a diffusion equation with highly varying coefficient. The method applies one solution of the Poisson equation at each step of iteration. In the present paper we extend the class of domains the iterative method is justified for.

The time-fractional phase transition problem, formulated in enthalpy form, is studied. This nonlinear problem with an unknown moving boundary includes, as an example, a mathematical model of one-phase Stefan problem with the latent heat accumulation memory. The posed problem is approximated by the backward Euler mesh scheme. The unique solvability of the mesh scheme is proved and a priori estimates for the solution are obtained. The properties of the mesh problem are studied, in particular, an estimate of movement rate for the mesh phase transition boundary is established. The proved estimate make it possible to localize the phase transition boundary and split the mesh scheme into the sum of a nonlinear problem of small algebraic dimension and a larger linear problem. This information can be used for further construction of efficient algorithms for implementing the mesh scheme. Several algorithms for implementing mesh scheme are briefly discussed.

A finite element method for a monolithic quasi-Lagrangian formulation of a fluid–porous structure interaction problem with a corrected balance of stresses on the fluid–structure interface is considered. Deformations of the elastic medium are not necessarily small and are modelled using Saint Venant–Kirchhoff (SVK) constitutive relation. The stability of the method is proved in a form of energy bound for the finite element solution.

A technique for constructing a sequence of difference schemes with the properties of a corrector and a predictor for integrating systems of the second-order ordinary differential equations is presented. The sequence of schemes begins with the explicit three-point Störmer method of the second order of approximation. Each subsequent scheme also implements the Störmer method corrected with additional terms calculated through the solution of the previous scheme. The stability of the resulting schemes and the increase in the order of convergence for the first of them are carefully substantiated. The results of calculations of the test problem are presented, confirming the increase in the order of accuracy of the constructed methods.