Pub Date : 2024-06-01DOI: 10.1007/s10092-024-00579-z
Yanping Chen, Yujing Yan, Xiaoli Li, Xuan Zhao
{"title":"Second-order energy-stable scheme and superconvergence for the finite difference method on non-uniform grids for the viscous Cahn–Hilliard equation","authors":"Yanping Chen, Yujing Yan, Xiaoli Li, Xuan Zhao","doi":"10.1007/s10092-024-00579-z","DOIUrl":"https://doi.org/10.1007/s10092-024-00579-z","url":null,"abstract":"","PeriodicalId":9522,"journal":{"name":"Calcolo","volume":null,"pages":null},"PeriodicalIF":1.7,"publicationDate":"2024-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141409380","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}
Pub Date : 2024-06-01DOI: 10.1007/s10092-024-00584-2
Claudio I. Correa, G. Gatica, Esteban Henríquez, R. Ruiz-Baier, Manuel Solano
{"title":"Banach spaces-based mixed finite element methods for the coupled Navier–Stokes and Poisson–Nernst–Planck equations","authors":"Claudio I. Correa, G. Gatica, Esteban Henríquez, R. Ruiz-Baier, Manuel Solano","doi":"10.1007/s10092-024-00584-2","DOIUrl":"https://doi.org/10.1007/s10092-024-00584-2","url":null,"abstract":"","PeriodicalId":9522,"journal":{"name":"Calcolo","volume":null,"pages":null},"PeriodicalIF":1.7,"publicationDate":"2024-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141394450","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}
Pub Date : 2024-06-01DOI: 10.1007/s10092-024-00594-0
Haining Wen
{"title":"Delay-dependent stability of predictor–corrector methods of Runge–Kutta type for stochastic delay differential equations","authors":"Haining Wen","doi":"10.1007/s10092-024-00594-0","DOIUrl":"https://doi.org/10.1007/s10092-024-00594-0","url":null,"abstract":"","PeriodicalId":9522,"journal":{"name":"Calcolo","volume":null,"pages":null},"PeriodicalIF":1.7,"publicationDate":"2024-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141396330","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}
Pub Date : 2024-05-13DOI: 10.1007/s10092-024-00577-1
Yanjun Zhang, Hanyu Li, Ling Tang
The nonlinear Kaczmarz method was recently proposed to solve the system of nonlinear equations. In this paper, we first discuss two greedy selection rules, i.e., the maximum residual and maximum distance rules, for the nonlinear Kaczmarz iteration. Then, based on them, two kinds of greedy randomized sampling methods are presented. Furthermore, we also devise four corresponding greedy randomized block methods, i.e., the multiple samples-based methods. The linear convergence in expectation of all the proposed methods is proved. Numerical results show that, in some applications, including brown almost linear function and generalized linear model, the greedy selection rules give faster convergence rates than the existing ones, and the block methods outperform the single sample-based ones.
{"title":"Greedy randomized sampling nonlinear Kaczmarz methods","authors":"Yanjun Zhang, Hanyu Li, Ling Tang","doi":"10.1007/s10092-024-00577-1","DOIUrl":"https://doi.org/10.1007/s10092-024-00577-1","url":null,"abstract":"<p>The nonlinear Kaczmarz method was recently proposed to solve the system of nonlinear equations. In this paper, we first discuss two greedy selection rules, i.e., the maximum residual and maximum distance rules, for the nonlinear Kaczmarz iteration. Then, based on them, two kinds of greedy randomized sampling methods are presented. Furthermore, we also devise four corresponding greedy randomized block methods, i.e., the multiple samples-based methods. The linear convergence in expectation of all the proposed methods is proved. Numerical results show that, in some applications, including brown almost linear function and generalized linear model, the greedy selection rules give faster convergence rates than the existing ones, and the block methods outperform the single sample-based ones.</p>","PeriodicalId":9522,"journal":{"name":"Calcolo","volume":null,"pages":null},"PeriodicalIF":1.7,"publicationDate":"2024-05-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140925864","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}
Pub Date : 2024-04-25DOI: 10.1007/s10092-024-00574-4
Xavier Blanc, Philippe Hoch, Clément Lasuen
{"title":"An asymptotic preserving scheme for the $$M_1$$ model on polygonal and conical meshes","authors":"Xavier Blanc, Philippe Hoch, Clément Lasuen","doi":"10.1007/s10092-024-00574-4","DOIUrl":"https://doi.org/10.1007/s10092-024-00574-4","url":null,"abstract":"","PeriodicalId":9522,"journal":{"name":"Calcolo","volume":null,"pages":null},"PeriodicalIF":1.7,"publicationDate":"2024-04-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140654631","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}
Pub Date : 2024-04-06DOI: 10.1007/s10092-024-00575-3
Abstract
A constructive numerical approximation of the two-dimensional unsteady stochastic Navier–Stokes equations of an incompressible fluid is proposed via a pseudo-compressibility technique involving a penalty parameter (varepsilon ). Space and time are discretized through a finite element approximation and an Euler method. The convergence analysis of the suggested numerical scheme is investigated throughout this paper. It is based on a local monotonicity property permitting the convergence toward the unique strong solution of the stochastic Navier–Stokes equations to occur within the originally introduced probability space.
{"title":"Numerical approximation of the stochastic Navier–Stokes equations through artificial compressibility","authors":"","doi":"10.1007/s10092-024-00575-3","DOIUrl":"https://doi.org/10.1007/s10092-024-00575-3","url":null,"abstract":"<h3>Abstract</h3> <p>A constructive numerical approximation of the two-dimensional unsteady stochastic Navier–Stokes equations of an incompressible fluid is proposed via a pseudo-compressibility technique involving a penalty parameter <span> <span>(varepsilon )</span> </span>. Space and time are discretized through a finite element approximation and an Euler method. The convergence analysis of the suggested numerical scheme is investigated throughout this paper. It is based on a local monotonicity property permitting the convergence toward the unique strong solution of the stochastic Navier–Stokes equations to occur within the originally introduced probability space.</p>","PeriodicalId":9522,"journal":{"name":"Calcolo","volume":null,"pages":null},"PeriodicalIF":1.7,"publicationDate":"2024-04-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140587009","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}
Pub Date : 2024-04-05DOI: 10.1007/s10092-024-00570-8
Abstract
Recent developments in the field of the radial basis function-finite difference (RBF-FD) framework have been focused on conditionally positive definite polyharmonic splines (PHS). Within this context, our research focuses on deriving analytical weights for the RBF-FD+polynomials method within the framework of PHS. We provide convergence analyses for various stencils. To validate the accuracy of our derived formulations, we conduct a series of computational experiments across a range of test problems.
{"title":"Weight calculation and convergence analysis of polyharmonic spline (PHS) with polynomials for different stencils","authors":"","doi":"10.1007/s10092-024-00570-8","DOIUrl":"https://doi.org/10.1007/s10092-024-00570-8","url":null,"abstract":"<h3>Abstract</h3> <p>Recent developments in the field of the radial basis function-finite difference (RBF-FD) framework have been focused on conditionally positive definite polyharmonic splines (PHS). Within this context, our research focuses on deriving analytical weights for the RBF-FD+polynomials method within the framework of PHS. We provide convergence analyses for various stencils. To validate the accuracy of our derived formulations, we conduct a series of computational experiments across a range of test problems.</p>","PeriodicalId":9522,"journal":{"name":"Calcolo","volume":null,"pages":null},"PeriodicalIF":1.7,"publicationDate":"2024-04-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140587209","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}
The proposed work discusses discrete collocation and discrete Galerkin methods for second kind Fredholm–Hammerstein integral equations on half line ([0,infty )) using Kumar and Sloan technique. In addition, the finite section approximation method is applied to transform the domain of integration from ([0, infty )) to ([0,alpha ],~ alpha >0). In contrast to previous studies in which the optimal order of convergence is achieved for projection methods, we attained superconvergence rates in uniform norm using piecewise polynomial basis function. Moreover, these superconvergence rates are further enhanced by using discrete multi-projection (collocation and Galerkin) methods. In order to support the provided theoretical framework, numerical examples are included as well.
{"title":"Discrete projection methods for Fredholm–Hammerstein integral equations using Kumar and Sloan technique","authors":"Ritu Nigam, Nilofar Nahid, Samiran Chakraborty, Gnaneshwar Nelakanti","doi":"10.1007/s10092-024-00573-5","DOIUrl":"https://doi.org/10.1007/s10092-024-00573-5","url":null,"abstract":"<p>The proposed work discusses discrete collocation and discrete Galerkin methods for second kind Fredholm–Hammerstein integral equations on half line <span>([0,infty ))</span> using Kumar and Sloan technique. In addition, the finite section approximation method is applied to transform the domain of integration from <span>([0, infty ))</span> to <span>([0,alpha ],~ alpha >0)</span>. In contrast to previous studies in which the optimal order of convergence is achieved for projection methods, we attained superconvergence rates in uniform norm using piecewise polynomial basis function. Moreover, these superconvergence rates are further enhanced by using discrete multi-projection (collocation and Galerkin) methods. In order to support the provided theoretical framework, numerical examples are included as well.</p>","PeriodicalId":9522,"journal":{"name":"Calcolo","volume":null,"pages":null},"PeriodicalIF":1.7,"publicationDate":"2024-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140316041","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}
Pub Date : 2024-03-23DOI: 10.1007/s10092-024-00572-6
Abstract
In this article, band structure calculations of two dimensional (2D) anisotropic photonic-crystal fibers (PhCFs) are considered. In 2D PhCFs, Maxwell’s equations for the transversal electric and magnetic mode become decoupled, but the difficulty, arising from the anisotropic permittivity ({{varvec{varepsilon }}}) and/or permeability ({{varvec{mu }}},) plaguing the frequency-domain finite difference method, especially the original Yee’s scheme, is our top concern. To resolve this difficulty, we re-establish the connection between the lowest order finite element method with the quasi-periodic condition and Yee’s scheme using 2D non-orthogonal mesh, whereby the decoupled Maxwell’s equations in 2D anisotropic PhCFs are readily discretized into a generalized eigenvalue problem (GEP). Moreover, we spell out the nullspace of the resulting GEP, if it exists, and explicitly construct the Moore–Penrose pseudoinverse of the singular coefficient matrix, whose smallest positive eigenvalues can be solved by the inverse Lanczos method. Extensive band structures of 2D PhCFs are calculated and benchmarked against reliable results to demonstrate the accuracy and efficiency of our method.
{"title":"An efficient and unified method for band structure calculations of 2D anisotropic photonic-crystal fibers","authors":"","doi":"10.1007/s10092-024-00572-6","DOIUrl":"https://doi.org/10.1007/s10092-024-00572-6","url":null,"abstract":"<h3>Abstract</h3> <p>In this article, band structure calculations of two dimensional (2D) <em>anisotropic</em> photonic-crystal fibers (PhCFs) are considered. In 2D PhCFs, Maxwell’s equations for the transversal electric and magnetic mode become decoupled, but the difficulty, arising from the anisotropic permittivity <span> <span>({{varvec{varepsilon }}})</span> </span> and/or permeability <span> <span>({{varvec{mu }}},)</span> </span> plaguing the frequency-domain finite difference method, especially the original Yee’s scheme, is our top concern. To resolve this difficulty, we re-establish the connection between the lowest order finite element method with the quasi-periodic condition and Yee’s scheme using 2D <em>non-orthogonal</em> mesh, whereby the decoupled Maxwell’s equations in 2D anisotropic PhCFs are readily discretized into a generalized eigenvalue problem (GEP). Moreover, we spell out the nullspace of the resulting GEP, if it exists, and explicitly construct the Moore–Penrose pseudoinverse of the singular coefficient matrix, whose smallest positive eigenvalues can be solved by the inverse Lanczos method. Extensive band structures of 2D PhCFs are calculated and benchmarked against reliable results to demonstrate the accuracy and efficiency of our method.</p>","PeriodicalId":9522,"journal":{"name":"Calcolo","volume":null,"pages":null},"PeriodicalIF":1.7,"publicationDate":"2024-03-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140202234","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}