Mathematical Modelling and Analysis https://mma.vgtu.lt/index.php/MMA <p>Mathematical Modelling and Analysis publishes original research on all areas of mathematical modelling and analysis.&nbsp;<a href="https://journals.vilniustech.lt/index.php/MMA/about">More information ...</a></p> Vilnius Gediminas Technical University en-US Mathematical Modelling and Analysis 1392-6292 <p>Authors who publish with this journal agree to the following terms</p> <ul> <li class="show">that this article contains no violation of any existing copyright or other third party right or any material of a libelous, confidential, or otherwise unlawful nature, and that I will indemnify and keep indemnified the Editor and THE PUBLISHER against all claims and expenses (including legal costs and expenses) arising from any breach of this warranty and the other warranties on my behalf in this agreement;</li> <li class="show">that I have obtained permission for and acknowledged the source of any illustrations, diagrams or other material included in the article of which I am not the copyright owner.</li> <li class="show">on behalf of any co-authors, I agree to this work being published in the above named journal, Open Access, and licenced under a Creative Commons Licence, 4.0 <a href="https://creativecommons.org/licenses/by/4.0/legalcode">https://creativecommons.org/licenses/by/4.0/legalcode</a>. This licence allows for the fullest distribution and re-use of the work for the benefit of scholarly information.</li> </ul> <p>For authors that are not copyright owners in the work (for example government employees), please <a href="mailto:%20journals@vilniustech.lt">contact VILNIUS TECH</a>to make alternative agreements.</p> Modelling for a new mechanical passive damper and its mathematical analysis https://mma.vgtu.lt/index.php/MMA/article/view/18963 <p>To improve ride comfort, an oil damper should restrain its damping force in a high-velocity range. A lot of dampers with such properties were developed. For example, the one controlling the flux of oil by a leaf valve is widely adapted and reasonable. However, it is difficult to represent its dynamics with a simple mathematical model, and the cost of a computational fluid dynamics is too expensive. To overcome the disadvantages, the first author in [15] developed the other mechanical oil damper with sub-pistons instead of the leaf valve, which enabled us to propose a simple mathematical model with linear ordinary differential equations, thanks to the simple mechanism to control the oil flow. In this article, we give a more detailed mathematical model for the damper taking the dynamic pressure resistance into account, which is represented by nonlinear ordinary differential equations. In addition, a numerical scheme for the model is also proposed and its mathematical analysis such as the validity of the numerical solutions is shown.</p> Dai Watanabe Shuji Yoshikawa Copyright (c) 2024 The Author(s). Published by Vilnius Gediminas Technical University. http://creativecommons.org/licenses/by/4.0 2024-10-11 2024-10-11 29 4 600–620 600–620 10.3846/mma.2024.18963 A class of explicit second derivative general linear methods for non-stiff ODEs https://mma.vgtu.lt/index.php/MMA/article/view/19325 <p>In this paper, we construct explicit second derivative general linear methods (SGLMs) with quadratic stability and a large region of absolute stability for the numerical solution of non-stiff ODEs. The methods are constructed in two different cases: SGLMs with <em>p </em>= <em>q </em>= <em>r </em>= <em>s </em>and SGLMs with <em>p </em>= <em>q </em>and <em>r </em>= <em>s </em>= 2 in which <em>p</em>, <em>q</em>, <em>r </em>and <em>s </em>are respectively the order, stage order, the number of external stages and the number of internal stages. Examples of the methods up to order five are given. The efficiency of the constructed methods is illustrated by applying them to some well-known non-stiff problems and comparing the obtained results with those of general linear methods of the same order and stage order.</p> Mohammad Sharifi Ali Abdi Michal Braś Gholamreza Hojjati Copyright (c) 2024 The Author(s). Published by Vilnius Gediminas Technical University. http://creativecommons.org/licenses/by/4.0 2024-10-11 2024-10-11 29 4 621–640 621–640 10.3846/mma.2024.19325 Micropolar fluid-thin elastic structure interaction: variational analysis https://mma.vgtu.lt/index.php/MMA/article/view/20053 <p>We consider the non-stationary flow of a micropolar fluid in a thin channel with an impervious wall and an elastic stiff wall, motivated by applications to blood flows through arteries. We assume that the elastic wall is composed of several layers with different elastic characteristics and that the domains occupied by the two media are infinite in one direction and the problem is periodic in the same direction. We provide a complete variational analysis of the two dimensional interaction between the micropolar fluid and the stratified elastic layer. For a suitable data regularity, we prove the existence, the uniqueness and the regularity of the solution to the variational problem associated to the physical system. Increasing the data regularity, we prove that the fluid pressure is unique, we obtain additional regularity for all the unknown functions and we show that the solution to the variational problem is solution for the physical system, as well.</p> Grigory Panasenko Laetitia Paoli Ruxandra Stavre Copyright (c) 2024 The Author(s). Published by Vilnius Gediminas Technical University. http://creativecommons.org/licenses/by/4.0 2024-10-11 2024-10-11 29 4 641–668 641–668 10.3846/mma.2024.20053 A new approach for Solving a nonlinear system of second-order BVPs https://mma.vgtu.lt/index.php/MMA/article/view/19217 <p>In this paper, we introduce a new approach based on the Reproducing Kernel Method (RKM) for solving a nonlinear system of second-order Boundary Value Problems (BVPs) without the Gram-Schmidt orthogonalization process. What motivates us to use the RKM without the Gram-Schmidt orthogonalization process is its easy implementation, elimination of the Gram-Schmidt process, fewer calculations, and high accuracy. Finally, the compatibility of numerical results and theorems demonstrates that the Present method is effective.</p> Taher Amoozad Saeid Abbasbandy Tofigh Allahviranloo Mohsen Rostamy Malkhalifeh Copyright (c) 2024 The Author(s). Published by Vilnius Gediminas Technical University. http://creativecommons.org/licenses/by/4.0 2024-10-11 2024-10-11 29 4 669–683 669–683 10.3846/mma.2024.19217 The BKM criterion to the 3D double-diffusive magneto convection systems involving planar components https://mma.vgtu.lt/index.php/MMA/article/view/19674 <p>In this paper, we investigate the BKM type blowup criterion applied to 3D double-diffusive magneto convection systems. Specifically, we demonstrate that a unique local strong solution does not experience blow-up at time <em>T</em>, given that <img src="/public/site/images/irena/Picture12.png" width="173" height="17">). To prove this, we employ the logarithmic Sobolev inequality in the Besov spaces with negative indices and a well-known commutator estimate established by Kato and Ponce. This result is the further improvement and extension of the previous works by O (2021) and Wu (2023).</p> Chol-Jun O Fan Wu Copyright (c) 2024 The Author(s). Published by Vilnius Gediminas Technical University. http://creativecommons.org/licenses/by/4.0 2024-10-11 2024-10-11 29 4 684–693 684–693 10.3846/mma.2024.19674 Averaged reaction for nonlinear boundary conditions on a grill-type Winkler foundation https://mma.vgtu.lt/index.php/MMA/article/view/20137 <p>We consider a homogenization problem for the elasticity operator posed in a bounded domain of the half-space, a part of its boundary being in contact with the plane. This surface is traction-free out of “small regions”, where we impose nonlinear Winkler-Robin boundary conditions containing “large reaction parameters”. Non-periodical distribution of these regions is allowed provided that they have the same area. We show the convergence of solutions towards those of the homogenized problems depending on the relations between the parameters distance, sizes, and reaction.</p> Delfina Gómez María-Eugenia Pérez-Martínez Copyright (c) 2024 The Author(s). Published by Vilnius Gediminas Technical University. http://creativecommons.org/licenses/by/4.0 2024-11-22 2024-11-22 29 4 694–713 694–713 10.3846/mma.2024.20137 A two-derivative time integrator for the Cahn-Hilliard equation https://mma.vgtu.lt/index.php/MMA/article/view/20646 <p>This paper presents a two-derivative energy-stable method for the Cahn-Hilliard equation. We use a fully implicit time discretization with the addition of two stabilization terms to maintain the energy stability. As far as we know, this is the first time an energy-stable multiderivative method has been developed for phase-field models. We present numerical results of the novel method to support our mathematical analysis. In addition, we perform numerical experiments of two multiderivative predictor-corrector methods of fourth and sixth-order accuracy, and we show numerically that all the methods are energy stable.</p> Eleni Theodosiou Carina Bringedal Jochen Schütz Copyright (c) 2024 The Author(s). Published by Vilnius Gediminas Technical University. http://creativecommons.org/licenses/by/4.0 2024-11-22 2024-11-22 29 4 714–730 714–730 10.3846/mma.2024.20646 Mirroring and nonlinear perturbation of a circuit's system with multiple attractors https://mma.vgtu.lt/index.php/MMA/article/view/21033 <p>We infix the duality-symmetric and the mirror symmetry conversion processes into a dynamical system representing an electric circuit diagram with three input (or output) as shown in Figure 2. Hence, a new non-linear variable order initial value problem is obtained and solved using the Haar wavelet numerical method (HWNM). Error, stability and entropy analyzes show the reliability of the method. Numerical simulations are then implemented and show for the new system, existence of various attractors’ types (point attractors (PAs), limit cycles, strange attractors (SAs), double attractor (DA), coexisting attractors (CoAs)) with their mirror reflections. Both are in a symmetrical structure in which they face each other, separated by a changing symmetry line and exhibiting similar properties. The circuit implementation using a Field Programmable Gate Array (FPGA) is performed and concur with the expected results.</p> Emile Franc Doungmo Goufo Copyright (c) 2024 The Author(s). Published by Vilnius Gediminas Technical University. http://creativecommons.org/licenses/by/4.0 2024-11-29 2024-11-29 29 4 731–752 731–752 10.3846/mma.2024.21033 Regions of existence and uniqueness for singular two-point boundary value problems https://mma.vgtu.lt/index.php/MMA/article/view/18638 <p>A monotone iterative technique with lower and upper solutions is presented to identify the regions of existence for the solutions of singular two-point boundary value problems</p> <p><img src="/public/site/images/irena/Picture13.png" width="338" height="54"></p> <p>without requiring the monotonicity conditions on <em>f</em>(<em>x,y</em>). Under an additional condition on <em>f</em>(<em>x,y</em>), uniqueness of the solution is also established. These existence and uniqueness results are constructive and complement the existing results. Four examples including some engineering problems are given to illustrate the applicability of the proposed approach.</p> Shih-Hsiang Chang Copyright (c) 2024 The Author(s). Published by Vilnius Gediminas Technical University. http://creativecommons.org/licenses/by/4.0 2024-11-29 2024-11-29 29 4 753–766 753–766 10.3846/mma.2024.18638 The inverse problem of determining profiles of electrophysical parameters in eddy-current structuroscopy using apriori information on multifrequency probing https://mma.vgtu.lt/index.php/MMA/article/view/20022 <p>Based on the proposed methodology, the essence of which is to identify the profiles of electrophysical parameters of planar objects of eddy-current testing by means of surrogate optimization in the active PCA-space of reduced dimensionality, the effectiveness of the approach is proved by modeling the process of measurement control using apriori accumulated information about an object, in particular, multifrequency probing. The particularity of these studies is the consideration of previously collected information not only on profile variations, but also on the effect of various object probing frequencies on the signal of the surface probe. The functions of the storage device and information carrier were performed by a neural network metamodel, characterized by a high computational efficiency. Numerical experiments have determined the accuracy indicators of the proposed improved method for determining the distributions of magnetic permeability and electrical conductivity along the subsurface layer of a metal object with changes in a microstructure. The analysis of the modeling results indicates a significant reduction in the level of computational resources required to solve the problem and an increase in the accuracy of profile identification.</p> Volodymyr Ya. Halchenko Ruslana Trembovetska Volodymyr Tychkov Nataliia Tychkova Copyright (c) 2024 The Author(s). Published by Vilnius Gediminas Technical University. http://creativecommons.org/licenses/by/4.0 2024-11-29 2024-11-29 29 4 767–780 767–780 10.3846/mma.2024.20022 An effective numerical algorithm for coupled systems of Emden-Fowler equations via shifted airfoil functions of the first kind https://mma.vgtu.lt/index.php/MMA/article/view/19540 <p>The present paper deals with the designation of a new efficient numerical scheme for numerical solution of a class of coupled systems of Emden-Fowler equations describing several processes in applied sciences and technology. This method combines the shifted version of airfoil functions of the first kind (AFFK) and quasilinearization technique. More specifically, the quasilinearization technique is first applied to the original problem and the shifted AFFK (SAFFK) collocation matrix technique is then constructed for obtaining the solution of resulting family of submodels. We derive the error bound and analyze the convergence properties of the SAFFK. Computational and experimental simulations are carried out to describe the applicability and accuracy of the new technique. To show the benefit of the new approach, the computed numerical outcomes are compared with those results obtained via the Bernoulli and Haar wavelets collocation methods. It is evident from the numerical illustrations that the new developed scheme is superior to the existing available ones. The elapsed CPU time of the proposed method is provided.</p> Mohammad Izadi Pradip Roul Copyright (c) 2024 The Author(s). Published by Vilnius Gediminas Technical University. http://creativecommons.org/licenses/by/4.0 2024-11-29 2024-11-29 29 4 781–800 781–800 10.3846/mma.2024.19540