Since 1978
Published in Sarov (Arzamas-16), Nizhegorodskaya oblast

RUSSIAN FEDERAL
NUCLEAR CENTER -
ALL-RUSSIAN RESEARCH INSTITUTE
OF EXPERIMENTAL PHYSICS
 
 Русский |  English
ABOUT EDITORIAL BOARD PUBLICATION ETHICS RULES FOR AUTHORS AUTHORS ARCHIVE MOST RECENT ISSUE IN NEXT ISSUE PAPER OF THE YEAR



Issue No 1, 2024

NUMERICAL SOLUTION OF THE RADIATIVE HEAT TRANSFER EQUATION IN 2D PROBLEMS WITH PHASE TRANSITIONS ON ADAPTIVE MESHES

A. M. Mustafin, N. N. Veselova, S. N. Lebedev
VANT. Ser.: Mat. Mod. Fiz. Proc. 2024. No 1. P. 4-17.

The paper discusses the development of a numerical technique for solving a 2D radiative heat transfer equation with account for phase transitions on adaptive difference meshes with a ROMB scheme, including the construction of the adaptive mesh and the efficiency testing of the developed technique. The problems of heat transfer during phase transitions, referred to as Stephan´s type problems, are solved by a once-through method based on the use of a difference analog of the δ-function. A ROMB scheme on an adaptively embedded mesh is obtained and an algorithm for constructing an adaptive difference mesh is described. Provided are numerical solutions for two problems: heat wave propagation in a heterogeneous region and a thermal beam impact on ice. The first clearly shows advantages of mesh adaptation near the front of the heat wave in 2D space. Both demonstrate gains from the use of mesh adaptation in combination with once-through calculation before applying traditional approaches to Stephan´s type problems.


Key words: Stephan´s problem, ROMB difference scheme, adaptive difference mesh, 2D radiative heat transfer equation, phase transition, heat wave.

STOCHASTIC METHOD OF NUMERICAL SIMULATION OF RADIATION BELTS OF THE EARTH

A. N. Zalialov, N. V. Ivanov, A. E. Shirokov
VANT. Ser.: Mat. Mod. Fiz. Proc. 2024. No 1. P. 18-29.

The paper describe a stochastic method of numerical simulation of the dynamics of relativistic electrons in radiation belts of the Earth that uses approximation of the leading center that, opposite other approaches, allows describing relaxation of relatively quick perturbations in the radiation belts of the Earth (for example, during geomagnetic storm) in the regions of low orbits. A stationary off-centered dipole is used as a physical–mathematical model of the geomagnetic field.
The solution of the kinetic equation in the Fokker - Plank approximation in this method is found using a Monte-Carlo method basing on the simulation of the motion of relativistic electrons along the geomagnetic line with the account for their interaction with the atmosphere. A method of splitting by physical processes into three stages is used to construct the trajectory. The motion of the electron under the effect of the magnetic field with the account for the deceleration force is computed at the first stage, electron scattering is simulated at the second stage, and its radial diffusion is simulated at the third stage.
The results of the numerical solution of the test problems is compared with the results of the analytical and difference solutions.


Key words: relativistic electrons, radiation belt of the Earth, Monte-Carlo method, radial diffusion.

SPECIFIC FEATURES OF DETONATION PROPAGATION IN SMALL CROSS-SECTION CHANNELS FOR THE PETN-BASED COMPOSITION

V. B. Titova, N. A. Volodina, M. O. Shirshova, E. V. Khaldeev, Y. M. Sustaeva
VANT. Ser.: Mat. Mod. Fiz. Proc. 2024. No 1. P. 30-36.

The works to study the processes of detonation wave propagation in the channels equipped with PETN-based (pentraerythritoltetranitrate) composition are currently underway to upgrade the systems of detonation initiation in experimental setups with explosives.
This paper studies specific features of the detonation wave propagation in the small-cross-section channels using model experiments with numerical simulation in the ideal and kinematic approximations. It shows the necessity to account for detonation kinetics in such computations. As a result of the computational analysis, we obtained a characteristic picture of the detonation propagation when the detonation wave follows the shadow zone, and the size of the shadow zone region is estimated after the detonation wave passes the turn of the channel.


Key words: detonation kinetics, PETN, channels interconnections, numerical simulation.

ACCOUNT FOR KINETIC EFFECTS IN APPROXIMATION OF RADIATING HEAT TRANSFER TO COMPUTE RADIATION TRANSFER PROBLEMS

D. A. Koshutin, A. A. Shestakov
VANT. Ser.: Mat. Mod. Fiz. Proc. 2024. No 1. P. 37-41.

