Based on finite-difference approximations in time and a bilinear finite-element approximation in spatial variables, numerical implementations of a new iterative method with boundary condition splitting are constructed for solving the Dirichlet initial-boundary value problem for the nonstationary Stokes system. The problem is considered in a strip with a periodicity condition along it. At each iteration step of the method, the original problem splits into two much simpler boundary value problems that can be stably numerically approximated. As a result, this approach can be used to construct new effective and stable numerical methods for solving the nonstationary Stokes problem. The velocity and pressure are approximated by identical bilinear finite elements, and there is no need to satisfy the well-known difficult-to-verify Ladyzhenskaya-Brezzi-Babuska condition, as is usually required when the problem is discretized as a whole. Numerical iterative methods are constructed that are first- and second-order accurate in time and second-order accurate in space in the max norm for both velocity and pressure. The numerical methods have fairly high convergence rates corresponding to those of the original iterative method at the differential level (the error decreases approximately 7 times per iteration step). Numerical results are presented that illustrate the capabilities of the methods developed.
