BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Martin Frank

Martin Frank appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAphysics.comp-phmath.APNumerical Analysismath-phmath.MPmath.OCphysics.med-phcond-mat.softcond-mat.stat-mechDistributed, Parallel, and Cluster Computingnucl-th

Catalog footprint

What is connected

19works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

19 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
preprint2016arXiv

An Asymptotic Preserving Two-Dimensional Staggered Grid Method for multiscale transport equations

We propose a two-dimensional asymptotic preserving scheme for linear transport equations with diffusive scalings. It is an extension of the time splitting developed by Jin, Pareschi and Toscani [SINUM,2000], but uses spatial discretizations on staggered grids, which preserves the discrete diffusion limit with a more compact stencil. The first novelty of this paper is that we propose a staggering in two dimensions that requires fewer unknowns than one could have naively expected. The second contribution of this paper is that we rigorously analyze the scheme of Jin, Pareschi, and Toscani [SINUM,2000] We show that the scheme is AP and obtain an explicit CFL condition, which couples a hyperbolic and a parabolic condition. This type of condition is common for asymptotic preserving schemes and guarantees uniform stability with respect to the mean free path. In addition, we obtain an upper bound on the relaxation parameter, which is the crucial parameter of the used time discretization. Several numerical examples are provided to verify the accuracy and asymptotic property of the scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Derivation and Analysis of Lattice Boltzmann Schemes for the Linearized Euler Equations

We derive Lattice Boltzmann (LBM) schemes to solve the Linearized Euler Equations in 1D, 2D, and 3D with the future goal of coupling them to an LBM scheme for Navier Stokes Equations and an Finite Volume scheme for Linearized Euler Equations. The derivation uses the analytical Maxwellian in a BGK model. In this way, we are able to obtain second-order schemes. In addition, we perform an $L^2$-stability analysis. Numerical results validate the approach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fractional Diffusion Limits of Non-Classical Transport Equations

We establish asymptotic diffusion limits of the non-classical transport equation derived in [E. W. Larsen, A generalized Boltzmann equation for non-classical particle transport, Joint international topical meeting on mathematics & computation and supercomputing in nuclear applications, 2007]. By introducing appropriate scaling parameters, the limits will be either regular or fractional diffusion equations depending on the tail behaviour of the path-length distribution. Our analysis uses the Fourier transform combined with a moment method. We conclude with remarks on the diffusion limit of the periodic Lorentz gas equation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Higher order mixed moment approximations for the Fokker-Planck equation in one space dimension

We study mixed-moment models (full zeroth moment, half higher moments) for a Fokker-Planck equation in one space dimension. Mixed-moment minimum-entropy models are known to overcome the zero net-flux problem of full-moment minimum entropy Mn models. Realizability theory for these mixed moments of arbitrary order is derived, as well as a new closure, which we refer to as Kershaw closures. They provide non-negative distribution functions combined with an analytical closure. Numerical tests are performed with standard first-order finite volume schemes and compared with a finite-difference Fokker-Planck scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sensitivity analysis for dose deposition in radiotherapy via a Fokker-Planck model

In this paper we study the sensitivities of electron dose calculations with respect to the stopping power and the transport coefficients. We focus on the application to radiotherapy simulations. We use a Fokker-Planck approximation to the Boltzmann transport equation. Equations for the sensitivities are derived by the adjoint method. The Fokker-Planck equation and its adjoint are solved numerically in slab geometry using the spherical harmonics expansion ($P_N$) and an HLL finite volume method. Our method is verified by comparison to finite difference approximations of the sensitivities. Finally, we present numerical results of the sensitivities for the normalized average dose deposition depth with respect to the stopping power and transport coefficients, demonstrating the increasing relative sensitivities as beam energy decreases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Non-Classical Boltzmann Equation, and Diffusion-Based Approximations to the Boltzmann Equation

We show that several diffusion-based approximations (classical diffusion or SP1, SP2, SP3) to the linear Boltzmann equation can (for an infinite, homogeneous medium) be represented exactly by a non-classical transport equation. As a consequence, we indicate a method to solve diffusion-based approximations to the Boltzmann equation via Monte Carlo, with only statistical errors - no truncation errors.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Perturbed, Entropy-Based Closure for Radiative Transfer

