@inbook{ef6d4e1030124b6bbbe357a9de814e9d,
title = "An Error-Based Low-Rank Correction for Pressure Schur Complement Preconditioners",
author = "R.S. Beddig and J. Behrens and S.L. Borne and K. Simon",
year = "2023",
doi = "10.1007/978-3-031-45158-4_5",
language = "English",
booktitle = "Lecture Notes in Computational Science and Engineering",
}
@article{9945020b2f94444dab6aacda3a21965f,
title = "Performance Assessment of the Cloud for Prototypical Instant Computing Approaches in Geoscientific Hazard Simulations",
abstract = "Computing forecasts of hazards, such as tsunamis, requires fast reaction times and high precision, which in turn demands for large computing facilities that are needed only in rare occasions. Cloud computing environments allow to configure largely scalable on-demand computing environments. In this study, we tested two of the major cloud computing environments for parallel scalability for relevant prototypical applications. These applications solve stationary and non-stationary partial differential equations by means of finite differences and finite elements. These test cases demonstrate the capacity of cloud computing environments to provide scalable computing power for typical tasks in geophysical applications. As a proof-of-concept example of an instant computing application for geohazards, we propose a workflow and prototypical implementation for tsunami forecasting in the cloud. We demonstrate that minimal on-site computing resources are necessary for such a forecasting environment. We conclude by outlining the additional steps necessary to implement an operational tsunami forecasting cloud service, considering availability and cost. Copyright {\textcopyright} 2022 Behrens, Schulz and Simon.",
keywords = "cloud computing, instant computing, natural hazard, parallel performance, tsunami",
author = "J. Behrens and Arne Schulz and K. Simon",
note = "Export Date: 20 April 2022 Correspondence Address: Behrens, J.; Department of Mathematics/CEN, Germany; email: joern.behrens@uni-hamburg.de Funding details: 603839 Funding details: Deutsche Forschungsgemeinschaft, DFG, 390683824 Funding details: Universit{\"a}t Hamburg, UH Funding text 1: Parts of this research were conducted in the framework of the ASTARTE project with funding from the European Unions Seventh Program for research, technological development and demonstration under grant agreement No. 603839. Additional funding was obtained by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany{\textquoteright}s Excellence Strategy—EXC 2037 {\textquoteright}CLICCS—Climate, Climatic Change, and Society{\textquoteright}—Project Number: 390683824, contribution to the Center for Earth System Research and Sustainability (CEN) of Universit{\"a}t Hamburg.",
year = "2022",
month = mar,
day = "9",
doi = "10.3389/feart.2022.762768",
language = "English",
volume = "10",
journal = "Frontiers in Earth Science",
issn = "2296-6463",
publisher = "Frontiers Media S. A.",
}
@article{725a4721ba754b13b543569b2ba773a3,
title = "Extending legacy climate models by adaptive mesh refinement for single-component tracer transport: a case study with ECHAM6-HAMMOZ (ECHAM6.3-HAM2.3-MOZ1.0)",
abstract = "The model error in climate models depends on mesh resolution, among other factors. While global refinement of the computational mesh is often not feasible computationally, adaptive mesh refinement (AMR) can be an option for spatially localized features. Creating a climate model with AMR has been prohibitive so far. We use AMR in one single-model component, namely the tracer transport scheme.Particularly, we integrate AMR into the tracer transport module of the atmospheric model ECHAM6 and test our implementation in several idealized scenarios and in a realistic application scenario (dust transport). To achieve this goal, we modify the flux-form semi-Lagrangian (FFSL) transport scheme in ECHAM6 such that we can use it on adaptive meshes while retaining all important properties (such as mass conservation) of the original FFSL implementation. Our proposed AMR scheme is dimensionally split and ensures that high-resolution information is always propagated on (locally) highly resolved meshes. We utilize a data structure that can accommodate an adaptive Gaussian grid.We demonstrate that our AMR scheme improves both accuracy and efficiency compared to the original FFSL scheme. More importantly, our approach improves the representation of transport processes in ECHAM6 for coarse-resolution simulations. Hence, this paper suggests that we can overcome the overhead of developing a fully adaptive Earth system model by integrating AMR into single components while leaving data structures of the dynamical core untouched. This enables studies to retain well-tested and complex legacy code of existing models while still improving the accuracy of specific components without sacrificing efficiency.",
author = "Yumeng Chen and Konrad Simon and J{\"o}rn Behrens",
year = "2021",
month = may,
day = "3",
doi = "10.5194/gmd-14-2289-2021",
language = "English",
volume = "14",
pages = "2289–2316",
journal = "Geoscientific Model Development",
issn = "1991-9603",
publisher = "Copernicus Publications",
number = "5",
}
@article{2911fda8d5e4482c95447fdc85e97132,
title = "Semi-Lagrangian Subgrid Reconstruction for Advection-Dominant Multiscale Problems with Rough Data",
abstract = "We introduce a new framework of numerical multiscale methods for advection-dominated problems motivated by climate sciences. Current numerical multiscale methods (MsFEM) work well on stationary elliptic problems but have difficulties when the model involves dominant lower order terms. Our idea to overcome the associated difficulties is a semi-Lagrangian based reconstruction of subgrid variability into a multiscale basis by solving many local inverse problems. Globally the method looks like a Eulerian method with multiscale stabilized basis. We show example runs in one and two dimensions and a comparison to standard methods to support our ideas and discuss possible extensions to other types of Galerkin methods, higher dimensions and nonlinear problems.",
keywords = "Advection–, Inverse problems, Multiscale finite elements, Multiscale simulation, Semi-Lagrangian, diffusion",
author = "Konrad Simon and J{\"o}rn Behrens",
year = "2021",
month = mar,
day = "28",
doi = "10.1007/s10915-021-01451-w",
language = "English",
volume = "87",
journal = "Journal of Scientific Computing",
issn = "0885-7474",
publisher = "Springer New York",
number = "2",
}
@inbook{2543a80d2a1b4d428a04ffbf58af1cb0,
title = "Massively Parallel Multiscale Simulations of the Feedback of Urban Canopies",
abstract = "Urban canopies consist of buildings and trees that are aligned along a street in the horizontal direction. These canopies in cities and forests modulate the local climate considerably in a complex way. Canopies constitute very fine subgrid features that actually have a significant impact on other components of earth system models but their feedbacks on larger scales are by now represented in rather heuristic ways. The problem in simulating their impact is twofold: First, their local modeling is delicate and, secondly, the numerical modeling of the scale interaction between fine and large scales is complicated since the fine scale structure is global. We will mostly focus on the second aspect.Multiscale finite element methods (MsFEM) in their classical form have been applied to various porous media problems but the situation in climate, and hence flow-dominated regimes is different from porous media applications. In order to study the effect of various parameters like the concentration of pollutants, or the dynamics of the background velocity and of the temperature in the atmospheric boundary layer, a semi-Lagrangian reconstruction based multiscale finite element framework (SLMsR) developed by [1, 2] for passive tracer transport modeled by an advection-diffusion equation with high-contrast oscillatory diffusion is applied.These methods are composed of two parts: a local-in-time semi-Lagrangian offline phase that pre-computes basis functions and an online phase that uses these basis functions to compute the solution on a coarse Eulerian simulation mesh. The overhead of pre-computing the basis functions in each coarse block can further be reduced by parallelization. The online phase is approximately as fast as a low resolution standard FEM but using the modified basis that carries subgrid information still allows to reveal the fine scale features of a highly resolved solution and is therefore accurate. This approach is studied in order to reveal the feedback of processes in the canopy layer on different scales present in climate simulation models and in particular on the atmospheric boundary layer.We will show the results of massively parallel simulations for passive tracer transport in an urban region using the new multiscale approach and compare them to classical approaches.",
author = "Heena Patel and Konrad Simon and J{\"o}rn Behrens",
year = "2021",
month = mar,
day = "3",
doi = "10.5194/egusphere-egu21-2507",
language = "English",
series = "EGU General Assembly ",
booktitle = "EGU General Assembly 2021",
note = "EGU General Assembly 2021 ; Conference date: 19-04-2021 Through 30-04-2021",
}
@inbook{cd8e7770388e47d3a35ef56307cc9a3e,
title = "Stable Multiscale Discretizations of L2-Differential Complexes",
abstract = "Global simulations over long time scales in climate sciences often require coarse grids due to computational constraints. This leaves dynamically important smaller scales unresolved. Thus the influence of small scale processes has to be taken care of by different means. State-of-the-art dynamical cores represent the influence of subscale processes typically via subscale parametrizations and often employ heuristic coupling of scales. This, however, unfortunately often lacks mathematical consistency. The aim of this work is to improve mathematical consistency of the upscaling process that transfers information from the subgrid to the coarse scales of the dynamical core and to largely extend the idea of adding subgrid correctors to basis functions for scalar and vector valued elements discretizing various function spaces.Discussing prototypically the issue of weighted Hodge decompositions I will show that standard techniques on coarse meshes fail to find good projections in all parts of a modified de Rham complex if rough data is involved and discuss an idea of how to construct multiscale finite element (MsFEM) correctors to scalar and vector valued finite elements and, further, how to construct stable multiscale element pairings using the theory of finite element exterior calculus (FEEC). This can be seen as a meta-framework that contains the construction of standard MsFEMs [Efendiev2009, Graham2012]. Application examples here comprise porous media, elasticity, and fluid flow as well as electromagnetism in fine-scale and high-contrast media. I will provide the necessary theoretical background in homological algebra and differential geometry, and discuss a scalable MPI based implementation technique suitable for large clusters. Several computational examples will be shown. I may, if time permits, discuss some ideas from homogenisation theory to attack the problem of a proof of accuracy.",
author = "Konrad Simon and J{\"o}rn Behrens",
year = "2020",
month = mar,
day = "10",
doi = "10.5194/egusphere-egu2020-14970",
language = "English",
booktitle = "EGU General Assembly 2020",
note = "EGU General Assembly 2020 ; Conference date: 04-05-2020 Through 08-05-2020",
}
@article{b2ff2f8660b3425c8b1b0d869f7c5967,
title = "Multiscale Finite Elements for Transient Advection-Diffusion Equations through Advection-Induced Coordinates",
abstract = "Long simulation times in climate science typically require coarse grids due to computational constraints. Nonetheless, unresolved subscale information significantly influences the prognostic variables and cannot be neglected for reliable long-term simulations. This is typically done via parametrizations, but their coupling to the coarse grid variables often involves simple heuristics. We explore a novel upscaling approach inspired by multiscale finite element methods. These methods are well established in porous media applications, where mostly stationary or quasi stationary situations prevail. In advection-dominated problems arising in climate simulations, the approach needs to be adjusted. We do so by performing coordinate transforms that make the effect of transport milder in the vicinity of coarse element boundaries. The idea of our method is quite general, and we demonstrate it as a proof-of-concept on a one-dimensional passive advection-diffusion equation with oscillatory background velocity and diffusion.",
keywords = "Advection-diffusion equation, Finite element method, Multiscale method, Subscale parametrization, Upscaling",
author = "Konrad Simon and J{\"o}rn Behrens",
year = "2020",
month = jan,
day = "1",
doi = "10.1137/18M117248X",
language = "English",
volume = "18",
pages = "543--571",
journal = "Multiscale Modeling Simulation",
publisher = "Society for Industrial and Applied Mathematics Publications",
number = "2",
}
@inbook{065e1c44ea96440f90d9b35721831ec1,
title = "A Semi-Lagrangian Multiscale Framework for Advection-Dominant Problems",
abstract = "We introduce a new parallelizable numerical multiscale method for advection-dominated problems as they often occur in engineering and geosciences. State of the art multiscale simulation methods work well in situations in which stationary and elliptic scenarios prevail but are prone to fail when the model involves dominant lower order terms which is common in applications. We suggest to overcome the associated difficulties through a reconstruction of subgrid variations into a modified basis by solving many independent (local) inverse problems that are constructed in a semi-Lagrangian step. Globally the method looks like a Eulerian method with multiscale stabilized basis. The method is extensible to other types of Galerkin methods, higher dimensions, nonlinear problems and can potentially work with real data. We provide examples inspired by tracer transport in climate systems in one and two dimensions and numerically compare our method to standard methods.",
author = "Konrad Simon and J{\"o}rn Behrens",
year = "2019",
doi = "10.1007/978-3-030-22747-0_30",
language = "English",
isbn = "978-3-030-22746-3",
series = "Lecture Notes in Computer Science",
publisher = "Springer International Publishing",
pages = "393--409",
editor = "Rodrigues, {Jo{\~a}o M. F.} and Cardoso, {Pedro J. S.} and J{\^a}nio Monteiro and Roberto Lam and Krzhizhanovskaya, {Valeria V.} and Lees, {Michael H.} and Dongarra, {Jack J.} and Sloot, {Peter M.A.}",
booktitle = "Computational Science -- ICCS 2019",
address = "Switzerland",
}
@inbook{f110a0b36a034587beba6353e48aad38,
title = "Comparison of various coupling methods for 1D diffusion equations with a analytical solution of two phase Stefan problem",
abstract = "We implemented loose and tight coupling methods to understand thermal diffusion between ocean and ice by means of a simplified one-dimensional model set-up proposed by Stefan. A Stefan problem is a prototypical two-phase model that can used to model, for example, melting and freezing of water due to the transfer of heat fluxes between the two phases. We discretized heat fluxes using low order derivatives for loose coupling and higher order derivatives for tight coupling while fluxes are computed at the (moving) interface. Compared to a known reference solution the tight coupling method exhibits a lower error when compared to the loose coupling discretization. However, further numerical tests are required to analyze these coupling methods.",
author = "Anusha Sunkisala and Konrad Simon and J{\"o}rn Behrens",
year = "2019",
language = "English",
isbn = "978-84-949194-5-9",
pages = "709–718",
booktitle = "COUPLED VIII : proceedings of the VIII International Conference on Computational Methods for Coupled Problems in Science and Engineering",
publisher = "CIMNE",
}
@article{5dcd135026df4fdb81075d7e18968653,
title = "A Hyperelastic Two-Scale Optimization Model for Shape Matching",
abstract = " We suggest a novel shape matching algorithm for three-dimensional surface meshes of disk or sphere topology. The method is based on the physical theory of nonlinear elasticity and can hence handle large rotations and deformations. Deformation boundary conditions that supplement the underlying equations are usually unknown. Given an initial guess, these are optimized such that the mechanical boundary forces that are responsible for the deformation are of a simple nature. We show a heuristic way to approximate the nonlinear optimization problem by a sequence of convex problems using finite elements. The deformation cost, i.e, the forces, is measured on a coarse scale while ICP-like matching is done on the fine scale. We demonstrate the plausibility of our algorithm on examples taken from different datasets. ",
keywords = "cs.CG, cs.CV, cs.GR",
author = "Konrad Simon and Sameer Sheorey and David Jacobs and Ronen Basri",
year = "2015",
month = jul,
day = "28",
language = "English",
journal = "arXiv.org",
publisher = "Cornell Univ. Library",
}