Top Banner
Modeling Feedback Control of Thermal Systems in COMSOL © Copyright 2015 COMSOL. COMSOL, COMSOL Multiphysics, Capture the Concept, COMSOL Desktop, COMSOL Server, and LiveLink are either registered trademarks or trademarks of COMSOL AB. All other trademarks are the property of their respective owners, and COMSOL AB and its subsidiaries and products are not affiliated with, endorsed by, sponsored by, or supported by those trademark owners. For a list of such trademark owners, see www.comsol.com/trademarks Lexi Carver Technical Marketing Engineer COMSOL Jon Ebert, PhD Director SC SOLUTIONS
43

0806ThermalSystems_webinar_2015

Apr 14, 2016

Download

Documents

ESTE DOCUMENTO ES INTRODUCTORIA AL ANALISIS DE SISTEMAS TERMICOS POR ELEMENTOS FINITOS
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: 0806ThermalSystems_webinar_2015

Modeling Feedback Control of Thermal Systems in COMSOL

© Copyright 2015 COMSOL. COMSOL, COMSOL Multiphysics, Capture the Concept, COMSOL Desktop, COMSOL Server, and LiveLink are either registered trademarks or trademarks of COMSOL AB. All other trademarks are the property of their respective owners, and COMSOL AB and its subsidiaries and products are not affiliated with, endorsed by, sponsored by, or supported by those trademark owners. For a list of such trademark owners, see www.comsol.com/trademarks

Lexi Carver Technical Marketing

Engineer COMSOL

Jon Ebert, PhD Director

SC SOLUTIONS

Page 2: 0806ThermalSystems_webinar_2015

Agenda• Introduction to COMSOL

Multiphysics® software• Simulating Feedback Control

of Thermal Systems• Demo: Rapid Thermal

Processing (RTP) System• Q&A• How To– Try COMSOL Multiphysics®– Contact Us

Prestressed plated metal layers in a model of a micromirror

Page 3: 0806ThermalSystems_webinar_2015

Why Simulate?• Conception and understanding

– Enables innovation

• Design and optimization– Achieve the highest possible

performance

• Testing and verification– Virtual testing is much faster

than testing physical prototypes

Displacement induced by thermal strains and applied pressure in a model of a capacitive pressure sensor

Page 4: 0806ThermalSystems_webinar_2015

Simulating with COMSOL Multiphysics®

• Electrical, mechanical, fluid, and chemical simulations

• Multiphysics – include and couple all relevant physical effects

• Single physics in one integrated environment

• Cross-disciplinary product development

Page 5: 0806ThermalSystems_webinar_2015

Model BuilderProvides instant access to any of the model settings• CAD/Geometry•Materials• Physics•Mesh• Solvers• Results

A Complete Simulation Environment

Graphics WindowUltrafast graphic presentation,stunning visualization

COMSOL Desktop®Straightforward to use, the Desktop gives insight and full control over the modeling process

Page 6: 0806ThermalSystems_webinar_2015

Product Suite – COMSOL 5.1

Page 7: 0806ThermalSystems_webinar_2015

Application Design Tools

Simulation ApplicationAny COMSOL model can be turned into an app with its own interface using the tools provided in the Application Builder

Application BuilderProvides all the tools needed to build and run simulation apps• Form Editor•Method Editor

Page 8: 0806ThermalSystems_webinar_2015

Run Applications

Simulation AppsThey can be run in a COMSOL® Client for Windows® and major web browsers

COMSOL Server™It’s the engine for running COMSOL apps and the hub for controlling their deployment, distribution, and use

Microsoft and Windows are either registered trademarks or trademarks of Microsoft Corporation in the United States and/or other countries.

Page 9: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 9August 2015

Feedback Control of Thermal Systemsin COMSOL

Jon EbertSC Solutions, Inc.

Email: [email protected]

Page 10: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 10August 2015

Overview

Feedback control of temperature is important in many manufacturing processes.

At SC we have developed methods for the model-based control of systems where we use physics-based models to develop high-performance temperature control.

This presentation will give an example of one method for designing a simple PID feedback controller (Q-Design).

We will present a simple example of feedback control within a COMSOL model.

A demo with 5-zone temperature control of a Rapid Thermal Processing (RTP) system will be shown.

Page 11: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 11August 2015

Closed-Loop Transfer Function

The Closed-Loop System

PC

Closed-Loop Response:

For a scalar system

Page 12: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 12August 2015

Proportional + Integral + Derivative Control (PID)

This is the most common type of feedback controller

The “gains”, Kp, Ki, and Kd (and td) must be selected by the engineer

There are systematic methods

Here we’ll outline ‘Q’-design

To limit high-frequency noise amplification, we usually low-pass filter the error driving the derivative control:

Page 13: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 13August 2015

Model-Based Control (Q-Design)

Incorporate a mathematical model of the system directly into the controller.

Often referred to as Q-parameterization or Youla parameterization.

References for Q-parameterization Control Design

For stable P, ALL stable controllers can be expressed in this form!

Control design becomes choice of transfer function Q

Page 14: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 14August 2015

A Model-Based Design (‘Q-Design’)

Closed-Loop Response:

If we select T and we know P, we can calculate C

For scalar (single input, single output)

Page 15: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 15August 2015

Selecting T (Closed-Loop Transfer Function)

Or, in time domain

