Published in Sarov (Arzamas-16), Nizhegorodskaya oblast
NUCLEAR CENTER -
ALL-RUSSIAN RESEARCH INSTITUTE
OF EXPERIMENTAL PHYSICS
Русский | English
Issue No 4, 2016
DSn METHOD WITH TVD RECONSTRUCTION AND SYNTHETIC P1SA ITERATION ACCELERATION FOR NUMERICAL SOLUTION OF THE TWO-DIMENSIONAL THERMAL RADIATION TRANSPORT EQUATION IN THE AXIALLY SYMMETRIC RZ GEOMETRY
A. D. Gadzhiev, D. A. Koshutin, A. A. ShestakovDevelopment of an efficient numerical solver for the thermal radiation transport equation is a challenging problem. We consider a new technique for numerical solution of the two-dimensional thermal radiation transport equation. The idea of the new approach is to use the discrete ordinate method with TVD reconstruction for solving the kinetic equation and a synthetic method based on a P1-approximation for accelerated convergence of iterations in joint solution of the transport and energy equations.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 4. P. 3-19.
Key words: radiation transport, TVD reconstruction, iteration method.
|PRACTICAL STABILITY CONDITIONS FOR A DIFFERENCE SCHEME FOR SOLVING THE HEAT TRANSFER EQUATION WITH DOUBLE THERMAL CONDUCTIVITY REMAPPING
Yu. A. Bondarenko, A. A. Gorbunov, B. P. TikhomirovIn some codes, the nonlinear heat transfer equation is solved numerically using implicit conditionally stable difference schemes with a fixed maximum permissible number of iterations with respect to the nonlinearity of thermal conductivity having the value of 1, 2 or 3. In this paper, the von Neumann-Fourier method is combined with the Newton numerical method for initial data stability analysis of a scheme with double remapping of the field of thermal conductivities. From the stability condition we derived a time step constraint suitable for practical simulations. It is shown that the dependence of the time step on the mesh spacing is not quadratic.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 4. P. 20-25.
Key words: radiant heat transfer, implicit difference scheme, thermal conductivity, initial data stability, von Neumann-Fourier method, iterative Newton method, time step, test simulations.
|STRUCTURE OF THREE-DIMENSIONAL STATIONARY ISOBARIC DOUBLE WAVES IN AN IDEAL INCOMPRESSIBLE FLUID. PART 2. SOLUTION OF REDUCED SYSTEM AND LOCAL CLASSIFICATION OF DOUBLE WAVES
V. E. ShemarulinStudies of the structure of isobaric three-dimensional stationary waves in an ideal incompressible fluid are coming to an end. In the first part of the work, the system of functional equations implicitly defining these waves was reduced to an equivalent system, which is more convenient for research. The present part shows that the solutions of the reduced system have a simple geometric structure and provides an explicit description of the diversity of local solutions of this system. As a result, all isobaric three-dimensional stationary double waves are proven to be clusters of flow regions of thee major types: shear, conical and tangential. Shear flows are the well known flows in parallel planes; for each plane, in parallel right lines. Conical flows are the flows along sets of tangent semiplanes of arbitrary convex conical surfaces; for each semiplane, in right lines parallel to the generatrix of the conical surface belonging to the semiplane. Tangential flows are the flows along sets of tangent surfaces (surfaces generated by right tangents of arbitrary space curves) governed by certain convexity conditions (ensuring that there are no intersections of tangent semiplanes); for each semiplane, in right lines parallel to the generatrix of the tangent surface (tangent to the curve) belonging to the semiplane. A similar explicit description of the structure of isobaric three-dimensional steady-state simple waves and two-dimensional non-steady-state isobaric flows has been given by the author before. Thus, the present work finalizes the local classification of isobaric flows of an ideal incompressible fluid for the case of three independent variables: three space variables or two space variables and time.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 4. P. 26-39.
Key words: isobaric flows, three-dimensional stationary double waves, local flow classification, shear, conical and tangential flows.
|COMPARISON OF TWO SWEEP PARALLELING METHODS FOR HYBRID COMPUTERS WITH GRAPHICS PROCESSING UNITS
A. A. Fedorov, A. N. BykovThe paper describes two methods (the parallel pipeline and the Yanenko method) for solving a system of linear algebraic equations with a tridiagonal matrix on GPU-accelerated parallel computes. Specific features of their implementation both on parallel computers having no accelerators and on hybrid computers are discussed. Arithmetic complexity of the Yanenko method is analyzed, and results of numerical scalability experiments are reported.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 4. P. 40-50.
Key words: system of linear algebraic equations with tridiagonal matrix, sweep method, parallel pipeline, Yanenko method, GPU-accelerated parallel computers, code RAMZES-KP.
|PROTON FAST IGNITION OF A CYLINDRICAL DT TARGET ENCLOSED IN A STATIONARY HEAT-INSULATED SHELL
K. V. Khishchenko, A. A. CharakhchianThe paper considers a two-dimensional axially-symmetric problem of fast ignition of a pre-compressed cylindrical target with a DT mixture enclosed in a stationary heat-insulated shell. The target is ignited on the end by a proton beam, the intensity of which does not depend on the radial coordinate. The difference from the one-dimensional problem is that the fusion alpha particles and the self-radiation of plasma freely emerge from the fuel through the side boundary with the shell. The threshold ignition energy for the mixture densities of 22 and 110 g/cm3 is shown to be about one tenth of that in the known case without the shell and with a beam radius much smaller than the target radius. As compared to the problem of interest, the quasi-one-dimensional model developed earlier understates the threshold ignition energy with respect to the radius by a factor of four for the given beam intensity.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 4. P. 51-59.
Key words: cylindrical ICF target, fast ignition, ignition energy.
|INVESTIGATION OF THE INFLUENCE OF AUSTENITIC STEEL SURFACING ON CRACK OPENING IN A PEARLITIC STEEL PIPELINE DN800
D. A. KuzminThe object of the study is a straight section of the main circulation pipeline from pearlitic steel with austenitic steel cladding. These steels have different linear expansion coefficients. The influence of the cladding on the opening of a through crack in the pipeline is studied. The size of crack opening and the shape of crack faces are determined for various crack lengths. A comparative analysis of the crack opening area for the pipeline with cladding and without it is made.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 4. P. 60-65.
Key words: main circulation pipeline, austenitic cladding, crack opening.
|A TECHNIQUE FOR CALCULATING RADIATION LOSSES IN THE INITIAL SECTION OF A BENT SINGLE-MODE OPTICAL FIBER
Yu. V. Malykh, V. V. ShubinA technique for calculating optical radiation losses in the initial section of a bent single-mode optical fiber as a function of the bending angle is presented. Known techniques and formulas for calculating radiation losses in bent optical fibers are discussed in terms of their applications. As far as the techniques considered cannot describe the losses in the initial bend section, we have developed a new approach. The resulting technique was used to calculate radiation losses as a function of the bending angle for an SMF-28 optical fiber. The calculated results are compared with experimental data.
VANT. Ser.: Mat. Mod. Fiz. Proc. 2016. No 4. P. 66-77.
Key words: bend of single-mode optical fiber, transition radiation losses, bending radiation losses, conformal mapping method, length and angle of bending.
|[ Back ]