Computers & Mathematics with Applications最新文献

{"title":"Integral transform technique for determining stress intensity factor in wave propagation through functionally graded piezoelectric-viscoelastic structure","authors":"Diksha , Soniya Chaudhary , Pawan Kumar Sharma , Qasem M. Al-Mdallal","doi":"10.1016/j.camwa.2025.03.021","DOIUrl":"10.1016/j.camwa.2025.03.021","url":null,"abstract":"<div><div>This study employs an integral transform approach for Love wave propagation in a rotating composite structure having an interfacial crack. The structure comprises an initially stressed functionally graded piezoelectric-viscoelastic half-space bonded to a piezoelectric-viscoelastic half-space, and is subjected to anti-plane mechanical loading and in-plane electrical loading. The study focuses on two material systems: the first material system consists of Epoxy-BNKLBT and Epoxy-KNLNTS, where BNKLBT stands for <span><math><mn>0.885</mn><mo>(</mo><msub><mrow><mtext>Bi</mtext></mrow><mrow><mn>0.5</mn></mrow></msub><msub><mrow><mtext>Na</mtext></mrow><mrow><mn>0.5</mn></mrow></msub><mo>)</mo><msub><mrow><mtext>TiO</mtext></mrow><mrow><mn>3</mn></mrow></msub><mo>−</mo><mn>0.05</mn><mo>(</mo><msub><mrow><mtext>Bi</mtext></mrow><mrow><mn>0.5</mn></mrow></msub><msub><mrow><mtext>K</mtext></mrow><mrow><mn>0.5</mn></mrow></msub><mo>)</mo><msub><mrow><mtext>TiO</mtext></mrow><mrow><mn>3</mn></mrow></msub><mo>−</mo><mn>0.015</mn><mo>(</mo><msub><mrow><mtext>Bi</mtext></mrow><mrow><mn>0.5</mn></mrow></msub><msub><mrow><mtext>Li</mtext></mrow><mrow><mn>0.5</mn></mrow></msub><mo>)</mo><msub><mrow><mtext>TiO</mtext></mrow><mrow><mn>3</mn></mrow></msub><mo>−</mo><mn>0.05</mn><msub><mrow><mtext>BaTiO</mtext></mrow><mrow><mn>3</mn></mrow></msub></math></span>, and KNLNTS represents <span><math><mo>(</mo><msub><mrow><mtext>K</mtext></mrow><mrow><mn>0.475</mn></mrow></msub><msub><mrow><mtext>Na</mtext></mrow><mrow><mn>0.475</mn></mrow></msub><msub><mrow><mtext>Li</mtext></mrow><mrow><mn>0.05</mn></mrow></msub><mo>)</mo><mo>(</mo><msub><mrow><mtext>Nb</mtext></mrow><mrow><mn>0.92</mn></mrow></msub><msub><mrow><mtext>Ta</mtext></mrow><mrow><mn>0.05</mn></mrow></msub><msub><mrow><mtext>Sb</mtext></mrow><mrow><mn>0.03</mn></mrow></msub><mo>)</mo><msub><mrow><mtext>O</mtext></mrow><mrow><mn>3</mn></mrow></msub></math></span>, doped with <span><math><mn>0.4</mn><mspace></mspace><mtext>wt%</mtext><mspace></mspace><msub><mrow><mtext>CeO</mtext></mrow><mrow><mn>2</mn></mrow></msub></math></span> and <span><math><mn>0.4</mn><mspace></mspace><mtext>wt%</mtext><mspace></mspace><msub><mrow><mtext>MnO</mtext></mrow><mrow><mn>2</mn></mrow></msub></math></span>. The second material system has Epoxy-BNKLBT and Epoxy-PZT7A, where PZT7A denotes Lead Zirconate Titanate. The viscoelastic materials are modeled to reflect their complex behavior under rotational and stress conditions. The Galilean transformation is applied to convert the Cartesian coordinate system into a moving reference frame aligned with the Love wave's propagation. Employing Bessel function properties, the system is converted into a set of double integral equations and subsequently reformulated into simultaneous Fredholm integral equations. Numerical solutions to these Fredholm integral equations are employed to compute the electric displacement intensity factor and the stress intensity factor near the interfaci","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"186 ","pages":"Pages 130-154"},"PeriodicalIF":2.9,"publicationDate":"2025-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143706299","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Semi-implicit Lax-Wendroff kinetic scheme for multi-scale phonon transport","authors":"Shuang Peng ,&nbsp;Songze Chen ,&nbsp;Hong Liang ,&nbsp;Chuang Zhang","doi":"10.1016/j.camwa.2025.03.019","DOIUrl":"10.1016/j.camwa.2025.03.019","url":null,"abstract":"<div><div>Fast and accurate predictions of the spatiotemporal distributions of temperature are crucial to the multi-scale thermal management and safe operation of microelectronic devices. To realize it, an efficient semi-implicit Lax-Wendroff kinetic scheme is developed for numerically solving the transient phonon Boltzmann transport equation (BTE) from the ballistic to diffusive regime. The biggest innovation of the present scheme is that the finite difference method is used to solve the phonon BTE for the reconstruction of the interfacial distribution function at the half-time step, where the second-order numerical schemes are used for both the temporal and spatial discretization. Consequently, the phonon scattering and migration are coupled together within one time step, and the evolution process of phonon distribution function follows the actual physical law even if the time step is much longer than the relaxation time. Numerical results show that the present scheme could accurately predict the steady/unsteady heat conduction in solid materials from the ballistic to diffusive regime, and its time step or cell size is not limited by the relaxation time or phonon mean free path. The present work could provide a useful tool for the efficient predictions of the macroscopic spatiotemporal distributions in the multi-scale thermal engineering.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"187 ","pages":"Pages 72-84"},"PeriodicalIF":2.9,"publicationDate":"2025-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143706248","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A non-convex and non-smooth weighted image denoising model","authors":"Huayu Fan ,&nbsp;Qiqi Feng ,&nbsp;Rui Chen ,&nbsp;Xiangyang Cao ,&nbsp;Zhi-Feng Pang","doi":"10.1016/j.camwa.2025.03.010","DOIUrl":"10.1016/j.camwa.2025.03.010","url":null,"abstract":"<div><div>In order to provide a more effective method to describe the local structure of the degraded image and to enhance the robustness of the denoising, we propose a non-convex total variational image denoising model that combines the non-convex log function with an adaptive weighted matrix within the total variation framework. In the proposed model, the weighted matrix is capable of effectively describing the primary direction of the edge structure, based on the coupling of the gradient operator of the denoising image and the diagonal matrix. As the proposed model is a non-convex and non-smooth optimisation problem, the iterative reweighted <span><math><msup><mrow><mi>ℓ</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> algorithm and alternating direction multiplier method are employed to decompose it into a number of readily solvable sub-problems. The results obtained from numerical experiments demonstrate that the proposed model is capable of effectively suppressing the noise while maintaining the local structure of the image.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"187 ","pages":"Pages 85-105"},"PeriodicalIF":2.9,"publicationDate":"2025-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143706252","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"H1− Galerkin mixed finite element method using tensor product of cubic B-splines for two-dimensional Kuramoto-Sivashinsky equation","authors":"L. Jones Tarcius Doss,&nbsp;V. Sindhujarani","doi":"10.1016/j.camwa.2025.03.009","DOIUrl":"10.1016/j.camwa.2025.03.009","url":null,"abstract":"<div><div>The two-dimensional <span><math><mo>(</mo><mn>2</mn><mi>D</mi><mo>)</mo></math></span> Kuramoto-Sivashinsky equation offers a robust framework for studying complex, chaotic, and nonlinear dynamics in various mathematical and physical contexts. Analyzing this model also provides insights into higher-dimensional spatio-temporal chaotic systems that are relevant to many fields. This article aims to solve the scalar form of the two-dimensional Kuramoto-Sivashinsky equation using the <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup><mo>−</mo></math></span> mixed Galerkin finite element method. By introducing an intermediate variable, the equation is transformed into a coupled system. This system is then approximated using the <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup><mo>−</mo></math></span> mixed Galerkin finite element method, with the tensor product of the cubic B-spline in <em>x</em> and <em>y</em> directions serving as the test and trial functions in both the semi-discrete and fully discrete schemes. In this approach, triangularization is avoided, thereby reducing the size of the stiffness matrix. In the fully discrete scheme, the time derivative is approximated using the backward Euler method. The suitable linearization method is used to simplify the nonlinear term in both schemes. The theoretical analysis yields optimal order error estimates for the scalar unknown and its flux in the <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span>, <span><math><msup><mrow><mi>L</mi></mrow><mrow><mo>∞</mo></mrow></msup></math></span>, and <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> norms, demonstrating the accuracy, efficiency, and stability of the proposed method. Additionally, three test problems are numerically analyzed to validate these theoretical results. The chaotic behavior of the equation is analyzed, in relation to the viscosity coefficient <em>γ</em>, and is also numerically investigated using the proposed method.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"188 ","pages":"Pages 19-39"},"PeriodicalIF":2.9,"publicationDate":"2025-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143687659","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Phase field lattice Boltzmann method for liquid-gas flows in complex geometries with efficient and consistent wetting boundary treatment","authors":"Dmytro Sashko ,&nbsp;Travis R. Mitchell ,&nbsp;Łukasz Łaniewski-Wołłk ,&nbsp;Christopher R. Leonardi","doi":"10.1016/j.camwa.2025.03.014","DOIUrl":"10.1016/j.camwa.2025.03.014","url":null,"abstract":"<div><div>This study investigates the application of wetting boundary conditions for modelling flows in complex curved geometries, such as rough fractures. It implements and analyses two common variants of the wetting boundary condition within the three-dimensional (3D) phase field lattice Boltzmann method. It provides a straightforward and novel extension of the geometrical approach to curved three-dimensional surfaces. It additionally implements surface-energy approach. A novel interpolation-based mitigation of the staircase approximation for curved boundaries is then developed and consistently applied to both wetting boundary conditions. The objectives of simplicity and parallel compute efficiency in implementation are emphasised. Through detailed validation on a series of 3D benchmark cases involving curved surfaces, such as droplet spread on a sphere, capillary intrusion, and droplet impact on a sphere, the behaviour of the wetting boundary conditions are validated and the differences between methods are highlighted. To demonstrate the applicability of the proposed approach in complex geometries with varying surface curvatures, two-phase flow through a synthetic rough fracture is presented. The suitability of the methods for complex simulations is also verified by comparing the computational performance between all investigated methods using this fracture flow test case. The present work thus contributes to the field of multiphase flow modelling with the lattice Boltzmann method in realistic applications where addressing the impact of complex geometries is essential.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"186 ","pages":"Pages 101-129"},"PeriodicalIF":2.9,"publicationDate":"2025-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143696762","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Mixed spectral element method combined with second-order time stepping schemes for a two-dimensional nonlinear fourth-order fractional diffusion equation","authors":"Jiarui Wang,&nbsp;Yining Yang,&nbsp;Hong Li,&nbsp;Yang Liu","doi":"10.1016/j.camwa.2025.03.015","DOIUrl":"10.1016/j.camwa.2025.03.015","url":null,"abstract":"<div><div>In this article, a mixed spectral element method combined with second-order time stepping schemes for solving a two-dimensional nonlinear fourth-order fractional diffusion equation is constructed. For formulating an efficient numerical scheme, an auxiliary function is introduced to transform the fourth-order fractional system into a low-order coupled system, then the time direction is discretized by second-order FBT-<em>θ</em> schemes, and the spatial direction is approximated using the Legendre mixed spectral element method (LMSEM). The stability and the optimal error estimate with <span><math><mi>O</mi><mo>(</mo><msup><mrow><mi>τ</mi></mrow><mrow><mn>2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>h</mi></mrow><mrow><mi>min</mi><mo>⁡</mo><mo>{</mo><mi>N</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>r</mi><mo>}</mo></mrow></msup><msup><mrow><mi>N</mi></mrow><mrow><mo>−</mo><mi>r</mi></mrow></msup><mo>)</mo></math></span> for the fully discrete scheme are derived, where <em>τ</em> stands for the time step size, <em>h</em> denotes the space step size, <em>N</em> indicates the degree of the polynomial, and <em>r</em> represents the order of Sobolev space. Finally, some numerical tests are carried out to verify the theory results and the effectiveness of the developed algorithm.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"188 ","pages":"Pages 1-18"},"PeriodicalIF":2.9,"publicationDate":"2025-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143687658","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Design and time-domain finite element analysis of a carpet thermal concentrator in metamaterials","authors":"Bin He ,&nbsp;Shouzhu Bao","doi":"10.1016/j.camwa.2025.03.016","DOIUrl":"10.1016/j.camwa.2025.03.016","url":null,"abstract":"<div><div>Traditional transform thermodynamic devices are designed from anisotropic materials which are difficult to fabricate. In this paper, we design and simulate a carpet thermal concentrator. Based on existing transformation thermodynamic techniques, we have derived the perfect parameters required for carpet heat concentrators. In order to eliminate the anisotropy of perfect parameters, we designed a heat concentration device for isotropic materials using the effective medium theory, and a posterior error analysis of the finite element discretization scheme in the metamaterial region is provided. Finally, we present the numerical simulation results to verify the correctness of the analysis and the performance of the designed heat concentration device.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"185 ","pages":"Pages 94-109"},"PeriodicalIF":2.9,"publicationDate":"2025-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143696089","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Fully consistent lowest-order finite element methods for generalised Stokes flows with variable viscosity","authors":"Felipe Galarce ,&nbsp;Douglas R.Q. Pacheco","doi":"10.1016/j.camwa.2025.03.013","DOIUrl":"10.1016/j.camwa.2025.03.013","url":null,"abstract":"<div><div>In finite element methods for incompressible flows, the most popular approach to allow equal-order velocity-pressure pairs are residual-based stabilisations. When using first-order elements, however, the viscous part of the residual cannot be approximated, which often degrades accuracy. For constant viscosity, or by assuming a Lipschitz condition on the viscosity field, we can construct stabilisation methods that fully approximate the residual, regardless of the polynomial order of the finite element spaces. This work analyses and tests two variants of such a fully consistent approach, with the generalised Stokes system as a model problem. We prove unique solvability and derive expressions for the stabilisation parameter, generalising some classical results for constant viscosity. Numerical results illustrate how our method completely eliminates the spurious pressure boundary layers typically induced by low-order PSPG-like stabilisations.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"188 ","pages":"Pages 40-49"},"PeriodicalIF":2.9,"publicationDate":"2025-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143687660","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A lattice-Boltzmann inspired finite volume solver for compressible flows","authors":"Jinhua Lu,&nbsp;Song Zhao,&nbsp;Pierre Boivin","doi":"10.1016/j.camwa.2025.03.007","DOIUrl":"10.1016/j.camwa.2025.03.007","url":null,"abstract":"<div><div>The lattice Boltzmann method (LBM) for compressible flow is characterized by good numerical stability and low dissipation, while the conventional finite volume solvers have intrinsic conversation and flexibility in using unstructured meshes for complex geometries. This paper proposes a strategy to combine the advantages of the two kinds of solvers by designing a finite volume solver to mimic the LBM algorithm. It assumes an ideal LBM that can recover all desired higher-order moments. Time-discretized moment equations with second-order temporal accuracy and physically consistent dissipation terms are derived from the ideal LBM. By solving the recovered moment equations, a finite volume solver that can be applied to nonuniform meshes naturally, enabling body-fitted mass-conserving simulations, is proposed. Numerical tests show that the proposed solver can achieve good numerical stability from subsonic to hypersonic flows, and low dissipation for a long-distance entropy spot convection. For the challenging direct simulations of acoustic waves, its dissipation can be significantly reduced compared with the Lax-Wendroff solver of the same second-order spatial and temporal accuracy, while only remaining higher than that of the LBM on coarse meshes. The analysis implies that approximations of third-order temporal accuracy are required to recover the low dissipation of LBM further.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"187 ","pages":"Pages 50-71"},"PeriodicalIF":2.9,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143666522","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A decoupled, convergent and fully linear algorithm for the Landau–Lifshitz–Gilbert equation with magnetoelastic effects","authors":"Hywel Normington ,&nbsp;Michele Ruggeri","doi":"10.1016/j.camwa.2025.03.008","DOIUrl":"10.1016/j.camwa.2025.03.008","url":null,"abstract":"<div><div>We consider the coupled system of the Landau–Lifshitz–Gilbert equation and the conservation of linear momentum law to describe magnetic processes in ferromagnetic materials including magnetoelastic effects in the small-strain regime. For this nonlinear system of time-dependent partial differential equations, we present a decoupled integrator based on first-order finite elements in space and an implicit one-step method in time. We prove unconditional convergence of the sequence of discrete approximations towards a weak solution of the system as the mesh size and the time-step size go to zero. Compared to previous numerical works on this problem, for our method, we prove a discrete energy law that mimics that of the continuous problem and, passing to the limit, yields an energy inequality satisfied by weak solutions. Moreover, our method does not employ a nodal projection to impose the unit length constraint on the discrete magnetisation, so that the stability of the method does not require weakly acute meshes. Furthermore, our integrator and its analysis hold for a more general setting, including body forces and traction, as well as a more general representation of the magnetostrain. Numerical experiments underpin the theory and showcase the applicability of the scheme for the simulation of the dynamical processes involving magnetoelastic materials at submicrometer length scales.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"187 ","pages":"Pages 1-29"},"PeriodicalIF":2.9,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143666523","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}