Publication activity of the Teaching staff of the Faculty of Mechanics and Mathematics for the last 5 years
№ 
Title of publication 
Output 
Article abstract 
Citation link (full name, article title, title, number and/or issue, journal volume, pages,

1. 
About the generalized DirichletNeumann problem for an elliptic equation of high order 
DOI:10.1063/1.5049007 
For the elliptic equation 2l th order with constant (and only) real coefficients considered boundary value problem of the job normal derivatives the (kj  1) order, j = 1, ⋯, l, where 1 ≤ k1< ⋯ 
Koshanov B., Soldatov A. . About the generalized DirichletNeumann problem for an elliptic equation of high order // AIP Conference Proceedings.  2018.  Vol.1997. DOI:10.1063/1.5049007

2 
New boundary value problems for fourthorder quasihyperbolic equations 
DOI:10.33048/SEMI.2019.16.098 
In this paper, we study the correctness in the spaces of S.L. Sobolev of new boundary value problems for quasihyperbolic differential equations u_{tttt} + Au = f(x; t) (A is an elliptic operator acting on spatial variables). For the proposed tasks theorems on the existence and uniqueness of solutions are proved, and examples of nonuniqueness are given. 
Kozhanov A.I., Koshanov B., Sultangazieva J. New boundary value problems for fourthorder quasihyperbolic equations // Siberian Electronic Mathematical Reports. – 2019.  Vol.16. – Р. 1410 – 1436. DOI:10.33048/SEMI.2019.16.098 
3 
On Fredholm solvability and on the index of the generalized Neumann problem for an elliptic equation 
DOI:10.1080/17476933.2021.1958797 
In this paper, we investigate the Fredholm solvability of the generalized Neumann problem for a highorder elliptic equation in the plane. The equivalence of the solvability conditions of the generalized Neumann problem to the complementarity condition (Shapiro–Lopatinsky condition) is proved. The formula for the index of the specified problem in the class (Formula presented.) is calculated. © 2021 Informa UK Limited, trading as Taylor & Francis Group. 
Koshanov B., Soldatov A. . On Fredholm solvability and on the index of the generalized Neumann problem for an elliptic equation // Complex Variables and Elliptic Equations. – 2021.  Р. 117. 
4 
A numerical simulation of air flow in the human respiratory system for various environmental conditions 
DOI: 10.1186/s12976020001338 
The functions of the nasal cavity are very important for maintaining the internal environment of the lungs since the inner walls of the nasal cavity control the temperature and saturation of the inhaled air with water vapor until the nasopharynx is reached. In this paper, threedimensional computational studies of airflow transport in the models of the nasal cavity were carried out for the usual inspiratory velocity in various environmental conditions. Threedimensional numerical results are compared with experimental data and calculations of other authors. Numerical results show that during normal breathing, the human nose copes with heat and relative moisture metabolism in order to balance the intraalveolar conditions. It is also shown in this paper that a normal nose can maintain balance even in extreme conditions, for example, in cold and hot weather. The nasal cavity accelerates heat transfer by narrowing the air passages and swirls from the nasal concha walls of the inner cavity. © 2020, The Author(s). 
Issakhov A., Zhandaulet Y., Issakhov A., Abylkassymova A. . A numerical simulation of air flow in the human respiratory system for various environmental conditions // Theoretical Biology and Medical Modelling.  2021.  Vol.18. DOI: 10.1186/s12976020001338 
4 
Solvability issues of a pseudo‐parabolic fractional order equation with a nonlinear boundary condition 
DOI: 10.3390/fractalfract5040134 
This paper is devoted to the fundamental problem of investigating the solvability of initial‐boundary value problems for a quasi‐linear pseudo‐parabolic equation of fractional order with a sufficiently smooth boundary. The difference between the studied problems is that the boundary conditions are set in the form of a nonlinear boundary condition with a fractional differentiation operator. The main result of this work is establishing the local or global solvability of stated problems, depending on the parameters of the equation. The Galerkin method is used to prove the existence of a quasi‐linear pseudo‐parabolic equation’s weak solution in a bounded domain. Using Sobolev embedding theorems, a priori estimates of the solution are obtained. A priori estimates and the Rellich–Kondrashov theorem are used to prove the existence of the desired solutions to the considered boundary value problems. The uniqueness of the weak generalized solutions of the initial boundary value problems is proved on the basis of the obtained a priori estimates and the application of the generalized Gronwall lemma. The need to consider and study such initial boundary value problems for a quasi‐linear pseudo‐parabolic equation follows from practical requirements, such as solving fractional differential equations that simulate physical processes that occur during the study of liquid filtration processes, etc. © 2021 by the authors. Licensee MDPI, Basel, Switzerland. 
Aitzhanov S.E., Berdyshev A.S., Bekenayeva K.S. Solvability issues of a pseudo‐parabolic fractional order equation with a nonlinear boundary condition // Fractal and Fractional.  2021.  Vol.5. DOI: 10.3390/fractalfract5040134

5 
Van der Corput lemmas for MittagLeffler functions. II. α–directions 
DOI 10.1016/j.bulsci.2021.103016 
The paper is devoted to study analogues of the van der Corput lemmas involving MittagLeffler functions. The generalisation is that we replace the exponential function with the MittagLefflertype function, to study oscillatory integrals appearing in the analysis of timefractional partial differential equations. More specifically, we study integral of the form Iα,β(λ)=∫REα,β(iαλϕ(x))ψ(x)dx, for the range 0[removed]0. This extends the variety of estimates obtained in the first part, where integrals with functions Eα,β(iλϕ(x)) have been studied. Several generalisations of the van der Corput lemmas are proved. As an application of the above results, the generalised RiemannLebesgue lemma, the Cauchy problem for the timefractional KleinGordon and timefractional Schrödinger equations are considered. © 2021 The Authors 
Ruzhansky M.,Torebek B.T. . Van der Corput lemmas for MittagLeffler functions. II. α–directions // Bulletin des Sciences Mathematiques. – 2021.  Vol.171. 
6 
On a model of the generation of turbulence 
DOI: 10.1016/j.chaos.2021.111099 
In this article, the nonlinear term of the NavierStokes equation is approximated to nonlinear convolutional expressions. At low values of viscosity, their values are close if the carrier of the convolution is of the same order of magnitude as the value of viscosity. It is expected that the dynamics of the thus obtained modified NavierStokes equation preserves the physical phenomena described by the NavierStokes equation. The dynamics of the modified NavierStokes equation is investigated in this work. © 2021 Elsevier Ltd 
Kanguzhin B.E. . On a model of the generation of turbulence // Chaos, Solitons and Fractals. – 2021.  Vol.150. DOI: 10.1016/j.chaos.2021.111099 
7 
An inverse problem for generalized KelvinVoigt equation with pLaplacian and damping term 
DOI 10.1088/13616420/ac1362 
In this paper, we consider the nonlinear inverse problem for generalized KelvinVoigt equations with the pLaplace diffusion and damping term, describing the motion of incompressible viscous fluids. We assume that the damping term in the momentum equation depends on whether its signal is positive or negative, which may realizes the presence of a source or a sink within the system. The investigated inverse problem consists of finding a coefficient f(t) of the righthand side of the momentum equation, a vector of velocity field v, and a pressure π. An additional information on a solution of the inverse problem is given as integral overdetermination condition. Under several assumptions on the exponents p, m, the coefficients μ, κ, γ, the dimension of the space d, and specified initial data, we prove the existence and uniqueness of the weak solution of the problem. © 2021 IOP Publishing Ltd. 
Antontsev S.N., Khompysh K. . An inverse problem for generalized KelvinVoigt equation with pLaplacian and damping term // Inverse Problems. – 2021.  Vol.37. DOI 10.1088/13616420/ac1362 
8 
Identification of the domain of the sturm–liouville operator on a star graph 
DOI 10.3390/sym13071210 
This article is devoted to the unique recovering of the domain of the Sturm–Liouville operator on a star graph. The domain of the Sturm–Liouville operator is uniquely identified from the set of spectra of a finite number of specially selected canonical problems. In the general case, the domain of the definition of the original operator can be specified by integrodifferential linear forms. In the case when the domain of the Sturm–Liouville operator on a star graph corresponds to the boundary value problem, it is sufficient to choose only finite parts of the spectra of canonical problems for a unique identification of the boundary form. Moreover, the above statement is valid only for a symmetric star graph. © 2021 by the authors. 
Kanguzhin B., Aimal Rasa G.H., Kaiyrbek Z. . Identification of the domain of the sturm–liouville operator on a star graph // Symmetry. – 2021.  Vol.13. DOI 10.3390/sym13071210 
9 
On the localization of the spectrum of some perturbations of a twodimensional harmonic oscillator 
DOI:10.1080/17476933.2021.1885386 
In this paper, we study the localization of the discrete spectrum of certain perturbations of a twodimensional harmonic oscillator. The convergence of the expansion of the source function in terms of the eigenfunctions of a twodimensional harmonic oscillator is investigated. A representation of Green's function of a twodimensional harmonic oscillator is obtained. The singularities of Green's function are highlighted. The wellposed definition of the maximal operator generated by a twodimensional harmonic oscillator on a specially extended domain of definition is given. Then, we describe everywhere solvable invertible restrictions of the maximal operator. We establish that the eigenvalues of a harmonic oscillator will also be the eigenvalues of wellposed restrictions. The results are supported by illustrative examples. © 2021 Informa UK Limited, trading as Taylor & Francis Group. 
Kanguzhin B.E., Fazullin Z. . On the localization of the spectrum of some perturbations of a twodimensional harmonic oscillator // Journal of Mathematical Sciences (United States). – 2021. .  Vol.66. – Р. 1194  1208. DOI:10.1080/17476933.2021.1885386 
10 
Operators Whose Resolvents Have Convolution Representations and their Spectral Analysis 
DOI 10.1007/s10958020051674 
In this paper, we study spectral decompositions with respect to a system of generalized eigenvectors of secondorder differential operators on an interval whose resolvents possess convolution representations. We obtain the convolution representation of resolvents of secondorder differential operators on an interval with integral boundary conditions. Then, using the convolution generated by the initial differential operator, we construct the Fourier transform. A connection between the convolution operation in the original functional space and the multiplication operation in the space of Fourier transforms is established. Finally, the problem on the convergence of spectral expansions generated by the original differential operator is studied. Examples of convolutions generated by operators are also presented. © 2021, Springer Science+Business Media, LLC, part of Springer Nature. 
Kanguzhin B.E.. Operators Whose Resolvents Have Convolution Representations and their Spectral Analysis // Journal of Mathematical Sciences (United States). – 2021. .  Vol.252. – Р. 384  398. DOI: 10.1007/s10958020051674 
11 
Identification of mathematical model of bacteria population under the antibiotic influence 
DOI 10.1515/jiip20170102 
This work is devoted to the identification of a mathematical model of bacteria population under the antibiotic influence, based on the solution of the corresponding inverse problems. These problems are solved by the gradient method, genetic algorithm and NelderMead method. Calculations are made using model and real data. © 2020 De Gruyter. 
Serovajsky S., Nurseitov D., Kabanikhin S., Azimov A., Ilin A., Islamov R.. Identification of mathematical model of bacteria population under the antibiotic influence // Journal of Inverse and IllPosed Problems. – 2021.  Vol.26. – Р. 384  398. DOI: 10.1515/jiip20170102. 
12 
On One Algorithm To Find a Solution to a Linear TwoPoint Boundary Value Problem 
DOI 10.1134/S1995080221030173 
Abstract: A twoparameter family of algorithms for finding an approximatesolution to a linear twopoint boundary value problem for a systemof ordinary differential equations is offered. The convergenceconditions for the algorithms are obtained. The necessary andsufficient coefficient conditions for the wellposedness ofconsidered problem are established. © 2021, Pleiades Publishing, Ltd. 
Temesheva S.M., Dzhumabaev D.S, Kabdrakhova S.S. . On One Algorithm To Find a Solution to a Linear TwoPoint Boundary Value Problem // Lobachevskii Journal of Mathematics. – 2021.  Vol.42. – Р. 606  612. DOI: 10.1134/S1995080221030173. 
13 
The classical KelvinVoigt problem for incompressible fluids with unknown nonconstant density: Existence, uniqueness and regularity 
DOI 10.1088/13616544/abe51e 
The classical KelvinVoigt equations for incompressible fluids with nonconstant density are investigated in this work. To the associated initialvalue problem endowed with zero Dirichlet conditions on the assumed Lipschitzcontinuous boundary, we prove the existence of weak solutions: velocity and density. We also prove the existence of a unique pressure. These results are valid for d ∈ {2, 3, 4}. In particular, if d ∈ {2, 3}, the regularity of the velocity and density is improved so that their uniqueness can be shown. In particular, the dependence of the regularity of the solutions on the smoothness of the given data of the problem is established. © 2021 IOP Publishing Ltd & London Mathematical Society. 
Antontsev S.N., De Oliveira H.B., Khompysh K. . The classical KelvinVoigt problem for incompressible fluids with unknown nonconstant density: Existence, uniqueness and regularity // Nonlinearity. – 2021.  Vol.34. – Р. 3083 – 3111. DOI: 10.1088/13616544/abe51e. 
14 
KelvinVoigt equations with anisotropic diffusion, relaxation and damping: Blowup and large time behavior 
DOI 10.3233/ASY201597 
A nonlinear initial and boundaryvalue problem for the KelvinVoigt equations with anisotropic diffusion, relaxation and absorption/damping terms is considered in this work. The global and local unique solvability of the problem was established in (J. Math. Anal. Appl. 473(2) (2019) 11221154). In the present work, we show how all the anisotropic exponents of nonlinearity and all anisotropic coefficients should interact with the problem data for the solutions of this problem display exponential and polynomial timedecays. We also establish the conditions for the solutions of this problem to blowup in a finite time in three different cases: problem without convection, full anisotropic problem, and the problem with isotropic relaxation. © 2021  IOS Press. All rights reserved. 
Antontsev S.N., De Oliveira H.B., Khompysh K. . KelvinVoigt equations with anisotropic diffusion, relaxation and damping: Blowup and large time behavior // Asymptotic Analysis. – 2021.  Vol.121. – Р. 125  157. DOI: 10.3233/ASY201597 
15 
Punctual Dimension of Algebraic Structures in Certain Classes 
DOI: 10.1134/S1995080221040089 
A nonlinear initial and boundaryvalue problem for the KelvinVoigt equations with anisotropic diffusion, relaxation and absorption/damping terms is considered in this work. The global and local unique solvability of the problem was established in (J. Math. Anal. Appl. 473(2) (2019) 11221154). In the present work, we show how all the anisotropic exponents of nonlinearity and all anisotropic coefficients should interact with the problem data for the solutions of this problem display exponential and polynomial timedecays. We also establish the conditions for the solutions of this problem to blowup in a finite time in three different cases: problem without convection, full anisotropic problem, and the problem with isotropic relaxation. © 2021  IOS Press. All rights reserved. 
Dorzhieva, M.V.,Issakhov, A.A., Kalmurzayev, B.S., Kornev, R.A., Kotov, M.V. . Punctual Dimension of Algebraic Structures in Certain Classes // Lobachevskii Journal of Mathematics. – 2021.  Vol. 42. – Р. 716  725. DOI: 10.1134/S1995080221040089 
16 
Noncommutative analogues of weak compactness criteria in symmetric spaces 
DOI 10.1007/s43036021001292 
The aim of this work is to study M. Nowak’s and K. M. Chong’s weak compactness criteria in noncommutative Orlicz spaces. In addition, we obtain a noncommutative analogue of Pełczyński’s property (V) in Orlicz spaces of τmeasurable operators. We also obtain Kolmogorov’s compactness criterion in terms of conditional expectations in noncommutative symmetric spaces. © 2021, This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply. 
Nessipbayev Y., Sukochev F., Tulenov K. . Noncommutative analogues of weak compactness criteria in symmetric spaces // Advances in Operator Theory. – 2021.  Vol. 6. DOI: 10.1007/s43036021001292 
17 
The boundedness of the Hilbert transformation from one rearrangement invariant Banach space into another and applications 
DOI 10.1016/j.bulsci.2020.102943 
In this paper, we study the boundedness of the Hilbert transformation in Lorentz function spaces, thereby complementing classical results of Boyd. We also characterize the optimal range of a triangular truncation operator in SchattenLorentz ideals. These results further entail sharp commutator estimates and applications to operator Lipschitz functions in SchattenLorentz ideals. © 2020 Elsevier Masson SAS 
Zanin D., Sukochev F., Tulenov K. . The boundedness of the Hilbert transformation from one rearrangement invariant Banach space into another and applications // Bulletin des Sciences Mathematiques. – 2021.  Vol. 167. DOI: 10.1016/j.bulsci.2020.102943 
18 
OPTIMAL REARRANGEMENTINVARIANT BANACH FUNCTION RANGE FOR THE HILBERT TRANSFORM 
DOI 10.32523/20779879202112290103 
We describe the optimal rearrangementinvariant Banach function range of the classical Hilbert transform acting on a rearrangementinvariant Banach function space. We also show the existence of such optimal range for the Lorentz and Marcinkiewicz spaces. © 2021. The L.N. Gumilyov Eurasian National University. 
Tulenov K. . OPTIMAL REARRANGEMENTINVARIANT BANACH FUNCTION RANGE FOR THE HILBERT TRANSFORM // Eurasian Mathematical Journal. – 2021.  Vol. 12. – Р. 01 – 23. DOI: 10.32523/20779879202112290103 
19 
Asymptotic behavior of the solution of a singularly perturbed general boundary value problem with boundary jumps 
DOI 10.1002/num.22719 
It is known that the study of boundary value and mixed problems for integrable linear equations encounters significant difficulties of a fundamental nature. Exceptions are problems with boundary conditions of a special type, which are often called integrable or linearizable. The purpose of this article is to study the asymptotic behaviors of solutions of singularly perturbed general boundary value problems with boundary jumps for higherorder equations. Using the Schlesinger–Birghof theorem, we constructed a fundamental system of solutions of a homogeneous perturbed equation of conditionally stable type in the critical case. Initial boundary functions are constructed based on the fundamental system of solutions. An analytical representation is found, the existence and uniqueness of a solution to this boundary value problem are proved. Asymptotic estimates of the solution and its derivatives are derived from the analytical representation of the solution of the given boundary value problem. The limit passage of solution of the perturbed problem to the solution of the unperturbed problem is proved. The conditions of the existence of jumps are found. The values of boundary jumps are determined. As a result, a class of boundary value problems is highlighted that has possessing of phenomenon of boundary jumps. © 2020 Wiley Periodicals LLC 
Nurgabyl D.N., Uaissov A.B. . Asymptotic behavior of the solution of a singularly perturbed general boundary value problem with boundary jumps // Numerical Methods for Partial Differential Equations. – 2021.  Vol. 37. – Р. 2375  2392. DOI: 10.1002/num.22719. 
20 
Integral Representations of Vector Functions Based on the Parametrix of FirstOrder Elliptic Systems 
DOI 10.1134/S0965542521030143 
Abstract: Generalized integrals are introduced with kernels depending on the difference of the arguments taken over a domain and a smooth contour, the boundary of this domain. These kernels arise as parametrixes of firstorder elliptic systems with variable coefficients. Using such integrals (with complex density over the domain and real density over the contour), representations of vector functions that are smooth in the closed domain are described. The Fredholmity of the representation obtained in the corresponding Banach spaces is established. © 2021, Pleiades Publishing, Ltd. 
Soldatov A.P., Otelbaev M.. Integral Representations of Vector Functions Based on the Parametrix of FirstOrder Elliptic Systems // Computational Mathematics and Mathematical Physics. – 2021.  Vol. 61. – Р. 964  973. DOI: 10.1002/num.22719. 
21 
On Green’s function of Cauchy–Dirichlet problem for hyperbolic equation in a quarter plane 
DOI 10.1186/s13661021015443 
The definition of a Green’s function of a Cauchy–Dirichlet problem for the hyperbolic equation in a quarter plane is given. Its existence and uniqueness have been proven. Representation of the Green’s function is given. It is shown that the Green’s function can be represented by the Riemann–Green function. © 2021, The Author(s). 
Sadybekov M., Derbissaly B. .On Green’s function of Cauchy–Dirichlet problem for hyperbolic equation in a quarter plane // Boundary Value Problems. – 2021.  Vol. 2021. DOI: 10.1186/s13661021015443. 
22 
Boundary Conditions of Volume Hyperbolic Potential in a Domain with Curvilinear Boundary 
DOI 10.1007/9783030692926_20 
A onedimensional volume hyperbolic potential in a domain with curvilinear boundaries is studied. As a kernel of the hyperbolic potential the fundamental solution of the Cauchy problem is chosen. It is wellknown that in this case the volume hyperbolic potential satisfies homogeneous initial conditions. The boundary conditions to which the hyperbolic potential satisfies at lateral boundaries of the domain are constructed. It is shown that the formulated initialboundary value problem has the unique classical solution. © 2021, Springer Nature Switzerland AG. 
Sadybekov M., Derbissaly B. . Boundary Conditions of Volume Hyperbolic Potential in a Domain with Curvilinear Boundary // Springer Proceedings in Mathematics and Statistics. – 2021.  Vol. 351.  Р. 257 – 271. DOI: 10.1007/9783030692926_20 
23 
Optimal Cubature Formulas on Classes of Periodic Functions in Several Variables 
DOI 10.1134/S0081543821010028 
Abstract: We establish sharp order estimates for the error of optimal cubature formulas on the Nikol’skii–Besov and Lizorkin–Triebel type spaces, $$B^{s\,\mathtt{m}}_{p\,q}(\mathbb T^m)$$ and $$L^{s\,\mathtt{m}}_{p\,q}(\mathbb T^m)$$, respectively, for a number of relations between the parameters $$s$$, $$p$$, $$q$$, and $$\mathtt{m}$$ ($$s=(s_1,\dots,s_n)\in\mathbb R^n_+$$, $$1\leq p,q\leq\infty$$, $$\mathtt{m}=(m_1,\dots,m_n)\in{\mathbb N}^n$$, $$m=m_1+\dots+m_n$$). Lower estimates are proved via Bakhvalov’s method. Upper estimates are based on Frolov’s cubature formulas. © 2021, Pleiades Publishing, Ltd. 
Bazarkhanov D.B. . Optimal Cubature Formulas on Classes of Periodic Functions in Several Variables // Proceedings of the Steklov Institute of Mathematics. – 2021.  Vol. 312.  Р. 16  36. DOI: 10.1134/S0081543821010028 
24 
Mixed convection in a channel with buoyancy force over backward and forward facing steps: The effects of inclination and geometry 
DOI 10.1016/j.csite.2021.101152 
This paper presents the computational results of heat transfer for a 2D laminar flow with different channel tilts with forward facing step and backward facing step, taking into account buoyancy forces for various bottom wall lengths. The inclination angle influence on the distribution of velocity and temperature is investigated. The validated numerical algorithm was used to the problem forward and backward facing steps with buoyancy force and at various tilt angles. From the obtained numerical results, it can be noticed that the length of the lower part of the channel has a very strong effect on the flow fluctuation and temperature distribution over the entire channel. It should be noticed that the tilt angle also has a very strong effect on the distribution of flow and temperature. Thus, taking into account the buoyancy force changes the shape of the main recirculation region, but at the same time, regardless of the different tilt angles, the number of vortices does not change, but only the size of the vortices changes. It should also be noticed that when the buoyancy force is taken into account, cooling occurs more efficiently in the middle of the channel. © 2021 
Issakhov A., Zhandaulet Y., Abylkassyomova A., Sakypbekova M., Issakhov A.. Mixed convection in a channel with buoyancy force over backward and forward facing steps: The effects of inclination and geometry // Case Studies in Thermal Engineering. – 2021.  Vol. 26. DOI: 10.1016/j.csite.2021.101152 
25 
Principal frequency of psubLaplacians for general vector fields 
DOI 10.4171/ZAA/1674 
In this paper, we prove the uniqueness and simplicity of the principal frequency (or the first eigenvalue) of the Dirichlet psubLaplacian for general vector fields. As a byproduct, we establish the Caccioppoli inequalities and also discuss the particular cases on the Grushin plane and on the Heisenberg group. © European Mathematical Society 
Ruzhansky M., Sabitbek B., Suragan D. . Principal frequency of psubLaplacians for general vector fields // Zeitschrift fur Analysis und ihre Anwendung. – 2021.  Vol. 40. – Р. 97 – 109. DOI: 10.4171/ZAA/1674 
26 
Geometric Hardy Inequalities on Starshaped Sets 
ISSN 09446532 
We present geometric Hardy inequalities on starshaped sets in Carnot groups. Also, we obtain geometric Hardy inequalities on halfspaces for general vector fields. © Heldermann Verlag 
Ruzhansky M., Sabitbek B., Suragan D. . Geometric Hardy Inequalities on Starshaped Sets s // Journal of Convex Analysis. – 2021.  Vol. 28.

27 
Sobolev, Hardy, Gagliardo–Nirenberg, and Caffarelli–Kohn–Nirenbergtype inequalities for some fractional derivatives

DOI 10.1007/s43037020000974 
In this paper, we show different inequalities for fractionalorder differential operators. In particular, the Sobolev, Hardy, Gagliardo–Nirenberg, and Caffarelli–Kohn–Nirenbergtype inequalities for the Caputo, Riemann–Liouville, and Hadamard derivatives are obtained. In addition, we show some applications of these inequalities. © 2020, Tusi Mathematical Research Group (TMRG). 
Kassymov A., Ruzhansky M., Tokmagambetov N., Torebek B.T.. Sobolev, Hardy, Gagliardo–Nirenberg, and Caffarelli–Kohn–Nirenbergtype inequalities for some fractional derivatives // Banach Journal of Mathematical Analysis. – 2021.  Vol. 15. DOI: 10.1007/s43037020000974 
28 
Cauchy type problems for fractional differential equations 
DOI 10.1080/10652469.2021.1900174 
While it is known that one can consider the Cauchy problem for evolution equations with Caputo derivatives, the situation for the initial value problems for the Riemann–Liouville derivatives is less understood. In this paper, we propose new type initial, inner, and innerboundary value problems for fractional differential equations with the Riemann–Liouville derivatives. The results on the existence and uniqueness are proved, and conditions on the solvability are found. The wellposedness of the new type of initial, inner, and innerboundary conditions is also discussed. Moreover, we give explicit formulas for the solutions. As an application fractional partial differential equations for general positive operators are studied. © 2021 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group. 
Karimov E., Ruzhansky M., Tokmagambetov N. . Sobolev, Hardy, Gagliardo–Nirenberg, and Caffarelli–Kohn–Nirenbergtype inequalities for some fractional derivatives // Integral Transforms and Special Functions. – 2021. DOI: 10.1080/10652469.2021.1900174 
29 
Direct and inverse problems for timefractional pseudoparabolic equations 
DOI 10.2989/16073606.2021.192832 
The purpose of this paper is to establish the solvability results to direct and inverse problems for timefractional pseudoparabolic equations with the selfadjoint operators. We are especially interested in proving existence and uniqueness of the solutions in the abstract setting of Hilbert spaces. © 2021 NISC (Pty) Ltd. 
Serikbaev D., Ruzhansky M., Tokmagambetov N., Torebek B.T. . Direct and inverse problems for timefractional pseudoparabolic equations // Quaestiones Mathematicae. – 2021. DOI: 2989/16073606.2021.192832 
30 
Fractional KleinGordon equation with singular mass 
DOI 10.1016/j.chaos.2020.110579 
We consider a spacefractional wave equation with a singular mass term depending on the position and prove that it is very weak wellposed. The uniqueness is proved in some appropriate sense. Moreover, we prove the consistency of the very weak solution with classical solutions when they exist. In order to study the behaviour of the very weak solution near the singularities of the coefficient, some numerical experiments are conducted where the appearance of a wall effect for the singular masses of the strength of δ2 is observed. © 2020 
Altybay A., Ruzhansky M., Tokmagambetov N., Sebih M.E. . Fractional KleinGordon equation with singular mass // Chaos, Solitons and Fractals. – 2021.  Vol. 143. DOI: 10.1016/j.chaos.2020.110579 
31 
Fractional Schrödinger Equation with Singular Potentials of Higher Order 
DOI 10.1016/S00344877(21)000161 
In this paper the spacefractional Schrödinger equations with singular potentials are studied. Delta like or even higherorder singularities are allowed. By using the regularising techniques, we introduce a family of ‘weakened’ solutions, calling them very weak solutions. The existence, uniqueness and consistency results are proved in an appropriate sense. Numerical simulations are done, and a particles accumulating effect is observed in the singular cases. From the mathematical point of view a “splitting of the strong singularity” phenomena is also observed. © 2021 Polish Scientific Publishers 
Altybay A., Ruzhansky M., Tokmagambetov N., Sebih M.E. . Fractional Schrödinger Equation with Singular Potentials of Higher Order // Reports on Mathematical Physics. – 2021.  Vol. 87. – Р. 129 – 144. DOI: 10.1016/S00344877(21)000161 
32. 
Electrodeposition of polymer electrolyte in nanostructured electrodes for enhanced electrochemical performance of thinfilm Liion microbatteries 
DOI: 10.1016/j.jpowsour.2016.11.078 (IF=7,467; Q1 in Chemistry, Physical). 
ABSTRACT : We report that electrodeposition of polymer electrolyte in nanostructured electrodes has a strong influence on the electrochemical properties of thinfilm Liion microbatteries. Electropolymerization of PMMAPEG (polymethyl methacrylatepolyethylene glycol) was carried out on both the anode (selfsupported titania nanotubes) and the cathode (porous LiNi_{0.5}Mn_{1.5}O_{4})_{ }by cyclic voltammetry and the resulting electrodeelectrolyte interface was examined by scanning electron microscopy. The electrochemical characterizations performed by galvanostatic experiments reveal that the capacity values obtained at different Crates are doubled when the electrodes are completely filled by the polymer electrolyte.

Salian, G.D., Lebouin, C., Demoulin, A., Lepihin, M.S., Maria, S., Galeyeva, A.K., Kurbatov, A.P., Djenizian, T. Electrodeposition of polymer electrolyte in nanostructured electrodes for enhanced electrochemical performance of thinfilm Liion microbatteries. //Journal of Power Sources. – 2017. – Vol.340. – P.242236 DOI: 10.1016/j.jpowsour.2016.11.078

33 
Numerical heat transfer simulation of solargeothermal hybrid source heat pump in Kazakhstan climates 
DOI:10.1016/j.renene.2019.06.102 (SJR 0,18 ; Q4, Procentile 15 (General Chemistry)). 
A numerical energy balance model has been proposed in this work for predicting the thermal performances of solargeothermal hybrid source heat pumps in winter climates of Kazakhstan. The numerical simulation was performed for the year round continental climate conditions. The energy balance model has been developed based on first law of thermodynamics. The proposed heat pump configuration is working in solar mode during sunshine hours and geothermal model during off sunshine hours. Moreover, the system is operating in solargeothermal hybrid model to meet the evaporator load during insufficient availability of solar and geothermal sources. The energy performance comparison between conventional geothermal source and solarground hybrid source heat pump configurations has been made. The influences of solar intensity, ambient temperature, heat pump operating temperatures are discussed. © 2018 Bulgarian Academy of Sciences.

@article{article, author = {Yerdesh, Yelnar and Abdulina, Zarina and Aliuly, Abdurashid and Belyayev, Yerzhan and Mohanraj, M. and Kaltayev, Aidarkhan}, year = {2019}, month = {06}, pages = {}, title = {Numerical simulation on solar collector and cascade heat pump combi water heating systems in Kazakhstan climates}, volume = {145}, journal = {Renewable Energy}, doi = {10.1016/j.renene.2019.06.102} } 
34 
Modelling of turbulence energy decay based on hybrid methods 
DOI: 10.1108/EC1120160395 0.313; Q3; Procentile 63 (General Engineering), 52 (Computational Theory and Mathematics), 47 (Computer Science Applications). 35 (Software))

