Vol 47, No 2 (2024)
Mathematics
The control problem for a heat conduction equation with Neumann boundary condition
Abstract
Previously, boundary control problems for a heat conduction equation with Dirichlet boundary condition were studied in a bounded domain. In this paper, we consider the boundary control problem for the heat conduction equation with Neumann boundary condition in a bounded one-dimensional domain. On the part of the border of the considered domain, the value of the solution with control parameter is given. Restrictions on the control are given in such a way that the average value of the solution in some part of the considered domain gets a given value. The studied initial boundary value problem is reduced to the Volterra integral equation of the first type using the method of separation of variables. It is known that the solution of Volterra’s integral equation of the first kind cannot always be shown to exist. In our work, the existence of a solution to the Volterra integral equation of the first kind is shown using the method of Laplace transform. For this, the necessary estimates for the kernel of the integral equation were found. Finally, the admissibility of the control function is proved.



Mathematical modeling
Mathematical Model of Van der Pol-Airy Fractional Oscillator
Abstract
The paper proposes a mathematical model of the nonlinear Van der Pol-Airy oscillator taking into account heredity. The nonlinearity of the oscillator is due to the dependence of the friction coefficient on the square of the displacement function, which is typical for the Van der Pol oscillator. Also, the natural frequency of oscillations is a function of time, which increases linearly as it increases. The latter is typical for the Airy oscillator. Heredity effects are introduced into the model equation through fractional derivatives in the Gerasimov-Caputo sense. They indicate that the oscillatory system may have memory effects that manifest themselves depending on its current state from previous ones. For the proposed mathematical model, a numerical algorithm was developed based on an explicit first-order finite-difference scheme. The numerical algorithm was implemented in a computer program in the Maple language, with the help of which the simulation results were visualized. Oscillograms and phase trajectories were constructed for various values of the model parameters. It is shown that a fractional mathematical model can have various oscillatory modes: from self-oscillatory, damped and chaotic. An interpretation of the simulation results is given.



The optimization problem for determining the functional dependence of the variable order of the fractional derivative of the Gerasimov-Caputo type
Abstract
When solving mathematical modeling problems, it is often necessary to turn to the theory of integral and differential calculus. This theory can be used to describe dynamic processes of various types. The use of fractional derivatives allows us to refine some models by taking into account the memory effect, which is expressed in the equations depending on the current state of the system from previous states. This effect is called non-locality and its intensity is determined by the value of the exponent in the fractional derivative. Classically, this value a a non-integer constant, but there are also generalizations for time-varying nonlocality and other functional dependencies. Fractional differential models are finding increasing application in the physical, mathematical, and technical sciences. However, given the nature of the modeled process, the selection of various parameters for such models must be carried out empirically. Model parameters are refined by iterating through values and comparing simulation results with experimental data representing the process. This process continues until the results begin to qualitatively approximate the data, which is a time-consuming process that inevitably leads to ideas about solving inverse problems. The purpose of this work is to demonstrate that it is possible to use methods of unconditional optimization to solve inverse problems and determine the type of functional dependence a(t). The direct problem is formulated as a Cauchy problem for a fractional differential equation, where the derivative is interpreted in the sense of Gerasimov-Caputo with a variable exponent a(t) for the fractional derivative. The direct problem is solved numerically using a non-local, implicit finite difference scheme. The inverse problem is defined as the problem of discrete minimization of the function a(t) based on experimental data. To solve this problem, we have chosen the Levenberg-Marquardt iterative method. Through test examples, we have shown that this method can be used for unconstrained optimization to determine the shape of the function a(t) and its optimal values in various models.



Information and computing technologies
Search signs changes in the state of the near-surface layer of the lithosphere based on the image analysis results reflecting the characteristics dynamics of the geoacoustic emission signal
Abstract
The lithosphere physical state, including its surface layers, in a certain way characterizes the likelihood of the risk of catastrophic seismic events, indicates the degree of threat to human life and the risk of economic damage. One of the directions for assessing the lithosphere state is the analyzing result the characteristics variability of the geoacoustic emission signal at the observation point. The studies results carried out at the Institute of Cosmophysical Research and Radio Wave Propagation Far Eastern Branch of the Russian Academy of Sciences are presented. This make it possible to identify the dynamic characteristics of the geoacoustic emission signal associated with changes in the near-surface of the lithosphere. Recognition and assessment of the characteristics variability of geoacoustic emission signals on a given time scale of observations was carried out using a neural network approach. A methodology has been developed for classifying observed anomalies in the dynamics of geoacoustic emission signal characteristics.



