Lecture 8 Turbocharger Modeling

Similar documents
GT RETROFIT AND RE-DESIGN OF A FOUR-STAGE CENTRIFUGAL COMPRESSOR

NODIA AND COMPANY. GATE SOLVED PAPER Chemical Engineering Plant Design and Economics. Copyright By NODIA & COMPANY

Partial Differential Equations of Fluid Dynamics

Riemannian Geometry, Key to Homework #1

BORGWARNER REPORTS SECOND QUARTER 2017 U.S. GAAP NET EARNINGS OF $1.00 PER DILUTED SHARE, OR $0.96 PER DILUTED SHARE EXCLUDING NON-COMPARABLE ITEMS

TECHNOLOGICAL PROGRESS AND ECONOMIC GROWTH. Dongpeng Liu Department of Economics Nanjing University

Case 3:17-md EMC Document Filed 01/18/19 Page 1 of 5. Exhibit 1 Class Member Benefits

"Pricing Exotic Options using Strong Convergence Properties

Kongsberg Automotive ASA. Fourth quarter 2015 February 12, 2016

Engineering 12 Physical systems analysis Laboratory 3: The Coupled Pendulums system

Value at Risk Ch.12. PAK Study Manual

Heterogeneous Firm, Financial Market Integration and International Risk Sharing

The Role of Firm-Level Productivity Growth for the Optimal Rate of Inflation

BORGWARNER REPORTS THIRD QUARTER 2018 U.S. GAAP NET EARNINGS OF $0.98 PER DILUTED SHARE, OR $1.00 PER DILUTED SHARE EXCLUDING NON-COMPARABLE ITEMS

MARCH 4-9, 2018 MACROECONOMIC MODELING OF ENERGY TRANSITIONS (ACCOUNTING FOR MONEY AND ENERGY)

Standard Motor Products, Inc. Q Investor Presentation

Compressed air energy savings in 4 moves

Rail Clamps Type BSZ. PINTSCH BUBENZER is certified according to DIN EN ISO 9001:2008. Safe Tried and Trusted High Performance Robust Easy Maintenance

ECON 815. A Basic New Keynesian Model II

Installed Pricing Information as of August, 2009

Dynamic Pricing in Ridesharing Platforms

Geology 3120 Powerpoint notes available online at:

DIN ISO 9001 : 2000 QA-Nr.: Safe Tried and Trusted High Performance Robust Easy Maintenance

Formula SAE Cost Report. Dalhousie University Formula SAE Michigan 2012 Car # 47

Foreign Competition and Banking Industry Dynamics: An Application to Mexico

Macroeconomics 2. Lecture 6 - New Keynesian Business Cycles March. Sciences Po

Monetary Macroeconomics Lecture 3. Mark Hayes

NO X CONFORMITY OF PRODUCTION PROCEDURE 2002 V07. Lloyd's Register of Shipping, 71 Fenchurch Street, London EC3M 4BS, UK

Standard Motor Products, Inc. Annual Meeting of Shareholders May 18 th, 2017

IN THE SEYCHELLES COURT OF APPEAL JUDGMENT OF THE COURT

JUNE 2018 COZY PARTS LIST

Short & Long Run impact of volatility on the effect monetary shocks

Chapter 8 Stocks and Their Valuation

Notice Inviting Tender

BorgWarner Inc. NEUTRAL ZACKS CONSENSUS ESTIMATES (BWA-NYSE)

Toolbox 7: Economic Feasibility Assessment Methods

IDENTIFYING BROAD AND NARROW FINANCIAL RISK FACTORS VIA CONVEX OPTIMIZATION: PART II

Investor Call Results H1 2017

Lightstream Scroll kW ULTRACOMPACT II COMPACT AIR-COOLED CHILLERS WITH SCROLL COMPRESSORS

Lars Davidson. Inlet boundary conditions for Hybrid LES-RANS

YAMAHA MOTOR CORPORATION, U.S.A. YAMALUBE ENGINE LUBRICATION LIMITED WARRANTY

ECONOMIC ANALYSIS AND LIFE CYCLE COSTING SECTION I

