Acesso aberto Acesso aberto  Acesso é fechado Acesso está concedido  Acesso é fechado Somente assinantes

Volume 56, Nº 9 (2016)

Article

Iteratively regularized methods for irregular nonlinear operator equations with a normally solvable derivative at the solution

Kokurin M.

Resumo

A group of iteratively regularized methods of Gauss–Newton type for solving irregular nonlinear equations with smooth operators in a Hilbert space under the condition of normal solvability of the derivative of the operator at the solution is considered. A priori and a posteriori methods for termination of iterations are studied, and estimates of the accuracy of approximations obtained are found. It is shown that, in the case of a priori termination, the accuracy of the approximation is proportional to the error in the input data. Under certain additional conditions, the same estimate is established for a posterior termination from the residual principle. These results generalize known similar estimates for linear equations with a normally solvable operator.

Computational Mathematics and Mathematical Physics. 2016;56(9):1523-1535
pages 1523-1535 views

An analysis of the matrix equation AX + X̄B = C

Ikramov K.

Resumo

The matrix equation AX + X̄B = C, where X̄B is obtained by the entry-wise conjugation of X, is examined. On the basis of analogy with the Sylvester matrix equation, special cases are distinguished where the former equation corresponds to normal and self-adjoint Sylvester equations. Efficient numerical algorithms are proposed for these special cases.

Computational Mathematics and Mathematical Physics. 2016;56(9):1536-1539
pages 1536-1539 views

Randomized projection method for estimating angular distributions of polarized radiation based on numerical statistical modeling

Mikhailov G., Tracheva N., Ukhinov S.

Resumo

To study the intensity of radiation transmitted through a layer of substance, a Monte Carlo algorithm is developed based on the expansion of the corresponding angular distribution density in terms of orthonormalized polynomials with a “Lambert” weight. The algorithm is optimized so as to simplify the computations as much as possible. Even a small effect of polarization and the deviation of the angular distribution from the Lambert one can be estimated rather accurately by applying the algorithm.

Computational Mathematics and Mathematical Physics. 2016;56(9):1540-1550
pages 1540-1550 views

Numerical continuation of solution at a singular point of high codimension for systems of nonlinear algebraic or transcendental equations

Krasnikov S., Kuznetsov E.

Resumo

Numerical continuation of solution through certain singular points of the curve of the set of solutions to a system of nonlinear algebraic or transcendental equations with a parameter is considered. Bifurcation points of codimension two and three are investigated. Algorithms and computer programs are developed that implement the procedure of discrete parametric continuation of the solution and find all branches at simple bifurcation points of codimension two and three. Corresponding theorems are proved, and each algorithm is rigorously justified. A novel algorithm for the estimation of errors of tangential vectors at simple bifurcation points of a finite codimension m is proposed. The operation of the computer programs is demonstrated by test examples, which allows one to estimate their efficiency and confirm the theoretical results.

Computational Mathematics and Mathematical Physics. 2016;56(9):1551-1564
pages 1551-1564 views

Differentiation of a functional in the problem of parametric coefficient optimization in the global electric circuit equation

Chernov A.

Resumo

For the problem of parametric optimization of the coefficient and the right-hand side of the linear global electric circuit equation, formulas for the first partial derivatives of an integral cost functional with respect to control parameters are obtained.

Computational Mathematics and Mathematical Physics. 2016;56(9):1565-1579
pages 1565-1579 views

Multiobjective optimization in a pseudometric objective space as applied to a general model of business activities

Khachaturov R.

Resumo

It is shown that finding the equivalence set for solving multiobjective discrete optimization problems is advantageous over finding the set of Pareto optimal decisions. An example of a set of key parameters characterizing the economic efficiency of a commercial firm is proposed, and a mathematical model of its activities is constructed. In contrast to the classical problem of finding the maximum profit for any business, this study deals with a multiobjective optimization problem. A method for solving inverse multiobjective problems in a multidimensional pseudometric space is proposed for finding the best project of firm’s activities. The solution of a particular problem of this type is presented.

Computational Mathematics and Mathematical Physics. 2016;56(9):1580-1590
pages 1580-1590 views

Approximate solution of the p-median minimization problem

Il’ev V., Il’eva S., Navrotskaya A.

Resumo

A version of the facility location problem (the well-known p-median minimization problem) and its generalization—the problem of minimizing a supermodular set function—is studied. These problems are NP-hard, and they are approximately solved by a gradient algorithm that is a discrete analog of the steepest descent algorithm. A priori bounds on the worst-case behavior of the gradient algorithm for the problems under consideration are obtained. As a consequence, a bound on the performance guarantee of the gradient algorithm for the p-median minimization problem in terms of the production and transportation cost matrix is obtained.

Computational Mathematics and Mathematical Physics. 2016;56(9):1591-1597
pages 1591-1597 views

Solvability of certain inverse problems for the nonstationary kinetic transport equation

Volkov N.

Resumo

Linear and nonlinear inverse problems for the nonstationary multispeed anisotropic kinetic transport equation are studied. Sufficient conditions for the existence and uniqueness of weak solutions to these problems in various function spaces are found. The proofs of the corresponding theorems imply that solutions of the inverse problems under study can be obtained by applying the method of successive approximations.

Computational Mathematics and Mathematical Physics. 2016;56(9):1598-1603
pages 1598-1603 views

Pseudo-spectral fourier method as applied to finding localized spherical soliton solutions of (3 + 1)-dimensional Klein–Gordon equations

Ekomasov E., Salimov R.

Resumo

