Model Order Reduction of Non-Linear Magnetostatic Problems Based on POD and DEI Methods

In the domain of numerical computation, Model Order Reduction approaches are more and more frequently applied in mechanics and have shown their efficiency in terms of reduction of computation time and memory storage requirements. One of these approaches, the Proper Orthogonal Decomposition (POD), can be very efficient in solving linear problems but encounters limitations in the non-linear case. In this paper, the Discret Empirical Interpolation Method coupled with the POD method is presented. This is an interesting alternative to reduce large-scale systems deriving from the discretization of non-linear magnetostatic problems coupled with an external electrical circuit

Model order reduction of large scale ODE systems : MOR for ANSYS versus ROM workbench

In this paper we compare the numerical results obtained by different model order reduction software tools, in order to test their scalability for relevant problems of the microelectronic-industry. MOR for ANSYS [2J is implemented in C++ and ROM Workbench [3J is a MATLAB code. The chosen benchmarks are large scale linear ODE systems, which arise from the finite element discretisation of electro-thermal MEMS models

Construction of high-order multirate Rosenbrock methods for stiff ODEs

Multirate time stepping is a numerical technique for efficiently solving large-scale ordinary differential equations (ODEs) with widely different time scales localized over the components. This technique enables one to use large time steps for slowly varying components, and small steps for rapidly varying ones. Multirate methods found in the literature are normally of low order, one or two. Focusing on stiff ODEs, in this paper we discuss multirate methods based on the higher-order, stiff Rosenbrock integrators. Special attention is paid to the treatment of the refinement interfaces with regard to the choice of the interpolant and the occurrence of order reduction. For stiff, linear systems containing a stiff source term, we propose modifications for the treatment of the source term which overcome order reduction originating from such terms and which we can implement in our multirate method

Trajectory Tracking Control Design for Large-Scale Linear Dynamical Systems With Applications to Soft Robotics

International audienceThis article presents new results to control process modeled through linear large-scale systems. Numerical methods are widely used to model physical systems, and the finite-element method is one of the most common methods. However, for this method to be precise, it requires a precise spatial mesh of the process. Large-scale dynamical systems arise from this spatial discretization. We propose a methodology to design an observer-based output feedback controller. First, a model reduction step is used to get a system of acceptable dimension. Based on this low-order system, two linear matrix inequality problems provide us, respectively, with the observer and controller gains. In both the cases, model and reduction errors are taken into account in the computations. This provides robustness with respect to the reduction step and guarantees the stability of the original large-scale system. Finally, the proposed method is applied to a physical setup-a soft robotics platform-to show its feasibility

Persistent Model Order Reduction for Complex Dynamical Systems Using Smooth Orthogonal Decomposition

Full-scale complex dynamic models are not effective for parametric studies due to the inherent constraints on available computational power and storage resources. A persistent reduced order model (ROM) that is robust, stable, and provides high-fidelity simulations for a relatively wide range of parameters and operating conditions can provide a solution to this problem. The fidelity of a new framework for persistent model order reduction of large and complex dynamical systems is investigated. The framework is validated using several numerical examples including a large linear system and two complex nonlinear systems with material and geometrical nonlinearities. While the framework is used for identifying the robust subspaces obtained from both proper and smooth orthogonal decompositions (POD and SOD, respectively), the results show that SOD outperforms POD in terms of stability, accuracy, and robustness

Efficient Numerical Solution of Large Scale Algebraic Matrix Equations in PDE Control and Model Order Reduction

Matrix Lyapunov and Riccati equations are an important tool in mathematical systems theory. They are the key ingredients in balancing based model order reduction techniques and linear quadratic regulator problems. For small and moderately sized problems these equations are solved by techniques with at least cubic complexity which prohibits their usage in large scale applications. Around the year 2000 solvers for large scale problems have been introduced. The basic idea there is to compute a low rank decomposition of the quadratic and dense solution matrix and in turn reduce the memory and computational complexity of the algorithms. In this thesis efficiency enhancing techniques for the low rank alternating directions implicit iteration based solution of large scale matrix equations are introduced and discussed. Also the applicability in the context of real world systems is demonstrated. The thesis is structured in seven central chapters. After the introduction chapter 2 introduces the basic concepts and notations needed as fundamental tools for the remainder of the thesis. The next chapter then introduces a collection of test examples spanning from easily scalable academic test systems to badly conditioned technical applications which are used to demonstrate the features of the solvers. Chapter four and five describe the basic solvers and the modifications taken to make them applicable to an even larger class of problems. The following two chapters treat the application of the solvers in the context of model order reduction and linear quadratic optimal control of PDEs. The final chapter then presents the extensive numerical testing undertaken with the solvers proposed in the prior chapters. Some conclusions and an appendix complete the thesis