MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design
|
|
- Abel Rice
- 5 years ago
- Views:
Transcription
1 MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design Cornell University, Fall 2009 Nicholas Zabaras Materials Process Design and Control Laboratory Sibley School of Mechanical and Aerospace Engineering 0 Rhodes Hall Cornell University Ithaca, NY
2 Mathematical preliminaries Let us consider a domain with interior and boundary. Usually the boundary is defined parametrically as follows: x(), s y() s, where s is the arc length along. For a function f(x,y) on the boundary, we can write: f ( s) = f( x( s), y( s)), s. We denote the primary (scalar) variable in 2D as u(x,y). Until further notice, we assume that all functions used here are smooth enough for the operations shown to be valid. 2
3 Mathematical preliminaries Recall that the gradient uxy (, ) is the rate of change of u as we move from (x,y) to nearby locations: uxy (, ) uxy (, ) uxy (, ) = i+ j x y Many times we will interpret this as a mathematical operation of the differential operator on uxy (, ). = i+ x y j The gradient uxy (, ) determines the rate of change of u at (x,y) in any direction The rate of change of u in a given direction t = cosθ i+ sinθ j is given as: duxy (, ) uxy (, ) uxy (, ) uxy (, ) t = i = cos θ + sin θ dt x y 3
4 Flux σ ( xy, ) σ ( xy, ) Similarly to the gradient uxy (, ) the flux is a vector-valued function or vector field. The flux () s at the boundary can be decomposed in normal and tangential components: σ () s σ n () s n = i () s σ () s = σ() s τ() s τ σ 4
5 Flux Consider an arbitrary part containing point ( x0, y0) of the domain. Consider the normal flux σ () s σ () s n n = i () s across the boundary The total flux across this boundary is ω. Σ = ω ω ω σ () sds n σ ( xy, ) Σ A ω If you divide ω by the area of the subregion, we can view the result as the amount of flowing into per unit area. σ ω 5
6 Divergence of the flux The limit of Σ ω / Aω σ n ( sds ) / Aω as decreases in size always ω containing point ( x0, y0) is called the divergence of the flux at : div P 0 σ ( x, y ) ω = 0 0 Taking as the square region shown, then: Σ ω =Δσ xδ y +Δσ yδx σ ( xy, ) Dividing by the area ΔxΔy and taking the limit Δx, Δy 0 : σ (, ) y ( xy, ) x xy σ divσ ( x, y) = + = i + j ( σ ) xi + σ y j σ x y x y So σ is the density of the net flux per unit area at a point ω 6
7 Gauss divergence theorem The total flux out of the region is then: Σ= σdxdy We can now write using an earlier expression: σdxdy σ nds This is the Gauss divergence theorem. Σ= = 7
8 Constitutive equation and balance law Constitutive equation: We consider problems with the following constitutive equation between the flux σ and the state variable u(x,y). σ ( x, y) = k( x, y) u( x, y) material modulus k( x, y) > 0in Conservation principle: Within any portion of the domain, the net flux across the boundary of that part must be equal to the total quantity produced by internal sources. ω σ nds = ω fdxdy f is source per unit area 8
9 Balance law Using the divergence theorem: ω σ nds = fdxdy ω ( σ f ) dxdy = 0 for all ω. ω Since ω is arbitrary, we conclude that the local form of the balance equation is: σ ( x, y) = f( x, y) To make things more interesting, we assume an additional source term proportional to u(x,y): σ ( xy, ) + bxyuxy (, ) (, ) = f( xy, ) 9
10 Conservation principle at interfaces Consider an interface Γ separating two materials with different modulus k and k 2. Take a thin strip of the body around a point on this interface. The balance law takes the form: s ( ) ( + ) ( ) 0 2 Σ= σ n+ σ n ds= s The local balance at the interface reduces to: σ ( ) ( ) () s σ + n () s n σ = () s n = 0, s Γ 0
11 Boundary conditions Let us assume that in the boundary Γ we apply natural boundary q 2 conditions: () s () s n () n s = () s = = ps ()(() us us ()), σ σ σ s 2 On the boundary Γ, we assume essential u boundary conditions: us () = us (), s
12 Summary of the problem of interest What we are given:, 2 and the interface Γ The source f = f( x, y) in i, i=,2 ki = ki( xy, ),( xy, ) i, i=,2 The material moduli bi = bi( xy, ),( xy, ) i, i=,2 The boundary coefficients ps (),() us in 2 or σ () s in 2 With this data we want to compute u(x,y) in : ( kxy (, ) uxy (, )) + bxyuxy (, ) (, ) = f( xy, ),( xy, ), i=,2 us () ks () = ps ()[() us us ()], s 2 or n us () ks () = σ (), s s n us () = us (), s 0, k u n = s Γ 2 i 2
13 A note regarding the boundary conditions The special case in which b=0 and when only natural us () BC of the form ks () = σ (), s s 2 is applied and when n 2, needs special attention: The solution u(x,y) can only be determined up to a constant For u to exist, the following compatibility condition needs to be satisfied: fdxdy = σds 3
14 Variational problem: Weak form We start by multiplying the residual r(x,y) by a sufficient smooth function w(x,y) and integrate over each domain in which rw is smooth and set the resulting weighted average equal to zero: rxy (, ) = ( kxy (, ) uxy (, )) + bxyuxy (, ) (, ) f( xy, ) ( kxy (, ) uxy (, )) bxyuxy (, ) (, ) f( xy, ) + wxydxdy (, ) + ( kxy (, ) uxy (, )) bxyuxy (, ) (, ) f( xy, ) + wxydxdy (, ) = 0 2 We need to integrate by parts the first terms in each of these integrals 4
15 2 Integration by parts gives: Variational problem ( kxy (, ) uxy (, )) bxyuxy (, ) (, ) f( xy, ) + wxydxdy (, ) + ( kxy (, ) uxy (, )) bxyuxy (, ) (, ) f( xy, ) + wxydxdy (, ) = 0 kxy (, ) uxy (, ) wxy (, ) bxyuxywxy (, ) (, ) (, ) f( xywxy, ) (, ) + dxdy+ kxy (, ) uxy (, ) wxy (, ) bxyuxywxy (, ) (, ) (, ) f( xywxy, ) (, ) + dxdy 2 ( wk u) dxdy ( wk u) dxdy = 0 2 The last 2 integrals can be simplified with use of the divergence theorem. 5
16 2 ( ) ( 2) Boundary of the domain Variational problem kxy (, ) uxy (, ) wxy (, ) bxyuxywxy (, ) (, ) (, ) f( xywxy, ) (, ) + dxdy+ kxy (, ) uxy (, ) wxy (, ) + bxyuxywxy (, ) (, ) (, ) f( xywxy, ) (, ) dxdy u k wds k n Boundary of the domain 2 u wds = 0 n We can separate integration over Γ in the last 2 terms: u u u u k wds k wds= k wds k wds n n n n ( ) ( ) ( ) Γ ( ) Γ 2 2 u u + ( k ) wds + ( k ) 2wds n n Γ Γ evaluated in region evaluated in region 2 6
17 Variational problem u u u u u u k wds k wds= k wds k wds+ ( k ) wds+ ( k ) 2wds n n n n n n ( ) ( 2) ( ) Γ ( 2) Γ Γ Γ evaluated in region evaluated in region 2 Note that the outward normal n to region is the negative of n 2 at each point on Γ. ( + ) ( ) u u ( + ) u ( ) u ( k ) wds+ ( k ) 2wds= ( k + k ) wds n n n n Γ Γ Γ evaluated in region evaluated in region 2 7
18 Variational problem The last integral is nothing else but: w σ () n s, which is 0! Finally, we can summarize our weak form as: 2 ( + ) ( ) u u ( + ) u ( ) u ( k ) wds+ ( k ) 2wds= ( k + k ) wds n n n n Γ Γ Γ evaluated in region evaluated in region 2 kxy (, ) uxy (, ) wxy (, ) bxyuxywxy (, ) (, ) (, ) f( xywxy, ) (, ) + dxdy+ kxy (, ) uxy (, ) wxy (, ) + bxyuxywxy (, ) (, ) (, ) f( xywxy, ) (, ) dxdy u u k wds k wds= 0 n n ( ) Γ ( 2) Γ u kxy (, ) uxy (, ) wxy (, ) bxyuxywxy (, ) (, ) (, ) f( xywxy, ) (, ) + dxdy k wds= 0 n 8
19 Variational problem u k( x, y) u( x, y) w( x, y) b( x, y) u( x, y) w( x, y) f ( x, y) w( x, y) + dxdy k wds = 0 n Substitution of the natural boundary condition us () ks () = ps ()[() us us ()], s Γq 2 gives: n kxy (, ) uxy (, ) wxy (, ) + bxyuxywxy (, ) (, ) (, ) f( xywxy, ) (, ) dxdy+ pu ( uwds ) = 0 for all admissible functions w( x, y) ( w = 0 on,sin ce we require u( s) = u ( s), s ) We write: Boundary of the whole domain 2 p( u u ) wds = puwds puwds = puwds γ wds γ
20 Space of admissible functions: H () Finally the weak form looks like: k( x, y) u( x, y) w( x, y) b( x, y) u( x, y) w( x, y) + dxdy + puwds = f wdxdy + γ wds for all admissible functions w( x, y), w = 0on 2 2 What is the class of admissible functions? x y 2 2 w w w dxdy < We indicate this space of functions as reflecting that first derivatives are square integrable indicates the domain over which these functions are defined. H ( ), 20
21 Summary of weak problem Find a function uxy (, ) H( ), such that us () = us (), s and the following holds: k( x, y) u( x, y) w( x, y) b( x, y) u( x, y) w( x, y) + dxdy + puwds = f wdxdy + γ wds for all functions w x y H with w = on (, ) ( ) We can repeat the above weak statement using matrix notation as follows (more appropriate for FEM implementation) Find uxy (, ) H( ), us () = us (), s, such that wxy (, ) H( ) withw= 0on the following holds: T ( w) k( x, y) u + b( x, y) u( x, y) w( x, y) dxdy + puwds = f wdxdy + γ wds 2 2 2x x2 x 2
22 Matrix form of the weak statement Find uxy (, ) H( ), us () = us (), s, such that wxy (, ) H( ) withw= 0on the following holds: T ( w) k( x, y) u + b( x, y) u( x, y) w( x, y) dxdy + puwds = f wdxdy + γ wds 2 2 2x x2 x The gradient operator here is defined as a column vector: u x x T w w =, u =, ( w) = u x y y y We will be interchanging the 2 notations here hoping that from the context you know which notation is implied. 22
23 Appendices Some extra slides are included here with: mathematics background, proof of equivalence of weak and strong formulations, Examples of deriving weak forms for 2D BVP, Including anisotropy (e.g. generalized Fourier s law), etc. 23
24 Equivalence between the weak and strong problems Consider the following simplified BVP and the corresponding weak form: ( qxy (, )) = f( xy, ),( xy, ) where : q = k us u, () () k s q n = q on Γ q, and u = u on Γu, with Γq Γu Γ( ) n The weak form is: Find uxy H w H withw on (, ) ( ): ( ) = 0 Γu : q wdxdy + fwdxdy = qwdγ Γ q We want to show that this weak form is equilavent to the strong problem (recover the PDE + natural BCs) 24
25 Equivalence between the weak and strong problems q wdxdy + fwdxdy = qwdγ Γ Integration by parts of the first term gives the following: q wq nd Γ w qd+ fwdxdy = qwdγ Γ Γ wq ndγ+ wq ndγ w q f d qwdγ= q ( ) 0 Γ Γ Γ u q q w q n q dγ+ wq ndγ w q f d= 0 w H ( ) with w = 0on Γ ( ) ( ) Γ Γ u q u 25
26 Equivalence between the weak and strong problems q wq n qdγ+ wq ndγ w q f d= 0 w H( ) withw= 0onΓ ( ) ( ) Γ Γ u Deriving the PDE first: Select w as follows: w = ψ( x) q f, withψ = 0on Γ and ψ > 0on ( ) The first equation above then results in: ψ q f d= q f = in ( ) u 26
27 Equivalence between the weak and strong problems Γ q wq n qdγ+ wq ndγ= 0 w H( ) withw= 0onΓ ( ) Γ u We next derive the natural BC: Select w as follows: w = ζ( x) q n q, with ζ = 0on Γ and ζ > 0on Γ ( ) The first equation above then results in: Γ q ζ q n q dγ= q n q= onγ ( ) u q u q 27
28 Generalized Fourier s law In general, the flux and temperature gradient are related with the anisotropic (generalized) Fourier s law: u q qx kxx kxy x [ D] u q = y kyx k = yy u 2x2conductivity matrix 2x y For isotropic materials k 0 [ D] = = k [ I] k 0 2 x 2 identity matrix For a 2D BVP with prescribed flux on Γ q and temperature on Γ u and a source term f, can you show that the weak form in matrix notation is: Find uxy H w H withw on (, ) ( ): ( ) 0 u : = Γ ( ) T w [ D] udxdy = qwdγ+ fwdxdy 2x2 2x Γq x2 28
29 Appendix: Green s theorem 0 If uxy (, ) C int egrable, then: ud= und Γ Γ Recall the one-dimensional version of this (where the boundary is only 2 points): du dx = un Γ = u ( x = L ) u ( x = 0) dx 29
30 Appendix: Divergence theorem 0 If qxy (, ) C int egrable, then: qd= q nd Γ Γ Note that: q x q y q = + x y Explicilty, we can write the divergence theorem as: q x x q d ( qxnx qyny) d Γ y + = + Γ y 30
31 Appendix: Green s theorem This is the analogue of integration by parts in D: w qd= wq ndγ w qd Γ The proof of this is simple if you notice that: w qx w qy ( wq) = ( wqx) + ( wqy) = qx + w + qy + w x y x x y y ( wq) = w q + w q This together with the divergence theorem give: w qd= ( wq) d w qd= wq ndγ w qd Γ 3
32 Appendix: Green s theorem in D w qd= wq ndγ w qd Γ The D analogue of Green s theorem is as follows: q w = Γ x w d wqxnd qxd x x Γ or more explicitly as: q w w w dx= wqn qdx= wq wq qdx ( ) x x Γ x x x= L x x= 0 x x x x Recall that we had used this compact form in D to treat boundary conditions on the left (x=0) and the right (x=l) end points in a unified way 32
33 Example problem on divergence theorem Consider the vector function q= ( qx. qy) = (3 x y+ y,3 x + y ), defined in the domain shown in the figure below. Demonstrate the validity of the divergence theorem, i.e. qd= q nd Γ Γ 33
34 Example problem on divergence theorem From q= ( qx. qy) = (3 x y+ y,3 x + y ), we compute the following: q x q y Using this, we compute qd as follows: 2 2 ( 6 ) (3 ) 6 3 q= + = xy + y = xy+ y x y 2 0.5x qd= (6xy + 3 y ) dydx = [3 x( 0.5 x) + ( 0.5 x) ] dx = y = 0.5x dx 34
35 Example problem on divergence theorem Also: Γ q ndγ= q n d Γ+ q n d Γ + q n d Γ AB (0, ) dx BC 5 5 CA (,0) dy (,2) dx (2) 5 Note that on segment BC: ds = dx and n = (, 2). Thus Γ q ndγ= (3 x+ y ) dx+ (3 x y + y ) + 2(3 x+ y ) ( dx) + ( 3 x y y )( dy) x dx q ndγ= 6 + (3 ( 0.5 ) ( 0.5 ) ) 2(3 ( 0.5 ) ) x x + x + x+ x dx Γ 0 q nd Γ= Γ y = 0.5x 2 2 q ndγ= =.5 ds = dx + dy = dx 2 2 dx 5 + = dx 2 2 Γ qd= q nd Γ Γ 35
36 Example: Deriving the weak form Consider isotropic heat conduction in the domain shown with the boundary conditions indicated. Provide a complete weak statement of the problem. ( kxy (, ) T( x, y)) = f( x, y) We start by multiplying with the weight function w and integrating over the domain: ( kxy (, ) T( x, y)) f( x, y) wxy (, ) dxdy = 0, sufficiently smooth w 36
37 Example: Deriving the weak form ( kxy (, ) T( x, y)) f( x, y) wxy (, ) dxdy = 0, sufficiently smooth w Integration by parts gives: T k T w fw dxdy k wds = 0, w H ( ), with w = 0 on Γ n Γ +Γ q h T The boundary term can be further simplified as: T k wds = ( q) wds + ( h( T T )) wds = qwds htwds + ht wds n Γ +Γ Γ Γ Γ Γ Γ q h q h q h h 37
38 Example: Deriving the weak form u k T w fw dxdy k wds = 0, w H ( ), withw= 0 onγ n Γ +Γ The final weak form is then: q h Find T ( x, y) H ( ), with T = T on ΓT such that + = + = Γ k T w fw dxdy htwds qwds ht wds, w H ( ), with w 0on Γ Γ Γ h q h In matrix form, we can also write the weak form as: T ( w) k T fw dxdy 2 x x x2 + htwds = qwds + ht wds Γ Γ Γ h q h w H ( ), with w = 0on Γ T T T 38
MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design
MAE4700/5700 Finite Element Analsis for Mechanical and Aerospace Design Cornell Universit, Fall 2009 Nicholas Zabaras Materials Process Design and Control Laborator Sible School of Mechanical and Aerospace
More informationMATH 332: Vector Analysis Summer 2005 Homework
MATH 332, (Vector Analysis), Summer 2005: Homework 1 Instructor: Ivan Avramidi MATH 332: Vector Analysis Summer 2005 Homework Set 1. (Scalar Product, Equation of a Plane, Vector Product) Sections: 1.9,
More informationPRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.
PRACTICE PROBLEMS Please let me know if you find any mistakes in the text so that i can fix them. 1.1. Let Show that f is C 1 and yet How is that possible? 1. Mixed partial derivatives f(x, y) = {xy x
More informationSummary of various integrals
ummary of various integrals Here s an arbitrary compilation of information about integrals Moisés made on a cold ecember night. 1 General things o not mix scalars and vectors! In particular ome integrals
More informationMath 234 Exam 3 Review Sheet
Math 234 Exam 3 Review Sheet Jim Brunner LIST OF TOPIS TO KNOW Vector Fields lairaut s Theorem & onservative Vector Fields url Divergence Area & Volume Integrals Using oordinate Transforms hanging the
More informationMathematics (Course B) Lent Term 2005 Examples Sheet 2
N12d Natural Sciences, Part IA Dr M. G. Worster Mathematics (Course B) Lent Term 2005 Examples Sheet 2 Please communicate any errors in this sheet to Dr Worster at M.G.Worster@damtp.cam.ac.uk. Note that
More informationJim Lambers MAT 280 Summer Semester Practice Final Exam Solution. dy + xz dz = x(t)y(t) dt. t 3 (4t 3 ) + e t2 (2t) + t 7 (3t 2 ) dt
Jim Lambers MAT 28 ummer emester 212-1 Practice Final Exam olution 1. Evaluate the line integral xy dx + e y dy + xz dz, where is given by r(t) t 4, t 2, t, t 1. olution From r (t) 4t, 2t, t 2, we obtain
More informationVector Calculus, Maths II
Section A Vector Calculus, Maths II REVISION (VECTORS) 1. Position vector of a point P(x, y, z) is given as + y and its magnitude by 2. The scalar components of a vector are its direction ratios, and represent
More informationPractice Problems for Exam 3 (Solutions) 1. Let F(x, y) = xyi+(y 3x)j, and let C be the curve r(t) = ti+(3t t 2 )j for 0 t 2. Compute F dr.
1. Let F(x, y) xyi+(y 3x)j, and let be the curve r(t) ti+(3t t 2 )j for t 2. ompute F dr. Solution. F dr b a 2 2 F(r(t)) r (t) dt t(3t t 2 ), 3t t 2 3t 1, 3 2t dt t 3 dt 1 2 4 t4 4. 2. Evaluate the line
More informationIntroduction and Vectors Lecture 1
1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum
More informationMath 234 Final Exam (with answers) Spring 2017
Math 234 Final Exam (with answers) pring 217 1. onsider the points A = (1, 2, 3), B = (1, 2, 2), and = (2, 1, 4). (a) [6 points] Find the area of the triangle formed by A, B, and. olution: One way to solve
More informationMathematical Notation Math Calculus & Analytic Geometry III
Name : Mathematical Notation Math 221 - alculus & Analytic Geometry III Use Word or WordPerect to recreate the ollowing documents. Each article is worth 10 points and can e printed and given to the instructor
More informationAP Calculus BC - Problem Solving Drill 19: Parametric Functions and Polar Functions
AP Calculus BC - Problem Solving Drill 19: Parametric Functions and Polar Functions Question No. 1 of 10 Instructions: (1) Read the problem and answer choices carefully () Work the problems on paper as
More informationMathematical Concepts & Notation
Mathematical Concepts & Notation Appendix A: Notation x, δx: a small change in x t : the partial derivative with respect to t holding the other variables fixed d : the time derivative of a quantity that
More informationMath 212-Lecture 20. P dx + Qdy = (Q x P y )da. C
15. Green s theorem Math 212-Lecture 2 A simple closed curve in plane is one curve, r(t) : t [a, b] such that r(a) = r(b), and there are no other intersections. The positive orientation is counterclockwise.
More informationCorrections to the First Printing: Chapter 6. This list includes corrections and clarifications through November 6, 1999.
June 2,2 1 Corrections to the First Printing: Chapter 6 This list includes corrections and clarifications through November 6, 1999. We should have mentioned that our treatment of differential forms, especially
More informationCreated by T. Madas LINE INTEGRALS. Created by T. Madas
LINE INTEGRALS LINE INTEGRALS IN 2 DIMENSIONAL CARTESIAN COORDINATES Question 1 Evaluate the integral ( x + 2y) dx, C where C is the path along the curve with equation y 2 = x + 1, from ( ) 0,1 to ( )
More information18.02 Multivariable Calculus Fall 2007
MIT OpenourseWare http://ocw.mit.edu 8.02 Multivariable alculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 8.02 Lecture 8. hange of variables.
More informationChapter 3 - Vector Calculus
Chapter 3 - Vector Calculus Gradient in Cartesian coordinate system f ( x, y, z,...) dr ( dx, dy, dz,...) Then, f f f f,,,... x y z f f f df dx dy dz... f dr x y z df 0 (constant f contour) f dr 0 or f
More informationExercise 1 (Formula for connection 1-forms) Using the first structure equation, show that
1 Stokes s Theorem Let D R 2 be a connected compact smooth domain, so that D is a smooth embedded circle. Given a smooth function f : D R, define fdx dy fdxdy, D where the left-hand side is the integral
More informationENGI 4430 Line Integrals; Green s Theorem Page 8.01
ENGI 4430 Line Integrals; Green s Theorem Page 8.01 8. Line Integrals Two applications of line integrals are treated here: the evaluation of work done on a particle as it travels along a curve in the presence
More informationMath 265H: Calculus III Practice Midterm II: Fall 2014
Name: Section #: Math 65H: alculus III Practice Midterm II: Fall 14 Instructions: This exam has 7 problems. The number of points awarded for each question is indicated in the problem. Answer each question
More informationMAC2313 Final A. (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative.
MAC2313 Final A (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative. ii. The vector field F = 5(x 2 + y 2 ) 3/2 x, y is radial. iii. All constant
More informationSections minutes. 5 to 10 problems, similar to homework problems. No calculators, no notes, no books, no phones. No green book needed.
MTH 34 Review for Exam 4 ections 16.1-16.8. 5 minutes. 5 to 1 problems, similar to homework problems. No calculators, no notes, no books, no phones. No green book needed. Review for Exam 4 (16.1) Line
More informationBTL What is the value of m if the vector is solenoidal. BTL What is the value of a, b, c if the vector may be irrotational.
VALLIAMMAI ENGINEERING OLLEGE SRM NAGAR, KATTANDKULATHUR Department of Mathematics MA65 - MATHEMATIS II QUESTION BANK - 6 UNIT - I VETOR ALULUS Part - A. Find, if at (, -, ). BTL-. Find the Directional
More informationMATH H53 : Final exam
MATH H53 : Final exam 11 May, 18 Name: You have 18 minutes to answer the questions. Use of calculators or any electronic items is not permitted. Answer the questions in the space provided. If you run out
More informationGradient Estimate of Mean Curvature Equations and Hessian Equations with Neumann Boundary Condition
of Mean Curvature Equations and Hessian Equations with Neumann Boundary Condition Xinan Ma NUS, Dec. 11, 2014 Four Kinds of Equations Laplace s equation: u = f(x); mean curvature equation: div( Du ) =
More information(b) Find the range of h(x, y) (5) Use the definition of continuity to explain whether or not the function f(x, y) is continuous at (0, 0)
eview Exam Math 43 Name Id ead each question carefully. Avoid simple mistakes. Put a box around the final answer to a question (use the back of the page if necessary). For full credit you must show your
More informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More informationDr. Allen Back. Dec. 3, 2014
Dr. Allen Back Dec. 3, 2014 forms are sums of wedge products of the basis 1-forms dx, dy, and dz. They are kinds of tensors generalizing ordinary scalar functions and vector fields. They have a skew-symmetry
More information29.3. Integral Vector Theorems. Introduction. Prerequisites. Learning Outcomes
Integral ector Theorems 9. Introduction arious theorems exist relating integrals involving vectors. Those involving line, surface and volume integrals are introduced here. They are the multivariable calculus
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationDivergence Theorem and Its Application in Characterizing
Divergence Theorem and Its Application in Characterizing Fluid Flow Let v be the velocity of flow of a fluid element and ρ(x, y, z, t) be the mass density of fluid at a point (x, y, z) at time t. Thus,
More informationName: Instructor: Lecture time: TA: Section time:
Math 222 Final May 11, 29 Name: Instructor: Lecture time: TA: Section time: INSTRUCTIONS READ THIS NOW This test has 1 problems on 16 pages worth a total of 2 points. Look over your test package right
More informationIdeas from Vector Calculus Kurt Bryan
Ideas from Vector Calculus Kurt Bryan Most of the facts I state below are for functions of two or three variables, but with noted exceptions all are true for functions of n variables..1 Tangent Line Approximation
More informationFinite-Elements Method 2
Finite-Elements Method 2 January 29, 2014 2 From Applied Numerical Analysis Gerald-Wheatley (2004), Chapter 9. Finite-Elements Method 3 Introduction Finite-element methods (FEM) are based on some mathematical
More informationGreen s, Divergence, Stokes: Statements and First Applications
Math 425 Notes 12: Green s, Divergence, tokes: tatements and First Applications The Theorems Theorem 1 (Divergence (planar version)). Let F be a vector field in the plane. Let be a nice region of the plane
More informationEELE 3331 Electromagnetic I Chapter 3. Vector Calculus. Islamic University of Gaza Electrical Engineering Department Dr.
EELE 3331 Electromagnetic I Chapter 3 Vector Calculus Islamic University of Gaza Electrical Engineering Department Dr. Talal Skaik 2012 1 Differential Length, Area, and Volume This chapter deals with integration
More informationx + ye z2 + ze y2, y + xe z2 + ze x2, z and where T is the
1.(8pts) Find F ds where F = x + ye z + ze y, y + xe z + ze x, z and where T is the T surface in the pictures. (The two pictures are two views of the same surface.) The boundary of T is the unit circle
More informationMath 31CH - Spring Final Exam
Math 3H - Spring 24 - Final Exam Problem. The parabolic cylinder y = x 2 (aligned along the z-axis) is cut by the planes y =, z = and z = y. Find the volume of the solid thus obtained. Solution:We calculate
More informationCalculus III 2004 Summer Practice Final 8/3/2004
.. Calculus III 4 ummer Practice Final 8/3/4. Compute the following limits if they exist: (a) lim (x,y) (,) e xy x+. cos x (b) lim x. (x,y) (,) x 4 +y 4 (a) ince lim (x,y) (,) exy and lim x + 6 in a (x,y)
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationMAT 211 Final Exam. Spring Jennings. Show your work!
MAT 211 Final Exam. pring 215. Jennings. how your work! Hessian D = f xx f yy (f xy ) 2 (for optimization). Polar coordinates x = r cos(θ), y = r sin(θ), da = r dr dθ. ylindrical coordinates x = r cos(θ),
More informationTwo-dimensional flow in a porous medium with general anisotropy
Two-dimensional flow in a porous medium with general anisotropy P.A. Tyvand & A.R.F. Storhaug Norwegian University of Life Sciences 143 Ås Norway peder.tyvand@umb.no 1 Darcy s law for flow in an isotropic
More information1 + f 2 x + f 2 y dy dx, where f(x, y) = 2 + 3x + 4y, is
1. The value of the double integral (a) 15 26 (b) 15 8 (c) 75 (d) 105 26 5 4 0 1 1 + f 2 x + f 2 y dy dx, where f(x, y) = 2 + 3x + 4y, is 2. What is the value of the double integral interchange the order
More informationLecture 10 Divergence, Gauss Law in Differential Form
Lecture 10 Divergence, Gauss Law in Differential Form ections: 3.4, 3.5, 3.6 Homework: ee homework file Properties of the Flux Integral: Recap flux is the net normal flow of the vector field F through
More informationMathematical Notation Math Calculus & Analytic Geometry III
Mathematical Notation Math 221 - alculus & Analytic Geometry III Use Word or WordPerect to recreate the ollowing documents. Each article is worth 10 points and should be emailed to the instructor at james@richland.edu.
More information4 Divergence theorem and its consequences
Tel Aviv University, 205/6 Analysis-IV 65 4 Divergence theorem and its consequences 4a Divergence and flux................. 65 4b Piecewise smooth case............... 67 4c Divergence of gradient: Laplacian........
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationVector and Tensor Calculus
Appendices 58 A Vector and Tensor Calculus In relativistic theory one often encounters vector and tensor expressions in both three- and four-dimensional form. The most important of these expressions are
More informationLecture 13 - Wednesday April 29th
Lecture 13 - Wednesday April 29th jacques@ucsdedu Key words: Systems of equations, Implicit differentiation Know how to do implicit differentiation, how to use implicit and inverse function theorems 131
More informationUniformity of the Universe
Outline Universe is homogenous and isotropic Spacetime metrics Friedmann-Walker-Robertson metric Number of numbers needed to specify a physical quantity. Energy-momentum tensor Energy-momentum tensor of
More informationSolution. This is a routine application of the chain rule.
EXAM 2 SOLUTIONS 1. If z = e r cos θ, r = st, θ = s 2 + t 2, find the partial derivatives dz ds chain rule. Write your answers entirely in terms of s and t. dz and dt using the Solution. This is a routine
More informationMATH 31CH SPRING 2017 MIDTERM 2 SOLUTIONS
MATH 3CH SPRING 207 MIDTERM 2 SOLUTIONS (20 pts). Let C be a smooth curve in R 2, in other words, a -dimensional manifold. Suppose that for each x C we choose a vector n( x) R 2 such that (i) 0 n( x) for
More informationMultivariable Calculus Lecture #12 Notes
Multivariable alculus Lecture #12 Notes In this lecture, we will develop a list of statements equivalent to a vector field being conservative and state and prove Green s heorem to help connect these facts.
More informationQ ( q(m, t 0 ) n) S t.
THE HEAT EQUATION The main equations that we will be dealing with are the heat equation, the wave equation, and the potential equation. We use simple physical principles to show how these equations are
More informationBrief Review of Vector Algebra
APPENDIX Brief Review of Vector Algebra A.0 Introduction Vector algebra is used extensively in computational mechanics. The student must thus understand the concepts associated with this subject. The current
More informationPractice problems **********************************************************
Practice problems I will not test spherical and cylindrical coordinates explicitly but these two coordinates can be used in the problems when you evaluate triple integrals. 1. Set up the integral without
More information1. (30 points) In the x-y plane, find and classify all local maxima, local minima, and saddle points of the function. f(x, y) = 3y 2 2y 3 3x 2 + 6xy.
APPM 35 FINAL EXAM FALL 13 INSTUTIONS: Electronic devices, books, and crib sheets are not permitted. Write your name and your instructor s name on the front of your bluebook. Work all problems. Show your
More informationMath 67. Rumbos Fall Solutions to Review Problems for Final Exam. (a) Use the triangle inequality to derive the inequality
Math 67. umbos Fall 8 Solutions to eview Problems for Final Exam. In this problem, u and v denote vectors in n. (a) Use the triangle inequality to derive the inequality Solution: Write v u v u for all
More informationC 3 C 4. R k C 1. (x,y)
16.4 1 16.4 Green s Theorem irculation Density (x,y + y) 3 (x+ x,y + y) 4 k 2 (x,y) 1 (x+ x,y) Suppose that F(x,y) M(x,y)i+N(x,y)j is the velocity field of a fluid flow in the plane and that the first
More informationES.182A Topic 44 Notes Jeremy Orloff
E.182A Topic 44 Notes Jeremy Orloff 44 urface integrals and flux Note: Much of these notes are taken directly from the upplementary Notes V8, V9 by Arthur Mattuck. urface integrals are another natural
More informationMultivariable Calculus Midterm 2 Solutions John Ross
Multivariable Calculus Midterm Solutions John Ross Problem.: False. The double integral is not the same as the iterated integral. In particular, we have shown in a HW problem (section 5., number 9) that
More informationSecond Order ODEs. Second Order ODEs. In general second order ODEs contain terms involving y, dy But here only consider equations of the form
Second Order ODEs Second Order ODEs In general second order ODEs contain terms involving y, dy But here only consider equations of the form A d2 y dx 2 + B dy dx + Cy = 0 dx, d2 y dx 2 and F(x). where
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
More informationChapter 2. First-Order Partial Differential Equations. Prof. D. C. Sanyal
BY Prof. D. C. Sanyal Retired Professor Of Mathematics University Of Kalyani West Bengal, India E-mail : dcs klyuniv@yahoo.com 1 Module-2: Quasi-Linear Equations of First Order 1. Introduction In this
More information18.1. Math 1920 November 29, ) Solution: In this function P = x 2 y and Q = 0, therefore Q. Converting to polar coordinates, this gives I =
Homework 1 elected olutions Math 19 November 9, 18 18.1 5) olution: In this function P = x y and Q =, therefore Q x P = x. We obtain the following integral: ( Q I = x ydx = x P ) da = x da. onverting to
More informationMIT 2.71/2.710 Optics 10/31/05 wk9-a-1. The spatial frequency domain
10/31/05 wk9-a-1 The spatial frequency domain Recall: plane wave propagation x path delay increases linearly with x λ z=0 θ E 0 x exp i2π sinθ + λ z i2π cosθ λ z plane of observation 10/31/05 wk9-a-2 Spatial
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationMath 265 (Butler) Practice Midterm III B (Solutions)
Math 265 (Butler) Practice Midterm III B (Solutions). Set up (but do not evaluate) an integral for the surface area of the surface f(x, y) x 2 y y over the region x, y 4. We have that the surface are is
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, < 1, 1 +, 1 < < 0, ψ() = 1, 0 < < 1, 0, > 1, so that it verifies ψ 0, ψ() = 0 if 1 and ψ()d = 1. Consider (ψ j ) j 1 constructed as
More informationFinal Exam. Monday March 19, 3:30-5:30pm MAT 21D, Temple, Winter 2018
Name: Student ID#: Section: Final Exam Monday March 19, 3:30-5:30pm MAT 21D, Temple, Winter 2018 Show your work on every problem. orrect answers with no supporting work will not receive full credit. Be
More informationMultiple Integrals and Vector Calculus: Synopsis
Multiple Integrals and Vector Calculus: Synopsis Hilary Term 28: 14 lectures. Steve Rawlings. 1. Vectors - recap of basic principles. Things which are (and are not) vectors. Differentiation and integration
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationDO NOT BEGIN THIS TEST UNTIL INSTRUCTED TO START
Math 265 Student name: KEY Final Exam Fall 23 Instructor & Section: This test is closed book and closed notes. A (graphing) calculator is allowed for this test but cannot also be a communication device
More informationD. Correct! This is the correct answer. It is found by dy/dx = (dy/dt)/(dx/dt).
Calculus II - Problem Solving Drill 4: Calculus for Parametric Equations Question No. of 0 Instructions: () Read the problem and answer choices carefully () Work the problems on paper as. Find dy/dx where
More informationSolutions to old Exam 3 problems
Solutions to old Exam 3 problems Hi students! I am putting this version of my review for the Final exam review here on the web site, place and time to be announced. Enjoy!! Best, Bill Meeks PS. There are
More information6. LIGHT SCATTERING 6.1 The first Born approximation
6. LIGHT SCATTERING 6.1 The first Born approximation In many situations, light interacts with inhomogeneous systems, in which case the generic light-matter interaction process is referred to as scattering
More information17th European Conference on Fracture 2-5 September,2008, Brno, Czech Republic. Thermal Fracture of a FGM/Homogeneous Bimaterial with Defects
-5 September,8, Brno, Czech Republic Thermal Fracture of a FGM/Homogeneous Bimaterial with Defects Vera Petrova, a, Siegfried Schmauder,b Voronezh State University, University Sq., Voronezh 3946, Russia
More information12. Stresses and Strains
12. Stresses and Strains Finite Element Method Differential Equation Weak Formulation Approximating Functions Weighted Residuals FEM - Formulation Classification of Problems Scalar Vector 1-D T(x) u(x)
More informationJoint Distributions: Part Two 1
Joint Distributions: Part Two 1 STA 256: Fall 2018 1 This slide show is an open-source document. See last slide for copyright information. 1 / 30 Overview 1 Independence 2 Conditional Distributions 3 Transformations
More informationMATHS 267 Answers to Stokes Practice Dr. Jones
MATH 267 Answers to tokes Practice Dr. Jones 1. Calculate the flux F d where is the hemisphere x2 + y 2 + z 2 1, z > and F (xz + e y2, yz, z 2 + 1). Note: the surface is open (doesn t include any of the
More informationG G. G. x = u cos v, y = f(u), z = u sin v. H. x = u + v, y = v, z = u v. 1 + g 2 x + g 2 y du dv
1. Matching. Fill in the appropriate letter. 1. ds for a surface z = g(x, y) A. r u r v du dv 2. ds for a surface r(u, v) B. r u r v du dv 3. ds for any surface C. G x G z, G y G z, 1 4. Unit normal N
More informationUnit 6 Line and Surface Integrals
Unit 6 Line and Surface Integrals In this unit, we consider line integrals and surface integrals and the relationships between them. We also discuss the three theorems Green s theorem, the divergence theorem
More informationReview problems for the final exam Calculus III Fall 2003
Review problems for the final exam alculus III Fall 2003 1. Perform the operations indicated with F (t) = 2t ı 5 j + t 2 k, G(t) = (1 t) ı + 1 t k, H(t) = sin(t) ı + e t j a) F (t) G(t) b) F (t) [ H(t)
More informationWORKSHEET #13 MATH 1260 FALL 2014
WORKSHEET #3 MATH 26 FALL 24 NOT DUE. Short answer: (a) Find the equation of the tangent plane to z = x 2 + y 2 at the point,, 2. z x (, ) = 2x = 2, z y (, ) = 2y = 2. So then the tangent plane equation
More informationf(p i )Area(T i ) F ( r(u, w) ) (r u r w ) da
MAH 55 Flux integrals Fall 16 1. Review 1.1. Surface integrals. Let be a surface in R. Let f : R be a function defined on. efine f ds = f(p i Area( i lim mesh(p as a limit of Riemann sums over sampled-partitions.
More informationMATH 311 Topics in Applied Mathematics I Lecture 36: Surface integrals (continued). Gauss theorem. Stokes theorem.
MATH 311 Topics in Applied Mathematics I Lecture 36: Surface integrals (continued). Gauss theorem. Stokes theorem. Surface integrals Let X : R 3 be a smooth parametrized surface, where R 2 is a bounded
More informationMATH 311 Topics in Applied Mathematics I Lecture 34: Conservative vector fields. Area of a surface. Surface integrals.
MATH 311 Topics in Applied Mathematics I Lecture 34: Conservative vector fields. Area of a surface. Surface integrals. Conservative vector fields Let R be an open region in R n such that any two points
More informationDiff. Eq. App.( ) Midterm 1 Solutions
Diff. Eq. App.(110.302) Midterm 1 Solutions Johns Hopkins University February 28, 2011 Problem 1.[3 15 = 45 points] Solve the following differential equations. (Hint: Identify the types of the equations
More informationGREEN S IDENTITIES AND GREEN S FUNCTIONS
GREEN S IENTITIES AN GREEN S FUNCTIONS Green s first identity First, recall the following theorem. Theorem: (ivergence Theorem) Let be a bounded solid region with a piecewise C 1 boundary surface. Let
More informationMath 220A - Fall 2002 Homework 5 Solutions
Math 0A - Fall 00 Homework 5 Solutions. Consider the initial-value problem for the hyperbolic equation u tt + u xt 0u xx 0 < x 0 u t (x, 0) ψ(x). Use energy methods to show that the domain of dependence
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More information16.2. Line Integrals
16. Line Integrals Review of line integrals: Work integral Rules: Fdr F d r = Mdx Ndy Pdz FT r'( t) ds r t since d '(s) and hence d ds '( ) r T r r ds T = Fr '( t) dt since r r'( ) dr d dt t dt dt does
More informationAn Introduction to Partial Differential Equations
An Introduction to Partial Differential Equations Ryan C. Trinity University Partial Differential Equations Lecture 1 Ordinary differential equations (ODEs) These are equations of the form where: F(x,y,y,y,y,...)
More informationSolutions to Sample Questions for Final Exam
olutions to ample Questions for Final Exam Find the points on the surface xy z 3 that are closest to the origin. We use the method of Lagrange Multipliers, with f(x, y, z) x + y + z for the square of the
More informationReview Sheet for the Final
Review Sheet for the Final Math 6-4 4 These problems are provided to help you study. The presence of a problem on this handout does not imply that there will be a similar problem on the test. And the absence
More informationGradient, Divergence and Curl in Curvilinear Coordinates
Gradient, Divergence and Curl in Curvilinear Coordinates Although cartesian orthogonal coordinates are very intuitive and easy to use, it is often found more convenient to work with other coordinate systems.
More information