Processing and preparation of observation data in the interests of highlighting the features of the dynamics of the characteristics of geoacoustic emission
Abstract
The lithospheric layer deformation under the action of seismic processes affects the characteristics of geoacoustic emission. The study of the geoacoustic emission dynamics is aimed at finding signs of pre-seismic events. There is a problem obtained for the high-quality processing of geoacoustic emission signals and the results classification. The study is aimed at finding the best combination of pre-processing and clustering tools for the pulse flow of geoacoustic emission to identify the features of the characteristics dynamics of such a signal. The processed signals were obtained during long-term measurements in the surface lithosphere layers of the seismically active region of the Kamchatka Peninsula. To identify the variability features of geoacoustic emission signals characteristics they are converted by sructurno-linguistic into a three-dimensional image. The images are processed, compared and clustered using convolutional neural networks of various architectures. The best result is assessed by three selected quality criteria. A technique has been developed for finding the best preprocessing and clustering result. The experimental result analisys are presented.



Physics
Generalization of the Results of Lidar Monitoring of the Temperature of the Middle Atmosphere over Tomsk
Abstract
The work carried out processing and analysis of long-term series of vertical temperature distribution in the middle atmosphere (altitude interval 15-60 km) over Tomsk, obtained on the basis of regular measurements of the lidar station of the Institute of Atmospheric Optics SB RAS for the period 2010-2023. During this period, about 1,000 individual night measurements were carried out and the corresponding total backscatter signals at a wavelength of 532 nm were accumulated and processed.. Based on the experimental data obtained, regional features of intra- and interannual variability of the thermal state of the middle atmosphere over Western Siberia were identified. For the warm period of each year (from May to September), a stable temperature distribution has been established with differences of up to several K in individual years. In spring and autumn, the difference between the average monthly temperature profiles increases to 5-10 K and in January reaches a maximum of 15 K. A fundamental difference between the vertical temperature distribution constructed from lidar data and the proposed CIRA-86 model has been revealed. From October to April, in the altitude range from 15 to 25 km, the lidar temperature profile is shifted from the model one in the negative direction, and above from 25 to 50 km, in the positive direction. The maximum negative shift in profiles is observed in December and is 16 K, and negative up to 15 K in January. A description is given of methods for analyzing a number of experimental data for the presence of values that differ significantly from the rest of the data and can lead to distortion of the results obtained.



Lidar studies of the vertical distribution of aerosol in the stratosphere over Тomsk in 2023
Abstract
The article presents studies of the intra-annual variability of the vertical-temporal structure of the background aerosol and its integral content in the stratosphere based on an analysis of the optical characteristics of the aerosol. The results were obtained at the lidar complex of the station for high-altitude atmospheric sounding of the IAO SB RAS for 2023. The sounding was carried out at night, the sounding altitude range was from 10 to 50 km. The optical characteristic — aerosol scattering ratio – is presented as a parameter describing the vertical stratification of the aerosol. According to monitoring data, as in previous years, a stable tendency for aerosol accumulation in the lower stratosphere during the cold period of the year with maximum content in January ( at altitudes of 10–30 km) has been established. In March, a decrease begins, and the vertical stratification of aerosol is estimated at an average value of in the 10-20 km layer, to virtual absence in June-July (). From September, the aerosol content in the lower stratosphere begins to increase to its maximum value in winter. In the upper stratosphere (30-50 km) there is virtually no background aerosol throughout the year. As a comparative analysis and control over the temporal dynamics of the aerosol filling of the stratosphere, the article provides parameter B - the integral aerosol backscatter coefficient for the period from 2017 to the beginning of 2024. In general, over the years, the maximum aerosol filling of the stratosphere occurs in the winter months and is centered relative to January. At the same time, it is unevenly distributed over the years and is marked by sharp and significant fluctuations.



