A SERVICE SYSTEM FOR THE 1-D OK COMPLEX

E. G. Voronov, M. I. Kaplunov, V. G. Podvalny, S. K. Rebrov, V. I. Filatov
Vant. Ser. Metodiki i Programmy Chislennogo Resheniya Zadach Matematicheskoy Fiziki. 1978. No 2. P. 3-12.

      The paper contains a general description of a program package service system for solving computational physics problems referred to as 1-D Complex (OK) where the application programs are written in FORTRAN. We define some non-FORTRAN elements, expand ing the application programming capabilities within OK Complex, along with their functional role and command language facilities for introducing these elements. The archive storage possibilities are briefly described. A general OK Complex operation scheme is presented.

V. E. Troshchiev, V. F. Yudintsev
Vant. Ser. Metodiki i Programmy Chislennogo Resheniya Zadach Matematicheskoy Fiziki. 1978. No 2. P. 13-16.

      An averaged continuous correction method is developed which is used to obtain a faster temperature convergence of simple iterations when solving a spectral thermal-radiation transport equation and a media energy equation.
      Two linear models are constructed for a nonlinear spectral transport equation which were used to investigate theoretically the convergence of both simple and correction iterations. The major study results include formulae for averaging spectral radiation runs, which cause the correction method to converge and to give a substantial speedup of the iteration convergence.

V. A. Bychenkov, V. V. Gadjieva
Vant. Ser. Metodiki i Programmy Chislennogo Resheniya Zadach Matematicheskoy Fiziki. 1978. No 2. P. 17-22.

      A computational medium model and numerical method for solving 2-D time-dependent problems, taking into account elastic/plastic properties and material destruction, are proposed.
       Several simple assumptions, introduced in the model, permit to simulate the flows of media having different strength properties (porous, crushed, fractured). The model studies two nontrivial parameters, that is pore and fracture specific volumes. The governing eguations are presented in a differential form.
       The numerical method is based on an explicit difference scheme approximating original differential eguations on rectangudar meshes in Lagrangian coordinates.
       An example is included which shows simultaneous burst of two charges in a destructable medium.

V. V. Bashurov, V. A. Svidinsky, T.F. Kryukova
Vant. Ser. Metodiki i Programmy Chislennogo Resheniya Zadach Matematicheskoy Fiziki. 1978. No 2. P. 23-29.

      A finite-difference scheme is considered for solving thermoelasticity equations based on differential equations written in mixed Eulerian-Lagrangian coordinates. An algorithm is given to evaluate contact boundaries and cuts, typical for thermoelasticity problems. The effect of some equation terms on the stress deviator is analyzed, which is then demonstrated using the simplest types of motion.
      The problem describing a hemisphere motion due to energy release within the medium is investigated. The interaction between the hemisphere and the rigid base is also studied. A numerical analysis of a problem illustrating the dependence of hemisphere average rebound velocity on its size and Poisson ratio is carried out.

V. A. Batalov, S. M. Bakhrakh, B. A. Klopov, N. P. Kovalyov, E. E. Meshkov, V. N. Mokhov, I. N. Pavlusha, A. V. Pevnitsky, V. A. Sarayev, V. P. Sevastyanov, V. I. Tarasov, Yu. V. Yanilkin, M. S. Samigulin
Vant. Ser. Metodiki i Programmy Chislennogo Resheniya Zadach Matematicheskoy Fiziki. 1978. No 2. P. 30-36.

      In order to study and numerically simulate the mechanics of an underground explosion a continuum model is considered, taking into account compressibility, elastic/plastic properties, initial porosity, destruction, heat conductivity, and turbulent effects. Appropriate algorithms and programs are written to solve both 1-D and 2-D (axisymmetric) problems. General laws, defining a containment cavity formation including, in particular, the medium parameters impact on the cavity size are estimated. Burst product and ground motion due to the explosion with nonspherical initial geometry are studied.
      A comparison between 2-D calculation and the U.S. Marvel experiment is made. Some laboratory experiment and 2-D calculations results are given, which show a channel collapse effect occurring in the case of explosion with such initial geometry.