{"title":"Regularity results and numerical solution by the discontinuous Galerkin method to semilinear parabolic initial boundary value problems with nonlinear Newton boundary conditions in a polygonal space-time cylinder","authors":"M. Balázsová, M. Feistauer, A. Sändig","doi":"10.1515/jnma-2021-0113","DOIUrl":"https://doi.org/10.1515/jnma-2021-0113","url":null,"abstract":"Abstract In this note we consider a parabolic evolution equation in a polygonal space-time cylinder. We show, that the elliptic part is given by a m-accretive mapping from Lq(Ω) → Lq(Ω). Therefore we can apply the theory of nonlinear semigroups in Banach spaces in order to get regularity results in time and space. The second part of the paper deals with the numerical solution of the problem. It is dedicated to the application of the space-time discontinuous Galerkin method (STDGM). It means that both in space as well as in time discontinuous piecewise polynomial approximations of the solution are used. We concentrate to the theoretical analysis of the error estimation.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86078351","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":"Fourier analysis of a time-simultaneous two-grid algorithm using a damped Jacobi waveform relaxation smoother for the one-dimensional heat equation","authors":"C. Lohmann, J. Dünnebacke, S. Turek","doi":"10.1515/jnma-2021-0045","DOIUrl":"https://doi.org/10.1515/jnma-2021-0045","url":null,"abstract":"Abstract In this work, the convergence behavior of a time-simultaneous two-grid algorithm for the one-dimensional heat equation is studied using Fourier arguments in space. The underlying linear system of equations is obtained by a finite element or finite difference approximation in space while the semi-discrete problem is discretized in time using the ϑ-scheme. The simultaneous treatment of all time instances leads to a global system of linear equations which provides the potential for a higher degree of parallelization of multigrid solvers due to the increased number of degrees of freedom per spatial unknown. It is shown that the all-at-once system based on an equidistant discretization in space and time stays well conditioned even if the number of blocked time-steps grows arbitrarily. Furthermore, mesh-independent convergence rates of the considered two-grid algorithm are proved by adopting classical Fourier arguments in space without assuming periodic boundary conditions. The rate of convergence with respect to the Euclidean norm does not deteriorate arbitrarily if the number of blocked time steps increases and, hence, underlines the potential of the solution algorithm under investigation. Numerical studies demonstrate why minimizing the spectral norm of the iteration matrix may be practically more relevant than improving the asymptotic rate of convergence.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-06-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76951897","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":"Diagonally implicit Runge-Kutta schemes: Discrete energy-balance laws and compactness properties","authors":"A. Salgado, I. Tomas","doi":"10.48550/arXiv.2205.13032","DOIUrl":"https://doi.org/10.48550/arXiv.2205.13032","url":null,"abstract":"Abstract We study diagonally implicit Runge-Kutta (DIRK) schemes when applied to abstract evolution problems that fit into the Gelfand-triple framework. We introduce novel stability notions that are well-suited to this setting and provide simple, necessary and sufficient, conditions to verify that a DIRK scheme is stable in our sense and in Bochner-type norms. We use several popular DIRK schemes in order to illustrate cases that satisfy the required structural stability properties and cases that do not. In addition, under some mild structural conditions on the problem we can guarantee compactness of families of discrete solutions with respect to time discretization.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-05-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"84938293","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}
Elena Bachini, Philip Brandner, Thomas Jankuhn, M. Nestler, S. Praetorius, A. Reusken, A. Voigt
{"title":"Diffusion of tangential tensor fields: numerical issues and influence of geometric properties","authors":"Elena Bachini, Philip Brandner, Thomas Jankuhn, M. Nestler, S. Praetorius, A. Reusken, A. Voigt","doi":"10.48550/arXiv.2205.12581","DOIUrl":"https://doi.org/10.48550/arXiv.2205.12581","url":null,"abstract":"Abstract We study the diffusion of tangential tensor-valued data on curved surfaces. For this purpose, several finite-element-based numerical methods are collected and used to solve a tangential surface n-tensor heat flow problem. These methods differ with respect to the surface representation used, the geometric information required, and the treatment of the tangentiality condition. We emphasize the importance of geometric properties and their increasing influence as the tensorial degree changes from n = 0 to n ≥ 1. A specific example is presented that illustrates how curvature drastically affects the behavior of the solution.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-05-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86812365","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 time-explicit weak Galerkin scheme for parabolic equations on polytopal partitions","authors":"Junping Wang, X. Ye, Shangyou Zhang","doi":"10.1515/jnma-2021-0128","DOIUrl":"https://doi.org/10.1515/jnma-2021-0128","url":null,"abstract":"Abstract In this paper a time-explicit weak Galerkin finite element method is introduced and analyzed for parabolic equations. The main idea relies on the inclusion of a stabilization term in the temporal direction in addition to the usual static stabilization in the weak Galerkin framework. Both semi-discrete and fully-discrete schemes in time are presented, as well as their stability and error analysis. Numerical results are reported for this new explicit weak Galerkin finite element method.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-05-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"83364144","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}
G. Gatica, Cristian Inzunza, R. Ruiz-Baier, Felipe Sandoval
{"title":"A posteriori error analysis of Banach spaces-based fully-mixed finite element methods for Boussinesq-type models","authors":"G. Gatica, Cristian Inzunza, R. Ruiz-Baier, Felipe Sandoval","doi":"10.1515/jnma-2021-0101","DOIUrl":"https://doi.org/10.1515/jnma-2021-0101","url":null,"abstract":"Abstract In this paper we consider Banach spaces-based fully-mixed variational formulations recently proposed for the Boussinesq and the Oberbeck–Boussinesq models, and develop reliable and efficient residual-based a posteriori error estimators for the 2D and 3D versions of the associated mixed finite element schemes. For the reliability analysis, we employ the global inf-sup condition for each sub-model, namely Navier–Stokes and heat equations in the case of Boussinesq, along with suitable Helmholtz decomposition in nonstandard Banach spaces, the approximation properties of the Raviart–Thomas and Clément interpolants, further regularity on the continuous solutions, and small data assumptions. In turn, the efficiency estimates follow from inverse inequalities and the localization technique through bubble functions in adequately defined local Lp spaces. Finally, several numerical results including natural convection in 3D differentially heated enclosures, are reported with the aim of confirming the theoretical properties of the estimators and illustrating the performance of the associated adaptive algorithm.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-04-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87041067","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":"How to prove optimal convergence rates for adaptive least-squares finite element methods","authors":"Philipp Bringmann","doi":"10.1515/jnma-2021-0116","DOIUrl":"https://doi.org/10.1515/jnma-2021-0116","url":null,"abstract":"Abstract The convergence analysis with rates for adaptive least-squares finite element methods (ALSFEMs) combines arguments from the a posteriori analysis of conforming and mixed finite element schemes. This paper provides an overview of the key arguments for the verification of the axioms of adaptivity for an ALSFEM for the solution of a linear model problem. The formulation at hand allows for the simultaneous analysis of first-order systems of the Poisson model problem, the Stokes equations, and the linear elasticity equations. Following [Carstensen and Park, SIAM J. Numer. Anal. 53(1), 2015], the adaptive algorithm is driven by an alternative residual-based error estimator with exact solve and includes a separate marking strategy for quasi-optimal data resolution of the right-hand side. This presentation covers conforming discretisations for an arbitrary polynomial degree and mixed homogeneous boundary conditions.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-04-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87369037","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 structure preserving front tracking finite element method for the Mullins–Sekerka problem","authors":"R. Nürnberg","doi":"10.1515/jnma-2021-0131","DOIUrl":"https://doi.org/10.1515/jnma-2021-0131","url":null,"abstract":"Abstract We introduce and analyse a fully discrete approximation for a mathematical model for the solidification and liquidation of materials of negligible specific heat. The model is a two-sided Mullins–Sekerka problem. The discretization uses finite elements in space and an independent parameterization of the moving free boundary. We prove unconditional stability and exact volume conservation for the introduced scheme. Several numerical simulations, including for nearly crystalline surface energies, demonstrate the practicality and accuracy of the presented numerical method.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-03-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"72998054","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":"Adaptive POD-DEIM correction for Turing pattern approximation in reaction–diffusion PDE systems","authors":"A. Alla, A. Monti, I. Sgura","doi":"10.48550/arXiv.2203.05998","DOIUrl":"https://doi.org/10.48550/arXiv.2203.05998","url":null,"abstract":"Abstract We investigate a suitable application of Model Order Reduction (MOR) techniques for the numerical approximation of Turing patterns, that are stationary solutions of reaction–diffusion PDE (RD-PDE) systems. We show that solutions of surrogate models built by classical Proper Orthogonal Decomposition (POD) exhibit an unstable error behaviour over the dimension of the reduced space. To overcome this drawback, first of all, we propose a POD-DEIM technique with a correction term that includes missing information in the reduced models. To improve the computational efficiency, we propose an adaptive version of this algorithm in time that accounts for the peculiar dynamics of the RD-PDE in presence of Turing instability. We show the effectiveness of the proposed methods in terms of accuracy and computational cost for a selection of RD systems, i.e., FitzHugh–Nagumo, Schnakenberg and the morphochemical DIB models, with increasing degree of nonlinearity and more structured patterns.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81964235","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":"Frontmatter","authors":"","doi":"10.1515/jnma-2022-frontmatter1","DOIUrl":"https://doi.org/10.1515/jnma-2022-frontmatter1","url":null,"abstract":"Article Frontmatter was published on March 1, 2022 in the journal Journal of Numerical Mathematics (volume 30, issue 1).","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138530246","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}