The purpose of this study is to present an exact and fastcalculated algorithm for the modelling of turbulent energy decay based on two different methods: finitedifference and spectral methods. Design/methodology/approach: The filtered threedimensional nonstationary Navier–Stokes equation is used for simulating the turbulent process. The problem is solved using hybrid methods, where the equation of motion is solved using finite difference methods in combination with cyclic pentadiagonal matrix, which allowed to reach high order of accuracy, and Poisson equation is solved using the spectral methods, which is proposed to speed up the solution procedure. For validation of the given algorithm, the turbulent characteristics were compared with the exact solution of the classical Taylor and Green vortex problem, showing good agreement. Findings: The proposed method shows high computational efficiency and good estimation quality. A numerical algorithm for solving nonstationary threedimensional Navier–Stokes equations for modelling of isotropic turbulence decay using hybrid methods was developed. The simulation’s turbulence characteristics show good agreements with analytical solution. The developed numerical algorithm can be used for simulation of turbulence decay with different values of viscosity. Originality/value: An efficient algorithm for simulation of turbulence processes depending on the properties of the viscosity was developed. © 2018, Emerald Publishing Limited. 
@article{article, author = {Abdibekova, Aigerim and Zhakebayev, Dauren and Abdigaliyeva, Akmaral and Zhubat, Kuanysh}, year = {2018}, month = {08}, pages = {}, title = {Modelling of turbulence energy decay based on hybrid methods}, volume = {35}, journal = {Engineering Computations}, doi = {10.1108/EC1120160395} } 
35 
Mathematical modeling of electron irradiation of oil 
DOI: 10.1108/EC1220160419 (SJR 0.313; Q3; Procentile 63 (General Engineering), 52 (Computational Theory and Mathematics), 47 (Computer Science Applications). 35 (Software))

Purpose: This paper aims to propose a mathematical model and numerical modeling to study the behavior of low conductive incompressible multicomponent hydrocarbon mixture in a channel under the influence of electron irradiation. In addition, it also aims to present additional mechanisms to study the radiation transfer and the separation of the mixture’s components. Design/methodology/approach: The threedimensional nonstationary Navier–Stokes equation is the basis for this model. The Adams–Bashforth scheme is used to solve the convective terms of the equation of motion using a fourthorder accuracy fivepoint elimination method, and the viscous terms are computed with the Crank–Nicolson method. The Poisson equation is solved with the matrix sweep method and the concentration and electron irradiation equations are solved with the Crank–Nicolson method too. Findings: It shows high computational efficiency and good estimation quality. On the basis of numerical results of mathematical model, the effect of the separation of mixture to fractions with various physical characteristics was obtained. The obtained results contribute to the improvement of technologies for obtaining highquality oil products through oil separation into light and heavy fractions. Mathematical model is approbated based on test problem, and has good agreement with the experimental data. Originality/value: The constructed mathematical model makes developing a methodology for conducting experimental studies of this phenomenon possible. © 2018, Emerald Publishing Limited. 
@article{article, author = {Kizbayev, Assylzhan and Zhakebayev, Dauren and Abdibekov, Ualikhan and Khikmetov, Askar}, year = {2018}, month = {08}, pages = {}, title = {Mathematical modeling of electron irradiation of oil}, volume = {35}, journal = {Engineering Computations}, doi = {10.1108/EC1220160419} } 
36 
Hyperbolicelliptic model for surface wave in a prestressed incompressible elastic halfspace 
DOI: 10.1016/j.mechrescom.2018.07.006 (SJR 0.681; Q3; Procentile 76 (Civil and Structural Engineering), 76 (Mechanical Engineering), 72 (Mechanics of Materials), 69 (Condensed Matter Physics), 66 (General Materials Science)) 
The paper aims at derivation of the asymptotic model for surface wave propagating in a prestressed incompressible elastic halfspace, subject to prescribed surface loading. The approach relies on the slowtime perturbation procedure, extending the previously known hyperbolicelliptic formulations for surface waves in compressible linearly elastic solids. Within the derived model, the decay away from the surface is governed by a pseudostatic elliptic equation, whereas wave propagation is described by a hyperbolic equation on the surface. The effect of prestress, namely, the principal Cauchy stress σ_{2}, is investigated. Finally, an illustrative example of the Lamb problem is considered, demonstrating the efficiency of the approach. 
@article{article, author = {Khajiyeva, Lelya and Prikazchikov, D. and Prikazchikova, Ludmila}, year = {2018}, month = {07}, pages = {}, title = {Hyperbolicelliptic model for surface wave in a prestressed incompressible elastic halfspace}, volume = {92}, journal = {Mechanics Research Communications}, doi = {10.1016/j.mechrescom.2018.07.006} } 
37 
The effect of gas and fluid flows on nonlinear lateral vibrations of rotating drill strings 
DOI: 10.1016/j.cnsns.2017.12.008 (SJR 1.159; Q1; Procentile 97 (Applied Mathematics), 96 (Numerical Analysis), 94 (Modeling and Simulation))

