Numerical Analysis
New submissions
[ showing up to 2000 entries per page: fewer | more ]
New submissions for Wed, 25 Dec 19
- [1] arXiv:1912.11129 [pdf, other]
-
Title: Uniqueness of an inverse source problem in experimental aeroacousticsSubjects: Numerical Analysis (math.NA)
This paper is concerned with the mathematical analysis of experimental methods for the estimation of the power of an uncorrelated, extended aeroacoustic source from measurements of correlations of pressure fluctuations. We formulate a continuous, infinite dimensional model describing these experimental techniques based on the convected Helmholtz equation in $\mathbb{R}^3$ or $\mathbb{R}^2$. As a main result we prove that an unknown, compactly supported source power function is uniquely determined by idealized, noise-free correlation measurements. Our framework further allows for a precise characterization of state-of-the-art source reconstruction methods and their interrelations.
- [2] arXiv:1912.11130 [pdf, other]
-
Title: Using trullekrul in pde2path -- anisotropic mesh adaptation for some Allen-Cahn models in 2D and 3DAuthors: Hannes UeckerComments: 10 pages, 6 figuresSubjects: Numerical Analysis (math.NA); Pattern Formation and Solitons (nlin.PS)
We describe by means of some examples how some functionality of the mesh adaptation package trullekrul can be used in pde2path.
- [3] arXiv:1912.11138 [pdf, other]
-
Title: Nonlinear Galerkin Model Reduction for Systems with Multiple Transport VelocitiesSubjects: Numerical Analysis (math.NA)
We propose a new model reduction framework for problems that exhibit transport phenomena. As in the moving finite element method (MFEM), our method employs time-dependent transformation operators and, especially, generalizes MFEM to arbitrary basis functions. The new framework is suitable to obtain a low-dimensional approximation with small errors even in situations where classical model order reduction techniques require much higher dimensions for a similar approximation quality. Analogously to the MFEM framework, the reduced model is designed to minimize the residual, which is also the basis for an a-posteriori error bound. Moreover, since the dependence of the transformation operators on the reduced state is nonlinear, the resulting reduced order model is obtained by projecting the original evolution equation onto a nonlinear manifold. Furthermore, for a special case, we show a connection between our approach and the method of freezing, which is also known as symmetry reduction. Besides the construction of the reduced order model, we also analyze the problem of finding optimal basis functions based on given data of the full order solution. Especially, we show that the corresponding minimization problem has a solution and reduces to the proper orthogonal decomposition of transformed data in a special case. Finally, we demonstrate the effectiveness of our method with several analytical and numerical examples.
- [4] arXiv:1912.11182 [pdf, ps, other]
-
Title: Analysis of adaptive BDF2 scheme for diffusion equationsComments: 20 pagesSubjects: Numerical Analysis (math.NA)
The variable two-step backward differentiation formula (BDF2) is revisited via a new theoretical framework using the positive semi-definiteness of BDF2 convolution kernels and a class of orthogonal convolution kernels. We prove that, if the adjacent time-step ratios $r_k:=\tau_k/\tau_{k-1}\le(3+\sqrt{17})/2\approx3.561$, the adaptive BDF2 time-stepping scheme for linear reaction-diffusion equations is unconditionally stable and (maybe, first-order) convergent in the $L^2$ norm. The second-order temporal convergence can be recovered if almost all of time-step ratios $r_k\le 1+\sqrt{2}$ or some high-order starting scheme is used. Specially, for linear dissipative diffusion problems, the stable BDF2 method preserves both the energy dissipation law (in the $H^1$ seminorm) and the $L^2$ norm monotonicity at the discrete levels. An example is included to support our analysis.
- [5] arXiv:1912.11207 [pdf, other]
-
Title: Implicit shock tracking using an optimization-based, $r$-adaptive, high-order discontinuous Galerkin methodComments: 35 pages, 20 figuresSubjects: Numerical Analysis (math.NA); Optimization and Control (math.OC)
A novel framework for resolving discontinuous solutions of conservation laws, e.g., contact lines, shock waves, and interfaces, using implicit tracking and a high-order discontinuous Galerkin (DG) discretization was introduced in [38]. Central to the framework is an optimization problem whose solution is a discontinuity-aligned mesh and the corresponding high-order approximation to the flow that does not require explicit meshing of the unknown discontinuity surface. The method was shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover optimal convergence rates even for problems with discontinuous solutions. This work extends the implicit tracking framework such that robustness is improved and convergence accelerated. In particular, we introduce an improved formulation of the central optimization problem and an associated sequential quadratic programming (SQP) solver. The new error-based objective function penalizes violation of the DG residual in an enriched test space and is shown to have excellent tracking properties. The SQP solver simultaneously converges the nodal coordinates of the mesh and DG solution to their optimal values and is equipped with a number of features to ensure robust, fast convergence: Levenberg-Marquardt approximation of the Hessian with weighted elliptic regularization, backtracking line search, and rigorous convergence criteria. We use the proposed method to solve a range of inviscid conservation laws of varying difficulty. We show the method is able to deliver accurate solutions on coarse, high-order meshes and the SQP solver is robust and usually able to drive the first-order optimality system to tight tolerances.
- [6] arXiv:1912.11293 [pdf, ps, other]
-
Title: Singular Value Decomposition in Sobolev Spaces: Part IIComments: 17 pages, 2 figures. Portions of this work previously appeared as arXiv:1809.11001 which was split during the refereeing for publicationSubjects: Numerical Analysis (math.NA); Functional Analysis (math.FA)
Under certain conditions, an element of a tensor product space can be identified with a compact operator and the singular value decomposition (SVD) applies to the latter. These conditions are not fulfilled in Sobolev spaces. In the previous part of this work (part I), we introduced some preliminary notions in the theory of tensor product spaces. We analyzed low-rank approximations in H1 and the error of the SVD performed in the ambient L2 space. In this work (part II), we continue by considering variants of the SVD in norms stronger than the L2-norm. Overall and, perhaps surprisingly, this leads to a more difficult control of the H1-error. We briefly consider an isometric embedding of H1 that allows direct application of the SVD to H1-functions. Finally, we provide a few numerical examples that support our theoretical findings.
- [7] arXiv:1912.11347 [pdf, ps, other]
-
Title: Robust Group Synchronization via Cycle-Edge Message PassingSubjects: Numerical Analysis (math.NA); Information Theory (cs.IT); Optimization and Control (math.OC); Probability (math.PR)
We propose a general framework for group synchronization with adversarial corruption and sufficiently small noise. Specifically, we apply a novel message passing procedure that uses cycle consistency information in order to estimate the corruption levels of group ratios and consequently infer the corrupted group ratios and solve the synchronization problem. We first explain why the group cycle consistency information is essential for effectively solving group synchronization problems. We then establish exact recovery and linear convergence guarantees for the proposed message passing procedure under a deterministic setting with adversarial corruption. These guarantees hold as long as the ratio of corrupted cycles per edge is bounded by a reasonable constant. We also establish the stability of the proposed procedure to sub-Gaussian noise. We further show that under a uniform corruption model, the recovery results are sharp in terms of an information-theoretic bound.
- [8] arXiv:1912.11375 [pdf, other]
-
Title: Diffuse optical tomography by simulated annealing via a spin HamiltonianSubjects: Numerical Analysis (math.NA); Optics (physics.optics)
The inverse problem of diffuse optical tomography is solved by the Markov-chain Monte Carlo. The Metropolis algorithm or single-component Metropolis-Hastings algorithm is used. The value of the unknown parameter is disretized and a spin Hamiltonian is introduced in the cost function. Then an initial random spin configuration is brought to a converged configuration by simulated annealing.
- [9] arXiv:1912.11382 [pdf, other]
-
Title: Sampling-free parametric model reduction for structured systemsSubjects: Numerical Analysis (math.NA); Systems and Control (eess.SY)
We consider the reduction of parametric families of linear dynamical systems having an affine parameter dependence that differ from one another by a low-rank variation in the state matrix. Usual approaches for parametric model reduction typically involve exploring the parameter space to isolate representative models on which to focus model reduction methodology, which are then combined in various ways in order to interpolate the response from these representative models. The initial exploration of the parameter space can be a forbiddingly expensive task. A different approach is proposed here that does not require any parameter sampling or exploration of the parameter space. Instead, we represent the system response in terms of four subsystems that are nonparametric. One may apply any one of a number of standard (nonparametric) model reduction strategies to reduce the subsystems independently, and then conjoin these reduced models with the underlying parameterized representation to obtain an overall parameterized response. Our approach has elements in common with the parameter mapping approach of Baur et al. [8], but offers greater flexibility and potentially greater control over accuracy. In particular, a data-driven variation of our approach is described that exercises this flexibility through the use of limited frequency-sampling of the underlying nonparametric models. The parametric structure of our system representation allows for a priori guarantees of system stability the resulting parametric reduced models, uniformly across all parameter values. Incorporation of system theoretic error bounds allow us to determine appropriate approximation orders for the nonparametric systems sufficient to yield uniformly high accuracy with respect to parameter variation.
- [10] arXiv:1912.11391 [pdf, ps, other]
-
Title: A Framework for Data-Driven Computational Dynamics Based on Nonlinear OptimizationComments: arXiv admin note: substantial text overlap with arXiv:1910.12736Subjects: Numerical Analysis (math.NA)
In this article, we present an extension of the formulation recently developed by the authors (A Framework for Data-Driven Computational Mechanics Based on Nonlinear Optimization, arXiv:1910.12736 [math.NA]) to the structural dynamics setting. Inspired by a structure-preserving family of variational integrators, our new formulation relies on a discrete balance equation that establishes the dynamic equilibrium. From this point of departure, we first derive an "exact" discrete-continuous nonlinear optimization problem that works directly with data sets. We then develop this formulation further into an "approximate" nonlinear optimization problem that relies on a general constitutive model. This underlying model can be identified from a data set in an offline phase. To showcase the advantages of our framework, we specialize our methodology to the case of a geometrically exact beam formulation that makes use of all elements of our approach. We investigate three numerical examples of increasing difficulty that demonstrate the excellent computational behavior of the proposed framework and motivate future research in this direction.
Cross-lists for Wed, 25 Dec 19
- [11] arXiv:1912.11110 (cross-list from math.ST) [pdf, ps, other]
-
Title: Kernel Embedding Linear ResponseComments: 28 pages, 14 figuresSubjects: Statistics Theory (math.ST); Numerical Analysis (math.NA)
In the paper, we study the problem of estimating linear response statistics under external perturbations using time series of unperturbed dynamics. A standard approach to this estimation problem is to employ the Fluctuation-Dissipation Theory, which, in turn, requires the knowledge of the functional form of the underlying unperturbed density that is not available in general. To overcome this issue, we consider a nonparametric density estimator formulated by the kernel embedding of distribution. To avoid the computational expense arises using radial type kernels, we consider the "Mercer-type" kernels constructed based on the classical orthogonal bases defined on non-compact domains, such as the Hermite and Laguerre polynomials. We establish the uniform convergence of the estimator, which justifies the use of the estimator for interpolation. Given a target function with a specific decaying property quantified by the available data, our framework allows one to choose the appropriate hypothesis space (an RKHS) that is "rich" enough for consistent estimation. For the linear response estimation, our study provides practical conditions for the well-posedness of both the estimator and the underlying response statistics. We offer a theoretical guarantee for the convergence of the estimator to the underlying actual linear response statistics. Finally, we provide a statistical error bound for the density estimation that accounts for the Monte-Carlo averaging over non-i.i.d time series and biases due to finite basis truncation. This error bound provides a mean to understand the feasibility as well as limitation of the kernel embedding with Mercer-type kernel. Numerically, we verify the effectiveness of the kernel embedding linear response estimator on two stochastic dynamics with known, yet, non-trivial equilibrium densities.
- [12] arXiv:1912.11292 (cross-list from cs.CE) [pdf, other]
-
Title: Computational framework for monolithic coupling for thin fluid flow in contact interfacesAuthors: Andrei G. Shvarts (1 and 2), Julien Vignollet (3), Vladislav A. Yastrebov (1) ((1) MINES ParisTech, PSL Research University, Centre des Matériaux, CNRS UMR 7633, France, (2) Glasgow Computational Engineering Centre, School of Engineering, University of Glasgow, United Kingdom, (3) Safran Tech, Safran Group, France)Comments: 35 pages, 14 figuresSubjects: Computational Engineering, Finance, and Science (cs.CE); Numerical Analysis (math.NA); Computational Physics (physics.comp-ph)
We develop a computational framework for simulating thin fluid flow in narrow interfaces between contacting solids, which is relevant for a range engineering, biological and geophysical applications. The treatment of this problem requires coupling between fluid and solid mechanics equations, further complicated by contact constraints and potentially complex geometrical features of contacting surfaces. We develop a monolithic finite-element framework for handling contact, thin incompressible viscous flow and fluid-induced tractions on the surface of the solid, suitable for both one- and two-way coupling approaches. Additionally, we consider fluid entrapment in "pools" delimited by contact patches and its pressurisation following a non-linear compressible constitutive law. Image analysis algorithms are adopted to identify the local status of each interface element (i.e. distinguish between contact, fluid flow and trapped fluid zones) within the Newton convergence loop. First, an application of the proposed framework for a problem with a model geometry is given, and the robustness is demonstrated by the DOF-wise and status-wise convergence. The full capability of the developed two-way coupling framework is demonstrated on a problem of a fluid flow in a contact interface between a solid with representative rough surface and a rigid flat. The evolution of the contact pressure, fluid flow pattern and the morphology of trapped fluid zones under increasing external load until the complete sealing of the interface is displayed. Finally, effective properties of flat-on-flat rough contact interfaces such as transmissivity and real contact area growth are calculated using the developed framework, showing qualitatively new results compared to the one-way coupling approximation.
- [13] arXiv:1912.11341 (cross-list from econ.GN) [pdf, other]
-
Title: Quantifying the Effects of the 2008 Recession using the Zillow DatasetSubjects: General Economics (econ.GN); Machine Learning (cs.LG); Numerical Analysis (math.NA)
This report explores the use of Zillow's housing metrics dataset to investigate the effects of the 2008 US subprime mortgage crisis on various US locales. We begin by exploring the causes of the recession and the metrics available to us in the dataset. We settle on using the Zillow Home Value Index (ZHVI) because it is seasonally adjusted and able to account for a variety of inventory factors. Then, we explore three methodologies for quantifying recession impact: (a) Principal Components Analysis, (b) Area Under Baseline, and (c) ARIMA modeling and Confidence Intervals. While PCA does not yield useable results, we ended up with six cities from both AUB and ARIMA analysis, the top 3 "losers" and "gainers" of the 2008 recession, as determined by each analysis. This gave us 12 cities in total. Finally, we tested the robustness of our analysis against three "common knowledge" metrics for the recession: geographic clustering, population trends, and unemployment rate. While we did find some overlap between the results of our analysis and geographic clustering, there was no positive regression outcome from comparing our methodologies to population trends and the unemployment rate.
Replacements for Wed, 25 Dec 19
- [14] arXiv:1807.02718 (replaced) [pdf, other]
-
Title: High-order, Dispersionless "Fast-Hybrid" Wave Equation Solver. Part I: $\mathcal{O}(1)$ Sampling Cost via Incident-Field Windowing and RecenteringComments: 32 pages, 8 figures, revised and extended manuscript (and now including direct comparisons to existing CQ and TDIE solver implementations) (Part I of II)Subjects: Numerical Analysis (math.NA); Computational Physics (physics.comp-ph)
- [15] arXiv:1808.10720 (replaced) [pdf, ps, other]
-
Title: An explicit P1 finite element scheme for Maxwell's equations with constant permittivity in a boundary neighborhoodComments: This version consists of 34 pages and contains several corrections and adjustments. Most of them were performed in order to conform to referee reports on an extended version of this article, submitted to a journal in May 2019Subjects: Numerical Analysis (math.NA)
- [16] arXiv:1912.08177 (replaced) [pdf, other]
-
Title: Lift & Learn: Physics-informed machine learning for large-scale nonlinear dynamical systemsSubjects: Numerical Analysis (math.NA); Machine Learning (cs.LG)
- [17] arXiv:1804.06245 (replaced) [pdf, other]
-
Title: 3d positive lattice walks and spherical trianglesSubjects: Combinatorics (math.CO); Numerical Analysis (math.NA); Probability (math.PR)
[ showing up to 2000 entries per page: fewer | more ]
Disable MathJax (What is MathJax?)
Links to: arXiv, form interface, find, math, recent, 1912, contact, help (Access key information)