Finite Element Models for Steady Flows of Viscous Incompressible Fluids

Size: px
Start display at page:

Download "Finite Element Models for Steady Flows of Viscous Incompressible Fluids"

Transcription

1 Finit Elmnt Modls for Stad Flows of Viscous Incomprssibl Fluids Rad: Chaptr 10 JN Rdd CONTENTS Govrning Equations of Flows of Incomprssibl Fluids Mid (Vlocit-Prssur) Finit Elmnt Modl Pnalt Function Mthod - algbraic problm Pnalt Finit Elmnt Modl of Viscous Incomprssibl Fluids Numrical Rsults Closur Viscous Incompr. Flows: 1

2 JN Rdd Govrning Equations of Flows of Viscous incomprssibl Fluids Equations of motion Dv σ + f = ρ Dt σ σ Dv + + f = ρ Dt σ σ Dv + + f = ρ Dt Matrial tim drivativ D D = + v = + v + v Dt t Dt t Consrvation of mass v v = 0 v + = 0 Constitutiv rlations σ = P I + τ, τ = 2µ D Viscous Incompr. Flows: 2

3 Govrning Equations of Flows of Viscous incomprssibl Fluids Kinmatics rlations 1 T D = v + ( v) 2 v v v v D =, D =, 2D = + Strss-vlocit-prssur rlations v v v v σ = 2µ P, σ = 2µ P, σ = µ + Boundar conditions ( v, t ) and ( v, t ) JN Rdd t = σ n + σ n, t = σ n + σ n Viscous Incompr. Flows: 3

4 Govrning Equations in Trms of Vlocitis and Prssur (2D) Diffrntial quations v v v σ σ ρ + v + v f = t 0 JN Rdd v v v σ σ ρ + v + v f = 0 t v v + = 0 Boundar conditions ( v, t) and ( v, t) v v v v σ = 2µ P, σ = 2µ P, σ = µ + Viscous Incompr. Flows: 4

5 WEAK FORMS OF THE EQUATIONS 0 0 v v v σ σ = w1 ρ + v + v f dd t Ω v v v w1 w 1 = ρw1 + v + v + σ + σ w1f dd w1t ds, t Ω v v v σ σ = w2 ρ + v + v f dd t Ω v v v w w = ρw + v + v + σ + σ wf t Ω dd w2 t ds, t = σ n + σ n, t = σ n + σ n w Γ w Γ v 1 v 2 0 v v = w3 + Ω dd w3 P JN Rdd Viscous Incompr. Flows: 5

6 JN Rdd Mid Finit Elmnt Modl for th stad-stat cas m m n v = v, v = v, P = P, j j j j j j j= 1 j= 1 j= M 0 0 v K + G(v) K K v F M 0 v + K K + G(v) K v = F P K K K P F 11 i j j j j i Kij = 2 + dd, Gij( v, v ) = i v + v dd 12 i j K, i, i ij = dd Kij = j dd Kij = j dd, K 22 ij i j i = + 2 j dd, K = K, K = K, K = K, K = ij ji ij ji ij ji ij F = f dd + t ds, F = f dd + t ds 1 2 i i i i i i 0 Viscous Incompr. Flows: 6

7 Mid Finit Elmnt Modl (continud) v = u, v = v v 1 1 P u 1 v 3 3 u 3 v 2 2 u 2 Quadratic (u,v); linar P Quadratic (u,v); linar P Linar (u,v); constant P Linar (u,v); constant P JN Rdd Viscous Incompr. Flows: 7

