Vol 28, No 1 (2024)

Differential Equations and Mathematical Physics

Dynamics of a thermal entanglement in the not-resonant three-qubit Tavis-Cummings model with Kerr nonlinearity

Bagrov A.R., Bashkirov E.K.

Abstract

In this article, we consider the dynamics of three identical qubits interacting not-resonantly with a thermal field of an ideal resonator with a Kerr medium. We have found the solutions of the Liouville quantum equation for the total density matrix of a system under consideration for the initial separable, biseparable, and genuine entangled states of the qubits and the thermal initial state of the resonator field. By averaging the total density matrix over the variables of the resonator field and the variables of one of the qubits, we found the reduced density matrix of the pair of remaining qubits. Two-qubit density matrices were used to calculate the qubit-qubit negativity. The results showed that detuning and Kerr nonlinearity can greatly enhance the amout of entanglement for initial separable state of a pair of qubits. It is also shown that detuning and a Kerr medium can inhibit the sudden death of entanglement.

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences. 2024;28(1):7-28
pages 7-28 views

Inverse problem for an integro-differential equation of hyperbolic type with additional information of a special form in a bounded domain

Safarov J.S.

Abstract

A one-dimensional inverse problem of determining the kernel of the integral
term of an integro-differential equation of hyperbolic type in a variablebounded
domain $x$ is considered. Firstly, the direct problem is investigated, for the regular part of which the Cauchy problem on the axis $x=0$ is obtained using the method of singularity extraction. Subsequently, an integral equation for the unknown function is derived by the d’Alembert formula.
For the direct problem, the inverse problem of determining the kernel entering the integral term of the equation is studied. To find it, an additional condition is specified in a special form. As a result, the inverse problem is reduced to an equivalent system of integral equations for unknown functions. The principle of contraction mappings in the space of continuous functions with weighted norms is applied to the obtained system.
For the given problem, a theorem of global unique solvability has been
proven, which is the main result of the study.

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences. 2024;28(1):29-44
pages 29-44 views

A modified Cauchy problem for an inhomogeneous equation of degenerate hyperbolic type of the second kind

Urinov A.K., Okboev A.B.

Abstract

In this study, a modified Cauchy problem was examined for an inhomogeneous equation of degenerate hyperbolic type of the second kind in a characteristic triangle. It is known that degenerate hyperbolic equations have a singularity, meaning that the well-posedness of the Cauchy problem with initial data on the line of parabolic degeneracy does not always hold for them. Therefore, in such cases, it is necessary to consider the problem with initial conditions in a modified form.
In present paper, modified Cauchy problems with initial conditions were formulated on the line of parabolic degeneracy for an inhomogeneous equation of degenerate hyperbolic type of the second kind. The considered problem is reduced to a modified Cauchy problem for a homogeneous equation and to a Cauchy problem for an inhomogeneous equation with zero initial conditions. The solutions of the modified Cauchy problem for a homogeneous equation are derived from the general solution of the considered equation. The explicit solutions of the modified Cauchy problem with homogeneous conditions for the inhomogeneous equation are found using the Riemann
method.
It is proven that the discovered solutions indeed satisfy the equation and the initial conditions.

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences. 2024;28(1):45-58
pages 45-58 views

Mechanics of Solids

Construction of elastic fields in the problem from the action of body forces of a cyclic nature

Ivanychev D.A., Levina E.Y.

Abstract

The paper presents a method for determining the stress-strain state of transversely isotropic bodies of revolution under the action of non-axisymmetric stationary volumetric forces. This problem involves the use of boundary state method definitions. The basis of the space of internal states is formed using fundamental polynomials. The polynomial is placed in any position of the displacement vector of the plane auxiliary state, and the spatial state is determined by the transition formulaes. The set of such states forms a finite-dimensional basis according to which, after orthogonalization, the desired state is expanded into Fourier series with the same coefficients. Series coefficients are scalar products of vectors of given and basic volumetric forces. Finally, the search for an elastic state is reduced to solving quadratures.
The solutions of problems of the theory of elasticity for a transversely isotropic circular cylinder from the action of volumetric forces given by various cyclic laws (sine and cosine) are analyzed. Recommendations are given for constructing the basis of internal states depending on the form of the function of given volumetric forces. The analysis of the series convergence and the estimation of the solution accuracy in graphical form are given.

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences. 2024;28(1):59-72
pages 59-72 views

Numerical method for structural and parametric identification of a mathematical model of incomplete inverse deformation of creep strain

Radchenko V.P., Zoteev V.E., Afanasieva E.A.

Abstract