Green Income Tax Incentives. By Colleen M. Berndt, CPA Meyers Brothers Kalicka, P.C. MBKCPA.COM

Dynamic Relative Valuation

Investor Call Annual Results 2016

Global low carbon energy investment (2015) Low carbon energy R&D investment ( ) Sectoral renewable energy targets worldwide

Draft Measure Information Template Cooling Tower Efficiency and Turndown

On VIX Futures in the rough Bergomi model

CAREY W. KING, PH.D. JUNE 14, 2018 A LONG-TERM GROWTH MODEL WITH ENDOGENOUS BIOPHYSICAL AND ECONOMIC STATES

Prize offered for the solution of a dynamic blocking problem

Non-Renewable Resources and the Sustainability of the Economic Growth under Bilateral Trade

International Economics B 9. Monopolistic competition and international trade: Firm Heterogeneity

2012 Results and Strategy Review

Standard Motor Products, Inc. Q Investor Presentation

EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS

Easystream kW PACKAGED SCROLL CHILLERS

Econ 110: Introduction to Economic Theory. 11th Class 2/14/11

Pressure relief valve, pilot operated, Types ZDB and Z2DB

Lecture 12: New Economic Geography

International Trade Lecture 14: Firm Heterogeneity Theory (I) Melitz (2003)

Flow Dividers Application and Selection Notes

Location, Productivity, and Trade

$150 Xtreme Pro Mod, Nitrous Pro Street, Street Outlaw, Xtreme Street, N/A 10.5, Factory Super Cars, Chevrolet Performance Stock, Mopar Hemi Stock

Notes on the EM Algorithm Michael Collins, September 24th 2005

What is Cyclical in Credit Cycles?

Crashcourse Interest Rate Models

Chapter II: Labour Market Policy

Project 1: Double Pendulum

UNIFORM NETWORK CODE TRANSPORTATION PRINCIPAL DOCUMENT SECTION Y CHARGING METHODOLOGIES. 1. The Gas Transmission Transportation Charging Methodology

SOLUTIONS TO CHAPTER 10 EXERCISES: HOPPER DESIGN

CONTINUOUS TIME PRICING AND TRADING: A REVIEW, WITH SOME EXTRA PIECES

Microeconomic Theory II Preliminary Examination Solutions

UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions.

Financial Econometrics Jeffrey R. Russell. Midterm 2014 Suggested Solutions. TA: B. B. Deng

Modeling Yields at the Zero Lower Bound: Are Shadow Rates the Solution?

MOCK TEST PAPER INTERMEDIATE (IPC): GROUP I PAPER 3: COST ACCOUNTING AND FINANCIAL MANAGEMENT

MATH SPEAK - TO BE UNDERSTOOD AND MEMORIZED

Standard Motor Products, Inc. Q Investor Presentation

Supplemental Material Optics formula and additional results

STRESS INTENSITY FACTOR CALCULATIONS FOR CRACKS EMANATING FROM BOLT HOLES IN A JET ENGINE COMPRESSOR DISC

EXAMINATION : MACROECONOMICS (MAC) ECONOMICS 1 (ECO101)

9.1 Principal Component Analysis for Portfolios

Career Progression and Formal versus on the Job Training

arxiv: v1 [math.pr] 15 Dec 2011

arxiv: v1 [math.oc] 23 Dec 2010

Atlas Copco LTV angle nutrunners The force of productivity

Real Time Evolution of Non-Gaussian Cumulants in the QCD Critical Regime

General Equilibrium Approach to Evaluate Real Option Value of Reserved Environments

RISK-NEUTRAL VALUATION AND STATE SPACE FRAMEWORK. JEL Codes: C51, C61, C63, and G13

MATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS

Sample Questions. MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.

Introduction to Reinforcement Learning. MAL Seminar

Missing Growth. Philippe Aghion (College de France & LSE) Brown University - March 2018

4 Reinforcement Learning Basic Algorithms

QUALITATIVE AND QUANTITAVE RISK EVALUATION ON THE BASIS OF MILITARY AVIATION EVENT ANALYSIS

