{"title":"Anisotropic Weakly Over-Penalised Symmetric Interior Penalty Method for the Stokes Equation","authors":"Hiroki Ishizaka","doi":"10.1007/s10915-024-02598-y","DOIUrl":"https://doi.org/10.1007/s10915-024-02598-y","url":null,"abstract":"<p>In this study, we investigate an anisotropic weakly over-penalised symmetric interior penalty method for the Stokes equation on convex domains. Our approach is a simple discontinuous Galerkin method similar to the Crouzeix–Raviart finite element method. As our primary contribution, we show a new proof for the consistency term, which allows us to obtain an estimate of the anisotropic consistency error. The key idea of the proof is to apply the relation between the Raviart–Thomas finite element space and a discontinuous space. While inf-sup stable schemes of the discontinuous Galerkin method on shape-regular mesh partitions have been widely discussed, our results show that the Stokes element satisfies the inf-sup condition on anisotropic meshes. Furthermore, we provide an error estimate in an energy norm on anisotropic meshes. In numerical experiments, we compare calculation results for standard and anisotropic mesh partitions.</p>","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":"160 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141512481","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":"Optimal Balanced-Norm Error Estimate of the LDG Method for Reaction–Diffusion Problems I: The One-Dimensional Case","authors":"Yao Cheng, Xuesong Wang, Martin Stynes","doi":"10.1007/s10915-024-02602-5","DOIUrl":"https://doi.org/10.1007/s10915-024-02602-5","url":null,"abstract":"<p>A singularly perturbed reaction–diffusion problem in 1D is solved numerically by a local discontinuous Galerkin (LDG) finite element method. For this type of problem the standard energy norm is too weak to capture the contribution of the boundary layer component of the true solution, so balanced norms have been used by many authors to give more satisfactory error bounds for solutions computed using various types of finite element method. But for the LDG method, up to now no optimal-order balanced-norm error estimate has been derived. In this paper, we consider an LDG method with central numerical flux on a Shishkin mesh. Using the superconvergence property of the local <span>(L^2)</span> projector and some local coupled projections around the two transition points of the mesh, we prove an optimal-order balanced-norm error estimate for the computed solution; that is, when piecewise polynomials of degree <i>k</i> are used on a Shishkin mesh with <i>N</i> mesh intervals, in the balanced norm we establish <span>(O((N^{-1}ln N)^{k+1}))</span> convergence when <i>k</i> is even and <span>(O((N^{-1}ln N)^{k}))</span> when <i>k</i> is odd. Numerical experiments confirm the sharpness of these error bounds.</p>","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":"210 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141512482","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":"Discretization of Non-uniform Rational B-Spline (NURBS) Models for Meshless Isogeometric Analysis","authors":"Urban Duh, Varun Shankar, Gregor Kosec","doi":"10.1007/s10915-024-02597-z","DOIUrl":"https://doi.org/10.1007/s10915-024-02597-z","url":null,"abstract":"<p>We present an algorithm for fast generation of quasi-uniform and variable-spacing nodes on domains whose boundaries are represented as computer-aided design (CAD) models, more specifically non-uniform rational B-splines (NURBS). This new algorithm enables the solution of partial differential equations within the volumes enclosed by these CAD models using (collocation-based) meshless numerical discretizations. Our hierarchical algorithm first generates quasi-uniform node sets directly on the NURBS surfaces representing the domain boundary, then uses the NURBS representation in conjunction with the surface nodes to generate nodes within the volume enclosed by the NURBS surface. We provide evidence for the quality of these node sets by analyzing them in terms of local regularity and separation distances. Finally, we demonstrate that these node sets are well-suited (both in terms of accuracy and numerical stability) for meshless radial basis function generated finite differences discretizations of the Poisson, Navier-Cauchy, and heat equations. Our algorithm constitutes an important step in bridging the field of node generation for meshless discretizations with isogeometric analysis.</p>","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":"28 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141512480","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":"The Splitting Characteristic Finite Difference Domain Decomposition Scheme for Solving Time-Fractional MIM Nonlinear Advection–Diffusion Equations","authors":"Zhongguo Zhou, Sihan Zhang, Wanshan Li","doi":"10.1007/s10915-024-02603-4","DOIUrl":"https://doi.org/10.1007/s10915-024-02603-4","url":null,"abstract":"<p>In this paper, we develop a new splitting characteristic finite difference scheme for solving the time-fractional mobile-immobile nonlinear advection–diffusion equation by combining non-overlapping block-divided domain decomposition method, the operator splitting technique and the characteristic finite difference method. Over each sub-domain, the solutions and fluxes along <i>x</i>-direction in the interiors of sub-domains are computed by the implicit characteristic finite difference method while the intermediate fluxes on the interfaces of sub-domains are computed by local multi-point weighted average from the approximate solutions at characteristic tracking points which are solved by the quadratic interpolation. Secondly, the solutions and fluxes along <i>y</i> direction in the interiors of sub-domains are computed lastly by the implicit characteristic difference method while the time fractional derivative is approximated by <i>L</i>1-format and the intermediate fluxes on the interfaces of sub-domains are computed by local multi-point weighted average from the approximate solutions at characteristic tracking points are solved by the quadratic interpolation. Applying Brouwer fixed point theorem, we prove strictly the existence and uniqueness of the proposed scheme. The conditional stability and convergence with <span>(Oleft( {varDelta t}+{varDelta t}^{2-alpha }+{h}^2+{H}^frac{5}{2}right) )</span> of the proposed scheme are given as well. Numerical experiments verify the theoretical results.</p>","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":"215 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141529162","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 Robust Randomized Indicator Method for Accurate Symmetric Eigenvalue Detection","authors":"Zhongyuan Chen, Jiguang Sun, Jianlin Xia","doi":"10.1007/s10915-024-02599-x","DOIUrl":"https://doi.org/10.1007/s10915-024-02599-x","url":null,"abstract":"<p>We propose a robust randomized indicator method for the reliable detection of eigenvalue existence within an interval for symmetric matrices <i>A</i>. An indicator tells the eigenvalue existence based on some statistical norm estimators for a spectral projector. Previous work on eigenvalue indicators relies on a threshold which is empirically chosen, thus often resulting in under or over detection. In this paper, we use rigorous statistical analysis to guide the design of a robust indicator. Multiple randomized estimators for a contour integral operator in terms of <i>A</i> are analyzed. In particular, when <i>A</i> has eigenvalues inside a given interval, we show that the failure probability (for the estimators to return very small estimates) is extremely low. This enables to design a robust rejection indicator based on the control of the failure probability. We also give a prototype framework to illustrate how the indicator method may be applied numerically for eigenvalue detection and may potentially serve as a new way to design randomized symmetric eigenvalue solvers. Unlike previous indicator methods that only detect eigenvalue existence, the framework also provides a way to find eigenvectors with little extra cost by reusing computations from indicator evaluations. Extensive numerical tests show the reliability of the eigenvalue detection in multiple aspects.</p>","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":"17 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141532366","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 Multiscale Finite Element Method for an Elliptic Distributed Optimal Control Problem with Rough Coefficients and Control Constraints","authors":"Susanne C. Brenner, José C. Garay, Li-yeng Sung","doi":"10.1007/s10915-024-02590-6","DOIUrl":"https://doi.org/10.1007/s10915-024-02590-6","url":null,"abstract":"<p>We construct and analyze a multiscale finite element method for an elliptic distributed optimal control problem with pointwise control constraints, where the state equation has rough coefficients. We show that the performance of the multiscale finite element method is similar to the performance of standard finite element methods for smooth problems and present corroborating numerical results.\u0000</p>","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":"7 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141512483","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":"On the Immersed Boundary Method with Time-Filter-SAV for Solving Fluid–Structure Interaction Problem","authors":"Qixing Chen, Li Cai, Feifei Jing, Pengfei Ma, Xiaoyu Luo, Hao Gao","doi":"10.1007/s10915-024-02591-5","DOIUrl":"https://doi.org/10.1007/s10915-024-02591-5","url":null,"abstract":"<p>In this work, the immersed boundary method with time filter and scalar auxiliary variable techniques is studied to solve the fluid–structure interaction problems. For the fluid flow, we first use the backward Euler differentiation formula in temporal discretization, we then employ the time filter technique to improve its convergence order, the scalar auxiliary variable strategy is visited to decouple the fluid equations and achieve fast solutions. We adopt the immersed boundary method to build the connection between the fluid and the structure, as well as characterize the deformations of the structure. We approximate the fluid–structure interaction model by the finite element method in space. The semi-discrete and fully-discrete implicit numerical schemes are proposed followed with unconditionally stability properties. We carry out several numerical experiments to validate the convergence behaviors and efficiency of the algorithms.\u0000</p>","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":"2013 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141530302","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":"Hybridizable Discontinuous Galerkin Methods for the Two-Dimensional Monge–Ampère Equation","authors":"Ngoc Cuong Nguyen, Jaime Peraire","doi":"10.1007/s10915-024-02604-3","DOIUrl":"https://doi.org/10.1007/s10915-024-02604-3","url":null,"abstract":"<p>We introduce two hybridizable discontinuous Galerkin (HDG) methods for numerically solving the two-dimensional Monge–Ampère equation. The first HDG method is devised to solve the nonlinear elliptic Monge–Ampère equation by using Newton’s method. The second HDG method is devised to solve a sequence of the Poisson equation until convergence to a fixed-point solution of the Monge–Ampère equation is reached. Numerical examples are presented to demonstrate the convergence and accuracy of the HDG methods. Furthermore, the HDG methods are applied to <i>r</i>-adaptive mesh generation by redistributing a given scalar density function via the optimal transport theory. This <i>r</i>-adaptivity methodology leads to the Monge–Ampère equation with a nonlinear Neumann boundary condition arising from the optimal transport of the density function to conform the resulting high-order mesh to the boundary. Hence, we extend the HDG methods to treat the nonlinear Neumann boundary condition. Numerical experiments are presented to illustrate the generation of <i>r</i>-adaptive high-order meshes on planar and curved domains.\u0000</p>","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":"61 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141529164","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}
Jean-Luc Guermond, Matthias Maier, Bojan Popov, Laura Saavedra, Ignacio Tomas
{"title":"First-Order Greedy Invariant-Domain Preserving Approximation for Hyperbolic Problems: Scalar Conservation Laws, and p-System","authors":"Jean-Luc Guermond, Matthias Maier, Bojan Popov, Laura Saavedra, Ignacio Tomas","doi":"10.1007/s10915-024-02592-4","DOIUrl":"https://doi.org/10.1007/s10915-024-02592-4","url":null,"abstract":"<p>The paper focuses on first-order invariant-domain preserving approximations of hyperbolic systems. We propose a new way to estimate the artificial viscosity that has to be added to make explicit, conservative, consistent numerical methods invariant-domain preserving and entropy inequality compliant. Instead of computing an upper bound on the maximum wave speed in Riemann problems, we estimate a minimum wave speed in the said Riemann problems such that the approximation satisfies predefined invariant-domain properties and predefined entropy inequalities. This technique eliminates non-essential fast waves from the construction of the artificial viscosity, while preserving pre-assigned invariant-domain properties and entropy inequalities.</p>","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":"14 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141512485","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":"Convergent Authalic Energy Minimization for Disk Area-Preserving Parameterizations","authors":"Shu-Yung Liu, Mei-Heng Yueh","doi":"10.1007/s10915-024-02594-2","DOIUrl":"https://doi.org/10.1007/s10915-024-02594-2","url":null,"abstract":"<p>An area-preserving parameterization of a surface is a bijective mapping that maps the surface onto a specified domain and preserves the local area. This paper formulates the computation of disk area-preserving parameterization as an improved optimization problem and develops a preconditioned nonlinear conjugate gradient method with guaranteed theoretical convergence for solving the problem. Numerical experiments indicate that our new approach has significantly improved area-preserving accuracy and computational efficiency compared to other state-of-the-art algorithms. Furthermore, we present an application of surface registration to illustrate the practical utility of area-preserving mappings as parameterizations of surfaces.</p>","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":"31 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141529165","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}