A new numerical method for parametric and structural identification of the physically nonlinear theory of reversibility of creep strain, valid within the first and second stages, has been developed. A series of stationary creep curves is used as basic experimental information. The problem is reduced to nonlinear regression analysis of determining estimates of random parameters based on time series of a sequence of observations of creep deformation at various constant stresses using difference equations. The obtained relationships between the coefficients of the difference equation and the parameters of nonlinear regression allow us to reduce the problem of estimating the coefficients of a linear-parametric discrete model. Corresponding iterative algorithms for refining parameter estimates with any given accuracy have been developed. Parametric and structural identification of the theory of incomplete reversibility of creep deformation has been carried out for steel EI736 (500 C) and alloys EI437A (700C), VZh98 (900 C), EP693 (700 C). Numerical values of model parameter estimates for these alloys are given. The adequacy of the constructed mathematical models was checked, and the relation between the calculated and experimental data was observed. Experimental data for all materials considered belong to the corresponding calculated confidence intervals for creep deformation, which indicates the reliability of the obtained estimates of the model parameters.

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences. 2024;28(1):73-95
pages 73-95 views

Mathematical models of nonlinear dynamics of functionally graded nano/micro/macroscale porous closed cylindrical Kirchhoff-Love shells

Yakovleva T.V., Krysko V.A.

Abstract

The article presents new mathematical models for the dynamics of nonlinear nano/micro/macro-scale functionally graded porous closed cylindrical shells. The Kirchhoff–Love hypothesis is chosen as the kinematic model for the shells. Geometric nonlinearity is considered according to the von Karman model. Nanoeffects are accounted for using by a modified moment theory of elasticity. Variational and differential equations, as well as boundary and initial conditions, are derived from Hamilton’s principle. A proof of the existence of a solution is conducted based on the theory of generalized solutions to differential equations (using methods of Hilbert spaces and variational methods).
As examples, nano/micro/macro-scale closed cylindrical shells are considered as systems with "almost" an infinite number of degrees of freedom subjected to banded transverse alternating loading. The Bubnov–Galerkin method in higher approximations is adopted as the method for reducing partial differential equations to the Cauchy problem. Its convergence is investigated.
The Cauchy problem is solved using Runge–Kutta methods of fourth to eighth order accuracy and the Newmark method. The application of several numerical methods at each stage of modeling is necessary to ensure the reliability of the obtained results. The study of complex oscillation characteristics of the closed cylindrical nano/micro/macro-scale shell is conducted using nonlinear dynamics methods, which involve constructing signals, phase portraits, applying Fourier analysis, and various wavelet transformations,
among which the Morlet wavelet proved to be the most informative.
An analysis of the type of chaotic oscillations is carried out based on the spectrum of Lyapunov exponents using the Sano–Sawada method and the dominant exponent through several methods: Kanca, Rosenstein, and Wolf. It is shown that the size-dependent parameter and the consideration of porosity have a significant impact on the nature of the oscillations of cylindrical shells. The phenomenon of hyper-chaos has been discovered.

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences. 2024;28(1):96-116
pages 96-116 views

Mathematical Modeling, Numerical Methods and Software Complexes

Implicit iterative scheme based on the pseudo--inversion algorithm and its application

Zhdanov A.I., Sidorov Y.V.

Abstract

A new version of the implicit iterative scheme is proposed for the implementation of which only matrix-vector computational procedures are required. This makes the proposed computational scheme potentially highly efficient for solving a wide class of high-dimensional problems on modern high-performance computing platforms, such as Nvidia Cuda. It is shown that the proposed algorithms can be used to solve ill-conditioned linear systems and least squares problems, as well as to construct iterative regularization algorithms. The results of computational experiments are presented, confirming the effectiveness of the proposed computational algorithms.

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences. 2024;28(1):117-129
pages 117-129 views

Hybrid numerical-analytical method for solving problems of salt ion transport in membrane systems with axial symmetry

Kazakovtseva E.V., Kovalenko A.V., Pismenskiy A.V., Urtenov M.K.

Abstract

The aim of this study is to develop a new hybrid numerical-analytical method for solving boundary value problems with axial symmetry, for example, with a rotating membrane disk, based on matching the asymptotic solution near the cation-exchange membrane (CEM) with the numerical solution in the rest of the region. For this, the following method is used:
1) a basic mathematical model for the transfer of salt ions in an electrochemical cell with a rotating cation-exchange membrane disk is proposed based on the general conservation laws represented by the Nernst-Planck-Poisson and Navier-Stokes equations with natural boundary and initial conditions. This model contains no fitting parameters or simplifying assumptions. However, the numerical solution of the corresponding boundary value problem presents significant computational difficulties for real solution concentrations and large jumps in the potential and angular velocity of the membrane disk rotation, associated with large concentration and potential gradients near the CEM in the quasi-equilibrium space charge region (SCR);
2) the solution region is divided into two parts, one of which is a small cation increase region (CIR) located near the CEM, and the remaining main part of the region (MPOR);
3) in the CIR, an analytical solution is found by the method of matching asymptotic solutions;
4) a simplified mathematical model is constructed in the MPOR, which differs from the basic mathematical model in such a boundary condition at the boundary with the CIR, which then allows us the solution of the corresponding boundary value problem to be matched with the solution in the CIR.
The main result is a hybrid numerical-analytical method that allows one to carry out a numerical analysis of the transfer of salt ions at real concentrations of a binary salt electrolyte solution in a wide range of changes in the potential jump and the angular velocity of the membrane disk.
Based on the results of the work, the following conclusion can be drawn, that the combination of the analytical (asymptotic) method of solving in the region of the boundary layer and the numerical solution in the rest of the region, with the exception of the boundary layer, with their subsequent splicing, makes it possible to construct an effective hybrid numerical-analytical method for solving the problems of salt ion transport in membrane systems with axial symmetry.

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences. 2024;28(1):130-151
pages 130-151 views

