Статьи журнала - Вестник Южно-Уральского государственного университета. Серия: Математическое моделирование и программирование
Все статьи: 729
Статья научная
In this paper, we prove some new results on operational second order differential equations of elliptic type with general Robin boundary conditions in a non-commutative framework. The study is performed when the second member belongs to a Sobolev space. Existence, uniqueness and optimal regularity of the classical solution are proved using interpolation theory and results on the class of operators with bounded imaginary powers. We also give an example to which our theory applies. This paper improves naturally the ones studied in the commutative case by M. Cheggag, A. Favini, R. Labbas, S. Maingot and A. Medeghri: in fact, introducing some operational commutator, we generalize the representation formula of the solution given in the commutative case and prove that this representation has the desired regularity.
Бесплатно
Non-Uniqueness of Solutions to Boundary Value Problems with Wentzell Condition
Статья научная
Recently, in the mathematical literature, theWentzell boundary condition is considered from two points of view. In the first case, let us call it classical one, this condition is an equation containing a linear combination of the values of the function and its derivatives on the boundary of the domain. Moreover, the function itself also satisfies the equation with an elliptic operator defined in the domain. In the second case, which we call neoclassical one, the Wentzell condition is an equation with the Laplace–Beltrami operator defined on the boundary of the domain understood as a smooth compact Riemannian manifold without boundary, and the external action is represented by the normal derivative of a function defined in the domain. The paper shows the non-uniqueness of solutions to boundary value problems with the Wentzell condition in the neoclassical sense both for the equation with the Laplacian and for the equation with the Bi-Laplacian given in the domain.
Бесплатно
Статья научная
Degenerate differential equations, as part of the differential-algebraic equations, the last few decades cause increasing interest among researchers, both because of the attractiveness of the considered theoretical questions, and by virtue of their applications. Currently, advanced methods developed in this area are used for system modelling and analysis of electrical and electronic circuits, chemical reaction simulations, optimization theory and automatic control, and many other areas. In this paper, the theory of normal forms of differential equations, originated in the works of Poincare and recently developed in the works of Arnold and his school, adapted to the simplest case of a degenerate differential equations. For this purpose we are using technique of Jordan chains, which was widely used in various problems of bifurcation theory. We study the normal forms of degenerate differential equations in the case of the existence of the maximal Jordan chain. Two and three dimensional spaces are studied in detail. Normal forms are the simplest representatives of the degenerate differential equations, which are equivalent to more complex ones. Therefore, normal forms should be considered as a model type of degenerate differential equations.
Бесплатно
Numerical analysis of fractional order integral dynamical models with piecewise continuous kernels
Статья научная
Volterra integral equations find their application in many areas, including mathematical physics, control theory, mechanics, electrical engineering, and in various industries. In particular, dynamic Volterra models with discontinuous kernels are effectively used in power engineering to determine the operating modes of energy storage devices, as well as to solve the problem of load balancing. This article proposes the numerical scheme for solution of the fractional order linear Volterra integral equations of the first kind with piecewise continuous kernels. The developed approach is based on a polynomial collocation method and effectively approximate such a weakly singular integrals. The efficiency of proposed numerical scheme is illustrated by two examples.
Бесплатно
Numerical investigation of the Boussinesq - Love mathematical models on geometrical graphs
Краткое сообщение
The article is devoted to the numerical investigation of the Boussinesq - Love mathematical models on geometrical graphs representing constructions made of thin elastic rods. The first paragraph describes the developed algorithm for numerical solution of the Boussinesq - Love equation with initial conditions and boundary conditions in the vertices. The block diagram of the algorithm is given and described. The result of computation experiment is given in the second paragraph.
Бесплатно
Краткое сообщение
The process of unsteady flow of incompressible viscoelastic fluid in a cylindrical tube of constant cross-section is considered. To describe the rheological properties of a viscoelastic fluid, the Kelvin-Voigt model is used and the mathematical model of this process is presented as an integro-differential partial differential equation. Within the framework of this model, the problem is to determine the pressure drop along the length of the pipe, which ensures the passage of a given flow rate of viscoelastic fluid through the pipe. This problem belongs to the class of inverse problems related to the recovery of the right parts of integro-differential equations. By replacing variables, the integro-differential equation is transformed into a third-order partial differential equation. First, a discrete analog of the problem is constructed using finite-difference approximations. To solve the resulting difference problem, we propose a special representation that allows splitting the problems into two mutually independent second-order difference problems. As a result, an explicit formula is obtained for determining the approximate value of the pressure drop along the length of the pipeline for each discrete value of the time variable. Based on the proposed computational algorithm, numerical experiments were performed for model problems.
Бесплатно
Numerical modelling of the dynamics of the galactic halos in the colliding galaxies
Статья научная
Based on parallel three-dimensional simulation of N-body and gas self-consistent dynamics, we study the behavior of hot coronal gas in the colliding galaxies with "live'' dark matter halos. We model a few scenarios of the galactic collisions including "bull-eye'' and non-central ones, and use different values of the initial velocities of the colliding galaxies. Taking into account the self-gravity, we demonstrate that the collision of gaseous and stellar components does not lead to the formation of a gaseous "protogalaxy'' observed in some numerical simulations. Also, we show that about sixty percent of hot halo gas is expelled into intergalactic space during the collision. Numerical simulations show that a considerable amount of gas (up to 70% for a bull-eye collisions) exchanges between two colliding galaxies.
Бесплатно
Numerical research of the Barenblatt - Zheltov - Kochina stochastic model
Краткое сообщение
At present, investigations of Sobolev-type models are actively developing. In the solution of applied problems the results allowing to get their numerical solutions are very significant. In the article the algorithm for numerical solving of the initial boundary value problem is developed. The problem describes the pressure distribution of the homogeneous fluid in the horizontal layer in the circle. The layer is opened by a vertical well of a small radius. In our research we suppose that random disturbing loads have an influence on the fluid. The problem was solved under two assumptions. Firstly, we suppose that an unstable fluid flow is axially symmetric, and secondly, that in initial moment the pressure in the layer is constant. After the process of the discretization we modify the original model to the Cauchy problem for the system of ordinary differential equations. For the numerical solution we use algorithms based on explicit one-step formulas of the Runge - Kutta type with the seventh-order accuracy and with the selection of the integration step. We also use the scheme of the eighth-order accuracy to evaluate the calculation accuracy on each steps of time. According to the results of this control, we choose the time-step. A lot of numerical experiments have shown high numerical efficiency of the algorithm that we use to solve the investigated initial-boundary problem.
Бесплатно
Numerical research of the mathematical model for traffic flow
Краткое сообщение
The problems of distribution of transport flows are currently relevant in connection with the increase in vehicles. In the 50s of the last century, the first macroscopic (hydrodynamic) models appeared, where the transport flow resembles the flow "motivated" compressible liquid. The scientific approach based on the Navier - Stokes system. The main idea of the scholars is consideration the hydrodynamic models on the grounds of interrelation between the transport flow and incompressible fluid. For modelling traffic flows we examine Oskolkov equation on the geometric graph, where the edge has two positive values corresponding to it "length" and "width". Certainly, in the context of mathematical model the values lk and bk are dimensionless, but for clarity it is convenient to imagine that lk is measured in linear metric units, for example, kilometers or miles, and bk is equal to the number of traffic lanes on the roadway in one direction. In terms of the Oskolkov model, we obtained a non-classical multipoint initial-final value condition. We will study such a model using the idea and methods of the Sobolean equation theory. These notes describe a numerical experiment based on the Galerkin method for the Oskolkov equation with a multipoint initial-final condition on the graph.
Бесплатно
Статья научная
The existing mathematical models of pressure swing adsorption (PSA) apply various assumptions regarding the mass and heat transfer mechanisms in the "gas mixture-adsorbent'' system. An increase in the number of assumptions leads to a simplification of the model, a decrease in the calculation time of one iteration in the model and, at the same time, a decrease in its accuracy. The simplification of the model is especially important in PSA processes, since the calculation of the model is carried out before the cyclic steady state and takes tens and even hundreds of cycles (iterations). Ensuring high accuracy of the PSA model and its minimum complexity is a contradictory requirement; therefore it is important to reasonably consider only those transfer mechanisms that are dominant in the model. The paper proposes a mathematical model of the PSA process, which takes into account the thermal effects of sorption, external and internal diffusion mechanisms of adsorptive transfer. A numerical research was carried out to determine the dominant transfer mechanism, and recommendations were proposed for using the preferred PSA model in terms of its accuracy and calculation time (for the processes of air oxygen enrichment and synthesis gas separation). It was found that to calculate PSA oxygen units with a capacity of less than 4 l/min at NTP, it is advisable to use an isothermal model, which saves at least 24,3% of the calculation time with a loss of accuracy of no more than 0,084 vol%. To calculate PSA hydrogen units, the use of an isothermal model is impractical even at the lowest productivity of 50 l/min at NTP. When the diameter of the adsorbent particles is less than 2 mm, it is advisable to use an external diffusion model, which saves at least 54,2% of the calculation time for oxygen units and at least 47,1% of the calculation time for hydrogen units with a slight loss of accuracy. At a gas flow velocity of more than 0,05 m/s, the model can ignore the diffusion in the gas. The research results can be used to calculate various PSA processes for separation of gas mixtures: rPSA, ultra rPSA, VSA, VPSA, and related processes.
Бесплатно
Numerical study of the SUSUPLUME air pollution model
Статья научная
In this paper, we propose a SUSUPLUME air pollution as a modern application of the classical Gaussian plume model. The presented model takes into account meteorological conditions and parameters of the pollution sources. The classical model is supplemented by the equations of motion of the center of mass of a single emission. A numerical study has shown that in stationary weather conditions the presented model qualitatively coincides with other known models. The results of calculating the concentrations of pollutants do not contradict the obtained values based on the official methodology for calculating the maximum concentrations of pollutants approved for usege in the territory of the Russian Federation. The SUSUPLUME model contains a number of identifiable parameters and it can be adapted to real conditions. The computational model consists of two blocks: a block for recording measurement information and a block for calculating the concentrations of pollutants. The measurement information registration unit has a low labor intensity (over a million registrations per second). The pollutant concentrations calculation block is laborious (400 points of calculations per second). Concentrations are calculated independently, it allows to use parallelization of the computational process in the future.
Бесплатно
Numerical study of the dynamics of air separation process by pressure swing adsorption
Статья научная
Using mathematical modelling and the finite element method, we carry out the calculation experiments to study the system connections and regularities of pressure swing adsorption process under the conditions of air separation and oxygen concentration (production). We study the influence of mode and construction variables on the dynamics and technological indicators of the effectiveness of this process. Namely, we study the influence of input variables (composition and temperature of atmospheric air, air pressure at the compressor outlet) on output variables (extraction degree, oxygen purity, unit capacity, etc.) of the studied object in a wide range of control variables (adsorption-desorption cycle time, pressure ratios at adsorption and desorption stages, and oxygen-enriched reverse air flow coefficient). Also, we study the influence of construction parameters (layer height, particle diameter and maximum adsorption volume of the adsorbent) on the amount of adsorption, which is equilibrium with the current concentration of the adsorptive in the gas mixture flow on the outer surface of the adsorbent granules, the value of the kinetic adsorption coefficient (the coefficient of external mass transfer of the adsorptive (mainly nitrogen) from the gas phase into the adsorbent). The results of calculation experiments allow to establish the most promising mode and construction parameters for the optimal design of oxygen enrichment systems by pressure swing adsorption with varying pressure.
Бесплатно
Observability of square membranes by Fourier series methods
Статья научная
Fourier series methods have been successfully applied in control theory for a long time. Some theorems, however, resisted this approach. Some years ago, Mehrenberger succeeded in establishing the boundary observability of vibrating rectangular membranes (and of analogous higher dimensional problems) by developing an ingenious generalization of Ingham's classical theorem on nonharmonic Fourier series. His method turn out to be useful for other applications as well. We improve Mehrenberger's approach by a shorter proof, and we improve and generalize some earlier applications.
Бесплатно
On a class of Sobolev-type equations
Статья научная
The article surveys the works of T.G. Sukacheva and her students studying the models of incompressible viscoelastic Kelvin-Voigt fluids in the framework of the theory of semilinear Sobolev-type equations. We focus on the unstable case because of greater generality. The idea is illustrated by an example: the non-stationary thermoconvection problem for the order 0 Oskolkov model. Firstly, we study the abstract Cauchy problem for a semilinear nonautonomous Sobolev-type equation. Then, we treat the corresponding initial-boundary value problem as its concrete realization. We prove the existence and uniqueness of a solution to the stated problem. The solution itself is a quasi-stationary semi-trajectory. We describe the extended phase space of the problem. Other problems of hydrodynamics can also be investigated in this way: for instance, the linearized Oskolkov model, Taylor''s problem, as well as some models describing the motion of an incompressible viscoelastic Kelvin-Voigt fluid in the magnetic field of the Earth.
Бесплатно
On a heat and mass transfer model for the locally inhomogeneous initial data
Краткое сообщение
We consider a model case of the problem of heat diffusion in a homogeneous body with a special initial state. The peculiarity of this initial state is its local inhomogeneity. That is, there is a closed domain Ω inside a body, the initial state is constant out of the domain. Mathematical modelling leads to the problem for a homogeneous multi-dimensional diffusion equation. We construct the boundary conditions on the boundary of the domain Ω, which can be characterized as "transparent" boundary conditions. We separately consider a special case - a model of redistribution of heat in a uniform linear rod, the side surface of which is insulated in the absence of (internal and external) sources of heat and of locally inhomogeneous initial state.
Бесплатно
Статья научная
We analyze, by means of singular potentials defined in terms of Dirac functions and their derivatives, a one dimensional symmetry breaking in quantum mechanics. From a mathematical point of view we use a technique of selfadjoint extensions applied to a symmetric differential operator with a domain containing smooth functions which vanish at two inner points of the real line. As is well known, the latter leads to a two-point boundary problem. We compute the resolvent of the corresponding extension and investigate its behavior in the case in which the inner points change their positions. The domain of these extensions can contain some functions with non differentiability or discontinuity at the points mentioned before. This fact can be interpreted as a presence of singular potentials like shifted Dirac delta functions and/or their first derivative centered at the same points. Then, we study the existence of broken-symmetry bound states. For some given entanglement boundary conditions we can show the existence of a ground state, which leads to a spontaneous symmetry breaking. We also prove that within a frame of Pontryagin spaces this type of symmetry breaking is saved if the distance between the mentioned above interior points tends to zero and then we can reformulate this result in terms of a larger Hilbert space.
Бесплатно
On a model of oscillations of a thin flat plate with a variety of mounts on opposite sides
Краткое сообщение
We consider a model case of stationary vibrations of a thin flat plate, one side of which is embedded, the opposite side is free, and the sides are freely leaned. In mathematical modeling there is a local boundary value problem for the biharmonic equation in a rectangular domain. Boundary conditions are given on all boundary of the domain. We show that the considered problem is self-adjoint. Herewith the problem is ill-posed. We show that the stability of solution to the problem is disturbed. Necessary and sufficient conditions of existence of the problem solution are found. Spaces of the ill-posedness of the considered problem are constructed.
Бесплатно
