Refine
Year of publication
Document Type
- In Proceedings (92)
- Article (57)
- ZIB-Report (31)
- Book chapter (10)
- Book (2)
- Doctoral Thesis (2)
- In Collection (2)
- Poster (2)
- Habilitation (1)
- Master's Thesis (1)
Is part of the Bibliography
- no (200)
Keywords
Institute
This book promotes the use of mathematical optimization and operations research methods in rail transportation. The editors assembled thirteen contributions from leading scholars to present a unified voice, standardize terminology, and assess the state-of-the-art.
There are three main clusters of articles, corresponding to the classical stages of the planning process: strategic, tactical, and operational. These three clusters are further subdivided into five parts which correspond to the main phases of the railway network planning process: network assessment, capacity planning, timetabling, resource planning, and operational planning. Individual chapters cover:
Simulation
Capacity Assessment
Network Design
Train Routing
Robust Timetabling
Event Scheduling
Track Allocation
Blocking
Shunting
Rolling Stock
Crew Scheduling
Dispatching
Delay Propagation
Back-reflector design in thin-film silicon solar cells by rigorous 3D light propagation modeling
(2014)
Hp-finite-elements for simulating electromagnetic fields in optical devices with rough textures
(2015)
Finite-element simulations of light propagation through subwavelength apertures in metal films
(2009)
Advanced finite-element methods for design and analysis of nano-optical structures: Applications
(2013)
Finite-element based electromagnetic field simulations: Benchmark results for isolated structures
(2013)
3D Finite-Element Simulations of Enhanced Light Transmission Through Arrays of Holes in Metal Films
(2009)
Decomposition of the high dimensional conformational space of bio-molecules into metastable subsets is used for data reduction of long molecular trajectories in order to facilitate chemical analysis and to improve convergence of simulations within these subsets. The metastability is identified by the Perron-cluster cluster analysis of a Markov process that generates the thermodynamic distribution. A necessary prerequisite of this analysis is the discretization of the conformational space. A combinatorial approach via discretization of each degree of freedom will end in the so called ''curse of dimension''. In the following paper we analyze Hybrid Monte Carlo simulations of small, drug-like biomolecules and focus on the dihedral degrees of freedom as indicators of conformational changes. To avoid the ''curse of dimension'', the projection of the underlying Markov operator on each dihedral is analyzed according to its metastability. In each decomposition step of a recursive procedure, those significant dihedrals, which indicate high metastability, are used for further decomposition. The procedure is introduced as part of a hierarchical protocol of simulations at different temperatures. The convergence of simulations within metastable subsets is used as an ''a posteriori'' criterion for a successful identification of metastability. All results are presented with the visualization program AmiraMol.
The paper is motivated by the need for a fast robust adaptive multigrid method to solve complex Helmholtz eigenvalue problems arising from the design of optical chips. A nonlinear multigrid method is developed, which can be regarded as an extension of a previous adaptive Rayleigh quotient minimization method for selfadjoint Helmholtz eigenproblems. Since the complex Helmholtz operator is just a compact nonselfadjoint perturbation of a selfadjoint operator, linear algebra techniques like Schur decomposition can be extended from the finite dimensional case. The efficiency of the derived adaptive nonlinear multigrid method is illustrated by computations for a technologically relevant integrated optics component containing Multi Quantum Well Layers.
{\bf Efficient eigenmode computation for the design of integrated optical chips.}The paper deals with adaptive multigrid methods for 2D Helmholtz eigenvalue problems arising in the design of integrated optical chips. Typical features of the technological problem are its geometric complexity, its multiscale structure, the possible occurrence of eigenvalue clusters, and the necessity of quite stringent required relative error tolerances. For reasons of sheer computational complexity, multigrid methods must be used to solve the discretized eigenvalue problems and adaptive grids must be automatically constructed to avoid an undesirable blow-up of the required number of nodes for these accuracies. In view of the problem specifications, an adaptive multigrid method based on Rayleigh quotient minimization, simultaneous eigenspace iteration, and conjugate gradient method as smoother is carefully selected. Its performance in the numerical simulation of a component of a rather recent optical chip (heterodyne receiver of HHI) is documented.
The paper presents a construction scheme of deriving transparent , i. e. reflection-free, boundary conditions for the numerical solution of Fresnel's equation (being formally equivalent to Schrödinger's equation). These boundary conditions appear to be of a nonlocal Cauchy type. As it turns out, each kind of linear implicit discretization induces its own discrete transparent boundary conditions.
Transparent boundary conditions for a wide-angle approximation of the one-way Helmholtz equation
(2000)
Transparent Boundary Conditions for a Wide-Angle Approximation of the One-Way Helmholtz Equation
(1999)
We present nonlocal discrete transparent boundary conditions for a fourth-order wide-angle approximation of the two-dimensional Helmholtz equation. The boundary conditions are exact in the sense that they supply the same discrete solution on a bounded interior domain as would be obtained by considering the problem on the entire unbounded domain with zero boundary conditions at infinity. The proposed algorithm results in an unconditionally stable propagation method. Numerical examples from optics illustrate the efficiency of our approach.
Boosting the photon-extraction efficiency of nanophotonic structures by deterministic microlenses
(2014)
Reliable numerical simulations of nano-optical structures are the key for design and pre-fabrication processes in diverse disciplines such as lithography in semiconductor industries, spectroscopy of biological molecules, optimization of quantum dot cavities for single-photon sources, and computation of atomic forces like the Casimir Effect.
For these purposes a variety of different methods, e.g. FDTD, FEM and RCWA, are in use. On the one hand, FDTD and FEM are investigated intensively - both in mathematics and numerical experiments - and their approximations and convergence properties are well known. On the other hand, there is a lack of these insights in RCWA. In spite of this, RCWA is commonly used to simulate a wide range of systems.
We review historical and modern contributions to convergence improvements with respect to RCWA from the early sugesstions to modern developments. We study the convergence rates of the open-source software S4 and analyze the algorithmic properties in detail.
Furthermore, we compare RCWA and FEM simulations for different classes of problems including 1D-binary gratings and 2D-periodic photonic crystals.
The introduction of the near-field quantity of optical chirality has emerged in various numerical and few experimental studies of local chirality enhancement due to its relation to the excitation rate of chiral molecules. This time-even pseudoscalar has been dismissed as being a higher-order version of helicity. Nevertheless, we revisit the derivation of the underlying conservation law and define optical chirality in media similar to. We identify the mechanism of chirality conversion by either inhomogeneous or anisotropic space to complement the conservation of optical chirality.
The conservation law of optical chirality in arbitrary space enables the extension of the concept of polarization to the near-field where no distiniguished propagation direction of light is present. We show that the connection of electromagnetic energy and optical chirality provide the ability to define a circular polarization basis in time-harmonic near-field analysis.
In order to illustrate our theory, we present electromagnetic field simulations of simple as well as more complex nanostructures. Results using the well-known far-field polarization concept are readily reproduced and extended from the point of view of chirality conversion.
Determining 2D photonic crystal geometries from reflectance spectra with a reduced basis method
(2016)
Reconstruction of photonic crystal geometries using a reduced basis method for nonlinear outputs
(2016)
Maxwell solvers based on the hp-adaptive finite element method allow for accurate geometrical modeling and high numerical accuracy. These features are indispensable for the optimization of optical properties or reconstruction of parameters through inverse processes. High computational complexity prohibits the evaluation of the solution for many parameters. We present a reduced basis method (RBM) for the time-harmonic electromagnetic scattering problem allowing to compute solutions for a parameter configuration orders of magnitude faster. The RBM allows to evaluate linear and nonlinear outputs of interest like Fourier transform or the enhancement of the electromagnetic field in milliseconds. We apply the RBM to compute light-scattering off two dimensional photonic crystal structures made of silicon and reconstruct geometrical parameters.
Reconstruction of photonic crystal geometries using a reduced basis method for nonlinear outputs
(2016)
Maxwell solvers based on the hp-adaptive finite element method allow for accurate geometrical modeling and high numerical accuracy. These features are indispensable for the optimization of optical properties or reconstruction of parameters through inverse processes. High computational complexity prohibits the evaluation of the solution for many parameters. We present a reduced basis method (RBM) for the time-harmonic electromagnetic scattering problem allowing to compute solutions for a parameter configuration orders of magnitude faster. The RBM allows to evaluate linear and nonlinear outputs of interest like Fourier transform or the enhancement of the electromagnetic field in milliseconds. We apply the RBM to compute light-scattering off two dimensional photonic crystal structures made of silicon and reconstruct geometrical parameters.
Optical 3D simulations in many-query and real-time contexts require new solution strategies. We study an adaptive, error controlled reduced basis method for solving parametrized time-harmonic optical scattering problems. Application fields are, among others, design and optimization problems of nano-optical devices as well as inverse problems for parameter reconstructions occuring e. g. in optical metrology. The reduced basis method presented here relies on a finite element modeling of the scattering problem with parametrization of materials, geometries and sources.
Model order reduction for the time-harmonic Maxwell equation applied to complex nanostructures
(2016)
Fields such as optical metrology and computational lithography require fast and efficient methods for solving the time-harmonic Maxwell's equation. Highly accurate geometrical modelling and numerical accuracy at low computational costs are a prerequisite for any simulation study of complex nano-structured photonic devices. We present a reduced basis method (RBM) for the time-harmonic electromagnetic scattering problem based on the hp-adaptive finite element solver JCMsuite capable of handling geometric and non-geometric parameter dependencies allowing for online evaluations in milliseconds. We apply the RBM to compute light-scattering at optical wavelengths of periodic arrays of fin field-effect transistors (FinFETs) where geometrical properties such as the width and height of the fin and gate can vary in a large range.
Optical 3D simulations in many-query and real-time contexts require new solution strategies. We study an adaptive, error controlled reduced
basis method for solving parametrized time-harmonic optical scattering problems. Application fields are, among others, design and optimization problems
of nano-optical devices as well as inverse problems for parameter reconstructions occurring e. g. in optical metrology. The reduced basis method pre-
sented here relies on a finite element modeling of the scattering problem with
parametrization of materials, geometries and sources.
Model order reduction for the time-harmonic Maxwell equation applied to complex nanostructures
(2016)
Fields such as optical metrology and computational lithography require fast and efficient methods for solving
the time-harmonic Maxwell’s equation. Highly accurate geometrical modeling and numerical accuracy atcomputational costs are a prerequisite for any simulation study of complex nano-structured photonic devices.
We present a reduced basis method (RBM) for the time-harmonic electromagnetic scattering problem based
on the hp-adaptive finite element solver JCMsuite capable of handling geometric and non-geometric parameter
dependencies allowing for online evaluations in milliseconds. We apply the RBM to compute light-scatteringoptical wavelengths off periodic arrays of fin field-effect transistors (FinFETs) where geometrical properties such
as the width and height of the fin and gate can vary in a large range.
Rigorous optical simulations of 3-dimensional nano-photonic structures are an important tool in the analysis and optimization of scattering properties of nano-photonic devices or parameter reconstruction. To construct geometrically accurate models of complex structured nano-photonic devices the finite element method (FEM) is ideally suited due to its flexibility in the geometrical modeling and superior convergence properties. Reduced order models such as the reduced basis method (RBM) allow to construct self-adaptive, error-controlled, very low dimensional approximations for input-output relationships which can be evaluated orders of magnitude faster than the full model. This is advantageous in applications requiring the solution of Maxwell's equations for multiple parameters or a single parameter but in real time. We present a reduced basis method for 3D Maxwell's equations based on the finite element method which allows variations of geometric as well as material and frequency parameters. We demonstrate accuracy and efficiency of the method for a light scattering problem exhibiting a resonance in the electric field.
In many experimentally realized applications, e.g. photonic crystals, solar cells and light-emitting diodes, nano-photonic systems are coupled to a thick substrate layer, which in certain cases has to be included as a part of the optical system. The finite element method (FEM) yields rigorous, high accuracy solutions of full 3D vectorial Maxwell's equations [1] and allows for great flexibility and accuracy in the geometrical modelling. Time-harmonic FEM solvers have been combined with Fourier methods in domain decomposition algorithms to compute coherent solutions of these coupled system. [2,3] The basic idea of a domain decomposition approach lies in a decomposition of the domain into smaller subdomains, separate calculations of the solutions and coupling of these solutions on adjacent subdomains.
In experiments light sources are often not perfectly monochromatic and hence a comparision to simulation results might only be justified if the simulation results, which include interference patterns in the substrate, are spectrally averaged.
In this contribution we present a scattering matrix domain decomposition algorithm for Maxwell's equations based on FEM. We study its convergence and advantages in the context of optical simulations of silicon thin film multi-junction solar cells. This allows for substrate light-trapping to be included in optical simulations and leads to a more realistic estimation of light path enhancement factors in thin-film devices near the band edge.
Adaptive sampling strategies for efficient parameter scans in nano-photonic device simulations
(2014)
Rigorous optical simulations are an important tool in optimizing scattering properties of nano-photonic devices and are used, for example, in solar cell optimization. The finite element method (FEM) yields rigorous, time-harmonic, high accuracy solutions of the full 3D vectorial Maxwell's equations [1] and furthermore allows for great flexibility and accuracy in the geometrical modeling of these often complex shaped 3D nano-structures. A major drawback of frequency domain methods is the limitation of single frequency evaluations. For example the accurate computation of the short circuit current density of an amorphous silicon / micro-crystalline multi-junction thin film solar cell may require the solution of Maxwell's equations for over a hundred different wavelengths if an equidistant sampling strategy is employed. Also in optical metrology, wavelength scans are frequently used to reconstruct unknown geometrical and material properties of optical systems numerically from measured
scatterometric data.
In our contribution we present several adaptive numerical integration and sampling routines and study their efficiency in the context of the determination of generation rate profiles of solar cells. We show that these strategies lead to a reduction in the computational effort without loss of accuracy. We discuss the employment of tangential information in a Hermite interpolation scheme to achieve similar accuracy on coarser grids. We explore the usability of these strategies for scatterometry and solar cell simulations.
Adaptive sampling strategies for effcient parameter scans in nano-photonic device simulations
(2014)
The aim of this paper is to develop fast methods for the solution of nonlinear Schrödinger type equations in fiber optics. Using the method of lines we have to solve a stiff system of ordinary differential equations where the eigenvalues of the Jacobian are close to the imaginary axis. This is usually done by a Split Step method. Here we consider the extrapolation of Split Step methods with adaptive order and step size control. For more complicated nonlinearities, in particular stimulated Raman scattering, Split Step methods are less efficient since symmetry is either destroyed or requires much additional effort. In this case we use implicit Runge Kutta formulas of Gauß type. The key point for the efficient implementation of these methods is that the system of nonlinear algebraic equations can be solved without setting up the Jacobian. The proposed methods are compared to other methods, in particular exponential integrators, the method of Marcuse, and the method of Blow and Wood.
Solving time-harmonic scattering problems based on the pole condition: Convergence of the PML method
(2001)
In this paper we study the PML method for Helmholtz-type scattering problems with radially symmetric potential. The PML method consists in surrounding the computational domain by a \textbf{P}erfectly \textbf{M}atched sponge \textbf{L}ayer. We prove that the approximate solution obtained by the PML method converges exponentially fast to the true solution in the computational domain as the thickness of the sponge layer tends to infinity. This is a generalization of results by Lassas and Somersalo based on boundary integral eqaution techniques. Here we use techniques based on the pole condition instead. This makes it possible to treat problems without an explicitly known fundamental solution.
We present a new efficient algorithm for the solution of direct time-harmonic scattering problems based on the Laplace transform. This method does not rely on an explicit knowledge of a Green function or a series representation of the solution, and it can be used for the solution of problems with radially symmetric potentials and problems with waveguides. The starting point is an alternative characterization of outgoing waves called \emph{pole condition}, which is equivalent to Sommerfeld's radiation condition for problems with radially symmetric potentials. We obtain a new representation formula, which can be used for a numerical evaluation of the exterior field in a postprocessing step. Based on previous theoretical studies, we discuss the numerical realization of our algorithm and compare its performance to the PML method.
The pole condition is a general concept for the theoretical analysis and the numerical solution of a variety of wave propagation problems. It says that the Laplace transform of the physical solution in radial direction has no poles in the lower complex half-plane. In the present paper we show that for the Helmholtz equation with a radially symmetric potential the pole condition is equivalent to Sommerfeld's radiation condition. Moreover, a new representation formula based on the pole condition is derived and used to prove existence, uniqueness and asymptotic properties of solutions. This lays the foundations of a promising new algorithm to solve time-harmonic scattering problems numerically and provides a new approach for analyzing existing algorithms such as the Perfectly Matched Layer (PML) method and the Bayliss-Gunzburger-Turkel (BGT) algorithm.