8 Pnalt Function Mthod-algbraic Problm: Find th minimum of th function F (; ) subjct to th constraint G(; ) = Lagrang multiplir mthod d = 0 F L (; ; ) F (; ) + G(; ) df d d d = 0 d + G(; @ @ = 0; G(; ) = 0 Viscous Incompr. Flows: 8

9 Pnalt Function Mthod-algbraic (continud) Pnalt function mthod F P (; ) = F (; ) + [G(; ) 0]2 2 df d d = 0 = + G(; )@G @ + G(; d + G(; + G(; = 0 JN Rdd Viscous Incompr. Flows: 9

10 df d d d = 0 µ = function d + G(; F P (; ) = F (; ) [G(; ) 0]2 + = 0; G(; ) = df d d = 0 = + G(; )@G @ + G(; @ + G(; + G(; = 0 JN Rdd Approimation of th Lagrang multiplid can b computd in th pnalt mthod from λ γ = γ G( γ, γ ) Viscous Incompr. Flows: 10

11 Pnalt Function Mthod-algbraic (An Eampl) F (; ) = ; G(; ) 2 = 0 Lagrang Multiplir Mthod = 0; = 0; 2 = 0 Pnalt Function Mthod = 0:5; = 1:0; = 3: (2 ) = 0; (2 ) = 0 Clarl, as! 1, w hav = ; = lim!1 = 0:5 = ; lim!1 = 1:0 = JN Rdd Viscous Incompr. Flows: 11

12 Pnalt Function Mthod-algbraic (Eampl - continud) Tabl: A comparison of th pnalt solution with th act for various valus of th pnalt paramtr G( ; ) λγ = γ G( γ, γ ) JN Rdd Viscous Incompr. Flows: 12

13 0 PENALTY FINITE ELEMENT FORMULATION for th Stad-Stat Cas Considr th wak forms v v v w2 w 2 0 = ρw2 + v + v + σ + σ w2f dd w2tds, t 0 v v v w1 w 1 = ρw1 + v + v + σ + σ w1f dd w1t ds, t Ω Ω = 3 + Ω w v v dd, t = σ n + σ n, t = σ n + σ n Γ Γ JN Rdd Viscous Incompr. Flows: 13

14 Pnalt Finit Elmnt Formulation (continud) Now suppos that th vlocit fild satisfis th constraint v v v v w w = = 0 Thn adding th thr wak statmnts, w obtain w v w v 1 2 w w v v w w v v 1 2 P w w f w f dd v v v w v v v w v 1 2 v dd w t w t ds 1 2 JN Rdd Viscous Incompr. Flows: 14

15 JN Rdd Pnalt Finit Elmnt Formulation - continud Thus, th wak form of th problm, subjctd to th constraint is w1 u w2 v w w 1 2 u v w f w f dd w t w t ds v v v v w v v w v v dd 1 2 Thus, th variational problm is 0 I( v, v) v subjct to th constraint v + = 0 Viscous Incompr. Flows: 15

16 JN Rdd Pnalt Finit Elmnt Formulation - continud Thn, th modifid wak form with th constraint is 0 v v I p( v, v) I( v, v) dd 2 w1 u w2 v w w 1 2 u v w f w f dd w t w t ds v v v v w v v w v v dd 1 2 w1 w 2 v v dd 2 Pnalt prssion Viscous Incompr. Flows: 16

17 Pnalt Finit Elmnt Formulation - continud Th wak form of th problm can b sparatd into th following two statmnts: w v w v v 1 1 w v v w f dd 1 v v w t ds w v v dd 1 1 w v 2 w v v 2 w v v w f dd 2 v v w t ds w v v dd 2 2 Ths statmnts form th basis of th pnalt FE Modl. JN Rdd Viscous Incompr. Flows: 17

18 JN Rdd Pnalt Finit Elmnt Formulation continud Altrnativl, th prssur (ngativ of th Lagrang multiplir) in th govrning quations can b rplacd b v v P W obtain v v v v v v f 0 v v v v v v f 0 Th wak forms of ths quations ar prcisl th sam as thos on th prvious slid. Viscous Incompr. Flows: 18

19 JN Rdd Pnalt Finit Elmnt Modl m v = v (, ), v = v (, ) j j j j j= 1 j= 1 Substitution into th wak forms (adding inrtia trms) ilds th quations M 0 v K + Gv ( ) K v F M v = K K + Gv ( ) v F 11 j j j K 2 i i i ij = + dd + dd K = dd + dd, K = K 12 i j i j ij ij ji 22 i j j j K 2 i i ij = + dd + dd F = f dd + t 1 i i m ds, F = f dd + t ds 2 i i i i Viscous Incompr. Flows: 19

20 Elmnts Usd for Pnalt FE Modl v 3 v 2 2 v 2 3 v 3 v 1 1 v 1 Nods with v and v JN Rdd Viscous Incompr. Flows: 20

21 Computational Aspcts of th Pnalt FEM Gnral form of th Pnalt FEM: ¹[K 1 ] + ½[K 2 ] + [K 3 ] f g = ff g Elmnt `locking': lim! 0 ¹[K 1 ] + ½[K 2 ] + [K 3 ] f g = ff g! [K 3 ]f g = ff g Choic of th pnalt paramtr: = 10 4 ¹ to = ¹ Rducd intgration of th pnalt trms JN Rdd Viscous Incompr. Flows: 21

22 Numrical Eampls Viscous fluid squzd btwn paralll plats V 0 Computational domain Linar lmnt Quadratic lmnt v = V0, v = 0 2b 2a v t = 0 = 0 t t = 0 = 0 V 0 (a) v 0, t = 0 = (b) JN Rdd Viscous Incompr. Flows: 22

23 Viscous fluid squzd btwn paralll plats: Vlocit fild Distanc, v at = 4 Analtical solution v at = 6 0.0, v Horizontal vlocit JN Rdd Viscous Incompr. Flows: 23

24 Viscous fluid squzd btwn paralll plats: Prssur fild Prssur, P Analtical solution ( = 2) FEM solution ( = ) Prssur, P Analtical solution ( = 0) FEM solution ( = ) Distanc, Distanc, JN Rdd Viscous Incompr. Flows: 24

25 LID-DRIVEN CAVITY FLOW v = 1, v = 0 Vr long v = 0, v = 0 a a 0.5 a v = 0, v = 0 v 0, v = 0 = JN Rdd Viscous Incompr. Flows: 25

26 1.0 LID-DRIVEN CAVITY FLOW Distanc along vrtical cntrlin, Linar lmnt (88 msh) Quadratic lmnt (44 msh) Linar lmnt (88 msh) Quadratic lmnt (44 msh) Circls dnot pnalt paramtr of 10^8 Squars dnot pnalt paramtr of 10^ JN Rdd Horizontal vlocit, u Viscous Incompr. Flows: 26

27 Wall-drivn cavit flow vlocit profils Distanc, Pnalt FEM 8 8 msh msh Horizontal vlocit, v (0.5,) JN Rdd Viscous Incompr. Flows: 27

28 Wall-Drivn Cavit rsults (continud) Distanc, msh of bilinar lmnts ( R = 500, β = 0.5, ε = 10 2 ) R = 10,000 R = 5,000 R = 1,000 R = 500 JN Rdd Horizontal vlocit, v (0.5,) Viscous Incompr. Flows: 28

29 Stramlins 4 R = 10 Prssur contours Dilatation contours JN Rdd FluidMch LSFEM - 29 Viscous Incompr. Flows: 29

30 SUMMARY Th following topics wr covrd in this lctur: Govrning quations of flows of incomprssibl fluids Mid (vlocit-prssur) finit lmnt modl Pnalt function mthod - algbraic problm Pnalt finit lmnt modl of viscous incomprssibl fluids Numrical rsults JN Rdd Viscous Incompr. Flows: 30

NONLINEAR ANALYSIS OF PLATE BENDING

NONLINEAR ANALYSIS OF PLATE BENDING NONLINEAR ANALYSIS OF PLATE BENDING CONTENTS Govrning Equations of th First-Ordr Shar Dformation thor (FSDT) Finit lmnt modls of FSDT Shar and mmbran locking Computr implmntation Strss calculation Numrical

More information

AS 5850 Finite Element Analysis

AS 5850 Finite Element Analysis AS 5850 Finit Elmnt Analysis Two-Dimnsional Linar Elasticity Instructor Prof. IIT Madras Equations of Plan Elasticity - 1 displacmnt fild strain- displacmnt rlations (infinitsimal strain) in matrix form

More information

OTHER TPOICS OF INTEREST (Convection BC, Axisymmetric problems, 3D FEM)

OTHER TPOICS OF INTEREST (Convection BC, Axisymmetric problems, 3D FEM) OTHER TPOICS OF INTEREST (Convction BC, Axisymmtric problms, 3D FEM) CONTENTS 2-D Problms with convction BC Typs of Axisymmtric Problms Axisymmtric Problms (2-D) 3-D Hat Transfr 3-D Elasticity Typical

More information

Nonlinear Bending of Strait Beams

Nonlinear Bending of Strait Beams Nonlinar Bnding of Strait Bams CONTENTS Th Eulr-Brnoulli bam thory Th Timoshnko bam thory Govrning Equations Wak Forms Finit lmnt modls Computr Implmntation: calculation of lmnt matrics Numrical ampls

More information

FEM FOR HEAT TRANSFER PROBLEMS دانشگاه صنعتي اصفهان- دانشكده مكانيك

FEM FOR HEAT TRANSFER PROBLEMS دانشگاه صنعتي اصفهان- دانشكده مكانيك FEM FOR HE RNSFER PROBLEMS 1 Fild problms Gnral orm o systm quations o D linar stady stat ild problms: For 1D problms: D D g Q y y (Hlmholtz quation) d D g Q d Fild problms Hat transr in D in h h ( D D

More information

Finite element discretization of Laplace and Poisson equations

Finite element discretization of Laplace and Poisson equations Finit lmnt discrtization of Laplac and Poisson quations Yashwanth Tummala Tutor: Prof S.Mittal 1 Outlin Finit Elmnt Mthod for 1D Introduction to Poisson s and Laplac s Equations Finit Elmnt Mthod for 2D-Discrtization

More information

COMPUTATIONAL NUCLEAR THERMAL HYDRAULICS

COMPUTATIONAL NUCLEAR THERMAL HYDRAULICS COMPUTTIONL NUCLER THERML HYDRULICS Cho, Hyoung Kyu Dpartmnt of Nuclar Enginring Soul National Univrsity CHPTER4. THE FINITE VOLUME METHOD FOR DIFFUSION PROBLEMS 2 Tabl of Contnts Chaptr 1 Chaptr 2 Chaptr

More information

1 Isoparametric Concept

1 Isoparametric Concept UNIVERSITY OF CALIFORNIA BERKELEY Dpartmnt of Civil Enginring Spring 06 Structural Enginring, Mchanics and Matrials Profssor: S. Govindj Nots on D isoparamtric lmnts Isoparamtric Concpt Th isoparamtric

More information

2008 AP Calculus BC Multiple Choice Exam

2008 AP Calculus BC Multiple Choice Exam 008 AP Multipl Choic Eam Nam 008 AP Calculus BC Multipl Choic Eam Sction No Calculator Activ AP Calculus 008 BC Multipl Choic. At tim t 0, a particl moving in th -plan is th acclration vctor of th particl

More information

CHAPTER 9. Interpolation functions for 2D elements. Numerical Integration. Modeling Considerations

CHAPTER 9. Interpolation functions for 2D elements. Numerical Integration. Modeling Considerations HAPTER 9 Intrpolation functions for D lmnts Numrical Intgration Modling onsidrations Pascal s triangl Dgr of Numbr of Elmnt with th complt trms in th nods polnomial polnomial Triangular Elmnts (Figur not

More information

Chapter 5. Introduction. Introduction. Introduction. Finite Element Modelling. Finite Element Modelling

Chapter 5. Introduction. Introduction. Introduction. Finite Element Modelling. Finite Element Modelling Chaptr 5 wo-dimnsional problms using Constant Strain riangls (CS) Lctur Nots Dr Mohd Andi Univrsiti Malasia Prlis EN7 Finit Elmnt Analsis Introction wo-dimnsional init lmnt ormulation ollows th stps usd

More information

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME 43 Introduction to Finit Elmnt Analysis Chaptr 3 Computr Implmntation of D FEM Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt

More information

Higher-Order Discrete Calculus Methods

Higher-Order Discrete Calculus Methods Highr-Ordr Discrt Calculus Mthods J. Blair Prot V. Subramanian Ralistic Practical, Cost-ctiv, Physically Accurat Paralll, Moving Msh, Complx Gomtry, Slid 1 Contxt Discrt Calculus Mthods Finit Dirnc Mimtic

More information

Derivation of Eigenvalue Matrix Equations

Derivation of Eigenvalue Matrix Equations Drivation of Eignvalu Matrix Equations h scalar wav quations ar φ φ η + ( k + 0ξ η β ) φ 0 x y x pq ε r r whr for E mod E, 1, y pq φ φ x 1 1 ε r nr (4 36) for E mod H,, 1 x η η ξ ξ n [ N ] { } i i i 1

More information

Sec 2.3 Modeling with First Order Equations

Sec 2.3 Modeling with First Order Equations Sc.3 Modling with First Ordr Equations Mathmatical modls charactriz physical systms, oftn using diffrntial quations. Modl Construction: Translating physical situation into mathmatical trms. Clarly stat

More information

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME Introduction to Finit Elmnt Analysis Chaptr 5 Two-Dimnsional Formulation Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt

More information

Section 11.6: Directional Derivatives and the Gradient Vector

Section 11.6: Directional Derivatives and the Gradient Vector Sction.6: Dirctional Drivativs and th Gradint Vctor Practic HW rom Stwart Ttbook not to hand in p. 778 # -4 p. 799 # 4-5 7 9 9 35 37 odd Th Dirctional Drivativ Rcall that a b Slop o th tangnt lin to th

More information

3 Finite Element Parametric Geometry

3 Finite Element Parametric Geometry 3 Finit Elmnt Paramtric Gomtry 3. Introduction Th intgral of a matrix is th matrix containing th intgral of ach and vry on of its original componnts. Practical finit lmnt analysis rquirs intgrating matrics,

More information

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values Dnamic Macroconomic Thor Prof. Thomas Lux Bifurcation Thor Bifurcation: qualitativ chang in th natur of th solution occurs if a paramtr passs through a critical point bifurcation or branch valu. Local

More information

( ) Differential Equations. Unit-7. Exact Differential Equations: M d x + N d y = 0. Verify the condition

( ) Differential Equations. Unit-7. Exact Differential Equations: M d x + N d y = 0. Verify the condition Diffrntial Equations Unit-7 Eat Diffrntial Equations: M d N d 0 Vrif th ondition M N Thn intgrat M d with rspt to as if wr onstants, thn intgrat th trms in N d whih do not ontain trms in and quat sum of

More information

Calculus II (MAC )

Calculus II (MAC ) Calculus II (MAC232-2) Tst 2 (25/6/25) Nam (PRINT): Plas show your work. An answr with no work rcivs no crdit. You may us th back of a pag if you nd mor spac for a problm. You may not us any calculators.

More information

VTU NOTES QUESTION PAPERS NEWS RESULTS FORUMS

VTU NOTES QUESTION PAPERS NEWS RESULTS FORUMS Diffrntial Equations Unit-7 Eat Diffrntial Equations: M d N d 0 Vrif th ondition M N Thn intgrat M d with rspt to as if wr onstants, thn intgrat th trms in N d whih do not ontain trms in and quat sum of

More information

Division of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):.

Division of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):. Division of Mchanics Lund Univrsity MULTIBODY DYNMICS Examination 7033 Nam (writ in block lttrs):. Id.-numbr: Writtn xamination with fiv tasks. Plas chck that all tasks ar includd. clan copy of th solutions

More information

(1) Then we could wave our hands over this and it would become:

(1) Then we could wave our hands over this and it would become: MAT* K285 Spring 28 Anthony Bnoit 4/17/28 Wk 12: Laplac Tranform Rading: Kohlr & Johnon, Chaptr 5 to p. 35 HW: 5.1: 3, 7, 1*, 19 5.2: 1, 5*, 13*, 19, 45* 5.3: 1, 11*, 19 * Pla writ-up th problm natly and

More information

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME 4 Introduction to Finit Elmnt Analysis Chaptr 4 Trusss, Bams and Frams Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt

More information

The Finite Element Method

The Finite Element Method Th Finit Elmnt Mthod Eulr-Brnoulli and Timoshnko Bams Rad: Chaptr 5 CONTENTS Eulr-Brnoulli bam thory Govrning Equations Finit lmnt modl Numrical ampls Timoshnko bam thory Govrning Equations Finit lmnt

More information

Construction of Mimetic Numerical Methods

Construction of Mimetic Numerical Methods Construction of Mimtic Numrical Mthods Blair Prot Thortical and Computational Fluid Dynamics Laboratory Dltars July 17, 013 Numrical Mthods Th Foundation on which CFD rsts. Rvolution Math: Accuracy Stability

More information

6.1 Integration by Parts and Present Value. Copyright Cengage Learning. All rights reserved.

6.1 Integration by Parts and Present Value. Copyright Cengage Learning. All rights reserved. 6.1 Intgration by Parts and Prsnt Valu Copyright Cngag Larning. All rights rsrvd. Warm-Up: Find f () 1. F() = ln(+1). F() = 3 3. F() =. F() = ln ( 1) 5. F() = 6. F() = - Objctivs, Day #1 Studnts will b

More information

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C Tchniqus of Intgration c Donald Kridr and Dwight Lahr In this sction w ar going to introduc th first approachs to valuating an indfinit intgral whos intgrand dos not hav an immdiat antidrivativ. W bgin

More information

Differentiation of Exponential Functions

Differentiation of Exponential Functions Calculus Modul C Diffrntiation of Eponntial Functions Copyright This publication Th Northrn Albrta Institut of Tchnology 007. All Rights Rsrvd. LAST REVISED March, 009 Introduction to Diffrntiation of

More information

Higher order derivatives

Higher order derivatives Robrto s Nots on Diffrntial Calculus Chaptr 4: Basic diffrntiation ruls Sction 7 Highr ordr drivativs What you nd to know alrady: Basic diffrntiation ruls. What you can larn hr: How to rpat th procss of

More information

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields Lctur 37 (Schrödingr Equation) Physics 6-01 Spring 018 Douglas Filds Rducd Mass OK, so th Bohr modl of th atom givs nrgy lvls: E n 1 k m n 4 But, this has on problm it was dvlopd assuming th acclration

More information

AP PHYSICS C: ELECTRICITY AND MAGNETISM 2015 SCORING GUIDELINES

AP PHYSICS C: ELECTRICITY AND MAGNETISM 2015 SCORING GUIDELINES AP PHYSICS C: ELECTRICITY AND MAGNETISM 2015 SCORING GUIDELINES Qustion 1 15 points total Distribution of points (a) i. For at last on arrow btwn th plats pointing downward from th positiv plats toward

More information

Abstract Interpretation: concrete and abstract semantics

Abstract Interpretation: concrete and abstract semantics Abstract Intrprtation: concrt and abstract smantics Concrt smantics W considr a vry tiny languag that manags arithmtic oprations on intgrs valus. Th (concrt) smantics of th languags cab b dfind by th funzcion

More information

DIFFERENTIAL EQUATION

DIFFERENTIAL EQUATION MD DIFFERENTIAL EQUATION Sllabus : Ordinar diffrntial quations, thir ordr and dgr. Formation of diffrntial quations. Solution of diffrntial quations b th mthod of sparation of variabls, solution of homognous

More information

GRAnada COde for the resolution of the adiabatic and non-adiabatic stellar oscillations. General scheme

GRAnada COde for the resolution of the adiabatic and non-adiabatic stellar oscillations. General scheme GRAnada COd for th rsolution of th adiabatic and non-adiabatic stllar oscillations A. Moa Instituto d Astrofísica d Andalucía, CSIC, Granada, Spain Gnral schm Equilibrium modls Standard intrior Connctin

More information

1997 AP Calculus AB: Section I, Part A

1997 AP Calculus AB: Section I, Part A 997 AP Calculus AB: Sction I, Part A 50 Minuts No Calculator Not: Unlss othrwis spcifid, th domain of a function f is assumd to b th st of all ral numbrs for which f () is a ral numbr.. (4 6 ) d= 4 6 6

More information

16. Electromagnetics and vector elements (draft, under construction)

16. Electromagnetics and vector elements (draft, under construction) 16. Elctromagntics (draft)... 1 16.1 Introduction... 1 16.2 Paramtric coordinats... 2 16.3 Edg Basd (Vctor) Finit Elmnts... 4 16.4 Whitny vctor lmnts... 5 16.5 Wak Form... 8 16.6 Vctor lmnt matrics...

More information

Hydrogen Atom and One Electron Ions

Hydrogen Atom and One Electron Ions Hydrogn Atom and On Elctron Ions Th Schrödingr quation for this two-body problm starts out th sam as th gnral two-body Schrödingr quation. First w sparat out th motion of th cntr of mass. Th intrnal potntial

More information

Inference Methods for Stochastic Volatility Models

Inference Methods for Stochastic Volatility Models Intrnational Mathmatical Forum, Vol 8, 03, no 8, 369-375 Infrnc Mthods for Stochastic Volatility Modls Maddalna Cavicchioli Cá Foscari Univrsity of Vnic Advancd School of Economics Cannargio 3, Vnic, Italy

More information

LEAST-SQUARES FINITE ELEMENT MODELS

LEAST-SQUARES FINITE ELEMENT MODELS LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows

More information

Elements of Statistical Thermodynamics

Elements of Statistical Thermodynamics 24 Elmnts of Statistical Thrmodynamics Statistical thrmodynamics is a branch of knowldg that has its own postulats and tchniqus. W do not attmpt to giv hr vn an introduction to th fild. In this chaptr,

More information

1973 AP Calculus AB: Section I

1973 AP Calculus AB: Section I 97 AP Calculus AB: Sction I 9 Minuts No Calculator Not: In this amination, ln dnots th natural logarithm of (that is, logarithm to th bas ).. ( ) d= + C 6 + C + C + C + C. If f ( ) = + + + and ( ), g=

More information

NEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA

NEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA NE APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA Mirca I CÎRNU Ph Dp o Mathmatics III Faculty o Applid Scincs Univrsity Polithnica o Bucharst Cirnumirca @yahoocom Abstract In a rcnt papr [] 5 th indinit intgrals

More information

Analysis of potential flow around two-dimensional body by finite element method

Analysis of potential flow around two-dimensional body by finite element method Vol. 7(2), pp. 9-22, May, 2015 DOI: 10.5897/JMER2014.0342 rticl Numbr: 20E80053033 ISSN 2141 2383 Copyright 2015 uthor(s) rtain th copyright of this articl http://www.acadmicjournals.org/jmer Journal of

More information

Southern Taiwan University

Southern Taiwan University Chaptr Ordinar Diffrntial Equations of th First Ordr and First Dgr Gnral form:., d +, d 0.a. f,.b I. Sparabl Diffrntial quations Form: d + d 0 C d d E 9 + 4 0 Solution: 9d + 4d 0 9 + 4 C E + d Solution:

More information

CHAPTER 2 LAGRANGIAN AND EULERIAN FINITE ELEMENTS IN ONE DIMENSION

CHAPTER 2 LAGRANGIAN AND EULERIAN FINITE ELEMENTS IN ONE DIMENSION CHAPTER 2 LAGRANGIAN AND EULERIAN FINITE ELEMENTS IN ONE DIMENSION by Td Blytschko Northwstrn Univrsity @ Copyright 1997 2.1 Introduction In this chaptr, th quations for on-dimnsional modls of nonlinar

More information

Direct Approach for Discrete Systems One-Dimensional Elements

Direct Approach for Discrete Systems One-Dimensional Elements CONTINUUM & FINITE ELEMENT METHOD Dirct Approach or Discrt Systms On-Dimnsional Elmnts Pro. Song Jin Par Mchanical Enginring, POSTECH Dirct Approach or Discrt Systms Dirct approach has th ollowing aturs:

More information

INTEGRATION BY PARTS

INTEGRATION BY PARTS Mathmatics Rvision Guids Intgration by Parts Pag of 7 MK HOME TUITION Mathmatics Rvision Guids Lvl: AS / A Lvl AQA : C Edcl: C OCR: C OCR MEI: C INTEGRATION BY PARTS Vrsion : Dat: --5 Eampls - 6 ar copyrightd

More information

Numerical methods for PDEs FEM implementation: element stiffness matrix, isoparametric mapping, assembling global stiffness matrix

Numerical methods for PDEs FEM implementation: element stiffness matrix, isoparametric mapping, assembling global stiffness matrix Platzhaltr für Bild, Bild auf Titlfoli hintr das Logo instzn Numrical mthods for PDEs FEM implmntation: lmnt stiffnss matrix, isoparamtric mapping, assmbling global stiffnss matrix Dr. Nomi Fridman Contnts

More information

MA 262, Spring 2018, Final exam Version 01 (Green)

MA 262, Spring 2018, Final exam Version 01 (Green) MA 262, Spring 218, Final xam Vrsion 1 (Grn) INSTRUCTIONS 1. Switch off your phon upon ntring th xam room. 2. Do not opn th xam booklt until you ar instructd to do so. 3. Bfor you opn th booklt, fill in

More information

ME469A Numerical Methods for Fluid Mechanics

ME469A Numerical Methods for Fluid Mechanics ME469A Numrical Mthods for Fluid Mchanics Handout #5 Gianluca Iaccarino Finit Volum Mthods Last tim w introducd th FV mthod as a discrtization tchniqu applid to th intgral form of th govrning quations

More information

High Energy Physics. Lecture 5 The Passage of Particles through Matter

High Energy Physics. Lecture 5 The Passage of Particles through Matter High Enrgy Physics Lctur 5 Th Passag of Particls through Mattr 1 Introduction In prvious lcturs w hav sn xampls of tracks lft by chargd particls in passing through mattr. Such tracks provid som of th most

More information

Chapter 6 Folding. Folding

Chapter 6 Folding. Folding Chaptr 6 Folding Wintr 1 Mokhtar Abolaz Folding Th folding transformation is usd to systmatically dtrmin th control circuits in DSP architctur whr multipl algorithm oprations ar tim-multiplxd to a singl

More information

VSMN30 FINITA ELEMENTMETODEN - DUGGA

VSMN30 FINITA ELEMENTMETODEN - DUGGA VSMN3 FINITA ELEMENTMETODEN - DUGGA 1-11-6 kl. 8.-1. Maximum points: 4, Rquird points to pass: Assistanc: CALFEM manual and calculator Problm 1 ( 8p ) 8 7 6 5 y 4 1. m x 1 3 1. m Th isotropic two-dimnsional

More information

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design MAE4700/5700 Finit Elmnt Analysis for Mchanical and Arospac Dsign Cornll Univrsity, Fall 2009 Nicholas Zabaras Matrials Procss Dsign and Control Laboratory Sibly School of Mchanical and Arospac Enginring

More information

Finite Element Analysis

Finite Element Analysis Finit Elmnt Analysis L4 D Shap Functions, an Gauss Quaratur FEA Formulation Dr. Wiong Wu EGR 54 Finit Elmnt Analysis Roamap for Dvlopmnt of FE Strong form: govrning DE an BCs EGR 54 Finit Elmnt Analysis

More information

u 3 = u 3 (x 1, x 2, x 3 )

u 3 = u 3 (x 1, x 2, x 3 ) Lctur 23: Curvilinar Coordinats (RHB 8.0 It is oftn convnint to work with variabls othr than th Cartsian coordinats x i ( = x, y, z. For xampl in Lctur 5 w mt sphrical polar and cylindrical polar coordinats.

More information

Maxwell s Equations on a Yee Grid

Maxwell s Equations on a Yee Grid Instructor Dr. Ramond Rumpf (95) 77 6958 rcrumpf@utp.du 57 Computational lctromagntics (CM) Lctur # Mawll s quations on a Y Grid Lctur Ths nots ma contain coprightd matrial obtaind undr fair us ruls. Distribution

More information

Thomas Whitham Sixth Form

Thomas Whitham Sixth Form Thomas Whitham Sith Form Pur Mathmatics Unit C Algbra Trigonomtr Gomtr Calculus Vctor gomtr Pag Algbra Molus functions graphs, quations an inqualitis Graph of f () Draw f () an rflct an part of th curv

More information

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the Copyright itutcom 005 Fr download & print from wwwitutcom Do not rproduc by othr mans Functions and graphs Powr functions Th graph of n y, for n Q (st of rational numbrs) y is a straight lin through th

More information

Integration by Parts

Integration by Parts Intgration by Parts Intgration by parts is a tchniqu primarily for valuating intgrals whos intgrand is th product of two functions whr substitution dosn t work. For ampl, sin d or d. Th rul is: u ( ) v'(

More information

Where k is either given or determined from the data and c is an arbitrary constant.

Where k is either given or determined from the data and c is an arbitrary constant. Exponntial growth and dcay applications W wish to solv an quation that has a drivativ. dy ky k > dx This quation says that th rat of chang of th function is proportional to th function. Th solution is

More information

nd the particular orthogonal trajectory from the family of orthogonal trajectories passing through point (0; 1).

nd the particular orthogonal trajectory from the family of orthogonal trajectories passing through point (0; 1). Eamn EDO. Givn th family of curvs y + C nd th particular orthogonal trajctory from th family of orthogonal trajctoris passing through point (0; ). Solution: In th rst plac, lt us calculat th di rntial

More information

INC 693, 481 Dynamics System and Modelling: The Language of Bound Graphs Dr.-Ing. Sudchai Boonto Assistant Professor

INC 693, 481 Dynamics System and Modelling: The Language of Bound Graphs Dr.-Ing. Sudchai Boonto Assistant Professor INC 693, 48 Dynamics Systm and Modlling: Th Languag o Bound Graphs Dr.-Ing. Sudchai Boonto Assistant Prossor Dpartmnt o Control Systm and Instrumntation Enginring King Mongkut s Unnivrsity o Tchnology

More information

INC 693, 481 Dynamics System and Modelling: Linear Graph Modeling II Dr.-Ing. Sudchai Boonto Assistant Professor

INC 693, 481 Dynamics System and Modelling: Linear Graph Modeling II Dr.-Ing. Sudchai Boonto Assistant Professor INC 69, 48 Dynamics Systm and Modlling: Linar Graph Modling II Dr.-Ing. Sudchai Boonto Assistant Profssor Dpartmnt of Control Systm and Instrumntation Enginring King Mongkut s Unnivrsity of Tchnology Thonuri

More information

Einstein Equations for Tetrad Fields

Einstein Equations for Tetrad Fields Apiron, Vol 13, No, Octobr 006 6 Einstin Equations for Ttrad Filds Ali Rıza ŞAHİN, R T L Istanbul (Turky) Evry mtric tnsor can b xprssd by th innr product of ttrad filds W prov that Einstin quations for

More information

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002 3.4 Forc Analysis of Linkas An undrstandin of forc analysis of linkas is rquird to: Dtrmin th raction forcs on pins, tc. as a consqunc of a spcifid motion (don t undrstimat th sinificanc of dynamic or

More information

surface of a dielectric-metal interface. It is commonly used today for discovering the ways in

surface of a dielectric-metal interface. It is commonly used today for discovering the ways in Surfac plasmon rsonanc is snsitiv mchanism for obsrving slight changs nar th surfac of a dilctric-mtal intrfac. It is commonl usd toda for discovring th was in which protins intract with thir nvironmnt,

More information

u x A j Stress in the Ocean

u x A j Stress in the Ocean Strss in t Ocan T tratmnt of strss and strain in fluids is comlicatd and somwat bond t sco of tis class. Tos rall intrstd sould look into tis rtr in Batclor Introduction to luid Dnamics givn as a rfrnc

More information

General Notes About 2007 AP Physics Scoring Guidelines

General Notes About 2007 AP Physics Scoring Guidelines AP PHYSICS C: ELECTRICITY AND MAGNETISM 2007 SCORING GUIDELINES Gnral Nots About 2007 AP Physics Scoring Guidlins 1. Th solutions contain th most common mthod of solving th fr-rspons qustions and th allocation

More information

1 General boundary conditions in diffusion

1 General boundary conditions in diffusion Gnral boundary conditions in diffusion Πρόβλημα 4.8 : Δίνεται μονοδιάτατη πλάκα πάχους, που το ένα άκρο της κρατιέται ε θερμοκραία T t και το άλλο ε θερμοκραία T 2 t. Αν η αρχική θερμοκραία της πλάκας

More information

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches. Subjct Chmistry Papr No and Titl Modul No and Titl Modul Tag 8/ Physical Spctroscopy / Brakdown of th Born-Oppnhimr approximation. Slction ruls for rotational-vibrational transitions. P, R branchs. CHE_P8_M

More information

SCALING OF SYNCHROTRON RADIATION WITH MULTIPOLE ORDER. J. C. Sprott

SCALING OF SYNCHROTRON RADIATION WITH MULTIPOLE ORDER. J. C. Sprott SCALING OF SYNCHROTRON RADIATION WITH MULTIPOLE ORDER J. C. Sprott PLP 821 Novmbr 1979 Plasma Studis Univrsity of Wisconsin Ths PLP Rports ar informal and prliminary and as such may contain rrors not yt

More information

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula 7. Intgration by Parts Each drivativ formula givs ris to a corrsponding intgral formula, as w v sn many tims. Th drivativ product rul yilds a vry usful intgration tchniqu calld intgration by parts. Starting

More information

Exercise 1. Sketch the graph of the following function. (x 2

Exercise 1. Sketch the graph of the following function. (x 2 Writtn tst: Fbruary 9th, 06 Exrcis. Sktch th graph of th following function fx = x + x, spcifying: domain, possibl asymptots, monotonicity, continuity, local and global maxima or minima, and non-drivability

More information

UNIT I PARTIAL DIFFERENTIAL EQUATIONS PART B. 3) Form the partial differential equation by eliminating the arbitrary functions

UNIT I PARTIAL DIFFERENTIAL EQUATIONS PART B. 3) Form the partial differential equation by eliminating the arbitrary functions UNIT I PARTIAL DIFFERENTIAL EQUATIONS PART B 1) Form th artial diffrntial quation b liminating th arbitrar functions f and g in z f ( x ) g( x ) ) Form th artial diffrntial quation b liminating th arbitrar

More information

SAFE HANDS & IIT-ian's PACE EDT-15 (JEE) SOLUTIONS

SAFE HANDS & IIT-ian's PACE EDT-15 (JEE) SOLUTIONS It is not possibl to find flu through biggr loop dirctly So w will find cofficint of mutual inductanc btwn two loops and thn find th flu through biggr loop Also rmmbr M = M ( ) ( ) EDT- (JEE) SOLUTIONS

More information

Broadband All-Angle Negative Refraction by Phononic Crystals

Broadband All-Angle Negative Refraction by Phononic Crystals Supplmntar Information Broadband All-Angl Ngativ Rfraction b Phononic Crstals Yang Fan Li, Fi Mng, Shiwi Zhou, Ming-Hui Lu and Xiaodong Huang 1 Optimization algorithm and procss Bfor th optimization procss,

More information

Finite Strain Elastic-Viscoplastic Model

Finite Strain Elastic-Viscoplastic Model Finit Strain Elastic-Viscoplastic Modl Pinksh Malhotra Mchanics of Solids,Brown Univrsity Introduction Th main goal of th projct is to modl finit strain rat-dpndnt plasticity using a modl compatibl for

More information

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0 unction Spacs Prrquisit: Sction 4.7, Coordinatization n this sction, w apply th tchniqus of Chaptr 4 to vctor spacs whos lmnts ar functions. Th vctor spacs P n and P ar familiar xampls of such spacs. Othr

More information

Rational Approximation for the one-dimensional Bratu Equation

Rational Approximation for the one-dimensional Bratu Equation Intrnational Journal of Enginring & Tchnology IJET-IJES Vol:3 o:05 5 Rational Approximation for th on-dimnsional Bratu Equation Moustafa Aly Soliman Chmical Enginring Dpartmnt, Th British Univrsity in

More information

2. Background Material

2. Background Material S. Blair Sptmbr 3, 003 4. Background Matrial Th rst of this cours dals with th gnration, modulation, propagation, and ction of optical radiation. As such, bic background in lctromagntics and optics nds

More information

First derivative analysis

First derivative analysis Robrto s Nots on Dirntial Calculus Chaptr 8: Graphical analysis Sction First drivativ analysis What you nd to know alrady: How to us drivativs to idntiy th critical valus o a unction and its trm points

More information

CIVL 8/ D Boundary Value Problems - Rectangular Elements 1/7

CIVL 8/ D Boundary Value Problems - Rectangular Elements 1/7 CIVL / -D Boundr Vlu Prolms - Rctngulr Elmnts / RECANGULAR ELEMENS - In som pplictions, it m mor dsirl to us n lmntl rprsnttion of th domin tht hs four sids, ithr rctngulr or qudriltrl in shp. Considr

More information

Dynamic Modelling of Hoisting Steel Wire Rope. Da-zhi CAO, Wen-zheng DU, Bao-zhu MA *

Dynamic Modelling of Hoisting Steel Wire Rope. Da-zhi CAO, Wen-zheng DU, Bao-zhu MA * 17 nd Intrnational Confrnc on Mchanical Control and Automation (ICMCA 17) ISBN: 978-1-6595-46-8 Dynamic Modlling of Hoisting Stl Wir Rop Da-zhi CAO, Wn-zhng DU, Bao-zhu MA * and Su-bing LIU Xi an High

More information

Discontinuous Galerkin and Mimetic Finite Difference Methods for Coupled Stokes-Darcy Flows on Polygonal and Polyhedral Grids

Discontinuous Galerkin and Mimetic Finite Difference Methods for Coupled Stokes-Darcy Flows on Polygonal and Polyhedral Grids Discontinuous Galrkin and Mimtic Finit Diffrnc Mthods for Coupld Stoks-Darcy Flows on Polygonal and Polyhdral Grids Konstantin Lipnikov Danail Vassilv Ivan Yotov Abstract W study locally mass consrvativ

More information

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers:

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers: APPM 6 Final 5 pts) Spring 4. 6 pts total) Th following parts ar not rlatd, justify your answrs: a) Considr th curv rprsntd by th paramtric quations, t and y t + for t. i) 6 pts) Writ down th corrsponding

More information

An Investigation on the Effect of the Coupled and Uncoupled Formulation on Transient Seepage by the Finite Element Method

An Investigation on the Effect of the Coupled and Uncoupled Formulation on Transient Seepage by the Finite Element Method Amrican Journal of Applid Scincs 4 (1): 95-956, 7 ISSN 1546-939 7 Scinc Publications An Invstigation on th Effct of th Coupld and Uncoupld Formulation on Transint Spag by th Finit Elmnt Mthod 1 Ahad Ouria,

More information

JOHNSON COUNTY COMMUNITY COLLEGE Calculus I (MATH 241) Final Review Fall 2016

JOHNSON COUNTY COMMUNITY COLLEGE Calculus I (MATH 241) Final Review Fall 2016 JOHNSON COUNTY COMMUNITY COLLEGE Calculus I (MATH ) Final Rviw Fall 06 Th Final Rviw is a starting point as you study for th final am. You should also study your ams and homwork. All topics listd in th

More information

Another view for a posteriori error estimates for variational inequalities of the second kind

Another view for a posteriori error estimates for variational inequalities of the second kind Accptd by Applid Numrical Mathmatics in July 2013 Anothr viw for a postriori rror stimats for variational inqualitis of th scond kind Fi Wang 1 and Wimin Han 2 Abstract. In this papr, w giv anothr viw

More information

Molecules and Covalent Bond

Molecules and Covalent Bond Molculs and ovalnt ond Qustion Papr 1 Lvl IGSE Subjct hmistry (0620/0971) Exam oard ambridg Intrnational Examinations (IE) Topic toms, lmnts and compounds Sub-Topic Molculs and covalnt bonds ooklt Qustion

More information

That is, we start with a general matrix: And end with a simpler matrix:

That is, we start with a general matrix: And end with a simpler matrix: DIAGON ALIZATION OF THE STR ESS TEN SOR INTRO DUCTIO N By th us of Cauchy s thorm w ar abl to rduc th numbr of strss componnts in th strss tnsor to only nin valus. An additional simplification of th strss

More information

Electromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology

Electromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology Elctromagntic scattring Graduat Cours Elctrical Enginring (Communications) 1 st Smstr, 1388-1389 Sharif Univrsity of Tchnology Contnts of lctur 8 Contnts of lctur 8: Scattring from small dilctric objcts

More information

Massachusetts Institute of Technology Department of Mechanical Engineering

Massachusetts Institute of Technology Department of Mechanical Engineering Massachustts Institut of Tchnolog Dpartmnt of Mchanical Enginring. Introduction to Robotics Mid-Trm Eamination Novmbr, 005 :0 pm 4:0 pm Clos-Book. Two shts of nots ar allowd. Show how ou arrivd at our

More information

dx equation it is called a second order differential equation.

dx equation it is called a second order differential equation. TOPI Diffrntial quations Mthods of thir intgration oncption of diffrntial quations An quation which spcifis a rlationship btwn a function, its argumnt and its drivativs of th first, scond, tc ordr is calld

More information

Pipe flow friction, small vs. big pipes

Pipe flow friction, small vs. big pipes Friction actor (t/0 t o pip) Friction small vs larg pips J. Chaurtt May 016 It is an intrsting act that riction is highr in small pips than largr pips or th sam vlocity o low and th sam lngth. Friction

More information

Thomas Whitham Sixth Form

Thomas Whitham Sixth Form Thomas Whitham Sith Form Pur Mathmatics Cor rvision gui Pag Algbra Moulus functions graphs, quations an inqualitis Graph of f () Draw f () an rflct an part of th curv blow th ais in th ais. f () f () f

More information

Text: WMM, Chapter 5. Sections , ,

Text: WMM, Chapter 5. Sections , , Lcturs 6 - Continuous Probabilit Distributions Tt: WMM, Chaptr 5. Sctions 6.-6.4, 6.6-6.8, 7.-7. In th prvious sction, w introduc som of th common probabilit distribution functions (PDFs) for discrt sampl

More information