Accurate Gradient Computation for Interface Problems
|
|
- Jennifer Charles
- 5 years ago
- Views:
Transcription
1 Accurate Gradient Computation for Interface Problems Nnenna Anoruo, Kevin Luna, William Reese, Justin Rivera Faculty Advsior: Zhilin Li North Carolina State University February 11, 016
2 Outline Problem : Solve ODE/PDE Interface Problems and Derivatives Piecewise nonconstant coecient problems Polar Coordinates (D) Spherical Coordinates (3D) Achieve second order accurate solutions and derivatives using FD method Approach Worked on nonconstant coecent interface problems and their derivatives Solved using Finite Dierence Method/IIM Worked on axis symmetric spherical and polar coordinate Work in progress Used the Immersed Interface Method 1D Stefan Problem 1D Moving Boundary Problem
3 Acheivements Derivatives on boundaries for regular problems Piecewise constant coecient 1D general problem Polar coordinates (D) Spherical coordinates (3D)
4 An Example of an Interface Problem Rubber band example u = Cδ(x α), 0 < x < 1, 0 < α < 1 u(0) = u(1) = 0 Consider an elastic material with two xed end points. The solution is continuous, but the derivative is not
5 Interfaces What is an interface? An interface is essentially a point (1D), a curve (D), or a surface (3D), that connects dierent states of a system Some examples include : ice and water, a string with dierent materials on it, the heating of an embedded metal block in another material with a dierent heat conductivity, etc. Simple model for Interface Problem (β(x)u x ) x = f(x) + vδ(x α), 0 < x < 1, 0 < α < 1 Where f(x) C (0, 1), and β(x) C 1 (0, α) C 1 (α, 1) β(x) = { β(x) 0 < x < α β(x) + α < x < 1 (1) Discontinous β(x) β(x) x
6 Mathematical Formulation of Rubber Band Problem in general Two equivalent ways to formulate the problem: On a single domain: (βu ) = Cδ(x α) 0 < x < 1 u(0) = u(1) = 0 Equivalently, splitting the domain so that the point source is avoided (βu x ) x = 0 x (0, α); (βu x ) x = 0 x (α, 1) [u] = 0, [βu x ] = C Note here [u] notation actually implies the following: [u] x=α = lim u(x) lim u(x) x α + x α
7 Immersed Interface Method (βu x ) x σu = f(x) + Cδ (x α), 0 < x < 1, 0 < α < 1 Peskin's method approximates the delta function approximation is 1st order accurate Immersed Interface Method reformulates the problem USING the jump condition (βu x ) x = f(x), x (0, α) (α, 1)
8 Finite dierence scheme at an irregular grid point We begin with the general FD stencil, γ 1 U j 1 + γ U j + γ 3 U j+1 = f j To nd the scheme at the irregular points, we want to minimize the local truncation error at α T j = γ 1 U j 1 + γ U j + γ 3 U j+1 f(x j ) C j We use method of undetermined coecients to solve for γ 1, γ, γ 3,and C j such that the local truncation error is minimized in magnitude by using the interface relations u + = w + u, u + x = β β + u x + v β +, u+ xx = β β + u xx
9 Finite dierence scheme at an irregular grid point T j = γ 1 u j 1 + γ u j + γ 3 u j+1 f(x j ) C j u(x j 1 ) = u + u x (x j 1 α) + (x j 1 α) u xx + O(h 3 ) u(x j ) = u + u x (x j α) + (x j α) u xx + O(h 3 ) u(x j+1 ) = u + + u + x (x j+1 α) + (x j+1 α) u + xx + O(h 3 ) = u + w + β β + u x (x j+1 α) + β (x j+1 α) β + u xx + O(h 3 )
10 FD at the irregular grid point, lll Taylor Expansion + interface relations + undetermined coecient method T j = γ j 1 u(x j 1 ) + γ j u(x j ) + γ j+1 u(x j+1 ) f j C j = γ j 1 (u + u x (x j 1 α) + (x j 1 α) u xx) + O(h 3 )) +γ j (u + u x (x j α) + (x j α) u β xx) + γ j+1 β + C(x j+1 α) +γ j+1 (u + w + β β + u x (x j+1 α) + β (x j 1 α) β + u xx) + O(h 3 )) = u (γ j 1 + γ j + γ j+1 ) + u x (γ j 1 α) + γ j (x j α) β +γ j+1 β + (x j+1 α)) + u (x j 1 α) xx(γ j 1 β (x j+1 α) γ j+1 β + ) f j C j x j α) + γ j +
11 Set-up equation The linear system for the coecients: γ j 1 + γ j + γ j+1 = 0 γ j 1 (x j 1 α) + γ j (x j α) + γ j+1 β γ j 1 (x j 1 α) + γ j (x j α) β + (x j+1 α) = 0 β (x j+1 α) + γ j+1 β + = β The correction term is ( ) C j = γ j+1 w + C β β + (x j+1 α)
12 Derivatives at Interface (Two Sided) To approximate u (x) we use the IIM on the approximated solutions of u(x) which were given by the IIM We collect the terms by u and its derivatives and set u equal to 1 and the other linear equations equal to 0 as follows: γ 1 + γ + γ 3 = 0 γ 1 (x j 1 α) + γ (x j α) + γ 3 β γ 1 (x j 1 α) and the correction term + γ (x j α) C j = γ 1 (w + β + (x j+1 α) = 1 β (x j+1 α) + γ 3 β + = 0 c β + (x j+1 α)). When this system is solved for γ 1, γ, and γ 3 we are left with the sided stencil for u (x) at an interface.
13 1-D Interface w/ Nonconstant Coecients In general 1D Interface problems where x (0, α) (α, 1) are expressed as: (β x u) x σu = f For our purposes we explored functions that have nite jumps at the interface α In all our examples we dened β to be: { 1 + cos(x) β(x) = if x α 1 + x if x > α and σ to be: σ(x) = { (1 + x) if x α ln( + x) if x > α Note: The appropriate Dirichlet boundary condition was chosen in each example according to the actual solution
14 1-D Interface w/ Non-constant Coecients Consider the following example: (10) ( ( cos(x) + 1)) sin(10x) 10 sin(x) cos(10x) (1 + x) cos(10x) if x < α f(x) = (30)( 30(x + 1) cos(30x) x sin(30x)) ln( + x) cos(30x) if x < α. The exact solution to the example can be shown as follows: { sin(10x) if x α u(x) = cos(30x) if x > α ()
15 Solution and Error Plot Still considering the same example: This gure shows the computed solution with n = 180 grid divisions. The actual solution has discontinuity that is connected by a plotter. This gure displays the error dierence of the estimated solution and actual solution at each grid point Error (difference) Plot
16 Error Plots on Log-log Scale Grid renement analysis of approximated solution with 15 renements starting at n = 0. Note the slope of the line is.1 implying that it is second order accurate. Grid renement analysis of the right derivative at the interface using one and two sided IIM methods. The slope of the lines are 1.08 and respectively. 1 These new results can be seen on Dr.Li's website: zhilin/
17 Polar coordinates We considered axis symmetric cases of interface problems in polar coordinates 1 r (rβu r) r = f(r) We generate a uniform grid r i = R 1 + i r, i = 0, 1,..., m, r = R R 1 m Standard dierence scheme for polar ( r i 1 + r i+ 1 1 r i 1 U i 1 r i ) U i + r i+ 1 U i+1 ( r) = f(r i )
18 Polar coordinates When the origin is in the boundary there is a singularity We use a staggered grid to address this issue r i = (i 1 ) r, r = R m 1, i = 1,...m. When i = 1 we have to use U 1 + U ( r) + 1 r 1 U r = f(r 1) For all other points on the grid we use the standard dierence scheme We can apply IIM at an interface if needed
19 Polar Coordinates with Axis Symmetry, Ex1 1 r (rβu r) r = f(r), 0 < r < 1, α = 3 { β = 1 r α β = β + = 100 r > α { f (sin(r)+r cos(r)) = f(r) = r r α f + = 9 cos(9r) r 81 sin(9r) r > α
20 Polar Coordinates with Axis Symmetry, Ex1 ctd. Exact Solution u(r) = { u = cos(r) u + = sin(9r) r α r > α Surface of Approximation
21 Solution and Error Plot Still considering the same example: This gure shows the computed solution with n = 160 grid divisions. The actual solution has discontinuity that is connected by a plotter. This gure displays the error dierence of the estimated solution and actual solution at each grid point
22 Error Plots on Log-log Scale Grid renement analysis of approximated solution with 8 renements starting at n = 0. Note the slope of the line is.1 implying that it is second order accurate. Grid renement analysis of the right derivative at the interface using one and two sided IIM methods. The slopes of the lines are 0.96 and 1.94 respectively.
23 Polar Coordinates with Axis Symmetry, Ex 1 r (rβu r) r = f(r), 0 < r < 1, α = 3 { β = 1 r α β = β + = 100 r > α f(r) = { f = sin(r)+r cos(r) r f + = r ( sin(r))+sin(r) r cos(r) r 3 r α r > α
24 Polar Coordinates with Axis Symmetry, Ex ctd. Exact Solution u(r) = Surface of Approximation { u = cos(r) u + = sin(r) r r α r > α
25 Solution and Error Plot Still considering the same example: This gure shows the computed solution with n = 160 grid divisions. The actual solution has discontinuity that is connected by a plotter Solution and Approximation plot for β (r u r ) r = rf(r) Approximation Exact Solution This gure displays the error dierence of the estimated solution and actual solution at each grid point Error (difference) Plot u(r) 0.6 U-u(x) r x
26 Error Plots on Log-log Scale Grid renement analysis of approximated solution with 8 renements starting at n = 0. Note the slope of the line is 1.99 implying that it is second order accurate Solution Error Plot on a log-log Scale for the IIM Immersed Interface Method Grid renement analysis of the right derivative at the interface using one and two sided IIM methods. The slopes of the lines are 0.9 and.06 respectively Two-Sided Method(IIM) One-sided method loglog Plot of U x + (α) Error log 10 ( U-u ) log 10 ( U x - (α)-ux - (α) ) log 10 (h) log 10 (h)
27 3D Spherical coordinates with Radial Symmetry We considered axis symmetric cases of interface problem in spherical coordinates 1 r (r βu r ) r = f(r) We generate the staggered grid r i = (i 1 ) r, r = R m 1, i = 1,...m. The standard dierence scheme for spherical coordinates (r i 1 ) U i 1 ( (r i 1 ) + (r i+ 1 ) ) U i + (r i+ 1 ) U i+1 ( r) = r f(r i )
28 Spherical coordinates For r i=1 the dierence scheme is ( ) r 1 U r + U1 + U r 1 ( r) U 0 4 = f(r 1) For all other points on the grid we use the standard dierence scheme We can apply IIM at an interface if needed
29 3D Spherical coordinates with Radial Symmetry, Ex1 1 r (r βu r ) r = f(r), 0 < r < 1, α = 3 { β = 1 r α β = β + = 100 r > α { f = 3r(4 sin(r 3 ) + 3r 3 cos(r 3 )) r α f(r) = f + = 6 cos(r ) 4r sin(r ) r > α
30 3D Spherical coordinates with Radial Symmetry, Ex1 ctd. Exact Solution u(r) = { u = cos(r 3 ) u + = sin(r ) r α r > α Surface of Approximation
31 Solution and Error Plot Still considering the same example: This gure shows the computed solution with n = 160 grid divisions. The actual solution has discontinuity that is connected by a plotter. This gure displays the error dierence of the estimated solution and actual solution at each grid point
32 Error Plots on Log-log Scale Grid renement analysis of approximated solution with 8 renements starting at n = 0. Note the slope of the line is 1.9 implying that it is second order accurate. Grid renement analysis of the right derivative at the interface using one and two sided IIM methods. The slopes of the lines are 0.99 and 1.99 respectively.
33 Work in Progress Stefan problem ( xed and moving boundaries) Parabolic interface problems Moving interface problems
34 Acknowledgements We would like to thank NCSU Department of Mathematics for hosting the REU program Faculty Advisor Dr. Li The NSF for grant DMS The NSA for grant H
ABSTRACT. coefficients. The main motivation is to get not only a second order accurate solution but
ABSTRACT CHEN, GUANYU. Accurate Gradient Computation for Elliptic Interface Problems with Discontinuous and Variable Coefficients. (Under the direction of Dr. Zhilin Li.) A new numerical method is proposed
More informationManoj Kumar 1, Pratibha Joshi 2. (Received 31 August 2011, accepted 26 December 2011)
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.14(2012) No.1,pp.11-22 A Mathematical Model and Numerical Solution of a One Dimensional Steady State Heat Conduction
More informationAN OVERVIEW OF THE IMMERSED INTERFACE METHOD AND ITS APPLICATIONS
. AN OVERVIEW OF THE IMMERSED INTERFACE METHOD AND ITS APPLICATIONS ZHILIN LI Abstract. Interface problems have many applications. Mathematically, interface problems usually lead to differential equations
More informationTaylor Series. Math114. March 1, Department of Mathematics, University of Kentucky. Math114 Lecture 18 1/ 13
Taylor Series Math114 Department of Mathematics, University of Kentucky March 1, 2017 Math114 Lecture 18 1/ 13 Given a function, can we find a power series representation? Math114 Lecture 18 2/ 13 Given
More informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for this
More informationUpscaling Wave Computations
Upscaling Wave Computations Xin Wang 2010 TRIP Annual Meeting 2 Outline 1 Motivation of Numerical Upscaling 2 Overview of Upscaling Methods 3 Future Plan 3 Wave Equations scalar variable density acoustic
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationMA1021 Calculus I B Term, Sign:
MA1021 Calculus I B Term, 2014 Final Exam Print Name: Sign: Write up your solutions neatly and show all your work. 1. (28 pts) Compute each of the following derivatives: You do not have to simplify your
More informationSection 10.7 Taylor series
Section 10.7 Taylor series 1. Common Maclaurin series 2. s and approximations with Taylor polynomials 3. Multiplication and division of power series Math 126 Enhanced 10.7 Taylor Series The University
More informationA cartesian grid finite volume method for the solution of the Poisson equation with variable coefficients and embedded interfaces
A cartesian grid finite volume method for the solution of the Poisson equation with variable coefficients and embedded interfaces M. Oevermann a, and R. Klein b, a Technische Universität Berlin, Institut
More informationANALYSIS AND NUMERICAL METHODS FOR SOME CRACK PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 2, Number 2-3, Pages 155 166 c 2011 Institute for Scientific Computing and Information ANALYSIS AND NUMERICAL METHODS FOR SOME
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More information1e N
Spectral schemes on triangular elements by Wilhelm Heinrichs and Birgit I. Loch Abstract The Poisson problem with homogeneous Dirichlet boundary conditions is considered on a triangle. The mapping between
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationOn the Solution of the Elliptic Interface Problems by Difference Potentials Method
On the Solution of the Elliptic Interface Problems by Difference Potentials Method Yekaterina Epshteyn and Michael Medvinsky Abstract Designing numerical methods with high-order accuracy for problems in
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More informationf (x) = k=0 f (0) = k=0 k=0 a k k(0) k 1 = a 1 a 1 = f (0). a k k(k 1)x k 2, k=2 a k k(k 1)(0) k 2 = 2a 2 a 2 = f (0) 2 a k k(k 1)(k 2)x k 3, k=3
1 M 13-Lecture Contents: 1) Taylor Polynomials 2) Taylor Series Centered at x a 3) Applications of Taylor Polynomials Taylor Series The previous section served as motivation and gave some useful expansion.
More informationA cartesian grid finite volume method for the solution of the Poisson equation with variable coefficients and embedded interfaces
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-495 Berlin-Dahlem Germany MICHAEL OEVERMANN, RUPERT KLEIN A cartesian grid finite volume method for the solution of the Poisson equation
More informationTheta Functions MAΘ National Convention 2018
Theta Functions MAΘ National Convention 018 Note that choice stands for None of the above answers is correct 1. How many times do y = x and y = x 3 intersect? A) 0 B) 1 C) D) 3. Let f(x) = x + 19 and g(x)
More informationPROPOSITION of PROJECTS necessary to get credit in the 'Introduction to Parallel Programing in 2014'
PROPOSITION of PROJECTS necessary to get credit in the 'Introduction to Parallel Programing in 2014' GENERALITIES The text may be written in Polish, or English. The text has to contain two parts: 1. Theoretical
More informationExamples of the Fourier Theorem (Sect. 10.3). The Fourier Theorem: Continuous case.
s of the Fourier Theorem (Sect. 1.3. The Fourier Theorem: Continuous case. : Using the Fourier Theorem. The Fourier Theorem: Piecewise continuous case. : Using the Fourier Theorem. The Fourier Theorem:
More informationMath 217 Practice Exam 1. Page Which of the following differential equations is exact?
Page 1 1. Which of the following differential equations is exact? (a) (3x x 3 )dx + (3x 3x 2 ) d = 0 (b) sin(x) dx + cos(x) d = 0 (c) x 2 x 2 = 0 (d) (1 + e x ) + (2 + xe x ) d = 0 CORRECT dx (e) e x dx
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationMinimum and maximum values *
OpenStax-CNX module: m17417 1 Minimum and maximum values * Sunil Kumar Singh This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 2.0 In general context, a
More informationReview: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function.
Taylor Series (Sect. 10.8) Review: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function. Review: Power series define functions Remarks:
More informationNotation Nodes are data points at which functional values are available or at which you wish to compute functional values At the nodes fx i
LECTURE 6 NUMERICAL DIFFERENTIATION To find discrete approximations to differentiation (since computers can only deal with functional values at discrete points) Uses of numerical differentiation To represent
More informationMATH 31B: MIDTERM 2 REVIEW. sin 2 x = 1 cos(2x) dx = x 2 sin(2x) 4. + C = x 2. dx = x sin(2x) + C = x sin x cos x
MATH 3B: MIDTERM REVIEW JOE HUGHES. Evaluate sin x and cos x. Solution: Recall the identities cos x = + cos(x) Using these formulas gives cos(x) sin x =. Trigonometric Integrals = x sin(x) sin x = cos(x)
More informationTaylor series - Solutions
Taylor series - Solutions. f(x) sin(x) sin(0) + x cos(0) + x x ( sin(0)) +!! ( cos(0)) + + x4 x5 (sin(0)) + 4! 5! 0 + x + 0 x x! + x5 5! x! + 0 + x5 (cos(0)) + x6 6! ( sin(0)) + x 7 7! + x9 9! 5! + 0 +
More informationExplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code
Eplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code V. Rutka A. Wiegmann November 25, 2005 Abstract The Eplicit Jump Immersed Interface method is a powerful tool to solve elliptic
More informationLecture 38 Insulated Boundary Conditions
Lecture 38 Insulated Boundary Conditions Insulation In many of the previous sections we have considered fixed boundary conditions, i.e. u(0) = a, u(l) = b. We implemented these simply by assigning u j
More informationMath 226 Calculus Spring 2016 Practice Exam 1. (1) (10 Points) Let the differentiable function y = f(x) have inverse function x = f 1 (y).
Math 6 Calculus Spring 016 Practice Exam 1 1) 10 Points) Let the differentiable function y = fx) have inverse function x = f 1 y). a) Write down the formula relating the derivatives f x) and f 1 ) y).
More informationDenition and some Properties of Generalized Elementary Functions of a Real Variable
Denition and some Properties of Generalized Elementary Functions of a Real Variable I. Introduction The term elementary function is very often mentioned in many math classes and in books, e.g. Calculus
More informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationby Ron Buckmire Mathematics Department Occidental College Los Angeles, CA , U.S.A.
A new nite-dierence scheme for singular dierential equations in cylinical or spherical coordinates by Ron Buckmire Mathematics Department Occidental College Los Angeles, CA 90041-3338, U.S.A. December
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
More informationMat1062: Introductory Numerical Methods for PDE
Mat1062: Introductory Numerical Methods for PDE Mary Pugh January 13, 2009 1 Ownership These notes are the joint property of Rob Almgren and Mary Pugh 2 Boundary Conditions We now want to discuss in detail
More informationContinuity, Intermediate Value Theorem (2.4)
Continuity, Intermediate Value Theorem (2.4) Xiannan Li Kansas State University January 29th, 2017 Intuitive Definition: A function f(x) is continuous at a if you can draw the graph of y = f(x) without
More informationProblem 1 Kaplan, p. 436: 2c,i
Mathematical Methods I - Fall 03 Homework Solutions Page Problem Kaplan, p 36: c,i Find the first three nonzero terms of the following Taylor series: c) ln + x) about x = 0 i) arctanh x about x = 0 The
More informationA Fourth Order Accurate Discretization for the Laplace and Heat Equations on Arbitrary Domains, with Applications to the Stefan Problem
A Fourth Order Accurate Discretization for the Laplace and Heat Equations on Arbitrary Domains, with Applications to the Stefan Problem Frédéric Gibou Ronald Fedkiw April 27, 2004 Abstract In this paper,
More informationComputer Problems for Fourier Series and Transforms
Computer Problems for Fourier Series and Transforms 1. Square waves are frequently used in electronics and signal processing. An example is shown below. 1 π < x < 0 1 0 < x < π y(x) = 1 π < x < 2π... and
More informationAn Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains
An Immersed Interace Method or the Incompressible Navier-Stokes Equations in Irregular Domains Duc-Vinh Le, Boo Cheong Khoo, Jaime Peraire Singapore-MIT Alliance Department o Mechanical Engineering, National
More informationHigh-order di erence potentials methods for 1D elliptic type models. Yekaterina Epshteyn and Spencer Phippen
High-order di erence potentials methods for 1D elliptic type models Yekaterina Epshteyn and Spencer Phippen Department of Mathematics, The University of Utah, Salt Lake City, UT, 84112 Abstract Numerical
More informationMath 250 Skills Assessment Test
Math 5 Skills Assessment Test Page Math 5 Skills Assessment Test The purpose of this test is purely diagnostic (before beginning your review, it will be helpful to assess both strengths and weaknesses).
More informationMTAEA Differentiation
School of Economics, Australian National University February 5, 2010 Basic Properties of the Derivative. Secant Tangent Applet l 3 l 2 l 1 a a 3 a 2 a 1 Figure: The derivative of f at a is the limiting
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationThe First Derivative and Second Derivative Test
The First Derivative and Second Derivative Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University April 9, 2018 Outline 1 Extremal Values 2
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationTaylor and Maclaurin Series. Approximating functions using Polynomials.
Taylor and Maclaurin Series Approximating functions using Polynomials. Approximating f x = e x near x = 0 In order to approximate the function f x = e x near x = 0, we can use the tangent line (The Linear
More informationProperties of Arithmetical Functions
Properties of Arithmetical Functions Zack Clark Math 336, Spring 206 Introduction Arithmetical functions, as dened by Delany [2], are the functions f(n) that take positive integers n to complex numbers.
More informationIncreasing and decreasing intervals *
OpenStax-CNX module: m15474 1 Increasing and decreasing intervals * Sunil Kumar Singh This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 2.0 A function is
More informationA NUMERICAL METHOD FOR SOLVING THE VARIABLE COEFFICIENT WAVE EQUATION WITH INTERFACE JUMP CONDITIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 6, Number, Pages 7 c 9 Institute for Scientific Computing and Information A NUMERICAL METHOD FOR SOLVING THE VARIABLE COEFFICIENT WAVE EQUATION
More informationDUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs
DUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs MATH 4354 Fall 2005 December 5, 2005 1 Duhamel s Principle
More informationEffective matrix-free preconditioning for the augmented immersed interface method
Effective matrix-free preconditioning for the augmented immersed interface method Jianlin Xia a, Zhilin Li b, Xin Ye a a Department of Mathematics, Purdue University, West Lafayette, IN 47907, USA. E-mail:
More informationJ.TAUSCH AND J.WHITE 1. Capacitance Extraction of 3-D Conductor Systems in. Dielectric Media with high Permittivity Ratios
JTAUSCH AND JWHITE Capacitance Extraction of -D Conductor Systems in Dielectric Media with high Permittivity Ratios Johannes Tausch and Jacob White Abstract The recent development of fast algorithms for
More informationStudy # 1 11, 15, 19
Goals: 1. Recognize Taylor Series. 2. Recognize the Maclaurin Series. 3. Derive Taylor series and Maclaurin series representations for known functions. Study 11.10 # 1 11, 15, 19 f (n) (c)(x c) n f(c)+
More informationMath 113 (Calculus 2) Exam 4
Math 3 (Calculus ) Exam 4 November 0 November, 009 Sections 0, 3 7 Name Student ID Section Instructor In some cases a series may be seen to converge or diverge for more than one reason. For such problems
More informationThe First Derivative and Second Derivative Test
The First Derivative and Second Derivative Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 2017 Outline Extremal Values The
More informationApplied Numerical Mathematics. High-order numerical schemes based on difference potentials for 2D elliptic problems with material interfaces
Applied Numerical Mathematics 111 (2017) 64 91 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum High-order numerical schemes based on difference potentials
More informationTaylor and Maclaurin Series. Approximating functions using Polynomials.
Taylor and Maclaurin Series Approximating functions using Polynomials. Approximating f x = e x near x = 0 In order to approximate the function f x = e x near x = 0, we can use the tangent line (The Linear
More informationLecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018
Lecture 10 Partial derivatives Dan Nichols nichols@math.umass.edu MATH 233, Spring 2018 University of Massachusetts February 27, 2018 Last time: functions of two variables f(x, y) x and y are the independent
More informationAdvanced Calculus Math 127B, Winter 2005 Solutions: Final. nx2 1 + n 2 x, g n(x) = n2 x
. Define f n, g n : [, ] R by f n (x) = Advanced Calculus Math 27B, Winter 25 Solutions: Final nx2 + n 2 x, g n(x) = n2 x 2 + n 2 x. 2 Show that the sequences (f n ), (g n ) converge pointwise on [, ],
More informationInstitute for Advanced Computer Studies. Department of Computer Science. On the Convergence of. Multipoint Iterations. G. W. Stewart y.
University of Maryland Institute for Advanced Computer Studies Department of Computer Science College Park TR{93{10 TR{3030 On the Convergence of Multipoint Iterations G. W. Stewart y February, 1993 Reviseed,
More information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationCalculus and linear algebra for biomedical engineering Week 11: The Riemann integral and its properties
Clculus nd liner lgebr for biomedicl engineering Week 11: The Riemnn integrl nd its properties Hrtmut Führ fuehr@mth.rwth-chen.de Lehrstuhl A für Mthemtik, RWTH Achen Jnury 9, 2009 Overview 1 Motivtion:
More informationNonconstant Coefficients
Chapter 7 Nonconstant Coefficients We return to second-order linear ODEs, but with nonconstant coefficients. That is, we consider (7.1) y + p(t)y + q(t)y = 0, with not both p(t) and q(t) constant. The
More informationNumerical solution of the 2-D Poisson equation on an irregular domain with Robin boundary conditions
Numerical solution of the 2-D Poisson equation on an irregular domain with Robin boundary conditions Z. Jomaa C. Macaskill August 8, 28 Abstract We describe a 2-D finite difference algorithm for inverting
More informationb n x n + b n 1 x n b 1 x + b 0
Math Partial Fractions Stewart 7.4 Integrating basic rational functions. For a function f(x), we have examined several algebraic methods for finding its indefinite integral (antiderivative) F (x) = f(x)
More information2 In this paper, we give a brief summary of the front tracking algorithm for axisymmetric ow. We report simulations of spherical shock refraction by a
A Fast Algorithm for Moving Interface Problems S. Dutta 1,J.Glimm 1 2,J.W.Grove 3,D.H.Sharp 4, and Y. Zhang 1 1 Department of Applied Mathematics and Statistics, University at Stony Brook, Stony Brook,
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More information8.7 Taylor s Inequality Math 2300 Section 005 Calculus II. f(x) = ln(1 + x) f(0) = 0
8.7 Taylor s Inequality Math 00 Section 005 Calculus II Name: ANSWER KEY Taylor s Inequality: If f (n+) is continuous and f (n+) < M between the center a and some point x, then f(x) T n (x) M x a n+ (n
More informationMATH 1231 MATHEMATICS 1B CALCULUS. Section 5: - Power Series and Taylor Series.
MATH 1231 MATHEMATICS 1B CALCULUS. Section 5: - Power Series and Taylor Series. The objective of this section is to become familiar with the theory and application of power series and Taylor series. By
More information2 Multidimensional Hyperbolic Problemss where A(u) =f u (u) B(u) =g u (u): (7.1.1c) Multidimensional nite dierence schemes can be simple extensions of
Chapter 7 Multidimensional Hyperbolic Problems 7.1 Split and Unsplit Dierence Methods Our study of multidimensional parabolic problems in Chapter 5 has laid most of the groundwork for our present task
More informationand lim lim 6. The Squeeze Theorem
Limits (day 3) Things we ll go over today 1. Limits of the form 0 0 (continued) 2. Limits of piecewise functions 3. Limits involving absolute values 4. Limits of compositions of functions 5. Limits similar
More informationReview of Section 1.1. Mathematical Models. Review of Section 1.1. Review of Section 1.1. Functions. Domain and range. Piecewise functions
Review of Section 1.1 Functions Mathematical Models Domain and range Piecewise functions January 19, 2017 Even and odd functions Increasing and decreasing functions Mathematical Models January 19, 2017
More informationMa 530 Power Series II
Ma 530 Power Series II Please note that there is material on power series at Visual Calculus. Some of this material was used as part of the presentation of the topics that follow. Operations on Power Series
More informationBi. lkent Calculus II Exams
Bi. lkent Calculus II Exams 988-208 Spring 208 Midterm I........ Spring 208 Midterm II....... 2 Spring 207 Midterm I........ 4 Spring 207 Midterm II....... 5 Spring 207 Final........... 7 Spring 206 Midterm
More informationNumerical Schemes from the Perspective of Consensus
Numerical Schemes from the Perspective of Consensus Exploring Connections between Agreement Problems and PDEs Department of Electrical Engineering and Computer Sciences University of California, Berkeley
More informationVoronoi cell nite dierence method for the diusion operator on arbitrary unstructured grids
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 003; 57:1 34 (DOI: 10.100/nme.664) Voronoi cell nite dierence method for the diusion operator on arbitrary unstructured
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
More informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More informationMath 106 Fall 2014 Exam 2.1 October 31, ln(x) x 3 dx = 1. 2 x 2 ln(x) + = 1 2 x 2 ln(x) + 1. = 1 2 x 2 ln(x) 1 4 x 2 + C
Math 6 Fall 4 Exam. October 3, 4. The following questions have to do with the integral (a) Evaluate dx. Use integration by parts (x 3 dx = ) ( dx = ) x3 x dx = x x () dx = x + x x dx = x + x 3 dx dx =
More informationProblem Set 1 October 30, 2017
1. e π can be calculated from e x = x n. But that s not a good approximation method. The n! reason is that π is not small compared to 1. If you want O(0.1) accuracy, then the last term you need to include
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationMATH 118, LECTURES 27 & 28: TAYLOR SERIES
MATH 8, LECTURES 7 & 8: TAYLOR SERIES Taylor Series Suppose we know that the power series a n (x c) n converges on some interval c R < x < c + R to the function f(x). That is to say, we have f(x) = a 0
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )
More information10. e tan 1 (y) 11. sin 3 x
MATH B FINAL REVIEW DISCLAIMER: WHAT FOLLOWS IS A LIST OF PROBLEMS, CONCEPTUAL QUESTIONS, TOPICS, AND SAMPLE PROBLEMS FROM THE TEXTBOOK WHICH COMPRISE A HEFTY BUT BY NO MEANS EXHAUSTIVE LIST OF MATERIAL
More informationAn Attempt To Understand Tilling Approach In proving The Littlewood Conjecture
An Attempt To Understand Tilling Approach In proving The Littlewood Conjecture 1 Final Report For Math 899- Dr. Cheung Amira Alkeswani Dr. Cheung I really appreciate your acceptance in adding me to your
More informationReview Smoothing Spatial Filters Sharpening Spatial Filters. Spatial Filtering. Dr. Praveen Sankaran. Department of ECE NIT Calicut.
Spatial Filtering Dr. Praveen Sankaran Department of ECE NIT Calicut January 7, 203 Outline 2 Linear Nonlinear 3 Spatial Domain Refers to the image plane itself. Direct manipulation of image pixels. Figure:
More informationMath Numerical Analysis
Math 541 - Numerical Analysis Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center San Diego State University
More informationCompact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems
Compact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems Daniel Ritter 1, Ulrich Rüde 1, Björn Gmeiner 1, Rochus Schmid 2 Copper Mountain, March 18th, 2013
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
More informationDRAFT - Math 101 Lecture Note - Dr. Said Algarni
2 Limits 2.1 The Tangent Problems The word tangent is derived from the Latin word tangens, which means touching. A tangent line to a curve is a line that touches the curve and a secant line is a line that
More information3 Green s functions in 2 and 3D
William J. Parnell: MT34032. Section 3: Green s functions in 2 and 3 57 3 Green s functions in 2 and 3 Unlike the one dimensional case where Green s functions can be found explicitly for a number of different
More informationChapter 5.8: Bessel s equation
Chapter 5.8: Bessel s equation Bessel s equation of order ν is: x 2 y + xy + (x 2 ν 2 )y = 0. It has a regular singular point at x = 0. When ν = 0,, 2,..., this equation comes up when separating variables
More informationSolution of the Two-Dimensional Steady State Heat Conduction using the Finite Volume Method
Ninth International Conference on Computational Fluid Dynamics (ICCFD9), Istanbul, Turkey, July 11-15, 2016 ICCFD9-0113 Solution of the Two-Dimensional Steady State Heat Conduction using the Finite Volume
More informationFall Math 3410 Name (Print): Solution KEY Practice Exam 2 - November 4 Time Limit: 50 Minutes
Fall 206 - Math 340 Name (Print): Solution KEY Practice Exam 2 - November 4 Time Limit: 50 Minutes This exam contains pages (including this cover page) and 5 problems. Check to see if any pages are missing.
More informationNumerical simulation of a droplet icing on a cold surface
Numerical simulation of a droplet icing on a cold surface Buccini Saverio Peccianti Francesco Supervisors: Prof. Alessandro Bottaro Prof. Dominique Legendre March 30, 2017 Università degli Studi di Genova
More information