A Hamiltonian Krylov-Schur-type method based on the symplectic Lanczos process

We discuss a Krylov-Schur like restarting technique applied within the symplectic Lanczos algorithm for the Hamiltonian eigenvalue problem. This allows to easily implement a purging and locking strategy in order to improve the convergence properties of the symplectic Lanczos algorithm. The Krylov-Schur-like restarting is based on the SR algorithm. Some ingredients of the latter need to be adapted to the structure of the symplectic Lanczos recursion. We demonstrate the efficiency of the new method for several Hamiltonian eigenproblems

Time integration and steady-state continuation for 2d lubrication equations

Lubrication equations allow to describe many structurin processes of thin liquid films. We develop and apply numerical tools suitable for their analysis employing a dynamical systems approach. In particular, we present a time integration algorithm based on exponential propagation and an algorithm for steady-state continuation. In both algorithms a Cayley transform is employed to overcome numerical problems resulting from scale separation in space and time. An adaptive time-step allows to study the dynamics close to hetero- or homoclinic connections. The developed framework is employed on the one hand to analyse different phases of the dewetting of a liquid film on a horizontal homogeneous substrate. On the other hand, we consider the depinning of drops pinned by a wettability defect. Time-stepping and path-following are used in both cases to analyse steady-state solutions and their bifurcations as well as dynamic processes on short and long time-scales. Both examples are treated for two- and three-dimensional physical settings and prove that the developed algorithms are reliable and efficient for 1d and 2d lubrication equations, respectively.Comment: 33 pages, 16 figure

KRYLOV SUBSPACE METHODS FOR SOLVING LARGE LYAPUNOV EQUATIONS

Published versio

Residual, restarting and Richardson iteration for the matrix exponential

A well-known problem in computing some matrix functions iteratively is a lack of a clear, commonly accepted residual notion. An important matrix function for which this is the case is the matrix exponential. Assume, the matrix exponential of a given matrix times a given vector has to be computed. We interpret the sought after vector as a value of a vector function satisfying the linear system of ordinary differential equations (ODE), whose coefficients form the given matrix. The residual is then defined with respect to the initial-value problem for this ODE system. The residual introduced in this way can be seen as a backward error. We show how the residual can efficiently be computed within several iterative methods for the matrix exponential. This completely resolves the question of reliable stopping criteria for these methods. Furthermore, we show that the residual concept can be used to construct new residual-based iterative methods. In particular, a variant of the Richardson method for the new residual appears to provide an efficient way to restart Krylov subspace methods for evaluating the matrix exponential.\u

Matrix-interpolation-based parametric model order reduction for multiconductor transmission lines with delays

A novel parametric model order reduction technique based on matrix interpolation for multiconductor transmission lines (MTLs) with delays having design parameter variations is proposed in this brief. Matrix interpolation overcomes the oversize problem caused by input-output system-level interpolation-based parametric macromodels. The reduced state-space matrices are obtained using a higher-order Krylov subspace-based model order reduction technique, which is more efficient in comparison to the Gramian-based parametric modeling in which the projection matrix is computed using a Cholesky factorization. The design space is divided into cells, and then the Krylov subspaces computed for each cell are merged and then truncated using an adaptive truncation algorithm with respect to their singular values to obtain a compact common projection matrix. The resulting reduced-order state-space matrices and the delays are interpolated using positive interpolation schemes, making it computationally cheap and accurate for repeated system evaluations under different design parameter settings. The proposed technique is successfully applied to RLC (R-resistor, L-inductor, C-capacitance) and MTL circuits with delays

Order reduction methods for solving large-scale differential matrix Riccati equations

We consider the numerical solution of large-scale symmetric differential matrix Riccati equations. Under certain hypotheses on the data, reduced order methods have recently arisen as a promising class of solution strategies, by forming low-rank approximations to the sought after solution at selected timesteps. We show that great computational and memory savings are obtained by a reduction process onto rational Krylov subspaces, as opposed to current approaches. By specifically addressing the solution of the reduced differential equation and reliable stopping criteria, we are able to obtain accurate final approximations at low memory and computational requirements. This is obtained by employing a two-phase strategy that separately enhances the accuracy of the algebraic approximation and the time integration. The new method allows us to numerically solve much larger problems than in the current literature. Numerical experiments on benchmark problems illustrate the effectiveness of the procedure with respect to existing solvers

Pressure Bifurcation Phenomenon on Supersonic Blowing Trailing Edges

Turbine blades operating in transonic-supersonic regime develop a complex shock wave system at the trailing edge, a phenomenon that leads to unfavorable pressure perturbations downstream and can interact with other turbine stages. Understanding the fluid behavior of the area adjacent to the trailing edge is essential in order to determine the parameters that have influence on these pressure fluctuations. Colder flow, bled from the high-pressure compressor, is often purged at the trailing edge to cool the thin blade edges, affecting the flow behavior and modulating the intensity and angle of the shock waves system. However, this purge flow can sometimes generate non-symmetrical configurations due to a pressure difference that is provoked by the injected flow. In this work, a combination of RANS simulations and global stability analysis is employed to explain the physical reasons of this flow bifurcation. Analyzing the features that naturally appear in the flow and become dominant for some value of the parameters involved in the problem, an anti-symmetrical global mode, related to the sudden geometrical expansion of the trailing edge slot, is identified as the main mechanism that forces the changes in the flow topology.Comment: Submitted to AIAA Journa