{"title":"Gaussian process learning of nonlinear dynamics","authors":"Dongwei Ye, Mengwu Guo","doi":"arxiv-2312.12193","DOIUrl":"https://doi.org/arxiv-2312.12193","url":null,"abstract":"One of the pivotal tasks in scientific machine learning is to represent\u0000underlying dynamical systems from time series data. Many methods for such\u0000dynamics learning explicitly require the derivatives of state data, which are\u0000not directly available and can be approximated conventionally by finite\u0000differences. However, the discrete approximations of time derivatives may\u0000result in a poor estimation when state data are scarce and/or corrupted by\u0000noise, thus compromising the predictiveness of the learned dynamical models. To\u0000overcome this technical hurdle, we propose a new method that learns nonlinear\u0000dynamics through a Bayesian inference of characterizing model parameters. This\u0000method leverages a Gaussian process representation of states, and constructs a\u0000likelihood function using the correlation between state data and their\u0000derivatives, yet prevents explicit evaluations of time derivatives. Through a\u0000Bayesian scheme, a probabilistic estimate of the model parameters is given by\u0000the posterior distribution, and thus a quantification is facilitated for\u0000uncertainties from noisy state data and the learning process. Specifically, we\u0000will discuss the applicability of the proposed method to two typical scenarios\u0000for dynamical systems: parameter identification and estimation with an affine\u0000structure of the system, and nonlinear parametric approximation without prior\u0000knowledge.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":"73 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138818017","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Plane Wave Discontinuous Galerkin methods for scattering by periodic structures","authors":"Armando Maria Monforte","doi":"arxiv-2312.12045","DOIUrl":"https://doi.org/arxiv-2312.12045","url":null,"abstract":"This thesis explores the application of Plane Wave Discontinuous Galerkin\u0000(PWDG) methods for the numerical simulation of electromagnetic scattering by\u0000periodic structures. Periodic structures play a pivotal role in various\u0000engineering and scientific applications, including antenna design, metamaterial\u0000characterization, and photonic crystal analysis. Understanding and accurately\u0000predicting the scattering behavior of electromagnetic waves from such\u0000structures is crucial in optimizing their performance and advancing\u0000technological advancements. The thesis commences with an overview of the theoretical foundations of\u0000electromagnetic scattering by periodic structures. This theoretical\u0000dissertation serves as the basis for formulating the PWDG method within the\u0000context of wave equation. The DtN operator is presented and it is used to\u0000derive a suitable boundary condition. The numerical implementation of PWDG methods is discussed in detail,\u0000emphasizing key aspects such as basis function selection and boundary\u0000conditions. The algorithm's efficiency is assessed through numerical\u0000experiments. We then present the DtN-PWDG method, which is discussed in detail and is used\u0000to derive numerical solutions of the scattering problem. A comparison with the\u0000finite element method (FEM) is presented. In conclusion, this thesis demonstrates that PWDG methods are a powerful tool\u0000for simulating electromagnetic scattering by periodic structures.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":"9 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138818636","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An exact divergence-free spectral method for incompressible and resistive magneto-hydrodynamic equations in two and three dimensions","authors":"Lechang Qin, Huiyuan Li, Zhiguo Yang","doi":"arxiv-2312.12218","DOIUrl":"https://doi.org/arxiv-2312.12218","url":null,"abstract":"In this paper, we present exact divergence-free spectral method for solving\u0000the incompressible and resistive magneto-hydrodynamic (MHD) equations in two\u0000and three dimensions, as well as the efficient solution algorithm and\u0000unconditionally energy-stable fully-discretized numerical schemes. We introduce\u0000new ideas of constructing two families of exact divergence-free vectorial\u0000spectral basis functions on domains diffeomorphic to squares or cubes. These\u0000bases are obtained with the help of orthogonality and derivative relation of\u0000generalised Jacobi polynomials, several de Rham complexes, as well as the\u0000property of contravariant Piola transformation. They are well-suited for\u0000discretizing the velocity and magnetic fields, respectively, thereby ensuring\u0000point-wise preservation of the incompressibility condition and the magnetic\u0000Gauss's law. With the aid of these bases, we propose a family of exact\u0000divergence-free implicit-explicit $k$-step backward differentiation formula\u0000(DF-BDF-$k$) fully-discretized schemes for the MHD system. These schemes\u0000naturally decouple the pressure field from the velocity field. Consequently,\u0000the stability of the space-time fully-discretized numerical schemes based on\u0000these bases are significantly enhanced. These schemes exhibit unconditional\u0000stability for $k=1,2$, and demonstrate exceptional stability and accuracy for\u0000$k=3,4$, verified with extensive numerical results for long time simulations\u0000using large time step sizes. Furthermore, we present efficient solution\u0000algorithms for these two decoupled equations for the velocity and magnetic\u0000fields, respectively, by exploiting the sparsity and structure of the resultant\u0000linear algebraic systems. Ample numerical examples in two and three dimensions\u0000are provided to demonstrate the distinctive accuracy, efficiency and stability\u0000of our proposed method.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":"33 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138818588","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Revisiting Diffusive Representations for Enhanced Numerical Approximation of Fractional Integrals","authors":"Renu Chaudhary, Kai Diethelm","doi":"arxiv-2312.11590","DOIUrl":"https://doi.org/arxiv-2312.11590","url":null,"abstract":"This study reexamines diffusive representations for fractional integrals with\u0000the goal of pioneering new variants of such representations. These variants aim\u0000to offer highly efficient numerical algorithms for the approximate computation\u0000of fractional integrals. The approach seamlessly aligns with established\u0000techniques used in addressing problems involving integer-order operators,\u0000contributing to a unified framework for numerical solutions.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":"35 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138818334","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An abstract framework for heterogeneous coupling: stability, approximation and applications","authors":"Silvia Bertoluzza, Erik Burman","doi":"arxiv-2312.11733","DOIUrl":"https://doi.org/arxiv-2312.11733","url":null,"abstract":"Introducing a coupling framework reminiscent of FETI methods, but here on\u0000abstract form, we establish conditions for stability and minimal requirements\u0000for well-posedness on the continuous level, as well as conditions on local\u0000solvers for the approximation of subproblems. We then discuss stability of the\u0000resulting Lagrange multiplier methods and show stability under a mesh\u0000conditions between the local discretizations and the mortar space. If this\u0000condition is not satisfied we show how a stabilization, acting only on the\u0000multiplier can be used to achieve stability. The design of preconditioners of\u0000the Schur complement system is discussed in the unstabilized case. Finally we\u0000discuss some applications that enter the framework.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":"111 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138817865","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Learned Regularization for Inverse Problems: Insights from a Spectral Model","authors":"Martin Burger, Samira Kabri","doi":"arxiv-2312.09845","DOIUrl":"https://doi.org/arxiv-2312.09845","url":null,"abstract":"The aim of this paper is to provide a theoretically founded investigation of\u0000state-of-the-art learning approaches for inverse problems. We give an extended\u0000definition of regularization methods and their convergence in terms of the\u0000underlying data distributions, which paves the way for future theoretical\u0000studies. Based on a simple spectral learning model previously introduced for\u0000supervised learning, we investigate some key properties of different learning\u0000paradigms for inverse problems, which can be formulated independently of\u0000specific architectures. In particular we investigate the regularization\u0000properties, bias, and critical dependence on training data distributions.\u0000Moreover, our framework allows to highlight and compare the specific behavior\u0000of the different paradigms in the infinite-dimensional limit.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":"38 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138716944","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Alina Chertock, Michael Herty, Arsen S. Iskhakov, Safa Janajra, Alexander Kurganov, Maria Lukacova-Medvidova
{"title":"New High-Order Numerical Methods for Hyperbolic Systems of Nonlinear PDEs with Uncertainties","authors":"Alina Chertock, Michael Herty, Arsen S. Iskhakov, Safa Janajra, Alexander Kurganov, Maria Lukacova-Medvidova","doi":"arxiv-2312.08280","DOIUrl":"https://doi.org/arxiv-2312.08280","url":null,"abstract":"In this paper, we develop new high-order numerical methods for hyperbolic\u0000systems of nonlinear partial differential equations (PDEs) with uncertainties.\u0000The new approach is realized in the semi-discrete finite-volume framework and\u0000it is based on fifth-order weighted essentially non-oscillatory (WENO)\u0000interpolations in (multidimensional) random space combined with second-order\u0000piecewise linear reconstruction in physical space. Compared with spectral\u0000approximations in the random space, the presented methods are essentially\u0000non-oscillatory as they do not suffer from the Gibbs phenomenon while still\u0000achieving a high-order accuracy. The new methods are tested on a number of\u0000numerical examples for both the Euler equations of gas dynamics and the\u0000Saint-Venant system of shallow-water equations. In the latter case, the methods\u0000are also proven to be well-balanced and positivity-preserving.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":"19 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138630327","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jane Shaw MacDonald, Yves Bourgault, Frithjof Lutscher
{"title":"A hybrid finite element method for moving-habitat models in two spatial dimensions","authors":"Jane Shaw MacDonald, Yves Bourgault, Frithjof Lutscher","doi":"arxiv-2312.07842","DOIUrl":"https://doi.org/arxiv-2312.07842","url":null,"abstract":"Moving-habitat models track the density of a population whose suitable\u0000habitat shifts as a consequence of climate change. Whereas most previous\u0000studies in this area consider 1-dimensional space, we derive and study a\u0000spatially 2-dimensional moving-habitat model via reaction-diffusion equations.\u0000The population inhabits the whole space. The suitable habitat is a bounded\u0000region where population growth is positive; the unbounded complement of its\u0000closure is unsuitable with negative growth. The interface between the two\u0000habitat types moves, depicting the movement of the suitable habitat poleward.\u0000Detailed modelling of individual movement behaviour induces a nonstandard\u0000discontinuity in the density across the interface. For the corresponding\u0000semi-discretised system we prove well-posedness for a constant shifting\u0000velocity before constructing an implicit-explicit hybrid finite element method.\u0000In this method, a Lagrange multiplier weakly imposes the jump discontinuity\u0000across the interface. For a stationary interface, we derive optimal a priori\u0000error estimates over a conformal mesh with nonconformal discretisation. We\u0000demonstrate with numerical convergence tests that these results hold for the\u0000moving interface. Finally, we demonstrate the strength of our hybrid finite\u0000element method with two biologically motivated cases, one for a domain with a\u0000curved boundary and the other for non-constant shifting velocity.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":"15 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138630640","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Hussein Albazzal, Alexei Lozinski, Roberta Tittarelli
{"title":"An a posteriori error estimate for a 0D/2D coupled model","authors":"Hussein Albazzal, Alexei Lozinski, Roberta Tittarelli","doi":"arxiv-2312.07959","DOIUrl":"https://doi.org/arxiv-2312.07959","url":null,"abstract":"This work is motivated by the need of efficient numerical simulations of gas\u0000flows in the serpentine channels used in proton-exchange membrane fuel cells.\u0000In particular, we consider the Poisson problem in a 2D domain composed of\u0000several long straight rectangular sections and of several bends corners. In\u0000order to speed up the resolution, we propose a 0D model in the rectangular\u0000parts of the channel and a Finite Element resolution in the bends. To find a\u0000good compromise between precision and time consuming, the challenge is double:\u0000how to choose a suitable position of the interface between the 0D and the 2D\u0000models and how to control the discretization error in the bends. We shall\u0000present an textit{a posteriori} error estimator based on an equilibrated flux\u0000reconstruction in the subdomains where the Finite Element method is applied.\u0000The estimates give a global upper bound on the error measured in the energy\u0000norm of the difference between the exact and approximate solutions on the whole\u0000domain. They are guaranteed, meaning that they feature no undetermined\u0000constants. (global) Lower bounds for the error are also derived. An adaptive\u0000algorithm is proposed to use smartly the estimator for aforementioned double\u0000challenge. A numerical validation of the estimator and the algorithm completes\u0000the work. end{abstract}","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":"260 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138630642","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Analysis of reconstruction of functions with rough edges from discrete Radon data in $mathbb R^2$","authors":"Alexander Katsevich","doi":"arxiv-2312.08259","DOIUrl":"https://doi.org/arxiv-2312.08259","url":null,"abstract":"We study the accuracy of reconstruction of a family of functions\u0000$f_epsilon(x)$, $xinmathbb R^2$, $epsilonto0$, from their discrete Radon\u0000transform data sampled with step size $O(epsilon)$. For each $epsilon>0$\u0000sufficiently small, the function $f_epsilon$ has a jump across a rough\u0000boundary $mathcal S_epsilon$, which is modeled by an $O(epsilon)$-size\u0000perturbation of a smooth boundary $mathcal S$. The function $H_0$, which\u0000describes the perturbation, is assumed to be of bounded variation. Let\u0000$f_epsilon^{text{rec}}$ denote the reconstruction, which is computed by\u0000interpolating discrete data and substituting it into a continuous inversion\u0000formula. We prove that\u0000$(f_epsilon^{text{rec}}-K_epsilon*f_epsilon)(x_0+epsiloncheck\u0000x)=O(epsilon^{1/2}ln(1/epsilon))$, where $x_0inmathcal S$ and $K_epsilon$\u0000is an easily computable kernel.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":"21 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2023-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138630321","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}