Model for rate of return to capital mathematical spiciness: ********** 10 stars (this appendix uses some advanced calculus) 1 Introduction

Earnings Release 1Q 2018 REVENUE GROWTH AND EXPANSION OF RETURN ON CAPITAL INVESTED

Case-study location: Battaglia Terme\Galzignano Terme

Kinetic Energy for a mass (m) at a velocity (V) (an Extensive Property) Frequently, one is interested in the change of KE from State 1 to State 2:

Application of Moment Expansion Method to Option Square Root Model

Transcription:

Lecture 8 Turbocharger Modeling System Modeling Institute for Dynamic Systems and Control (IDSC) Camillo Balerna Dr. Guillaume Ducard Camillo Balerna 07/11/2017 1

Lecture Overview Turbocharger modeling Turbine & Compressor Causality diagram Inputs, outputs Maps and operation Source: https://auto.howstuffworks.com/turbo2.htm Turbocharged internal combustion engines Engine power modeling Naturally aspirated (NA) Vs turbocharged (TC) engines Benefits of turbocharging F1 electrified turbocharger Camillo Balerna 07/11/2017 2

Turbocharger Θ tc dω tc dt = T t T c + T ext Source: https://auto.howstuffworks.com/turbo2.htm Camillo Balerna 07/11/2017 3

Turbine Power pressure T position Source: https://auto.howstuffworks.com/turbo2.htm Camillo Balerna 07/11/2017 4

Turbine Causality Diagram Inputs Outputs Inputs θ 3 : Temperature before the turbine [K] p 3 : pressure before the turbine [Pa] p 4 : pressure after the turbine [Pa] ω t : Turbine speed [rad/s] u vng : Variable nozzle geometry control input [ ] Π t = p 3 p 4 = p bef,t p aft,t Camillo Balerna 07/11/2017 5

Turbine Causality Diagram Inputs Outputs Outputs θ 4 : Temperature of the flow exiting the turbine [K] t : Mass flow through the turbine [Kg/s] T t : Torque generated by the turbine [Nm] Camillo Balerna 07/11/2017 6

Turbine Outputs o Temperature ofthe flow exiting the turbine θ 4 = θ 3 1 κ κ 1 η t 1 Π t o Mass Flow through the turbine Π t = p 3 p 4 = p bef,t p aft,t t = p 3 p ref,0 θ ref,0 θ 3 μሶ t o Torque produced by the turbine T t = P t ω t = t c p θ 3 ω t 1 Π t 1 κ κ η t Camillo Balerna 07/11/2017 7

Turbine Outputs derivation Open system de dt = ሶ H in ሶ H out ሶ W t + ሶ Q o Turbine does not store energy over time de dt = 0 o Turbine is assumed to be adiabatic no heat transfer ሶ Q = 0 P t = ሶ W t = ሶ H in ሶ H out = t c p θ 3 θ 4 Isentropic relation θ 3 θ 4,is = p 3 p 4 κ 1 κ κ 1 κ = Π t Turbine exit temperature θ 4 = θ 3 1 κ κ 1 η t 1 Π t Isentropic efficiency η t = θ 3 θ 4 θ 3 θ 4,is Turbine power produced P t = 1 κ κ t c p θ 3 1 Π t η t Camillo Balerna 07/11/2017 8

Turbine Outputs o Temperature ofthe flow exiting the turbine θ t = θ 3 1 κ κ 1 η t 1 Π t o Mass Flow through the turbine Π t = p 3 p 4 = p bef,t p aft,t t = p 3 p ref,0 θ ref,0 θ 3 μሶ t o Torque produced by the turbine T t = P t ω t = t c p θ 3 ω t 1 Π t 1 κ κ η t Camillo Balerna 07/11/2017 9

Turbine Efficiency Map Input Output c us = 2 c p θ 3 1 Π t (1 κ)/κ c us ǁ = r t ω t c us Since the turbine efficiency mainly depends on the angle of incidence of the inflowing gas, the turbine blade speed ratio c us ǁ is used as variable. Camillo Balerna 07/11/2017 10

