Some contributions to the model order reduction of large scale non-linear electric circuits

Model order Reduction (MOR) has become an ubiquitous technique in the simulation of large-scale dynamical systems (i.e. 10⁴ and more equations). One technique for non-linear MOR is the trajectory piecewise-linear approach (TPWL) [1]. TPWL approximates a non-linear differential system by a weighted sum of linear systems which have a significantly reduced number of equations. One open question is which weight representations provide physical meaning of the weighted sum [2]. In this article we propose two representations. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

No Borel Connections for the Unsplitting Relations

We prove that there is no Borel connection for non‐trivial pairs of unsplitting relations. This was conjectured in [3].     

On the numerical solution of large-scale sparse discrete-time Riccati equations

We discuss the numerical solution of large-scale discrete-time algebraic Riccati equations (DAREs) as they arise, e.g., in fully discretized linear-quadratic optimal control problems for parabolic partial differential equations (PDEs). We employ variants of Newton??s method that allow to compute an approximate low-rank factor of the solution of the DARE. The principal computation in the Newton iteration is the numerical solution of a Stein (aka discrete Lyapunov) equation in each step. For this purpose, we present a low-rank Smith method as well as a low-rank alternating-direction-implicit (ADI) iteration to compute low-rank approximations to solutions of Stein equations arising in this context. Numerical results are given to verify the efficiency and accuracy of the proposed algorithms.     

On the influence of design parameters and nonlinearities in vehicle suspensions on road deformation

Among others, two main objectives of modern vehicle design are road friendliness and ride comfort. Both aspects are strongly related since the dynamical tire forces depend on the vertical acceleration of the vehicle. In order to investigate the influence of design and operation parameters, different car models are considered which move with constant velocity on a rippled road. First, a linear half car model is examined and the influence of different design parameters is discussed. Second, nonlinear suspensions with Coulomb friction due to sealings as well as with bilinear shock absorbers are taken into account. The vertical dynamics of the vehicle model and the dynamic tire forces between vehicle and road are calculated using analytical methods. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Missing Point Estimation for steady aerodynamic applications

Model reduction in aerodynamic applications usually makes use of Proper Orthogonal Decomposition (POD). In this work a POD-based method, called Missing Point Estimation (MPE), will be applied to steady state flows with variation of the flow parameter angle of attack. The basic idea of MPE is to select a subset of the computational grid points (control volumes) and limit the governing equations to these. Subsequently, the remaining equations are projected onto the POD subspace. This approach has the advantage that the nonlinear right hand side of the governing equations has to be evaluated only for some selected points (control volumes) which makes the MPE faster than classical Galerkin projection. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Meeting infinitely many cells of a partition once

We investigate several versions of a cardinal characteristic defined by Frankiewicz. Vojtáš showed , and Blass showed . We show that all the versions coincide and that is greater than or equal to the splitting number. We prove the consistency of and of . Received: 2 October 1996 / Revised version: 22 May 1997     

On numerical methods for discrete least-squares approximation by trigonometric polynomials

Fast, efficient and reliable algorithms for discrete least-squares approximation of a real-valued function given at arbitrary distinct nodes in by trigonometric polynomials are presented. The algorithms are based on schemes for the solution of inverse unitary eigenproblems and require only arithmetic operations as compared to operations needed for algorithms that ignore the structure of the problem. An algorithm which solves this problem with real-valued data and real-valued solution using only real arithmetic is given. Numerical examples are presented that show that the proposed algorithms produce consistently accurate results that are often better than those obtained by general QR decomposition methods for the least-squares problem.

     

Measuring Lineup Difficulty By Matching Distance Metrics With Subject Choices in Crowd-Sourced Data

Graphics play a crucial role in statistical analysis and data mining. Being able to quantify structure in data that is visible in plots, and how people read the structure from plots is an ongoing challenge. The lineup protocol provides a formal framework for data plots, making inference possible. The data plot is treated like a test statistic, and lineup protocol acts like a comparison with the sampling distribution of the nulls. This article describes metrics for describing structure in data plots and evaluates them in relation to the choices that human readers made during several large Amazon Turk studies using lineups. The metrics that were more specific to the plot types tended to better match subject choices, than generic metrics. The process that we followed to evaluate metrics will be useful for general development of numerically measuring structure in plots, and also in future experiments on lineups for choosing blocks of pictures. Supplementary materials for this article are available online.     

Framed! Reproducing and Revisiting 150-Year-Old Charts

The Statistical Atlases published by the Census Bureau in the late 1800s utilized a number of novel methods for displaying data. In this paper, we examine the use of framed spine and mosaic plots used in two plates of the Statistical Atlas of 1870. We use forensic statistics to recreate the data using available census information, and then use that data to create framed charts using modern plotting methods. We then examine the effectiveness of the framed charts compared to other alternatives with a user study. The data and code for this study are available online.     

Phase Field Modeling of Brittle and Ductile Fracture

The phase field modeling of brittle fracture was a topic of intense research in the last few years and is now well-established. We refer to the work [1-3], where a thermodynamically consistent framework was developed. The main advantage is that the phase-field-type diffusive crack approach is a smooth continuum formulation which avoids the modeling of discontinuities and can be implemented in a straightforward manner by multi-field finite element methods. Therefore complex crack patterns including branching can be resolved easily. In this paper, we extend the recently outlined phase field model of brittle crack propagation [1-3] towards the analysis of ductile fracture in elastic-plastic solids. In particular, we propose a formulation that is able to predict the brittle-to-ductile failure mode transition under dynamic loading that was first observed in experiments by Kalthoff and Winkler [4]. To this end, we outline a new thermodynamically consistent framework for phase field models of crack propagation in ductile elastic-plastic solids under dynamic loading, develop an incremental variational principle and consider its robust numerical implementation by a multi-field finite element method. The performance of the proposed phase field formulation of fracture is demonstrated by means of the numerical simulation of the classical Kalthoff-Winkler experiment that shows the dynamic failure mode transition. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)