In this work we develop nonlinear mathematical models describing coupled lateral vibrations of a rotating drill string under the effect of external supersonic gas and internal fluid flows. An axial compressive load and a torque also affect the drill string. The mathematical models are derived by the use of Novozhilov's nonlinear theory of elasticity with implementation of Hamilton's variation principle. expressions for the gas flow pressure are determined according to the piston theory. The fluid flow is considered as added mass inside the curved tube of the drill string. Using an algorithm developed in the Mathematica computation program on the basis of the Galerkin approach and the stiffness switching method the numerical solution of the obtained approximate differential equations is found. Influences of the external loads, drill string angular speed of rotation, parameters of the gas and fluid flows on the drill string vibrations are shown. 
@article{article, author = {Khajiyeva, Lelya and Kudaibergenov, Askar and Kudaibergenov, Askat}, year = {2017}, month = {12}, pages = {}, title = {The effect of gas and fluid flows on nonlinear lateral vibrations of rotating drill strings}, volume = {59}, journal = {Communications in Nonlinear Science and Numerical Simulation}, doi = {10.1016/j.cnsns.2017.12.008} } 
38 
Detailed Comparative Analysis of Interaction of a Supersonic Flow with a Transverse Gas Jet at High Pressure Ratios 
DOI: 10.1134/S1063784219100049 (SJR 0.297; Q4; Procentile 23 (Physics and Astronomy (miscellaneous))

Abstract: The interaction of a 3D supersonic turbulent gas flow with a transverse sonic jet injected from the wall has been studied in detail both numerically and experimentally. However, the main drawback of such studies is the lack of detailed description of formation and propagation of vortex structures for moderate and large parameters n (ratio of pressure in the jet to pressure in the flow). Analysis performed in this study is aimed at revelation and detailed explanation of mechanisms of formation of vortices behind the injected sonic jet in a supersonic oncoming flow depending on n for improving the effectiveness of mixing of the jet with the flow. As initial equations, we have used 3D Favreaveraged Navier–Stokes equations closed by the k–ω model of turbulence; these equations are solved using the algorithm based on the essentially nonoscillatory scheme of the third approximation order. We have demonstrated the presence of the following vortex structures known from a number of theoretical publications: two oppositely rotating vortices in front of the jet, a horseshoe vortex; and two pairs formed in the mixing zone between the jet and the flow (one in the wake behind the jet and the other on the lateral line of the jet). We have determined the pressure ratios for which extra pairs of vortices appear (one pair emerges at the Mach disk edge as a result of interaction of a retarded flow of the jet behind the Mach disk with a highvelocity ascending flow behind the barrel and the other pair is formed due to the interaction of the ascending jet flow with the incoming main gas flow). As a result of comparative analysis, the pressure ratios for which a clear pattern of additional horn vortices is observed near the wall in the region behind the jet, have been determined. The dependence of the slope of the bow shock on the pressure ratio has been plotted. It is found that the pressure distribution at the wall in front of the jet in the symmetry plane is in satisfactory agreement with experimental data. © 2019, Pleiades Publishing, Ltd. 
@article{article, author = {Beketaeva, A. and Bruel, Pascal and Naimanova, A.}, year = {2019}, month = {10}, pages = {14301440}, title = {Detailed Comparative Analysis of Interaction of a Supersonic Flow with a Transverse Gas Jet at High Pressure Ratios}, volume = {64}, journal = {Technical Physics}, doi = {10.1134/S1063784219100049} } 
39 
Numerical study of the interaction of vortex ring and background turbulence 
DOI: 10.1134/S0869864319060040 (SJR 0.473; Q3; Procentile 43 (Aerospace Engineering) 39 (Energy Engineering and Power Technology) 35 (Modeling and Simulation)

The article presents the results of direct numerical simulation of a turbulent vortex ring with a moderate Reynolds number, interacting with the field of external turbulent fluctuations. The ring is formed by buoyancy forces from a sphericallyshaped volume of elevated temperature. We consider the interaction of the ring with the field of turbulent fluctuations (of temperature and velocity) located in the form of a horizontal layer in front of the ring. The effect of separation of vortices by the sign of vorticity in a layer of fluctuations (with respect to the rotation of the ring) during the passage of the ring through this field has been found. This effect first causes the ring to slow down when the fluctuations pass outside the ring, and then to accelerate with decreasing radius when the fluctuations pass through its center. Due to buoyancy effects in the fluctuations layer, the separation of hot and cold air occurs. During the interaction with the ring, this leads to the accumulation of the reduced temperature in the core of the ring, while the increased temperature is accumulated near the axis. This temperature distribution results in a baroclinic vorticity generation leading to the reduction of the ring radius. © 2019, M.Yu. Hrebtov, M.S. Bobrov, D.B. Zhakebaev, and K.K. Karzhaubaev. 
@article{article, author = {Hrebtov, Michael and Bobrov, M. and Zhakebaev, D. and Karzhaubaev, K.}, year = {2019}, month = {11}, pages = {821836}, title = {Numerical study of the interaction of vortex ring and background turbulence}, volume = {26}, journal = {Thermophysics and Aeromechanics}, doi = {10.1134/S0869864319060040} } 
40 
Numerical mass transfer studies in case of convective flows occurrence in isothermal ternary gas mixtures 
DOI: 10.1615/ComputThermalScien.2018024724 (SJR 0.245; Q3; Procentile 35 (Energy Engineering and Power Technology), 34 (Computational Mathematics), 31 (Fluid Flow and Transfer Processes), 21 (Surfaces and Interfaces))

Isothermal diffusion mixing in multicomponent gas mixtures at various pressures is studied numerically in He + Ar – N_{2} and CH_{4} + Ar – N_{2} . It is shown that in systems where the diffusion coefficients differ significantly from each other, nonlinear distributions of the concentrations of the components arise with increasing pressure, which leads to a nonmonotonous density distribution of the gas mixture in the computational domain. The resulting nonlinearity is the cause of the formation of convective structured flows. The transition from the diffusion to the convective regime is characterized by a significant increase in the average mixing rate. The dynamics of the development of convective currents and their structure is investigated. © 2019 by Begell House, Inc. 
@article{article, author = {Kossov, Vladimir and Fedorenko, Olga and Zhakebayev, Dauren}, year = {2018}, month = {01}, pages = {}, title = {NUMERICAL MASS TRANSFER STUDIES IN CASE OF CONVECTIVE FLOWS OCCURRENCE IN ISOTHERMAL TERNARY GAS MIXTURES}, volume = {11}, journal = {Computational Thermal Sciences: An International Journal}, doi = {10.1615/ComputThermalScien.2018024724} } 
41 
Peculiarities of the rise of structured formations at the boundary of the change of the regimes “diffusionconcentration convection” at an isothermal mixing of a binary mixture equally diluted by the third component 
DOI: 10.1134/S0869864319010049 (SJR 0.245; Q3; Procentile 43 (Aerospace Engineering) 39 (Energy Engineering and Power Technology) 35 (Modeling and Simulation)) 
The peculiarities of the change of the regimes “diffusiongravitational concentration convection” have been studied numerically. It is shown that at a certain mixture composition and pressure, substantially nonlinear distributions of the concentrations of components arise at the expense of the difference in the diffusion coefficients of components, which lead to the corresponding distribution of the gaseous mixture density. This is a reason for the formation of structured formations and the rise of convective instability in the mixtures under study. The time of the rise of convective flows and of the mean velocity of the components transfer has been estimated. Obtained results are compared with experimental data. © 2019, V.N. Kossov, O.V. Fedorenko, D.B. Zhakebaev, and A.P. Kizbaev. 
@article{article, author = {Kossov, V. and Fedorenko, O. and Zhakebaev, D. and Kizbaev, A.}, year = {2019}, month = {01}, pages = {2735}, title = {Peculiarities of the rise of structured formations at the boundary of the change of the regimes “diffusionconcentration convection” at an isothermal mixing of a binary mixture equally diluted by the third component}, volume = {26}, journal = {Thermophysics and Aeromechanics}, doi = {10.1134/S0869864319010049} } 
42 
Features of Multicomponent Mass Transfer in Gas Mixtures Containing Hydrocarbon Components 
DOI: 10.1002/ceat.201800604 (SJR 0.403; Q3; Procentile 61 (Industrial and Manufacturing Engineering) 55 (General Chemical Engineering) 53 (General Chemistry)) 
The features characterizing the transition from a diffusion process to the concentration gravitational convection in multicomponent gas mixtures containing hydrocarbon components are analyzed. It is demonstrated that this transition is possible due to the difference in interdiffusion coefficients of components, initial concentration of components, and the pressure of the experiment. A mathematical model for describing the transition from the diffusion process to the concentration gravitational convection in multicomponent gas mixtures is considered. The key feature indicating the transfer from the diffusion regime to the convective one is the curvature of an isoline of the heavy component of the gas mixture. © 2019 WILEYVCH Verlag GmbH & Co. KGaA, Weinheim

@article{article, author = {Kossov, Vladimir and Fedorenko, Olga and Zhakebayev, Dauren}, year = {2019}, month = {02}, pages = {}, title = {Features of multicomponent mass transfer at the boundary of change in the “diffusion ‐ concentration gravitational convection” regime in gas mixtures containing hydrocarbon components}, volume = {42}, journal = {Chemical Engineering & Technology}, doi = {10.1002/ceat.201800604} } 
43 
Mixed LagrangianEulerian Simulation of Interaction between a Shockwave and a Cloud of Water Droplets 
DOI: 10.1134/S1810232820020071 (SJR 0.457; Q3; Procentile 50 (Energy Engineering and Power Technology), 49 (Modeling and Simulation), 48 (Environmental Engineering), 41 (Condensed Matter Physics))

Abstract: We present a numerical model suitable for simulation of shock orblast waves passing through a cloud of water droplets. The model takesinto account the droplet breakup, radiation, and evaporation effects.The gas phase (a mixture of air and water vapor) is solved within anEulerian framework with a set of compressible transport equations. Thedisperse phase (water droplets) is represented by a number of Lagrangianparcels of a specified size and mass distribution. The model has beenverified with experimental data. The results show that for large(millimetersized) droplets with high Weber numbers, the breakup modelis the most important part for accurate representation of thewavedroplet interaction phenomena, while for very fine droplets(1 μmor less), the evaporation effects are the strongest in the shockwavemitigation process. The radiative heat flux increase due to high dropletemissivity is found to be possible for an intermediate size waterdroplets (about 10–20 μm)in the case of continuous heat release. © 2020, Pleiades Publishing, Ltd.

@article{article, author = {Hrebtov, Michael and Bobrov, M. and Zhakebayev, Dauren and Karzhaubayev, K.}, year = {2020}, month = {04}, pages = {254263}, title = {Mixed LagrangianEulerian Simulation of Interaction between a Shockwave and a Cloud of Water Droplets}, volume = {29}, journal = {Journal of Engineering Thermophysics}, doi = {10.1134/S1810232820020071} } 
44 
Laminar to turbulent flow transition inside the boundary layer adjacent to isothermal wall of natural convection flow in a cubical cavity 
DOI: 10.1016/j.ijheatmasstransfer.2020.120822 (SJR 1.713; Q1; Procentile 98 (Fluid Flow and Transfer Processes) 95(Mechanical Engineering) 94(Condensed Matter Physics))

We investigate threedimensional natural convection flow in an airfilled, differentially heated cubical cavity. The vertical wall on the left is heated and the vertical wall on the right is cooled, with the remaining four walls being adiabatic. We performed direct numerical simulations of the natural convection flow using discrete unified gaskinetic scheme (DUGKS), with an improved implementation of boundary conditions. Thin boundary layers are developed along the two isothermal walls. The laminar to turbulent flow transition inside the boundary layers is studied in this paper. The simulations are conducted at three Rayleigh numbers of 1.5×10^{9}, 1.0×10^{10}, 1.0×10^{11} using nonuniform grids with resolution up to 320^{3}. The Prandtl number is fixed at 0.71. We provide a detailed analysis of the transition from laminar to turbulent flow inside the vertical boundary layers and its influence on the rate of heat transfer. Time traces of temperature and velocity, timeaveraged flow field, statistics of fluctuation fields are presented to illustrate distinct behaviors in the laminar and turbulent thermal boundary layer, as well as to determine the transition location at different Ra numbers. The average Nusselt numbers for different Ra numbers are compiled and compared to previous results. A guideline of the resolution requirement is suggested based on the Ra scaling of laminar thermal boundary layer. © 2020 Elsevier Ltd 
@article{article, author = {Wen, Xin and Wang, LianPing and Guo, Zhaoli and Zhakebayev, Dauren}, year = {2021}, month = {03}, pages = {120822}, title = {Laminar to turbulent flow transition inside the boundary layer adjacent to isothermal wall of natural convection flow in a cubical cavity}, volume = {167}, journal = {International Journal of Heat and Mass Transfer}, doi = {10.1016/j.ijheatmasstransfer.2020.120822} } 
45 
The Verhulstlike equations: Integrable OΔE and ODE with chaotic behavior 
DOI: 10.3390/sym11121446 (SJR 0.385; Q2; Procentile 90 (General Mathematics) 77 (Computer Science miscellaneous) 66 (Physics and Astronomy (miscellaneous)) 56 (Chemistry (miscellaneous)) 
In this paper, we study various variants of Verhulstlike ordinary differential equations (ODE) and ordinary difference equations (OΔE). Usually Verhulst ODE serves as an example of a deterministic system and discrete logistic equation is a classic example of a simple system with very complicated (chaotic) behavior. In our paper we present examples of deterministic discretization and chaotic continualization. Continualization procedure is based on Pade approximants. To correctly characterize the dynamics of obtained ODE we measured such characteristic parameters of chaotic dynamical systems as the Lyapunov exponents and the Lyapunov dimensions. Discretization and continualization lead to a change in the symmetry of the mathematical model (i.e., group properties of the original ODE and OΔE). This aspect of the problem is the aim of further research. © 2019 by the authors. 
@article{article, author = {Andrianov, Igor and Starushenko, Galina and Kvitka, Sergiy and Khajiyeva, Lelya}, year = {2019}, month = {11}, pages = {1446}, title = {The VerhulstLike Equations: Integrable OΔE and ODE with Chaotic Behavior}, volume = {11}, journal = {Symmetry}, doi = {10.3390/sym11121446} } 
46 
On the dynamics of drilling 
DOI: 10.1016/j.ijengsci.2019.103184 (SJR 2.731; Q1; Procentile 99 (General Engineering) 97 (Mechanical Engineering) 96 (Mechanics of Materials) 92 (General Materials Science)

A multiparameter approach for analysing drilling dynamics is developed. It is oriented to better understanding peculiarities of the process and optimising existing formulations. The procedure is specified for bending vibrations of a rotating elastic beam prestressed simultaneously by an axial compressive force and torque. The appropriately normalized ratio of the bending stiffness to the axial force magnitude is assumed to be small. Two other independent dimensionless parameters correspond to the vibration frequency and the relation between force and torque magnitudes. A classification, similar to that in thin shell theory, is established for free bending vibrations. The associated shortened equations are derived. The asymptotic results are validated by comparison with the dispersion curves calculated from the original full equation of motion. An example of a boundaryvalue problem for a drill string of a finite length is also presented. © 2019 
@article{article, author = {Kaplunov, J. and Khajiyeva, Lelya and Martyniuk, M. and Sergaliyev, A.S.}, year = {2020}, month = {01}, pages = {103184}, title = {On the dynamics of drilling}, volume = {146}, journal = {International Journal of Engineering Science}, doi = {10.1016/j.ijengsci.2019.103184} } 
47 
Development of an information system for storing digitized works of the Almaty Academgorodok research institutes 
DOI: 10.47086/pims.613534 (SJR 0.190; Q3; Procentile 43 (General Engineering))

The present article describes the architecture of the integrated distributed information system created for storing digitized works of employees of Almaty Akademgorodok research institutes (Kazakhstan) and providing access to them using Web technology. Comparative analysis of two data storage systems for storing digitized works, Ceph and GlusterFS, is provided. The description of the software part of the information system is provided which consists of four subsystems: repository of digital objects, subsystem for managing current research information, subsystem of integration of distributed information resources, subsystem of access to distributed information resources based on Web technologies. The relation between the subsystems and their integration is described. The work defines the requirements to the repository of digital objects. The requirements for the repository of digital objects are defined; a comparative analysis of open source software used for these purposes is made. © 2019 Author(s). 
@article{article, author = {TEMİRBEKOV, Nurlan and Baigereyev, Dossan and TEMİRBEKOV, Almas and OMİRZHANOVA, Bakytzhan}, year = {2021}, month = {07}, pages = {}, title = {Development of an Information System for Storing Digitized Works of the Almaty Academgorodok Research Institutes}, journal = {Proceedings of International Mathematical Sciences}, doi = {10.47086/pims.613534} } 
48 
[Architecture of the distributed information system of the Almaty Academgorodok] [Architektura rozproszonego systemu informacyjnego w Kampusie Almaty] 
DOI: 10.15199/48.2019.05.19 (SJR 0.190; Q3; Procentile 22 (Electrical and Electronic Engineering)) 
The paper describes the architecture of an integrated distributed information system, which allows to preserve the result of the intellectual activity of Kazakhstan Engineering Technological University and a number of research institutes, located in the Academgorodok of Almaty, Kazakhstan. A description of the implementation of information system subsystems is presented. © 2019, Wydawnictwo SIGMANOT. All rights reserved. 
@article{article, author = {TEMIRBEKOV, Nurlan and Baigereyev, Dossan}, year = {2019}, month = {05}, pages = {7780}, title = {Architecture of the distributed information system of the Almaty Academgorodok}, volume = {1}, journal = {PRZEGLĄD ELEKTROTECHNICZNY}, doi = {10.15199/48.2019.05.19} } 
49 
Development of a Distributed Information System of the Almaty Academgorodok 
DOI:10.15199/48.2019.05.19 (SJR 0.190; Q3; Procentile 43 (General Engineering)

The present article describes the architecture of an integrated distributed information system used to store and manage digitized works of employees of research institutes of the Almaty Academgorodok. The Ceph open source software object storage network is used as data storage. Testing of the POSIXcompatible CephFS file system abstraction is performed. The software part of the information system consists of four subsystems: repository of digital objects, subsystem for managing current research information, subsystem of integration of distributed information resources, subsystem of access to distributed information resources based on web technologies. The description of the software part of the information system is provided. Integration between the subsystems of the information system is performed. © 2020 School of Science, IHU. All Rights Reserved. 
@article{article, author = {TEMIRBEKOV, Nurlan and Baigereyev, Dossan}, year = {2019}, month = {05}, pages = {7780}, title = {Architecture of the distributed information system of the Almaty Academgorodok}, volume = {1}, journal = {PRZEGLĄD ELEKTROTECHNICZNY}, doi = {10.15199/48.2019.05.19} } 
50 
Parallel CUDA implementation of a numerical algorithm for solving the NavierStokes equations using the pressure uniqueness condition 
DOI: 10.1063/5.0041039 (SJR 0.177; Q4; Procentile 17 (General Physics and Astronomy)) 
In this paper, we study numerical methods for solving the NavierStokes equations in doubly connected domains. Two methods for solving the problem are considered. The first method is based on constructing a difference problem in variables of the stream function and the vortex of velocity using the uniqueness condition for pressure. The numerical solution of the elliptic equation for stream functions is found as the sum of the solutions of two simple problems of an elliptic type. One problem is with homogeneous boundary conditions, and the other is with a homogeneous equation. An alternative approach to solving the problem is the fictitious domain method with the continuation of the least coefficient. This method does not require satisfying the pressure uniqueness condition, and is simple to implement. An important direction in the development of numerical simulation methods is the study of approximate methods for solving problems of mathematical physics in complex multidimensional areas. To solve many applied problems in irregular areas, the fictitious domain method is widely used, which is characterized by a high degree of automation of programming. The main idea of the fictitious domain method is that the problem is solved not in the original complex domain, but in some other, simpler domain. This allows to create software immediately for a fairly wide class of problems with arbitrary computational domains. The possibilities of applying the fictitious domain method to the problems of hydrodynamics in the variables "stream function, vortex of velocity"are considered in many works. In this paper, we study a numerical method for solving the NavierStokes equations in doubly connected domains. A computational finite difference algorithm for solving an auxiliary problem of the fictitious domain method has been developed. The results of numerical modeling of the twodimensional NavierStokes equations by the fictitious domain method with continuation by the lowest coefficient are presented. For this problem, a parallel algorithm was developed using the CUDA architecture, which was tested on various grid dimensions. © 2021 Author(s). 
@inproceedings{inproceedings, author = {Temirbekov, Almas and Baigereyev, Dossan and Temirbekov, Nurlan and Urmashev, Baidaulet and Amantayeva, Aidana}, year = {2021}, month = {02}, pages = {020063}, title = {Parallel CUDA implementation of a numerical algorithm for solving the NavierStokes equations using the pressure uniqueness condition}, volume = {2325}, journal = {AIP Conference Proceedings}, doi = {10.1063/5.0041039} } 
51 
Information systems of integrated machine learning modules on the example of a verbal robot 
DOI: 10.32014/2019.2518170X.173 (SJR 0.323; Q3; Procentile 40 (Geology) 37 (Geotechnical Engineering and Engineering Geology) 
In this work, information systems of integrated machine learning modules have been performed using the verbal robot as an example. Hardware components have been developed, logical components that have been assembled and /or developed to implement an automated verbal robot system, module tracking has been also performed that can track human faces in real time through the OpenCV library and automated services on Jetson TX1 SoC for maneuvering a mobile robot chassis. © 2019, National Academy of Sciences of the Republic of Kazakhstan. All rights reserved. 
@article{article, author = {Kalimoldayev, Maksat and Akhmetzhanov, Maxat and Kunelbayev, Murat}, year = {2019}, month = {12}, pages = {215222}, title = {INFORMATION SYSTEMS OF INTEGRATED MACHINE LEARNING MODULES ON THE EXAMPLE OF A VERBAL ROBOT}, volume = {6}, journal = {NEWS of National Academy of Sciences of the Republic of Kazakhstan}, doi = {10.32014/2019.2518170X.173} } 
52 
Implementation of covering algorithm for the robot with parallel structure 
DOI: 10.32014/2019.2518170X.142 (SJR 0.323; Q3; Procentile 40 (Geology) 37 (Geotechnical Engineering and Engineering Geology) 
The paper considers a 3RPR robot with a parallel structure. One of the main tasks in robotics is to determine the working area of the robot. Algorithms for solving systems of this type are given. The properties and accuracy estimates of the obtained approximations are proved. As an approach to determining the work area, the method of nonuniform coatings was used in this work, which allows one to determine the external and internal approximation of the set of solutions of the system with a given accuracy. © 2019, National Academy of Sciences of the Republic of Kazakhstan. All rights reserved. 
@article{article, author = {Kalimoldayev, M. and Akhmetzhanov, M. and Mukanova, Balgaisha and Azimova, Dinara}, year = {2019}, month = {10}, pages = {197201}, title = {IMPLEMENTATION OF COVERING ALGORITHM FOR THE ROBOT WITH PARALLEL STRUCTURE}, volume = {5}, journal = {NEWS of National Academy of Sciences of the Republic of Kazakhstan}, doi = {10.32014/2019.2518170X.142} } 
53 
NearResonant Regimes of a Moving Load on a PreStressed Incompressible Elastic HalfSpace 
DOI: 10.2478/ama20210005 (SJR 0.314; Q3; Procentile 45 (Mechanical Engineering) 41 (Control and Systems Engineering) 
The article is concerned with the analysis of the problem for a concentrated line load moving at a constant speed along the surface of a prestressed, incompressible, isotropic elastic halfspace, within the framework of the planestrain assumption. The focus is on the nearcritical regimes, when the speed of the load is close to that of the surface wave. Both steadystate and transient regimes are considered. Implementation of the hyperbolicelliptic asymptotic formulation for the surface wave field allows explicit approximate solution for displacement components expressed in terms of the elementary functions, highlighting the resonant nature of the surface wave. Numerical illustrations of the solutions are presented for several material models. 
@article{article, author = {Kudaibergenov, Askar and Kudaibergenov, Askat and Prikazchikov, D.}, year = {2021}, month = {03}, pages = {3036}, title = {NearResonant Regimes of a Moving Load on a PreStressed Incompressible Elastic HalfSpace}, volume = {15}, journal = {Acta Mechanica et Automatica}, doi = {10.2478/ama20210005} } 
54 
Differential evolution algorithm of solving an inverse problem for the spatial Solow mathematical model 
DOI: 10.1515/jiip20200108 (SJR 0.498; Q3; Procentile 54 (Applied Mathematics)) 
The differential evolution algorithm is applied to solve the optimization problem to reconstruct the production function (inverse problem) for the spatial Solow mathematical model using additional measurements of the gross domestic product for the fixed points. Since the inverse problem is illposed the regularized differential evolution is applied. For getting the optimized solution of the inverse problem the differential evolution algorithm is paralleled to 32 kernels. Numerical results for different technological levels and errors in measured data are presented and discussed. © 2020 Walter de Gruyter GmbH, Berlin/Boston 2020. 
@article{article, author = {Kabanikhin, S. and Krivorotko, Olga and Bektemessov, Zh and Bektemesov, Maktagali and Zhang, Shuhua}, year = {2020}, month = {11}, pages = {761774}, title = {Differential evolution algorithm of solving an inverse problem for the spatial Solow mathematical model}, volume = {28}, journal = {Journal of Inverse and Illposed Problems}, doi = {10.1515/jiip20200108} } 
55 
Numerical Solution of the Inverse Pharmacokinetic Problem for the Three Compartment Model 
DOI:10.51301/vest.su.2020.v138.i2.136 (SJR 0.19; Q3; Procentile 43 (Engineering:

This article considers the numerical solution of the inverse pharmacokinetics problem for a threecompartment linear model. First, the article presents some reviews of the pharmacokinetics problem and the threecompartment model. The following describes the formulation of the pharmacokinetics problem for a threecompartment linear model. The direct problem is the Cauchy problem for systems of ordinary differential equations. Solving the direct problem analytically, we find the concentration for the first compartment, since it is the object of the study. The formulation of the inverse problem is reduced to a nonlinear operator equation. For the inverse problem, seven coefficients concentration for the first compartment should be found for some additional information of a given concentration. The inverse problem is reduced to minimizing the objective functional. For the numerical solution, an adaptive search method is used genetic algorithm. The numerical results of this problem are given. © 2020. School of Science, IHU. All rights reserved. 
@article{article, author = {Urmashev, B.A. and Kasenov, S.Y. and Tursynbay, А.Т and Temirbekov, A.N. and Amantayeva, A.B. and Sagimbayeva, L.A.}, year = {2020}, month = {01}, pages = {772779}, title = {Solutions of the inverse problem of pharmacokinetics for the threecompartment model}, volume = {138}, journal = {Vestnik KazNRTU}, doi = {10.51301/vest.su.2020.v138.i2.136} } 
56 
Numerical algorithm for solving the inverse problem for the helmholtz equation 
DOI: 10.1007/9783030122034_20 (SJR 0.16; Q3; Procentile 30 (Mathematics: 18 (Computer Science:

In this paper we consider acoustic equation. The equation by separation of variables is reduced to a boundary value problem for the Helmholtz equation. We consider problem for the Helmholtz equation. We reduce the solution of the operator equation to the problem of minimizing the functional. And we build numerical algorithm for solving the inverse problem. At the end of the article is given the numerical calculations of this problem. © Springer Nature Switzerland AG 2019. 
@inbook{inbook, author = {Shishlenin, Maxim and Kasenov, S. and Askerbekova, Zh}, year = {2019}, month = {01}, pages = {197207}, title = {Numerical Algorithm for Solving the Inverse Problem for the Helmholtz Equation: Methods and Protocols}, isbn = {9781493990856}, doi = {10.1007/9783030122034_20} } 
57 
Numerical Solution of the Inverse Pharmacokinetic Problem for the Three Compartment Model 
DOI: 10.51301/vest.su.2020.v138.i2.136 (SJR 0.19; Q4; Procentile 43 (Engineering:

This article considers the numerical solution of the inverse pharmacokinetics problem for a threecompartment linear model. First, the article presents some reviews of the pharmacokinetics problem and the threecompartment model. The following describes the formulation of the pharmacokinetics problem for a threecompartment linear model. The direct problem is the Cauchy problem for systems of ordinary differential equations. Solving the direct problem analytically, we find the concentration for the first compartment, since it is the object of the study. The formulation of the inverse problem is reduced to a nonlinear operator equation. For the inverse problem, seven coefficients concentration for the first compartment should be found for some additional information of a given concentration. The inverse problem is reduced to minimizing the objective functional. For the numerical solution, an adaptive search method is used genetic algorithm. The numerical results of this problem are given. © 2020 School of Science, IHU. All Rights Reserved. 
@article{article, author = {Urmashev, B.A. and Kasenov, S.Y. and Tursynbay, А.Т and Temirbekov, A.N. and Amantayeva, A.B. and Sagimbayeva, L.A.}, year = {2020}, month = {01}, pages = {772779}, title = {Solutions of the inverse problem of pharmacokinetics for the threecompartment model}, volume = {138}, journal = {Vestnik KazNRTU}, doi = {10.51301/vest.su.2020.v138.i2.136} } 
58 
Determination of the main parameters of the photovoltaic solar module 
DOI: 10.12720/sgce.9.1.162169 (SJR 0.203; Q4; Procentile 25 (Energy: General: Energy), 20 (Environmental Science: General Environmental Science), 19 (Earth and Planetary Sciences: General Earth and Planetary Sciences) 
This article deals with the determination of the main operating parameters of a photovoltaic solar module. In laboratory tests, the study of the dependence of current, voltage and power on time and density of solar radiation, as well as monitoring of environmental parameters: temperature and humidity of the outside air. Analysis of the test results shows that a photoelectric module with an installed capacity of 800 W and a total battery capacity of 800 Ah provides the electric power industry with a daily consumption of 2.0.. 2.2 kWh. The discharge time of the battery varies from 11.7 to 3.5 hours when the average electric load of the consumer changes from 300 to 1000 watts. © 2020 The Authors, published by EDP Sciences. 
@article{article, author = {Urmashev, Baydaulet and Kunelbayev, Murat and Temirbekov, Almas and Kassenov, Syrym and Zhaksylykova, Zhadra and Amenova, Farida}, year = {2020}, month = {01}, pages = {162169}, title = {Determination of the main parameters of the photovoltaic solar module}, journal = {International Journal of Smart Grid and Clean Energy}, doi = {10.12720/sgce.9.1.162169} } 
59. 
Numerical heat transfer simulation of solargeothermal hybrid source heat pump in Kazakhstan climates 
DOI: http://bcc.bas.bg/BCC_Volumes/Volume_50_Special_G_2018/50G_PD_713.56.pdf (SJR 0,18 ; Q4, Procentile 15 (General Chemistry)). 
ABSTRACT : A numerical energy balance model has been proposed in this work for predicting the thermal performances of solargeothermal hybrid source heat pumps in winter climates of Kazakhstan. The numerical simulation was performed for the year round continental climate conditions. The energy balance model has been developed based on first law of thermodynamics. The proposed heat pump configuration is working in solar mode during sunshine hours and geothermal model during off sunshine hours. Moreover, the system is operating in solargeothermal hybrid model to meet the evaporator load during insufficient availability of solar and geothermal sources. The energy performance comparison between conventional geothermal source and solarground hybrid source heat pump configurations has been made. The influences of solar intensity, ambient temperature, heat pump operating temperatures are discussed. © 2018 Bulgarian Academy of Sciences.

Saktashova G., Aliuly A., Belyayev Ye., Mohanraj M., Singh R.M. ,Numerical heat transfer simulation of solargeothermal hybrid source heat pump in Kazakhstan climates. //Bulgarian Chemical Communications. – 2018. – Vol.50. – P.713 DOI: http://bcc.bas.bg/BCC_Volumes/Volume_50_Special_G_2018/50G_PD_713.56.pdf

60. 
Research and developments on solar assisted compression heat pump systems – A comprehensive review (Part A: Modeling and modifications) 
DOI: https://doi.org/10.1016/j.rser.2017.08.022 (SJR 3,52 ; Q1, Procentile 97 (Renewable Energy, Sustainability and the Environment)). 
This paper presents a comprehensive review on research and developments on solar assisted compression heat pump systems, mostly reported during the last two decades. The first part of this paper provides a detailed description of the past efforts on various system configurations, system modeling, enhancement of system performance, modifications in compression heat pump cycles and environment friendly refrigerant options for solar assisted compression heat pump systems. The economical and environmental impacts of solar assisted compression heat pump systems are also described. Further research needs in this important field of solar assisted compression heat pump are listed. The outcomes of this review confirm that, there is a lot of research scope for improving the performance of the system and reducing its initial cost to make it competitive in the global market. The information presented in this paper is shall be highly beneficial for the active researchers working presently on solar assisted compression heat pumps. 
M.Mohanraj, Ye.Belyayev, S.Jayaraj, A.Kaltayev , Numerical heat transfer simulation of solargeothermal hybrid source heat pump in Kazakhstan climates. //Renewable and Sustainable Energy Reviews. – 2018. – Vol.83. – P.90123 
61. 
Dynamics of orthogonal mechanism of vibrating table in view of friction 
DOI: 10.1007/9783319607023_27 (SJR 0,16; Q4, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials)). 
The article deals with dynamics of orthogonal mechanism of vibrating table in view of friction. The emphasis is placed on the influence of friction on the oscillatory motion of the mechanism at its interaction with nonideal energy source. As a result of analysis of numerical solutions of nonlinear motion equation has been determined that the coefficient of sliding friction affects the average angular speed of the driving link, and at that its growth in the first half cycle, when the vibrating table moves upward will cause reduction of the average angular speed, but in the second halfcycle, conversely, will cause increase of this value of angular speed. A coefficient of rolling friction influences on the vibration amplitude of the angular speed of the driving link: with its increasing the amplitude of vibrations increases. Therefore, sliding friction may affect the coefficient of nonuniformity of rotation in two ways, i.e. one half cycle of the mechanism rotation can reduce it and increase it in the other half cycle, and only rolling friction destabilizing effect has been found. © Springer International Publishing AG 2018. 
Iskakov Z.,Bissembayev K.,Jamalov N., Dynamics of orthogonal mechanism of vibrating table in view of friction. //Mechanisms and Machine Science. – 2018. – Vol.52. – P.261268 DOI: 10.1007/9783319607023_27 
62. 
The nonlinear vibrations of orthogonal mechanism of vibrating table in view of friction 
DOI: 10.5194/ms93072018 (SJR 0,31; Q2, Procentile 52(Engineering Industrial and Manufacturing Engineering) 50(Engineering Mechanical Engineering) 48(Engineering Civil and Structural Engineering) 48(Chemical Engineering Fluid Flow and Transfer Processes) 46(Engineering Mechanics of Materials) 45(Engineering Control and Systems Engineering)). 
In dynamics of the orthogonal mechanism of vibrating table main attention is paid to the effect of friction on vibrational proceß when it is interacting with nonideal energy source. As a result of analysis of numerical solutions of nonlinear equation of motion, it is established that in general, for the period of rotation, both sliding friction and rotational friction increase the difference between the maximum and minimum angular velocity values of the driving link (engine shaft) when its average value is constant; both increase the coefficient of nonuniformity of rotation. When there is a cargo of significant maß in the working link, the effect of sliding friction impacts more on the amplitude and frequency of the angular velocity oscillations of the driving link (engine shaft) and on the frequency of the coordinate changes, as well as on the frequencies and maximum and minimum values of the projections of the velocity and acceleration of the hinge linking the work platform to coupler. Similarly, the larger of the load maß, the stronger the effect of rotational friction on the amplitude and frequency of angular velocity oscillations of the driving link (engine shaft). © Author(s) 2018. 
Iskakov Z.,Bissembayev K.,Jamalov N., The nonlinear vibrations of orthogonal mechanism of vibrating table in view of friction. //Mechanical Sciences. – 2018. – Vol.9. – P.307325 DOI: 10.5194/ms93072018 
63. 
Analysis of manipulator grains for overloading largesized goods in the sewing industry 
(SJR 0,2 ; Q3, Procentile 17(Materials Science Polymers and Plastics) 17(Business, Management and Accounting General Business, Management and Accounting) 16(Engineering Industrial and Manufacturing Engineering) ). 
The paper examines the gripper scheme of the manipulator for reloading bulky goods in the garment industry. Grips with a circuit of rods and with cable drives are studied. The results of static calculations for the gripper with pneumatic drive and calculations for strength and rigidity for the threephalanx adaptive gripper are given. In the garment industry for the reloading of largesized cargoes, the use of robot manipulators is complicated by the lack of reliable grippers. To fill this problem, two grasping schemes are proposed in this paper. Some calculations of these two grasping schemes developed by the authors are given. The drive of the first gripper is pneumatic, and the second one is based on flexible traction elements. © 2018 Ivanovo State University of Chemistry and Technology. All rights reserved. 
Usenbekov Z.H.,Temirbekov E.S.,Kaimov S.T.,Karasaev B.A.,Kandidat M., Analysis of manipulator grains for overloading largesized goods in the sewing industry. //Izvestiya Vysshikh Uchebnykh Zavedenii, Seriya Teknologiya Tekstil'noi Promyshlennosti. – 2018. – Vol.378. – P.332336 DOI: 
64. 
Numerical heat transfer simulation of solargeothermal hybrid source heat pump in Kazakhstan climates 
DOI: http://www.bcc.bas.bg/BCC_Volumes/Volume_50_Special_G_2018/50G_PD_713.56.pdf (SJR 0,18; Q4, Procentile 15 (General Chemistry)). 
A numerical energy balance model has been proposed in this work for predicting the thermal performances of solargeothermal hybrid source heat pumps in winter climates of Kazakhstan. The numerical simulation was performed for the year round continental climate conditions. The energy balance model has been developed based on first law of thermodynamics. The proposed heat pump configuration is working in solar mode during sunshine hours and geothermal model during off sunshine hours. Moreover, the system is operating in solargeothermal hybrid model to meet the evaporator load during insufficient availability of solar and geothermal sources. The energy performance comparison between conventional geothermal source and solarground hybrid source heat pump configurations has been made. The influences of solar intensity, ambient temperature, heat pump operating temperatures are discussed. © 2018 Bulgarian Academy of Sciences. 
Saktashova G.,Aliuly A.,Belyayev Ye.,Mohanraj M.,Singh R.M., Numerical heat transfer simulation of solargeothermal hybrid source heat pump in Kazakhstan climates. //Bulgarian Chemical Communications. – 2018. – Vol.50. – P.713 DOI: http://www.bcc.bas.bg/BCC_Volumes/Volume_50_Special_G_2018/50G_PD_713.56.pdf 
65. 
Identification of productive horizons of hydrocarbon deposits along spatial models of the structure and geomechanical state of the upper part of the earth's crust 
DOI: 0.5593/sgem2018/1.4/S06.062 (SJR 0,22; Q, Procentile 17(Earth and Planetary Sciences Geotechnical Engineering and Engineering Geology) 15(Earth and Planetary Sciences Geology) ). 
Modern ideas about the formation of hydrocarbon deposits are based on the results of experimental and theoretical studies of the lithosphere dynamics. They are based on the established facts of the "instability of fluid systems in the sedimentary cover and foundation, as well as the migration of fluids through the permeable zones resulting from the varying stresses in the earth's crust". It follows that the geological environment is non equilibrium. Forces of external influence and thermodynamic conditions inside the crust promote the extrusion of fluids from the lower horizons to the upper horizons and ensure their movement through the zones of decompaction to the places of discharge, which are often the nearfault zones limited by fluid. The developed technology makes it possible to identify horizons in the context of the geological environment with which hydrocarbon accumulations can be associated, visualize their position in space and, thus, create a basis for the design of the location of wells and solve several problems: a) Obtain a clear spatial representation of the properties of the geological environment. In particular, the provisions of reservoirs supplying hydrocarbon migration channels, detection of new zones of increased permeability; b) According to the distribution of the parameters of the stressstrain state, it is possible to obtain a spatial representation of the direction of migration of the fluid flows; c) The notion of the morphology of productive horizons provides a basis for designing wells and allows conclusions about the possible location of new deposits; d) The proposed methodology can be used as one of the stages of prospecting and exploration of hydrocarbon deposits on the projected areas in order to identify productive horizons in the crustal section and to assess the possible direction of fluid flows in the block under study;. © SGEM2018. 
Zhantayev Z.,Iskakov B.,Fremd A., Bibossinov A., Identification of productive horizons of hydrocarbon deposits along spatial models of the structure and geomechanical state of the upper part of the earth's crust. //International Multidisciplinary Scientific GeoConference Surveying Geology and Mining Ecology Management, SGEM. – 2018. – Vol.18. – P.467474 DOI: 0.5593/sgem2018/1.4/S06.062 
66. 
GPSderived velocity fields of northern tien shan from permanent stations in Kazakhstan 
DOI: 10.1109/IGARSS.2018.8518103 (SJR 0,26; Q, Procentile 34(Earth and Planetary Sciences General Earth and Planetary Sciences) 22(Computer Science Computer Science Applications)). 
This study presents and analyzed the data on GPS observations from 10 permanent stations located in the Northern Tien Shan, Kazakhstan. Analyzed data covers the 6 years, for the period of 20092015, which is acceptable for the geodynamic analysis. This investigation based on wellestablished processing method with GAMIT/GLOBK software package. The velocity field and motion rates of GPS sites were analyzed in the study. All stations showed northwest velocity field with the rate of 2 7 mm per year. Motion rates of GPS sites across the Northern Tien Shan decrease northwestward in gradual way. The speeds and directions of the vector fields in Southeastern Kazakhstan do not contrary with modern concepts of regional geodynamics and crustal deformation rates. © 2018 IEEE 
Zhantayev Z.,Kaldybayev A.,Bibossinov A., Vilyaev A.,Turgumbayev A.,Nurakynov S ,GPSderived velocity fields of northern tien shan from permanent stations in Kazakhstan. //International Geoscience and Remote Sensing Symposium (IGARSS). – 2018. – Vol.2018July. – P.31893191 DOI:10.1109/IGARSS.2018.8518103 
67. 
Determination of dissipative parameters of crank press 
DOI: 10.32014/2018.2518170X.40 (SJR 0,32; Q3, Procentile 40(Earth and Planetary SciencesGeology) 37(Earth and Planetary Sciences  Geotechnical Engineering and Engineering Geology) ). 
The paper considers the determination of the dissipation coefficients of the crank press's assemblies (units). Slidecrank mechanism of the press performs stamping of various parts. In such a case, there are significant dynamic loads in the assemblies and connections of the press. Work features of the crank press are connected with shock cyclic loads. Dynamic research of crank presses, at present, is a critical task. Dynamic model of the crank press is represented as an oscillatory system with many degrees of freedom. The dynamic model of the crank press consists of lumped masses connected by elasticdissipative elements. To calculate the dynamic model of the crank press, it is required to determine dissipation coefficients of individual sections of the drive shaft. Determination of the dissipation coefficients is of great complexity. In the case of vibrations of elastic systems, energy is dissipated in machine paits and mechanical energy is transferred into thermal. Losses of mechanical energy are caused by forces of inelastic resistance, called dissipative forces. The amplitudes of the resonance oscillations, the conditions for the excitation of parametric oscillations and selfoscillations depend on the dissipation coefficients. The determination of dissipative coefficients is based on available information on the dissipation coefficients of individual elastic dissipative elements of the system. An example is given for determining the coefficients of dissipation of the crank press. The dynamic model of the crank press was compiled on the program complex SimulationX. The obtained results show that the dynamics of oscillations of the drive shaft of the crank press changes significantly after adjustment of the coefficients of dissipation. © 2018, National Academy of Sciences of the Republic of Kazakhstan. All rights reserved. 
Tuleshov А.K.,Jomartov А.А., Determination of dissipative parameters of crank press. //News of the National Academy of Sciences of the Republic of Kazakhstan, Series of Geology and Technical Sciences. – 2018. – Vol.6. – P.102109 DOI: 10.32014/2018.2518170X.40 
68. 
Vector method for kinetostatic analysis of planar linkages 
DOI: 10.1007/s404300181022y (SJR 0,43; Q2, Procentile 74(Mathematics Applied Mathematics) 73(Engineering General Engineering) 71(Engineering Automotive Engineering) 68(Engineering Aerospace Engineering) 65(Engineering Industrial and Manufacturing Engineering) 64(Engineering Mechanical Engineering) ). 
The article reviews the vector method of solving one of the important problems of the dynamic analysis of planar linkages—the kinetostatic analysis, which is to determine the reactions in kinematic pairs and to balance the moment (force) on the driving medium with a given law of motion of the mechanism. This problem is of great practical value when choosing the electric motor and calculating bearings in kinematic pairs. Proposed quality criteria for an optimal kinematic synthesis of linkage mechanism depend only on the geometrical dimensions and are independent of the forces and moments applied to the linkages. © 2018, The Brazilian Society of Mechanical Sciences and Engineering. 
Jomartov A.,Tuleshov A.,Vector method for kinetostatic analysis of planar linkages. //Journal of the Brazilian Society of Mechanical Sciences and Engineering. – 2018. – Vol.40. DOI: 10.1007/s404300181022y 
69. 
The model of a transition region with smoothness conditions of the second order 
DOI: 10.1063/1.5049032 (SJR 0,18; Procentile 17(Physics and Astronomy General Physics and Astronomy)). 
We consider, a problem of determining the position of the conjugate points of a combined trajectory, ensuring the conditions of tangency, continuity, and equality of the radius of curvature. When the arcs of a circular and elliptical trajectory are conjugated at the junction where, there is a common tangent, a jump in the centrifugal force occurs due to the inequality of the radii of curvature, which will lead to a shock / 12 /. © 2018 Author(s). 
Bostanov B.O.,Temirbekov E.S.,,Matin D.T , The model of a transition region with smoothness conditions of the second ordert. //AIP Conference Proceedings. – 2018. – Vol.1997. DOI: 10.1063/1.5049032 
70. 
Noise resistance processing of speech signals 
DOI: 10.1109/EIConRus.2018.8317301 (SJR 0,87; Procentile 68). 
This article presents a new superior noise resistance method for processing speech signals. The proposed method of processing is based on adaptive sequential filtering. The presented results of processing speech signals demonstrate an increase in the quantitative parameters of noise resistance and reduction of errors in the detection of speech signals. The reliability of the obtained processing results is confirmed by cluster analysis. © 2018 IEEE. 
Zhumasheva Z.T..Altay Y.A.,Kremlev A.S, Noise resistance processing of speech signals. //Proceedings of the 2018 IEEE Conference of Russian Young Researchers in Electrical and Electronic Engineering, ElConRus 2018. – 2018.– P.11801185 DOI: 10.1109/EIConRus.2018.8317301 
71. 
Synthesis of Reconfigurable Positioning Parallel Manipulator of a Class RoboMech 
DOI: 10.1109/REMAR.2018.8449852 (SJR 1,95; Procentile 87). 
In this paper the methods of structural and dimensional synthesis of reconfigurable positioning parallel manipulator of a class RoboMech are presented. These parallel manipulators are formed by connecting the output point to the base by three closed kinematic chains: passive, active and negative. Parameters of synthesis of the active and negative closed kinematic chains are determined on the basis of the formulated Chebyshev and leastsquare approximations. © 2018 IEEE. 
Baigunchekov Z.,Ibrayev S.,Izmambetov M.,Zhumasheva Z.,Baigunchekov T., Naurushev B.,Mustafa A.,Sagyntay M. Synthesis of Reconfigurable Positioning Parallel Manipulator of a Class RoboMech. //2018 International Conference on Reconfigurable Mechanisms and Robots, ReMAR 2018  Proceedings. – 2018. DOI: 10.1109/REMAR.2018.8449852 
72. 
Thermal Performance Simulation of a Heat Pump Assisted Solar Desalination System for Kazakhstan Climatic Conditions 
DOI: 110.1080/01457632.2018.1451246 (SJR 0,18;Q4, Procentile 15 (General Chemistry)). 
A numerical model has been proposed in this research work for predicting the energy performance of a heat pump assisted regenerative solar still using phase change heat storage materials such as paraffin wax and hydrate salt. The numerical model is based on energy and mass balance. The performance was simulated for Fort Sherchenko town in Kazakhstan. The performance comparison between the conventional solar still and heat pump assisted regenerative solar still with and without phase change materials are presented. Moreover, the influence of two different phase change materials on thermal performance of a heat pump assisted regenerative solar still are described. The numerical simulation results showed that heat pump assisted regenerative solar still is more energy efficient and produces better productivity when compared to the conventional simple solar still. The use of hydrate salt phase change heat storage material was found to be a good option for enhancing the performance of a heat pump assisted regenerative solar still. © 2018, © 2018 Taylor & Francis Group, LLC. 
Belyayev, Y., Mohanraj, M., Jayaraj, S., Kaltayev, A., Thermal Performance Simulation of a Heat Pump Assisted Solar Desalination System for Kazakhstan Climatic Conditions. //Heat Transfer Engineerings. – 2019. – Vol.40. – P.10601072 DOI: 10.1080/01457632.2018.1451246 
73. 
Numerical simulation of a heat pump assisted solar dryer for continental climates 
DOI: 10.1016/j.renene.2019.04.119 (SJR 1,83;Q1, Procentile 88 (Renewable Energy, Sustainability and the Environment)). 
A numerical model has been proposed in this work for predicting the energy performance of a heat pump assisted solar dryer under continental climates. The model is based on energy and mass balance. The energy performance comparison between heat pump dryer, solar dryer and heat pump assisted solar dryer are presented. The simulation was performed for four different climatic conditions of Almaty city in Kazakhstan. The simulation results showed that, the heat pump assisted solar dryer is more energy efficient when compared to the conventional solar dryers. It is also confirmed that, the conventional solar dryers are not suitable for continental climatic conditions having low ambient temperatures. The heat pump dryer reduces the initial moisture content of banana (on wet basis) from about 74% to the final moisture content (on wet basis) of about 19% in 21 h. Similarly, the solar dryer reduces the initial moisture content (wet basis) from about 74% to the final moisture content (wet basis) of about 20% in 35 h. The specific moisture extraction rate and coefficient of performance of a heat pump assisted solar dryer are estimated to be about 0.6 kg/kWh and 2.72, respectively. © 2019 Elsevier Ltd 
Kuan M.,Shakir Y.,Mohanraj M.,Belyayev Y.,Jayaraj S.,Kaltayev A.,Numerical simulation of a heat pump assisted solar dryer for continental climates. //Renewable Energy. – 2019. – Vol.143. – P.214225 DOI: 10.1016/j.renene.2019.04.119 
74. 
Stabilization of the movement of a small spacecraft in a geostationary orbit 
DOI: 0.32014/2019.2518170X.162 (SJR 0,32;Q3 Procentile 40(Earth and Planetary SciencesGeology) 37(Earth and Planetary Sciences  Geotechnical Engineering and Engineering Geology)). 
A system for controlling the movement of a small spacecraft, which will later be part of a formation located in the geostationary orbit for the remote sensing of the Earth is being considered. The article assesses the effect of acting moments on the small spacecraft in the geostationary orbit and justifies the choice of these moments in the motion simulation. A motion control system for small spacecraft in the geostationary orbit using a Pcontroller (a controller with a proportional control law) is proposed. The results of the simulation to verify the correctness of the proposed control law are given. 
Rakisheva, Z.B., Nakasuka, S., Doszhan, N.S., Ibrayev, G.E.,Stabilization of the movement of a small spacecraft in a geostationary orbit. //News of the National Academy of Sciences of the Republic of Kazakhstan, Series of Geology and Technical Sciences. – 2019. – Vol.6. – P.112120 DOI:10.32014/2019.2518170X.162 
75. 
Optimization Issues in the Problem of Small Satellite Attitude Determination and Control 
DOI: 10.1007/9783030105013_14 (SJR 0,52;Q2 Procentile 17(Mathematics  Control and Optimization)). 
The problems of synthesis of attitude determination and control system of small satellites regarding the influence of external perturbations due to their small mass and restrictions in using highprecision actuators, due to the limitations of energy budget and construction, are considered. The solutions to this problem are proposed involving the development of highaccuracy algorithms for satellite attitude determination and control, using the minimum set of sensors, various types of actuators, and optimization principles. 
Rakisheva, Z., Sukhenko, A., Kaliyeva, N., Optimization Issues in the Problem of Small Satellite Attitude Determination and Control. //Springer Optimization and Its Applications,. – 2019. – Vol.144. – P.373393 
76. 
Spatial variations in the caspian sea wave climate in 2002–2013 from satellite altimetry 
DOI: 10.3176/earth.2019.16 (SJR 0,34;Q3, Procentile 44(Earth and Planetary Sciences General Earth and Planetary Sciences) 38(Environmental Science  Water Science and Technology)). 
The core properties of the wave climate and its changes in the Caspian Sea are established in terms of the annual mean significant wave height and its regional changes in 2002–2013 based on the outcome of the satellite altimetry mission JASON1. Remotely estimated wave heights are validated against properties of the empirical distribution of instrumentally measured wave heights in the southern Caspian Sea and monthly averages of visually observed wave heights at three locations. A correction for systematic differences leads to very good correspondence between monthly averaged in situ and satellite data with a typical root mean square difference of 0.06 m. The average significant wave height in the Caspian Sea is 0.5–0.7 m in the northern basin of the sea, around 1.2 m in large parts of the central and southern basins and reaches up to 1.8 m in the northern segment of the central basin. The basinwide average wave intensity varied insignificantly in the range of 1.02–1.14 m in 2002–2013. These estimates overestimate the wave heights by about 30% because low wave conditions are ignored. Substantial and statistically significant changes in the wave height occurred in certain areas. The wave height decreased by 0.019 ± 0.007 m/yr in the eastern segment of the central basin and by 0.04 ± 0.04 m/yr in the western segment of the southern basin. These changes can be explained by an increase in the frequency of westerly winds at the expence of southerly winds. Both basinwide and regional extreme wave heights exhibit large interannual variations but do not show any significant trend. The patterns of changes in mean and extreme wave height are different. The average wave height has increased while the extreme wave height has decreased in the eastern segment of the southern basin. © 2019 Authors. 
Kudryavtseva N.,Kussembayeva K.,Rakisheva Z.B.,Soomere T., Spatial variations in the caspian sea wave climate in 2002–2013 from satellite altimetry. //Estonian Journal of Earth Sciences. – 2019. – Vol.68. – P.225240 DOI: 10.3176/earth.2019.16 
77. 
Optimal structural synthesis of agricultural legged robot with minimal damage on soil 
DOI: 10.1051/e3sconf/201913501027 (SJR 0,2; Procentile 25(Energy General Energy) 20(Environmental Science General Environmental Science) 19(Earth and Planetary Sciences General Earth and Planetary Sciences)). 
Optimal structural synthesis of agricultural legged robot is carried out, that causes minimal damage on soil and provide the most favorable conditions for plant growth. A rational structure of a legged robot with orthogonal propel based on kinematic decoupling of the motion is proposed. Most of traditionally used walking robots have universal "insect type" structure with multiple actuators to be synchronized which result in complex control. The alternative design is realizing a concept of functional independence of the actuators when each actuator is responsible for specified purpose: the main actuators are responsible for rectilinear translational motion of cabine/chassis whereas another group of actuators are responsible for adaptation purposes and anothers participate in turning/maneuvering. This allows to carry out the cabine/hull shifting and turning with a minimum number of actuators and simplified control. A new kinematic equivalent scheme of turning mechanism was proposed in order to optimize turning modes of the robot. The proposed planar model allows to determine the optimal parameters of the robot by applying multicriteria synthesis methods of parallel robot manipulators based on isotropy, maneuverability and other criteria. © The Authors, published by EDP Sciences, 2019. 
Ibrayev S.,Jamalov N.,Ibrayeva A.,Mukhambetkaliyeva G., Optimal structural synthesis of agricultural legged robot with minimal damage on soil. //E3S Web of Conferences. – 2019. – Vol.135. DOI: 10.1051/e3sconf/201913501027 
78. 
[A comprehensive study to evaluate absolute uptake of carbon dioxide adsorption onto composite adsorbent] [Une étude complète visant à mesurer l'absorption absolue de dioxyde de carbone sur un adsorbant composite] 
DOI: 10.1016/j.ijrefrig.2019.01.014 (SJR 1.5; Q1 , Procentile 90(Engineering Building and Construction) 88(Engineering Mechanical Engineering)). 
In this study, new composite adsorbent with enhanced thermal conductivity and adsorption capacity was synthesized and analyzed comprehensively for the development of compact CO 2 based adsorption cooling system. The consolidated composite was prepared employing activated carbon, graphene nanoplatelets and hydroxyl cellulose as a parent adsorbent, thermal conductivity enhancer, and binder, respectively. The surface area and pore volume of the composite were found to be 1778 ± 13 m 2 g −1 and 1.014 cm 3 g −1 , respectively. In addition, the composite showed 233% higher thermal conductivity compared to the parent activated carbon. Adsorption characteristics of CO 2 were measured at temperature ranging from 20 to 70 °C and pressures up to 5 MPa. Absolute uptake was evaluated from excess adsorption based on the following two methods: (i) the adsorbed phase volume is equal to the pore volume of the adsorbent; and (ii) the adsorbed phase volume is almost zero under low pressure and/or high temperature conditions. Furthermore, the averaging of above two methods was also taken for avoiding these two extreme assumptions. Obtained absolute adsorption uptake data were fitted with modified DubininAstakhov and Tóth models. Results indicated good approximation between data points and models. The average isosteric heats of adsorption estimated using modified DA and Tóth model were found to be 19.742 kJ mol −1 and 19.023 kJ mol −1 , respectively. The obtained characteristics of composite adsorbent are prerequisites for designing compact CO 2 based adsorption cooling systems. © 2019 Elsevier Ltd and IIR

Berdenova B.,Pal A.,Muttakin M.,Mitra S.,Thu K.,Saha B.B.,Kaltayev A., A comprehensive study to evaluate absolute uptake of carbon dioxide adsorption onto composite adsorbent. //International Journal of Refrigerations. – 2019. – Vol.100. – P.131140 DOI: 10.1016/j.ijrefrig.2019.01.014 
79. 
Grasps of robot manipulator when overloading solid highradioactive elements and their calculation 
DOI: 10.1007/9783030033200_34 (SJR 0,16;Q4, Procentile\ 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
In this paper considers the problem of simulating two grasping schemes of the robot manipulator, created by the authors for the overload of solid highly radioactive waste of fuel elements. Two schemes, previously proposed by the authors as grasps, are investigated. One  rod, the second  innovative with flexible actuators. For the grasps considered, some results are given: Static calculations for flat grasping models; 3D models of grasp; the strength and rigidity of one 3D model are calculated. © Springer Nature Switzerland AG 2019 
Temirbekov E.S.,Kaiym T.T.,Ceccarelli M., Bostanov B.O.,Kaimov S.T.,Kaimov A.T. Grasps of robot manipulator when overloading solid highradioactive elements and their calculation. //Mechanisms and Machine Sciences. – 2019. – Vol.68. – P.316323 DOI: 10.1007/9783030033200_34 
80. 
Combined trajectory of continuous curvature 
DOI: 10.1007/9783030033200_2 (SJR 0,16;Q4, Procentile\ 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
The problem of creating a smooth combined trajectory and determining the position of the connection points ensuring the conditions of tangency, continuity and equality of the radius of curvature are considered. When the circular and conical arcs are conjoined at the junction point, there are nonsmoothness, there is a jump in the centrifugal force due to the inequality of the radii of curvature, which will lead to a strike. A method is used to insert a transitional section, the model of which is a conic, in order to exclude the undesirable impact effect. There are known methods of continuous, smooth connection of trajectories, which include contours, splines. The condition of equality of the radii of curvature is not fulfilled at such conjugations of these curves at the junction point. The mathematical condition of unstressed connection of trajectories is established. You can define the final connection point by specifying the starting point of the connection and using the condition of the unstressed connection. The proposed method allows to design complex technical forms and create on their basis new models of a smooth trajectory from conical arcs. © Springer Nature Switzerland AG 2019 
Temirbekov E.S.,Bostanov B.O.,Dudkin M.V.,Kaimov S.T.,Kaimov A.T., Combined trajectory of continuous curvature. //Mechanisms and Machine Sciences. – 2019. – Vol.68. – P.1219 DOI: 10.1007/9783030033200_2 
81. 
A gripper mechanism to automate overload process for fuel elements 
DOI: 10.1007/9783030003654_15 (SJR 0,16;Q4 Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials)). 
The paper reviews new and innovative ways to overload fuel elements of the intermediate container of solid high level radioactive waste into the cavity of the main container made from a strong block of rock, the implementation of which is planned to be performed with the use of innovative design gripper of the gripper of the manipulator of an industrial robot. Also reported main requirements for the functioning of a gripper and calculations. © 2019, Springer Nature Switzerland AG. 
Kaimov S.,Kaimov A.T.,Ceccarelli M., Kaiym T.,Kaimova G.,Jomartov A.,Temirbekov E., A gripper mechanism to automate overload process for fuel elements. //Mechanisms and Machine Science. – 2019. – Vol.66. – P.118128 DOI: 10.1007/9783030003654_15 
82. 
Computer modeling of residual stresses and strains at arc welding by modulated current 
DOI: 10.1007/9783319956305_265 (SJR 0,15.Q4; Procentile 17(Engineering Aerospace Engineering) 17(Engineering Mechanical Engineering) 15(Engineering Automotive Engineering) 13(Chemical Engineering Fluid Flow and Transfer Processes) ). 
The article reviews the investigations of forming residual stresses and strains at arc welding by the modulated current. The maintenance of a favorable thermal cycle in a zone (that is close to a joint) is very important at welding many metals and their fusions. By changing a thermal cycle it is possible to reduce a grain value, to achieve a minimal speed of increasing elastoplastic strains in a zone of high temperature, to control a cooling rate and thus to receive favorable structures in a zone near to a joint. The numerical modeling of thermodeformation processes of thin plates at arc welding by the modulated current is suggested in the paper. The complex of investigations depends on various technological parameters. The comparison of the received calculation data with experimental data has shown the acceptability of the positions and suppositions accepted in the design model. © 2019, Springer Nature Switzerland AG. 
Nurguzhin M.,Danenova G.,Akhmetzhanov T., Computer modeling of residual stresses and strains at arc welding by modulated current. //Lecture Notes in Mechanical Engineering. – 2019. – Vol.0. – P.24532460 DOI: 10.1007/9783319956305_265 
83. 
Optimization Issues in the Problem of Small Satellite Attitude Determination and Control 
DOI: 10.1007/9783030105013_14 (SJR 0,52;Q2, Procentile 17(Mathematics Control and Optimization)). 
The problems of synthesis of attitude determination and control system of small satellites regarding the influence of external perturbations due to their small mass and restrictions in using highprecision actuators, due to the limitations of energy budget and construction, are considered. The solutions to this problem are proposed involving the development of highaccuracy algorithms for satellite attitude determination and control, using the minimum set of sensors, various types of actuators, and optimization principles. © 2019, Springer Nature Switzerland AG. 
Rakisheva Z.,Sukhenko A.,Kaliyeva N., Optimization Issues in the Problem of Small Satellite Attitude Determination and ControlDetermination of dissipative parameters of crank press. //Springer Optimization and Its Applicationss. – 2019. – Vol.144. – P.373393 DOI: 10.1007/9783030105013_14 
84. 
Designing of a Crank press on the Basis of High Class Planar Linkages 
DOI: 10.1007/9783030201319_298 (SJR 0,16;Q4, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
Currently, crank presses based on planar linkages received prevalence for presswork. One of the major problems of these presses is skewness of slide caused by eccentric application of deforming force and total linear elastic deformation of the press links and stamps under load. The skewness of the slide of the crank press reduces the accuracy of punching and causes its jamming and lead to its breakage. To eliminate these problems, it is proposed to use highclass planar linkages to design crank presses in the work. Schemes of planar linkages of class IV for designing of new crank presses were obtained. A new crank press was synthesized on the basis of planar linkages of class IV and its prototype was made. The tests of the prototype of the new crank press were carried out. © 2019, Springer Nature Switzerland AG. 
Tuleshov А.K.,Jomartov А.А., Kuatova M.,Designing of a Crank press on the Basis of High Class Planar Linkagess. //Mechanisms and Machine Science. – 2019. – Vol.73. – P.30273036 DOI: 10.1007/9783030201319_298 
85. 
Modeling dynamics of planetary gears of crank press on simulationx 
DOI: 10.1007/9783030033200_5 (SJR 0,16;Q4, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
The paper investigates the dynamics of a planetary gear of crank press. Currently, dynamic research of planetary gear of crank presses with account for elasticity of links, friction and backlashes, is a critical task. Dynamics of planetary gear of crank presses is well simulated on a software package SimulationX, with the help of universal basic structures. The paper shows that the dynamics of planetary mechanism with account for backlashes essentially differs from the dynamics of the ideal mechanism. The maximum acceleration of the carrier, when taking into account the backlashes in bearings of gears of the planetary gears of the crank press, increases by 30%. © Springer Nature Switzerland AG 2019 
Tuleshov А.K.,Jomartov А.А., Modeling dynamics of planetary gears of crank press on simulationx. //Mechanisms and Machine Science. – 2019. – Vol.68. – P.4148 DOI: 10.1007/9783030033200_5

86. 
Analysis of the oscillating motion of a solid body on vibrating bearers 
DOI: 10.3390/MACHINES7030058 (SJR 0,39;Q2 Procentile 83(Computer Science Computer Science (miscellaneous)) 82(Mathematics Control and Optimization) 77(Engineering Mechanical Engineering) 76(Engineering Industrial and Manufacturing Engineering) 72(Engineering Electrical and Electronic Engineering) 71(Engineering Control and Systems Engineering) ). 
This article considers the oscillation of a solid body on kinematic foundations, the main elements of which are rolling bearers bounded by highorder surfaces of rotation at horizontal displacement of the foundation. Equations of motion of the vibroprotected body have been obtained. It is ascertained that the obtained equations of motion are highly nonlinear differential equations. Stationary and transitional modes of the oscillatory process of the system have been investigated. It is determined that several stationary regimes of the oscillatory process exist. Equations of motion have been investigated also by quantitative methods. In this paper the cumulative curves in the phase plane are plotted, a qualitative analysis for singular points and a study of them for stability are performed. In the Hayashi plane a cumulative curve of a body protected against vibration forms a closed path which does not tend to the stability of a singular point. This means that the vibration amplitude of a body protected against vibration does not remain constant in a steady state, but changes periodically. © 2019 by the authors. 
Bissembayev K.,Jomartov A.,Tuleshov A.,Dikambay T., Analysis of the oscillating motion of a solid body on vibrating bearers. //Machines. – 2019. – Vol.6. – №58 DOI: 10.3390/MACHINES7030058 
87. 
Dynamic analyses of a clutch of crank press 
DOI: 10.32014/2019.2518170X.132 (SJR 0,32;Q3, Procentile 40(Earth and Planetary SciencesGeology) 37(Earth and Planetary Sciences  Geotechnical Engineering and Engineering Geology)). 
The paper studies the dynamic of a clutch of crank press. At present, the dynamic study of clutch of the crank presses, with account of interaction with other blocks, is a priority. The crank press contains movable parts and assemblies, the mass of which is from one hundred kilograms to several tons. These parts and assemblies are connected cyclically by clutch of crank press with high speeds and they are subject to large dyjnamic loads. To simu late and analyze the movement of crank press with clutch, a software package: SimulationX is used. SimulationX is a software package for modeling and analyzing the dynamics and kinematics of cars, industrial equipment, electric, pneumatic and hydraulic drives, hybrid engines, etc. As a result of dynamic calculation, important dynamic para meters of the crank press clutch and working slide are determined. It is shown that dynamic loads sharply increase almost in all blocks of the crank press when the clutch is switched on. © 2019, National Academy of Sciences of the Republic of Kazakhstan. All rights reserved. 
Tuleshov A.K.,Jomartov A.A.,Kuatova M.J., Halicioglu R., Dynamic analyses of a clutch of crank press. //News of the National Academy of Sciences of the Republic of Kazakhstan, Series of Geology and Technical Sciences. – 2019. – Vol.5. – P.114123 DOI: 10.32014/2019.2518170X.132 
88. 
Creating the form of a complicated treadmill 
DOI: 10.1007/9783030201319_38 (SJR 0,16; Q4 , Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials)). 
The problem of creating a smooth complex trajectory and determining the position of the points of connection, providing the conditions of tangency, continuity and equality of the radius of curvature, is considered. When conjugation trajectory in the form of arcs of curves, there are nonsmoothness at the junction, a centrifugal force jump occurs due to the inequality of the radii of curvature, which will lead to the impact. Known methods and methods for continuous, smooth connection of trajectories, which include lines, splines, as well as specially designed transition curves, such as a majorant curve, clothoid, elastic line, lemniscate and velocity curve. However, the proposed methods are less accurate and approximate. With such conjugations, the condition of equality of the radii of curvature of these curves at the junction point is not satisfied. To eliminate the undesirable effect of impact, the method of inserting a transition section is used, the model of which is a secondorder curve (conic). An original method is proposed for analytically determining a smooth (smooth second order) transition region in the form of a conic. A mathematical condition has been established that ensures a smooth connection without a jump in the radius of curvature. By setting the start point of connection and using mathematical condition for smooth connection can determine the finish point of connection. The process of determining the smoothness of conjugation is proposed to be modeled by a link mechanism. The proposed method makes it possible to construct complex technical forms and to create on their basis new models of a smooth trajectory from arcs of highorder curves with a conical transition section. © 2019, Springer Nature Switzerland AG. 
Temirbekov Y.S.,Bostanov B.O.,Karassayev B.A.,Tukeshova G.A Creating the form of a complicated treadmill. //Mechanisms and Machine Science. – 2019. – Vol.73. – P.377386 DOI: 10.1007/9783030201319_38 
89. 
Mechanicsmathematical model of conjugation of a part of a trajectory with conditions of continuity, touch and smoothness 
DOI: 10.1007/9783030122034_8 (SJR 0,16; Q4 , Procentile 30(Mathematics General Mathematics) 18(Computer Science General Computer Science)). 
Development process of combined trajectories, in places of joining conic arcs, undesirable intermittent effects inevitably arise due to the secondorder nonsmoothness. A second order tangency is considered taking into account the curvature and the equality condition of the arcs curvature radii to be joined at the conjugation points. A kinematic method for determining joints on the basis of a rocker mechanism is given, which ensures smooth joints. © Springer Nature Switzerland AG 2019. 
Bostanov B.O.,Temirbekov E.S., Dudkin M.V.,Kim A.I., Mechanicsmathematical model of conjugation of a part of a trajectory with conditions of continuity, touch and smoothness. //Communications in Computer and Information Science. – 2019. – Vol.998. – P.7181 DOI: 10.1007/9783030122034_8 
90. 
Grasps of robot manipulator when overloading solid highradioactive elements and their calculation 
DOI: 10.1007/9783030033200_34 (SJR 0,16;Q4, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
In this paper considers the problem of simulating two grasping schemes of the robot manipulator, created by the authors for the overload of solid highly radioactive waste of fuel elements. Two schemes, previously proposed by the authors as grasps, are investigated. One  rod, the second  innovative with flexible actuators. For the grasps considered, some results are given: Static calculations for flat grasping models; 3D models of grasp; the strength and rigidity of one 3D model are calculated. © Springer Nature Switzerland AG 2019 
Temirbekov E.S.,Kaiym T.T.,Ceccarelli M., Bostanov B.O.,Kaimov S.T.,Kaimov A.T.,Grasps of robot manipulator when overloading solid highradioactive elements and their calculation. //Mechanisms and Machine Science. – 2019. – Vol.68. – P.316323 DOI: 10.1007/9783030033200_34 
91. 
Combined trajectory of continuous curvature 
DOI: 10.1007/9783030033200_2 (SJR 0,16;Q4, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
The problem of creating a smooth combined trajectory and determining the position of the connection points ensuring the conditions of tangency, continuity and equality of the radius of curvature are considered. When the circular and conical arcs are conjoined at the junction point, there are nonsmoothness, there is a jump in the centrifugal force due to the inequality of the radii of curvature, which will lead to a strike. A method is used to insert a transitional section, the model of which is a conic, in order to exclude the undesirable impact effect. There are known methods of continuous, smooth connection of trajectories, which include contours, splines. The condition of equality of the radii of curvature is not fulfilled at such conjugations of these curves at the junction point. The mathematical condition of unstressed connection of trajectories is established. You can define the final connection point by specifying the starting point of the connection and using the condition of the unstressed connection. The proposed method allows to design complex technical forms and create on their basis new models of a smooth trajectory from conical arcs. © Springer Nature Switzerland AG 2019 
Temirbekov E.S.,Bostanov B.O.,Dudkin M.V.,Kaimov S.T.,Kaimov A.T., Combined trajectory of continuous curvature. //Mechanisms and Machine Science. – 2019. – Vol.68. – P.1219 DOI: 10.1007/9783030033200_2 
92. 
A gripper mechanism to automate overload process for fuel elements 
DOI: 10.1007/9783030003654_15 (SJR 0,16;Q4, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
The paper reviews new and innovative ways to overload fuel elements of the intermediate container of solid high level radioactive waste into the cavity of the main container made from a strong block of rock, the implementation of which is planned to be performed with the use of innovative design gripper of the gripper of the manipulator of an industrial robot. Also reported main requirements for the functioning of a gripper and calculations. © 2019, Springer Nature Switzerland AG. 
Kaimov S.,Kaimov A.T.,Ceccarelli M.Kaiym T.,Kaimova G.,Jomartov A.,Temirbekov E., A gripper mechanism to automate overload process for fuel elements. //Mechanisms and Machine Science. – 2019. – Vol.66. – P.118128 DOI: 10.1007/9783030003654_15 
93. 
Parallel manipulator of a class RoboMech for generation of horizontal trajectories family 
DOI: 10.1007/9783030201319_137 (SJR 0,16;Q4, Procentile 23 (Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials)). 
In this paper the methods of structural and dimensional synthesis of parallel manipulator of a class RoboMech for reproducing the series of horizontal trajectories are developed. This parallel manipulator is formed by connecting the output point to the base using two closing kinematic chains: one passive and one active ones. The synthesis parameters of the active closing kinematic chain of RPR type are determined on the basis of the leastsquare approximation. © 2019, Springer Nature Switzerland AG. 
Baigunchekov Z.,Izmambetov M.,Zhumasheva Z.,Baigunchekov T.,Mustafa A.,Parallel manipulator of a class RoboMech for generation of horizontal trajectories family. //Mechanisms and Machine Science. – 2019. – Vol.73. – P.13951402 DOI: 10.1007/9783030201319_137

94. 
Numerical simulation on solar collector and cascade heat pump combi water heating systems in Kazakhstan climates 
DOI: https://doi.org/10.1016/j.renene.2019.06.102 (SJR 0,83;Q1, Procentile 88 (Renewable Energy, Sustainability and the Environment) ). 
At low ambient temperatures, the heating capacity and coefficient of performance of a single stage vapour compression heat pump cycle is significantly getting reduced. The two stage cascade heat pump cycle operating with two different refrigerants provides a sustainable solution to lift the condenser temperature above 343 K. In this work, a numerical simulation model was developed for predicting the performance of a solar collector and two stage cascade heat pump combi water heating systems under Kazakhstan climatic conditions. The numerical simulation was performed for winter climatic conditions using nine refrigerant pairs such as, R32/R290, R32/R1234yf, R32/R134a, R410A/R290, R410A/R1234yf, R410A/R134a, R744/R290, R744/R1234yf and R744/R134a. The influences of solar irradiation, ambient temperature and condenser temperature are discussed. The solar collector and two stage heat pump combi heating system has about 30% energy savings when compared to the conventional two stage cascade heat pump without integration solar collectors. The R32/R290 refrigerant pair is having maximum coefficient of performance of 2.4 at 323 K condensing temperature and 263 K evaporating temperature. The refrigerant R744/R290 pair is identified as an environment friendly sustainable option in terms of its global warming impact for two stage cascade heat pump applications. 
Yerdesh, Y., Abdulina, Z., Aliuly, A..Mohanraj, M., Kaltayev, A, Numerical simulation on solar collector and cascade heat pump combi water heating systems in Kazakhstan climates. //Renewable Energy. – 2020. – Vol.145. – P.12221234 
95. 
A Study of Secular Perturbations of TranslationalRotational Motion in a Nonstationary TwoBody Problem Using Computer Algebra 
DOI: 10.1134/S0965542520010054 (SJR 0,53;Q2, Procentile 35(Mathematics  Computational Mathematics) ). 
A nonstationary twobody problem is considered such that one of the bodies has a spherically symmetric density distribution and is central, while the other one is a satellite with axisymmetric dynamical structure, shape, and variable oblateness. Newton’s interaction force is characterized by an approximate expression of the force function up to the second harmonic. The body masses vary isotropically at different rates. Equations of motion of the satellite in a relative system of coordinates are derived. The problem is studied by the methods of perturbation theory. Equations of secular perturbations of the translationalrotational motion of the satellite in analogues of Delaunay–Andoyer osculating elements are deduced. All necessary symbolic computations are performed using the Wolfram Mathematica computer algebra system. © 2020, Pleiades Publishing, Ltd. 
Bizhanova, S.B., Minglibayev, M.Z., Prokopenya, A.N.,A Study of Secular Perturbations of TranslationalRotational Motion in a Nonstationary TwoBody Problem Using Computer Algebra. //Mechanisms and Machine ScienceComputational Mathematics and Mathematical Physics,. – 2020. – Vol.60. – P.2635 DOI: 10.1134/S0965542520010054 
96. 
Numerical calculation of the pressure drop and saturation of twophase flow through porous medium 
DOI: https://doi.org/10.2118/202570MS (SJR 0, Procentile ). 
Twophase flow through a fibrous porous medium is numerically simulated. The goal of the study is to investigate the effect of numerical slip length on variation of pressure drop and saturation. The simulation is based on the numerical solution of incompressible NavierStokes equations for twophase flow around a circular cylinder with radii in the ranges from 0.1 to 0.4. Variations of pressure drop and saturation in time are numerically calculated on simple and snap meshes for different values of surface tension (σ), viscosity ratio (M = μnw/μw), and mesh sizes. Also, numerical results using simple and snap meshes are compared. Numerical calculations are performed using the interFoam solver in OpenFOAM® finite volume library. The presented results show that solutions of saturation and pressure drop in time vary with mesh refinement, so the numerical slip length significantly affects the flow. For a coarser mesh the influence of numerical slip length is less than for a finer mesh, also as bigger the value of the viscosity ratio or surface tension as lower the influence of numerical slip length. The results show that the contact line motion significantly depends not only on the microstructure of the porous media, but also on the fluid properties (viscosity ratio, surface tension and etc.). Copyright 2020, Society of Petroleum Engineers. 
Akasheva, Z., Assilbekov, B., Kudaikulov, A., Bolysbek, D.,Numerical calculation of the pressure drop and saturation of twophase flow through porous medium. //Society of Petroleum Engineers  SPE Annual Caspian Technical Conference 2020, CTC 2020. – 2020. SPE Annual Caspian Technical Conference 2020, CTC 2020 
97. 
Experimental and Theoretical Studies of the Efficiency of Autonomous Multistory Wind Power Plants 
DOI: 10.1134/S1063784220010168 (SJR 0,3;Q3, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials)23(Physics and Astronomy  Physics and Astronomy (miscellaneous)) ). 
The compact multistory wind power plants developed by authors for the first time in the world have been studied. Their efficiency due to the autonomy, compactness, and use of draft effect, which occurs between stories, has been shown. Thanks to this, compact multistory wind power plants will have some advantages in comparison with conventional wind power plants and separated winddriven generators with the same power. © 2020, Pleiades Publishing, Ltd. 
Kunakbaev, T., Tanasheva, N.K., Dyusembaeva, A.N., Shaimerdenova, K.M., Sagitzhanov, B.M.,Experimental and Theoretical Studies of the Efficiency of Autonomous Multistory Wind Power Plants. //Technical Physics. – 2020. – Vol.65. – P.3740 DOI: 10.1134/S1063784220010168 
98. 
Computeraided creation of coordinate grid for blasted rock block 
DOI: 10.25018/023614932020804051 (SJR 0,33;Q2, Procentile 24(Engineering Industrial and Manufacturing Engineering) 21(Environmental Science Ecology) 19(Earth and Planetary Sciences Geotechnical Engineering and Engineering Geology) 17(Earth and Planetary Sciences Geology) 9(Earth and Planetary Sciences Geochemistry and Petrology) ). 
The notions of ‘coordinate grid of blast block’ and ‘coordinate grid of blasted block’ are defined. The computeraided determination of the coordinate grid points for a blasted rock block is described. The starting point is using relationships of the initial parameters of a blast rock block such as height and the final pattern of the determined points in the blasted rock block. The governing factors of a blasted block coordinate grid are the brokenrock disintegration height and width dependent on physical and mechanical properties of rocks, physicochemical characteristics of explosives and on blasting pattern design. Considerable influence on arrangement of the blasted block coordinate grid points is exerted by the degree of fragmentation of rocks in horizontal and inclined layers. It is found that the degree of fragmentation tends to increase from the rear to the periphery of the broken rock pile, both horizontally and vertically. Based on the integration of all determinants of a broken rock pile after blasting, the analytical method is newly developed to determine the points and create the coordinate grid of a blasted block. The new computer program for the automatic creation of a blasted block coordinate grid takes into account the influences exerted by physical and mechanical properties of rocks, physicochemical characteristics of explosives and blasting pattern design. The program allows fast and sufficiently accurate position fixing of a blasted block coordinate grid points at different parameters and conditions of blasting. The resultant coordinate grids of blasted blocks are the tool for locating different elements inside brokenrock disintegration. © 2020 Publishing house Mining book. All rights reserved. 
Rakishev, B.R., Rakisheva, Z.B., Orynbay, A.A,Computeraided creation of coordinate grid for blasted rock block. //Mining Informational and Analytical Bulletin. – 2019. – Vol.2020(8). – P.4051 DOI: 10.25018/023614932020804051 
99. 
Automated determination of internal points of the coordinate grid of the blasted rock mass 
DOI: 10.1051/e3sconf/202016800015 (SJR 0,2; Procentile 25(Energy General Energy) 20(Environmental Science General Environmental Science) 19(Earth and Planetary Sciences General Earth and Planetary Sciences) ). 
An automated method for determining the internal points of the coordinate grid of the blasted rock mass is described. It is based on the method of determining the nodal points of the coordinate grid of the blasted rock mass, which is based on taking into account the dependencies that connect the initial parameters of the blasting rock mass with the final location of the fixed points of the blasted rock mass. The determining factors are the specific height and width of the collapse, the coefficient of loosening of the rocks. The method of analysis of experimental and industrial mass explosions in quarries, methods of analytical geometry, matrix theory and linear algebra are used. For the first time in mining, an analytical method has been developed for determining the internal points of the coordinate grid of an exploded block. It includes the established functions of the movement of nodal points, components of the vectors of movement of nodal and internal points of the coordinate grid. The established dependencies allow one to determine the displacements of any point inside the coordinate grid of the blasted block from the initial coordinates of the nodal and internal points. © The Authors, published by EDP Sciences, 2020. 
Rakishev, B., Rakisheva, Z., Auezova, A., Orynbay,, Automated determination of internal points of the coordinate grid of the blasted rock mass. //E3S Web of Conferences DOI: 10.1051/e3sconf/202016800015 
100. 
Digital hierarchical model of lumpiness of blasted rock mass 
DOI: 10.1080/25726668.2020.1838775 (SJR 0,56;Q2, Procentile 54(Earth and Planetary Sciences Geology) 48(Earth and Planetary Sciences Geotechnical Engineering and Engineering Geology) ). 
The article describes the developed new digital hierarchical model of lumpiness of blasted rocks. It is based on a combination of classes of pieces of rock by size, taken as hierarchical levels of lumpiness of rocks, and the percentage of pieces in classes, adopted for its numerical characteristics. 7 classes of rocks by size are considered, each of which is characterized by the content of its pieces. The key component of the proposed model  the granulometric composition of the blasted rocks is determined by the block mass of the rocks, the size of the zones of intense crushing, the volume of crushed rocks due to the action of stress waves and reflected waves (I stage of explosion), detonation products (swelling effect of the explosion) (II stage of explosion), and collisions of large pieces when moving (III stage of explosion). The regularities of the formation of these explosion results are established depending on the different combination of the physicomechanical properties of the rocks of the massif, the chemophysical characteristics of the explosive used, the parameters of the explosives and the explosive method. © 2020 Institute of Materials, Minerals and Mining and The AusIMM Published by Taylor & Francis on behalf of the Institute and The AusIMM. 
Rakishev, B., Rakisheva, Z.B., Auezova, A.M., Orynbay, A.A., Digital hierarchical model of lumpiness of blasted rock mass. //Mining Technology: Transactions of the Institute of Mining and Metallurgy. – 2020. – Vol.128(4).P.228237 DOI: 10.1080/25726668.2020.1838775 
101. 
The Effect of Applied Pressure Function on Thermoelastic Problem in the Dry Friction Clutches 
DOI: 10.1007/s11668020010314 (SJR 0,25;Q3, Procentile 36(Engineering Safety, Risk, Reliability and Quality) 34(Engineering Mechanical Engineering) 29(Materials Science General Materials Science) 29(Engineering Mechanics of Materials)). 
The main purpose of this paper is to investigate deeply the effect the contact pressure function on the thermomechanical behavior of the friction clutch system during the slipping time (heating stage). The other purpose is to explore theoretically the complex interaction among the contact pressure, sliding speed and frictional characteristics of frictional facings to specify the magnitude and distribution of the frictional heating generation on the contact surfaces of the dry friction clutch under different applied pressure. It was developed a numerical code based on finite element method (ANSYS/APDL 2019) to determine accurately the contact pressure, temperature and frictional heat generated on contact surfaces of the friction clutch disc that has two effective frictional faces at any instant during the slipping period. It was found a significant effect of the magnitude and variation of applied pressure during the heating phase on the surface temperatures, contact pressure and frictional heat generated. Where, the hot spot can be appeared when the applied pressure is constant. Under such circumstances, high amount of temperature and contact pressure focused on a small zone of the nominal contact area. This phenomenon is considered one of the main reasons for the early failure of the contacting surfaces of friction clutch. © 2020, ASM International. 
Stojanovic, N., Abdullah, O.I., Rakisheva, Z.B., Lattieff, F.A., Hashim, E.T., The Effect of Applied Pressure Function on Thermoelastic Problem in the Dry Friction Clutches. //Journal of Failure Analysis and Prevention. – 2020. – Vol.20(6). – P.21452152 DOI: 10.1007/s11668020010314 
102. 
Experimental and Theoretical Studies of the Efficiency of Autonomous Multistory Wind Power Plants 
DOI: 10.1134/S1063784220010168 (SJR 0,3;Q3, Procentile 23(Physics and Astronomy Physics and Astronomy (miscellaneous)) ). 
The compact multistory wind power plants developed by authors for the first time in the world have been studied. Their efficiency due to the autonomy, compactness, and use of draft effect, which occurs between stories, has been shown. Thanks to this, compact multistory wind power plants will have some advantages in comparison with conventional wind power plants and separated winddriven generators with the same power. © 2020, Pleiades Publishing, Ltd. 
Kunakbaev, T., Tanasheva, N.K., Dyusembaeva, A.N., Shaimerdenova, K.M., Sagitzhanov, B.M.,Experimental and Theoretical Studies of the Efficiency of Autonomous Multistory Wind Power Plants. //Technical Physics. – 2020. – Vol.65. – P.3740 DOI: 10.1134/S1063784220010168 
103. 
Synthesis of fourbar linkage with adjustable crank length for multipath generation 
DOI:10.18178/ijmerr.9.4.489495 (SJR 0,19;Q4, Procentile 26(Engineering Mechanical Engineering) 22(Engineering Control and Systems Engineering) 14(Computer Science Artificial Intelligence) ). 
Synthesis of planar mechanism with adjustable crank length for generating multiple paths is presented. Leastsquare approximation problem is considered which allows carrying out approximate synthesis with unlimited number of desired coupler point positions and with unlimited number of prescribed trajectories. By reducing the task to synthesis of twoelement link with variable binary link length, which is called RPRmodule, the analytical solution is obtained to determine not only constant design parameters (mechanism link lengths) but the adjusting parameter values as well. Thus the number of design variables for nonlinear optimization (applied to find the remaining parameters) will be decreased significantly. The applied method is exemplified by synthesis of the mechanism for variable straight line generation, where the required height of the endeffector is adjusted by adjusting the crank length. Combined with random search technique the method allows to find all local minimums of the optimized goal function and thus allows to take full advantage from the considered mechanism structure during design. © 2020 by the authors. 
Ibrayev S.,Jomartov A.,Tuleshov A.,Jamalov N.,Ibrayev A.,Mukhambetkalieva G.,Aidasheva G.,Kamal A.,Synthesis of fourbar linkage with adjustable crank length for multipath generation. //International Journal of Mechanical Engineering and Robotics Research. – 2020. – Vol.9(4). – P.489495 DOI: 10.18178/ijmerr.9.4.489495 
104. 
Designing of the stephenson ii sixlink linkage actuator for servo mechanical press 
DOI: 10.24247/ijmperdapr202053 (SJR 0; Procentile 22(Engineering Mechanical Engineering) 22(Engineering Aerospace Engineering) 17(Chemical Engineerin Fluid Flow and Transfer Processes) ). 
Currently, servo mechanical presses with an actuator based on a slidercrank mechanism are widely used for stamping details. One of the major problems of these servo mechanical presses is inclination of slide caused by eccentric application of deforming force and total linear elastic deformation of the press links and stamps under load. The inclination of the slide of the servo mechanical press reduces the accuracy of punching and causes its jamming and lead to its breakage. To eliminate these problems, it is proposed to use Stephenson II sixbar linkage as an actuator of servo mechanical press. The scheme structure of the actuator of servo mechanical press with two connecting rods and one crank based on the Stephenson II sixbar linkage is obtained. Synthesis of the Stephenson II sixbar linkage was carried out. As a result of synthesis, a new actuator of servo mechanical press on the base of the Stephenson II sixbar linkage was obtained and its prototype was made. The prototype of the new actuator of servo mechanical press based on the Stephenson II sixbar linkage was tested. Testing of the prototype of the new actuator of servo mechanical press showed a good distribution of the efforts applied during the presswork and a better tolerance of the eccentric load. © TJPRC Pvt. Ltd. 
Jomartov, A., Tuleshov, A., Jamalov, N., Kuatova, K., Kaimov, A., Designing of the stephenson ii sixlink linkage actuator for servo mechanical press. //International Journal of Mechanical and Production Engineering Research and Development. – 2020. – Vol.10(2). – P.501512 DOI: 10.24247/ijmperdapr202053 
105. 
Optimal synthesis of planar linkages 
DOI: 10.32014/2020.2518170X.21 (SJR 0,32;Q3, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials)40(Earth and Planetary Sciences Geology) 37(Earth and Planetary Sciences Geotechnical Engineering and Engineering Geology) ). 
This paper investigates the optimal synthesis of planar linkages. The main idea of this paper is to find the initial approximations based on the use of Burmester points for function generator linkages, path generator linkages, motion generator linkages. The results of the numerical synthesis of the linkages depend on the choice of the initial approximations. A more flexible method to the search for initial approximations is the method based on the use of Burmester points. This method allows the determination of the initial approximations analytically for three, four or five by established initial data of synthesis. In this case, the problem is reduced to determining the solutions of polynomials, respectively the second, third and fourth degree. The method consists in that the synthesized linkage is conditionally divided into initial kinematic chains and closing kinematic chains, and Burmester points are determined for each chain. After the choice of initial approximations, an objective function is formed according to the output criteria, depending on the synthesis parameters, using the Chebyshevsky (best) or quadratic approximation problems. The synthesis parameters of planar linkages are determined from objective function minimum. According to this method, a program for the synthesis of planar linkages has been developed. An example is included to demonstrate the method. © 2020, National Academy of Sciences of the Republic of Kazakhstan. All rights reserved. 
Tuleshov, А.K., Jomartov, А.А., Ibrayev, S., Jamalov, N.K., Halicioglu, R., A gripper mechanism to automate overload process for fuel elements. //News of the National Academy of Sciences of the Republic of Kazakhstan, Series of Geology and Technical Sciences. – 2020. – Vol.1(439). – P.172180 DOI: 10.32014/2020.2518170X.21 
106. 
Synthesis of fourbar linkage with adjustable crank length for multipath generation 
DOI: 10.18178/ijmerr.9.4.489495 (SJR 0,19;Q4, Procentile 26(Engineering Mechanical Engineering) 22(Engineering Control and Systems Engineering) 14(Computer Science Artificial Intelligence) ). 
Synthesis of planar mechanism with adjustable crank length for generating multiple paths is presented. Leastsquare approximation problem is considered which allows carrying out approximate synthesis with unlimited number of desired coupler point positions and with unlimited number of prescribed trajectories. By reducing the task to synthesis of twoelement link with variable binary link length, which is called RPRmodule, the analytical solution is obtained to determine not only constant design parameters (mechanism link lengths) but the adjusting parameter values as well. Thus the number of design variables for nonlinear optimization (applied to find the remaining parameters) will be decreased significantly. The applied method is exemplified by synthesis of the mechanism for variable straight line generation, where the required height of the endeffector is adjusted by adjusting the crank length. Combined with random search technique the method allows to find all local minimums of the optimized goal function and thus allows to take full advantage from the considered mechanism structure during design. © 2020 by the authors. 
Ibrayev S.,Jomartov A.,Tuleshov A.,Jamalov N.,Ibrayev A.,Mukhambetkalieva G.,Aidasheva G.,Kamal A.,Synthesis of fourbar linkage with adjustable crank length for multipath generation. //International Journal of Mechanical Engineering and Robotics Research. – 2020. – Vol.9(4). – P.489495 DOI: 10.18178/ijmerr.9.4.489495 
107. 
GNSSBased Attitude Determination TechniquesA Comprehensive Literature Survey 
DOI: 10.1109/ACCESS.2020.2970083 (SJR 0,59;Q1, Procentile 87(Engineering General Engineering) 81(Computer Science General Computer Science) 70(Materials Science General Materials Science) ). 
GNSSbased Attitude Determination (AD) of a mobile object using the readings of the Global Navigation Satellite Systems (GNSS) is an active area of research. Numerous attitude determination methods have been developed lately by making use of various sensors. However, the last two decades have witnessed an accelerated growth in research related to GNSSbased navigational equipment as a reliable and competitive device for determining the attitude of any outdoor moving object using data demodulated from GNSS signals. Because of constantly increasing number of GNSSbased AD methods, algorithms, and techniques, introduced in scientific papers worldwide, the problem of choosing an appropriate approach, that is optimal for the given application, operational environment, and limited financial funding becomes quite a challenging task. The work presents an extensive literature survey of the methods mentioned above which are classified in many different categories. The main aim of this survey is to help researchers and developers in the field of GNSS applications to understand pros and cons of the current state of the art methods and their computational efficiency, the scope of use and accuracy of the angular determination. © 2020 IEEE. 
Raskaliyev, A., Patel, S.H., Sobh, T.M., Ibrayev, A., GNSSBased Attitude Determination TechniquesA Comprehensive Literature Survey. //IEEE Access. – 2020. – Vol.8. – P.24873–24886, 8972427 DOI: 10.1109/ACCESS.2020.2970083 
108. 
Some approaches to assessing the quality of masking noise interference of spatial noise generators 
DOI: (SJR 0,15;Q4, Procentile 36(Computer Science General Computer Science) 16(Mathematics Theoretical Computer Science) ). 
The article discusses the characteristics of spatial electromagnetic noise generators and the formation of a broadband noise signal. It also describes a number of known methods and methods for assessing the quality of masking noise interference and their differences. Different approaches to measuring masking noise when evaluating its quality are proposed. The first method is based on the measurement of the instantaneous values of the amplitudes of the noise signal and the calculation of the entropy coefficient based on this method. The second method involves searching for correlation of masking noise signals of noise generators in different frequency subbands. The third approach is to use statistical and (or) graphical methods (tests) for randomness. The completeness and objectivity of assessing the quality of masking noise interference from spatial noise generators will be achieved by combining all the methods. © 2005 – ongoing JATIT & LLS. 
Smailov, N., Batyrgaliyev, A., Seilova, N., Kuttybaeva, A., Ibrayev, A., Some approaches to assessing the quality of masking noise interference of spatial noise generators. //Journal of Theoretical and Applied Information Technology. – 2020. – Vol.98. – P.35553574 DOI: 
109. 
Numerical simulation on solar collector and cascade heat pump combi water heating systems in Kazakhstan climates 
DOI: 10.1016/j.renene.2019.06.102 (SJR 1,83;Q1, Procentile 88 (Renewable Energy, Sustainability and the Environment) ). 
At low ambient temperatures, the heating capacity and coefficient of performance of a single stage vapour compression heat pump cycle is significantly getting reduced. The two stage cascade heat pump cycle operating with two different refrigerants provides a sustainable solution to lift the condenser temperature above 343 K. In this work, a numerical simulation model was developed for predicting the performance of a solar collector and two stage cascade heat pump combi water heating systems under Kazakhstan climatic conditions. The numerical simulation was performed for winter climatic conditions using nine refrigerant pairs such as, R32/R290, R32/R1234yf, R32/R134a, R410A/R290, R410A/R1234yf, R410A/R134a, R744/R290, R744/R1234yf and R744/R134a. The influences of solar irradiation, ambient temperature and condenser temperature are discussed. The solar collector and two stage heat pump combi heating system has about 30% energy savings when compared to the conventional two stage cascade heat pump without integration solar collectors. The R32/R290 refrigerant pair is having maximum coefficient of performance of 2.4 at 323 K condensing temperature and 263 K evaporating temperature. The refrigerant R744/R290 pair is identified as an environment friendly sustainable option in terms of its global warming impact for two stage cascade heat pump applications. © 2019 
Yerdesh Y.,Abdulina Z.,Aliuly A.,Belyayev Y.,Mohanraj M.,Kaltayev A Numerical simulation on solar collector and cascade heat pump combi water heating systems in Kazakhstan climates. //Renewable Energy. – 2019. – Vol.145. – P.12221234 DOI: 10.1016/j.renene.2019.06.102 
110. 
GPU Accelerated Modeling of InSitu Leaching Process and Streamline Based Reactive Transport Simulation 
DOI: 10.1016/j.procs.2020.11.016 (SJR 0,33; Procentile 68(Computer Science General Computer Science) ). 
In present paper GPU acceleration modeling and Streamline based reactive transport simulation for uranium InSitu leaching process is studied. InSitu leaching (ISL) is a method of selective dissolution of mineral inside of rock by injection of leaching solution to the layer through the wells and chemical interaction of leaching solution with ore. Changing of flow rates on wells over time leads to increasing computational time and simulation of InSitu leaching problem becomes resourceintensive problem. Acceleration of the hydrodynamic of the ISL process is carried out on GPU; parallelization of reactive mass transport due to reaction of leaching solute with mineral is performed by with multiThreading and streamlines simulation. © 2020 Elsevier B.V.. All rights reserved. 
Tungatarova, M.S., Kurmanseiit, M.B., Shayakhmetov, N.M., GPU Accelerated Modeling of InSitu Leaching Process and Streamline Based Reactive Transport Simulation. //Procedia Computer Science. – 2020. – Vol.178. – P.145152 DOI: 10.1016/j.procs.2020.11.016 
111. 
Synthesis of fourbar linkage with adjustable crank length for multipath generation 
DOI: 10.18178/ijmerr.9.4.489495 (SJR 0,19;Q4, Procentile 26(Engineering Mechanical Engineering) 22(Engineering Control and Systems Engineering) 14(Computer Science Artificial Intelligence) ). 
Synthesis of planar mechanism with adjustable crank length for generating multiple paths is presented. Leastsquare approximation problem is considered which allows carrying out approximate synthesis with unlimited number of desired coupler point positions and with unlimited number of prescribed trajectories. By reducing the task to synthesis of twoelement link with variable binary link length, which is called RPRmodule, the analytical solution is obtained to determine not only constant design parameters (mechanism link lengths) but the adjusting parameter values as well. Thus the number of design variables for nonlinear optimization (applied to find the remaining parameters) will be decreased significantly. The applied method is exemplified by synthesis of the mechanism for variable straight line generation, where the required height of the endeffector is adjusted by adjusting the crank length. Combined with random search technique the method allows to find all local minimums of the optimized goal function and thus allows to take full advantage from the considered mechanism structure during design. © 2020 by the authors. 
Ibrayev S.,Jomartov A.,Tuleshov A.,Jamalov N.,Ibrayev A.,Mukhambetkalieva G.,Aidasheva G.,Kamal A.,Synthesis of fourbar linkage with adjustable crank length for multipath generation. //International Journal of Mechanical Engineering and Robotics Research. – 2020. – Vol.9(4). – P.489495 DOI: 10.18178/ijmerr.9.4.489495 
112. 
Designing of the stephenson ii sixlink linkage actuator for servo mechanical press 
DOI: 10.24247/ijmperdapr202053 (SJR 0; Procentile 22(Engineering Mechanical Engineering) 22(Engineering Aerospace Engineering) 17(Chemical Engineering Fluid Flow and Transfer Processes) ). 
Currently, servo mechanical presses with an actuator based on a slidercrank mechanism are widely used for stamping details. One of the major problems of these servo mechanical presses is inclination of slide caused by eccentric application of deforming force and total linear elastic deformation of the press links and stamps under load. The inclination of the slide of the servo mechanical press reduces the accuracy of punching and causes its jamming and lead to its breakage. To eliminate these problems, it is proposed to use Stephenson II sixbar linkage as an actuator of servo mechanical press. The scheme structure of the actuator of servo mechanical press with two connecting rods and one crank based on the Stephenson II sixbar linkage is obtained. Synthesis of the Stephenson II sixbar linkage was carried out. As a result of synthesis, a new actuator of servo mechanical press on the base of the Stephenson II sixbar linkage was obtained and its prototype was made. The prototype of the new actuator of servo mechanical press based on the Stephenson II sixbar linkage was tested. Testing of the prototype of the new actuator of servo mechanical press showed a good distribution of the efforts applied during the presswork and a better tolerance of the eccentric load. © TJPRC Pvt. Ltd. 
Jomartov, A., Tuleshov, A., Jamalov, N., Kuatova, K., Kaimov, A., Designing of the stephenson ii sixlink linkage actuator for servo mechanical press. //International Journal of Mechanical and Production Engineering Research and Development. – 2020. – Vol.10(2). – P.501512 DOI: 10.24247/ijmperdapr202053 
113. 
Dynamic Model of a Crank Press in the Process of Braking 
DOI: 10.1007/9783030300364_12 (SJR 0,16;Q4, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
The paper studies the dynamic of a brake of crank press. At present, the dynamic research of brake of the crank presses, with account of interaction with other links, is a priority. The crank press contains movable parts and links, the mass of which is from one hundred kilograms to several tons. These parts and links are cyclically stopped when braking with a crank press almost instantaneously, and they are subject to high dynamic loads. To simulate and analyze the movement of crank press with brake, a software package: SimulationX is used. SimulationX is a software package for modeling and analyzing the dynamics and kinematics of cars, industrial equipment, electric, pneumatic and hydraulic drives, hybrid engines, etc. As a result of dynamic calculation, important dynamic parameters of the crank press brake and working ram are determined. It is shown that dynamic loads sharply increase almost in all links of the crank press when the brake is switched on. © 2020, Springer Nature Switzerland AG. 
Jomartov, A., Tuleshov, A., Kuatova, M., Dynamic Model of a Crank Press in the Process of Braking. //Mechanisms and Machine Science. – 2020. – Vol.78. – P.141150 DOI: 10.1007/9783030300364_12 
114. 
Structurally parametric synthesis and position analysis of a robomech class parallel manipulator with two endeffectors 
DOI: (SJR 0,26;Q3, Procentile 49(Mathematics Applied Mathematics) ). 
In this paper, the methods of structurally parametric synthesis and position analysis of a RoboMech class parallel manipulator with two endeffectors are presented. This parallel manipulator is formed by connecting the two moving output objects with the fixed base by two passive, one active and two negative closing kinematic chains. Geometrical parameters of the active and negative closing kinematic chains are determined by the Chebyshev and leastsquare approximations. Position analysis is made on base of the conditional generalized coordinates method. © 2020 International Association of Engineers. 
Baigunchekov Z.бNaurushev B.,Zhumasheva Z.,Mustafa A.,Kairov R., Amanov B, Structurally parametric synthesis and position analysis of a robomech class parallel manipulator with two endeffectors. //IAENG International Journal of Applied Mathematics. – 2020. – Vol.5(177). – P.111 DOI: 
115. 
The First Type of Singularity of a 3PRRS Parallel Manipulator 
DOI: 10.1007/9783030489892_38 (SJR 0,16;Q4, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
In this paper, the first type of singularity of a 3PRRS parallel manipulator is considered. This type of singularity appears when the determinant of the Jacobian matrix of the generalized velocities goes to zero. In this case, the parallel manipulator loses one or more degrees of freedom (DOF). From the analysis of the degeneracy of the generalized velocities Jacobian matrix, the conditions of the first type singular configurations are determined. 
Baigunchekov, Z., Laribi, M.A., Izmambetov, M., Zhumasheva, Z., Kaiyrov, R..,The First Type of Singularity of a 3PRRS Parallel Manipulator. //Mechanisms and Machine Science. – 2020. – Vol.84. – P.356363 
116. 
Performance analysis of crushed gravel sand heat storage and biomass evaporatorassisted single slope solar still 
DOI:10.1007/s1135602115487w (SJR 0,85;Q2, Procentile 79 (Environmental SciencePollution) 76(Environmental ScienceHealth, Toxicology and Mutagenesis) 67(Environmental ScienceEnvironmental Chemistry) ). 
In this research work, the productivity, energy, exergy, and economic and enviroeconomic performance in crushed gravel sand heat storage and biomass evaporatorassisted solar still (CGSBSS) have been investigated and compared the results with conventional solar still (CSS) under the similar climatic conditions of Coimbatore City during the year 2019. The heat accumulated in crushed gravel sand and biomass evaporator have been used to preheat the inlet saline water and air vapor before entering into the solar still. This results in enhanced air vapor mixture temperature and evaporative heat transfer rate of CGSBSS significantly. The productivity, energy, and exergy efficiencies in CGSBSS were improved by 34.6%, 34.4%, and 35%, respectively when compared to CSS. In economic analysis, the payback period (PBP) in both CGSBSS and CSS was estimated to be about 4.7 months and 3.9 months, respectively. Furthermore, in enviroeconomic analysis, the CO2 emission estimated in CGSBSS and CSS was about 16.63 tons and 8.18 tons, respectively during its lifetime of 10 years. 
Ramasamy Dhivagar, Murugesan Mohanraj,Yerzhan Belyayev ,Performance analysis of crushed gravel sand heat storage and biomass evaporatorassisted single slope solar still. //Environmental Science and Pollution Research. – 2021. DOI: 10.1007/s1135602115487w

117. 
Walking Robot Leg Design Based on Translatory StraightLine Generator 
DOI:10.1007/9783030583804_32 (SJR 0,21;Q3, Procentile 37(Engineering Mechanical Engineering) 32(Engineering Mechanics of Materials) 30(Computer Science Computer Science Applications) 27(Mathematics Modeling and Simulation) ). 
The planar straightline generating linkages are of great interest for walking robot propel, especially to reduce power consumption and simplify control. Optimal design of sixlink leg linkage is proposed with unlimited foot adaptation on terrain irregularities due to rectilinear and translatory motion of the output link, referred to as shinlink. The analytical method of synthesis is proposed based on leastsquare approximation. Due to minimizing directly the deviation from the desired output motion, the method does not suffer from so called branching defect and allows synthesizing the mechanism with desired transmission angle. Multicriteria optimization attaining both the best accuracy and transmission angle is carried out and leglinkage with maximal value of duration ratio of support and transference phases of the leg step cycle is presented. © 2021, CISM International Centre for Mechanical Sciences. 
Ibrayev S.,Jamalov N.,Tuleshov A.,Jomartov A.,Ibrayev A.,Kamal A.,Ibrayeva A.,Bissembayev K., Walking Robot Leg Design Based on Translatory StraightLine Generator. //CISM International Centre for Mechanical Sciences, Courses and Lectures. – 2021. – Vol.601. – P.264271 DOI: 10.1007/9783030583804_32 
118. 
Dynamic Model of Servo Mechanical Press 
DOI:10.1007/9783030583804_21 (SJR 0,21;Q3, Procentile 37(Engineering Mechanical Engineering) 32(Engineering Mechanics of Materials) 30(Computer Science Computer Science Applications) 27(Mathematics Modeling and Simulation)). 
Servo mechanical press is a mechatronic system, consisting of a reducer, actuator and a servomotor with a controller. Control of the servo mechanical press is carried out with a help of the controller depending of technological process of pressing. Reverse kinematic analysis of the actuator of the servo mechanical press is presented in the work to determine the law of motion of a crank, reproducing the specified law of motion of the servo mechanical press slide. Obtained law of motion of the actuator’s crank of the servo mechanical press is necessary to choose the servomotor. Considering that servo mechanical press works by complex programmable motion large dynamic loads occur in its nodal points and links. The work proposes to simulate the dynamics of the servo mechanical press on the software complex SimulationX. Dynamic model of the servo mechanical press is composed with account for elasticdissipative characteristics of the links and the parameters of the servomotor, for different laws of motion of the slide. A dynamic model of a 50ton servo mechanical press is developed using the SimulationX software package. Calculation was carried out using real data of the existing 50 ton servo mechanical press. As a result of the calculation of the dynamic model of the servo mechanical press, the following data were obtained: motion, velocity, acceleration of the slide, torque effect on the crank, angular rate, torque effect, power of the servomotor. 
Jomartov, A., Tuleshov, A., Jamalov, N.,Temirbekov, Y., Bostanov, B.,Dynamic Model of Servo Mechanical Press. //CISM International Centre for Mechanical Sciences, Courses and Lectures. – 2021. – Vol.601. – P.170178 
119. 
Unsteady Resonant Oscillations of a Gyroscopic Rigid Rotor with Nonlinear Damping and Nonlinear Rigidity of the Elastic Support 
DOI:10.1007/9783030835941_9 (SJR 0,16;Q4, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
The article is concerned with the effect of linear and cubic nonlinear damping of an elastic bearing on forced resonant vibrations of a gyroscopic vertical rigid rotor taking into account nonlinear stiffness of the cubic nature of the bearing material. It is confirmed that nonlinear cubic damping of the support can suppress not only the maximum amplitude, but also the amplitudes of forced unsteady oscillations behind the rotation speed corresponding to the maximum amplitude and the variation of its values in time along the main curve, around its mean values. It shifts the speed of rotation of the amplitude maximum, with rigid and soft nonlinear elastic characteristics of the support material downwards and upwards, respectively. It is shown that with a “slow” increase in the shaft rotation speed, an increase in the absolute value of the angular acceleration is accompanied by a shift of the amplitude peak towards high speeds, with a “slow” decrease in the shaft rotation speed – towards low speeds with a decrease in the amplitude of oscillations. It is shown that during the rotor takeoff run, the maximum amplitude for the case with a rigid nonlinear elasticity characteristic of the support material is greater than the same value for the case with a soft nonlinear elasticity characteristic of the support material, and conversely, during the rotor rundown for similar cases. 
Iskakov, Z., Jamalov, N., Bissembayev, K., Unsteady Resonant Oscillations of a Gyroscopic Rigid Rotor with Nonlinear Damping and Nonlinear Rigidity of the Elastic Support. //Mechanisms and Machine Science. – 2021. – Vol.85. – P.8393 
120. 
Nonisothermal pore change model predicting CO2 adsorption onto consolidated activated carbon 
DOI:https://doi.org/10.1016/j.ijheatmasstransfer.2021.121480 (SJR 1,71;Q1, Procentile 98(Chemical Engineering Fluid Flow and Transfer Processes) 95(Engineering Mechanical Engineering) 94(Physics and Astronomy Condensed Matter Physics) ). 
Accurate simulation and detailed description of the dynamics of the adsorption process play a significant role in forecasting the performance of new materials when used in various adsorption systems, like cooling/heating. The activated carbon (AC) consolidation allows improving the heat transfer rate inside the adsorption/desorption bed and compacting the systems. There are numerous mathematical models in literature for gas adsorption onto granular AC. But for consolidated AC, because of the absence of macropores, most assumptions that work well for granular AC may lead to significant discrepancies. Therefore, the present research proposes a new mathematical model for gas (CO2) adsorption onto consolidated adsorbent, a nonisothermal pore change model. The model takes into account the porosity and permeability changes due to the adsorption. The validation of the developed model is performed via comparison with the results obtained experimentally and numerically using an isothermal model. The effective Knudsen diffusion coefficient for the working pair is evaluated from porosity data. The rate of adsorption or mass transfer coefficient is estimated using the van't Hoff plot. The study results could be applied in the development of waste heatdriven cooling systems employing consolidated composite material as the adsorbent. The proposed mathematical model is also applicable for many other working pairs. 
Berdenova, B., Pal, A., Saha, B.B., Kaltayev, A.,Nonisothermal pore change model predicting CO2 adsorption onto consolidated activated carbon. //International Journal of Heat and Mass Transfe. – 2021. – Vol.177. – 121480 DOI: https://doi.org/10.1016/j.ijheatmasstransfer.2021.121480 
121. 
The use of the linear form of dynamical equations of the satellite attitude control system for its analysis and synthesis 
DOI:https://doi.org/10.15632/jtampl/129071 (SJR 0,29;Q3, Procentile 46(Engineering Mechanical Engineering) ). 
At present, the methods based on using linearized dynamical equations are applied for syn thesis of an attitude control system of a satellite with nonlinear dynamics. Linearized equa tions describe the satellite dynamics approximately, which is the main their disadvantage. This article shows that basing on the angular momentum theorem, the nonlinear dynamical equations of the satellite attitude control system can be represented in the form of linear differential equations with variable coefficients, which makes it possible to use engineering methods of stability analysis and analysis of transient quality in the process of synthesis of the satellite attitude control system. 
Moldabekov, M., Sukhenko, A., Shapovalova, D., Yelubayev, S.,The use of the linear form of dynamical equations of the satellite attitude control system for its analysis and synthesis. //Journal of Theoretical and Applied Mechanics (Poland). – 2021. – Vol.59(1). – P.109120 
122. 
Multiparametric dynamic analysis of a rolling bearings system 
DOI:10.5545/svjme.2021.7178 (SJR 0,27;Q3, Procentile 47(Engineering Mechanical Engineering) 42(Engineering Mechanics of Materials) ). 
A method for calculating amplitudes and constructing frequency characteristics of forced and selfexcited vibrations of a rotorfluidfoundation system on rolling bearings with a nonlinear characteristic based on the method of complex amplitudes and harmonic balance has been developed. Nonlinear equations of motion of the rotorfluidfoundation system are derived, and analytical methods of their solution are presented. Frequencies of fundamental and ultraharmonic resonances are determined. The intervals between selfoscillation frequencies are estimated. The dependence of amplitudes on the amount of fluid in the rotor cavity, the mass of the foundation, linear imbalance, the value of the stiffness coefficient, and the damping coefficient is shown. © 2021 Journal of Mechanical Engineering. 
Kydyrbekuly, A., Ibrayev, G.G.A., Ospan, T., Nikonov, A.,Multiparametric dynamic analysis of a rolling bearings system. //Strojniski Vestnik/Journal of Mechanical Engineeringe. – 2021. – Vol.67(9). – P.421432 DOI: 10.5545/svjme.2021.7178 
123. 
Modeling the Separation Process in Vertical Rotor Systems 
DOI:10.1007/9783030835941_11 (SJR 0,16;Q4, Procentile 23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
In this paper, we study and analyze the features of the separation process in a centrifugal force field, i.e. centrifugation process in vertical rotor systems. The main parameters that determine the time of separation of particles are revealed, and the optimal modes are indicated both from a constructive and from an economic point of view. Special cases of a fixed rotor are considered. Nonlinear differential equations of motion of a suspension particle are obtained, which do not have an exact solution. The study is carried out by analytical and numerical methods. The dependences of the slope angles of the tubes on the angular velocity of rotation of the rotor, sedimentation curves that allow one to estimate the time of deposition of particles, as well as the effect of the dispersed composition on the separation process as a whole, are obtained. The results of the study of this work allow us to determine with sufficient accuracy all the necessary characteristics working process of separation and sedimentation, and also allow in certain cases to exclude experimental work. © 2022, The Author(s), under exclusive license to Springer Nature Switzerland AG. 
Kydyrbekuly, A.B., Ibrayev, G.E, Modeling the Separation Process in Vertical Rotor Systems. //Mechanisms and Machine Science. – 2021. – Vol.85. – P.104113. DOI: 10.1007/9783030835941_11 
124. 
A Comparison of Machine Learning Algorithms in Predicting Lithofacies: Case Studies from Norway and Kazakhstan 
DOI:10.3390/en14071896 (SJR 0,6;Q2, Procentile 49(Mathematics Applied Mathematics)23(Engineering Mechanical Engineering) 18(Engineering Mechanics of Materials) ). 
Defining distinctive areas of the physical properties of rocks plays an important role in reservoir evaluation and hydrocarbon production as core data are challenging to obtain from all wells. In this work, we study the evaluation of lithofacies values using the machine learning algorithms in the determination of classification from various well log data of Kazakhstan and Norway. We also use the wavelettransformed data in machine learning algorithms to identify geological properties from the well log data. Numerical results are presented for the multiple oil and gas reservoir data which contain more than 90 released wells from Norway and 10 wells from the Kazakhstan field. We have compared the the machine learning algorithms including KNN, Decision Tree, Random Forest, XGBoost, and LightGBM. The evaluation of the model score is conducted by using metrics such as accuracy, Hamming loss, and penalty matrix. In addition, the influence of the dataset features on the prediction is investigated using the machine learning algorithms. The result of research shows that the Random Forest model has the best score among considered algorithms. In addition, the results are consistent with outcome of the SHapley Additive exPlanations (SHAP) framework. Copyright: © 2021 by the authors. 
Merembayev, T., Kurmangaliyev, D., Bekbauov, B., Amanbek, Y., A Comparison of Machine Learning Algorithms in Predicting Lithofacies: Case Studies from Norway and Kazakhstan. //Energies. – 2021. – Vol.14(7). – P.1896 DOI: 10.3390/en14071896 
125. 
Kinematic synthesis method and eccentricity effects of a Stephenson mechanism 
DOI:10.5194/ms1212021 (SJR 0,31;Q2, Procentile 52(Engineering Industrial and Manufacturing Engineering) 50(Engineering Mechanical Engineering) 48(Engineering Civil and Structural Engineering) 48(Chemical Engineering Fluid Flow and Transfer Processes) 46(Engineering Mechanics of Materials) 45(Engineering Control and Systems Engineering) ). 
When implementing the technological process on crank presses, it is necessary to provide a predetermined working cycle of the slider motion: fast lifting, dwell, and slow lowering. The cycle cannot be realized without controlling the motor. In addition, using controllable motors increases the manufacturing cost. Due to the geometric and kinematic capabilities of the mechanism, changing the kinematics of the working link is the best choice. Thanks to the use of the Stephenson II mechanism, the slider skew is eliminated due to the parallel connecting rods and the increased area of slider contact. This study presents a numerical method for kinematic synthesis of the Stephenson mechanism that has kinematic advantages. The method is based on mean square deviation which is the minimizing of an objective function. Thanks to the proposed synthesis method, approximate dwell movement can be performed when the slider is on the bottom dead center. In this study, values of the crank length and parallel connecting rods' lengths, angular coordinates of the crank and connecting rods, and the eccentricity of the guide slider relative to the crank rotation axis were obtained. It is observed that eccentricity affects the lower forward and higher backward speed of the slider. The kinematic results of the slider movement are comparatively presented in this article. © 2021 Copernicus GmbH. All rights reserved. 
Tuleshov, A., Halicioglu, R., Shadymanova, A., Kuatova, M., Kinematic synthesis method and eccentricity effects of a Stephenson mechanism. //Mechanical Sciences. – 2021. – Vol.12(1). – P.18 DOI: 10.5194/ms1212021 
126. 
Analytical method for determination of internal forces of mechanisms and manipulators 
DOI:10.3390/robotics7030053 (SJR 0,39; Q3, Procentile 70(Control and Optimization) 67(Mechanical Engineering) 53(Artificial Intelligence) ). 
This paper presents a theory for the analytical determination of internal forces in the links of planar linkage mechanisms and manipulators with statically determinate structures, considering the distributed dynamic loads. Linkage mechanisms and manipulators were divided into elements and joints. Discrete models were created for both the elements and the entire mechanism. The dynamic equations of equilibrium for the discrete model of the elements and the hinged and rigid joints, under the action of longitudinal and transverse distributed dynamic trapezoidal loads, were derived. In the dynamic equations of the equilibrium of the discrete model of the elements and joints, the connections between the components of the force vector in the calculated crosssections and the geometric, physical, and kinematic characteristics of the element were established for its planeparallel motion. According to the developed technique, programs were created in the Maple system, and animations of the motion of the mechanisms were produced. The links were constructed with the intensity of transverse and longitudinaldistributed dynamic loads, bending moments, and shearing and normal forces, depending on the kinematic characteristics of the links. © 2018 by the authors. 
Utenov, M., Sobh, T., Baigunchekov, Z., Zhilkibayeva, S., Patel, S., Analytical method for determination of internal forces of mechanisms and manipulators. //Robotics,. – 2018. – Vol.7(3). – P.53DOI: 10.3390/robotics7030053 
127. 
Synthesis of Reconfigurable Positioning Parallel Manipulator of a Class RoboMech

DOI:10.1109/REMAR.2018.8449852 (SJR 0,18;). 
In this paper the methods of structural and dimensional synthesis of reconfigurable positioning parallel manipulator of a class RoboMech are presented. These parallel manipulators are formed by connecting the output point to the base by three closed kinematic chains: passive, active and negative. Parameters of synthesis of the active and negative closed kinematic chains are determined on the basis of the formulated Chebyshev and leastsquare approximations. © 2018 IEEE. 
Baigunchekov Z.,Ibrayev S.,Izmambetov M.,Zhumasheva Z.,Baigunchekov T.,Naurushev B.,Mustafa A.,Sagyntay M., Synthesis of Reconfigurable Positioning Parallel Manipulator of a Class RoboMech. //2018 International Conference on Reconfigurable Mechanisms and Robots, ReMAR 2018  Proceedings. – 2018. – 8449852DOI: 10.1109/REMAR.2018.8449852 
128. 
Synthesis of cartesian manipulator of a class RoboMech 
DOI:10.1007/9783030003654_9 (SJR 0,172; Q4, Procentile 25(Mechanical Engineering) 19(Mechanics of Materials) ). 
In this paper the methods of structural and dimensional synthesis of Cartesian manipulators of a class RoboMech reproducing the series of orthogonal trajectories are developed. The Cartesian manipulator of a class RoboMech is formed by connecting the output point to the base by three closing kinematic chains: one passive and two active ones. Parameters of synthesis of active closing kinematic chains are determined on the basis of the formulated problem of leastsquare approximation. © 2019, Springer Nature Switzerland AG. 
Baigunchekov Z.,Ibrayev S.,Izmambetov M.,Baigunchekov T.,Naurushev B.,Mustafa A.,Synthesis of cartesian manipulator of a class RoboMech. //RoboticsMechanisms and Machine Science,. – 2019. – Vol.66. – P.6976 DOI: 10.1007/9783030003654_9 
129. 
Parallel manipulator of a class RoboMech for generation of horizontal trajectories family

DOI:10.1007/9783030201319_137 (SJR 0,172; Q4, Procentile 23(Mechanical Engineering) 18(Mechanics of Materials) ). 
In this paper the methods of structural and dimensional synthesis of parallel manipulator of a class RoboMech for reproducing the series of horizontal trajectories are developed. This parallel manipulator is formed by connecting the output point to the base using two closing kinematic chains: one passive and one active ones. The synthesis parameters of the active closing kinematic chain of RPR type are determined on the basis of the leastsquare approximation. © 2019, Springer Nature Switzerland AG. 
Baigunchekov, Z., Izmambetov, M., Zhumasheva, Z., Baigunchekov, T., Mustafa, A.,Parallel manipulator of a class RoboMech for generation of horizontal trajectories family. //Mechanisms and Machine Science,. – 2019. – Vol.73. – P.13951402 DOI: 10.1007/9783030201319_137 
130. 
Geometry and Inverse Kinematics of 3PRRS Type Parallel Manipulator 
DOI:10.1007/9783030196486_2 (SJR 0; Q4, Procentile 27(Control and OptimizationGeneral Computer Science) 22(Control and Systems Engineering)). 
In this paper the methods of structural synthesis and inverse kinematics of 3PRRS type parallel manipulator (PM) are developed. This PM is formed by connecting of a moving platform with a fixed base by three passive closing kinematic chains of PRRS type. Constant and variable parameters characterizing the geometry of links and relative motions of elements of joints are defined. Inverse kinematics of the PM is solved of the basis of solution of the loopclosure matrix equations of the legs. © Springer Nature Switzerland AG 2020. 
Baigunchekov Z.,Laribi M.A.,Mustafa A.,Kaiyrov R.,Amanov B.,Kassinov A., Geometry and Inverse Kinematics of 3PRRS Type Parallel Manipulator. //Advances in Intelligent Systems and Computing,. – 2019. – Vol.980. – P.1218 DOI: 10.1007/9783030196486_2 
131. 
Structurally Parametric Synthesis of a RoboMech Class Parallel Manipulator with Three DOF

DOI:10.1007/9783030489892_40 (SJR 0,172; Q4, Procentile 23(Mechanical Engineering) 18(Mechanics of Materials)). 
This paper presents the methods of structuralparametric synthesis of a parallel manipulator with three DOF (Degree of Freedom) working in a cylindrical coordinate system. This parallel manipulator belongs to the RoboMech class because it works under the setting laws of motions of the endeffector and actuators, which simplifies the control system and improves its dynamics. Parallel manipulators of a RoboMech class work with certain structural schemes and geometrical parameters of their links. The considered parallel manipulator is formed by connecting the output object to a base using one passive and two active closing kinematic chains. The geometrical parameters of their links are determined on the base of the leastsquare approximations. © The Editor(s) (if applicable) and The Author(s), under exclusive license to Springer Nature Switzerland AG 2020. 
Baigunchekov, Z., Tarek, S., Patel, S., Mustafa, A.,Structurally Parametric Synthesis of a RoboMech Class Parallel Manipulator with Three DOF. //Mechanisms and Machine Science. – 2020. – Vol.84. – P.371379 DOI: 10.1007/9783030489892_40 
132. 
The First Type of Singularity of a 3PRRS Parallel Manipulator, 2020, 84, стр. 356–363

DOI:10.1007/9783030489892_38 (SJR 0,172; Q4, Procentile 23(Mechanical Engineering) 18(Mechanics of Materials)) 
In this paper, the first type of singularity of a 3PRRS parallel manipulator is considered. This type of singularity appears when the determinant of the Jacobian matrix of the generalized velocities goes to zero. In this case, the parallel manipulator loses one or more degrees of freedom (DOF). From the analysis of the degeneracy of the generalized velocities Jacobian matrix, the conditions of the first type singular configurations are determined. © The Editor(s) (if applicable) and The Author(s), under exclusive license to Springer Nature Switzerland AG 2020. 
Baigunchekov, Z., Laribi, M.A., Izmambetov, M., Zhumasheva, Z., Kaiyrov, R.,The First Type of Singularity of a 3PRRS Parallel Manipulator. //Mechanisms and Machine Science. – 2020. – Vol.84. – P.356363 DOI: 10.1007/9783030489892_38 
133. 
Inverse Kinematics of a 3PRPS Type Parallel Manipulator

DOI:10.1007/9783030489892_39 (SJR 0,172; Q4, Procentile 23(Mechanical Engineering) 18(Mechanics of Materials)) 
In this paper, geometry and inverse kinematics of a 3PRPS type parallel manipulator are studied. This parallel manipulator is formed by connecting a moving platform with a base by three passive closing kinematic chains of a PRPS type. Constant and variable parameters characterizing the geometry of links and relative movements of elements of kinematic pairs, respectively, are defined, and the matrices of binary links and kinematic pairs are derived. On the base of these matrices of binary links and kinematic pairs, the inverse kinematics problem is solved and numerical results are presented. © The Editor(s) (if applicable) and The Author(s), under exclusive license to Springer Nature Switzerland AG 2020. 
Baigunchekov, Z., Zeghloul, S., Kassinov, A.,Inverse Kinematics of a 3PRPS Type Parallel Manipulator. //Mechanisms and Machine Science. – 2020. – Vol.84. – P.364370 DOI: 10.1007/9783030489892_39 
134. 
Structurally parametric synthesis and position analysis of a robomech class parallel manipulator with two endeffectors 
DOI: (SJR 0,259; Q3, Procentile 49(Applied Mathematics)) 
In this paper, the methods of structurally parametric synthesis and position analysis of a RoboMech class parallel manipulator with two endeffectors are presented. This parallel manipulator is formed by connecting the two moving output objects with the fixed base by two passive, one active and two negative closing kinematic chains. Geometrical parameters of the active and negative closing kinematic chains are determined by the Chebyshev and leastsquare approximations. Position analysis is made on base of the conditional generalized coordinates method. © 2020 International Association of Engineers. 
Baigunchekov Z.,Naurushev B.,Zhumasheva Z.,Mustafa A.,Kairov R., Amanov B,Structurally parametric synthesis and position analysis of a robomech class parallel manipulator with two endeffectors. //IAENG International Journal of Applied Mathematics. – 2020. – Vol.5. – P.111 DOI: 
135. 
A robomech class parallel manipulator with three degrees of freedom2020, 3(7105), стр. 44–56

DOI:10.15587/17294061.2020.203131 (SJR 0,268; Q3, Procentile 56(Applied Mathematics) 54(Management of Technology and Innovation) 52(Industrial and Manufacturing Engineering) 49(Industrial and Manufacturing Engineering) 46(Industrial and Manufacturing Engineering) 46(Energy Engineering and Power Technology) 45(Computer Science Applications) 44(Control and Systems Engineering)) 
This paper presents the methods of structuralparametric synthesis and kinematic analysis of a parallel manipulator with three degrees of freedom working in a cylindrical coordinate system. This parallel manipulator belongs to a RoboMech class because it works under the set laws of motions of the endeffector and actuators, which simplifies the control system and improves its dynamics. Parallel manipulators of a RoboMech class work with certain structural schemes and geometrical parameters of their links. The considered parallel manipulator is formed by connecting the output point to a base using one passive and two active closing kinematic chains (CKC). Passive CKC have zero degree of freedom and it does not impose a geometrical constraint on the movement of the output point, so the geometrical parameters of the links of the passive CKC are freely varied. Active CKCs have active kinematic pairs and they impose geometrical constraints on the movement of the output point. The geometrical parameters of the links of the active CKCs are determined on the basis of the approximation problems of the Chebyshev and leastsquare approximations. For this, the equations of geometrical constraints are derived in the forms of functions of weighted differences, which are presented in the forms of generalized (Chebyshev) polynomials. This leads to linear iterative problems. The direct and inverse problems of the kinematics of the investigated parallel manipulator are solved. In the direct kinematics problem, the coordinates of the output point are determined by the given position of the input links. In the inverse kinematics problem, the positions of the input links are determined by the coordinates of the output point. The direct and inverse problems of the kinematics of the investigated parallel manipulator are reduced to solving problems on the positions of Sylvester dyads. Numerical results of structuralparametric synthesis and kinematic analysis of the considered parallel manipulator are presented. The numerical results of the kinematic analysis show that the maximum deviation of the movement of the output point from the orthogonal trajectories is 1.65 % 
Baigunchekov, Z., Mustafa, A., Sobh, T., Patel, S., Utenov, M.,A robomech class parallel manipulator with three degrees of freedom. //A robomech class parallel manipulator with three degrees of freedom. – 2020. – Vol.3(7105). – P.4456 DOI: 10.15587/17294061.2020.203131 
136. 
Direct kinematics of a 3PRRS type parallel manipulator

DOI:10.18178/ijmerr.9.7.967972 (SJR 0,187; Q4, Procentile 23(Mechanical Engineering) 18(Mechanics of Materials)) 
This paper is devoted to solving the direct kinematics of a novel 3PRRS type parallel manipulator with sixdegreesoffreedom, where P, R, and S are prismatic, revolute and spherical kinematic pairs respectively. This parallel manipulator is formed by connecting a moving platform with a fixed platform (base) through three closing kinematic chains of a PRRS type in which the prismatic kinematic pairs and their adjacent revolute kinematic pairs are active and they are located on a fixed platform. The constant and variable parameters of the considered parallel manipulator characterizing its geometry and kinematics respectively are determined. In the direct kinematics, the positions of the moving platform are determined by the known constant parameters of the links and the given variable parameters of the active kinematic pairs. An analysis of the obtained equations of the direct kinematics showed that the variable parameters of the active prismatic kinematic pairs are set free, and these equations are reduced to a 16 –order polynomial equation with the passive kinematic pairs variables. Numerical examples of the considered parallel manipulator’s direct kinematics are presented, and the results showed that the direct kinematics equations have four solutions corresponding to the four assemblies of the parallel manipulator. 
Baigunchekov, Z.Z., Kaiyrov, R.A.,Direct kinematics of a 3PRRS type parallel manipulator. //International Journal of Mechanical Engineering and Robotics Research. – 2020. – Vol.9(7). – P.967972 
137. 
Inverse Kinematics and Workspace of a 3PRRS Type Parallel Manipulator

DOI:10.1007/9783030752712_8 (SJR 0,159; Q4, Procentile 23(Mechanical Engineering) 18(Mechanics of Materials)) 
In this paper, methods of workspace analysis of a 3PRRS type parallel manipulator are described. The equations of spheres and circles on these spheres, along which the center of the moving platform can move, are derived, and it is shown that the total reachable area of these spheres is the workspace of the considered parallel manipulator. Numerical examples of defining the workspace of the 3PRRS type parallel manipulator are presented. © 2021, The Author(s), under exclusive license to Springer Nature Switzerland AG. 
Baigunchekov, Z., Laribi, M.A., Kaiyrov, R., Zholdassov, E.,Inverse Kinematics and Workspace of a 3PRRS Type Parallel Manipulator. //Mechanisms and Machine Science. – 2021. – Vol.103. – P.7178 DOI: 10.1007/9783030752712_8 
138. 
Kinematic synthesis and analysis of the robomech class parallel manipulator with two grippers, 2021, 10(3), 99

DOI:10.3390/robotics10030099 (SJR 0,390; Q2, Procentile 67(Mechanical Engineering) 70(Control and Optimization) 53(Artificial Intelligence)) 
In this paper, methods of kinematic synthesis and analysis of the RoboMech class parallel manipulator (PM) with two grippers (end effectors) are presented. This PM is formed by connecting two output objects (grippers) with a base using two passive and one negative closing kinematic chains (CKCs). A PM with two end effectors can be used for reloading operations of stamped products between two adjacent main technologies in a cold stamping line. Passive CKCs represent two serial manipulators with two degrees of freedom, and negative CKC is a three‐joined link with three negative degrees of freedom. A negative CKC imposes three geometric constraints on the movements of the two output objects. Geometric parameters of the negative CKC are determined on the basis of the problems of the Chebyshev and least‐square approximations. Problems of positions and analogues of velocities and accelerations of the PM with two end effectors have been solved. © 2021 by the authors. Licensee MDPI, Basel, Switzerland. 
Baigunchekov, Z., Laribi, M.A., Mustafa, A., Kassinov, A.,Kinematic synthesis and analysis of the robomech class parallel manipulator with two grippers . //Robotics. – 2021. – Vol.10(3). – P.99 DOI: 10.3390/robotics10030099 
139. 
Structuralparametric synthesis of the robomech class parallel mechanism with two sliders

DOI:10.3390/app11219831 (SJR 0,435; Q2, Procentile 71(General Engineering) 62(Instrumentation) 59(Fluid Flow and Transfer Processes) 56(Computer Science Applications) 51(General Materials Science) 50(Process Chemistry and Technology)) 
This paper addresses the structuralparametric synthesis and kinematic analysis of the RoboMech class of parallel mechanisms (PM) having two sliders. The proposed methods allow the synthesis of a PM with its structure and geometric parameters of the links to obtain the given laws of motions of the input and output links (sliders). The paper outlines a possible application of the proposed approach to design a PM for a cold stamping technological line. The proposed PM is formed by connecting two sliders (input and output objects) using one passive and one negative closing kinematic chain (CKC). The passive CKC does not impose a geometric constraint on the movements of the sliders and the geometric parameters of its links are varied to satisfy the geometric constraint of the negative CKC. The negative CKC imposes one geometric constraint on the movements of the sliders and its geometric parameters are determined on the basis of the Chebyshev and leastsquare approximations. Problems of positions and analogues of velocities and accelerations of the considered PM are solved to demonstrate the feasibility and effectiveness of the proposed formulations and case of study. © 2021 by the authors. Licensee MDPI, Basel, Switzerland. 
Baigunchekov, Z., Laribi, M.A., Carbone, G.,Mustafa A.,Amanov, B., Zholdassov, Y..,Structuralparametric synthesis of the robomech class parallel mechanism with two sliders. //Applied Sciences (Switzerland). – 2021. – Vol.11(21). – P.9831 DOI: 10.3390/app11219831 
140. 
Nonisothermal pore change model predicting CO2 adsorption onto consolidated activated carbon

DOI:10.1016/j.ijheatmasstransfer.2021.121480 (SJR 1,713; Q1, Procentile 98(Fluid Flow and Transfer Processes) 95(Mechanical Engineering) 94(Condensed Matter Physics)) 
Accurate simulation and detailed description of the dynamics of the adsorption process play a significant role in forecasting the performance of new materials when used in various adsorption systems, like cooling/heating. The activated carbon (AC) consolidation allows improving the heat transfer rate inside the adsorption/desorption bed and compacting the systems. There are numerous mathematical models in literature for gas adsorption onto granular AC. But for consolidated AC, because of the absence of macropores, most assumptions that work well for granular AC may lead to significant discrepancies. Therefore, the present research proposes a new mathematical model for gas (CO2) adsorption onto consolidated adsorbent, a nonisothermal pore change model. The model takes into account the porosity and permeability changes due to the adsorption. The validation of the developed model is performed via comparison with the results obtained experimentally and numerically using an isothermal model. The effective Knudsen diffusion coefficient for the working pair is evaluated from porosity data. The rate of adsorption or mass transfer coefficient is estimated using the van't Hoff plot. The study results could be applied in the development of waste heatdriven cooling systems employing consolidated composite material as the adsorbent. The proposed mathematical model is also applicable for many other working pairs. © 2021 Elsevier Ltd 
Berdenova, B., Pal, A., Saha, B.B., Kaltayev, A.,Nonisothermal pore change model predicting CO2 adsorption onto consolidated activated carbon . //International Journal of Heat and Mass Transfer. – 2021. – Vol.177. – P.121480 DOI: 10.1016/j.ijheatmasstransfer.2021.121480 
141. 
Multiparametric dynamic analysis of a rolling bearings system

DOI:10.5545/svjme.2021.7178 (SJR 0,271; Q3, Procentile 47(Mechanical Engineering) 42(Mechanics of Materials)) 
A method for calculating amplitudes and constructing frequency characteristics of forced and selfexcited vibrations of a rotorfluidfoundation system on rolling bearings with a nonlinear characteristic based on the method of complex amplitudes and harmonic balance has been developed. Nonlinear equations of motion of the rotorfluidfoundation system are derived, and analytical methods of their solution are presented. Frequencies of fundamental and ultraharmonic resonances are determined. The intervals between selfoscillation frequencies are estimated. The dependence of amplitudes on the amount of fluid in the rotor cavity, the mass of the foundation, linear imbalance, the value of the stiffness coefficient, and the damping coefficient is shown. © 2021 Journal of Mechanical Engineering. 
Kydyrbekuly, A., Ibrayev, G.G.A., Ospan, T., Nikonov, A.,Multiparametric dynamic analysis of a rolling bearings system. //Strojniski Vestnik/Journal of Mechanical Engineering. – 2021. – Vol.67(9). – P.421432 DOI: 10.5545/svjme.2021.7178 
142. 
Modeling the Separation Process in Vertical Rotor Systems

DOI:10.1007/9783030835941_11 (SJR 0,159; Q4, Procentile 23(Mechanical Engineering) 18(Mechanics of Materials)) 
In this paper, we study and analyze the features of the separation process in a centrifugal force field, i.e. centrifugation process in vertical rotor systems. The main parameters that determine the time of separation of particles are revealed, and the optimal modes are indicated both from a constructive and from an economic point of view. Special cases of a fixed rotor are considered. Nonlinear differential equations of motion of a suspension particle are obtained, which do not have an exact solution. The study is carried out by analytical and numerical methods. The dependences of the slope angles of the tubes on the angular velocity of rotation of the rotor, sedimentation curves that allow one to estimate the time of deposition of particles, as well as the effect of the dispersed composition on the separation process as a whole, are obtained. The results of the study of this work allow us to determine with sufficient accuracy all the necessary characteristics working process of separation and sedimentation, and also allow in certain cases to exclude experimental work. © 2022, The Author(s), under exclusive license to Springer Nature Switzerland AG. 
Kydyrbekuly, A.B., Ibrayev, G.E.,Modeling the Separation Process in Vertical Rotor Systems. //Mechanisms and Machine Science. – 2021. – Vol.85. – P.104113 DOI: 10.1007/9783030835941_11 
143. 
Structuralparametric synthesis of the robomech class parallel mechanism with two sliders

DOI:10.3390/app11219831 (SJR 0,435; Q2, Procentile 71(General Engineering) 62(Instrumentation) 59(Fluid Flow and Transfer Processes) 56(Computer Science Applications) 51(General Materials Science) 50(Process Chemistry and Technology)) 
This paper addresses the structuralparametric synthesis and kinematic analysis of the RoboMech class of parallel mechanisms (PM) having two sliders. The proposed methods allow the synthesis of a PM with its structure and geometric parameters of the links to obtain the given laws of motions of the input and output links (sliders). The paper outlines a possible application of the proposed approach to design a PM for a cold stamping technological line. The proposed PM is formed by connecting two sliders (input and output objects) using one passive and one negative closing kinematic chain (CKC). The passive CKC does not impose a geometric constraint on the movements of the sliders and the geometric parameters of its links are varied to satisfy the geometric constraint of the negative CKC. The negative CKC imposes one geometric constraint on the movements of the sliders and its geometric parameters are determined on the basis of the Chebyshev and leastsquare approximations. Problems of positions and analogues of velocities and accelerations of the considered PM are solved to demonstrate the feasibility and effectiveness of the proposed formulations and case of study. © 2021 by the authors. Licensee MDPI, Basel, Switzerland. 
Baigunchekov, Z., Laribi, M.A., Carbone, G.,Mustafa A.,Amanov, B., Zholdassov, Y.,Structuralparametric synthesis of the robomech class parallel mechanism with two sliders. //Applied Sciences (Switzerland). – 2021. – Vol.11(21). – P.9831 DOI: 10.3390/app11219831 
144. 
Experimental investigation of an asphalt concrete deformation under cyclic loadingNews of the National Academy of Sciences of the Republic of Kazakhstan 
DOI: (SJR 0,323; Q3, Procentile 40(Geology) 37(Geotechnical Engineering and Engineering Geology) ) 
In the paper results of experimental investigation of an asphalt concrete deformation under cyclic loading on the scheme of direct tension are presented. Two series of the conventional hot dense asphalt concrete were tested. Durations of loading and rest periods were equal to 600 seconds. A stress has been applied to the tested sample quickly (within 1 second) and kept constant for the following 600 seconds. Then the stress was removed quickly and the sample was free of stress for the following 600 seconds. In the next cycle a stress greater than the preceding one was applied to the sample. In the first and second series 10 an 12 samples of the asphalt concrete were tested. In the first of them cyclic stresses were equal to 0.041; 0.111; 0.148; 0.183 and 0.219 MPa, in the second one stresses were equal to 0.744; 1.448; 2.232 and 2.976 MPa. Test temperature was 2224°C. The tests were carried out in an equipment constructed especially. The sample strain was measured by means of two clock typed indicators and was recorded in a video camera. Results of the tests showed that in each cycle maximum creep strain of an asphalt concrete occurs in the end of loading period and maximum recovered strain occurs in the end of relax period. Both these strains increase with the growth of cycle number, but the increase rate of the first one is higher than of the second one. In the considered test conditions (durations of loading and relax periods equal to 600 seconds, values of applied stresses from 0.041 MPa till 3.0 MPa, temperature T=2224°C) an asphalt concrete is a very plastic body: the total amount of plastic strain is equal to 8087% from the strain achieved in the end of loading period in each cycle. © 2018, National Academy of Sciences of the Republic of Kazakhstan. All rights reserved. 
Iskakbayev, A., Teltayev, B.B., Rossi, C.O., Yensebayeva, G.M..,Experimental investigation of an asphalt concrete deformation under cyclic loading News of the National Academy of Sciences of the Republic of Kazakhstan. //Series of Geology and Technical Sciences. – 2018. – Vol.2(428). – P.104111 DOI: 
145. 
A new simple damage accumulation model for predicting of an asphalt concrete cyclic strength

DOI:10.32014/2018.2518170X.8 (SJR 0,323; Q3, Procentile 40(Geology) 37(Geotechnical Engineering and Engineering Geology) ) 
Paper shows part of series test results for finegrained dense asphalt concrete, prepared with the use of oxidated bitumen to failure at cyclic creep. Load and relax duration for each cycle were equal to 10 and 60 s respectively. Test temperature: 22 °C. It has been determined that creep strain under load and recovered strain at relax increase with the growth of cyclic number. The value of the recovered strain in the end of relax period decreases (from 52 to 20%) till half of the cyclic number to failure according to exponential dependence with the growth of cyclic number and further it remains constant till sample failure. Using the curve of asphalt concrete longterm strength, constructed according to the test results of more than 110 samples to failure under creep scheme at the stresses from 0.05 to 0.31 MPa, based on Bailey’s criterion, the expression was obtained for determination of cyclic strength for the asphalt concrete. © National Academy of Sciences of the Republic of Kazakhstan, 2018. 
Iskakbayev, A.I., Teltayev, B.B., Oliviero Rossi, C., Estayev, K.,A new simple damage accumulation model for predicting of an asphalt concrete cyclic strength. //News of the National Academy of Sciences of the Republic of Kazakhstan, Series of Geology and Technical Sciences. – 2018. – Vol.5(431). – P.3847 DOI: 10.32014/2018.2518170X.8 
146. 
Computer modeling of creep for hereditary materials by abel’s kernel

DOI:10.32014/2018.2518170X.36 (SJR 0,323; Q3, Procentile 40(Geology) 37(Geotechnical Engineering and Engineering Geology) ) 
The work is devoted to the computer modeling of creep procedure for hereditary materials. Creep procedure is described by nonlinear integral equation of Yu. N. Rabotnov, and creep kernel is represented by Abel’s kernel. New efficient method has been proposed for determining of parameters (α, δ) for Abel’s kernel. Bisection method is used for obtaining of parameter α. Algorithm and relevant software have been developed for calculating of parameters α и δ. Efficient methods and relevant software have been developed for calculation of values of instantaneous strain and creep strain for hereditary materials. High accuracy of modeling of creep procedure has been shown by means of developed methods and software based on materials Nylon 6 and glassreinforced plastic TC 8/3250. © 2018, National Academy of Sciences of the Republic of Kazakhstan. All rights reserved. 
Iskakbayev, A.I., Teltayev, B.B., Yensebayeva, G.M., Kutimov, K.S.,Computer modeling of creep for hereditary materials by abel’s kernel . //News of the National Academy of Sciences of the Republic of Kazakhstan Series of Geology and Technical Sciences. – 2018. – Vol.6(432). – P.6676 DOI: 10.32014/2018.2518170X.36 
147. 
Determination of nonlinear creep parameters for hereditary materials

DOI:10.3390/app8050760 (SJR 0,435; Q2, Procentile 71(General Engineering) 62(Instrumentation) 59(Fluid Flow and Transfer Processes) 56(Computer Science Applications) 51(General Materials Science) 50(Process Chemistry and Technology)) 
This work proposes an effective algorithm for description of nonlinear deformation of hereditary materials based on Rabotnov's method of isochronous creep curves. The notions have been introduced for experimental and model rheological parameters and similarity coefficients of isochronous curves. It has been shown how using them, one can find instantaneous strains at various stress levels for description of nonlinear deformation of hereditary materials at creep. Relevant equations have been determined from the nonlinear integral equation of Yu. N. Rabotnov for the application cases of Rabotnov's fractional exponential kernel and Abel's kernel for nonlinear deformation of hereditary materials at creep. The improved methods have been given for determination of creep parameters α, ε0, δ, β, and λ. By processing and using test results for material Nylon 6 and glassreinforced plastic TC 8/3250, the process has been shown for sequential implementation of the developed methods for description of linear and nonlinear deformation of these materials at creep. From the results of the experimental investigation performed by the authors of this paper, it has been determined that finegrained, dense asphalt concrete at the temperature of 20 ± 2 °C and stresses up to 0.183 MPa at direct tension is deformed considerably in a nonlinear way. It has been shown in an illustrative way by construction of isochronous creep curves at various load durations and curves of experimental rheological parameter at various stresses. Nonlinear deformation of asphalt concrete at creep is adequately described by the proposed methods. © 2018 by the authors. 
Iskakbayev, A., Teltayev, B., Rossi, C.O., Yensebayeva, G.,Determination of nonlinear creep parameters for hereditary materials. //Applied Sciences (Switzerland). – 2018. – Vol.8(5). – P.760 DOI: 10.3390/app8050760 
148. 
Mechanical characteristics of asphalt concrete at different loading rates

DOI:10.32014/2019.2518170X.156 (SJR 0,323; Q3, Procentile 40(Geology) 37(Geotechnical Engineering and Engineering Geology) ) 
The results for experimental determination of characteristics have been given and analyzed in this article for deformation and failure of an asphalt concrete at two loading rates (0.058 and 0.652 MPa/s). Hot finegrained dense asphalt concrete of type B prepared with viscous bitumen of grade BND 100/130 has been selected which is traditionally used in road construction. The tests have been performed at the temperature of 2224°C in a specially invented and assembled device according to the scheme of direct tension. The asphalt concrete samples had a shape of rectangular beam with dimensions 5x5x15 cm. It is found that loading rate effects greatly the characteristics of deformation and failure of an asphalt concrete: failure time and specific work of deformation are decreased in 8 times approximately, failure strain is decreased in 1.5 times and as much as that the strength is increased (failure at tension) at the loading rate increase in 11 times from 0.058 MPa to 0.652 MPa/s. From the moment of loading to the moment of failure the asphalt concrete is deformed nonlinearly. The rate of nonlinearity is increased with the load increase. © 2019, National Academy of Sciences of the Republic of Kazakhstan. All rights reserved. 
Iskakbayev, A.I., Teltayev, B.B., Yensebayeva, G.M.,Mechanical characteristics of asphalt concrete at different loading rates. //News of the National Academy of Sciences of the Republic of Kazakhstan, Series of Geology and Technical Sciences. – 2019. – Vol.6(438). – P.6469 DOI: 10.32014/2019.2518170X.156

149. 
Impact of loading rate on asphalt concrete deformation and failure 
DOI:10.18720/MCE.100.8 (SJR 0,4; Q3, Procentile 81(Building and Construction) 79(Civil and Structural Engineering) ) 
The results for experimental determination of characteristics have been given and analyzed in this article for deformation and failure of an asphalt concrete at eleven loading rates from 0.000563 MPa/s to 0.652 MPa/s differing in 1158 times. A hot finegrained dense asphalt concrete of type B prepared with a viscous bitumen of grade BND 100/130 which is traditionally used in road construction has been selected for the research. The tests have been performed at the temperature of 22–24 °C in a specially invented and assembled device according to the scheme of direct tension. The asphalt concrete samples had a shape of rectangular beam with dimensions 5×5×15 cm. It is found that from the moment of loading to the moment of failure the asphalt concrete is deformed nonlinearly. The rate of nonlinearity is increased with the load increase. Loading rate effects greatly the characteristics of deformation and failure of the asphalt concrete: failure time, specific work of deformation and failure deformation are decreased in 242, 160 and 3 times respectively at the loading rate increase in 1158 (nearly 1200) times from 0.000563 MPa/s to 0.652 MPa/s and the strength is increased in 5 times. Dependences for characteristics of the asphalt concrete failure (failure time, failure deformation, specific work of deformation and strength) on a loading rate are described with a high accuracy by power functions. © Iskakbayev, A.,Teltayev, B.B.,Rossi, C.O.,Yensebayeva, G.,Abu, B., Kutimov, K.S., 2020. 
Iskakbayev, A., Teltayev, B., Rossi, C.O., Yensebayeva, G.Abu B.,Kutimov K.,,Impact of loading rate on asphalt concrete deformation and failure. //Magazine of Civil Engineering. – 2020. – Vol.100(8). – P.10008 DOI: 10.18720/MCE.100.8 
150. 
Modeling of hereditary materials relaxation by Abel Kernel 
DOI:10.32014/2020.2518170X.154 (SJR 0,323; Q3, Procentile 40(Geology) 37(Geotechnical Engineering and Engineering Geology) ) 

Zhurinov, M.Zh., Iskakbayev, A.I., Teltayev, B.B., Kutimov, K.S.,Modeling of hereditary materials relaxation by Abel Kernel. //News of the National Academy of Sciences of the Republic of Kazakhstan, Series of Geology and Technical Sciences. – 2020. – Vol.6(444). – P.254260 DOI: 10.32014/2020.2518170X.154 
151. 
Longterm strength of asphalt concrete and its applications2020, 244, 118325 
DOI:10.1016/j.conbuildmat.2020.118325 (SJR 1,662; Q1, Procentile 95(Civil and Structural Engineering) 94(Geotechnical Engineering and Engineering GeologyBuilding and Construction) 87(General Materials Science) ) 
An international journal dedicated to the investigation and innovative use of materials in construction and repair. Construction and Building Materials provides an international forum for the dissemination of innovative and original research and development in the field of construction and building materials and their application in new works and repair practice. The journal publishes a wide range of innovative research and application papers which describe laboratory and to a limited extent numerical investigations or report on full scale projects. Multipart papers are discouraged. Construction and Building Materials also publishes detailed case studies and some incisive review articles that contribute new understandings. We are focusing on construction materials papers and we exclude papers on structural engineering, geotechnics and unbound highway layers. The construction materials and technology covered include: cement, concrete reinforcement, bricks and mortars, additives, corrosion technology, ceramics, timber, steel, polymers, glass fibres, recycled materials, bamboo, rammed earth, nonconventional building materials, bituminous materials and railway material applications. The scope of Construction and Building Materials includes, but is not restricted to, materials, NDT and monitoring aspects of new works and repair and maintenance of the following: bridges, highrise buildings, dams, civil engineering structures, silos, highway pavements, tunnels, water containment structures, sewers, roofing, housing, coastal defences and railways. 
Iskakbayev, A.I., Teltayev, B.B., Yestayev, K.Z., Abu, B.D.,Longterm strength of asphalt concrete and its applications. //Construction and Building Materials. – 2020. – Vol.244. – P.118325 DOI: 10.1016/j.conbuildmat.2020.118325 
152. 
Research and developments on solar assisted compression heat pump systems – A comprehensive review (Part A: Modeling and modifications) 
DOI:10.1016/j.rser.2017.08.022 (SJR 3,522; Q1, Procentile 97(Renewable Energy, Sustainability and the Environment) ) 
This paper presents a comprehensive review on research and developments on solar assisted compression heat pump systems, mostly reported during the last two decades. The first part of this paper provides a detailed description of the past efforts on various system configurations, system modeling, enhancement of system performance, modifications in compression heat pump cycles and environment friendly refrigerant options for solar assisted compression heat pump systems. The economical and environmental impacts of solar assisted compression heat pump systems are also described. Further research needs in this important field of solar assisted compression heat pump are listed. The outcomes of this review confirm that, there is a lot of research scope for improving the performance of the system and reducing its initial cost to make it competitive in the global market. The information presented in this paper is shall be highly beneficial for the active researchers working presently on solar assisted compression heat pumps. © 2017 
Mohanraj, M., Belyayev, Y., Jayaraj, S., Kaltayev, A.,Research and developments on solar assisted compression heat pump systems – A comprehensive review (Part A: Modeling and modifications). //Renewable and Sustainable Energy Reviews,. – 2018. – Vol.83. – P.90123 DOI: 10.1016/j.rser.2017.08.022 
153. 
Research and developments on solar assisted compression heat pump systems – A comprehensive review (PartB: Applications) 
DOI:10.1016/j.rser.2017.08.086 (SJR 3,522; Q1, Procentile 97(Renewable Energy, Sustainability and the Environment)) 
The second part of this review presents a comprehensive review on applications of solar assisted compression heat pump systems and its associated research potentials. The reported investigations on solar assisted heat pumps are categorized into five major groups as follows: (a) drying, (b) room space heating; (c) agricultural green house space heating, (d) water heating and (e) desalination applications. Most of the reported investigations in open literature are reviewed. The study will then presents the major limitations of solar assisted compression heat pumps and further research tasks in various applications. This paper concludes that solar assisted compression heat pump system is the promising equipment used for heating applications. The information presented in the second part of the review is beneficial to the researchers, energy experts and decision makers in the field. © 2017 Elsevier Ltd 
Mohanraj, M., Belyayev, Y., Jayaraj, S., Kaltayev, A.,Research and developments on solar assisted compression heat pump systems – A comprehensive review (PartB: Applications). //Renewable and Sustainable Energy Reviews. – 2018. – Vol.83. – P.124155DOI: 10.1016/j.rser.2017.08.086 
154. 
Quantitative model of the formation mechanism of the rollfront uranium deposits 
DOI:10.18321/ectj724 (SJR 0,181; Q3, Procentile 35(General Chemical Engineering) 34(General Chemistry) 29(General Materials Science) 23(Condensed Matter Physics) ) 
The rollfront type deposits are crescent shaped accumulation of mineralization including uranium, selenium, molybdenum in reduced permeable sandstones. It generally forms within a geochemical barrier between mostly reduced and predominantly oxidized environments. Redox reactions between oxidant and reductant creates favorable conditions for uranium precipitation, while constant flow of oxidant continuously dissolves uranium minerals thereby creating a reactive transport. Several previous works had either focused on the characteristics of the rollfront type deposits, or on the description of chemical and geological processes involved in their genesis. Based on these previous works, authors aimed to mimic laboratory experiments numerically by reactive flow and numerical simulation. Data from one particular experiment was used to determine reaction rates between reactants to produce a model of reactive transport and chemical processes involved in the formation of rollfront type deposits. The resulting model was used to identify the causes of crescent like formations and to determine main mechanisms influencing rollfront evolution. A better understanding and simulation of the mechanism involved in the formation of rollfront type deposits and their properties would contribute to decreased exploration and production costs of commodities trapped within such accumulations. The results of this work can be used to model other deposits formed through infiltration and subsequent precipitation of various minerals at the redox interface. © 2018 AlFarabi Kazakh National University. 
Aizhulov, D.Y., Shayakhmetov, N.M., Kaltayev, A.,Quantitative model of the formation mechanism of the rollfront uranium deposits. //Eurasian ChemicoTechnological Journal. – 2018. – Vol.20(3). – P.213221DOI: 10.18321/ectj724 
155. 
A novel hybrid approach for insitu determining the thermal properties of subsurface layers around borehole heat exchanger 
DOI: 10.1016/j.ijheatmasstransfer.2018.05.107 (SJR 1,624; Q1, Procentile 98(Fluid Flow and Transfer Processes) 95(Mechanical Engineering) 94(Condensed Matter Physics) ) 
Performance of shallow geothermal systems such as borehole thermal energy storage (BTES) and ground source heat pump (GSHP) mainly depends on the thermal properties of the subsurface and proper design of borehole heat exchangers (BHE). This paper introduces a novel hybrid approach for measuring the effectiveness of BHEs and surrounding subsurface thermal properties, which combines traditional thermal response test (TRT) with the borehole temperature relaxation method (BTR), based on two dimensional radial conductive heat transfer. The new method allows for: (1) evaluation of how convective heat loss at groundwater layers influence estimation of subsurface thermal properties; (2) examination of nonuniform heat transfer through a BHE to stratified subsurface layers; and, (3) calculation of depthdependency of thermal properties of unsaturated subsurface layers. The hybrid approach was tested using a 50 m Utype BHE, the results of which indicated that convective heat transfer at the groundwater level altered the real value of effective thermal conductivity from 0.45 to 1.56 W/m K. The nonuniformity of heat transfer along the BHE was confirmed by calculations that showed subsurface thermal conductivities were depth dependent, varying between 0.34 and 0.61 W/m K. © 2018 Elsevier Ltd 
Akhmetov, B., Georgiev, A., Popov, R.,Turtaeva Z,.Kaltayev, A., Ding, Y,A novel hybrid approach for insitu determining the thermal properties of subsurface layers around borehole heat exchanger. //International Journal of Heat and Mass Transfer. – 2018. – Vol.126. – P.11381149 DOI: 10.1016/j.ijheatmasstransfer.2018.05.107 
156. 
A comprehensive study to evaluate absolute uptake of carbon dioxide adsorption onto composite adsorbent  Une étude complète visant à mesurer l'absorption absolue de dioxyde de carbone sur un adsorbant composite 
DOI:10.1016/j.ijrefrig.2019.01.014 (SJR 1,543; Q1, Procentile 90(Building and Construction) 88(Mechanical Engineering)) 
In this study, new composite adsorbent with enhanced thermal conductivity and adsorption capacity was synthesized and analyzed comprehensively for the development of compact CO 2 based adsorption cooling system. The consolidated composite was prepared employing activated carbon, graphene nanoplatelets and hydroxyl cellulose as a parent adsorbent, thermal conductivity enhancer, and binder, respectively. The surface area and pore volume of the composite were found to be 1778 ± 13 m 2 g −1 and 1.014 cm 3 g −1 , respectively. In addition, the composite showed 233% higher thermal conductivity compared to the parent activated carbon. Adsorption characteristics of CO 2 were measured at temperature ranging from 20 to 70 °C and pressures up to 5 MPa. Absolute uptake was evaluated from excess adsorption based on the following two methods: (i) the adsorbed phase volume is equal to the pore volume of the adsorbent; and (ii) the adsorbed phase volume is almost zero under low pressure and/or high temperature conditions. Furthermore, the averaging of above two methods was also taken for avoiding these two extreme assumptions. Obtained absolute adsorption uptake data were fitted with modified DubininAstakhov and Tóth models. Results indicated good approximation between data points and models. The average isosteric heats of adsorption estimated using modified DA and Tóth model were found to be 19.742 kJ mol −1 and 19.023 kJ mol −1 , respectively. The obtained characteristics of composite adsorbent are prerequisites for designing compact CO 2 based adsorption cooling systems. © 2019 Elsevier Ltd and IIR 
Berdenova, B., Pal, A., Muttakin, M.,Mitra S.,Thu K.,Saha, B.B., Kaltayev, ,A comprehensive study to evaluate absolute uptake of carbon dioxide adsorption onto composite adsorbent. //International Journal of Refrigeration. – 2019. – Vol.100. – P.131140 DOI: 10.1016/j.ijrefrig.2019.01.014 
157. 
Thermal Performance Simulation of a Heat Pump Assisted Solar Desalination System for Kazakhstan Climatic Conditions 
DOI:10.1080/01457632.2018.1451246 (SJR 0.818; Q2, Procentile 68(Fluid Flow and Transfer Processes) 67(Fluid Flow and Transfer Processes) 58(Condensed Matter Physics) ) 
A numerical model has been proposed in this research work for predicting the energy performance of a heat pump assisted regenerative solar still using phase change heat storage materials such as paraffin wax and hydrate salt. The numerical model is based on energy and mass balance. The performance was simulated for Fort Sherchenko town in Kazakhstan. The performance comparison between the conventional solar still and heat pump assisted regenerative solar still with and without phase change materials are presented. Moreover, the influence of two different phase change materials on thermal performance of a heat pump assisted regenerative solar still are described. The numerical simulation results showed that heat pump assisted regenerative solar still is more energy efficient and produces better productivity when compared to the conventional simple solar still. The use of hydrate salt phase change heat storage material was found to be a good option for enhancing the performance of a heat pump assisted regenerative solar still. © 2018, © 2018 Taylor & Francis Group, LLC. 
Belyayev, Y., Mohanraj, M., Jayaraj, S., Kaltayev, A.,Thermal Performance Simulation of a Heat Pump Assisted Solar Desalination System for Kazakhstan Climatic Conditions. //Heat Transfer Engineering. – 2019. – Vol.40(12). – P.10601072 DOI: 10.1080/01457632.2018.1451246 
158. 
Numerical simulation of a heat pump assisted solar dryer for continental climates 
DOI: 10.1016/j.renene.2019.04.119 (SJR 2,052; Q1, Procentile 88(Renewable Energy, Sustainability and the Environment)) 
A numerical model has been proposed in this work for predicting the energy performance of a heat pump assisted solar dryer under continental climates. The model is based on energy and mass balance. The energy performance comparison between heat pump dryer, solar dryer and heat pump assisted solar dryer are presented. The simulation was performed for four different climatic conditions of Almaty city in Kazakhstan. The simulation results showed that, the heat pump assisted solar dryer is more energy efficient when compared to the conventional solar dryers. It is also confirmed that, the conventional solar dryers are not suitable for continental climatic conditions having low ambient temperatures. The heat pump dryer reduces the initial moisture content of banana (on wet basis) from about 74% to the final moisture content (on wet basis) of about 19% in 21 h. Similarly, the solar dryer reduces the initial moisture content (wet basis) from about 74% to the final moisture content (wet basis) of about 20% in 35 h. The specific moisture extraction rate and coefficient of performance of a heat pump assisted solar dryer are estimated to be about 0.6 kg/kWh and 2.72, respectively. © 2019 Elsevier Ltd 
Kuan M.,Shakir Y.,Mohanraj M.,Belyayev Y. Jayaraj S.,Kaltayev A.,A.,Numerical simulation of a heat pump assisted solar dryer for continental climates. //Renewable Energy. – 2019. – Vol.143. – P.214225 DOI: 10.1016/j.renene.2019.04.119 
159. 
Numerical study of integrated latent heat thermal energy storage devices using nanoparticleenhanced phase change materials 
DOI: 10.1016/j.solener.2019.10.015 (SJR 1,537; Q1, Procentile 87(General Materials Science) 83(Renewable Energy, Sustainability and the Environment)) 
Two sequentially integrated LHTES devices based on paraffin waxes (PW), PWL and PWH with different phase change temperature ranges are numerically studied using Comsol Multiphysics for efficient thermal energy storage (TES). Thermal properties of the PWs are characterized and aluminum oxide nanoparticles (nanoAl2O3) are dispersed into the PWs to improve their heat transfer ability. According to the laser flash apparatus results, when the nanoAl2O3 composes 4 wt% of the mass of the PWL, its thermal diffusivity can be enhanced up to 40%. The same amount of the nanoAl2O3 improves the thermal diffusivity of the PWH approximately by 25%. Further characterization studies show that the incorporation of the nanoAl2O3 does not significantly change the specific heat capacity, latent heat of melting and cooling of the PCMs, but improves the heat transfer efficiency of the PCMs. Measured thermal properties of the PCMs are considered as input data in the numerical simulation of operating regimes of the devices. The full charging time of the integrated LHTES devices is reduced by 57 min and 106 min when the nanoAl2O3 composed 2 wt% and 4 wt% of the mass of the PCMs respectively. Likewise, the full discharging time of the integrated devices is decreased by 32 min and 74 min by the addition of the nanoAl2O3. Such reductions lead to improved charging and discharging efficiency of the LHTES devices. Moreover, simulation results show that the total amount of the stored energy in the devices fairly approximates the differential scanning calorimetry (DSC) results. © 2019 International Solar Energy Society 
Akhmetov, B., Navarro, M.E., Seitov, A.,Kaltayev A.,Bakenov, Z., Ding, Y.,Numerical study of integrated latent heat thermal energy storage devices using nanoparticleenhanced phase change materials. //Solar Energy. – 2019. – Vol.194. – P.724741 DOI: 10.1016/j.solener.2019.10.015 
160. 
Numerical simulation on solar collector and cascade heat pump combi water heating systems in Kazakhstan climates 
DOI: 10.1016/j.renene.2019.06.102 (SJR 2,052; Q1, Procentile 88(Renewable Energy, Sustainability and the Environment)) 
At low ambient temperatures, the heating capacity and coefficient of performance of a single stage vapour compression heat pump cycle is significantly getting reduced. The two stage cascade heat pump cycle operating with two different refrigerants provides a sustainable solution to lift the condenser temperature above 343 K. In this work, a numerical simulation model was developed for predicting the performance of a solar collector and two stage cascade heat pump combi water heating systems under Kazakhstan climatic conditions. The numerical simulation was performed for winter climatic conditions using nine refrigerant pairs such as, R32/R290, R32/R1234yf, R32/R134a, R410A/R290, R410A/R1234yf, R410A/R134a, R744/R290, R744/R1234yf and R744/R134a. The influences of solar irradiation, ambient temperature and condenser temperature are discussed. The solar collector and two stage heat pump combi heating system has about 30% energy savings when compared to the conventional two stage cascade heat pump without integration solar collectors. The R32/R290 refrigerant pair is having maximum coefficient of performance of 2.4 at 323 K condensing temperature and 263 K evaporating temperature. The refrigerant R744/R290 pair is identified as an environment friendly sustainable option in terms of its global warming impact for two stage cascade heat pump applications. © 2019 
Yerdesh, Y., Abdulina, Z., Aliuly, A.,Belyayev Y.,Mohanraj, M., Kaltayev, A.,Numerical simulation on solar collector and cascade heat pump combi water heating systems in Kazakhstan climates. //Renewable Energy. – 2020. – Vol.145. – P.12221234 DOI: 10.1016/j.renene.2019.06.102 
161. 
A GPUaccelerated Simulator of Turbulent Reacting Flows 
DOI: 10.1080/10618562.2020.1787996 (SJR 0,470; Q2, Procentile 60(Computational Mechanics) 57(Aerospace Engineering) 53(Mechanical Engineering) 51(Energy Engineering and Power Technology) 49(Mechanics of Materials) 43(Condensed Matter Physics)) 
A new computational methodology is developed for large eddy simulation (LES) of turbulent reacting flows using graphic processing units (GPUs). The LES is based on the filtered density function (FDF) of the scalarcomposition in conjunction with a discontinuous Galerkin (DG) discretisation scheme on a structured rectangular mesh. This hybrid solver is developed in a manner suitable for GPU computing. The simulator, as devised, is shown to be of the order of 200 times faster than the serial CPUbased calculations; facilitating the use of FDF for practical applications. The consistency and the accuracy of the methodology are demonstrated by simulations of a temporally developing mixing layer, under both nonreacting and reacting conditions. © 2020, © 2020 Informa UK Limited, trading as Taylor & Francis Group. 
Inkarbekov, M., Aitzhan, A., Kaltayev, A., Sammak, S.,A GPUaccelerated Simulator of Turbulent Reacting Flows. //International Journal of Computational Fluid Dynamics. – 2020. – P.381396 DOI: 10.1080/10618562.2020.1787996 
162. 
Control of elastic shear waves by periodic geometric transformation: cloaking, high reflectivity and anomalous resonancesa 
DOI: 10.1016/j.jmps.2019.103816 (SJR 1.857; Q1, Procentile 93(Mechanical Engineering)) 
A doubly periodic geometric transformation is applied to the problem of out of plane shear wave propagation in a doubly periodic perforated elastic medium. The technique leads to the design of a system of radially anisotropic and inhomogeneous shells surrounding the void inclusions, that can be tuned to give the desired filtering properties. For a regular transformation, the transformed elastic system displays the same dispersion properties than the original homogeneous one, but for overlapping and unfolding transformations new filtering properties can be obtained, which include anomalous resonances at zero and finite frequencies. Lowfrequency homogenisation reveals how it is possible to tune the phase and group velocity in the longwave limit at any value, or to obtain a zero frequency band gap for Neumann boundary conditions. The dispersion properties of the medium are studied both semianalytically by the multipole expansion and numerically by the finite element methods. Several applications are shown, including the transmission problem throughout a grating of void inclusions and an interface in a waveguide, where the capability of the proposed model is quantitatively demonstrated by computing the transmitted power flow. Finally, we gave a demonstration of defect modes in a waveguide and of tuning the transformation for the research of Dirac points. 
MEIRBEKOVA B., BRUN M. JOURNAL OF THE MECHANICS AND PHYSICS OF SOLIDSDOI: 10.1016/j.jmps.2019.103816

163. 
GMRES based numerical simulation and parallel implementation of multicomponent multiphase flow in porous media 
DOI: 10.1080/23311916.2020.1785189 (SJR 0.312; Q2, Procentile 67 (Engineering: General Engineering), 62 (Computer Science: General Computer Science), 52 (Chemical Engineering: General Chemical Engineering)) 
This article considered the numerical simulation of multicomponent multiphase flow in porous media. The resulting system of nonlinear equations linearized by the NewtonRaphson method and solved with the iterative Generalized minimal residual method (GMRES) algorithm. To achieve better convergence, we used the ILU(0) preconditioner to the GMRES algorithm. As a result, we used a completely implicit scheme called the NewtonILU0GMRES algorithm to solve the problem of interest. Based on the obtained sequential algorithm, we implemented a parallel algorithm using Message Passing Interface (MPI) technology. Additionally, we made comparisons between the parallel program of the presented algorithm and the parallel program using the readymade Portable Extensible Toolkit for Scientific Computation (PETSc) library. We developed an MPI parallel algorithm and tested it on the MVS10P supercomputer of the Interdepartmental Supercomputer Center of the Russian Academy of Sciences 
Mukhambetzhanov S.T. Lebedev D.V. Kassymbek N.M. Imankulov T.S. Matkerim B., AkhmedZaki D.Z.

164. 
Cogeneration Plants with Solar Radiation Concentrators 
DOI: 10.1134/S0040601520100079 (SJR 0.602; Q2, 43 (Engineering: General Engineering) 
Abstract: Results from experimental studies of a solar cogeneration system with linear photovoltaic modules of a fundamentally new design are presented. The Ʌshaped frontal walls are installed facetoface at an angle to each other and mutually shield their own thermal radiation, which decreases the radiation heat losses by 27% compared with linear photovoltaic modules of the known designs. The photocurrent generated by cooled solar cells is directed to a system for charging chemical batteries and the thermal energy released is transmitted to the unconsumed intermediate heattransfer fluid and then, through the surface of coil pipes of countercurrent heat exchangers, to the consumed process water of the outer circulation circuit. The further transportation of thermal energy to the storage system occurs by natural circulation of the consumed process water through the temperature gradient formed by the control system over the height between the heat source, the heat exchanger, and the heat receiver, an insulated container (a heat accumulator). For the first time, efficient controlled transportation of heat has been implemented without using a circulation pump owing to the excess thermal energy released during the conversion of solar energy by the solar cells and a photoselective film installed in the focal spot of the optical concentrator. Thus, a possibility of increasing the temperature of the heattransfer fluids at the cogeneration system outlet has been offered. A twocircuit circulation system allows for separation of unconsumed heattransfer fluids (antifreezing solutions) and the consumed fluid (the process water) by the pressure in the channels and installation of a linear countercurrent heat exchanger that performs the functions of a supporting platform’s mechanical axis along the rotational axis of the optical concentrator. 
Nesterenkov P.A. Nesterenkov A.G. Temirbekov A.N

165. 
Drying of Vegetable Products in Mobile Solar Dryer with Movable Shelving 
DOI: 10.1134/S1810232821010112 (SJR 0.457; Q2, 50 (Engineering: General Engineering) 
Abstract: The article outlines the relevance of production of dried products using a mobile solar dryer plant with mobile shelving. For intensification of the drying process, the convection of the drying agent flow in the proposed design of drying chamber was studied via numerical solution of the Navier–Stokes equations. As a result, a graphical interpretation of isolines of moving stream of drying agent in a chamber was obtained. Varying dimensionless parameters of the drying agent enabled complete coverage of all zones of the dryer chamber. The motion of shelves due to gravitational forces allowed uniform drying of materials. Through multivariate experiments, the influence of the temperaturetime regimes of the dryer chamber on the particle size was investigated. The optimal drying conditions for cut fruit are presented. 
Urmashev B.A. Khazimov K.M., Temirbekov A.N. Tursynbay A.T. Torzhenova T.V. Khazimov M.Z.