Turbine Efficiency Map Variable Geometry Source: https://www.dieselnet.com/tech/air_turbo_vgt.php Camillo Balerna 07/11/2017 11

Turbine Outputs o Temperature ofthe flow exiting the turbine θ t = θ 3 1 κ κ 1 η t 1 Π t o Mass Flow through the turbine Π t = p 3 p 4 = p bef,t p aft,t t = p 3 p ref,0 θ ref,0 θ 3 μሶ t o Torque produced by the turbine T t = P t ω t = t c p θ 3 ω t 1 Π t 1 κ κ η t Camillo Balerna 07/11/2017 12

Turbine Mass Flow Map Input Output For control purposes, the mass flow behaviour of fluid-dynamic turbines can be modeled quite well as orifice compressible flow through a valve. If the turbine is a Variable Nozzle Turbine (VNT) or Variable Geometry Turbine (VGT), the mass flow and its maximum value depend on the nozzle position (as it is for the compressible flow through a valve). Camillo Balerna 07/11/2017 13

Turbine Mass Flow Map Fixed Geometry Source: https://www.dieselnet.com/tech/air_turbo_vgt.php Camillo Balerna 07/11/2017 14

Turbine Mass Flow Map Variable Geometry Source: https://www.dieselnet.com/tech/air_turbo_vgt.php Camillo Balerna 07/11/2017 15

Turbine Variable Geometry Turbine (VGT) At Low Engine Speed low mass flow low pressure low turbine power Narrow inlet area Better incidence Increased efficiency Increased power At High Engine Speed high mass flow high pressure Wide inlet area Avoid choke! Source: https://www.intmarketing.org/en/automotive/113-variable-turbine-geometry.html Camillo Balerna 07/11/2017 16

Compressor Source: https://auto.howstuffworks.com/turbo2.htm Camillo Balerna 07/11/2017 17

Compressor pressure C Power position Source: https://auto.howstuffworks.com/turbo2.htm Camillo Balerna 07/11/2017 18

Compressor Causality Diagram Inputs Outputs Inputs p 1 : Pressure before the compressor [Pa] p 2 : Pressure after the compressor [Pa] θ 1 : Temperature before the compressor [K] ω c : Compressor speed [rad/s] Π c = p 2 p 1 = p aft,c p bef,c Camillo Balerna 07/11/2017 19

Compressor Causality Diagram Inputs Outputs Outputs θ c : Temperature of the flow exiting the compressor [K] c : Mass flow through the compressor [Kg/s] T c : Torque absorbed by the compressor [Nm] Camillo Balerna 07/11/2017 20

Compressor Outputs o Temperature ofthe flow exiting the compressor θ c = θ 1 + Π c κ 1 κ 1 θ 1 η c o Mass Flow through the compressor Π c = p 2 p 1 = p aft,c p bef,c c = p 1 p ref,0 θ ref,0 θ 1 μሶ c o Torque absorbed by the compressor T c = P c ω c = c c p θ 1 ω c Π c κ 1 κ 1 1 η c Camillo Balerna 07/11/2017 21

Compressor Outputs derivation Open system de dt = ሶ H in ሶ H out ሶ W c + ሶ Q o Compressor does not store energy over time de dt = 0 o Compressor is assumed to be adiabatic no heat transfer ሶ Q = 0 P c = ሶ W c = Hሶ out Hሶ in = c c p θ 2 θ 1 Isentropic relation θ 2,is θ 1 = p 2 p 1 κ 1 κ Isentropic efficiency η c = θ 2,is θ 1 θ 2 θ 1 κ 1 κ = Π c Compressor exit temperature θ 2 = θ 1 + Π c κ 1 κ 1 θ 1 η c Compressor power absorbed P c = κ 1 c c p θ 1 Π κ c 1 1 η c Camillo Balerna 07/11/2017 22

