Published in Sarov (Arzamas-16), Nizhegorodskaya oblast
NUCLEAR CENTER -
ALL-RUSSIAN RESEARCH INSTITUTE
OF EXPERIMENTAL PHYSICS
Русский | English
Issue No 2, 2016
INVESTIGATION OF THE ITERATION CONVERGENCE RATE OF THE VDEPF METHOD USED FOR JOINTLY SOLVING THE ENERGY AND RADIATION TRANSPORT EQUATIONS
A. A. ShestakovThe application of implicit schemes for the mathematical modeling of the thermal radiation transport requires studying the convergence problem of the iterative process of solving a system of photon and energy transport equations. The paper considered some modifications to the method of identifying a diagonal element in the corrector form for iteratively solving a nonlinear system consisting of the photon transport equation and energy transport equation. The iteration convergence rate is theoretically investigated.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 2. P. 3-16.
Key words: radiation transport, iteration convergence rate.
|NONSTATIONARY PROBLEMS FOR TESTING ELASTO-PLASTICITY CODES
A. A. Krayukhin, V. A. Svidinskii, A. L. Stadnik, Yu. V. YanilkinThe paper describes three nonstationary test problems for the verification of elasto-plasticity codes. The first problem consists in finding the velocity of a moving longitudinal wave in elastic media of different configurations, such as an unbounded medium, a plate, and a bar. The second one is the problem of accelerating a system of elastic components by a pressure pulse, or an impactor. Its analytical solution has been found for the system of three components in the form of bars, or plates. In the third problem a plane elasto-plastic wave propagates in an unbounded medium. The problem has been analytically studied in details. Results of numerical simulations with the EGAK code are presented for the problems above.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 2. P. 17-30.
Key words: verification, test problems, d’Alembert method, longitudinal wave, the EGAK code, elasto-plastic wave.
|APPLICATION OF THE ALGORITHM OF CORRECTING A GRID WITH RESPECT TO THE SURFACE OF REVOLUTION
O. V. UshakovaThe paper considers the application of the algorithm of correcting a grid with respect to the surface of revolution generated by the revolution through 180° around its axis of a plane curve generator consisting of straight-line segments, arcs of circles, and ellipses. The algorithm is used during the global reconstruction, or optimization of structured grids.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 2. P. 31-37.
Key words: grids, boundary nodes, correction, surfaces of revolution.
|3D NUMERICAL SIMULATION OF THE VORTEX RING GENERATION DUE TO A RISING THERMAL IN ATMOSPHERE
O. G. Sin’kova, V. P. Statsenko, Yu. V. YanilkinThe generation of a vortex ring due to a rising atmospheric thermal (light cloud) in the earth gravity field is studied numerically. 3D simulations were carried out with the EGAK code. Results of these simulations are compared with the approximate analytical solutions for the motion of a buoyant vortex ring verified using the existing experimental data.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 2. P. 38-47.
Key words: buoyant volume, field of gravity, vortex ring, turbulence, numerical simulation.
|DIRECT NUMERICAL SIMULATION OF NEAR-WALL TURBULENCE IN PLANE CHANNEL IN A WIDE RANGE OF REYNOLDS NUMBERS
D. G. Asfandiyarov, S. A. Finogenov, V. M. GolovizninThe paper presents direct numerical simulation results for the near-wall turbulence in a plane channel obtained using computational grids of 4 x 106 to 34 x 106 cells for Reynolds numbers 5 600, 13 750, 21 900. Simulations were carried out using the modified CABARET scheme. The choice of geometrical sizes for the computational domain and computational grids is justified based on the analysis of correlation dependences and spectral characteristics of flow. The simulation results were used to obtain standard parameters of the turbulent flow in plane channel, namely, the averaged velocity profile, components of turbulent stress tensor, contributions from the molecular and turbulent components of viscosity tensor, and the channel resistance coefficient. It is shown that they are in a good agreement with the data calculated for a similar problem using the pseudo-spectral scheme. All simulations were carried out on supercomputer "Lomonosov".
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 2. P. 48-58.
Key words: direct numerical simulation, scheme CABARET, turbulent flow in plane channel, viscous incompressible fluid.
|NUMERICAL SIMULATION OF IMPACTS OF STEEL IMPACTORS AND ALUMINUM TARGETS AT A VELOCITY OF 6.2 KM/S USING SPH CODE OF LOGOS SOFTWARE PACKAGE
I. V. Kovalenko, N. A. Zhavoronkova, E. I. Lipenkova, Yu. N. BukharevThe paper describes several problems of a high-speed impact and presents results of the 3D simulation for these problems using the SPH code of LOGOS software package developed at FSUE "RFNC-VNIIEF" (Sarov, N.Novgorod region). The work objective was to obtain verification data for the conditions of interaction at a velocity of 6.2 km/s between a disc-shaped target made of aluminum alloy AMg-6 and an impactor (a sphere, or a disc) made of steel EP-637 by comparing results of simulations with the earlier published data of experiments and results of simulations using other codes. Stability of the numerical simulation results obtained for the problems of impacting an Al disk by a steel spherical impactor and their satisfactory agreement with the available experimental and earlier calculated data have been demonstrated. For the shock-wave phase of a planar impact of two discs, there have been found noticeable oscillations of pressure and mass velocity relative to the expected analytical values, which can be lowered, in particular, by optimizing the artificial viscosity parameters in the computational algorithm.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 2. P. 59-65.
Key words: numerical simulation, high-speed impact, SPH code, LOGOS software package, impactor, target.
|GEOMETRICAL INTEGRATION OF 3D ISOBARIC FLOW EQUATIONS OF AN IDEAL INCOMPRESSIBLE FLUID
V. E. ShemarulinA geometrical method is used to integrate equations describing 3D isobaric flows of an ideal incompressible fluid. An equivalent system of external equations is matched with the system of differential equations. It is shown that the original differential system integrability is conditioned by expansibility of the corresponding external forms and full integrability of the distribution associated with these forms. As a result, a general solution to the 3D isobaric flow equations has been found and the geometrical cause of integrability of these equations has been identified. The same is valid for 2D equations also, however, the integration of 2D equations is trivial and not considered in the paper. Some other authors have already found the general solution to isobaric equations using traditional methods and, therefore, the integration method described here is of interest from viewpoint of methodology.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 2. P. 66-74.
Key words: ideal incompressible fluid, isobaric flows, differential forms, associated distribution, full integrability, a general solution to isobaric flow equations.
|MATHEMATICAL MODEL OF URANIUM HEXAFLUORIDE DESUBLIMATION
A. A. Orlov, A. F. Tsimbalyuk, R. V. MalyuginThe paper presents a nonstationary mathematical model of the UF6 desublimation in a vertical storage vessel regarding the convection and UF6 desublimation on a bottom wall, as well as the ellipticity of the upper and lower walls of the vessel. Simulation results are presented for the linear desublimation rate, thickness of a solid layer, temperature distribution in the vessel wall-desublimate layer system, and dynamics of filling out the 2,5 m3 cylindrical vessel with solid UF6 as functions of time.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 2. P. 75-83.
Key words: uranium hexafluoride, desublimation, mathematical model, heat-and-mass transport, nonstationary process, gas dynamics, Stephen problem, a mobile grid, SIMPLE algorithm, iteration-interpolation method.
|[ Back ]