1 / 64

# Introduction to Model Order Reduction - PowerPoint PPT Presentation

Introduction to Model Order Reduction. Luca Daniel Massachusetts Institute of Technology [email protected] http://onigo.mit.edu/~dluca/2006PisaMOR. www.rle.mit.edu/cpg. Course Outline. Numerical Simulation Quick intro to PDE Solvers Quick intro to ODE Solvers Model Order reduction

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.

## PowerPoint Slideshow about ' Introduction to Model Order Reduction' - sawyer-herrera

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

### Introduction to Model Order Reduction

Luca Daniel

Massachusetts Institute of Technology

[email protected]

http://onigo.mit.edu/~dluca/2006PisaMOR

www.rle.mit.edu/cpg

Numerical Simulation

Quick intro to PDE Solvers

Quick intro to ODE Solvers

Model Order reduction

Linear systems

Common engineering practice

Optimal techniques in terms of model accuracy

Efficient techniques in terms of time and memory

Non-Linear Systems

Parameterized Model Order Reduction

Linear Systems

Non-Linear Systems

Today

• Introduction

• Reduction of weakly nonlinear systems (Volterra Series)

• Reduction of strongly nonlinear systems Trajectory PieceWise Linear (TPWL) and Polynomial (PWP)

• with moment matching

• with Truncated Balance Realization

• Inlet density disturbance d:

DISCRETIZE

NonLinear Reduction

PDE Field Solvers or Circuit Simulators

Non-Linear analog components

e.g. MEMs, VCO, LNA

Projection Framework

Equation Testing

for nonlinear dynamic systems

in

• Substitute:

• Reduced system:

• A problem:

q=10

n=104

n=104

q=10

• Using VTF(Ux) is too expensive!

• In what basis should we project?

• No simple frequency domain insight

• No eigenmodes

• No Krylov subspace

• How do you represent the projection?

• New problem for nonlinear

• Analysis of linearized models [Ma88]

• Sampling of time-simulation data + principal components analysis [Sirovich87]

• Nonlinear balancing [Scherpen93],

• Time-derivatives [Gunupudi99], …

for NonLinear Systems

• Representation of F(x) using a polynomial (e.g. Taylor’s expansions, Volterra Series) [Phillips00]*****

• Representation of F(x) using several linearizations (Trajectory Piece-Wise Linear TPWL) [Rewienski01]

• Representation of F(x) with several polynomials (PWP PieceWise Polynomial) [Dong03]

III.2 – Reduction of Weakly Non Linear Systems

Luca Daniel

Thanks to Joel Phillips, Cadence Berkeley Labs

• Approximate nonreduced model:

MOR for nonlinear dynamic systems

(J. Phillips, Y. Chen, F. Wang):

• Linear, quadratic reduced order models:

Differential equation becomes

To match first few terms in functional series expansion, only need first few polynomial terms

Each term is a -dimensional

tensor, represented as an matrix

Polynomial Approximations

[note : not static models; frequency effects retained]

where

etc.

Identity for Kronecker products

Project tensors

Gives reduced model

Projection of polynomial terms

etc.

Projection procedure produces reduced model in same polynomial form

Key tensor components are compressed to lower dimensionality

Kronecker forms provide convenient general notation

How to get projection spaces V?

Reduced polynomial models

Introduce variational parameter polynomial form

Consider parametrized system

Expand response in powers of

Substitute into differential equation

Compare and collect terms for each power of

Variational Analysis

Coupled set of differential equations polynomial form

Each set computes response at one order

Each set is linear in its own state variable

First set is just linearized model

Variational Analysis

etc.

Krylov subspace for projection

Second order equation set driven by input

Key insight : recall is well-approximated by thus

Model Reduction

Generate second Krylov space using

Compute models by projection

each equation order set separately, or

project original polynomial system using UNION of

number of matched moments of order-i given by size of appropriate Krylov subspace

Model Reduction

potentially large models

intrinsic in polynomial descriptions (e.g. Volterra series)