Formation of capillary-gravity waves in the flow under the influence of a system consisting of two vortices
Abstract
Capillary-gravity waves significantly change the general circulation of the water surface of the World Ocean: attenuation and collapse of longer waves, gas exchange and mixing in the upper layer, kinematics of surface suspensions, and, therefore, require additional study. In the framework of a two-dimensional problem, the paper considers surface capillary-gravity waves without taking into account wind effects, but taking into account the isobaric approximation. The approach proposed by Keldysh is chosen as a quantitative basis. The object of perturbation located at a finite depth is a vortex dipole with constant curl. Two asymptotic solutions were obtained on the basis of analytical calculations: the first solution describes the profile of capillary-gravity waves located behind the vortex dipole; the second solution is the profile of capillary-gravity waves in front of the perturbation source. It is shown that the capillary component of the wave dominates in the formation of waves in front of the perturbation source, and the gravitational component dominates behind the obstacle. The relations between isobaric and barotropic effects on the free surface are qualitatively analyzed. The case for the gravitational component is considered; for this purpose, the technique of representing the wave profile using a Taylor polynomial was used.



Instruments and measurement methods
Experience with the «Triton-M» georadar device in the Petropavlovsk-Kamchatsky city (Kamchatka)
Abstract
Determining the possibilities of using the ground penetrating radar method, which is a method of shallow geophysics, in seismically active areas will increase the information content and reliability of the results, especially in engineering and geological surveys to ensure seismic-resistant construction in geodynamically active areas, which include Kamchatka. To optimize the ground penetrating radar method, it is necessary to compare and analyze the results of observations in various geological conditions. Since the ground penetrating radar method has been intensively developing in recent years, theoretical research in the field of data processing must be illustrated by numerous examples of the practical use of ground penetrating radars at various research sites. The purpose of the research is to create a scientific and methodological basis for the georadar method in relation to Kamchatka, including methods for processing and interpreting data, taking into account practical experience of application at various sites. The object of the research is the geological features and soil conditions of the territory of the Petropavlovsk-Kamchatsky city, requiring clarification and detailed study. The article describes the experience and features of working with the «Triton-M» georadar device in the Petropavlovsk-Kamchatsky city. Some results of ground penetrating radar profiling are characterized. Recommended values of measurement parameters during work are given. From a practical point of view, the most important parameters of the ground penetrating radar are the probing depth and the depth resolution. The «Triton-M» ground penetrating radar device is applicable for solving geological problems in the territory of the Petropavlovsk-Kamchatsky city, while optimally selected values of measurement parameters minimize the recording of interference and provide more accurate useful information about the probed environment. For accurate interpretation of radargrams, it is necessary to use information about the geological structure of the studied territory.



Experience with the «PYTHON-3» georadar device in the Petropavlovsk-Kamchatsky city (Kamchatka)
Abstract
Determining the possibilities of using the ground penetrating radar method, which is a method of shallow geophysics, in seismically active areas will increase the information content and reliability of the results, especially in engineering and geological surveys to ensure seismic-resistant construction in geodynamically active areas, which include Kamchatka. To optimize the ground penetrating radar method, it is necessary to compare and analyze the results of observations in various geological conditions. Since the ground penetrating radar method has been intensively developing in recent years, theoretical research in the field of data processing must be illustrated by numerous examples of the practical use of ground penetrating radars at various research sites. The purpose of the research is to create a scientific and methodological basis for the georadar method in relation to Kamchatka, including methods for processing and interpreting data, taking into account practical experience of application at various sites. The object of the research is the geological features and soil conditions of the territory of the Petropavlovsk-Kamchatsky city, requiring clarification and detailed study. The article describes the experience and features of working with the «PYTHON-3» georadar device in the Petropavlovsk-Kamchatsky city. Some results of ground penetrating radar profiling are characterized. Recommended values of measurement parameters during work are given. Georadar «PYTON-3» - geophysical equipment for radar subsurface probing. This georadar is informative for research at great depths when solving geological problems in the Petropavlovsk-Kamchatsky city. It is worth considering that the information content of the upper part of the probed environment is lost (smaller geological structures may not be detected). For an accurate interpretation of the resulting radargrams, it is necessary to use information about the geological structure of the territory under study.