Compressor Outputs o Temperature ofthe flow exiting the compressor θ c = θ 1 + Π c κ 1 κ 1 θ 1 η c o Mass Flow through the compressor Π c = p 2 p 1 = p aft,c p bef,c c = p 1 p ref,0 θ ref,0 θ 1 μሶ c o Torque absorbed by the compressor T c = P c ω c = c c p θ 1 ω c Π c κ 1 κ 1 1 η c Camillo Balerna 07/11/2017 23

Compressor Mass Flow & Efficiency Map Input Output Camillo Balerna 07/11/2017 24

Compressor Mass Flow & Efficiency Map Source: http://www.enginelabs.com/engine-tech/poweradders/understanding-compressor-maps-sizing-a-turbocharger/ Camillo Balerna 07/11/2017 25

Compressor Operational Limits Fluid-dynamic instabilities destroy the regular flow pattern possible back-flow Maximum speed allowed to avoid mechanical damages centrifugal forces Behaviour at zero (or very low) speed blocking orifice Flow reaches sonics conditions choked orifice Camillo Balerna 07/11/2017 26

Formula 1 Turbocharged Engine Source: https://sport.sky.it/formula1/2017/03/21/formula-1--il-dizionario--power-unit-ed-elettronica.html Camillo Balerna 07/11/2017 27

Internal Combustion Engine Engine Power can be approximated as following: P engine = P comb,fuel + P fric + P pump Assume constant engine speed ω e o Engine Power coming from the fuel combustion: P comb,fuel = e comb P f = e comb H l fuel k 1 fuel o Engine Power coming from the pistons mechanical friction: P fric k 2 o Engine Power coming from the gas exchange: P pump = p intake p exhaust V d ω e 4π k 3 p intake p exhaust Camillo Balerna 07/11/2017 28

Internal Combustion Engine o Engine Power, neglecting the friction and assuming p intake = p exhaust : P engine k 1 fuel + k 3 p intake p exhaust k 1 fuel o Air to Fuel Ratio is defined as following: λ AF = air 1 fuel σ 0 λ AF = 1 fuel = air σ 0 o Engine Air Mass Flow isapproximated as following: air = p intake R air θ intake ω e 4π V d λ vol k 4 p intake P engine fuel air p intake NA 1 bar 100 kw TC 4 bar 400 kw Camillo Balerna 07/11/2017 29

Turbocharged Engine o Engine Power, neglecting the friction and for a specific P engine k 1 fuel : fuel + k 3 p intake p exhaust NA 1 bar TC 4 bar NA 1 bar TC 2 bar C intake exhaust T ω e = 10 000 rpm V d = 1. 6 L P pump,na 0 kw P pump,tc 26 kw Camillo Balerna 07/11/2017 30

Engine Response P engine P engine,req Naturally Aspirated Slower increase time Camillo Balerna 07/11/2017 31

Formula 1 Engine Response How and how fast is the engine power response of a conventional turbocharger compared to an electrified turbocharger (e.g. F1)? Turbocharger 1) Press throttle pedal 2) More fuel injected 3) Exhaust temperature increases 4) More power extracted by the turbine 5) Turbocharger speed increases 6) Compressor mass flow increases 7) More air more fuel can be injected 8) More Engine power Electrified Turbocharger 1) Press throttle pedal 2) More fuel injected & MGU-H positive 3) Turbocharger speed increases faster ( MGU-H & exhaust temperature) 4) Compressor mass flow increases 5) More air more fuel can be injected 6) More Engine power Camillo Balerna 07/11/2017 32

Formula 1 Engine Response P engine P engine,req Steeper increase time Camillo Balerna 07/11/2017 33

IDSC Open Lab 2017 Camillo Balerna 07/11/2017 34

IDSC Open Lab 2017 Formula 1 Power Unit Efficient control algorithms are designed for the hybrid electric propulsion system of the Formula 1 car, in order to achieve the fastest possible lap-time. (Presentations in English or German) Mauro Salazar, maurosalazar@idsc.mavt.ethz.ch Camillo Balerna, balernac@idsc.mavt.ethz.ch ML K37.1 http://www.idsc.ethz.ch/research-guzzellaonder/research-projects/formula1.html Camillo Balerna 07/11/2017 35