EXTENDED MATHEMATICAL MODEL OF THE INVERSE PROBLEM OF NUCLEAR GAMMA-RESONANCE: RELIABILITY AND INFORMATIVE OF APPLICATION

Nemtsova O.M., Konygin G.N., Veselkov I.S.

Abstract

Research the properties of iron-based solid solutions by Mössbauer spectroscopy has the problem of interpreting the results of processing experimental data within the traditional mathematical model. Since the disordered solid solutions, for example, as a result of mechanical activation, are consisted of a set of the different local atomic configurations, the corresponding Mössbauer spectra contain a large number of the shifted relative to each other spectral components with close values of the hyperfine interaction parameters. The magnitude and sign of these shift are determined by many factors: the quantitative distribution of atoms of each type in the coordination spheres, the symmetry of their distribution relative to the quantization axis, the possible local shift relative to the average statistical position in the crystallographic structure, etc. In the mathematical model, as a rule, it’s not possible to taken into account all these effects of the shift by analytically.
The proposed extended mathematical model for describing the Mössbauer spectra of solid solutions makes it possible to take into account the shifts in the spectral components by using Gaussian normal distribution function, as a function of statistical set of local distortions. The width of the Gaussian distribution makes it possible to estimate the degree of local distortions of the crystal lattice that arise due to differences in the sizes of atoms of the mixed components, local distortions of the structure and
symmetry of the environment of the resonant atom.
The inverse problem of nuclear gamma-resonance is formulated by the Fredholm integral equation of the first kind and is an ill-posed problem with a priori constraints on the desired solution. The introduction of two Gaussian functions with a priori unknown linewidths into the kernel of the integral equation leads to the problem of solving the equation by classical methods. Algorithm for obtaining a reliable solution based on the Tikhonov regularization method with correction of the parameters of the kernel of the
integral equation is proposed in this paper. On the examples of the study of real objects, the reliability and informative application of the extended mathematical model of the inverse problem of nuclear gamma-resonance is proved.

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences. 2024;28(1):152-170
pages 152-170 views

Short Communications

The solution to a boundary value problem for a third-order equation with variable coefficients

Apakov Y.P., Umarov R.A.

Abstract

In the paper we consider the second boundary value problem in a rectangular domain for an inhomogeneous partial differential equation of the third
order with multiple characteristics with variable coefficients. The uniqueness
of the solution of the posed problem is proved by the method of energy integrals. The uniqueness theorem is proved. A counter-example is constructed
in the case of violation of the conditions of the uniqueness theorem. By the
method of separation of variables, the solution of the problem is sought as
a product of two functions
X (x) and Y (y). In order to determine Y (y),
we generate an ordinary differential equation of the second order with two
boundary conditions at the boundaries of the segment
[0, q]. For this problem, the eigenvalues and the corresponding eigenfunctions, when n = 0 and
n N, are found. For determining X (x), we generate an ordinary differential
equation of the third order with three boundary conditions at the boundaries
of the segment
[0, p]. A new function, which makes the boundary conditions
homogeneous, is introduced. The solution of the given problem is constructed
by means of the Green function. For
n = 0 and for n N, the Green function
is individually constructed. It is verified that the obtained Green’s functions
satisfy the boundary conditions and properties of the Green function. The
solution
X (x) is written by virtue of the constructed Green functions. After
some transformations, the Fredholm integral equation of the second kind is
generated and its solution is written in terms of the resolvent. Estimates
for the resolvent and the Green function are obtained. The uniform convergence of both solution and its possible partial derivatives are shown under the conditions on the given functions. The convergence of the third order derivative of the solution with respect to the variable x is proved using the Cauchy-Bunyakovsky and Bessel inequalities. When justifying the uniform
convergence of the solution, the absence of a “small denominator” is proved.

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences. 2024;28(1):171-185
pages 171-185 views

Analysis for two-component composite coatings in the production of electronic components by computer vision methods

Nenashev M.V., Rachmanin O.S., Kiyashchenko V.V.

Abstract

A method for analyzing two-component composite coatings in the production of electronic components is presented, based on a system of mathematically grounded image processing algorithms. This method allows for the determination of the specific surface area, total material area, and contact boundaries, ensuring high accuracy and result stability. The obtained results can be successfully integrated into industrial processes for material quality assessment and production control. Within the scope of the study, an information-measurement image processing system has been developed, minimizing error accumulation at each stage and ensuring high precision in determining material characteristics. Examples of successful method application are presented, highlighting its effectiveness and prospects in various areas, including industrial production of electronic components. The obtained results serve as a basis for further research and refinement of methods for analyzing composite materials.

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences. 2024;28(1):186-198
pages 186-198 views