Practical for weak nonlinearities needing only few terms in functional series

example: RF datapath

strong nonlinearities included in bias conditions

Many local optimizations, extensions are possible

two-sided projectors

multiple reduction steps

Computational Considerations

-20 order

-60

-80

0

50

100

150

200

MHz

RF mixer example

dB

Projection + polynomial series + variational analysis = rigorous nonlinear reduction procedure

Compact Kronecker notation is very useful

Is there a way to reduce number of terms in higher-order Krylov space?

Smaller models faster

Diagnose when models based on first order system are failing?

Follow-on reduction?

Summary weakly nonlinear MOR

Model Order Reduciton rigorous nonlinear reduction procedure

for NonLinear Systems

• Representation of F(x) using a polynomial (e.g. Taylor’s expansions, Volterra Series) [Phillips00]*****

• Representation of F(x) using several linearizations (Trajectory Piece-Wise Linear TPWL) [Rewienski01]

• Representation of F(x) with several polynomials (PWP PieceWise Polynomial) [Dong03]

### III.3 - Reduction of Strongly NonLinear Systems via Trajectory PieceWise Linear Method

Luca Daniel

Massachusetts Institute of Technology

Thanks to Michał Rewieński, Jacob White and Brad Bond

www.rle.mit.edu/cpg

Piecewise-linear representation Trajectory PieceWise Linear Method

• Linearizations around xi , i=0,…, s-1

• Weighted combination of the models:

• Project linearized models:

Background – TPWL [Reiwenski01] Trajectory PieceWise Linear Method

x8

x1

x7

x6

x5

x2

x4

x14

x10

xa

xb

x0

x11

x9

x13

x3

x15

Model i only valid near xi

x2

#linearizations =#samplesn

n = 104

#samples = 100

10010000 = LARGE

Use collection of linear models

x1

Background – TPWL: Picking Linearization Points Trajectory PieceWise Linear Method

Use training trajectories to pick linearization points

y(t)

x2

State Space

Time Domain Simulation

t

x1

Linearization at current state xi

Quasi-piecewise-linear MOR Trajectory PieceWise Linear Method– computing weights

• For i=0,…,(s-1) compute:

• Compute

• For i=0,…,(s-1) compute:

• Normalize wi.

= Trajectory PieceWise Linear Method

Model from linearization 1

Model from linearization k

Model from linearization 2

K2

A1

K1

A2

Ak

Kk

Background – TPWLReduction of the Linearized Systems

= Trajectory PieceWise Linear Method

Background – TPWLConstructing the projection matrix V

Use moments from EACH linear model to construct V

Background – TPWL: Weighting / Simulation Trajectory PieceWise Linear Method[Riewinski01,Tiwary05,Dong05]

Linearization 3

x2

Use weighting functions to combine linear models during simulation

also well

approximated

C – poorly

approximated

Well approximated

Current state

Linearization 2

x1

Linearization 1

Test example of a nonlinear circuit: Trajectory PieceWise Linear Method

Computational results – circuit example Trajectory PieceWise Linear Method

• Input:

training

input

testing

input

Computational results - circuit example, SSS Trajectory PieceWise Linear Method

DISCRETIZE Trajectory PieceWise Linear Method

Example of a microswitch

Computational results - microswitch example Trajectory PieceWise Linear Method

• Input:

training

input

testing

input

Computational results - microswitch example Trajectory PieceWise Linear Method

Computational results – microswitch example Trajectory PieceWise Linear Method

training

input

testing

input

Simulation with the piecewise-linear model Trajectory PieceWise Linear Method

Generate the Trajectory PieceWise Linear Methodreduced linear model

Use the reduced linear system

for simulation

x2

x1

x0

A

x3

Fast approximate simulation

x2

x1

Computational results - model order reduction Trajectory PieceWise Linear Method

*Matlab implementation

Computational results – OpAmp example Trajectory PieceWise Linear Method

OpAmp:

• 70 MOSFETs

• 13 resistors

• 9 capacitors

