Communications on Applied Mathematics and Computation最新文献

{"title":"A Well-Balanced Active Flux Method for the Shallow Water Equations with Wetting and Drying","authors":"Wasilij Barsukow, J. P. Berberich","doi":"10.48550/arXiv.2212.02426","DOIUrl":"https://doi.org/10.48550/arXiv.2212.02426","url":null,"abstract":"Active Flux is a third order accurate numerical method which evolves cell averages and point values at cell interfaces independently. It naturally uses a continuous reconstruction, but is stable when applied to hyperbolic problems. In this work, the Active Flux method is extended for the first time to a nonlinear hyperbolic system of balance laws, namely, to the shallow water equations with bottom topography. We demonstrate how to achieve an Active Flux method that is well-balanced, positivity preserving, and allows for dry states in one spatial dimension. Because of the continuous reconstruction all these properties are achieved using new approaches. To maintain third order accuracy, we also propose a novel high-order approximate evolution operator for the update of the point values. A variety of test problems demonstrates the good performance of the method even in presence of shocks.","PeriodicalId":29916,"journal":{"name":"Communications on Applied Mathematics and Computation","volume":"67 1","pages":"1-46"},"PeriodicalIF":1.6,"publicationDate":"2022-12-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78649495","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Exponentially Convergent Multiscale Finite Element Method","authors":"Yifan Chen, T. Hou, Yixuan Wang","doi":"10.48550/arXiv.2212.00823","DOIUrl":"https://doi.org/10.48550/arXiv.2212.00823","url":null,"abstract":"We provide a concise review of the exponentially convergent multiscale finite element method (ExpMsFEM) for efficient model reduction of PDEs in heterogeneous media without scale separation and in high-frequency wave propagation. The ExpMsFEM is built on the non-overlapped domain decomposition in the classical MsFEM while enriching the approximation space systematically to achieve a nearly exponential convergence rate regarding the number of basis functions. Unlike most generalizations of the MsFEM in the literature, the ExpMsFEM does not rely on any partition of unity functions. In general, it is necessary to use function representations dependent on the right-hand side to break the algebraic Kolmogorov n -width barrier to achieve exponential convergence. Indeed, there are online and offline parts in the function representation provided by the ExpMsFEM. The online part depends on the right-hand side locally and can be computed in parallel efficiently. The offline part contains basis functions that are used in the Galerkin method to assemble the stiffness matrix; they are all independent of the right-hand side, so the stiffness matrix can be used repeatedly in multi-query scenarios.","PeriodicalId":29916,"journal":{"name":"Communications on Applied Mathematics and Computation","volume":"1 1","pages":"1-17"},"PeriodicalIF":1.6,"publicationDate":"2022-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90316726","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Preface to the Focused Issue on High-Order Numerical Methods for Evolutionary PDEs","authors":"A. Hidalgo, M. Dumbser, E. Toro","doi":"10.1007/s42967-022-00229-7","DOIUrl":"https://doi.org/10.1007/s42967-022-00229-7","url":null,"abstract":"","PeriodicalId":29916,"journal":{"name":"Communications on Applied Mathematics and Computation","volume":"4 1","pages":"529-531"},"PeriodicalIF":1.6,"publicationDate":"2022-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78788968","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Preface to the Focused Issue in Honor of Professor Tong Zhang on the Occasion of His 90th Birthday","authors":"Jiequan Li, Wancheng Sheng, Chi-Wang Shu, P. Zhang, Yuxi Zheng","doi":"10.1007/s42967-022-00230-0","DOIUrl":"https://doi.org/10.1007/s42967-022-00230-0","url":null,"abstract":"","PeriodicalId":29916,"journal":{"name":"Communications on Applied Mathematics and Computation","volume":"1 1","pages":"965 - 966"},"PeriodicalIF":1.6,"publicationDate":"2022-11-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86675649","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Techniques, Tricks, and Algorithms for Efficient GPU-Based Processing of Higher Order Hyperbolic PDEs","authors":"Sethupathy S., D. Balsara, Deepak Bhoriya, Harish Kumar","doi":"10.1007/s42967-022-00235-9","DOIUrl":"https://doi.org/10.1007/s42967-022-00235-9","url":null,"abstract":"","PeriodicalId":29916,"journal":{"name":"Communications on Applied Mathematics and Computation","volume":"14 1","pages":"1-49"},"PeriodicalIF":1.6,"publicationDate":"2022-11-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78782111","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Characterizations and Properties of Dual Matrix Star Orders","authors":"Hongxing Wang, P. Huang","doi":"10.1007/s42967-023-00255-z","DOIUrl":"https://doi.org/10.1007/s42967-023-00255-z","url":null,"abstract":"","PeriodicalId":29916,"journal":{"name":"Communications on Applied Mathematics and Computation","volume":"78 1","pages":"1-24"},"PeriodicalIF":1.6,"publicationDate":"2022-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88565049","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Bound-Preserving Discontinuous Galerkin Methods with Modified Patankar Time Integrations for Chemical Reacting Flows","authors":"Fangyao Zhu, Juntao Huang, Yang Yang","doi":"10.48550/arXiv.2211.05600","DOIUrl":"https://doi.org/10.48550/arXiv.2211.05600","url":null,"abstract":"In this paper, we develop bound-preserving discontinuous Galerkin (DG) methods for chemical reactive flows. There are several difficulties in constructing suitable numerical schemes. First of all, the density and internal energy are positive, and the mass fraction of each species is between 0 and 1. Second, due to the rapid reaction rate, the system may contain stiff sources, and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes. To obtain physically relevant numerical approximations, we apply the bound-preserving technique to the DG methods. Though traditional positivity-preserving techniques can successfully yield positive density, internal energy, and mass fractions, they may not enforce the upper bound 1 of the mass fractions. To solve this problem, we need to (i) make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density; (ii) choose conservative time integrations, such that the summation of the mass fractions is preserved. With the above two conditions, the positive mass fractions have summation 1, and then, they are all between 0 and 1. For time discretization, we apply the modified Runge-Kutta/multi-step Patankar methods, which are explicit for the flux while implicit for the source. Such methods can handle stiff sources with relatively large time steps, preserve the positivity of the target variables, and keep the summation of the mass fractions to be 1. Finally, it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations. The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces, while Patankar methods can keep the positivity of the pre-selected point values of the target variables. To match the degree of freedom, we use $$Q^k$$ Q k polynomials on rectangular meshes for problems in two space dimensions. To evolve in time, we first read the polynomials at the Gaussian points. Then, suitable slope limiters can be applied to enforce the positivity of the solutions at those points, which can be preserved by the Patankar methods, leading to positive updated numerical cell averages. In addition, we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux. Numerical examples are given to demonstrate the good performance of the proposed schemes.","PeriodicalId":29916,"journal":{"name":"Communications on Applied Mathematics and Computation","volume":"47 1","pages":"1-28"},"PeriodicalIF":1.6,"publicationDate":"2022-11-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"84145003","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical Investigations on the Resonance Errors of Multiscale Discontinuous Galerkin Methods for One-Dimensional Stationary Schrödinger Equation","authors":"Bo Dong, Wei Wang","doi":"10.48550/arXiv.2211.04900","DOIUrl":"https://doi.org/10.48550/arXiv.2211.04900","url":null,"abstract":"In this paper, numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin (DG) methods (Dong et al. in J Sci Comput 66: 321–345, 2016; Dong and Wang in J Comput Appl Math 380: 1–11, 2020) for a one-dimensional stationary Schrödinger equation. Previous work showed that penalty parameters were required to be positive in error analysis, but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes. In this work, by performing extensive numerical experiments, we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods, and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.","PeriodicalId":29916,"journal":{"name":"Communications on Applied Mathematics and Computation","volume":"216 1","pages":"1-14"},"PeriodicalIF":1.6,"publicationDate":"2022-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91388700","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Norms of Dual Complex Vectors and Dual Complex Matrices","authors":"Xin He Miao, Zhenghai Huang","doi":"10.1007/s42967-022-00215-z","DOIUrl":"https://doi.org/10.1007/s42967-022-00215-z","url":null,"abstract":"","PeriodicalId":29916,"journal":{"name":"Communications on Applied Mathematics and Computation","volume":"28 1","pages":"1484 - 1508"},"PeriodicalIF":1.6,"publicationDate":"2022-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"74343877","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Existence of Two Limit Cycles in Zeeman’s Class 30 for 3D Lotka-Volterra Competitive System","authors":"Yaoqi Li","doi":"10.1007/s42967-022-00220-2","DOIUrl":"https://doi.org/10.1007/s42967-022-00220-2","url":null,"abstract":"","PeriodicalId":29916,"journal":{"name":"Communications on Applied Mathematics and Computation","volume":"9 1","pages":"1584-1590"},"PeriodicalIF":1.6,"publicationDate":"2022-11-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"82743894","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}