Mathematics for Sustainable Energies
The last years have brought a fundamental political change towards energy management in Germany: As soon as possible regenerative energies on all scales are to replace nuclear as well as fossil fuel energies. Electromobility is changing traffic profiles, in particular urban traffic. Decentralized compact power stations that save one portion of the inevitable thermodynamic degree of efficiency by consuming the heat at the generation site will have to be established and optimally linked. Apart from these individual topics, the overall topic of energy efficiency has gained a dominant importance. All of these topics induce numerous new quantitative problems of high complexity, which are the source of new mathematical challenges.
Photovoltaics is one of the fastest growing regenerative energy sectors. In the focus of present interest are demands of higher efficiency with lower material usage and production costs.
In the Action Plan for Electromobility 2020, Berlin sets the goal of becoming Europe's leading location for electromobility. As for the technological basis, the key words are fuel cells, novel batteries or plug-in hybrids. Of course, urban traffic will be the topic of main interest for Berlin.
Decentralized power station networks. In the years to come, an increasing number of buildings will produce energy in a regenerative way, e.g., via solar cell roofs and walls. Whenever the produced electrical energy is not needed at the homes, the locally available electrical power will have to be fed into the communal electrical network systems. This gives rise to multiple problems of energy distribution and optimal network management, problems of utmost mathematical and computational complexity.
Projects
Dr. Christian Schröder
Description
One of the bottlenecks of current procedures for the generation and distribution of green (wind or solar) energy is the accurate and timely simulation of processes in the ocean and atmosphere that can be used in short term planning and real time control of energy systems. A particular difficulty is the real time construction of physically plausible model initializations and 'controls/inputs' to bring simulations into coherence with available observations when observation locations and observations are coming in at variable times and locations.The currently best approach for fixed observation times and locations are variational data assimilation techniques. These methods use a four dimensional model that is adapted to the incoming observations using a combination of different filtering techniques and numerical integration of the dynamical system. In order to make these methods efficient in real time data assimilation they have to be combined with appropriate model order reduction methods. A major difficulty in these techniques is the combination of approximate transfer functions and approximate initial and boundary conditions as well as the construction of guaranteed error estimates and the capturing of essential features of the original model. The so-called representer approach formulates the data assimilation problem as the numerical solution of a large-scale nonlinear optimal control problem and incorporates the assimilation of the model to the observations, via an extended ensemble Kalman filter, and the adaptation of the initial data in one approach. Adding further assumptions and linearization this optimization problem usually reduces to a linear quadratic optimal control problem which is solved via the solution of a boundary value problem with Hamiltonian structure.
Website
X Close projectProf. Dr. Volker Mehrmann
Description
In the field of sustainable energies, microbial cell factories such as yeasts and cyanobacteria are receiving increasing interest due to their potential to produce biofuels. A major question is how the metabolism of these microorganisms is coordinated in a dynamic environment such that the correct macromolecules are synthesized at the right time in order to enable growth and survival. Recent mathematical modeling approaches have made it possible to study this problem using an optimal dynamic resource allocation formalism such as dynamic enzyme-cost flux balance analysis (deFBA). The goal of this project is to study the mathematical properties properties of the underlying optimal control problem involving differential-algebraic constraints and to develop efficient numerical solution strategies.Website
X Close projectProf. Dr. Alexander Mielke
Description
The aim of the project D-SE2 is to find adequate spatially resolved PDE models for the electrothermal description of organic semiconductor devices describing self-heating and thermal switching phenomena. Moreover, the project intends to investigate their analytical properties, derive suitable numerical approximation schemes, and provide simulation results which can help to optimize large-area organic light emitting diodes.Click here for more information
Website
X Close projectProf. Dr. Peter Karl Friz
Dr. Clemens Guhlke
Dr. Manuel Landstorfer
Description
Currently lithium-ion batteries are the most promising storage devices to store and convert chemical energy into electrical energy. An important class of modern lithium batteries contain electrodes that consist of many nano-particles. During the charging process of a battery, lithium is reversibly stored in the ensemble of the nano-particles and the particles undergo a phase transition from a Li-rich to a Li-poor phase. For this type of batteries a successful mathematical model was developed in the previous ECMath project SE8, based on a stochastic mean field interacting particle system. The new project focuses on modeling, analysis and simulations of extreme conditions in battery operation like fast charging, mostly full/empty discharge states, mechanical stresses within the electrode. The aim of the project is to achieve deeper understanding of the behavior of lithium-ion batteries in extreme conditions.Website
X Close projectProf. Dr. Volker Mehrmann
Prof. Dr. Caren Tischendorf
Description
In the project the stability of power networks and power network models is analyzed. The classical way of modeling a power network is via a large differential-algebraic system of network equations (DAE). Modifications of the power network by adding extra power lines into the network grid or by removing some power lines can be interpreted as low rank perturbations of matrices and matrix pencils that linearize the DAE system mentioned above. In the project, the influence of these perturbation on the stability of the network is analyzed.Website
X Close projectDr. Matthias Liero
Description
Organic materials lead to innovative electronic components with fine-tuned properties and promise more sustainable, eco-friendly electronic technologies. The potential for greater sustainability extends across the entire life cycle of organic electronics, beginning with the use of materials that are synthesized rather than mined from the earth, over low temperature production of devices, and ending with potentially biodegradable or recyclable devices. The aim of the project is to develop a thermodynamically correct energy-drift-diffusion model for organic semiconductor devices and its discretization and implementation in a simulation tool. First, the transport of charge carriers in the isothermal case will be described on the basis of a drift-diffusion model, taking the distinctive features of organic materials into account. Second, the model will be extended in a thermodynamically consistent way to include the self-heating and the resulting feedback as well. In both points, the aspects of modeling, analysis, numerics, and simulation are considered. There are several mathematical challenges regarding a drift-diffusion description of organic devices: In organic semiconductors, the energy levels are Gaussian-distributed with disorder parameter σ such that the densities of electrons and holes are described by the Gauss-Fermi integrals. This leads to a generalized Einstein relation and results in a nonlinear diffusion enhancement in the relation between drift and diffusion current. Moreover, the mobility functions for organic semiconductor materials with Gaussian disorder are increasing with respect to temperature, carrier density, and electrical field strength. The outcome of the project is a fundamental building block for a more efficient multi-scale and multi-physics description and simulation of organic devices.Website
X Close projectProf. Dr. Gitta Kutyniok
Prof. Dr. Barbara Wagner
Dr. Arne Roggensack
Description
The project SE4 aims to develop and study mathematical models in order to understand, functionalize and optimize modern nanostructured materials. Such materials are fundamental for the design of next generation thin-film solar cells as well as batteries for the production and storage of sustainable energy, respectively. Besides the mathematical modeling, the main goals of this research project are the analysis of the developed phase field systems and the construction of numerical algorithms that efficiently capture the material properties and, in particular, their anisotropic nature. More information...Website
X Close projectDescription
The estimation of the temperature and airflow distribution in buildings via the location of sensor networks is a nonlinear and multiscale problem where dynamics and measurements are (in general) stochastically perturbed. The goal here is to reliably outline temperature and airflow in certain areas by placing sensors on prescribed admissible locations while optimizing several criteria. A trustworthy estimation, provided that closed-loop controllers are in place, becomes a main step in reducing the energy consumption of such control systems.Website
X Close projectDescription
Solar energy is mostly harvested by means of photovoltaic (PV) or concentrating photovoltaic (CPV) solar cells. The efficiency of CPV is higher (at least twice) than the traditional PV but significantly more expensive. To reduce costs, optical condensers (e.g., a Fresnel lens) to concentrate solar light on each CPV cell are used. Moreover, since the energy production is maximized when the panels are perpendicular to the light beam, mechanical tracking systems that move the array of solar panels based on the position of the sun. But these tracking system increases costs, requires power and are error-prone. The goal of this project is the optimal design and control of steerers and concentrators for PV or CPV using electrowetting (EW) and electrowetting-on-dielectric (EWOD).Website
X Close projectDescription
The performance of optoelectronic devices, such as photovoltaic cells, is critically influenced by the transport of excitonic energy because the majority of photo-excitations occur in the bulk of the crystal from where the energy has to be transported to the interfaces with the electrodes, where charge generation often only occurs. In organic semiconductors, e.g. molecular crystals, polymer chains or dendrimers, the excitons are strongly localized, and the energy transport is normally modeled in terms of excitons diffusively hopping between sites. The present proposal aims at an improved understanding of the excitonic energy transport in organic semiconductors, which is relevant for the characterization of organic solar cells, on a microscopic basis, with emphasis on the role of the electron-phonon coupling. Using mixed quantum-classical dynamics schemes, the electronic degrees of freedom (excitons) are to be treated quantum-mechanically while the nuclear motions (phonons, molecular vibrations) are treated classically. To this end, stochastic surface-hopping algorithms shall be applied and further developed.Website
X Close projectProf. Dr. Frank Schmidt
Dr. Martin Hammerschmidt
Description
Artificial photosynthesis and water splitting, i.e. the sustainable production of chemical fuels like hydrogen and carbohydrates from water and carbon dioxide, has the potential to store the abundance of solar energy that reaches the earth in chemical bonds. Fundamental in this process is the conversion of electromagnetic energy. In photoelectrochemical water splitting semiconductor materials are employed to generate electron hole pairs with sufficient energy to drive the electrochemical reactions. In this project we investigate the use of metallic nanoparticles to excite plasmonic resonances by means of numerical simulations. These resonances localize electromagnetic nearfields which is beneficial for the electrochemical reactions. We develop electromagnetic models and numerical methods to facilite in depth analysis of these processes in close contact with our collaboration partners within the ECMath and the joint lab ``Berlin Joint Lab for Optical Simulation for renewable Energy research'' (BerOSE) between the ZIB, FU and HZB.Website
X Close projectProf. Dr. Volker Mehrmann
Dr. Matthias Voigt
Description
In this project we will study the modeling of power networks by employing the port-Hamiltonian framework. Energy based modeling with port-Hamiltonian descriptor systems has many advantages, e. g., it accounts for the physical interpretation of its variables, it is best suited for the modular structure of the network, since coupled port-Hamiltonian systems form again a port-Hamiltonian system and it encodes these properties in algebraic and geometric properties that simplify Galerkin type model reduction, stability analysis, and also efficient discretization techniques. To improve the predictions that one obtains from such models we suggest to employ data assimilation and state estimation techniques by incorporating the measurement data. These would allow to take the uncertainty in the measurements and the presence of unmodeled dynamics as well as data and modeling errors into account. The improved predictions can then be used to control the network such that (the expected value of) the load is kept as constant as possible. To control the network we propose to use techniques of model predictive control (MPC) which solve a sequence of finite horizon optimal control problems. The method uses predictions of the state and computes a local optimal control which is then used for the model simulation in the next iteration. This framework is very flexible, since it allows control in real time and the incorporation of nonlinear dynamics and/or inequality constraints. It has already been used successfully within other areas of energy network control. Our new ansatz will also incorporate the stochastic effects into the model predictive control framework using data assimilation. Our vision is to develop numerical methods for network operators that allows the incorporation of model uncertainities for improving simulation and control of power networks.Website
X Close projectProf. Dr. Vladimir Spokoiny
Description
The project aims at developing numerical methods for the solution of complex optimal control problems arising in energy production, storage, and trading on energy markets. As a first step, we implement a Monte-Carlo approach to a hydro-electricity production and storage problem coupled with a stochastic model of the electricity market. Further we develop algorithms for pricing of complex energy derivatives based on the dual martingale approach.Website
X Close projectProf. Dr. Vladimir Spokoiny
Description
One of the main goals in this project is a systematic numerical treatment of generic optimal decision problems in real-life applications that encounter in energy markets by incorporating ``deep Learning'', a recent concept for data analysis and prediction. Also it is intended to include principles of deep learning in methods for forecasting and estimating price distribution processes in a systematic way.Website
X Close projectProf. Dr. Peter Karl Friz
Dr Mario Maurelli
Description
The aim of the project is to better understand and to give simulations for a successful model for the charging and discharging of lithium-ion batteries, which are currently the most promising storage devices to store and convert chemical energy into electrical energy and vice versa. The model exhibits phase transition under different small parameter regimes and gives rise to hysteresis. We study these phenomena using the interpretation of the model as a stochastic particle system, with the goal of providing stability bounds, fast simulations, improvement of the model itself and optimization of the device. More information...Website
X Close projectDescription
Many computational models are stochastic and the model output needs require some sort of sampling. Besides this intrinsic stochasticity, the models usually depend on a number of uncertain parameters. We develop a multi-level adaptive sparse grid strategy to address this parametric uncertainty, where the sampling effort is adjusted to the level of the sparse grid. This methodology is applied to stochastic simulation models of charge transport, as they appear in photovoltaics and photocatalysis.Website
X Close projectProf. Dr. Fredi Tröltzsch
Description
The project D-SE9 focuses on the analysis and efficient numerical solution of optimal control problems for nonlinear evolution equations with Maxwell's evolution equations as a challenging benchmark example. In particular, we aim at developing low rank approximation techniques for the solution of forward-backward optimality systems that arise whenever optimal control problems for evolution equations are considered. In this project, we thus merge existing expertise in optimal control and low rank matrix and tensor approximation.Website
X Close projectDescription
In the project D-SE10 we aspire to recover higher order tensors from a relatively small number of measurements using low rank assumptions. As straight forward generalizations of the matrix recovery techniques to the problem of tensor recovery are often either infeasible or impossible, the focus of this project is twofold. First, to investigate those generalizations that might still be feasible in a tensor setting in particular Riemannian methods on low rank tensor manifolds, and second, to apply and specialize existing techniques from tensor product approximation like the ALS to the tensor recovery and completion settings.Website
X Close projectDescription
Photocatalysis is a key application in the field of femtochemistry where chemical reaction dynamics is controlled by temporally shaped femtosecond laser pulses, with the target to promote specific product channels while suppressing competing undesired channels, e.g. pollutants. The optimal shaping of the laser pulse requires a detailed insight into the underlying reaction mechanisms at the atomic or molec- ular level that can often only be obtained by theoretical modelling and computer simulations of the quantum mechanical equations of motion. For catalytic system, this boils down to the iterated integration of the dissipative Liouville–von–Neumann (LvN) equation for reduced quantum mechanical density matrices, which represents the computational bottleneck for theoretical modelling, as the size of the matrices grows quadratically with the number of quantum states involved. The aim of this project is to study model order reduction (MOR) of LvN-based models to beat the curse of dimensionality in the simulation and (optimal) control of photocatalytic processes. In the setting of first-order perturbation theory, the laser field in these models is linearly coupled to the density matrix, which leads to a time- inhomogeneous bilinear system of equations of motion. MOR of bilinear systems has recently been a field of intense research. The downside of many available methods is their lack of structure preservation, most importantly, asymptotic stability of fixed points. An alternative that is in the focus of this project is MOR based on balancing the controllable and observable subspace of the system. Even though the identification of the essential subspace requires the solution of large-scale Lyapunov equations, which limits the applicability of the method to systems of moderate size (up to 100,000 DOFs), it has proven powerful for linear control systems in terms of computable error bounds and structure preservation. Whether these results carry over to the bilinear case is still open. Goals: Extending existing approaches to MOR from the linear to the bilinear case (required for LvN-based models) Developing, implementing and testing numerical methods for the solution of large-scale generalized Sylvester and Lyapunov equations Exploring structure preservation of MOR approaches Applying MOR approaches with optimal control of open quantum systems Identifying relevant photochemical benchmark systems to test various MOR / OC approachesWebsite
X Close projectDescription
The aim of this project is to develop fast linear solvers for heterogeneous saddle point problems as appearing during the iterative solution of non-smooth optimization problems, e.g., in the context of phase field models. While the development of nonlinear solvers for phase field models has reached a certain maturity, the existing solvers for the linear subproblems are still unsatisfactory. As a first step, we focus on two-phase models. Later, we plan to extend these solvers for multiphase systems.Matheon-C-SE12.php" target="_blank">Website
X Close projectPD Dr. René Henrion
Prof. Dr. Dietmar Hömberg
Prof. Dr. Reinhold Schneider
Dr. Thomas Petzold