Search Results
On the algorithmic solution of optimization problems subject to probabilistic/robust (probust) constraints
We present an adaptive grid refinement algorithm to solve probabilistic optimization problems with infinitely many random constraints. Using a bilevel approach, we iteratively aggregate inequalities that provide most information not in a geometric but in a probabilistic sense. This conceptual idea, for which a convergence proof is provided, is then adapted to an implementable algorithm. The efficiency of our approach when compared to naive methods based on uniform grid refinement is illustrated for a numerical test example as well as for a water reservoir problem with joint probabilistic filling level constraints.
Analysis and simulations for a phase-field fracture model at finite strains based on modified invariants
Phase-field models have already been proven to predict complex fracture patterns for brittle fracture at small strains. In this paper we discuss a model for phase-field fracture at finite deformations in more detail. Among the identification of crack location and projection of crack growth the numerical stability is one of the main challenges in solid mechanics. Here we present a phase-field model at finite strains, which takes into account the anisotropy of damage by applying an anisotropic split of the modified invariants of the right Cauchy-Green strain tensor. We introduce a suitable weak notion of solution that also allows for a spatial and temporal discretization of the model. In this framework we study the existence of solutions and we show that the time-discrete solutions converge in a weak sense to a solution of the time-continuous formulation of the model. Numerical examples in two and three space dimensions illustrate the range of validity of the analytical results.
Experimental proof of Joule heating-induced switched-back regions in OLEDs
Organic light-emitting diodes (OLEDs) have become a major pixel technology in the display sector, with products spanning the entire range of current panel sizes. The ability to freely scale the active area to large and random surfaces paired with flexible substrates provides additional application scenarios for OLEDs in the general lighting, automotive, and signage sectors. These applications require higher brightness and, thus, current density operation compared to the specifications needed for general displays. As extended transparent electrodes pose a significant ohmic resistance, OLEDs suffering from Joule self-heating exhibit spatial inhomogeneities in electrical potential, current density, and hence luminance. In this article, we provide experimental proof of the theoretical prediction that OLEDs will display regions of decreasing luminance with increasing driving current. With a two-dimensional OLED model, we can conclude that these regions are switched back locally in voltage as well as current due to insufficient lateral thermal coupling. Experimentally, we demonstrate this effect in lab-scale devices and derive that it becomes more severe with increasing pixel size, which implies its significance for large-area, high-brightness use cases of OLEDs. Equally, these non-linear switching effects cannot be ignored with respect to the long-term operation and stability of OLEDs; in particular, they might be important for the understanding of sudden-death scenarios. © 2020, The Author(s).
Modelling and simulation of flame cutting for steel plates with solid phases and melting
The goal of this work is to describe in detail a quasi-stationary state model which can be used to deeply understand the distribution of the heat in a steel plate and the changes in the solid phases of the steel and into liquid phase during the flame cutting process. We use a 3D-model similar to previous works from Thiébaud (J. Mater. Process. Technol. 214(2):304–310, 2014) and expand it to consider phases changes, in particular, austenite formation and melting of material. Experimental data is used to validate the model and study its capabilities. Parameters defining the shape of the volumetric heat source and the power density are calibrated to achieve good agreement with temperature measurements. Similarities and differences with other models from literature are discussed. © 2020, The Author(s).
Mode competition in broad-ridge-waveguide lasers
The lateral brightness achievable with high-power GaAs-based laser diodes having long and broad waveguides is commonly regarded to be limited by the onset of higher-order lateral modes. For the study of the lateral-mode competition two complementary simulation tools are applied, representing different classes of approximations. The first tool bases on a completely incoherent superposition of mode intensities and disregards longitudinal effects like spatial hole burning, whereas the second tool relies on a simplified carrier transport and current flow. Both tools yield agreeing power-current characteristics that fit the data measured for 5-23 µm wide ridges. Also, a similarly good qualitative conformance of the near and far fields is found. However, the threshold of individual modes, the partition of power between them at a given current, and details of the near and far fields show differences. These differences are the consequence of a high sensitivity of the mode competition to details of the models and of the device structure. Nevertheless, it can be concluded concordantly that the brightness rises with increasing ridge width irrespective of the onset of more and more lateral modes. The lateral brightness W mm-1at 10 MW cm-2 power density on the front facet of the investigated laser with widest ridge (23 µm) is comparable with best values known from much wider broad-area lasers. In addition, we show that one of the simulation tools is able to predict beam steering and coherent beam coupling without introducing any phenomenological coupling coefficient or asymmetries. © 2020 The Author(s). Published by IOP Publishing Ltd.
Quantitative Heat-Kernel Estimates for Diffusions with Distributional Drift
[For Abstract, see PDF]
On periodic solutions for one-phase and two-phase problems of the Navier–Stokes equations
This paper is devoted to proving the existence of time-periodic solutions of one-phase or two-phase problems for the Navier–Stokes equations with small periodic external forces when the reference domain is close to a ball. Since our problems are formulated in time-dependent unknown domains, the problems are reduced to quasilinear systems of parabolic equations with non-homogeneous boundary conditions or transmission conditions in fixed domains by using the so-called Hanzawa transform. We separate solutions into the stationary part and the oscillatory part. The linearized equations for the stationary part have eigen-value 0, which is avoided by changing the equations with the help of the necessary conditions for the existence of solutions to the original problems. To treat the oscillatory part, we establish the maximal Lp–Lq regularity theorem of the periodic solutions for the system of parabolic equations with non-homogeneous boundary conditions or transmission conditions, which is obtained by the systematic use of R-solvers developed in Shibata (Diff Int Eqns 27(3–4):313–368, 2014; On the R-bounded solution operators in the study of free boundary problem for the Navier–Stokes equations. In: Shibata Y, Suzuki Y (eds) Springer proceedings in mathematics & statistics, vol. 183, Mathematical Fluid Dynamics, Present and Future, Tokyo, Japan, November 2014, pp 203–285, 2016; Comm Pure Appl Anal 17(4): 1681–1721. https://doi.org/10.3934/cpaa.2018081, 2018; R boundedness, maximal regularity and free boundary problems for the Navier Stokes equations, Preprint 1905.12900v1 [math.AP] 30 May 2019) to the resolvent problem for the linearized equations and the transference theorem obtained in Eiter et al. (R-solvers and their application to periodic Lp estimates, Preprint in 2019) for the Lp boundedness of operator-valued Fourier multipliers. These approaches are the novelty of this paper. © 2020, The Author(s).
Convergence bounds for empirical nonlinear least-squares
We consider best approximation problems in a nonlinear subset ℳ of a Banach space of functions (𝒱,∥•∥). The norm is assumed to be a generalization of the L 2-norm for which only a weighted Monte Carlo estimate ∥•∥n can be computed. The objective is to obtain an approximation v ∈ ℳ of an unknown function u ∈ 𝒱 by minimizing the empirical norm ∥u − v∥n. We consider this problem for general nonlinear subsets and establish error bounds for the empirical best approximation error. Our results are based on a restricted isometry property (RIP) which holds in probability and is independent of the specified nonlinear least squares setting. Several model classes are examined and the analytical statements about the RIP are compared to existing sample complexity bounds from the literature. We find that for well-studied model classes our general bound is weaker but exhibits many of the same properties as these specialized bounds. Notably, we demonstrate the advantage of an optimal sampling density (as known for linear spaces) for sets of functions with sparse representations.
A discussion of the cell voltage during discharge of an intercalation electrode for various C-rates based on non-equilibrium thermodynamics and numerical simulations
In this work we discuss the modeling procedure and validation of a non-porous intercalation half-cell during galvanostatic discharge. The modeling is based on continuum thermodynamics with non-equilibrium processes in the active intercalation particle, the electrolyte, and the common interface where the intercalation reaction Li+ + e- ↔ Li occurs. The model is in detail investigated and discussed in terms of scalings of the non-equilibrium parameters, i.e. the diffusion coefficients DA and DE of the active phase and the electrolyte, conductivity sA and sE of both phases, and the exchange current density e0L, with numerical solutions of the underlying PDE system. The current density i as well as all non-equilibrium parameters are scaled s with respect to the 1-C current density iC A of the intercalation electrode. We compute then numerically the cell voltage E as function of the capacity Q and the C-rate Ch. Within a hierarchy of approximations we provide computations of E(Q) for various scalings of the diffusion coefficients, the conductivities and the exchange current density. For the later we provide finally a discussion for possible concentration dependencies. © The Author(s) 2019. Published by ECS.
Beyond just “flattening the curve”: Optimal control of epidemics with purely non-pharmaceutical interventions
When effective medical treatment and vaccination are not available, non-pharmaceutical interventions such as social distancing, home quarantine and far-reaching shutdown of public life are the only available strategies to prevent the spread of epidemics. Based on an extended SEIR (susceptible-exposed-infectious-recovered) model and continuous-time optimal control theory, we compute the optimal non-pharmaceutical intervention strategy for the case that a vaccine is never found and complete containment (eradication of the epidemic) is impossible. In this case, the optimal control must meet competing requirements: First, the minimization of disease-related deaths, and, second, the establishment of a sufficient degree of natural immunity at the end of the measures, in order to exclude a second wave. Moreover, the socio-economic costs of the intervention shall be kept at a minimum. The numerically computed optimal control strategy is a single-intervention scenario that goes beyond heuristically motivated interventions and simple “flattening of the curve”. Careful analysis of the computed control strategy reveals, however, that the obtained solution is in fact a tightrope walk close to the stability boundary of the system, where socio-economic costs and the risk of a new outbreak must be constantly balanced against one another. The model system is calibrated to reproduce the initial exponential growth phase of the COVID-19 pandemic in Germany. © 2020, The Author(s).