One of the most difficult problems of radiation gas dynamics is solution of a multi-dimensional equation of radiative heat transfer, the difficulty of which mostly comes from the large-dimensional character of the phase space. Because of this, they implement different approximations that simplify the radiation transfer problem. These approximations comprise, for example, approximation of radiative heat transfer. Here, complex enough problems of selecting an effective iteration method and construction of a monotonous scheme of upgraded accuracy for the kinetic equation are eliminated. When we implement approximation of radiative heat transfer, different correction factors are used. The need in accurate correction factors in approximation of the radiative heat transfer stimulated development of new mathematical models.
In 2016, M. Yu. Kozmanov, D. A. Koshutin and A. A. Shestakov suggested replacing a system of 2D transfer equations with the equations of quasi-thermal conductivity, where the solutions are obtained by parabolic equations. This is a generalization of the quasi-transfer method offered by M. Yu. Kozmanov and N. G. Karlykhanov. To solve a 2D system of parabolic equations of quasi-heat transfer, they used the method splitting by directions and the difference scheme of the second order of approximation. In the approximation of quasi-heat transfer, a parabolic system is solved, for which you can construct a monotonous difference scheme of the second order of approximation.
The paper describes algorithms proposed in 2016, but not previously published.
The ideas of quasi-thermal conductivity were further developed in the works by N. G. Karlykhanov, where averaging of multi-group equation of transfer by the spectrum of photons was additionally applied.


Key words: radiative heat transfer, account for kinetic effects.

VERIFICATION OF THE EQUATIONS OF STATE FIR THE EXPLOSION PRODUCTS OF THE PLASTICIZED HMX USING THE EXPERIMENTAL RESULTS ON LINERS ACCELERATION

V. B. Titova, N. A. Volodina, M. O. Shirshova, M. N. Kityukhina, E. N. Bogdanov, A. A. Stanovov
VANT. Ser.: Mat. Mod. Fiz. Proc. 2024. No 1. P. 42-47.

The work describes the experimental results where continuous registration was carried out for the motion of liners made of aluminum and the tungsten-nickel-iron alloy thrown with the explosion products of the plasticized HMX. The computational analysis is provided for the experimental results on the research on the propelling properties of the plasticized HMX in order to verify the equations of state of its explosion products. Verification of the equations of state of the explosion products was carried out for three forms of the equations of state relying of the experimental data. Involvement of modern experimental research methods allowed for more high-quality verification of the equations of state of the explosion products.


Key words: explosion products, detonation, equation of state, kinetics, numerical simulation.

50 YEARS OF FLECK PROBLEMS

A. A. Shestakov
VANT. Ser.: Mat. Mod. Fiz. Proc. 2024. No 1. P. 48-56.

In 1971, specialists on radiation transfer from Livermore National Lab (LLNA) of the USA suggested problems that with time became classical tests for radiation transfer. In all problems, a flat layer is heated with the radiation flux that corresponds to the temperature of 1 keV. The basic advantages of these problems are:
      – the simplicity of the setup of the boundary and initial conditions that allows computing these problems in different approximations;
      – simple analytical formulas for spectral runs;
      – a wide range of implementation of spectral coefficient of absorption;
      – simulation of disruptors of the spectral coefficient разрывов of absorption in space and energy spectrum;
      – the simplest equation of state of the material.
The paper is devoted to the fifty-year anniversary of Fleck tests in the method of numerical solution of heat radiation transfer. The paper provides the setups and numerical computations of two modifications of these problems.


Key words: heat radiation transfer, Fleck problems.

NUMERICAL SIMULATION ON THE BASIS OF LINEARIZED SYSTEM OF GAS-DYNAMIC EQUATIONS OF RISING SWIRL FLOWS WITH THE ACCOUNT FOR CORIOLIS FORCES

A. A. Bugaenko, I. Yu. Krutova
VANT. Ser.: Mat. Mod. Fiz. Proc. 2024. No 1. P. 57-63.

A system of gas-dynamic equations is linearized with the account for Coriolis forces using an accurate solution (uniform rest). Origination and development of rising swirl flows are numerically simulated using the gas blow up the pipe. The results of numerical computations of the gas-dynamic parameters for average by intensity cyclones are described.


Key words: a system of gas-dynamic equations, Coriolis force, linearization, numerical simulation, gas blow.

ALGORITHMS AND SPECIFIC FEATURES OF PROGRAM REALIZATION OF SEMI-EMPIRICAL MODEL OF WIDE-RANGE "ROSA-MFI" EQUATION OF STATE IN "URS-OF" PACKAGE

A. S. Danilov, D. G. Gordeev, O. N. Shumilina, I. N. Arapov
VANT. Ser.: Mat. Mod. Fiz. Proc. 2024. No 1. P. 64-77.

The paper briefly describes algorithms that realize the computation of thermodynamic functions by semi-empirical model of wide-range ROSA-MFI equations of state and specific features of their program realization in the library of programs of UFS-OF package. The ROSA-MFI model is the development of the earlier created ROSA-MI model and it additionally accounts for the changes in the thermodynamic functions at a solid - a liquid phase transition (melting). Possibilities of the model in the description of thermodynamic functions in the region of melting are shown using the example of the copper equation of state.
Program realization of the ROSA-MFI model is done in Fortran 90 language using vectorized structures in multi-architectural reference format. Using specially developed pro-processor, equivalent transformations of this format onto programs are done that account for the specific features of not-vectorized and vectorized computations and the computations using the CUDA architecture.


Key words: URS-OF program library, algorithms of the ROSA-MFI model, algorithms to compute thermodynamic functions at phase transitions, melting, evaporation, multi-version reference file, parallelization efficiency.

 
 
© FSUE "RFNC-VNIIEF", 2000-2024