Researcher profile

Kolja Brix

Kolja Brix contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
2topics
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

5 published item(s)

preprint2015arXiv

Designing Illumination Lenses and Mirrors by the Numerical Solution of Monge-Ampère Equations

We consider the inverse refractor and the inverse reflector problem. The task is to design a free-form lens or a free-form mirror that, when illuminated by a point light source, produces a given illumination pattern on a target. Both problems can be modeled by strongly nonlinear second-order partial differential equations of Monge-Ampère type. In [Math. Models Methods Appl. Sci. 25 (2015), pp. 803--837, DOI: 10.1142/S0218202515500190] the authors have proposed a B-spline collocation method which has been applied to the inverse reflector problem. Now this approach is extended to the inverse refractor problem. We explain in depth the collocation method and how to handle boundary conditions and constraints. The paper concludes with numerical results of refracting and reflecting optical surfaces and their verification via ray tracing.

preprint2014arXiv

Multilevel Preconditioning of Discontinuous-Galerkin Spectral Element Methods, Part I: Geometrically Conforming Meshes

This paper is concerned with the design, analysis and implementation of preconditioning concepts for spectral Discontinuous Galerkin discretizations of elliptic boundary value problems. While presently known techniques realize a growth of the condition numbers that is logarithmic in the polynomial degrees when all degrees are equal and quadratic otherwise, our main objective is to realize full robustness with respect to arbitrarily large locally varying polynomial degrees degrees, i.e., under mild grading constraints condition numbers stay uniformly bounded with respect to the mesh size and variable degrees. The conceptual foundation of the envisaged preconditioners is the auxiliary space method. The main conceptual ingredients that will be shown in this framework to yield "optimal" preconditioners in the above sense are Legendre-Gauss-Lobatto grids in connection with certain associated anisotropic nested dyadic grids as well as specially adapted wavelet preconditioners for the resulting low order auxiliary problems. Moreover, the preconditioners have a modular form that facilitates somewhat simplified partial realizations. One of the components can, for instance, be conveniently combined with domain decomposition, at the expense though of a logarithmic growth of condition numbers. Our analysis is complemented by quantitative experimental studies of the main components.

preprint2014arXiv

Solving the Monge-Ampère Equations for the Inverse Reflector Problem

The inverse reflector problem arises in geometrical nonimaging optics: Given a light source and a target, the question is how to design a reflecting free-form surface such that a desired light density distribution is generated on the target, e.g., a projected image on a screen. This optical problem can mathematically be understood as a problem of optimal transport and equivalently be expressed by a secondary boundary value problem of the Monge-Ampère equation, which consists of a highly nonlinear partial differential equation of second order and constraints. In our approach the Monge-Ampère equation is numerically solved using a collocation method based on tensor-product B-splines, in which nested iteration techniques are applied to ensure the convergence of the nonlinear solver and to speed up the calculation. In the numerical method special care has to be taken for the constraint: It enters the discrete problem formulation via a Picard-type iteration. Numerical results are presented as well for benchmark problems for the standard Monge-Ampère equation as for the inverse reflector problem for various images. The designed reflector surfaces are validated by a forward simulation using ray tracing.

preprint2013arXiv

Legendre-Gauss-Lobatto grids and associated nested dyadic grids

Legendre-Gauss-Lobatto (LGL) grids play a pivotal role in nodal spectral methods for the numerical solution of partial differential equations. They not only provide efficient high-order quadrature rules, but give also rise to norm equivalences that could eventually lead to efficient preconditioning techniques in high-order methods. Unfortunately, a serious obstruction to fully exploiting the potential of such concepts is the fact that LGL grids of different degree are not nested. This affects, on the one hand, the choice and analysis of suitable auxiliary spaces, when applying the auxiliary space method as a principal preconditioning paradigm, and, on the other hand, the efficient solution of the auxiliary problems. As a central remedy, we consider certain nested hierarchies of dyadic grids of locally comparable mesh size, that are in a certain sense properly associated with the LGL grids. Their actual suitability requires a subtle analysis of such grids which, in turn, relies on a number of refined properties of LGL grids. The central objective of this paper is to derive just these properties. This requires first revisiting properties of close relatives to LGL grids which are subsequently used to develop a refined analysis of LGL grids. These results allow us then to derive the relevant properties of the associated dyadic grids.

preprint2012arXiv

Robust Preconditioners for DG-Discretizations with Arbitrary Polynomial Degrees

Discontinuous Galerkin (DG) methods offer an enormous flexibility regarding local grid refinement and variation of polynomial degrees for a variety of different problem classes. With a focus on diffusion problems, we consider DG discretizations for elliptic boundary value problems, in particular the efficient solution of the linear systems of equations that arise from the Symmetric Interior Penalty DG method. We announce a multi-stage preconditioner which produces uniformly bounded condition numbers and aims at supporting the full flexibility of DG methods under mild grading conditions. The constructions and proofs are detailed in an upcoming series of papers by the authors. Our preconditioner is based on the concept of the auxiliary space method and techniques from spectral element methods such as Legendre-Gauß-Lobatto grids. The presentation for the case of geometrically conforming meshes is complemented by numerical studies that shed some light on constants arising in four basic estimates used in the second stage.