We derive a hierarchy of closures based on perturbations of well-known entropy-based closures; we therefore refer to them as perturbed entropy-based models. Our derivation reveals final equations containing an additional convective and diffusive term which are added to the flux term of the standard closure. We present numerical simulations for the simplest member of the hierarchy, the perturbed M1 or PM1 model, in one spatial dimension. Simulations are performed using a Runge-Kutta discontinuous Galerkin method with special limiters that guarantee the realizability of the moment variables and the positivity of the material temperature. Improvements to the standard M1 model are observed in cases where unphysical shocks develop in the M1 model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Optimal radiotherapy treatment planning using minimum entropy models

We study the problem of finding an optimal radiotherapy treatment plan. A time-dependent Boltzmann particle transport model is used to model the interaction between radiative particles with tissue. This model allows for the modeling of inhomogeneities in the body and allows for anisotropic sources modeling distributed radiation---as in brachytherapy---and external beam sources---as in teletherapy. We study two optimization problems: minimizing the deviation from a spatially-dependent prescribed dose through a quadratic tracking functional; and minimizing the survival of tumor cells through the use of the linear-quadratic model of radiobiological cell response. For each problem, we derive the optimality systems. In order to solve the state and adjoint equations, we use the minimum entropy approximation; the advantages of this method are discussed. Numerical results are then presented.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Adaptive Finite Element Simulation of the Time-dependent Simplified PN Equations

The steady-state simplified $P_N$ approximation to the radiative transport equation has been successfully applied to many problems involving radiation. Recently, time-dependent simplified $P_N$ equations have been derived by an asymptotic analysis similar to the asymptotic derivation of the steady-state $SP_N$ equations \cite{FraKlaLarYas07}. In this paper, we present computational results for the time-dependent $SP_N$ equations in two dimensions, obtained by using an adaptive finite element approach. Several numerical comparisons with other existing models are shown.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Deterministic Partial Differential Equation Model for Dose Calculation in Electron Radiotherapy

Treatment with high energy ionizing radiation is one of the main methods in modern cancer therapy that is in clinical use. During the last decades, two main approaches to dose calculation were used, Monte Carlo simulations and semi-empirical models based on Fermi-Eyges theory. A third way to dose calculation has only recently attracted attention in the medical physics community. This approach is based on the deterministic kinetic equations of radiative transfer. Starting from these, we derive a macroscopic partial differential equation model for electron transport in tissue. This model involves an angular closure in the phase space. It is exact for the free-streaming and the isotropic regime. We solve it numerically by a newly developed HLLC scheme based on [BerCharDub], that exactly preserves key properties of the analytical solution on the discrete level. Several numerical results for test cases from the medical physics literature are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Diffusive Corrections to Pn Approximations

In this paper, we investigate moment methods from a general point of view using an operator notation. This theoretical approach lets us explore the moment closure problem in more detail. This gives rise to a new idea, proposed in [Levermore2005, Levermore2009], of how to improve the well-known Pn approximations. We systematically develop a diffusive correction to the Pn equations from the operator formulation - the so-called Dn approximation. We validate the new approach with numerical examples in one and two dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Optimal treatment planning governed by kinetic equations

In this paper we study a problem in radiotherapy treatment planning, where the evolution of the radiation field is governed by a deterministic Boltzmann transport equation. We show existence, uniqueness and regularity of solutions to an optimal dose distribution problem constrained by the Boltzmann Continuous Slowing-Down equation in an appropriate function space. The main new difficulty is the treatment of the stopping power term. Furthermore, we characterize optimal controls for problems governed by this transport equation.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Brillouin Scattering Study of Propylene Carbonate: An Evaluation of Phenomenological and Mode Coupling Analyses

Brillouin scattering spectra of the molecular glassformer propylene carbonate (PC) in the temperature range 140 K to 350 K were analyzed using both the phenomenological Cole-Davidson memory function and a hybrid memory function consisting of the Cole-Davidson function plus a power-law term representing the critical decay part of the fast beta relaxation. The spectra were also analyzed using the extended two-correlator schematic MCT model recently employed by Gotze and Voigtmann to analyze depolarized light backscattering,dielectric, and neutron-scattering spectra of PC [W.Gotze and Th. Voigtman, Phys. Rev. E61, 4133 (2000). We assess the ability of the phenomenological and MCT fits, each with three free fitting parameters, to simultaneously describe the spectra and give reasonable values for the alpha-relaxation time tau-alpha

0 citations0 reviewsNo code linkedOpen access