Approximation of self-avoiding inextensible curves
|
|
- Sheena Simpson
- 5 years ago
- Views:
Transcription
1 Approximation of self-avoiding inextensible curves Sören Bartels Department of Applied Mathematics University of Freiburg, Germany Joint with Philipp Reiter (U Duisburg-Essen) & Johannes Riege (U Freiburg) 8th Workshop on Numerical Methods for Evolution Equations Heraklion, Crete, September 23-24, / 15
2 Motivation Energy for large bending deformations of thin elastic sheets: I bend (y) = 1 II 2 dx f y dx, I = id 2 ω with fundamental forms I, II : ω R 2 2 associated with y : ω R 3 flat isometries: angle/length relations preserved nonlinear constraint implies linearity: II 2 = D 2 y 2 Kirchhoff s bending model derived by Friesecke, James & Müller 02 Nanotechnology simulation (B., Bonito & Nochetto 15): bilayer bending ω Problem: Large deformations ok but injectivity not guaranteed 2 / 15
3 Curves Simpler setting: inextensible curves describing elastic rods I bend (u) = 1 u (x) 2 dx, u (x) 2 = 1 2 curves parametrized by arc-length well posed minimization problems on H 2 (I; R 3 ) derivation from hyperelasticity: Müller & Mora 04 I I u EL eq s: stationary u H 2 (I; R 3 ) with u (x) 2 = 1 characterized by (u, v ) = u (x) v (x) dx = 0 for admissible v H 2 (I; R 3 ) with u v = 0 in I. I 3 / 15
4 Iteration B 12: Practical minimization via gradient flow ( t u, v) + (u, v ) = 0, u(0) = u 0, u (t, x) 2 = 1 for admissible v H 2 (I; R 3 ) with v u (t, ) = 0 Observation: relation t u u = 0 implies constraint preservation, u (t) 2 u (0) 2 = t 0 d ds u (s) 2 ds = 2 t 0 t u u ds = 0 Algorithm: given u n 1 compute d t u n with [d t u n ] [u n 1 ] = 0 and (d t u n, v) + ([u n 1 + τd t u n ], v ) = 0 for v with v [u n 1 ] = 0 and update u n = u n 1 + τd t u n Unconditional stability and controlled constraint violation: n [u n ] 2 = [u n 1 ] 2 +τ 2 [d t u n ] 2 =... = [u 0 ] 2 +τ 2 [d t u l ] 2 = 1+O(τ) l=1 4 / 15
5 Discretization Use H 2 conforming piecewise cubic C 1 curves DOFs: positions and tangents good approximation properties flexibility for boundary conditions u h (z) u (z) h Impose constraint/orthogonality at nodes z, i.e., u h(z) 2 = 1, d t u h(z) u h(z) = 0 constraint satisfied in limit h 0 if H 2 bounds available 5 / 15
6 Self-avoidance Maddocks & Gonzales 99: Injectivity enforced via potential TP(u) = 2 q 1 q r q( ) dx dy u(y), u(x) I I with r = radius of circle tangential at u(y) and intersecting in u(x) u(x) l(y) α r u(y) α u(x) u(y) u(y) r 0 as u(x) u(y) u(x) Using u (y) = 1 and sin α = u (y) [u(y) u(x)] / u(y) u(x) r ( u(y), u(x) ) = and 1/r u (x) /2 for x y u(y) u(x) 2 2 u (y) [u(y) u(x)] 6 / 15
7 Properties Tangent-point potential TP(u) = 1 q I I u (y) [u(y) u(x)] q u(y) u(x) 2q dx dy Sobolev Slobodeckii spaces defined via seminorm, 0 < s < 1, ( f (x) f (y) p ) 1/p [f ] W s,p (I;R 3 ) = x y 1+sp dx dy I I Class of embedded and arc-length parametrized curves: { } C q = u W 2 1/q,q (I; R 3 ) : u embedded and u (x) = 1 Basic properties (Blatt 13, Blatt & Reiter 15) for q > 2: u C q implies 1 bil(u) = Lip(u 1 ) < TP finite on C q and TP(u) M implies bil(u) C M,q TP is knot energy: TP(u k ) if pointw. lim u k C q 7 / 15
8 Variation Blatt & Reiter 15: TP is C 1 on embedded W 2 1/q,q curves with where M(u; v, w) = A(u; v, w) = with δtp(u) = M(u; u, ϕ) + M(u; ϕ, u) 2A(u; u, ϕ) I I I I Φ(x, y) (v (y) ( w(x) w(y) )) dx dy Ψ(x, y) ( v(x) v(y) ) (w(x) w(y) ) dx dy Φ(x, y) = u (y) ( u(x) u(y) ) q 2 u(x) u(y) 2q ( u (y) ( u(x) u(y) )), Ψ(x, y) = u (y) ( u(x) u(y) ) q u(x) u(y) 2q+2 Note: A(u;, ) positive semidefinite and M(u; u, u) 0 8 / 15
9 Approximation (A) Removal of singular diagonal u (y) [u(y) u(x)] q TP ε (u) = u(y) u(x) 2q dx dy leads to error TP(u) TP ε (u) (B) Use of quadrature TP ε,h (u) = x y ε { c δ,q ε 2 2/q δ bil(u) 2q u H 1 u H 2 x y ε c q ε bil(u) 2q u 2 L q u W 2,q Q h [ u (y) [u(y) u(x)] q u(y) u(x) 2q ] dx dy with node average on rectangles I j I k gives TPε (u) TP ε,h (u) { c q h 1/2 ε ( (q+1) u L ) u H 2 c q h ( u W 3, + 1 ) u W 3, 9 / 15
10 Proofs With S ε = {(x, y) I I : x y < ε} we have q ( TP(u) TP ε (u) ) u (y) ( u(y + z) u(y) zu (y) ) q S ε u(y + z) u(y) 2q dz dy 1 ( bil(u) 2q 0 u (y + ϑz) u (y) ) dϑ q z q S ε z 2q dz dy 1 u bil(u) 2q (y + ϑz) u (y) q S ε 0 z q dϑ dz dy 1 u bil(u) 2q ϑ q 1 (y + z) u (y) q 0 S ϑε z q dz dy dϑ u bil(u) 2q (y + z) u (y) q ( ) z q dz dy... For quadrature estimate show that integrand is in C 0,1/2 S ε 10 / 15
11 Gradient flow Weighted sum of bending energy and self-avoidance potential leads to gradient flow dynamics I tot (u) = I bend (u) + ϱtp ε (u) ( t u, v) = δi tot (u)[v] = B(u, v) 2ϱM sym ε (u; u, v) + 2ϱA ε (u; u, v) with 2M sym ε (u; v, w) = M ε (u; v, w) + M ε (u; w, v) Semi-implicit discretization (d t u n, v) + B(u n, v) + 2ϱM sym ε (u n 1 ; u n, v) = 2ϱA ε (u n 1 ; u n 1, v) constrained linear system of equations fully populated system matrices parameter ϱ defines length scale 11 / 15
12 Experiment I With I = [0, 2π] define planar circle cos(y) u(y) = sin(y). 0 Relative error quantities: δ h TP = TPε,h (I h u) TP(u) /TP(u), δ h M = Mε,h (I h u; I h u, I h u) M(u; u, u) /M(u; u, u), Parameters: q = 3, ε = 2h, uniform partitions h j δ j TP δ j M 2 1 / / / / / / 15
13 Experiment II With Ĩ = [0, 1] define trefoil knot (2 + cos(6πỹ)) cos(4πỹ) ũ(ỹ) = (2 + cos(6πỹ)) sin(4πỹ), sin(6πỹ) and discrete transformation gives arc-length parametrized u : I R 3 Movie TP eps,h (u h ) (M=20) -"- (M=40) -"- (M=80) I tot h (u h ) (M=20) -"- (M=40) -"- (M=80) Stationary states 1/h 20, 40, τ = h/4, ε = 2h, ϱ = h 1/ / 15
14 Experiment III With Ĩ = [0, 1] define planar spiral 20 cos(2π(2ỹ + 1))/(2π(2ỹ + 1)) u(ỹ) = 20 sin(2π(2ỹ + 1))/(2π(2ỹ + 1)) 0 and discrete transformation gives arc-length parametrized u : I R 3 Movie TP eps,h (u h ) (M=20) -"- (M=40) -"- (M=80) I tot h (u h ) (M=20) -"- (M=40) -"- (M=80) τ = h/4, ε = 2h, ϱ = h 1/ / 15
15 Summary Numerical scheme for self-avoiding elastic curves Error bounds for approximation and discretization Good experimental stability properties Future aspects Stability analysis for time-stepping Treatment of surfaces, efficiency More information 15 / 15
Existence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
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 informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationSample Final Questions: Solutions Math 21B, Winter y ( y 1)(1 + y)) = A y + B
Sample Final Questions: Solutions Math 2B, Winter 23. Evaluate the following integrals: tan a) y y dy; b) x dx; c) 3 x 2 + x dx. a) We use partial fractions: y y 3 = y y ) + y)) = A y + B y + C y +. Putting
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 informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationThere are some trigonometric identities given on the last page.
MA 114 Calculus II Fall 2015 Exam 4 December 15, 2015 Name: Section: Last 4 digits of student ID #: No books or notes may be used. Turn off all your electronic devices and do not wear ear-plugs during
More informationEnergy scaling law for a single disclination in a thin elastic shee
Energy scaling law for a single disclination in a thin elastic sheet 7 December, 2015 Overview 1 Introduction: Energy focusing in thin elastic sheets 2 3 4 Energy focusing in thin elastic sheets Experimental
More informationName: SOLUTIONS Date: 11/9/2017. M20550 Calculus III Tutorial Worksheet 8
Name: SOLUTIONS Date: /9/7 M55 alculus III Tutorial Worksheet 8. ompute R da where R is the region bounded by x + xy + y 8 using the change of variables given by x u + v and y v. Solution: We know R is
More informationMATH 280 Multivariate Calculus Fall Integrating a vector field over a curve
MATH 280 Multivariate alculus Fall 2012 Definition Integrating a vector field over a curve We are given a vector field F and an oriented curve in the domain of F as shown in the figure on the left below.
More informationFact Sheet Functional Analysis
Fact Sheet Functional Analysis Literature: Hackbusch, W.: Theorie und Numerik elliptischer Differentialgleichungen. Teubner, 986. Knabner, P., Angermann, L.: Numerik partieller Differentialgleichungen.
More information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informatione x2 dxdy, e x2 da, e x2 x 3 dx = e
STS26-4 Calculus II: The fourth exam Dec 15, 214 Please show all your work! Answers without supporting work will be not given credit. Write answers in spaces provided. You have 1 hour and 2minutes to complete
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 informationMultiple Choice. Compute the Jacobian, (u, v), of the coordinate transformation x = u2 v 4, y = uv. (a) 2u 2 + 4v 4 (b) xu yv (c) 3u 2 + 7v 6
.(5pts) y = uv. ompute the Jacobian, Multiple hoice (x, y) (u, v), of the coordinate transformation x = u v 4, (a) u + 4v 4 (b) xu yv (c) u + 7v 6 (d) u (e) u v uv 4 Solution. u v 4v u = u + 4v 4..(5pts)
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 informationThere are five problems. Solve four of the five problems. Each problem is worth 25 points. A sheet of convenient formulae is provided.
Preliminary Examination (Solutions): Partial Differential Equations, 1 AM - 1 PM, Jan. 18, 16, oom Discovery Learning Center (DLC) Bechtel Collaboratory. Student ID: There are five problems. Solve four
More informationMATH 280 Multivariate Calculus Fall Integration over a curve
dr dr y MATH 28 Multivariate Calculus Fall 211 Integration over a curve Given a curve C in the plane or in space, we can (conceptually) break it into small pieces each of which has a length ds. In some
More informationMath 210, Final Exam, Spring 2012 Problem 1 Solution. (a) Find an equation of the plane passing through the tips of u, v, and w.
Math, Final Exam, Spring Problem Solution. Consider three position vectors (tails are the origin): u,, v 4,, w,, (a) Find an equation of the plane passing through the tips of u, v, and w. (b) Find an equation
More information16.2 Line Integrals. Lukas Geyer. M273, Fall Montana State University. Lukas Geyer (MSU) 16.2 Line Integrals M273, Fall / 21
16.2 Line Integrals Lukas Geyer Montana State University M273, Fall 211 Lukas Geyer (MSU) 16.2 Line Integrals M273, Fall 211 1 / 21 Scalar Line Integrals Definition f (x) ds = lim { s i } N f (P i ) s
More informationENGI Parametric Vector Functions Page 5-01
ENGI 3425 5. Parametric Vector Functions Page 5-01 5. Parametric Vector Functions Contents: 5.1 Arc Length (Cartesian parametric and plane polar) 5.2 Surfaces of Revolution 5.3 Area under a Parametric
More informationMcGill University April 20, Advanced Calculus for Engineers
McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student
More informationGeometry and Motion Selected answers to Sections A and C Dwight Barkley 2016
MA34 Geometry and Motion Selected answers to Sections A and C Dwight Barkley 26 Example Sheet d n+ = d n cot θ n r θ n r = Θθ n i. 2. 3. 4. Possible answers include: and with opposite orientation: 5..
More informationStokes s Theorem 17.2
Stokes s Theorem 17.2 6 December 213 Stokes s Theorem is the generalization of Green s Theorem to surfaces not just flat surfaces (regions in R 2 ). Relate a double integral over a surface with a line
More informationHyperbolic Geometry on Geometric Surfaces
Mathematics Seminar, 15 September 2010 Outline Introduction Hyperbolic geometry Abstract surfaces The hemisphere model as a geometric surface The Poincaré disk model as a geometric surface Conclusion Introduction
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More information(a) The points (3, 1, 2) and ( 1, 3, 4) are the endpoints of a diameter of a sphere.
MATH 4 FINAL EXAM REVIEW QUESTIONS Problem. a) The points,, ) and,, 4) are the endpoints of a diameter of a sphere. i) Determine the center and radius of the sphere. ii) Find an equation for the sphere.
More informationSharp Korn inequalities in thin domains: The rst and a half Korn inequality
Sarp Korn inequalities in tin domains: Te rst and a alf Korn inequality Davit Harutyunyan (University of Uta) joint wit Yury Grabovsky (Temple University) SIAM, Analysis of Partial Dierential Equations,
More informationA distance space on Cayley s ruled surface
A distance space on Cayley s ruled surface Hans Havlicek Joint work with Johannes Gmainer DIFFERENTIALGEOMETRIE UND GEOMETRISCHE STRUKTUREN Basic notions Let P 3 (K) be the 3-dimensional projective space
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 information7a3 2. (c) πa 3 (d) πa 3 (e) πa3
1.(6pts) Find the integral x, y, z d S where H is the part of the upper hemisphere of H x 2 + y 2 + z 2 = a 2 above the plane z = a and the normal points up. ( 2 π ) Useful Facts: cos = 1 and ds = ±a sin
More information1. If the line l has symmetric equations. = y 3 = z+2 find a vector equation for the line l that contains the point (2, 1, 3) and is parallel to l.
. If the line l has symmetric equations MA 6 PRACTICE PROBLEMS x = y = z+ 7, find a vector equation for the line l that contains the point (,, ) and is parallel to l. r = ( + t) i t j + ( + 7t) k B. r
More informationMath 127C, Spring 2006 Final Exam Solutions. x 2 ), g(y 1, y 2 ) = ( y 1 y 2, y1 2 + y2) 2. (g f) (0) = g (f(0))f (0).
Math 27C, Spring 26 Final Exam Solutions. Define f : R 2 R 2 and g : R 2 R 2 by f(x, x 2 (sin x 2 x, e x x 2, g(y, y 2 ( y y 2, y 2 + y2 2. Use the chain rule to compute the matrix of (g f (,. By the chain
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationMath 11 Fall 2018 Practice Final Exam
Math 11 Fall 218 Practice Final Exam Disclaimer: This practice exam should give you an idea of the sort of questions we may ask on the actual exam. Since the practice exam (like the real exam) is not long
More informationIndian Institute of Technology Bombay
Indian Institute of Technology Bombay MA 5 Calculus Autumn SRG/AA/VDS/KSK Solutions and Marking Scheme for the Mid-Sem Exam Date: September, Weightage: 3 % Time: 8.3 AM.3 AM Max. Marks:. (i) Consider the
More informationGG303 Lecture 15 10/6/09 1 FINITE STRAIN AND INFINITESIMAL STRAIN
GG303 Lecture 5 0609 FINITE STRAIN AND INFINITESIMAL STRAIN I Main Topics on infinitesimal strain A The finite strain tensor [E] B Deformation paths for finite strain C Infinitesimal strain and the infinitesimal
More informationA Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term
A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term Peter Sternberg In collaboration with Dmitry Golovaty (Akron) and Raghav Venkatraman (Indiana) Department of Mathematics
More informationQuestion. [The volume of a cone of radius r and height h is 1 3 πr2 h and the curved surface area is πrl where l is the slant height of the cone.
Q1 An experiment is conducted using the conical filter which is held with its axis vertical as shown. The filter has a radius of 10cm and semi-vertical angle 30. Chemical solution flows from the filter
More informationBilayer Plates: From Model Reduction to Γ-Convergent Finite Element Approximations and More
Bilayer Plates: From Model Reduction to Γ-Convergent Finite Element Approximations and More Department of Mathematics Texas A&M University Joint work with: Soeren Bartels, University of Freiburg (Germany)
More informationII. DIFFERENTIABLE MANIFOLDS. Washington Mio CENTER FOR APPLIED VISION AND IMAGING SCIENCES
II. DIFFERENTIABLE MANIFOLDS Washington Mio Anuj Srivastava and Xiuwen Liu (Illustrations by D. Badlyans) CENTER FOR APPLIED VISION AND IMAGING SCIENCES Florida State University WHY MANIFOLDS? Non-linearity
More informationRule ST1 (Symmetry). α β = β α for 1-forms α and β. Like the exterior product, the symmetric tensor product is also linear in each slot :
2. Metrics as Symmetric Tensors So far we have studied exterior products of 1-forms, which obey the rule called skew symmetry: α β = β α. There is another operation for forming something called the symmetric
More informationExact Solutions of the Einstein Equations
Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.
More informationProblem Solving 1: Line Integrals and Surface Integrals
A. Line Integrals MASSACHUSETTS INSTITUTE OF TECHNOLOY Department of Physics Problem Solving 1: Line Integrals and Surface Integrals The line integral of a scalar function f ( xyz),, along a path C is
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 information49. Green s Theorem. The following table will help you plan your calculation accordingly. C is a simple closed loop 0 Use Green s Theorem
49. Green s Theorem Let F(x, y) = M(x, y), N(x, y) be a vector field in, and suppose is a path that starts and ends at the same point such that it does not cross itself. Such a path is called a simple
More informationCS 468. Differential Geometry for Computer Science. Lecture 17 Surface Deformation
CS 468 Differential Geometry for Computer Science Lecture 17 Surface Deformation Outline Fundamental theorem of surface geometry. Some terminology: embeddings, isometries, deformations. Curvature flows
More informationMATH 2083 FINAL EXAM REVIEW The final exam will be on Wednesday, May 4 from 10:00am-12:00pm.
MATH 2083 FINAL EXAM REVIEW The final exam will be on Wednesday, May 4 from 10:00am-12:00pm. Bring a calculator and something to write with. Also, you will be allowed to bring in one 8.5 11 sheet of paper
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationDG APPROACH TO LARGE BENDING DEFORMATIONS WITH ISOMETRY CONSTRAINT
DG APPROACH TO LARGE BENDING DEFORMATIONS WITH ISOMETRY CONSTRAINT Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA joint with Soeren Bartels, University
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
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 informationLent 2019 VECTOR CALCULUS EXAMPLE SHEET 1 G. Taylor
Lent 29 ECTOR CALCULU EXAMPLE HEET G. Taylor. (a) The curve defined parametrically by x(t) = (a cos 3 t, a sin 3 t) with t 2π is called an astroid. ketch it, and find its length. (b) The curve defined
More information14.1. Multiple Integration. Iterated Integrals and Area in the Plane. Iterated Integrals. Iterated Integrals. MAC2313 Calculus III - Chapter 14
14 Multiple Integration 14.1 Iterated Integrals and Area in the Plane Objectives Evaluate an iterated integral. Use an iterated integral to find the area of a plane region. Copyright Cengage Learning.
More informationMidterm Solution
18303 Midterm Solution Problem 1: SLP with mixed boundary conditions Consider the following regular) Sturm-Liouville eigenvalue problem consisting in finding scalars λ and functions v : [0, b] R b > 0),
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationPrint Your Name: Your Section:
Print Your Name: Your Section: Mathematics 1c. Practice Final Solutions This exam has ten questions. J. Marsden You may take four hours; there is no credit for overtime work No aids (including notes, books,
More informationMath 350 Solutions for Final Exam Page 1. Problem 1. (10 points) (a) Compute the line integral. F ds C. z dx + y dy + x dz C
Math 35 Solutions for Final Exam Page Problem. ( points) (a) ompute the line integral F ds for the path c(t) = (t 2, t 3, t) with t and the vector field F (x, y, z) = xi + zj + xk. (b) ompute the line
More informationMA 441 Advanced Engineering Mathematics I Assignments - Spring 2014
MA 441 Advanced Engineering Mathematics I Assignments - Spring 2014 Dr. E. Jacobs The main texts for this course are Calculus by James Stewart and Fundamentals of Differential Equations by Nagle, Saff
More informationKinematics (2) - Motion in Three Dimensions
Kinematics (2) - Motion in Three Dimensions 1. Introduction Kinematics is a branch of mechanics which describes the motion of objects without consideration of the circumstances leading to the motion. 2.
More informationSolutions to Practice Exam 2
Solutions to Practice Eam Problem : For each of the following, set up (but do not evaluate) iterated integrals or quotients of iterated integral to give the indicated quantities: Problem a: The average
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationMath 221 Examination 2 Several Variable Calculus
Math Examination Spring Instructions These problems should be viewed as essa questions. Before making a calculation, ou should explain in words what our strateg is. Please write our solutions on our own
More informationSolutions for the Practice Final - Math 23B, 2016
olutions for the Practice Final - Math B, 6 a. True. The area of a surface is given by the expression d, and since we have a parametrization φ x, y x, y, f x, y with φ, this expands as d T x T y da xy
More informationPoincaré Models of Hyperbolic Geometry
Chapter 9 Poincaré Models of Hyperbolic Geometry 9.1 The Poincaré Upper Half Plane Model The next model of the hyperbolic plane that we will consider is also due to Henri Poincaré. We will be using the
More informationNumerical solution of an Eikonal equation on a Graph
Numerical solution of an Eikonal equation on a Graph Adriano Festa Sapienza" Univ. of Roma Math. Dep. oint work with Fabio Camilli ( Sapienza University, Rome) London. September 2011 SADCO Workshop Outline
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 informationComplex Analysis MATH 6300 Fall 2013 Homework 4
Complex Analysis MATH 6300 Fall 2013 Homework 4 Due Wednesday, December 11 at 5 PM Note that to get full credit on any problem in this class, you must solve the problems in an efficient and elegant manner,
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 informationConstrained Optimization in Two Variables
in Two Variables James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 17, 216 Outline 1 2 What Does the Lagrange Multiplier Mean? Let
More informationOne side of each sheet is blank and may be used as scratch paper.
Math 244 Spring 2017 (Practice) Final 5/11/2017 Time Limit: 2 hours Name: No calculators or notes are allowed. One side of each sheet is blank and may be used as scratch paper. heck your answers whenever
More informationMaximum-norm a posteriori estimates for discontinuous Galerkin methods
Maximum-norm a posteriori estimates for discontinuous Galerkin methods Emmanuil Georgoulis Department of Mathematics, University of Leicester, UK Based on joint work with Alan Demlow (Kentucky, USA) DG
More informationDeposited on: 19 July 2012
Scardia, L. 6 he nonlinear bending-torsion theory for curved rods as Gamma-limit of three-dimensional elasticity. Asymptotic Analysis, 47 3-4. pp. 37-343. ISSN 9-734 http://eprints.gla.ac.uk/6763 eposited
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationLecture 3: Asymptotic Methods for the Reduced Wave Equation
Lecture 3: Asymptotic Methods for the Reduced Wave Equation Joseph B. Keller 1 The Reduced Wave Equation Let v (t,) satisfy the wave equation v 1 2 v c 2 = 0, (1) (X) t2 where c(x) is the propagation speed
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationContents. Prologue Introduction. Classical Approximation... 19
Contents Prologue........................................................................ 15 1 Introduction. Classical Approximation.................................. 19 1.1 Introduction................................................................
More information3. The linear 3-D elasticity mathematical model
3. The linear 3-D elasticity mathematical model In Chapter we examined some fundamental conditions that should be satisfied in the modeling of all deformable solids and structures. The study of truss structures
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
More informationx+ y = 50 Dividing both sides by 2 : ( ) dx By (7.2), x = 25m gives maximum area. Substituting this value into (*):
Solutions 7(b 1 Complete solutions to Exercise 7(b 1. Since the perimeter 100 we have x+ y 100 [ ] ( x+ y 50 ividing both sides by : y 50 x * The area A xy, substituting y 50 x gives: A x( 50 x A 50x x
More informationME 509, Spring 2016, Final Exam, Solutions
ME 509, Spring 2016, Final Exam, Solutions 05/03/2016 DON T BEGIN UNTIL YOU RE TOLD TO! Instructions: This exam is to be done independently in 120 minutes. You may use 2 pieces of letter-sized (8.5 11
More informationMath 185 Fall 2015, Sample Final Exam Solutions
Math 185 Fall 2015, Sample Final Exam Solutions Nikhil Srivastava December 12, 2015 1. True or false: (a) If f is analytic in the annulus A = {z : 1 < z < 2} then there exist functions g and h such that
More informationSolutions to Practice Test 3
Solutions to Practice Test 3. (a) Find the equation for the plane containing the points (,, ), (, 2, ), and (,, 3). (b) Find the area of the triangle with vertices (,, ), (, 2, ), and (,, 3). Answer: (a)
More informationOn the Hardy constant of non-convex planar domains: the case of the quadrilateral
On the Hardy constant of non-convex planar domains: the case of the quadrilateral Gerassimos Barbatis University of Athens joint work with Achilles Tertikas University of Crete Heraklion, Hardy constant
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 informationWithout fully opening the exam, check that you have pages 1 through 12.
Name: Section: Recitation Instructor: INSTRUCTIONS Fill in your name, etc. on this first page. Without fully opening the exam, check that you have pages 1 through 12. Show all your work on the standard
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 informationQuiz 4 (Discussion Session) Phys 1302W.400 Spring 2018
Quiz 4 (Discussion ession) Phys 1302W.400 pring 2018 This group quiz consists of one problem that, together with the individual problems on Friday, will determine your grade for quiz 4. For the group problem,
More informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationWilson loops at strong coupling for curved contours with cusps. Contour with two cusps, formed by segments of two circles
1 Wilson loops at strong coupling for curved contours with cusps based on H.D. 1509.00222 Motivation and introduction Contour with two cusps, formed by segments of two circles Cusp anomalous dimension
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationCANONICAL EQUATIONS. Application to the study of the equilibrium of flexible filaments and brachistochrone curves. By A.
Équations canoniques. Application a la recherche de l équilibre des fils flexibles et des courbes brachystochrones, Mem. Acad. Sci de Toulouse (8) 7 (885), 545-570. CANONICAL EQUATIONS Application to the
More informationSTRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD
Shape optimization 1 G. Allaire STRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD G. ALLAIRE CMAP, Ecole Polytechnique with F. JOUVE and A.-M. TOADER 1. Introduction 2. Setting of the problem 3. Shape differentiation
More informationMAE 101A. Homework 7 - Solutions 3/12/2018
MAE 101A Homework 7 - Solutions 3/12/2018 Munson 6.31: The stream function for a two-dimensional, nonviscous, incompressible flow field is given by the expression ψ = 2(x y) where the stream function has
More information1 Integration in many variables.
MA2 athaye Notes on Integration. Integration in many variables.. Basic efinition. The integration in one variable was developed along these lines:. I f(x) dx, where I is any interval on the real line was
More informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More information