Or, in time domain

Page 16: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 16August 2015

Solving for Model-Based Controller Gains

where and select

Resulting controller:

Our PID form:

Tuning PID

It turns out that if you pick T as we have done for the second order P, the controller is in the form of a PID controller.

Assume the system (P) is a second order system

Page 17: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 17August 2015

PID Example (w=1, z=1)

Page 18: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 18August 2015

PID Example (w=2, z=1)

Page 19: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 19August 2015

Implementing PID in COMSOL (An Example)

One dimensional heat transfer in a plate

Governing Equations:

Measure y on this surface

Heat surface at x=0 and measure temperature at surface x=L (y=T2)

Page 20: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 20August 2015

Build the Model in COMSOL

Page 21: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 21August 2015

Model Parameters

Page 22: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 22August 2015

Apply Boundary Conditions

Control input

Page 23: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 23August 2015

Create Mesh and Run the Simulation

This is not a second order system, but a much higher order PDE.

When we discretize the PDE using finite elements, we create a system of ODE’s.

For the mesh used here, there are 101 differential equations, or degrees of freedom (DOF).

Still, we will look for a PID controller based on a second order system.

Page 24: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 24August 2015

Open-Loop Response

We step the input flux from 0 to qmax*uctrl

Page 25: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 25August 2015

Measure Model Response Properties from OL Step

Time where y crosses 63% of final value

Page 26: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 26August 2015

Finding Approximate Time Constants

Measure the rate of change in y.Maximum is approximately 18[K/s]

Page 27: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 27August 2015

Approximate 2nd Order Model Parameters

We can use these to compute controller parameters.

We also need to pick the closed-loop speed, w, and damping, z.

Approximate 2nd order model parameters for plate model:

Page 28: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 28August 2015

Implementing PID in COMSOL

Proportional control:

Integral control:

Derivative control:

We solve these differential equations in COMSOL.

Page 29: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 29August 2015

Implementing PID in COMSOLAdd “Global ODEs and DAEs” to Model

Enter the Equations

State variableEquation

Page 30: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 30August 2015

Implementing PID in COMSOLDefine T2 as the plate temperature at x=L

Define control variables

The control: limit to 0<= uctrl <=1

Page 31: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 31August 2015

Define the Reference

Define the reference using an interpolation table.

Page 32: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 32August 2015

Controller Parameters

Tuning PID

Page 33: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 33August 2015

Closed-Loop Response

Works pretty well.

The open-loop response took almost 100s to reach the setpoint, the closed loop gets there in 12s.

Sensor

ref

Actuator face

Page 34: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 34August 2015

More Advanced Model-based Control

(MIMO) Multi-input (u)/Multi-output (y) require more complicated methods.

Each operation is a matrix operation…

The design principle is the same, but the difficulty of control design is balancing the tradeoffs between bandwidth, noise, and actuator saturation.

PC

Page 35: 0806ThermalSystems_webinar_2015

Copyright © 2015, SC Solutions, Inc. All Rights Reserved 35August 2015

Summary

A methodology for tuning a PID controller has been presented.

A model for the system is built and an approximate second order model is assumed for choosing the PID parameters.

A desired closed-loop transfer function (T) is selected. Here we selected it so the resulting controller has the same form as the PID.

The PID gains were derived that give the desired closed-loop response.

An example of implementing PID in COMSOL was presented.

Page 36: 0806ThermalSystems_webinar_2015

Poll QuestionHow do you model your feedback control law in your application?• In the same software that I use for studying

thermal systems.• In a different software that I interface.• In COMSOL using the Global ODEs and DAEs

interface.• In COMSOL with LiveLink™ for MATLAB®.• I haven't done it yet and want to learn how.

Page 37: 0806ThermalSystems_webinar_2015

DemoClosed-Loop Temperature Control of a Rapid

Thermal Processing (RTP) System

Wafer (200mm dia.)

Lamps (5)

FeedbackControllerS

Rapid Thermal Processor (RTP)

-+ Lamp

power commands

Wafer temperature sensors

Temperature reference

Cold walls600°C

1000°C

Page 38: 0806ThermalSystems_webinar_2015

Conclusion• A model-based method of tuning a PID

controller has been presented (Q-Design).

• For many thermal systems an open-loop step response can be used to characterize parameters (time-constants and gain) needed to tune the PID controller.

• COMSOL provides a simple ODE interface that allows quick implementation of PID control.

Page 39: 0806ThermalSystems_webinar_2015

Q&A Session

Page 40: 0806ThermalSystems_webinar_2015

Product Suite – COMSOL® 5.1

Page 41: 0806ThermalSystems_webinar_2015

Try COMSOL Multiphysics®North America Mansfield, MA San Antonio, TX Carlsbad, CA Pittsburgh, PA Boulder, CO Danbury, CT Beaverton, OR Tampa, FL Miami, FL

Europe Stuttgart, Germany Edinburgh, United Kingdom Zoetermeer, Netherlands Zürich, Switzerland Helsinki, Finland

Free hands-on workshops REGISTER TODAYwww.comsol.com/events

Page 42: 0806ThermalSystems_webinar_2015
Page 43: 0806ThermalSystems_webinar_2015

Contact Us• Questions?

www.comsol.com/contact

• www.comsol.com– User Stories– Videos– Application Gallery– Discussion Forum– Blog– Product News