BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Martin Frank

Martin Frank contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.comp-phDistributed, Parallel, and Cluster Computingmath-phmath.APmath.MP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2022arXiv

Constructing relaxation systems for lattice Boltzmann methods

We present the first top-down ansatz for constructing lattice Boltzmann methods (LBM) in d dimensions. In particular, we construct a relaxation system (RS) for a given scalar, linear, d-dimensional advection-diffusion equation. Subsequently, the RS is linked to a d-dimensional discrete velocity Boltzmann model (DVBM) on the zeroth and first energy shell. Algebraic characterizations of the equilibrium, the moment space, and the collision operator are carried out. Further, a closed equation form of the RS expresses the added relaxation terms as prefactored higher order derivatives of the conserved quantity. Here, a generalized (2d+1)x(2d+1) RS is linked to a DdQ(2d+1) DVBM which, upon complete discretization, yields an LBM with second order accuracy in space and time. A rigorous convergence result for arbitrary scaling of the RS, the DVBM and conclusively also for the final LBM is proven. The top-down constructed LBM is numerically tested on multiple GPUs with smooth and non-smooth initial data in d=3 dimensions for several grid-normalized non-dimensional numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Neural network-based, structure-preserving entropy closures for the Boltzmann moment system

This work presents neural network based minimal entropy closures for the moment system of the Boltzmann equation, that preserve the inherent structure of the system of partial differential equations, such as entropy dissipation and hyperbolicity. The described method embeds convexity of the moment to entropy map in the neural network approximation to preserve the structure of the minimal entropy closure. Two techniques are used to implement the methods. The first approach approximates the map between moments and the minimal entropy of the moment system and is convex by design. The second approach approximates the map between moments and Lagrange multipliers of the dual of the minimal entropy optimization problem, which present the gradients of the entropy with respect to the moments, and is enforced to be monotonic by introduction of a penalty function. We derive an error bound for the generalization gap of convex neural networks which are trained in Sobolev norm and use the results to construct data sampling methods for neural network training. Numerical experiments are conducted, which show that neural network-based entropy closures provide a significant speedup for kinetic solvers while maintaining a sufficient level of accuracy. The code for the described implementations can be found in the Github repositories.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A low-rank method for two-dimensional time-dependent radiation transport calculations

The low-rank approximation is a complexity reduction technique to approximate a tensor or a matrix with a reduced rank, which has been applied to the simulation of high dimensional problems to reduce the memory required and computational cost. In this work, a dynamical low-rank approximation method is developed for the time-dependent radiation transport equation in 1-D and 2-D Cartesian geometries. Using a finite volume discretization in space and a spherical harmonics basis in angle, we construct a system that evolves on a low-rank manifold via an operator splitting approach. Numerical results on five test problems demonstrate that the low-rank solution requires less memory than solving the full rank equations with the same accuracy. It is furthermore shown that the low-rank algorithm can obtain high-fidelity results at a moderate extra cost by increasing the number of basis functions while keeping the rank fixed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Intrusive acceleration strategies for Uncertainty Quantification for hyperbolic systems of conservation laws

Methods for quantifying the effects of uncertainties in hyperbolic problems can be divided into intrusive and non-intrusive techniques. Non-intrusive methods allow the usage of a given deterministic solver in a black-box manner, while being embarrassingly parallel. However, avoiding intrusive modifications of a given solver takes away the ability to use several inherently intrusive numerical acceleration tools. Moreover, intrusive methods are expected to reach a given accuracy with a smaller number of unknowns compared to non-intrusive techniques. This effect is amplified in settings with high dimensional uncertainty. A downside of intrusive methods is however the need to guarantee hyperbolicity of the resulting moment system. In contrast to stochastic-Galerkin (SG), the Intrusive Polynomial Moment (IPM) method is able to maintain hyperbolicity at the cost of solving an optimization problem in every spatial cell and every time step. In this work, we propose several acceleration techniques for intrusive methods and study their advantages and shortcomings compared to the non-intrusive Stochastic Collocation method. When solving steady problems with IPM, the numerical costs arising from repeatedly solving the IPM optimization problem can be reduced by using concepts from PDE-constrained optimization. Additionally, we propose an adaptive implementation and efficient parallelization strategy of the IPM method. The effectiveness of the proposed adaptations is demonstrated for multi-dimensional uncertainties in fluid dynamics applications, resulting in the observation of requiring a smaller number of unknowns to achieve a given accuracy when using intrusive methods. Furthermore, using the proposed acceleration techniques, our implementation reaches a given accuracy faster than Stochastic Collocation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Massively Parallel Stencil Strategies for Radiation Transport Moment Model Simulations

The radiation transport equation is a mesoscopic equation in high dimensional phase space. Moment methods approximate it via a system of partial differential equations in traditional space-time. One challenge is the high computational intensity due to large vector sizes (1600 components for P39) in each spatial grid point. In this work, we extend the calculable domain size in 3D simulations considerably, by implementing the StaRMAP methodology within the massively parallel HPC framework NAStJA, which is designed to use current supercomputers efficiently. We apply several optimization techniques, including a new memory layout and explicit SIMD vectorization. We showcase a simulation with 200 billion degrees of freedom, and argue how the implementations can be extended and used in many scientific domains.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Structured Approach to the Construction of Stable Linear Lattice-Boltzmann Collision Operators

We introduce a structured approach to the construction of linear BGK-type collision operators ensuring that the resulting Lattice-Boltzmann methods are stable with respect to a weighted $L^2$-norm. The results hold for particular boundary conditions including periodic, bounce-back, and bounce-back with flipping of sign. This construction uses the equivalent moment-space definition of BGK-type collision operators and the notion of stability structures as guiding principle for the choice of the equilibrium moments for those moments influencing the error term only but not the order of consistency. The presented structured approach is then applied to the 3D isothermal linearized Euler equations with non-vanishing background velocity. Finally, convergence results in the strong discrete $L^\infty$-norm highlight the suitability of the structured approach introduced in this manuscript.

0 citations0 reviewsNo code linkedOpen access