Computational Methods for Domains with! Complex Boundaries-I!

Size: px
Start display at page:

Download "Computational Methods for Domains with! Complex Boundaries-I!"

Transcription

1 Computational Methods or Domains with Comple Boundaries-I Grétar Trggvason Spring For most engineering problems it is necessar to deal with comple geometries, consisting o arbitraril curved and oriented boundaries Outline Grid Generation How to deal with irregular domains Overview o various strategies Boundar-itted coordinates - Navier-Stokes equations in vorticit orm - Navier-Stokes equations in primitive orm Grid generation or bod-itted coordinates - Algebraic methods - Dierential methods From: Overview Various Strategies or Comple Geometries and to concentrate grid points in speciic regions Staircasing Originall, rectangular staircasing on rectangular Cartesian grids was used to represent comple boundaries This was ollowed b bod itted grids, the grids are still structured but grid lines are not straight. In commercial codes, unstructured grids have now mostl replaced bod itted grids Regular structured grids are still o interest, particularl when coupled with Adaptive Mesh Reinement (AMR and immersed boundaries Approimate a curved boundar b a the nearest grid lines

2 Unstructured Grids Unstructured versus structured grids structured grids: an ordered laout o grid points. unstructured grids: an arbitrar laout o grid points. Inormation about the laout must be provided Overview Finite Dierence methods on bod itted grids are generall derived b mapping the equations Finite Volume methods are generall derived b using arbitraril shaped control volumes. Boundar-Fitted Coordinates Boundar-Fitted Coordinates Coordinate mapping: transorm the domain into a simpler (usuall rectangular domain. Boundaries are aligned with a constant coordinate line, thus simpliing the treatment o boundar conditions The mathematical equations become more complicated Boundar-Fitted Coordinates (, (, (, (, (, (, First consider the D case: ( (( ( d d d d d d d d d d Boundar-Fitted Coordinates For the irst derivative the change o variables is straightorward using the chain rule: For the second derivative the derivation becomes considerabl more comple:

3 Boundar-Fitted Coordinates The second derivative is given b d d d d d d d d d d d d d d d d d d d d Where we have used the epression or the ist derivative or the inal step. However, since the equations will be discretized in the new grid sstem, it is important to end up with terms like /, not /. d d d d d d d d d d d d d d + d d d d d d d d d d d d + d d d d d d d d d d + d d d d Boundar-Fitted Coordinates To do so, we look at the second derivative in the new sstem d d d d d d d d d d d d d d d d d d d d d d d d + d d d d d d d d d d d d + d d d d d d d d + d d d d d d Solving or the original derivative (which is the one we need to transorm we get: d d d d d d d d d d d d d d d d d d B the chain rule we have d d d d d d Boundar-Fitted Coordinates Oten we need the derivatives o the transormation itsel: d d d d For the second derivative we dierentiate the above: d d d d d d d d d d + d d d d d d d d d d + d d d d Giving: Boundar-Fitted Coordinates D: First Derivatives Change o variables (, ( (,, (, (, The equations will be discretized in the new grid sstem (,. Thereore, it is important to end up with terms like /, not /. d d d d d d d d d d d d d d d d d d Boundar-Fitted Coordinates We want to derive epressions or /, / in the mapped coordinate sstem. Using the chain rule, as we did or the D case: + + Boundar-Fitted Coordinates Solving or the derivatives + + Subtracting + +

4 ,, Solving or the original derivatives ields: where is the acobian. Boundar-Fitted Coordinates A short-hand notation: Boundar-Fitted Coordinates ( ( Rewriting in short-hand notation where is the acobian. Boundar-Fitted Coordinates These relations can also be written in conservative orm: Boundar-Fitted Coordinates Since: + And similarl or the other equation The second derivatives is ound b repeated application o the rules or the irst derivative Similarl D: Second Derivatives Boundar-Fitted Coordinates Adding and ields an epression or the aplacian: + Boundar-Fitted Coordinates

5 + Boundar-Fitted Coordinates ( ( ( ( + where Boundar-Fitted Coordinates + q q q + q + q + + q q + q + q + q q Boundar-Fitted Coordinates Epanding the derivatives ields ( q q + q + ( + ( where [ ( + q + q ] [ ( + + q q ] + + Derivation o i hence similarl Boundar-Fitted Coordinates + ( ( Boundar-Fitted Coordinates Putting them together, it can be shown that (prove it [ q ( + q ( ] q [ q ( + q ( ] q Boundar-Fitted Coordinates We also have, or an unction and g g g ( g g

6 Boundar-Fitted Coordinates A comple domain can be mapped into a rectangular domain where all grid lines are straight. The equations must, however, be rewritten in the new domain. (, (, (, (, Thus: ( ( And more comple epressions or the higher derivatives Vorticit-Stream Function Formulation Vorticit-Stream Function Formulation Vorticit-Stream Function Formulation The Navier-Stokes equations in vorticit orm are: ω t + ψ ω ψ ω ν ω ψ ω Using the transormation relations obtained earlier, The Navier-Stokes equations in vorticit orm become: ω t + (ψ ω ψ ω ν ( q ω q ω + q ω + ( ω + ( ω ( q ψ q ψ + q ψ + ( ψ + ( ψ ω q + q + q + Boundar Conditions Inlow Vorticit-Stream Function Formulation Δ Δ : Inlow N : Outlow, M : No slip ψ Q ψ Q ud vd Outlow ( dψ ψ ψ ower wall Vorticit-Stream Function Formulation Stream unction: ψ Vorticit: (no-slip ψ (, ψ (, + ψ (, + ψ (, + HOT Using that + ω(, ψ (, We have: ω(, + [ ψ(, ψ(, ]

7 Upper wall M Vorticit-Stream Function Formulation Stream unction: ψ Q Vorticit: M Q (ud vd dψ ψ M ψ ω(,m + [ ψ(,m ψ(,m ] Inlet low Vorticit-Stream Function Formulation ( Considering a ull-developed parabolic proile u(, C( 6Q + d C C Q C( 6 6Q u u (, ( ω 6Q Q( ( d Q Inlet low Vorticit-Stream Function Formulation ( Considering a ull-developed parabolic proile and assume that M 6Q u(, M M ψ, Q M M Q ω M (, 6 Outlow Vorticit-Stream Function Formulation ( N Tpicall, assuming straight streamlines ψ n I is normal to the outlow boundar, this ields ψ ψ I not, then a proper transormation is needed or n

Boundary-Fitted Coordinates!

Boundary-Fitted Coordinates! Computational Fluid Dnamics http:wwwndedu~gtrggvacfdcourse Computational Fluid Dnamics Computational Methods or Domains with Comple BoundariesI Grétar Trggvason Spring For most engineering problems it

More information

Topic02_PDE. 8/29/2006 topic02_pde 1. Computational Fluid Dynamics (AE/ME 339) MAE Dept., UMR

Topic02_PDE. 8/29/2006 topic02_pde 1. Computational Fluid Dynamics (AE/ME 339) MAE Dept., UMR MEAE 9 Computational Fluid Dnamics Topic0_ 89006 topic0_ Partial Dierential Equations () (CLW: 7., 7., 7.4) s can be linear or nonlinear Order : Determined b the order o the highest derivative. Linear,

More information

AE/ME 339. K. M. Isaac Professor of Aerospace Engineering. December 21, 2001 topic13_grid_generation 1

AE/ME 339. K. M. Isaac Professor of Aerospace Engineering. December 21, 2001 topic13_grid_generation 1 AE/ME 339 Professor of Aerospace Engineering December 21, 2001 topic13_grid_generation 1 The basic idea behind grid generation is the creation of the transformation laws between the phsical space and the

More information

SIMULATION OF FLOW OF FIBRE SUSPENSIONS USING AN IRBF-BCF BASED MULTISCALE APPROACH

SIMULATION OF FLOW OF FIBRE SUSPENSIONS USING AN IRBF-BCF BASED MULTISCALE APPROACH Eleventh International Conerence on CFD in the Minerals and Process Industries CSIRO Melbourne Australia 7-9 December 25 SIMULATION OF FLOW OF FIBRE SUSPENSIONS USING AN IRBF-BCF BASED MULTISCALE APPROACH

More information

RATIONAL FUNCTIONS. Finding Asymptotes..347 The Domain Finding Intercepts Graphing Rational Functions

RATIONAL FUNCTIONS. Finding Asymptotes..347 The Domain Finding Intercepts Graphing Rational Functions RATIONAL FUNCTIONS Finding Asymptotes..347 The Domain....350 Finding Intercepts.....35 Graphing Rational Functions... 35 345 Objectives The ollowing is a list o objectives or this section o the workbook.

More information

9.1 The Square Root Function

9.1 The Square Root Function Section 9.1 The Square Root Function 869 9.1 The Square Root Function In this section we turn our attention to the square root unction, the unction deined b the equation () =. (1) We begin the section

More information

. This is the Basic Chain Rule. x dt y dt z dt Chain Rule in this context.

. This is the Basic Chain Rule. x dt y dt z dt Chain Rule in this context. Math 18.0A Gradients, Chain Rule, Implicit Dierentiation, igher Order Derivatives These notes ocus on our things: (a) the application o gradients to ind normal vectors to curves suraces; (b) the generaliation

More information

An Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains

An Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains An Immersed Interace Method or the Incompressible Navier-Stokes Equations in Irregular Domains Duc-Vinh Le, Boo Cheong Khoo, Jaime Peraire Singapore-MIT Alliance Department o Mechanical Engineering, National

More information

y2 = 0. Show that u = e2xsin(2y) satisfies Laplace's equation.

y2 = 0. Show that u = e2xsin(2y) satisfies Laplace's equation. Review 1 1) State the largest possible domain o deinition or the unction (, ) = 3 - ) Determine the largest set o points in the -plane on which (, ) = sin-1( - ) deines a continuous unction 3) Find the

More information

A Fourier Transform Model in Excel #1

A Fourier Transform Model in Excel #1 A Fourier Transorm Model in Ecel # -This is a tutorial about the implementation o a Fourier transorm in Ecel. This irst part goes over adjustments in the general Fourier transorm ormula to be applicable

More information

Controlling the Heat Flux Distribution by Changing the Thickness of Heated Wall

Controlling the Heat Flux Distribution by Changing the Thickness of Heated Wall J. Basic. Appl. Sci. Res., 2(7)7270-7275, 2012 2012, TextRoad Publication ISSN 2090-4304 Journal o Basic and Applied Scientiic Research www.textroad.com Controlling the Heat Flux Distribution by Changing

More information

Review of Prerequisite Skills for Unit # 2 (Derivatives) U2L2: Sec.2.1 The Derivative Function

Review of Prerequisite Skills for Unit # 2 (Derivatives) U2L2: Sec.2.1 The Derivative Function UL1: Review o Prerequisite Skills or Unit # (Derivatives) Working with the properties o exponents Simpliying radical expressions Finding the slopes o parallel and perpendicular lines Simpliying rational

More information

Differential Equations

Differential Equations LOCUS Dierential Equations CONCEPT NOTES 0. Motiation 0. Soling Dierential Equations LOCUS Dierential Equations Section - MOTIVATION A dierential equation can simpl be said to be an equation inoling deriaties

More information

Mathematical Preliminaries. Developed for the Members of Azera Global By: Joseph D. Fournier B.Sc.E.E., M.Sc.E.E.

Mathematical Preliminaries. Developed for the Members of Azera Global By: Joseph D. Fournier B.Sc.E.E., M.Sc.E.E. Mathematical Preliminaries Developed or the Members o Azera Global B: Joseph D. Fournier B.Sc.E.E., M.Sc.E.E. Outline Chapter One, Sets: Slides: 3-27 Chapter Two, Introduction to unctions: Slides: 28-36

More information

Lecture Outline. Basics of Spatial Filtering Smoothing Spatial Filters. Sharpening Spatial Filters

Lecture Outline. Basics of Spatial Filtering Smoothing Spatial Filters. Sharpening Spatial Filters 1 Lecture Outline Basics o Spatial Filtering Smoothing Spatial Filters Averaging ilters Order-Statistics ilters Sharpening Spatial Filters Laplacian ilters High-boost ilters Gradient Masks Combining Spatial

More information

8. THEOREM If the partial derivatives f x. and f y exist near (a, b) and are continuous at (a, b), then f is differentiable at (a, b).

8. THEOREM If the partial derivatives f x. and f y exist near (a, b) and are continuous at (a, b), then f is differentiable at (a, b). 8. THEOREM I the partial derivatives and eist near (a b) and are continuous at (a b) then is dierentiable at (a b). For a dierentiable unction o two variables z= ( ) we deine the dierentials d and d to

More information

Curve Sketching. The process of curve sketching can be performed in the following steps:

Curve Sketching. The process of curve sketching can be performed in the following steps: Curve Sketching So ar you have learned how to ind st and nd derivatives o unctions and use these derivatives to determine where a unction is:. Increasing/decreasing. Relative extrema 3. Concavity 4. Points

More information

3. Several Random Variables

3. Several Random Variables . Several Random Variables. Two Random Variables. Conditional Probabilit--Revisited. Statistical Independence.4 Correlation between Random Variables. Densit unction o the Sum o Two Random Variables. Probabilit

More information

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is:

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is: Noll 6 CT Notes : Page Compute Tomograph Notes Part Challenges with Projection X-ra Sstems The equation that governs the image intensit in projection imaging is: z I I ep μ z Projection -ra sstems are

More information

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is:

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is: Noll 3 CT Notes : Page Compute Tomograph Notes Part Challenges with Projection X-ra Sstems The equation that governs the image intensit in projection imaging is: z I I ep µ z Projection -ra sstems are

More information

New Functions from Old Functions

New Functions from Old Functions .3 New Functions rom Old Functions In this section we start with the basic unctions we discussed in Section. and obtain new unctions b shiting, stretching, and relecting their graphs. We also show how

More information

CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a

CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a Chapter 4 CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS Convective heat transer analysis o nanoluid lowing inside a straight tube o circular cross-section under laminar and turbulent conditions

More information

Lecture : Feedback Linearization

Lecture : Feedback Linearization ecture : Feedbac inearization Niola Misovic, dipl ing and Pro Zoran Vuic June 29 Summary: This document ollows the lectures on eedbac linearization tought at the University o Zagreb, Faculty o Electrical

More information

Particular Solutions

Particular Solutions Particular Solutions Our eamples so far in this section have involved some constant of integration, K. We now move on to see particular solutions, where we know some boundar conditions and we substitute

More information

Categories and Natural Transformations

Categories and Natural Transformations Categories and Natural Transormations Ethan Jerzak 17 August 2007 1 Introduction The motivation or studying Category Theory is to ormalise the underlying similarities between a broad range o mathematical

More information

Physics 5153 Classical Mechanics. Solution by Quadrature-1

Physics 5153 Classical Mechanics. Solution by Quadrature-1 October 14, 003 11:47:49 1 Introduction Physics 5153 Classical Mechanics Solution by Quadrature In the previous lectures, we have reduced the number o eective degrees o reedom that are needed to solve

More information

Lecture 8 Optimization

Lecture 8 Optimization 4/9/015 Lecture 8 Optimization EE 4386/5301 Computational Methods in EE Spring 015 Optimization 1 Outline Introduction 1D Optimization Parabolic interpolation Golden section search Newton s method Multidimensional

More information

Solving Partial Differential Equations Numerically. Miklós Bergou with: Gary Miller, David Cardoze, Todd Phillips, Mark Olah

Solving Partial Differential Equations Numerically. Miklós Bergou with: Gary Miller, David Cardoze, Todd Phillips, Mark Olah Solving Partial Dierential Equations Numerically Miklós Bergou with: Gary Miller, David Cardoze, Todd Phillips, Mark Olah Overview What are partial dierential equations? How do we solve them? (Example)

More information

Chapter 7: Boundary Layer Theory

Chapter 7: Boundary Layer Theory Proessor Fred Stern Fall 17 1 Chapter 7: Boundar Laer Theor 7.1. Introduction: Boundar laer lows: Eternal lows around streamlined bodies at high Re have viscous (shear and no-slip) eects conined close

More information

Differentiation. The main problem of differential calculus deals with finding the slope of the tangent line at a point on a curve.

Differentiation. The main problem of differential calculus deals with finding the slope of the tangent line at a point on a curve. Dierentiation The main problem o dierential calculus deals with inding the slope o the tangent line at a point on a curve. deinition() : The slope o a curve at a point p is the slope, i it eists, o the

More information

! " k x 2k$1 # $ k x 2k. " # p $ 1! px! p " p 1 # !"#$%&'"()'*"+$",&-('./&-/. !"#$%&'()"*#%+!'",' -./#")'.,&'+.0#.1)2,'!%)2%! !"#$%&'"%(")*$+&#,*$,#

!  k x 2k$1 # $ k x 2k.  # p $ 1! px! p  p 1 # !#$%&'()'*+$,&-('./&-/. !#$%&'()*#%+!',' -./#)'.,&'+.0#.1)2,'!%)2%! !#$%&'%()*$+&#,*$,# "#$%&'()"*#%+'",' -./#")'.,&'+.0#.1)2,' %)2% "#$%&'"()'*"+$",&-('./&-/. Taylor Series o a unction at x a is " # a k " # " x a# k k0 k It is a Power Series centered at a. Maclaurin Series o a unction is

More information

CH.7. PLANE LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics

CH.7. PLANE LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics CH.7. PLANE LINEAR ELASTICITY Multimedia Course on Continuum Mechanics Overview Plane Linear Elasticit Theor Plane Stress Simplifing Hpothesis Strain Field Constitutive Equation Displacement Field The

More information

Estimation of the Vehicle Sideslip Angle by Means of the State Dependent Riccati Equation Technique

Estimation of the Vehicle Sideslip Angle by Means of the State Dependent Riccati Equation Technique Proceedings o the World Congress on Engineering 7 Vol II WCE 7, Jul 5-7, 7, London, U.K. Estimation o the Vehicle Sideslip Angle b Means o the State Dependent Riccati Equation Technique S. Strano, M. Terzo

More information

Syllabus Objective: 2.9 The student will sketch the graph of a polynomial, radical, or rational function.

Syllabus Objective: 2.9 The student will sketch the graph of a polynomial, radical, or rational function. Precalculus Notes: Unit Polynomial Functions Syllabus Objective:.9 The student will sketch the graph o a polynomial, radical, or rational unction. Polynomial Function: a unction that can be written in

More information

A Wavelet Collocation Method for Optimal Control. Problems

A Wavelet Collocation Method for Optimal Control. Problems A Wavelet ollocation ethod or Optimal ontrol Problems Ran Dai * and John E. ochran Jr. Abstract A Haar wavelet technique is discussed as a method or discretizing the nonlinear sstem equations or optimal

More information

CHAPTER-III CONVECTION IN A POROUS MEDIUM WITH EFFECT OF MAGNETIC FIELD, VARIABLE FLUID PROPERTIES AND VARYING WALL TEMPERATURE

CHAPTER-III CONVECTION IN A POROUS MEDIUM WITH EFFECT OF MAGNETIC FIELD, VARIABLE FLUID PROPERTIES AND VARYING WALL TEMPERATURE CHAPER-III CONVECION IN A POROUS MEDIUM WIH EFFEC OF MAGNEIC FIELD, VARIABLE FLUID PROPERIES AND VARYING WALL EMPERAURE 3.1. INRODUCION Heat transer studies in porous media ind applications in several

More information

Environmental Process Analysis, 2: Dynamic Steady States

Environmental Process Analysis, 2: Dynamic Steady States Environmental Process Analysis, 2: Dynamic Steady States T. Torgersen Department o Marine Sciences, University o Connecticut, Groton, CT 06340-6097, Thomas.Torgersen@uconn.edu B. Branco Department o Marine

More information

On the Development of Implicit Solvers for Time-Dependent Systems

On the Development of Implicit Solvers for Time-Dependent Systems School o something FACULTY School OF OTHER o Computing On the Development o Implicit Solvers or Time-Dependent Systems Peter Jimack School o Computing, University o Leeds In collaboration with: P.H. Gaskell,

More information

Constantine, Algeria. Received Accepted Keywords: Copper nanoparticles; heat transfer; circular cylinder; steady regime.

Constantine, Algeria. Received Accepted Keywords: Copper nanoparticles; heat transfer; circular cylinder; steady regime. Metallurgical and Materials Engineering Association o Metallurgical Engineers o Serbia AMES Scientiic paper UDC: 669.245 NUMERICAL INVESTIGATION OF FLUID FLOW AND HEAT TRANSFER AROUND A CIRCULAR CYLINDER

More information

Roberto s Notes on Differential Calculus Chapter 8: Graphical analysis Section 1. Extreme points

Roberto s Notes on Differential Calculus Chapter 8: Graphical analysis Section 1. Extreme points Roberto s Notes on Dierential Calculus Chapter 8: Graphical analysis Section 1 Extreme points What you need to know already: How to solve basic algebraic and trigonometric equations. All basic techniques

More information

Increasing and Decreasing Functions and the First Derivative Test. Increasing and Decreasing Functions. Video

Increasing and Decreasing Functions and the First Derivative Test. Increasing and Decreasing Functions. Video SECTION and Decreasing Functions and the First Derivative Test 79 Section and Decreasing Functions and the First Derivative Test Determine intervals on which a unction is increasing or decreasing Appl

More information

Analog Computing Technique

Analog Computing Technique Analog Computing Technique by obert Paz Chapter Programming Principles and Techniques. Analog Computers and Simulation An analog computer can be used to solve various types o problems. It solves them in

More information

APPENDIX 1 ERROR ESTIMATION

APPENDIX 1 ERROR ESTIMATION 1 APPENDIX 1 ERROR ESTIMATION Measurements are always subject to some uncertainties no matter how modern and expensive equipment is used or how careully the measurements are perormed These uncertainties

More information

8.4 Inverse Functions

8.4 Inverse Functions Section 8. Inverse Functions 803 8. Inverse Functions As we saw in the last section, in order to solve application problems involving eponential unctions, we will need to be able to solve eponential equations

More information

Separation of Variables in Cartesian Coordinates

Separation of Variables in Cartesian Coordinates Lecture 9 Separation of Variables in Cartesian Coordinates Phs 3750 Overview and Motivation: Toda we begin a more in-depth loo at the 3D wave euation. We introduce a techniue for finding solutions to partial

More information

Feedback Linearization

Feedback Linearization Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)

More information

Computers & Fluids. Improving lattice Boltzmann simulation of moving particles in a viscous flow using local grid refinement

Computers & Fluids. Improving lattice Boltzmann simulation of moving particles in a viscous flow using local grid refinement Computers & Fluids 00 (204) 2 Computers & Fluids Improving lattice Boltzmann simulation o moving particles in a viscous low using local grid reinement Songying Chen a,, Cheng Peng b, Yihua Teng c, Lian-Ping

More information

Functions of Several Variables

Functions of Several Variables Chapter 1 Functions of Several Variables 1.1 Introduction A real valued function of n variables is a function f : R, where the domain is a subset of R n. So: for each ( 1,,..., n ) in, the value of f is

More information

(C) The rationals and the reals as linearly ordered sets. Contents. 1 The characterizing results

(C) The rationals and the reals as linearly ordered sets. Contents. 1 The characterizing results (C) The rationals and the reals as linearly ordered sets We know that both Q and R are something special. When we think about about either o these we usually view it as a ield, or at least some kind o

More information

Global Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance

Global Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam) Abstract. A reormulation o the planetary

More information

Physics 2B Chapter 17 Notes - First Law of Thermo Spring 2018

Physics 2B Chapter 17 Notes - First Law of Thermo Spring 2018 Internal Energy o a Gas Work Done by a Gas Special Processes The First Law o Thermodynamics p Diagrams The First Law o Thermodynamics is all about the energy o a gas: how much energy does the gas possess,

More information

( x) f = where P and Q are polynomials.

( x) f = where P and Q are polynomials. 9.8 Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm ( ) ( ) ( ) P where P and Q are polynomials. Q An eample o a simple rational

More information

Available online at ScienceDirect. Energy Procedia 83 (2015 ) Václav Dvo ák a *, Tomáš Vít a

Available online at   ScienceDirect. Energy Procedia 83 (2015 ) Václav Dvo ák a *, Tomáš Vít a Available online at www.sciencedirect.com ScienceDirect Energy Procedia 83 (205 ) 34 349 7th International Conerence on Sustainability in Energy and Buildings Numerical investigation o counter low plate

More information

Math Review and Lessons in Calculus

Math Review and Lessons in Calculus Math Review and Lessons in Calculus Agenda Rules o Eponents Functions Inverses Limits Calculus Rules o Eponents 0 Zero Eponent Rule a * b ab Product Rule * 3 5 a / b a-b Quotient Rule 5 / 3 -a / a Negative

More information

ENERGY ANALYSIS: CLOSED SYSTEM

ENERGY ANALYSIS: CLOSED SYSTEM ENERGY ANALYSIS: CLOSED SYSTEM A closed system can exchange energy with its surroundings through heat and work transer. In other words, work and heat are the orms that energy can be transerred across the

More information

Numerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods

Numerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods Numerical Methods - Lecture 1 Numerical Methods Lecture. Analysis o errors in numerical methods Numerical Methods - Lecture Why represent numbers in loating point ormat? Eample 1. How a number 56.78 can

More information

1 General description of the OpenFOAM suite

1 General description of the OpenFOAM suite Santiago Márquez Damián-Final Work-Computational Fluid Dynamics Description and utilization o interfoam multiphase solver 1 General description o the OpenFOAM suite OpenFOAM Open Field Operation and Manipulation

More information

MHF 4U Unit 7: Combining Functions May 29, Review Solutions

MHF 4U Unit 7: Combining Functions May 29, Review Solutions MHF 4U Unit 7: Combining Functions May 9, 008. Review Solutions Use the ollowing unctions to answer questions 5, ( ) g( ), h( ) sin, w {(, ), (3, ), (4, 7)}, r, and l ) log ( ) + (, ) Determine: a) + w

More information

6.1 The Linear Elastic Model

6.1 The Linear Elastic Model Linear lasticit The simplest constitutive law or solid materials is the linear elastic law, which assumes a linear relationship between stress and engineering strain. This assumption turns out to be an

More information

A NOVEL METHOD OF INTERPOLATION AND EXTRAPOLATION OF FUNCTIONS BY A LINEAR INITIAL VALUE PROBLEM

A NOVEL METHOD OF INTERPOLATION AND EXTRAPOLATION OF FUNCTIONS BY A LINEAR INITIAL VALUE PROBLEM A OVEL METHOD OF ITERPOLATIO AD EXTRAPOLATIO OF FUCTIOS BY A LIEAR IITIAL VALUE PROBLEM Michael Shatalov Sensor Science and Technolog o CSIR Manuacturing and Materials, P.O.Bo 395, Pretoria, CSIR and Department

More information

Enhancement Using Local Histogram

Enhancement Using Local Histogram Enhancement Using Local Histogram Used to enhance details over small portions o the image. Deine a square or rectangular neighborhood hose center moves rom piel to piel. Compute local histogram based on

More information

The Effect of Internal Obstructions in Naturally Ventilated Greenhouse Applications

The Effect of Internal Obstructions in Naturally Ventilated Greenhouse Applications HEFAT27 5 th International Conerence on Heat Transer, Fluid Mechanics and Thermodynamics 1-4 July 27, Sun City, South Arica Paper number: KS2 The Eect o Internal Obstructions in Naturally Ventilated Greenhouse

More information

Chapter 2. Basic concepts of probability. Summary. 2.1 Axiomatic foundation of probability theory

Chapter 2. Basic concepts of probability. Summary. 2.1 Axiomatic foundation of probability theory Chapter Basic concepts o probability Demetris Koutsoyiannis Department o Water Resources and Environmental Engineering aculty o Civil Engineering, National Technical University o Athens, Greece Summary

More information

0,0 B 5,0 C 0, 4 3,5. y x. Recitation Worksheet 1A. 1. Plot these points in the xy plane: A

0,0 B 5,0 C 0, 4 3,5. y x. Recitation Worksheet 1A. 1. Plot these points in the xy plane: A Math 13 Recitation Worksheet 1A 1 Plot these points in the y plane: A 0,0 B 5,0 C 0, 4 D 3,5 Without using a calculator, sketch a graph o each o these in the y plane: A y B 3 Consider the unction a Evaluate

More information

Periodic Structures in FDTD

Periodic Structures in FDTD EE 5303 Electromagnetic Analsis Using Finite Difference Time Domain Lecture #19 Periodic Structures in FDTD Lecture 19 These notes ma contain coprighted material obtained under fair use rules. Distribution

More information

Chapter 7: Boundary Layer Theory

Chapter 7: Boundary Layer Theory Proessor Fred Stern Fall 14 1 Chapter 7: Boundar Laer Theor 7.1. Introduction: Boundar laer lows: Eternal lows around streamlined bodies at high Re have viscous (shear and no-slip) eects conined close

More information

COMBINED EFFECTS OF RADIATION AND HEAT GENERATION ON MHD NATURAL CONVECTION FLOW ALONG A VERTICAL FLAT PLATE IN PRESENCE OF HEAT CONDUCTION

COMBINED EFFECTS OF RADIATION AND HEAT GENERATION ON MHD NATURAL CONVECTION FLOW ALONG A VERTICAL FLAT PLATE IN PRESENCE OF HEAT CONDUCTION BRAC University Journal, vol.vi, no., 9, pp 11- COMBINED EFFECTS OF RADIATION AND HEAT GENERATION ON MHD NATURAL CONVECTION FLOW ALONG A VERTICAL FLAT PLATE IN PRESENCE OF HEAT CONDUCTION Mohammad Mokaddes

More information

Definition: Let f(x) be a function of one variable with continuous derivatives of all orders at a the point x 0, then the series.

Definition: Let f(x) be a function of one variable with continuous derivatives of all orders at a the point x 0, then the series. 2.4 Local properties o unctions o several variables In this section we will learn how to address three kinds o problems which are o great importance in the ield o applied mathematics: how to obtain the

More information

ROAD MAP... D-1: Aerodynamics of 3-D Wings D-2: Boundary Layer and Viscous Effects D-3: XFLR (Aerodynamics Analysis Tool)

ROAD MAP... D-1: Aerodynamics of 3-D Wings D-2: Boundary Layer and Viscous Effects D-3: XFLR (Aerodynamics Analysis Tool) AE301 Aerodynamics I UNIT D: Applied Aerodynamics ROAD MAP... D-1: Aerodynamics o 3-D Wings D-2: Boundary Layer and Viscous Eects D-3: XFLR (Aerodynamics Analysis Tool) AE301 Aerodynamics I : List o Subjects

More information

Optimal Control. with. Aerospace Applications. James M. Longuski. Jose J. Guzman. John E. Prussing

Optimal Control. with. Aerospace Applications. James M. Longuski. Jose J. Guzman. John E. Prussing Optimal Control with Aerospace Applications by James M. Longuski Jose J. Guzman John E. Prussing Published jointly by Microcosm Press and Springer 2014 Copyright Springer Science+Business Media New York

More information

Numerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective

Numerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective Numerical Solution o Ordinary Dierential Equations in Fluctuationlessness Theorem Perspective NEJLA ALTAY Bahçeşehir University Faculty o Arts and Sciences Beşiktaş, İstanbul TÜRKİYE TURKEY METİN DEMİRALP

More information

Description of a One-Dimensional Numerical Model of an Active Magnetic Regenerator Refrigerator

Description of a One-Dimensional Numerical Model of an Active Magnetic Regenerator Refrigerator This is a 1D model o an active magnetic regenerative rerigerator (AMRR) that was developed in MATLAB. The model uses cycle inputs such as the luid mass low and magnetic ield proiles, luid and regenerator

More information

Objectives. By the time the student is finished with this section of the workbook, he/she should be able

Objectives. By the time the student is finished with this section of the workbook, he/she should be able FUNCTIONS Quadratic Functions......8 Absolute Value Functions.....48 Translations o Functions..57 Radical Functions...61 Eponential Functions...7 Logarithmic Functions......8 Cubic Functions......91 Piece-Wise

More information

Basic properties of limits

Basic properties of limits Roberto s Notes on Dierential Calculus Chapter : Limits and continuity Section Basic properties o its What you need to know already: The basic concepts, notation and terminology related to its. What you

More information

Infinitesimal Rotations

Infinitesimal Rotations Universit of Connecticut DigitalCommons@UConn Chemistr Education Materials Department of Chemistr Januar 007 Infinitesimal Rotations Carl W. David Universit of Connecticut, Carl.David@uconn.edu Follow

More information

Strain Transformation and Rosette Gage Theory

Strain Transformation and Rosette Gage Theory Strain Transformation and Rosette Gage Theor It is often desired to measure the full state of strain on the surface of a part, that is to measure not onl the two etensional strains, and, but also the shear

More information

Time Series Analysis for Quality Improvement: a Soft Computing Approach

Time Series Analysis for Quality Improvement: a Soft Computing Approach ESANN'4 proceedings - European Smposium on Artiicial Neural Networks Bruges (Belgium), 8-3 April 4, d-side publi., ISBN -9337-4-8, pp. 19-114 ime Series Analsis or Qualit Improvement: a Sot Computing Approach

More information

Iterative Methods for Stokes/Darcy Coupling

Iterative Methods for Stokes/Darcy Coupling Iterative Methods or Stokes/Darcy Coupling Marco Discacciati Ecole Polytechnique Fédérale de Lausanne Institut d Analyse et Calcul Scientiique Chair o Modelling and Scientiic Computing marco.discacciati@epl.ch

More information

Chapter 3: Image Enhancement in the. Office room : 841

Chapter 3: Image Enhancement in the.   Office room : 841 Chapter 3: Image Enhancement in the Spatial Domain Lecturer: Jianbing Shen Email : shenjianbing@bit.edu.cn Oice room : 841 http://cs.bit.edu.cn/shenjianbing cn/shenjianbing Principle Objective o Enhancement

More information

The Plane Stress Problem

The Plane Stress Problem . 4 The Plane Stress Problem 4 Chapter 4: THE PLANE STRESS PROBLEM 4 TABLE OF CONTENTS Page 4.. INTRODUCTION 4 3 4... Plate in Plane Stress............... 4 3 4... Mathematical Model.............. 4 4

More information

Global Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance

Global Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam Abstract. A reormulation o the planetary

More information

Probabilistic Optimisation applied to Spacecraft Rendezvous on Keplerian Orbits

Probabilistic Optimisation applied to Spacecraft Rendezvous on Keplerian Orbits Probabilistic Optimisation applied to pacecrat Rendezvous on Keplerian Orbits Grégory aive a, Massimiliano Vasile b a Université de Liège, Faculté des ciences Appliquées, Belgium b Dipartimento di Ingegneria

More information

Last modified Spring 2016

Last modified Spring 2016 Math 00 Final Review Questions In problems 6, perform the indicated operations and simplif if necessar.. 8 6 8. 7 6. ( i) ( 4 i) 4. (8 i). ( 9 i)( 7 i) 6. ( i)( i) In problems 7-, solve the following applications.

More information

y,z the subscript y, z indicating that the variables y and z are kept constant. The second partial differential with respect to x is written x 2 y,z

y,z the subscript y, z indicating that the variables y and z are kept constant. The second partial differential with respect to x is written x 2 y,z 8 Partial dierentials I a unction depends on more than one variable, its rate o change with respect to one o the variables can be determined keeping the others ied The dierential is then a partial dierential

More information

Math 0210 Common Final Review Questions (2 5 i)(2 5 i )

Math 0210 Common Final Review Questions (2 5 i)(2 5 i ) Math 0 Common Final Review Questions In problems 1 6, perform the indicated operations and simplif if necessar. 1. ( 8)(4) ( )(9) 4 7 4 6( ). 18 6 8. ( i) ( 1 4 i ) 4. (8 i ). ( 9 i)( 7 i) 6. ( i)( i )

More information

NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION

NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION by Asterios Pantokratoras School o Engineering, Democritus University o Thrace, 67100

More information

UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX METHOD

UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX METHOD Proceedings of the 3rd ASME/JSME Joint Fluids Engineering Conference Jul 8-23, 999, San Francisco, California FEDSM99-8 UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX

More information

Department of Physics and Astronomy 2 nd Year Laboratory. L2 Light Scattering

Department of Physics and Astronomy 2 nd Year Laboratory. L2 Light Scattering nd ear laborator script L Light Scattering Department o Phsics and Astronom nd Year Laborator L Light Scattering Scientiic aims and objectives To determine the densit o nano-spheres o polstrene suspended

More information

Shai Avidan Tel Aviv University

Shai Avidan Tel Aviv University Image Editing in the Gradient Domain Shai Avidan Tel Aviv Universit Slide Credits (partial list) Rick Szeliski Steve Seitz Alosha Eros Yacov Hel-Or Marc Levo Bill Freeman Fredo Durand Slvain Paris Image

More information

NONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION

NONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION NONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION Steven Ball Science Applications International Corporation Columbia, MD email: sball@nmtedu Steve Schaer Department o Mathematics

More information

CHAPTER 1: INTRODUCTION. 1.1 Inverse Theory: What It Is and What It Does

CHAPTER 1: INTRODUCTION. 1.1 Inverse Theory: What It Is and What It Does Geosciences 567: CHAPTER (RR/GZ) CHAPTER : INTRODUCTION Inverse Theory: What It Is and What It Does Inverse theory, at least as I choose to deine it, is the ine art o estimating model parameters rom data

More information

This is a problem of calculus of variations with an equality constraint. This problem is similar to the problem of minimizing a function

This is a problem of calculus of variations with an equality constraint. This problem is similar to the problem of minimizing a function 3-7 ISOPERIMETRIC PROBLEMS Isoperimetric implies constant perimeter. This is a problem of calculus of variations ith an equalit constraint in the integral form. ( V ) d ; W G( ) d This problem is similar

More information

Basic mathematics of economic models. 3. Maximization

Basic mathematics of economic models. 3. Maximization John Riley 1 January 16 Basic mathematics o economic models 3 Maimization 31 Single variable maimization 1 3 Multi variable maimization 6 33 Concave unctions 9 34 Maimization with non-negativity constraints

More information

Mechanical Engineering Research Journal BUOYANT FLOW OF NANOFLUID FOR HEAT-MASS TRANSFER THROUGH A THIN LAYER

Mechanical Engineering Research Journal BUOYANT FLOW OF NANOFLUID FOR HEAT-MASS TRANSFER THROUGH A THIN LAYER Dept. o Mech. Eng. CUET Published Online March 2015 (http://www.cuet.ac.bd/merj/index.html) Mechanical Engineering Research Journal Vol. 9, pp. 712, 2013 M E R J ISSN: 1990-5491 BUOYANT FLOW OF NANOFLUID

More information

CHAPTER 4 Reactor Statics. Table of Contents

CHAPTER 4 Reactor Statics. Table of Contents CHAPTER 4 Reactor Statics Prepared by Dr. Benjamin Rouben, & Consulting, Adjunct Proessor, McMaster University & University o Ontario Institute o Technology (UOIT) and Dr. Eleodor Nichita, Associate Proessor,

More information

Computational Fluid Dynamics (CFD, CHD)*

Computational Fluid Dynamics (CFD, CHD)* 1 / 1 Computational Fluid Dnamics (CFD, CHD)* PDE (Shocks 1st); Part I: Basics, Part II: Vorticit Fields Rubin H Landau Sall Haerer, Producer-Director Based on A Surve of Computational Phsics b Landau,

More information

OE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004

OE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004 OE465 Vaclav Matousek October 13, 004 1 Dredge Vermelding Pumps onderdeel and Slurry organisatie Transport OE465 Vaclav Matousek October 13, 004 Dredge Vermelding Pumps onderdeel and Slurry organisatie

More information

Tangent Line Approximations

Tangent Line Approximations 60_009.qd //0 :8 PM Page SECTION.9 Dierentials Section.9 EXPLORATION Tangent Line Approimation Use a graphing utilit to graph. In the same viewing window, graph the tangent line to the graph o at the point,.

More information

1036: Probability & Statistics

1036: Probability & Statistics 1036: Probabilit & Statistics Lecture 4 Mathematical pectation Prob. & Stat. Lecture04 - mathematical epectation cwliu@twins.ee.nctu.edu.tw 4-1 Mean o a Random Variable Let be a random variable with probabilit

More information

Fluid Mechanics II. Newton s second law applied to a control volume

Fluid Mechanics II. Newton s second law applied to a control volume Fluid Mechanics II Stead flow momentum equation Newton s second law applied to a control volume Fluids, either in a static or dnamic motion state, impose forces on immersed bodies and confining boundaries.

More information