TOWARDS GRADIENT BASED AERODYNAMIC OPTIMIZATION OF WIND TURBINE BLADES USING OVERSET GRIDS

Size: px
Start display at page:

Download "TOWARDS GRADIENT BASED AERODYNAMIC OPTIMIZATION OF WIND TURBINE BLADES USING OVERSET GRIDS"

Transcription

1 TOWARDS GRADIENT BASED AERODYNAMIC OPTIMIZATION OF WIND TURBINE BLADES USING OVERSET GRIDS S. H. Jongsm E. T. A. vn de Weide H. W. M. Hoeijmkers Overset symposium Deprtment of mechnicl engineering 1

2 OUTLINE Motivtion Optimiztion procedure Prmeteriztion method Flow model Hole cutting Grdient computtion Test cse nd results 2

3 MOTIVATION (1 Conventionl pproch in design of lde: Design 2D erofoil sections Comine sections to otin lde geometry Incorporte effects 3D flow fetures very difficult even for experienced designers Current pproch: Simultion sed design of lde shpe in 3D 3D flow fetures tken into ccount Design expertise required for defining ojective function nd constrints for optimiztion method 3

4 GRADIENT BASED OPTIMIZATION PROCEDURE 4

5 FLOW MODEL (1 Consider rigid isolted rotor: i.e. effects of presence of tower nd ground plne re neglected effect of lde deformtion is neglected Assume rotor plne perpendiculr to stedy wind: flow is periodic Stedy flow round single lde in rotting frme i.e. effect of unstedy inflow nd yw is not considered Currently solve Euler equtions Lter: RANS 5

6 FLOW MODEL (2 Cell centered finite volume method 2 nd Order sptil ccurcy Integrtion to stedy stte y: Newton s method Runge-Kutt time stepping Use PETSc to solve the resulting system of liner equtions Re c = M = 0.2 = 2.3 o O-grid nd Crtesin ckground grid 6

7 PARAMETERIZATION METHOD Requirements: Spn complete design spce Limit numer of design vriles Choice: Non-uniform rtionl sis spline (NURBS surfce Flexile Comptile with CAD Design vriles: (Selection of coordintes of control points Weights (if desired 7

8 FLOW DOMAIN DISCRETIZATION Requirements: Account for chnge of shpe fter ech design itertion Not computtionlly intensive Fully utomtic Choice: Multi-lock overset grids Use hyperolic grid genertion for field grids 8

9 GRID PREPROCESSING (1 Approch Evlution of surfce integrls Zipper grid (Chn 2009 Elimintion cells inside geometricl entities Ry csting Block connectivity Implicit hole cutting (Lee 2008 No sequentil ottleneck Resonly efficient prllel implementtion Chn (2009, Enhncements to the Hyrid Mesh Approch to Surfce Lods Integrtion on Overset Structured Grids Lee (2008, On Overset Grids Connectivity nd Vortex Trcking in Rotorcrft CFD (PhD thesis 9

10 GRID PREPROCESSING (2 ZIPPER GRID (1 Closed non-overlpping surfce grid required for: Evlution of surfce integrls Ry-csting procedure Achieved y ppliction of zipper grids: remove overlpping prt of surfce grids connect neighoring grids y creting tringles in etween Works for prolems considered so fr 10

11 GRID PREPROCESSING (3 ZIPPER GRID (2 11

12 GRID PREPROCESSING (4 RAY-CASTING Approch: Choose p nd csting direction Find possile mtches using tree serch Check possile mtches with ccurte check Accurte check: Split qudrilterl in 2 tringles Trnslte tringle to origin: Project p' nd tringle on fce perpendiculr to csting direction Use trnsformtion to rycentric coordintes: Solve: Ry crosses tringle if: 12

13 GRID PREPROCESSING (5 IMPLICIT HOLE CUTTING (1 Concept Use dul grid for hole cutting Identify cells tht reside in region covered y multiple locks In region of overlp: Consider user defined priority of locks Consider index distnce to wll of cells Consider cell volume Use cells with the smllest volume to solve governing equtions (field cell Cell with lrger volume gets sttus: hole cell Specify fringe cells t interfce etween hole nd field cells (stencil dependent 13

14 GRID PREPROCESSING (6 IMPLICIT HOLE CUTTING (2 Approch Construct lternting digitl tree for ll cells in ech lock Perform tree serch to look for cells for which the ounding ox overlps prticulr cell Use tri-liner trnsformtion to identify overlp for potentil cndidtes from tree serch Use criteri on lock priority, index distnce nd cell volume to identify field nd hole cells Identify fringe cells 14

15 GRADIENT COMPUTATION (1 INTRODUCTION (1 Requirements: Accurte Otined efficiently Methods: Finite difference Complex step finite difference Discrete djoint eqution method 15

16 GRADIENT COMPUTATION (2 DISCRETE ADJOINT EQUATION METHOD (1 Design vrile α Computtionl mesh X = X (α Flow vriles u = u(α,x Ojective function I = I(α,u, X Residul of flow solution R = R(α,u, X Totl derivtive of ojective function: 16

17 GRADIENT COMPUTATION (3 DISCRETE ADJOINT EQUATION METHOD (2 Totl derivtive of residul: Rewrite to: Sustitute: 17

18 GRADIENT COMPUTATION (4 DISCRETE ADJOINT EQUATION METHOD (3 Define djoint vector: Sustitute to otin new expression for totl derivtive: Advntge: Only one flow solution needed for ll derivtives, once djoint vector is found y solving: 18

19 GRADIENT COMPUTATION (5 DUAL NUMBERS (1 (Fike et l Choose h = 1 Fike & Alonso (2011, The Development of Hyper-Dul Numers for Exct Second- Derivtive Clcultions, AIAA pper

20 20 GRADIENT COMPUTATION (6 DUAL NUMBERS (2 / / ( / 1/( ( /( ( / 1/ d 1/(c Division : ( ( ( Multipliction : ( ( ( ( Addition : 2 2 c d c c d c d c c d c c d c d c d c d c = = = = = ] ln( exp[( ( / ln( ln( exp( exp( exp( 1] ( [tn tn( tn( cos( sin( sin( sin( cos( Mth functions : cos( d (c 2 d c = = = = = = Computtionl rules

21 GRADIENT COMPUTATION (7 CONSIDERATIONS Limit implementtion time y reusing existing functions Use of C templte functions Use grph coloring to limit numer of function evlutions Existence of fringe cells nd hlo cells requires specil ttention Tke effect of oundry conditions into ccount explicitly Determine corresponding mtrix index of donors Computtion djoint vector Use first order Jcoin for construction ILU(k preconditioning mtrix Use PETSc to solve system of liner equtions y mens of GMRES itertions 21

22 GRADIENT COMPUTATION (8 HANDLING FRINGE CELLS Construction Jcoin mtrix: Fringe cells re represented y row in the mtrix Mke sure ech fringe cell hs explicit representtion y field cells Required to fcilitte implicit solution method Updte solution fter prtil convergence GMRES itertions For consistency, fringe cells re updted using solution of field cells 22

23 OPTIMIZATION (1 APPROACH Choose design vriles Specify ojective function nd constrint functions Specify ounds on design vriles Evlute ojective function nd constrint functions Compute grdients Provide vlue of ojective function, constrint functions nd grdients to grdient sed optimiztion lgorithm: SNOPT (sprse non-liner optimizer (Stnford University (Gill et l BFGS qusi Newton method Provides new design vriles o Repet procedure Gill, Murry & Sunders (2005, SNOPT: An SQP Algorithm for Lrge-Scle Constrined Optimiztion, SIAM Review, Vol. 47, No

24 OPTIMIZATION (2 PROBLEM DEFINITION (1 Inverse(/reverse design: Ojective: Find irfoil geometry which minimizes (c p c p_trget 2 Approch: Red input file with c p dt nd corresponding x-coordinte Interpolte input dt for comprison of c p current design for different x-coordinte center of fce Sum squre of difference Minimize resulting sum x 24

25 OPTIMIZATION (3 PROBLEM DEFINITION (2 Constrint functions: 25

26 OPTIMIZATION (4 SETUP Initil design: NACA 0012 t 0 ngle of ttck 13 control points 36 design vriles Flow conditions: M = 0.3 Grid dimensions: Grid n i n j O-grid Bckground Fr field: At 50 chord lengths 26

27 OPTIMIZATION (5 RESULTS Trget nd initil c p -distriution Finl nd initil geometry Fit Originl dt point 27

28 CONCLUDING REMARKS Prmeteriztion using NURBS surfce Discretiztion using hyperolic grid genertion method Fully utomtic hole cutting procedure Derivtives computed using: Discrete djoint eqution method Dul numer rithmetic Method successfully pplied for solving 2D optimiztion prolems on overset grids Fit of c p distriution Drg minimiztion in trnsonic flow 28

29 FURTURE WORK Testing the method for simple 3D optimiztion prolem Applying the method for full 3D wind turine lde in co-rotting frme of reference Use RANS equtions Extension to time-periodic optimiztion vi time-spectrl method 29

30

6.3 Volumes. Just as area is always positive, so is volume and our attitudes towards finding it.

6.3 Volumes. Just as area is always positive, so is volume and our attitudes towards finding it. 6.3 Volumes Just s re is lwys positive, so is volume nd our ttitudes towrds finding it. Let s review how to find the volume of regulr geometric prism, tht is, 3-dimensionl oject with two regulr fces seprted

More information

Lecture 10 Evolutionary Computation: Evolution strategies and genetic programming

Lecture 10 Evolutionary Computation: Evolution strategies and genetic programming Lecture 10 Evolutionry Computtion: Evolution strtegies nd genetic progrmming Evolution strtegies Genetic progrmming Summry Negnevitsky, Person Eduction, 2011 1 Evolution Strtegies Another pproch to simulting

More information

Stained Glass Design. Teaching Goals:

Stained Glass Design. Teaching Goals: Stined Glss Design Time required 45-90 minutes Teching Gols: 1. Students pply grphic methods to design vrious shpes on the plne.. Students pply geometric trnsformtions of grphs of functions in order to

More information

GENERATING ORTHOIMAGES FOR CLOSE-RANGE OBJECTS BY AUTOMATICALLY DETECTING BREAKLINES

GENERATING ORTHOIMAGES FOR CLOSE-RANGE OBJECTS BY AUTOMATICALLY DETECTING BREAKLINES GENEATING OTHOIMAGES FO CLOSE-ANGE OBJECTS BY AUTOMATICALLY DETECTING BEAKLINES Efstrtios Stylinidis 1, Lzros Sechidis 1, Petros Ptis 1, Spiros Sptls 2 Aristotle University of Thessloniki 1 Deprtment of

More information

Math 35 Review Sheet, Spring 2014

Math 35 Review Sheet, Spring 2014 Mth 35 Review heet, pring 2014 For the finl exm, do ny 12 of the 15 questions in 3 hours. They re worth 8 points ech, mking 96, with 4 more points for netness! Put ll your work nd nswers in the provided

More information

2 Computing all Intersections of a Set of Segments Line Segment Intersection

2 Computing all Intersections of a Set of Segments Line Segment Intersection 15-451/651: Design & Anlysis of Algorithms Novemer 14, 2016 Lecture #21 Sweep-Line nd Segment Intersection lst chnged: Novemer 8, 2017 1 Preliminries The sweep-line prdigm is very powerful lgorithmic design

More information

8.2 Areas in the Plane

8.2 Areas in the Plane 39 Chpter 8 Applictions of Definite Integrls 8. Ares in the Plne Wht ou will lern out... Are Between Curves Are Enclosed Intersecting Curves Boundries with Chnging Functions Integrting with Respect to

More information

a < a+ x < a+2 x < < a+n x = b, n A i n f(x i ) x. i=1 i=1

a < a+ x < a+2 x < < a+n x = b, n A i n f(x i ) x. i=1 i=1 Mth 33 Volume Stewrt 5.2 Geometry of integrls. In this section, we will lern how to compute volumes using integrls defined by slice nlysis. First, we recll from Clculus I how to compute res. Given the

More information

B. Definition: The volume of a solid of known integrable cross-section area A(x) from x = a

B. Definition: The volume of a solid of known integrable cross-section area A(x) from x = a Mth 176 Clculus Sec. 6.: Volume I. Volume By Slicing A. Introduction We will e trying to find the volume of solid shped using the sum of cross section res times width. We will e driving towrd developing

More information

Before We Begin. Introduction to Spatial Domain Filtering. Introduction to Digital Image Processing. Overview (1): Administrative Details (1):

Before We Begin. Introduction to Spatial Domain Filtering. Introduction to Digital Image Processing. Overview (1): Administrative Details (1): Overview (): Before We Begin Administrtive detils Review some questions to consider Winter 2006 Imge Enhncement in the Sptil Domin: Bsics of Sptil Filtering, Smoothing Sptil Filters, Order Sttistics Filters

More information

Explicit Decoupled Group Iterative Method for the Triangle Element Solution of 2D Helmholtz Equations

Explicit Decoupled Group Iterative Method for the Triangle Element Solution of 2D Helmholtz Equations Interntionl Mthemticl Forum, Vol. 12, 2017, no. 16, 771-779 HIKARI Ltd, www.m-hikri.com https://doi.org/10.12988/imf.2017.7654 Explicit Decoupled Group Itertive Method for the Tringle Element Solution

More information

Math 17 - Review. Review for Chapter 12

Math 17 - Review. Review for Chapter 12 Mth 17 - eview Ying Wu eview for hpter 12 1. Given prmetric plnr curve x = f(t), y = g(t), where t b, how to eliminte the prmeter? (Use substitutions, or use trigonometry identities, etc). How to prmeterize

More information

Algorithm Design (5) Text Search

Algorithm Design (5) Text Search Algorithm Design (5) Text Serch Tkshi Chikym School of Engineering The University of Tokyo Text Serch Find sustring tht mtches the given key string in text dt of lrge mount Key string: chr x[m] Text Dt:

More information

Optimal aeroacoustic shape design using approximation modeling

Optimal aeroacoustic shape design using approximation modeling Center for Turbulence Reserch Annul Reserch Briefs 22 21 Optiml erocoustic shpe design using pproximtion modeling By Alison L. Mrsden, Meng Wng AND Petros Koumoutskos 1. Introduction Reduction of noise

More information

Slides for Data Mining by I. H. Witten and E. Frank

Slides for Data Mining by I. H. Witten and E. Frank Slides for Dt Mining y I. H. Witten nd E. Frnk Simplicity first Simple lgorithms often work very well! There re mny kinds of simple structure, eg: One ttriute does ll the work All ttriutes contriute eqully

More information

ZZ - Advanced Math Review 2017

ZZ - Advanced Math Review 2017 ZZ - Advnced Mth Review Mtrix Multipliction Given! nd! find the sum of the elements of the product BA First, rewrite the mtrices in the correct order to multiply The product is BA hs order x since B is

More information

PARALLEL AND DISTRIBUTED COMPUTING

PARALLEL AND DISTRIBUTED COMPUTING PARALLEL AND DISTRIBUTED COMPUTING 2009/2010 1 st Semester Teste Jnury 9, 2010 Durtion: 2h00 - No extr mteril llowed. This includes notes, scrtch pper, clcultor, etc. - Give your nswers in the ville spce

More information

Approximation by NURBS with free knots

Approximation by NURBS with free knots pproximtion by NURBS with free knots M Rndrinrivony G Brunnett echnicl University of Chemnitz Fculty of Computer Science Computer Grphics nd Visuliztion Strße der Ntionen 6 97 Chemnitz Germny Emil: mhrvo@informtiktu-chemnitzde

More information

Geometric transformations

Geometric transformations Geometric trnsformtions Computer Grphics Some slides re bsed on Shy Shlom slides from TAU mn n n m m T A,,,,,, 2 1 2 22 12 1 21 11 Rows become columns nd columns become rows nm n n m m A,,,,,, 1 1 2 22

More information

A dual of the rectangle-segmentation problem for binary matrices

A dual of the rectangle-segmentation problem for binary matrices A dul of the rectngle-segmenttion prolem for inry mtrices Thoms Klinowski Astrct We consider the prolem to decompose inry mtrix into smll numer of inry mtrices whose -entries form rectngle. We show tht

More information

If f(x, y) is a surface that lies above r(t), we can think about the area between the surface and the curve.

If f(x, y) is a surface that lies above r(t), we can think about the area between the surface and the curve. Line Integrls The ide of line integrl is very similr to tht of single integrls. If the function f(x) is bove the x-xis on the intervl [, b], then the integrl of f(x) over [, b] is the re under f over the

More information

N-Level Math (4045) Formula List. *Formulas highlighted in yellow are found in the formula list of the exam paper. 1km 2 =1000m 1000m

N-Level Math (4045) Formula List. *Formulas highlighted in yellow are found in the formula list of the exam paper. 1km 2 =1000m 1000m *Formul highlighted in yellow re found in the formul lit of the em pper. Unit Converion Are m =cm cm km =m m = m = cm Volume m =cm cm cm 6 = cm km/h m/ itre =cm (ince mg=cm ) 6 Finncil Mth Percentge Incree

More information

USING HOUGH TRANSFORM IN LINE EXTRACTION

USING HOUGH TRANSFORM IN LINE EXTRACTION Stylinidis, Efstrtios USING HOUGH TRANSFORM IN LINE EXTRACTION Efstrtios STYLIANIDIS, Petros PATIAS The Aristotle University of Thessloniki, Deprtment of Cdstre Photogrmmetry nd Crtogrphy Univ. Box 473,

More information

Network Interconnection: Bridging CS 571 Fall Kenneth L. Calvert All rights reserved

Network Interconnection: Bridging CS 571 Fall Kenneth L. Calvert All rights reserved Network Interconnection: Bridging CS 57 Fll 6 6 Kenneth L. Clvert All rights reserved The Prolem We know how to uild (rodcst) LANs Wnt to connect severl LANs together to overcome scling limits Recll: speed

More information

Complete Coverage Path Planning of Mobile Robot Based on Dynamic Programming Algorithm Peng Zhou, Zhong-min Wang, Zhen-nan Li, Yang Li

Complete Coverage Path Planning of Mobile Robot Based on Dynamic Programming Algorithm Peng Zhou, Zhong-min Wang, Zhen-nan Li, Yang Li 2nd Interntionl Conference on Electronic & Mechnicl Engineering nd Informtion Technology (EMEIT-212) Complete Coverge Pth Plnning of Mobile Robot Bsed on Dynmic Progrmming Algorithm Peng Zhou, Zhong-min

More information

Summer Review Packet For Algebra 2 CP/Honors

Summer Review Packet For Algebra 2 CP/Honors Summer Review Pcket For Alger CP/Honors Nme Current Course Mth Techer Introduction Alger uilds on topics studied from oth Alger nd Geometr. Certin topics re sufficientl involved tht the cll for some review

More information

Lecture 7: Integration Techniques

Lecture 7: Integration Techniques Lecture 7: Integrtion Techniques Antiderivtives nd Indefinite Integrls. In differentil clculus, we were interested in the derivtive of given rel-vlued function, whether it ws lgeric, eponentil or logrithmic.

More information

Fig.1. Let a source of monochromatic light be incident on a slit of finite width a, as shown in Fig. 1.

Fig.1. Let a source of monochromatic light be incident on a slit of finite width a, as shown in Fig. 1. Answer on Question #5692, Physics, Optics Stte slient fetures of single slit Frunhofer diffrction pttern. The slit is verticl nd illuminted by point source. Also, obtin n expression for intensity distribution

More information

Ray surface intersections

Ray surface intersections Ry surfce intersections Some primitives Finite primitives: polygons spheres, cylinders, cones prts of generl qudrics Infinite primitives: plnes infinite cylinders nd cones generl qudrics A finite primitive

More information

Tree Structured Symmetrical Systems of Linear Equations and their Graphical Solution

Tree Structured Symmetrical Systems of Linear Equations and their Graphical Solution Proceedings of the World Congress on Engineering nd Computer Science 4 Vol I WCECS 4, -4 October, 4, Sn Frncisco, USA Tree Structured Symmetricl Systems of Liner Equtions nd their Grphicl Solution Jime

More information

MTH 146 Conics Supplement

MTH 146 Conics Supplement 105- Review of Conics MTH 146 Conics Supplement In this section we review conics If ou ne more detils thn re present in the notes, r through section 105 of the ook Definition: A prol is the set of points

More information

Math 142, Exam 1 Information.

Math 142, Exam 1 Information. Mth 14, Exm 1 Informtion. 9/14/10, LC 41, 9:30-10:45. Exm 1 will be bsed on: Sections 7.1-7.5. The corresponding ssigned homework problems (see http://www.mth.sc.edu/ boyln/sccourses/14f10/14.html) At

More information

Section 9.2 Hyperbolas

Section 9.2 Hyperbolas Section 9. Hperols 597 Section 9. Hperols In the lst section, we lerned tht plnets hve pproimtel ellipticl orits round the sun. When n oject like comet is moving quickl, it is le to escpe the grvittionl

More information

4452 Mathematical Modeling Lecture 4: Lagrange Multipliers

4452 Mathematical Modeling Lecture 4: Lagrange Multipliers Mth Modeling Lecture 4: Lgrnge Multipliers Pge 4452 Mthemticl Modeling Lecture 4: Lgrnge Multipliers Lgrnge multipliers re high powered mthemticl technique to find the mximum nd minimum of multidimensionl

More information

Graphing Conic Sections

Graphing Conic Sections Grphing Conic Sections Definition of Circle Set of ll points in plne tht re n equl distnce, clled the rdius, from fixed point in tht plne, clled the center. Grphing Circle (x h) 2 + (y k) 2 = r 2 where

More information

A modified directed search domain algorithm for multiobjective engineering and design optimization

A modified directed search domain algorithm for multiobjective engineering and design optimization Struct Multidisc Optim DOI.7/s58-3-946- RESEARCH PAPER A modified directed serch domin lgorithm for multiojective engineering nd design optimiztion Tohid Erfni Sergey V. Utyuzhnikov Brin Kolo Received:

More information

Linear analysis of texture property relationships using process-based representations of Rodrigues space

Linear analysis of texture property relationships using process-based representations of Rodrigues space Act Mterili 55 (7) 573 587 www.ctmt-journls.com Liner nlysis of texture property reltionships using process-sed representtions of Rodrigues spce Veer Sundrrghvn, Nichols Zrs * Mterils Process Design nd

More information

Agilent Mass Hunter Software

Agilent Mass Hunter Software Agilent Mss Hunter Softwre Quick Strt Guide Use this guide to get strted with the Mss Hunter softwre. Wht is Mss Hunter Softwre? Mss Hunter is n integrl prt of Agilent TOF softwre (version A.02.00). Mss

More information

Objective: Students will understand what it means to describe, graph and write the equation of a parabola. Parabolas

Objective: Students will understand what it means to describe, graph and write the equation of a parabola. Parabolas Pge 1 of 8 Ojective: Students will understnd wht it mens to descrie, grph nd write the eqution of prol. Prols Prol: collection of ll points P in plne tht re the sme distnce from fixed point, the focus

More information

Optimization of Air Bearing Slider Design

Optimization of Air Bearing Slider Design Proceedings of TC2005 orld Tribology Congress III Proceedings of TC2005 September 2-6, orld 2005, Tribology shington, Congress D.C., III SA September 2-6, 2005, shington, D.C., SA Optimiztion of Air Bering

More information

Chapter44. Polygons and solids. Contents: A Polygons B Triangles C Quadrilaterals D Solids E Constructing solids

Chapter44. Polygons and solids. Contents: A Polygons B Triangles C Quadrilaterals D Solids E Constructing solids Chpter44 Polygons nd solids Contents: A Polygons B Tringles C Qudrilterls D Solids E Constructing solids 74 POLYGONS AND SOLIDS (Chpter 4) Opening prolem Things to think out: c Wht different shpes cn you

More information

A multiview 3D modeling system based on stereo vision techniques

A multiview 3D modeling system based on stereo vision techniques Mchine Vision nd Applictions (2005) 16: 148 156 Digitl Oject Identifier (DOI) 10.1007/s00138-004-0165-2 Mchine Vision nd Applictions A multiview 3D modeling system sed on stereo vision techniques Soon-Yong

More information

Applications of the Definite Integral ( Areas and Volumes)

Applications of the Definite Integral ( Areas and Volumes) Mth1242 Project II Nme: Applictions of the Definite Integrl ( Ares nd Volumes) In this project, we explore some pplictions of the definite integrl. We use integrls to find the re etween the grphs of two

More information

Topic 3: 2D Transformations 9/10/2016. Today s Topics. Transformations. Lets start out simple. Points as Homogeneous 2D Point Coords

Topic 3: 2D Transformations 9/10/2016. Today s Topics. Transformations. Lets start out simple. Points as Homogeneous 2D Point Coords Tody s Topics 3. Trnsformtions in 2D 4. Coordinte-free geometry 5. (curves & surfces) Topic 3: 2D Trnsformtions 6. Trnsformtions in 3D Simple Trnsformtions Homogeneous coordintes Homogeneous 2D trnsformtions

More information

Unsteady Adjoint Method for the Optimal Control of Advection and Burger s Equations using High-Order Spectral Difference Method

Unsteady Adjoint Method for the Optimal Control of Advection and Burger s Equations using High-Order Spectral Difference Method 49th AIAA Aerospce Sciences Meeting including the New Horizons Forum nd Aerospce Eposi 4-7 Jnur 2, Orlndo, Florid AIAA 2-24 Unsted Adjoint Method for the Optiml Control of Advection nd Burger s Equtions

More information

Parallel Square and Cube Computations

Parallel Square and Cube Computations Prllel Squre nd Cube Computtions Albert A. Liddicot nd Michel J. Flynn Computer Systems Lbortory, Deprtment of Electricl Engineering Stnford University Gtes Building 5 Serr Mll, Stnford, CA 945, USA liddicot@stnford.edu

More information

Chapter 2 Sensitivity Analysis: Differential Calculus of Models

Chapter 2 Sensitivity Analysis: Differential Calculus of Models Chpter 2 Sensitivity Anlysis: Differentil Clculus of Models Abstrct Models in remote sensing nd in science nd engineering, in generl re, essentilly, functions of discrete model input prmeters, nd/or functionls

More information

CS-C3100 Computer Graphics, Fall 2016 Ray Casting II Intersection Extravaganza

CS-C3100 Computer Graphics, Fall 2016 Ray Casting II Intersection Extravaganza CS-C3100 Computer Grphics, Fll 2016 Ry Csting II Intersection Extrvgnz Henrik Wnn Jensen Jkko Lehtinen with lots of mteril from Frédo Durnd CS-C3100 Fll 2016 Lehtinen 1 Pinholes in Nture Flickr user Picture

More information

Analysis of Computed Diffraction Pattern Diagram for Measuring Yarn Twist Angle

Analysis of Computed Diffraction Pattern Diagram for Measuring Yarn Twist Angle Textiles nd Light ndustril Science nd Technology (TLST) Volume 3, 2014 DO: 10.14355/tlist.2014.0301.01 http://www.tlist-journl.org Anlysis of Computed Diffrction Pttern Digrm for Mesuring Yrn Twist Angle

More information

Ranking of Hexagonal Fuzzy Numbers for Solving Multi Objective Fuzzy Linear Programming Problem

Ranking of Hexagonal Fuzzy Numbers for Solving Multi Objective Fuzzy Linear Programming Problem Interntionl Journl of Computer pplictions 097 8887 Volume 8 No 8 Decemer 0 nking of egonl Fuzzy Numers Solving ulti Ojective Fuzzy Liner Progrmming Prolem jrjeswri. P Deprtment of themtics Chikknn Government

More information

COLOUR IMAGE MATCHING FOR DTM GENERATION AND HOUSE EXTRACTION

COLOUR IMAGE MATCHING FOR DTM GENERATION AND HOUSE EXTRACTION Hee Ju Prk OLOUR IMAGE MATHING FOR DTM GENERATION AND HOUSE EXTRATION Hee Ju PARK, Petr ZINMMERMANN * Swiss Federl Institute of Technology, Zuric Switzerlnd Institute for Geodesy nd Photogrmmetry heeju@ns.shingu-c.c.kr

More information

Section 5.3 : Finding Area Between Curves

Section 5.3 : Finding Area Between Curves MATH 9 Section 5. : Finding Are Between Curves Importnt: In this section we will lern just how to set up the integrls to find re etween curves. The finl nswer for ech emple in this hndout is given for

More information

LU Decomposition. Mechanical Engineering Majors. Authors: Autar Kaw

LU Decomposition. Mechanical Engineering Majors. Authors: Autar Kaw LU Decomposition Mechnicl Engineering Mjors Authors: Autr Kw Trnsforming Numericl Methods Eduction for STEM Undergrdutes // LU Decomposition LU Decomposition LU Decomposition is nother method to solve

More information

CS481: Bioinformatics Algorithms

CS481: Bioinformatics Algorithms CS481: Bioinformtics Algorithms Cn Alkn EA509 clkn@cs.ilkent.edu.tr http://www.cs.ilkent.edu.tr/~clkn/teching/cs481/ EXACT STRING MATCHING Fingerprint ide Assume: We cn compute fingerprint f(p) of P in

More information

Grade 7/8 Math Circles Geometric Arithmetic October 31, 2012

Grade 7/8 Math Circles Geometric Arithmetic October 31, 2012 Fculty of Mthemtics Wterloo, Ontrio N2L 3G1 Grde 7/8 Mth Circles Geometric Arithmetic Octoer 31, 2012 Centre for Eduction in Mthemtics nd Computing Ancient Greece hs given irth to some of the most importnt

More information

Text mining: bag of words representation and beyond it

Text mining: bag of words representation and beyond it Text mining: bg of words representtion nd beyond it Jsmink Dobš Fculty of Orgniztion nd Informtics University of Zgreb 1 Outline Definition of text mining Vector spce model or Bg of words representtion

More information

Accuracy in Warping ECT LV Surfaces to CT Angiography Coronary Vessels

Accuracy in Warping ECT LV Surfaces to CT Angiography Coronary Vessels Accurcy in Wrping ECT LV Surfces to CT Angiogrphy Coronry Vessels BC Lee, JN Kritzmn, JR Corbett, EP Ficro * University of Michign Helth System, Ann Arbor, MI Disclosure: * Receive softwre roylties from

More information

AML710 CAD LECTURE 16 SURFACES. 1. Analytical Surfaces 2. Synthetic Surfaces

AML710 CAD LECTURE 16 SURFACES. 1. Analytical Surfaces 2. Synthetic Surfaces AML7 CAD LECTURE 6 SURFACES. Anlticl Surfces. Snthetic Surfces Surfce Representtion From CAD/CAM point of view surfces re s importnt s curves nd solids. We need to hve n ide of curves for surfce cretion.

More information

P(r)dr = probability of generating a random number in the interval dr near r. For this probability idea to make sense we must have

P(r)dr = probability of generating a random number in the interval dr near r. For this probability idea to make sense we must have Rndom Numers nd Monte Crlo Methods Rndom Numer Methods The integrtion methods discussed so fr ll re sed upon mking polynomil pproximtions to the integrnd. Another clss of numericl methods relies upon using

More information

Lecture 7: Building 3D Models (Part 1) Prof Emmanuel Agu. Computer Science Dept. Worcester Polytechnic Institute (WPI)

Lecture 7: Building 3D Models (Part 1) Prof Emmanuel Agu. Computer Science Dept. Worcester Polytechnic Institute (WPI) Computer Grphics (CS 4731) Lecture 7: Building 3D Models (Prt 1) Prof Emmnuel Agu Computer Science Dept. Worcester Polytechnic Institute (WPI) Stndrd d Unit itvectors Define y i j 1,0,0 0,1,0 k i k 0,0,1

More information

Optics and Optical design Problems

Optics and Optical design Problems Optics nd Opticl design 0 Problems Sven-Görn Pettersson / Cord Arnold 0-09-06 4:3 This mteril is tken from severl sources. Some problems re from the book Våglär och Optik by Görn Jönsson nd Elisbeth Nilsson.

More information

From Indexing Data Structures to de Bruijn Graphs

From Indexing Data Structures to de Bruijn Graphs From Indexing Dt Structures to de Bruijn Grphs Bstien Czux, Thierry Lecroq, Eric Rivls LIRMM & IBC, Montpellier - LITIS Rouen June 1, 201 Czux, Lecroq, Rivls (LIRMM) Generlized Suffix Tree & DBG June 1,

More information

ONU Calculus I Math 1631

ONU Calculus I Math 1631 ONU Clculus I Mth 1631 2013-2014 Syllus Mrs. Trudy Thompson tthompson@lcchs.edu Text: Clculus 8 th Edition, Anton, Bivens nd Dvis Prerequisites: C or etter in Pre-Clc nd techer s permission This course

More information

this grammar generates the following language: Because this symbol will also be used in a later step, it receives the

this grammar generates the following language: Because this symbol will also be used in a later step, it receives the LR() nlysis Drwcks of LR(). Look-hed symols s eplined efore, concerning LR(), it is possile to consult the net set to determine, in the reduction sttes, for which symols it would e possile to perform reductions.

More information

4-1 NAME DATE PERIOD. Study Guide. Parallel Lines and Planes P Q, O Q. Sample answers: A J, A F, and D E

4-1 NAME DATE PERIOD. Study Guide. Parallel Lines and Planes P Q, O Q. Sample answers: A J, A F, and D E 4-1 NAME DATE PERIOD Pges 142 147 Prllel Lines nd Plnes When plnes do not intersect, they re sid to e prllel. Also, when lines in the sme plne do not intersect, they re prllel. But when lines re not in

More information

Tool Vendor Perspectives SysML Thus Far

Tool Vendor Perspectives SysML Thus Far Frontiers 2008 Pnel Georgi Tec, 05-13-08 Tool Vendor Perspectives SysML Thus Fr Hns-Peter Hoffmnn, Ph.D Chief Systems Methodologist Telelogic, Systems & Softwre Modeling Business Unit Peter.Hoffmnn@telelogic.com

More information

A PDE based approach to multi-domain partitioning and quadrilateral meshing

A PDE based approach to multi-domain partitioning and quadrilateral meshing A PDE sed pproch to multi-domin prtitioning nd qudrilterl meshing Nicols Kowlski 1, Frnck Ledoux 1, nd Pscl Frey 2 1 CEA, DAM, DIF, F-91297 Arpjon, Frnce kowlskin@ocre.ce.fr, frnck.ledoux@ce.fr 2 UPMC

More information

)

) Chpter Five /SOLUTIONS Since the speed ws between nd mph during this five minute period, the fuel efficienc during this period is between 5 mpg nd 8 mpg. So the fuel used during this period is between

More information

Lecture 5: Spatial Analysis Algorithms

Lecture 5: Spatial Analysis Algorithms Lecture 5: Sptil Algorithms GEOG 49: Advnced GIS Sptil Anlsis Algorithms Bsis of much of GIS nlsis tod Mnipultion of mp coordintes Bsed on Eucliden coordinte geometr http://stronom.swin.edu.u/~pbourke/geometr/

More information

Area & Volume. Chapter 6.1 & 6.2 September 25, y = 1! x 2. Back to Area:

Area & Volume. Chapter 6.1 & 6.2 September 25, y = 1! x 2. Back to Area: Bck to Are: Are & Volume Chpter 6. & 6. Septemer 5, 6 We cn clculte the re etween the x-xis nd continuous function f on the intervl [,] using the definite integrl:! f x = lim$ f x * i )%x n i= Where fx

More information

1 Quad-Edge Construction Operators

1 Quad-Edge Construction Operators CS48: Computer Grphics Hndout # Geometric Modeling Originl Hndout #5 Stnford University Tuesdy, 8 December 99 Originl Lecture #5: 9 November 99 Topics: Mnipultions with Qud-Edge Dt Structures Scribe: Mike

More information

Conic Sections Parabola Objective: Define conic section, parabola, draw a parabola, standard equations and their graphs

Conic Sections Parabola Objective: Define conic section, parabola, draw a parabola, standard equations and their graphs Conic Sections Prol Ojective: Define conic section, prol, drw prol, stndrd equtions nd their grphs The curves creted y intersecting doule npped right circulr cone with plne re clled conic sections. If

More information

Solving Problems by Searching. CS 486/686: Introduction to Artificial Intelligence

Solving Problems by Searching. CS 486/686: Introduction to Artificial Intelligence Solving Prolems y Serching CS 486/686: Introduction to Artificil Intelligence 1 Introduction Serch ws one of the first topics studied in AI - Newell nd Simon (1961) Generl Prolem Solver Centrl component

More information

Pythagoras theorem and trigonometry (2)

Pythagoras theorem and trigonometry (2) HPTR 10 Pythgors theorem nd trigonometry (2) 31 HPTR Liner equtions In hpter 19, Pythgors theorem nd trigonometry were used to find the lengths of sides nd the sizes of ngles in right-ngled tringles. These

More information

Geometry Subsystem Design

Geometry Subsystem Design Geometr Susstem Design Ln-D Vn ( 范倫達 ), Ph. D. Deprtment of Computer Science Ntionl Chio Tung Universit Hisnchu, Tiwn Fll, 206 206/0/4 Outline Geometr Susstem Introduction to Shding Algorithms Proposed

More information

Fig.25: the Role of LEX

Fig.25: the Role of LEX The Lnguge for Specifying Lexicl Anlyzer We shll now study how to uild lexicl nlyzer from specifiction of tokens in the form of list of regulr expressions The discussion centers round the design of n existing

More information

Approximation of Two-Dimensional Rectangle Packing

Approximation of Two-Dimensional Rectangle Packing pproximtion of Two-imensionl Rectngle Pcking Pinhong hen, Yn hen, Mudit Goel, Freddy Mng S70 Project Report, Spring 1999. My 18, 1999 1 Introduction 1-d in pcking nd -d in pcking re clssic NP-complete

More information

LETKF compared to 4DVAR for assimilation of surface pressure observations in IFS

LETKF compared to 4DVAR for assimilation of surface pressure observations in IFS LETKF compred to 4DVAR for ssimiltion of surfce pressure oservtions in IFS Pu Escrià, Mssimo Bonvit, Mts Hmrud, Lrs Isksen nd Pul Poli Interntionl Conference on Ensemle Methods in Geophysicl Sciences Toulouse,

More information

What are suffix trees?

What are suffix trees? Suffix Trees 1 Wht re suffix trees? Allow lgorithm designers to store very lrge mount of informtion out strings while still keeping within liner spce Allow users to serch for new strings in the originl

More information

9.1 apply the distance and midpoint formulas

9.1 apply the distance and midpoint formulas 9.1 pply the distnce nd midpoint formuls DISTANCE FORMULA MIDPOINT FORMULA To find the midpoint between two points x, y nd x y 1 1,, we Exmple 1: Find the distnce between the two points. Then, find the

More information

Solving Problems by Searching. CS 486/686: Introduction to Artificial Intelligence Winter 2016

Solving Problems by Searching. CS 486/686: Introduction to Artificial Intelligence Winter 2016 Solving Prolems y Serching CS 486/686: Introduction to Artificil Intelligence Winter 2016 1 Introduction Serch ws one of the first topics studied in AI - Newell nd Simon (1961) Generl Prolem Solver Centrl

More information

Section 10.4 Hyperbolas

Section 10.4 Hyperbolas 66 Section 10.4 Hyperbols Objective : Definition of hyperbol & hyperbols centered t (0, 0). The third type of conic we will study is the hyperbol. It is defined in the sme mnner tht we defined the prbol

More information

Tilt-Sensing with Kionix MEMS Accelerometers

Tilt-Sensing with Kionix MEMS Accelerometers Tilt-Sensing with Kionix MEMS Accelerometers Introduction Tilt/Inclintion sensing is common ppliction for low-g ccelerometers. This ppliction note describes how to use Kionix MEMS low-g ccelerometers to

More information

Solutions to Math 41 Final Exam December 12, 2011

Solutions to Math 41 Final Exam December 12, 2011 Solutions to Mth Finl Em December,. ( points) Find ech of the following its, with justifiction. If there is n infinite it, then eplin whether it is or. ( ) / ln() () (5 points) First we compute the it:

More information

6.2 Volumes of Revolution: The Disk Method

6.2 Volumes of Revolution: The Disk Method mth ppliction: volumes by disks: volume prt ii 6 6 Volumes of Revolution: The Disk Method One of the simplest pplictions of integrtion (Theorem 6) nd the ccumultion process is to determine so-clled volumes

More information

Illumination and Shading

Illumination and Shading Illumintion nd hding In order to produce relistic imges, we must simulte the ppernce of surfces under vrious lighting conditions. Illumintion models: given the illumintion incident t point on surfce, wht

More information

MATH 2530: WORKSHEET 7. x 2 y dz dy dx =

MATH 2530: WORKSHEET 7. x 2 y dz dy dx = MATH 253: WORKSHT 7 () Wrm-up: () Review: polr coordintes, integrls involving polr coordintes, triple Riemnn sums, triple integrls, the pplictions of triple integrls (especilly to volume), nd cylindricl

More information

Computing offsets of freeform curves using quadratic trigonometric splines

Computing offsets of freeform curves using quadratic trigonometric splines Computing offsets of freeform curves using qudrtic trigonometric splines JIULONG GU, JAE-DEUK YUN, YOONG-HO JUNG*, TAE-GYEONG KIM,JEONG-WOON LEE, BONG-JUN KIM School of Mechnicl Engineering Pusn Ntionl

More information

50 AMC LECTURES Lecture 2 Analytic Geometry Distance and Lines. can be calculated by the following formula:

50 AMC LECTURES Lecture 2 Analytic Geometry Distance and Lines. can be calculated by the following formula: 5 AMC LECTURES Lecture Anlytic Geometry Distnce nd Lines BASIC KNOWLEDGE. Distnce formul The distnce (d) between two points P ( x, y) nd P ( x, y) cn be clculted by the following formul: d ( x y () x )

More information

Unit #9 : Definite Integral Properties, Fundamental Theorem of Calculus

Unit #9 : Definite Integral Properties, Fundamental Theorem of Calculus Unit #9 : Definite Integrl Properties, Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl

More information

such that the S i cover S, or equivalently S

such that the S i cover S, or equivalently S MATH 55 Triple Integrls Fll 16 1. Definition Given solid in spce, prtition of consists of finite set of solis = { 1,, n } such tht the i cover, or equivlently n i. Furthermore, for ech i, intersects i

More information

Integration. October 25, 2016

Integration. October 25, 2016 Integrtion October 5, 6 Introduction We hve lerned in previous chpter on how to do the differentition. It is conventionl in mthemtics tht we re supposed to lern bout the integrtion s well. As you my hve

More information

15. 3D-Reconstruction from Vanishing Points

15. 3D-Reconstruction from Vanishing Points 15. 3D-Reconstruction from Vnishing Points Christin B.U. Perwss 1 nd Jon Lseny 2 1 Cvendish Lortory, Cmridge 2 C. U. Engineering Deprtment, Cmridge 15.1 Introduction 3D-reconstruction is currently n ctive

More information

Introduction Transformation formulae Polar graphs Standard curves Polar equations Test GRAPHS INU0114/514 (MATHS 1)

Introduction Transformation formulae Polar graphs Standard curves Polar equations Test GRAPHS INU0114/514 (MATHS 1) POLAR EQUATIONS AND GRAPHS GEOMETRY INU4/54 (MATHS ) Dr Adrin Jnnett MIMA CMth FRAS Polr equtions nd grphs / 6 Adrin Jnnett Objectives The purpose of this presenttion is to cover the following topics:

More information

Digital Design. Chapter 6: Optimizations and Tradeoffs

Digital Design. Chapter 6: Optimizations and Tradeoffs Digitl Design Chpter 6: Optimiztions nd Trdeoffs Slides to ccompny the tetbook Digitl Design, with RTL Design, VHDL, nd Verilog, 2nd Edition, by Frnk Vhid, John Wiley nd Sons Publishers, 2. http://www.ddvhid.com

More information

Loop Pipelining in Hardware-Software Partitioning

Loop Pipelining in Hardware-Software Partitioning Loop Pipelining in Hrdwre-Softwre Prtitioning Jinhwn Jeon nd Kioung Choi School of Electricl Engineering Seoul Ntionl Universit Seoul, Kore 151-742 Tel: +82-2-880-5457 F: +82-2-887-6575 e-mil: {jeonjinh,kchoi}@popp.snu.c.kr

More information

CHAPTER III IMAGE DEWARPING (CALIBRATION) PROCEDURE

CHAPTER III IMAGE DEWARPING (CALIBRATION) PROCEDURE CHAPTER III IMAGE DEWARPING (CALIBRATION) PROCEDURE 3.1 Scheimpflug Configurtion nd Perspective Distortion Scheimpflug criterion were found out to be the best lyout configurtion for Stereoscopic PIV, becuse

More information

Compiler-Assisted Cache Replacement

Compiler-Assisted Cache Replacement LCPC 3 Formulting The Prolem of Compiler-Assisted Cche Replcement Hongo Yng LCPC 3 Agend Bckground: Memory hierrchy, ISA with cche hints Prolem definition: How should compiler give cche hint to minimize

More information

Chapter Spline Method of Interpolation More Examples Electrical Engineering

Chapter Spline Method of Interpolation More Examples Electrical Engineering Chpter. Spline Method of Interpoltion More Exmples Electricl Engineering Exmple Thermistors re used to mesure the temperture of bodies. Thermistors re bsed on mterils chnge in resistnce with temperture.

More information

Qubit allocation for quantum circuit compilers

Qubit allocation for quantum circuit compilers Quit lloction for quntum circuit compilers Nov. 10, 2017 JIQ 2017 Mrcos Yukio Sirichi Sylvin Collnge Vinícius Fernndes dos Sntos Fernndo Mgno Quintão Pereir Compilers for quntum computing The first genertion

More information