


Volume 56, Nº 9 (2016)
- Ano: 2016
- Artigos: 15
- URL: https://journal-vniispk.ru/0965-5425/issue/view/11098
Article
Iteratively regularized methods for irregular nonlinear operator equations with a normally solvable derivative at the solution
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.



An analysis of the matrix equation AX + X̄B = C
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.



Randomized projection method for estimating angular distributions of polarized radiation based on numerical statistical modeling
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.



Numerical continuation of solution at a singular point of high codimension for systems of nonlinear algebraic or transcendental equations
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.



Differentiation of a functional in the problem of parametric coefficient optimization in the global electric circuit equation
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.



Multiobjective optimization in a pseudometric objective space as applied to a general model of business activities
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.



Approximate solution of the p-median minimization problem
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.



Solvability of certain inverse problems for the nonstationary kinetic transport equation
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.



Pseudo-spectral fourier method as applied to finding localized spherical soliton solutions of (3 + 1)-dimensional Klein–Gordon equations
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.



Simulation of collisionless ultrarelativistic electron–proton plasma dynamics in a self-consistent electromagnetic field
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.



Numerical study of wave propagation in porous media with the use of the grid-characteristic method
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).



On the equivalence of the electromagnetic problem of diffraction by an inhomogeneous bounded dielectric body to a volume singular integro-differential equation
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.



Generation of a longitudinal current by a transverse electromagnetic field in collisional degenerate plasma
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.



Efficient parallel shock-capturing method for aerodynamics simulations on body-unfitted cartesian grids
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.



MHD three-dimensional flow of viscoelastic fluid over an exponentially stretching surface with variable thermal conductivity
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.


