Journal of Computational and Applied Mathematics最新文献

{"title":"Reconstructing algebraic functions from a nonconforming exponential weighted enriched finite element","authors":"Francesco Dell’Accio ,&nbsp;Allal Guessab ,&nbsp;Gradimir V. Milovanović ,&nbsp;Federico Nudo","doi":"10.1016/j.cam.2025.116603","DOIUrl":"10.1016/j.cam.2025.116603","url":null,"abstract":"<div><div>The reconstruction of functions is a fundamental task in various applications, ranging from computer graphics to remote sensing. This paper addresses the challenge of function reconstruction in scenarios where, instead of pointwise function evaluations, only a set of integrals over specific lines is available. We propose a novel method based on a one-parameter family of weighted finite elements that incorporates exponential Hermite weight functions within bounded domains. Numerical experiments demonstrate that the proposed approach significantly improves computational efficiency and accuracy in function reconstruction compared to the classical Crouzeix–Raviart finite element.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"466 ","pages":"Article 116603"},"PeriodicalIF":2.1,"publicationDate":"2025-02-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143527179","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 method for solving ill-conditioned separable nonlinear least squares problems and its application","authors":"Luyao Wang ,&nbsp;Guolin Liu ,&nbsp;Yang Chen ,&nbsp;Huadong Ma","doi":"10.1016/j.cam.2025.116624","DOIUrl":"10.1016/j.cam.2025.116624","url":null,"abstract":"<div><div>For ill-conditioned separable nonlinear least squares problems, the LM (Levenberg-Marquardt) iteration method based on the VP (Variable Projection) algorithm and SVD (Singular Value Decomposition) for solving nonlinear parameters is explored in this paper, and a ridge estimation is proposed based on the TSVD (Truncated Singular Value Decomposition) and MSVD (Modified Singular Value Decomposition) methods to solve linear parameters. It is proved that the LM method based on SVD and the improved algorithm based on TSVD and MSVD can overcome the ill-conditioning of the matrix to a certain extent and enhance the reliability of parameter estimation through Mackey-Glass time series simulation and height anomaly fitting experiments. TSVD, MSVD and improved algorithm are compared and analyzed in terms of the accuracy and stability of parameter estimation and the curve fitting effect in the Mackey-Glass time series simulation experiment. The height anomaly fitting experiment verifies the feasibility and applicability of the algorithm in such practical problems. The experimental results show that the improved algorithm based on TSVD and MSVD enables the parameter estimation to be more stable, and the parameter values calculated by the algorithm bring about a better goodness of fit and fitting effect of the model.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"466 ","pages":"Article 116624"},"PeriodicalIF":2.1,"publicationDate":"2025-02-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143551037","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 new double relaxed inertial viscosity-type algorithm for split equilibrium problems and its application to detecting osteoporosis health problems","authors":"Watcharaporn Yajai ,&nbsp;Wongthawat Liawrungrueang ,&nbsp;Watcharaporn Cholamjiak","doi":"10.1016/j.cam.2025.116602","DOIUrl":"10.1016/j.cam.2025.116602","url":null,"abstract":"<div><div>This article proposes a new double relaxed inertial viscosity-type algorithm for solving split equilibrium problems; this algorithm is flexible to use by the relaxed extrapolation parameters in real numbers. We obtain a strong convergence theorem under suitable assumptions of two bifunctions in real Hilbert spaces. In data classification, we apply our proposed algorithm for finding optimal output weight in an extreme learning machine. The primary features of the Osteoporosis dataset from Harvard Dataverse are used for the algorithm’s experiments. The right extrapolation parameters show that the algorithm converges faster than the standard algorithm. The comparison with the existing algorithm is demonstrated for our proposed algorithm’s high efficiency. Finally, our algorithm’s accuracy and loss plots are presented, obtaining our good-fitting model.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"466 ","pages":"Article 116602"},"PeriodicalIF":2.1,"publicationDate":"2025-02-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143519490","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 compact scheme for the Munk boundary-layer equation in one dimension","authors":"M. Ben-Artzi ,&nbsp;J.-P. Croisille ,&nbsp;D. Fishelov","doi":"10.1016/j.cam.2025.116595","DOIUrl":"10.1016/j.cam.2025.116595","url":null,"abstract":"<div><div>In this paper, we introduce a two-scale compact finite difference scheme for the equation <span><span><span>(MK-1D)</span><span><math><mfenced><mrow><mtable><mtr><mtd></mtd><mtd><mo>−</mo><mi>β</mi><mfrac><mrow><mi>d</mi></mrow><mrow><mi>d</mi><mi>x</mi></mrow></mfrac><mi>u</mi><mo>+</mo><mi>ɛ</mi><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mi>d</mi></mrow><mrow><mi>d</mi><mi>x</mi></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mrow><mn>4</mn></mrow></msup><mi>u</mi><mo>=</mo><mi>f</mi><mo>,</mo><mspace></mspace><mi>x</mi><mo>∈</mo><mrow><mo>(</mo><mi>a</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd></mtd><mtd><mi>u</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><mi>u</mi><mrow><mo>(</mo><mi>b</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mi>u</mi></mrow><mrow><mo>′</mo></mrow></msup><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mi>u</mi></mrow><mrow><mo>′</mo></mrow></msup><mrow><mo>(</mo><mi>b</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>.</mo></mtd></mtr></mtable></mrow></mfenced></math></span></span></span>This equation serves as a model for the nonlinear barotropic equation (NB) governing oceanic flows. <span><span><span>(NB)</span><span><math><mrow><msub><mrow><mi>∂</mi></mrow><mrow><mi>t</mi></mrow></msub><mi>Δ</mi><mi>ψ</mi><mo>+</mo><msup><mrow><mo>∇</mo></mrow><mrow><mo>⊥</mo></mrow></msup><mi>ψ</mi><mo>.</mo><mo>∇</mo><mi>Δ</mi><mi>ψ</mi><mo>+</mo><mi>β</mi><msub><mrow><mi>∂</mi></mrow><mrow><mi>x</mi></mrow></msub><mi>ψ</mi><mo>=</mo><mfrac><mrow><mn>1</mn></mrow><mrow><mi>H</mi></mrow></mfrac><msub><mrow><mrow><mo>(</mo><mo>∇</mo><mo>×</mo><mi>τ</mi><mo>)</mo></mrow></mrow><mrow><mi>v</mi></mrow></msub><mo>−</mo><mi>μ</mi><mi>Δ</mi><mi>ψ</mi><mo>+</mo><mi>ɛ</mi><msup><mrow><mi>Δ</mi></mrow><mrow><mn>2</mn></mrow></msup><mi>ψ</mi><mo>,</mo></mrow></math></span></span></span>where <span><math><mrow><mi>ψ</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow></math></span> and <span><math><mi>τ</mi></math></span> are the streamfunction and the wind stress tensor, respectively. This equation encodes the <em>western boundary layer problem</em> (Ghil et al. 2008) for the potential vorticity <span><math><mi>ψ</mi></math></span>, which corresponds to the sharp contrast between the gyres flow in the oceanic circulation at mid-latitude and the strong western boundary currents. Numerical results for Equation (MK-1D) show that, with this two-scale scheme, high order accuracy is preserved for <span><math><mi>u</mi></math></span> and <span><math><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mi>d</mi></mrow><mrow><mi>d</mi><mi>x</mi></mrow></mfrac></mrow><mo>)</mo></mrow><mi>u</mi></mrow></math></span> both in the boundary layer and in the central zone of the domain. The test cases are taken from Chekroun et al. 2020.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"466 ","pages":"Article 116595"},"PeriodicalIF":2.1,"publicationDate":"2025-02-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143551036","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":"Auto-stabilized weak Galerkin finite element methods on polytopal meshes without convexity constraints","authors":"Chunmei Wang","doi":"10.1016/j.cam.2025.116572","DOIUrl":"10.1016/j.cam.2025.116572","url":null,"abstract":"<div><div>This paper introduces an auto-stabilized weak Galerkin (WG) finite element method with a built-in stabilizer for Poisson equations. By utilizing bubble functions as a key analytical tool, our method extends to both convex and non-convex elements in finite element partitions, marking a significant advancement over existing stabilizer-free WG methods. It overcomes the restrictive conditions of previous approaches and is applicable in any dimension <span><math><mi>d</mi></math></span>, offering substantial advantages. The proposed method maintains a simple, symmetric, and positive definite structure. These benefits are evidenced by optimal order error estimates in both discrete <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> and <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> norms, highlighting the effectiveness and accuracy of our WG method for practical applications.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"466 ","pages":"Article 116572"},"PeriodicalIF":2.1,"publicationDate":"2025-02-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143519489","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":"Non-reflecting boundary conditions for solving Maxwell’s equations in waveguides","authors":"Denise Aregba-Driollet ,&nbsp;Patrick Lacoste ,&nbsp;Corentin Prigent","doi":"10.1016/j.cam.2025.116598","DOIUrl":"10.1016/j.cam.2025.116598","url":null,"abstract":"<div><div>We present here a method for solving the electromagnetic wave–matter interaction, on the one hand in closed waveguides, single or multi-port, and on the other hand with infinite materials periodic in two directions. In this presentation, the guides are coaxial or rectangular, but extension to any regular guide section is straightforward. This paper presents a method for solving Maxwell’s equations in real or virtual waveguides, using a decomposition of the DtN operator, and also introduces an adapted basis of the solution space based on a variational formulation. Finally, 2D axisymmetric and 3D finite elements are used. The text presents a number of resolutions of guided plane-wave scattering problems in a material, confirming the validity of the method and its generality.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"466 ","pages":"Article 116598"},"PeriodicalIF":2.1,"publicationDate":"2025-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143519590","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":"Effects of higher order discretization on the predictions of steady turbulent airflow in human respiratory system","authors":"A. Lancmanová,&nbsp;T. Bodnár","doi":"10.1016/j.cam.2025.116593","DOIUrl":"10.1016/j.cam.2025.116593","url":null,"abstract":"<div><div>This study aims to evaluate the efficiency and accuracy of two distinct numerical schemes, the first order upwind scheme and the second-order central scheme with Van Leer limiter for convective terms, in predicting steady turbulent airflow in human airways for different flow rates ( <span><math><mrow><mi>Q</mi><mo>=</mo><mn>15</mn><mo>,</mo><mn>30</mn><mo>,</mo><mn>60</mn><mspace></mspace><mrow><mo>[</mo><mstyle><mi>ℓ</mi></mstyle><mspace></mspace><msup><mrow><mstyle><mi>m</mi><mi>i</mi><mi>n</mi></mstyle></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo>]</mo></mrow></mrow></math></span> ). A 3D model of the human respiratory system, extending from the oral cavity to the trachea and incorporating a simplified lung model (with branches from the 4th to 7th generations merged into larger funnels), is used to simulate steady inhalation and exhalation, excluding nasal tract respiration.</div><div>Through a series of simulations and subsequent analysis of the results, this study mutually compares the performance of the two numerical schemes in terms of their assumed accuracy, convergence, and computational efficiency with reference to available experimental data.</div><div>This study has shown significant differences in numerical predictions of respiratory flows obtained by first and second order schemes, pointing out that the second order approximation not only offers a potentially more accurate results, but also lower computational cost.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"466 ","pages":"Article 116593"},"PeriodicalIF":2.1,"publicationDate":"2025-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143527178","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 symmetric ADMM-type algorithm for robust tensor completion problems using a regularized SCAD-Schatten-p model with application in color image and video recovery","authors":"Zhechen Zhang ,&nbsp;Sanyang Liu ,&nbsp;Lixia Liu ,&nbsp;Zhiping Lin","doi":"10.1016/j.cam.2025.116604","DOIUrl":"10.1016/j.cam.2025.116604","url":null,"abstract":"<div><div>In this work, we design an efficient method for tackling the robust tensor completion problem. Tensor nuclear norm (TNN) is a conventional approach to solve the robust tensor completion problem. However, TNN may yield suboptimal solutions because the tensor rank is non-convex. With this purpose, we introduce an innovative tensor rank approximation that combines the Schatten-p norm with the Smoothly Clipped Absolute Deviation function. Additionally, we incorporate the Total Variation technique into the model to preserve the local smoothing properties of the image. To address the resulting model, we formulate a symmetric Alternating Direction Method of Multipliers algorithm (ADMM). Under some mild conditions, we validate that the solution produced by the algorithm converges to the KKT point of the model. Comprehensive experiments indicate the excellent performance of the proposed approach.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"466 ","pages":"Article 116604"},"PeriodicalIF":2.1,"publicationDate":"2025-02-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143512123","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":"Fractional multiplicative Newton-type inequalities for multiplicative s-convex positive functions with application","authors":"Abdelghani Lakhdari ,&nbsp;Djaber Chemseddine Benchettah ,&nbsp;Badreddine Meftah","doi":"10.1016/j.cam.2025.116600","DOIUrl":"10.1016/j.cam.2025.116600","url":null,"abstract":"<div><div>In this paper, we begin by demonstrating a novel fractional multiplicative integral identity applicable to multiplicative differentiable functions. Leveraging this identity, we derive fractional Newton-type inequalities for multiplicative <span><math><mi>s</mi></math></span>-convex functions. Graphical representations accompany illustrative examples to validate the accuracy of our findings. Additionally, we provide applications of our results to special means.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"465 ","pages":"Article 116600"},"PeriodicalIF":2.1,"publicationDate":"2025-02-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143508326","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":"Boundary value problems modeling moisture transport in soils","authors":"Vasyl Marynets ,&nbsp;Kateryna Marynets ,&nbsp;Oksana Kohutych","doi":"10.1016/j.cam.2025.116597","DOIUrl":"10.1016/j.cam.2025.116597","url":null,"abstract":"<div><div>To model the moisture transport in soil and to better understand physics underneath, we study a boundary value problem for a nonlinear hyperbolic PDE. Using a constructive method for approximation of solutions of the problem, we derive sufficient conditions for existence and uniqueness of its regular solutions and show that these solutions satisfy the sign-preserving inequalities. Additionally, we prove a comparison theorem and a theorem about differential inequalities, and derive an posteriori error of the method. Theoretical results are validated on an illustrative numerical example.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"465 ","pages":"Article 116597"},"PeriodicalIF":2.1,"publicationDate":"2025-02-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143508573","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}