Nonlinear Klein–Gordon equations with fractional power and logarithmic potentials and with a variation in the φ4 potential are found for which the existence of long-lived stable spherically symmetric solutions in the form of pulsons is numerically established. Their mean oscillation amplitude and the frequency of the fast oscillation mode do not vary in the course of the numerical simulation. It is shown that the stability of these pulsons is explained by the presence of a potential well.

Computational Mathematics and Mathematical Physics. 2016;56(9):1604-1610
pages 1604-1610 views

Simulation of collisionless ultrarelativistic electron–proton plasma dynamics in a self-consistent electromagnetic field

Ginzburg S., Dyachenko V., Orlov Y., Fimin N., Chechetkin V.

Resumo

The evolution of a collisionless electron–proton plasma in the self-consistent approximation is investigated. The plasma is assumed to move initially as a whole in a vacuum with the Lorentz factor. The behavior of the dynamical system is analyzed by applying a three-dimensional model based on the Vlasov–Maxwell equations with allowance for retarded potentials. It is shown that the analysis of the solution to the problem is not valid in the “center-of-mass frame” of the plasmoid (since it cannot be correctly defined for a relativistic plasma interacting via an electromagnetic field) and the transition to a laboratory frame of reference is required. In the course of problem solving, a chaotic electromagnetic field is generated by the plasma particles. As a result, the particle distribution functions in the phase space change substantially and differ from their Maxwell–Juttner form. Computations show that the kinetic energies of the electron and proton components and the energy of the self-consistent electromagnetic field become identical. A tendency to the isotropization of the particle momentum distribution in the direction of the initial plasmoid motion is observed.

Computational Mathematics and Mathematical Physics. 2016;56(9):1611-1619
pages 1611-1619 views

Numerical study of wave propagation in porous media with the use of the grid-characteristic method

Kvasov I., Leviant V., Petrov I.

Resumo

Elastic wave propagation in a porous medium is numerically studied by applying the grid-characteristic method. On the basis of direct measurements of reflected and transmitted wave amplitudes, the reflection and decay coefficients are investigated as depending on the degree of porosity (percentage of the pore volume) and on the type of the filling substance (solid, liquid, or nothing). The reflection and decay coefficients are shown to be closely related to the porosity of the medium, which can be used in geological applications (estimation of porosity) and engineering applications (acoustic response attenuation).

Computational Mathematics and Mathematical Physics. 2016;56(9):1620-1630
pages 1620-1630 views

On the equivalence of the electromagnetic problem of diffraction by an inhomogeneous bounded dielectric body to a volume singular integro-differential equation

Smirnov Y.

Resumo

The paper is concerned with the smoothness of the solutions to the volume singular integrodifferential equations for the electric field to which the problem of electromagnetic-wave diffraction by a local inhomogeneous bounded dielectric body is reduced. The basic tool of the study is the method of pseudo-differential operators in Sobolev spaces. The theory of elliptic boundary problems and field-matching problems is also applied. It is proven that, for smooth data of the problem, the solution from the space of square-summable functions is continuous up to the boundaries and smooth inside and outside of the body. The results on the smoothness of the solutions to the volume singular integro-differential equation for the electric field make it possible to resolve the issues on the equivalence of the boundary value problem and the equation.

Computational Mathematics and Mathematical Physics. 2016;56(9):1631-1640
pages 1631-1640 views

Generation of a longitudinal current by a transverse electromagnetic field in collisional degenerate plasma

Latyshev A., Yushkanov A.

Resumo

From the Vlasov–Boltzmann kinetic equation for a collisional degenerate plasma, the electron distribution function is constructed in the quadratic approximation in the electric field strength. A formula for calculating the electric current is derived. It is shown that nonlinearity leads to the rise of a longitudinal electric current directed along the wave vector. The longitudinal current is orthogonal to the known transverse classical current obtained in the linear analysis. When the collision frequency tends to zero, all results obtained for a collisional plasma pass into the corresponding results for a collisionless plasma. The case of small wavenumbers is considered. It is shown that, when the collision frequency tends to zero, the expression for the current passes into the corresponding expression for the current in a collisionless plasma. Graphic analysis of the real and imaginary parts of the current density is performed. The dependence of the electromagnetic field oscillation frequency and electron–plasma-particle collision frequency on the wavenumber is studied.

Computational Mathematics and Mathematical Physics. 2016;56(9):1641-1650
pages 1641-1650 views

Efficient parallel shock-capturing method for aerodynamics simulations on body-unfitted cartesian grids

Menshov I., Pavlukhin P.

Resumo

For problems with complex geometry, a numerical method is proposed for solving the three-dimensional nonstationary Euler equations on Cartesian grids with the use of hybrid computing systems. The baseline numerical scheme, a method for implementing internal boundary conditions on body-unfitted grids, and an iterative matrix-free LU-SGS method for solving the discretized equations are described. An efficient software implementation of the numerical algorithm on a multiprocessor hybrid CPU/GPU computing system is considered. Results of test computations are presented.

Computational Mathematics and Mathematical Physics. 2016;56(9):1651-1664
pages 1651-1664 views

MHD three-dimensional flow of viscoelastic fluid over an exponentially stretching surface with variable thermal conductivity

Alsaedi A., Hayat T., Muhammad T., Shehzad S.

Resumo

This study models the magnetohydrodynamic (MHD) three-dimensional boundary layer flow of viscoelastic fluid. The flow is due to the exponentially stretching surface. The heat transfer analysis is performed through prescribed surface temperature (PST) and prescribed surface heat flux (PHF). The thermal conductivity is taken temperature dependent. Series solutions of velocities and temperatures are constructed. Graphical results for PST and PHF cases are plotted and analyzed. Numerical values of skin-friction coefficients and Nusselt numbers are presented and discussed.

Computational Mathematics and Mathematical Physics. 2016;56(9):1665-1678
pages 1665-1678 views