• 51 nodes

training input

testing

input

Computational results – OpAmp example Trajectory PieceWise Linear Method

OpAmp:

• 70 MOSFETs

• 13 resistors

• 9 capacitors

• 51 nodes

Computational results – OpAmp example, SSS Trajectory PieceWise Linear Method

• Density disturbance d:

Computational results inlet– shock propagation in 1D

• Burgers equation:

• How should we pick training inputs?

• How “often” should we compute linearized models along the training trajectory?

• What are the errors of representing a nonlinear systems as a trajectory piecewise linear system?

• What are “the best” reduced order bases for the trajectory PWL model?

Conclusions inlet

• TPWL macromodels can be accurate for a number of strongly nonlinear dynamical systems, outperforming reduced models based on single-state polynomial expansions.

• TPWL models use simple, cost-efficient representation of system’s nonlinearity, resulting in very short simulation times.

• TPWL models may be extracted very efficiently by using a fast approximate simulator.

• Introduction

• Reduction of weakly nonlinear systems (Volterra Series)

• Reduction of strongly nonlinear systems (TPWL) Trajectory Pieace Wise Linear and Polynomial

• with moment matching

• with Truncated Balance Realization

### III.4 Reduction of Strongly NonLinear Systems inletvia Trajectory-PieceWise Linear (TPWL) + Truncated Balance Realization (TBR)

Luca Daniel

Massachusetts Institute of Technology

Thanks to Dmitry Vasilyev, Michał Rewieński, Jacob White

Fast

Don’t guarantee accuracy

Balanced-truncation methods

Expensive (~n3)

Guarantee accuracy

Obtaining projection basis

For example,

V=W=

colspan(A-1B, A-2B, … , A-q B)

Let’s now use this!

Our Approach: inlet

• Use TPWL to handle nonlinearity

• Before we used Krylov-subspace linear reduction (less accurate)

• Here we use TBR for projection matrices W and V

x0

x2

x1

xn

Hankel operator inlet

TBR reduction

u

y

LTI SYSTEM

t

t

Past

input

Future

output

X (state)

P (controllability)

Which states are easier to reach?

Q (observability)

Which states produces more output?

TBR algorithm includes into projection basis

most controllable and most observable states

Numerical results – inletRLC transmission line

TBR-based TPWL beat

Krylov-based

4-th order TBR TPWL reaches the limit of TPWL representation

Error in transient

||yr – y||2

Order of the reduced model

Krylov reduction (Wi , Vi):

Ai = WiTAVi

Bi = WiTB

Ci = CVi

Initial Model:

(A B C), n

Intermediate Model:

(Ai Bi Ci), ni

TBR reduction (Wt , Vt):

Ar = WtTAVt

Br = WtTB

Cr = CVt

Reduced Model:

(Ar Br Cr), q

Micromachined device example algorithm

FD model

non-symmetric

indefinite Jacobian

Cost of Krylov-TBR almost equals Krylov

*Matlab implementation

Comparing accuracies of Krylov TPWL method and TBR-based TPWL algorithm

5x reduction in order – 125x improvement in efficiency

*Testing input equal to training input

Future work TPWL algorithm

• Are TBR-based TPWL models valid for unstable linearizations?

• What about systems having the following form (i.e. circuits with nonlinear capacitors):

Conclusions TPWL algorithm

• TBR-based linear reduction procedure to generate TPWL reduced models

• Order reduced 5 times while maintaining comparable accuracy with Krylov TPWL method (efficiency improved 125 times!)

• Combined Krylov-TBR reduction allows to extract TPWL models at low cost

Course Outline TPWL algorithm

Numerical Simulation

Quick intro to PDE Solvers

Quick intro to ODE Solvers

Model Order reduction

Linear systems

Common engineering practice

Optimal techniques in terms of model accuracy

Efficient techniques in terms of time and memory

Non-Linear Systems

Parameterized Model Order Reduction

Linear Systems

Non-Linear Systems

Today

Tomorrow