Research reports

Years: 2024 2023 2022 2021 2020 2019 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2004 2003 2002 2001 2000 1999 1998 1997 1996 1995 1994 1993 1992 1991

Adaptive Anisotropic Petrov-Galerkin Methods for First Order Transport Equations

by W. Dahmen and G. Kutyniok and W. Lim and Ch. Schwab and G. Welper

(Report number 2016-03)

Abstract
This paper builds on recent developments of adaptive methods for linear transport equations based on certain stable variational formulations of Petrov-Galerkin type. The key issues can be summarized as follows. The variational formulations allow us to employ meshes with cells of arbitrary aspect ratios. We develop a refinement scheme generating highly anisotropic partitions that is inspired by shearlet systems. We establish approximation rates for \(N\)-term approximations from corresponding piecewise polynomials for certain compact cartoon classes of functions. In contrast to earlier results in a curvelet or shearlet context the cartoon classes are concisely defined through certain characteristic parameters and the dependence of the approximation rates on these parameters is made explicit here. The approximation rate results serve then as a benchmark for subsequent applications to adaptive Galerkin solvers for transport equations. We outline a new class of directionally adaptive, Petrov-Galerkin discretizations for such equations. In numerical experiments, the new algorithms track \(C^2\)-curved shear layers and discontinuities stably and accurately, and realize essentially optimal rates. Finally, we treat parameter dependent transport problems, which arise in kinetic models as well as in radiative transfer. In heterogeneous media these problems feature propagation of singularities along curved characteristics precluding, in particular, fast marching methods based on ray-tracing. Since now the solutions are functions of spatial variables and parameters one has to address the curse of dimensionality. We show computationally, for a model parametric transport problem in heterogeneous media in \(2+1\) dimension, that sparse tensorization of the presently proposed spatial directionally adaptive scheme with hierarchic collocation in ordinate space based on a stable variational formulation high-dimensional phase space, the curse of dimensionality can be removed when approximating averaged bulk quantities.

Keywords: Linear transport problems, $L_2$-stable Petrov-Galerkin formulations, $\delta$-proximality, adaptive refinements, anisotropic discretizations, best $N$-term approximation

@Techreport{DKLSW16_640,
  author = {W. Dahmen and G. Kutyniok and W. Lim and Ch. Schwab and G. Welper},
  title = {Adaptive Anisotropic Petrov-Galerkin Methods for First Order Transport Equations},
  institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
  number = {2016-03},
  address = {Switzerland},
  url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2016/2016-03.pdf },
  year = {2016}
}

Disclaimer
© Copyright for documents on this server remains with the authors. Copies of these documents made by electronic or mechanical means including information storage and retrieval systems, may only be employed for personal use. The administrators respectfully request that authors inform them when any paper is published to avoid copyright infringement. Note that unauthorised copying of copyright material is illegal and may lead to prosecution. Neither the administrators nor the Seminar for Applied Mathematics (SAM) accept any liability in this respect. The most recent version of a SAM report may differ in formatting and style from published journal version. Do reference the published version if possible (see SAM Publications).