{"about":{"site":"https://codewithpapers.app","non_affiliation":"Code with Papers and Syntology are not affiliated with, endorsed by, or sponsored by Papers with Code, Meta, or the pwc-archive mirror.","licence":"CC BY-SA 4.0","licence_url":"https://creativecommons.org/licenses/by-sa/4.0/legalcode","attribution":"https://codewithpapers.app/attribution","modified":"archive material modified by Syntology; see the attribution page"},"url":"/paper/constructing-stable-high-order-finite","title":"Constructing stable, high-order finite-difference operators on point clouds over complex geometries","arxiv_id":"2409.00809","date":"2024-09-01","proceeding":null,"authors":["Jason Hicken","Ge Yan","Sharanjeet Kaur"],"abstract":"High-order difference operators with the summation-by-parts (SBP) property can be used to build stable discretizations of hyperbolic conservation laws; however, most high-order SBP operators require a conforming, high-order mesh for the domain of interest. To circumvent this requirement, we present an algorithm for building high-order, diagonal-norm, first-derivative SBP operators on point clouds over level-set geometries. The algorithm is not mesh-free, since it uses a Cartesian cut-cell mesh to define the sparsity pattern of the operators and to provide intermediate quadrature rules; however, the mesh is generated automatically and can be discarded once the SBP operators have been constructed. Using this temporary mesh, we construct local, cell-based SBP difference operators that are assembled into global SBP operators. We identify conditions for the existence of a positive-definite diagonal mass matrix, and we compute the diagonal norm by solving a sparse system of linear inequalities using an interior-point algorithm. We also describe an artificial dissipation operator that complements the first-derivative operators when solving hyperbolic problems, although the dissipation is not required for stability. The numerical results confirm the conditions under which a diagonal norm exists and study the distribution of the norm's entries. In addition, the results verify the accuracy and stability of the point-cloud SBP operators using the linear advection equation.","url_abs":"https://arxiv.org/abs/2409.00809v2","url_pdf":"https://arxiv.org/pdf/2409.00809v2.pdf","source":{"archive":"pwc-archive (Hugging Face), CC BY-SA 4.0","snapshot":"2025-07-28","licence_url":"https://creativecommons.org/licenses/by-sa/4.0/legalcode","row_kind":"links_only","authors_date_abstract":"arXiv metadata, CC0 1.0 (https://info.arxiv.org/help/license), from the Kaggle arXiv metadata snapshot of 2026-09-12"},"code_links":[{"paper_slug":"constructing-stable-high-order-finite","repo_url":"https://github.com/jehicken/cloudsbp.jl","is_official":1,"mentioned_in_paper":1,"mentioned_in_github":1,"framework":"none","reach":null}],"tasks":[],"methods":[],"datasets_introduced":[],"methods_introduced":[],"results":[],"syntology":{"atlas_url":null,"mcp":null,"developers":"https://syntology.ai/developers"},"arxiv_metadata":null,"syntology_extracted_results":null}