NUMERICAL GAS-DYNAMIC SIMULATION OF STRONG EXPLOSIONS IN EXPONENTIAL ATMOSPHERE

V.A. Zhmaylo, V.V. Zmushko, F.A. Pletenev, I.D. Sofronov
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 3-9.

      Numerical research on 2D problem of single explosion and 3D problem of two explosions at the same height in exponential atmosphere is described. The research was carried out within das-dynamic models with numerical approaches and within a modified method of infinitely thin shock layer. The computational results produced with different models both for 2D and 3D problems are compared.

N.Yu. Bakaev, V.N. Mikhailov
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 10-14.

      A new approach to efficient numerical solution of stiff systems of non-linear differential equations is suggested. A principle of multiplication and identification of the main part of the mapping that predetermines this stiff system formulated in this work lies in the basis of this approach.

A.V. Zabrodin, N.N. Chentsov
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 15-17.

      A classical problem of two flying-away flat pistons separated with a volume of highly-compressed ideal gas with γ isentrope index is under consideration. Approximated formulas to find final energy of light pistons and gas velocity distribution as functions of Eulerian and Lagrangian coordinates are suggested; their relation with Sedov similarity solutions for flat cylindrical and spherical flows is shown. Computational difficulties for the ballistic problem with difference methods at γ ≈ 1 are analyzed; the volume of the model is estimated from below.

S.B. Serov
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 18-21.

      Algorithms of numerical solution of 3D kinetic equation written in spherical and cylindrical coordinate systems and in a fixed basis momentum space Ω are described.

V.I. Ilin
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 22-27.

      The article is the finalizing one in a series of three devoted to methods of incomplete factorization to solve mesh linear equations. It discusses the issues of theoretical substantiation of methods, in particular the issues of no degeneracy of matrices, of the iteration processes convergence, of the evaluation of the iteration convergence velocity. The analysis of experimental data is provided and practical recommendations on implementation of various algorithms are given.

Yu.A. Skvortsov
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 28-31.

      A method of numerical simulation of non-unidimensional non-stationary flows of incompressible non-viscous gas in the air-gas channel of power units based on implementation of direct replacement of local approximating functions into the integral form of the preservation equations for the element of the computational domain is described. The performance of the method is shown using the computational example of 2D non-stationary gas flow in axisymmetric ring channel of variable area.

O.B. Khayrullina
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 32-39.

      A finite-difference method to construct stationary subsonic flows of the ideal gas in axisymmetric simply-connected channels of complex geometries is developed. The method is based on implementation of nonorthogonal optimum curvilinear meshes and allows calculating vortex flows with closed-loop lines of the flow. The results of numerical computations are provided to demonstrate the performance of the method. It is shown that in case of a large variation range of the Mach number when calculating subsonic flows we need to account for the compressibility of the gas.

A.M. Afonin, A.Yu. Pautkin, A.S. Roshal
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 40-43.

      A new efficient method for calculations on quasi-stationary electromagnetic field created by a bunch of charged particles that occupies a relatively small part of the area under investigation is proposed. To do the computations of the potentials and the fields, we used the method of partition followed by Fourier series expansion by two coordinates and implementation of spline-approximation by the third one. An analytical solution of auxiliary tasks followed by joining numerical and analytical solutions at the boundary of the bunch is used in the areas not occupied by the bunch. The method has high accuracy and is cost-effective as it does not use finite-difference approximation of derivatives and the computational grid is introduced only in the area of the bunch. The results of comparison between the analytical and numerical solutions for several benchmark problems are provided.

B.L. Voronin, V.A. Zhmaylo, A.G. Kozub, S.I. Skrypnik, I.D. Sofronov
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 44-50.

      The computational results on the problem of two strong explosions of the same yield simultaneously in the exponential atmospheres at the same height at a certain distance from each other are described. 2D and 3D options of the problem (depending on the distance) are considered. The computations were carried out using RAMZES software package based on Eulerian-Lagrangian computational method for multidimensional non-stationary gas-dynamic problems. A distribution pattern of gas-dynamic parameters in the area covered with a shock wave was produced, as well as its development with time. There was noted the preservation of the 3D character of the perturbed area shape in the respective computations at all the time range of the computation starting with the collision moment of the shock waves spread from the blasts.

A.V. Zabrodin, L.A. Pliner
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 51-55.

      The explosive in the problem under consideration is initiated by the contact with the impactor that flies at a constant velocity. Specific features of the on-going processes are reproduced in the numerical experiment. A formula for approximated estimation from above of the velocity of piston flying that gives good agreement with computational data is produced.

N.Yu. Bakaev
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 56-61.

      A stability theory of additive difference schemes in Banach spaces is described. The produced estimations of stability are the generalization of the results produced by the author on the stability of non-split difference schemes. And the components of the initial operator splitting are not supposed to be commutating, but we need them to be close to the commutating in the sense specified in the paper. The results of the work can be applied to the additive schemes both with constant and with alternating transition operator. Stability estimations for additive schemes with perturbed operators are constructed.

A.F. Nikiforov, V.G. Novikov, S.K. Trukhanov, V.B. Uvarov
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 62-74.

      The tables of pressure, internal energy and average ionization degree are provided for aluminum in a wide range of temperature and density variation (4 eV ≤4 keV; 0.01 ≤ ρ≤ 100 g/cm³) when calculated using the modified Hartree-Fock-Slatter model. The produced data were analyzed, Hugoniots for uniform and porous aluminum were constructed, and the results were compared with experimental data and computation results by other authors. The work is the first one in the series of efforts on functional completion of the THERMOS data base developed in the M.V. Keldysh Institute of Applied Mathematics of the Academy of Science of the USSR.

A.A. Solonenkov
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 75-78.

      General description of the central part of the software of the film preparation system at the non-homogeneous computer system is provided. Graphic data are collected within centralized system output. The archive structure of the system is described. A microfilming device is used for the output of vector images on the film.

L.V. Nesterenko, D.M. Obuvalin, S.V. Sivolgin
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 79-83.

      A technology to work with external storage located at the same conceptual level with widely spread imperative programming language is described. From the user’s point of view it has independent properties from the language used and from the operational environment. LINGO complex comprises four basic components: a compiler of declarative statements; pre-compiler(s) of the access operators of the programming language(s); an executive system and an administrator. User’s interface of LINGO system is represented by data access language, which is submerged into the user’s programming language and is processed by the pre-compiler. The executive system and the administrator are hidden from the level of the applied programming. They are meant for the data access, dynamic connection of the declarers with the copies of data, program operation synchronization with the data being separated.

I.A. Bashkirtsev
VANT. Ser.: Mat. Mod. Fiz. Rroc. 1990. No 3. P. 84-89.

      A problem of gas seepage into the vacuum in case of spherical and cylindrical symmetry is considered. The convergence rate of characteristic series was experimentally studied there. Acceleration methods for convergence of series were implemented; the values of gas dynamic parameters were produced with high accuracy both for the case of the cavity implosion and for the case of the gas seepage from the ball (cyllinder). The provided results can be used as references when solving gas dynamic problems with difference methods.