Grid adaptation for modeling trapezoidal difractors by using pseudo-spectral methods for solving Maxwell equations
|
|
- Alban Nash
- 6 years ago
- Views:
Transcription
1 Keldysh Institute Publication search Keldysh Institute preprints Preprint No., Zaitsev N.A., Sofronov I.L. Grid adaptation for modeling trapeoidal difractors by using pseudo-spectral methods for solving Maxwell equations ecommended form of bibliographic references: Zaitsev N.A., Sofronov I.L. Grid adaptation for modeling trapeoidal difractors by using pseudo-spectral methods for solving Maxwell equations. Keldysh Institute preprints,, No., 6 p. UL:
2 О р д е н а Л е н и н а ИНСТИТУТ ПРИКЛАДНОЙ МАТЕМАТИКИ имени М.В.Келдыша Р о с с и й с к о й а к а д е м и и н а у к N.A. Zaitsev, I.L. Sofronov Grid adaptation for modeling trapeoidal difractors by using pseudo-spectral methods for solving Maxwell equations Москва
3 N. A. Zaitsev, I. L. Sofronov, Grid adaptation for modeling trapeoidal difractors by using pseudo-spectral methods for solving Maxwell equations Abstract. We propose a series of successive one-dimensional coordinate transformations to construct an adaptive two-dimensional grid for the problem of diffraction on a grating with a trapeoidal tooth. Special attention is paid to the smoothness of the grid and grid refinement in domains of large solution gradients. Key words. Adaptive grid, Maxwell equation, difraction. Н. А. Зайцев, И. Л. Софронов, Построение адаптивной сетки для трапециевидных структур при моделировании дифракции псевдоспектральным методом решения уравнений Максвелла Аннотация. В работе предложен ряд последовательных одномерных преобразований координат для построения адаптивной двумерной сетки для задачи дифракции на решетке с трапециевидным зубом. Особое внимание уделено гладкости получаемой сетки и сгущению узлов к местам наибольших градиентов решения. Ключевые слова: адаптивная сетка, уравнение Максвелла, дифракция. Contents I. Preliminaries... II. Problem formulation... 4 III. Generation of adaptive -grid... 5 IV. Generation of adaptive x-grid... 7 V. Numerical examples... VI. eferences... 6 Работа выполнена при финансовой поддержке РФФИ, гранты 4 и -567.
4 I. Preliminaries We solve D diffraction problem described by Maxwell equations c c D H t H t E () with the dispersion law t D( r, t) E( r, t) ( ) E( r, t ) d, r ( x, y, ), P p ( t) p exp ( p ip) t exp ( p ip) t i. The scalar p-polariation case of (), normal incidence, reads: D H x c t D H c t x H E E c t x y y y x () Diffraction is considered on a grating with trapeoidal tooth. The pseudo-spectral method [] used for the solution of this problem requires high-quality grid in this periodical domain. We propose step-by-step method of generating the grid and make analysis of quality of the grid by considering various examples of values of geometric parameters.
5 4 II. Problem formulation Suppose we need to construct computational grid in D rectangular computational domain of width W and height H for a periodic with respect to x problem governed by (). The physical domain is the single spacing of a grating (W coincides with the period of the problem), see Figure. The height of the tooth of the grating is h g and its width is w g. The angle of the slope of sides of the tooth is. Let us introduce the following notations: x min, x max are the horiontal ends of the computational domain, W xmax xmin ; min, max are the vertical ends of the computational domain, H max min ; a is the base of the tooth of the grating; b is the top of the tooth of the grating; x( ), xl() are the right and left supporting lines, see Figure. For a b these lines coincide with the sides of the tooth. One has to continue the lines to the bottom and the top boundaries of the computational domain. They must have vertical tangents at the boundaries in order to use simpler transparent boundary conditions []. α wg hg x Figure. The general view of the computational domain.
6 5 III. Generation of adaptive -grid We will consider the symmetric case, i.e. when the tooth of the grating is an isosceles trapeoid. In this case xmax W /, xmin xmax ; max H /, min max ; h /, ; x ( ) x ( ) ; b g a b, x b wg hg x ( ) w cot( ) h / a g g L ( ) cot( ) /. In order to construct adaptive grid and calculate derivatives x,, such that ( ), x x(, ). use a one-to-one map / and / x we We construct the map as a superposition of simpler maps. Let us start from the ( ) transformation. max b x L() x () h g a min x min x max Figure. Supporting lines.
7 6 In -direction we have 4 supporting points: min, a, b, max. The grid in - direction should be Chebyshev s type grid and should have small grid spacing near the points a and b. First we normalie the domain in order to use sine function as a condensing map. So we need to map min on, a on, b on and max on. It is done by means of the polynomial function of fifth order: ( ) c c, c 5 c, [,] () where c5 max b/, c b c5. for max / b.7 ; and by cc ( ) c arctan( c ),, [,] ( c ) where c and c are chosen so that carctan( c) b, c arctan( c) max for max / b.7. (4) Now we concentrate grid points to the interfaces the following function: a and b by means of G G ( ) sin /, cos, [,] ; (5) where coefficient G is the approximate ratio of grid step sies near the interfaces to sies of the grid steps without this transformation. Usually one needs additional transformation in order to redistribute grid points between the inner region b and the outer region b. It can be done making use of the following function: where ( ) d d, d d, [,] (6) d, d d, here characteries a coarsening from the centre and is equal to the value of ratio ()/ () (the derivative d/ d at the middle point d ). If - is equal to
8 7 grid without transform (6) is equidistant then the ratio of grid steps at the middle point and at the outer boundary is equal to If then instead of coarsening the transform concentrates grid points to the center. However, if then map (6) is not one-to-one. In this case we can use the following map: where ( ) d sin( )/, d cos( ), [,] (7) d. Sometimes this transformation is not aggressive enough. Therefore we can use the following instead of (7): where t t( ), t, t ( t) t d sin( t)/, d cos( t), t [,], t( ) d sin( )/, t d cos( ), [,] d. (8) Now one can construct -grid using the Chebyshev s points cos ( N k)/( N ), k,..., N k and calculate derivative d/ d with the spectral accuracy, see []. IV. Generation of adaptive x-grid The essential difference between x-grid and -grid is that the problem is periodic in x- direction. Therefore all our maps must be periodic along lines const. We start the constructing of the x-grid from the eliminating of -dependence of the grid lines. For this purpose we use transform x( y, ) y asin( y/ ymax ), y[ xmax, xmax ], x a sin( y/ y ), x a cos( y/ y )/ y max y max max (9)
9 8 where x( ) y dx ( )/, d a a, sin( y / y ) sin( y / y ) max max max y x, y h /. g max () x and dx / We need also functions () a polynomial of the second order: d. The simplest way to continue x is using min a x ( ) a a a, dx / d a a, [, ], x ( ) x ( ) a, a a, a x ( ) a a, b a min a a a b a ( a min ) x( ) b b b, dx / d b b, [ b, max ], x ( ) x ( ) b, b b b a max b a ( b max ) One can use also polynomials of the third order: b b b, b x ( ) b b. min a x ( ) a a a a, dx / d a a a, [, ], x ( ) x ( ) a, a a, b a min b a ( a min ) a a min amin, a x( a) a a aa aa, b () x ( ) b b b b, dx / d b b b, [, max ], x( b) x( ) b a, b b max, ( ) and of the fourth order: b a b max max max b b b b b b b, b x ( ) b b b. 4 dx x ( ) 4, a a a a a a a a 4 a4, [ min, a], d x ( ) x ( ) a, a a ( ), a a 6 a, 4 b a 4 a min min 4 min b a ( a min ) a a min a min 4 a4min, a x( a) a a aa aa a4a 4, ()
10 9 4 dx x ( ) 4, b b b b b b b b 4 b4, [ b, max ], d x ( ) x ( ) b, b b ( ), b b 6 b, 4 b a 4 b max max 4 max b a ( b max ) b b max b max 4 b4 max, b x( b) b b bb b b b4b 4. () To provide higher smoothness one can need even the sixth order. Formulas are derived as follows. We introduce auxiliary map s min for [ min, a ]. Then x ( s) a a4s a5s a6s, dx / d 4a4s a5s 6 a6s, [ min, a] (4) where the coefficients are obtained from solution of system a a s a s a s x ( ) 4 a 5 a 6 a a 4 5 x ( ) ( ) b x a s a a a sa a6sa 4 4 a 5 a 6 a a4sa a5sa a6sa 6a s a s 5a s 5 b a (5) with and s a a min. Similarly, for b max [, ] s max x ( s) b b4s b5s b6s, dx / d 4b4 s b5 s 6 b6s, [ b, max ] (6) where the coefficients are solution of the system b 5 b 6 b b b b s b s b s x ( ) 4 5 x ( ) ( ) b x b s a b b sb b6s b 4 4 b 5 b 6 b b4sb b5s b b6sb 6b s b s 5b s 5 b a (7)
11 with sbb max. One can also use the following infinitely smooth continuation of sides of the trapeoid to the boundaries: C ( ) xinf C t xinf, min, a ; min a dx x( b) x( ) a erf yt ( ) ; d ( ) b a dx x ( ) ( ) x a d d where the values of the integral have to be computed numerically. The difference between the variants of x() is presented in Figure. x() and xl(), current order= a (8) -4 = a x ( a ) order x = order x = order x = 4 order x = 6 order x = = min Figure. Variants of x () obtained using polynomials of various orders in the region min a [, ] (-axis is vertical) The dependence of x() on parameter C xinf is as follows: if Cxinf the curve x() (red line in Figure 4) tends to the line BED, see Figure 4. In this figure the segment BC is the linear extrapolation of the right side of the trapeoid, the segment ED is the line joining midpoints of sides BC and AC of triangle ABC. If
12 Cxinf the curve x() tends to the segment BD. Numerical experiments show that the optimal value of C xinf is about. One can see that for all the values of parameter C xinf point D is the end point of line x(). If we would like to move the end from point D we can use the following function: min a C ( ) x inf C t x inf, min, a. min a The right side of the trapeoid B =a E =min A D C Figure 4. Lay-out of x () For the interval [ b, max ] an infinitely smooth continuation of the sides of the trapeoid is achieved by the following formulas: C ( ) xinf C t xinf, b, max ; max b dx x( b) x( ) a erf yt ( ) ; d ( ) b a dx x ( ) ( ) x b d. d b (9)
13 The normaliation of the domain y [ ymin, ymax ] we do by the same way as for - axis. Therefore, we us the map ymax y ( ) max d y d sin, y cos, [,] () where d y y max / for. ymax / y 5 (see ()); and the map y( ) ymax d exp d cos t dt, y d exp d cos, [,] () where parameters d and d must meet the following conditions d exp d cos d y y, for ymax / y. or ymax / y 5. max d exp d cos d y max The following map concentrates grid points to the interfaces x x() and x xl( ) : ( ) G sin /, G cos, [,] () x Approximately G x is the ratio of grid step sie near the interfaces to grid steps sie without this transformation. Usually one needs an additional transformation in order to redistribute grid points between the inner region x x() and the outer region x x(). It can be done by using the following function: x where t( ) d sin( )/, t d cos( ), [,], ( t) t d sin( t)/, d cos( t), t [,], t t t x d ; x ()
14 characteries the value of ratio ()/ (). If x-grid without transform () is x equidistant then the ratio of grid steps at the middle point and at the outer boundary is equal to x. However, for big values, if x, map () is not one-to-one. In this case we use the following map: ( ) d exp d cos t dt, d exp d cos, [,] (4) where parameters d and d must meet the following conditions: log, d x d Now we can construct equidistant -grid exp d cos d m/ N, m,..., N m x x ( Nx is excluded because of the periodicity) and calculate derivative d/ d with the spectral accuracy, see [].. V. Numerical examples In short notation the presented transform has the following form:,,, x x x x x x x x y,, y y,,. (5) Thus for a function u( x, ) we have:
15 4 u u x x, u u d u x x, d (6) where x x x x dy d d d d d d,, x x y d d d d d d d x da da dx ( ) sin( / max ), y y. d d d sin( y / y ) The corresponding Matlab code which computes the grid and implements formulas (6) has been written. For numerical tests we take the following function max x u( x, ) exp sin x max max (7) for which u u u x, exp cos x x x. (8) max max max max In order to check the spectral accuracy of the presented method and its implementation we take the following setup: H=, W=, h, w 8, 7 ; g g (9) and the following values of parameters: Cx inf, G., Gx., _ rarefaction.5, x_ rarefaction 8. () Table shows obtained accuracy of computed values of u/ and u/ x for various numbers of grid points.
16 5 Table. Error of derivatives versus discretiation parameters N x N Error of u/ Error of u/ x It s clear that the method has the spectral accuracy. Nevertheless, it seems that too many points needed to achieve the desired accuracy. The reason is that even a very smooth in x, -space function can has very big gradient in, -space. Let us investigate sensitivity to separate parameters of grid. Graph of function u( x, ) (7) in Cartesian coordinates is presented in Figure 5. Graph of the same function in parametric coordinates and for the following values of the parameters H, W, Gx, x, h 6, w 5, G, g g () is presented in Figure 6. Here hg H and wg W/ so the normaliation doesn t influence almost. Similar graph for Gx. is shown in Figure 7. Here u(, ) has much bigger gradient. Change of parameter G has similar effect: graph of u(, ) for parameters Gx and G. is shown in Figure 8 (other parameters are not changed). The next Figure shows the picture when we condense grid points to the interfaces in both directions: Gx., G.. Graph of u(, ) without any condensation but with coarsening x 8 is presented in Figure. The same picture with x and x /8 (concentration) is presented in Figure. The picture with x 8 and x /8 is presented in Figure. Graph of u(, ) without any condensation and coarsening but for the real positions of the interfaces ( hg, wg 8) is presented in Figure. One can see very strong deformation of the original function (see Figure 6). Graph of u(, ) when all the parameters are active is shown in Figure 4. One can see that the coarsening transform reduces the gradient near the boundaries but generates additional gradient near the center. Next set of figures shows graphs of correspondent grids.
17 6 In Figure 5 and Figure 6 -grid and its grid steps are shown when there is no coarsening in direction. Figure 7 and Figure 8 show the same picture for coarsening.5. Figure 9 and Figure show x-grid grid and its grid steps when there is no coarsening in x direction. Figure Figure show the same picture for coarsening x 8. The final grid when all the parameters are active is presented in Figure. VI. eferences [] Зайцев Н.А., Софронов И.Л. Схемы высокой точности для трёхмерных уравнений Максвелла в средах Лоренца, базирующиеся на неявных схемах переменных направлений // Препринты ИПМ им.м.в.келдыша с. UL: [] Софронов И.Л., Воскобойникова О.И., Зайцев Н.А. Прозрачные граничные условия на неравномерных сетках для двумерных уравнений Максвелла с дисперсиейсхемы высокой точности для трёхмерных уравнений Максвелла в средах Лоренца, базирующиеся на неявных схемах переменных направлений // Препринты ИПМ им.м.в.келдыша с. UL: [] Lloyd N. Trefethen, Spectral Methods in MATLAB, SIAM, Philadelphia,
18 x -5 Figure 5. Function u( x, ) Figure 6. Function u(, ) for parameters of grid transform H, h 6, W, w 5, G, G,,. g g x x
19 Figure 7. Function u(, ) for parameters of grid transform H, h 6, W, w 5, G., G,,. g g x x Figure 8. Function u(, ) for parameters of grid transform H, h 6, W, w 5, G, G.,,. g g x x
20 Figure 9. Function u(, ) for parameters of grid transform H, h 6, W, w 5, G., G.,,. g g x x Figure. Function u(, ) for parameters of grid transform H, h 6, W, w 5, G, G, 8,. g g x x
21 Figure. Function u(, ) for parameters of grid transform H, h 6, W, w 5, G, G,,.5 g g x x Figure. Function u(, ) for parameters of grid transform H, h 6, W, w 5, G, G, 8,.5 g g x x
22 Figure. Function u(, ) for parameters of grid transform H, h, W, w 8, G, G,,. g g x x Figure 4. Function u(, ) for parameters of grid transform H, h, W, w 8, G., G., 8,.5 g g x x
23 6 -grid, K_intervals= Figure 5. -grid when there is no coarsening in direction, 6 Grid steps: h min =.476 h max = h Figure 6. Grid steps () h when there is no coarsening in direction,
24 6 -grid, K_intervals= Figure 7. -grid for coarsening.5 6 Grid steps: h min =.594 h max = h Figure 8. Grid steps () h for coarsening.5
25 4 6 x-grid, Kx_intervals= Figure 9. x-grid when there is no coarsening in x direction, x 4 Grid steps: hx min =.854 hx max = h x x Figure. Grid steps () x h x when there is no coarsening in x direction, x
26 5 6 x-grid, Kx_intervals= Figure. x-grid for coarsening x 8 4 Grid steps: hx min =.6975 hx max = h x x Figure. Grid steps () x h x for coarsening 8 x
27 6 6 D grid, Kx_intervals=64 K_intervals= Figure. Final grid for H, h, W, w 8, G., G., 8,.5 g x x g
Study of the accuracy of active magnetic damping algorithm
Keldysh Institute Publication search Keldysh Institute preprints Preprint No 6 ISSN 7-898 (Print) ISSN 7-9 (Online) Ovchinnikov MYu Roldugin DS Tkachev SS Study of the accuracy of active magnetic damping
More informationMAT137 Calculus! Lecture 6
MAT137 Calculus! Lecture 6 Today: 3.2 Differentiation Rules; 3.3 Derivatives of higher order. 3.4 Related rates 3.5 Chain Rule 3.6 Derivative of Trig. Functions Next: 3.7 Implicit Differentiation 4.10
More informationMath-Net.Ru All Russian mathematical portal
Math-Net.Ru All Russian mathematical portal L. D. Pustyl nikov, M. V. Deryabin, Nonequilibrium gas, entropy and generalized billiards, Keldysh Institute preprints, 2014, 003 Use of the all-russian mathematical
More informationNational Quali cations
H 08 X747/76/ National Quali cations Mathematics Paper (Non-Calculator) THURSDAY, MAY 9:00 AM 0:0 AM Total marks 60 Attempt ALL questions. You may NOT use a calculator. Full credit will be given only to
More informationMath-Net.Ru All Russian mathematical portal
Math-Net.Ru All Russian mathematical portal A. B. Batkhin, Families of symmetric periodic solutions of the Generalized Hill s Problem, Keldysh Institute preprints, 2013, 060 Use of the all-russian mathematical
More informationchapter 1 vector geometry solutions V Consider the parallelogram shown alongside. Which of the following statements are true?
chapter vector geometry solutions V. Exercise A. For the shape shown, find a single vector which is equal to a)!!! " AB + BC AC b)! AD!!! " + DB AB c)! AC + CD AD d)! BC + CD!!! " + DA BA e) CD!!! " "
More informationAP Calculus Testbank (Chapter 9) (Mr. Surowski)
AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series
More informationDifferential equations for the radial limits in Z. solutions of a discrete integrable system. 2 + of the. Keldysh Institute Publication search
Keldysh Institute Publication search Keldysh Institute preprints Preprint No. 214, 2018 ISSN 2071-2898 (Print) ISSN 2071-2901 (Online) Aptekarev A.I., Kozhan R. Differential equations for the radial limits
More informationCO-ORDINATE GEOMETRY. 1. Find the points on the y axis whose distances from the points (6, 7) and (4,-3) are in the. ratio 1:2.
UNIT- CO-ORDINATE GEOMETRY Mathematics is the tool specially suited for dealing with abstract concepts of any ind and there is no limit to its power in this field.. Find the points on the y axis whose
More information2.2 The derivative as a Function
2.2 The derivative as a Function Recall: The derivative of a function f at a fixed number a: f a f a+h f(a) = lim h 0 h Definition (Derivative of f) For any number x, the derivative of f is f x f x+h f(x)
More informationChapter 10 Conics, Parametric Equations, and Polar Coordinates Conics and Calculus
Chapter 10 Conics, Parametric Equations, and Polar Coordinates 10.1 Conics and Calculus 1. Parabola A parabola is the set of all points x, y ( ) that are equidistant from a fixed line and a fixed point
More informationTopic 3 Part 1 [449 marks]
Topic 3 Part [449 marks] a. Find all values of x for 0. x such that sin( x ) = 0. b. Find n n+ x sin( x )dx, showing that it takes different integer values when n is even and when n is odd. c. Evaluate
More information6x 2 8x + 5 ) = 12x 8
Example. If f(x) = x 3 4x + 5x + 1, then f (x) = 6x 8x + 5 Observation: f (x) is also a differentiable function... d dx ( f (x) ) = d dx ( 6x 8x + 5 ) = 1x 8 The derivative of f (x) is called the second
More information(1) Recap of Differential Calculus and Integral Calculus (2) Preview of Calculus in three dimensional space (3) Tools for Calculus 3
Math 127 Introduction and Review (1) Recap of Differential Calculus and Integral Calculus (2) Preview of Calculus in three dimensional space (3) Tools for Calculus 3 MATH 127 Introduction to Calculus III
More informationSection I 10 marks (pages 2 5) Attempt Questions 1 10 Allow about 15 minutes for this section
017 HIGHER SCHOOL CERTIFICATE EXAMINATION Mathematics General Instructions Reading time 5 minutes Working time hours Write using black pen NESA approved calculators may be used A reference sheet is provided
More informationTABLE OF CONTENTS 2 CHAPTER 1
TABLE OF CONTENTS CHAPTER 1 Quadratics CHAPTER Functions 3 CHAPTER 3 Coordinate Geometry 3 CHAPTER 4 Circular Measure 4 CHAPTER 5 Trigonometry 4 CHAPTER 6 Vectors 5 CHAPTER 7 Series 6 CHAPTER 8 Differentiation
More information. State whether the triangle is an isosceles triangle, a right triangle, neither of these, or both.
PRECALCULUS MIDTERM PRACTICE TEST (2008-2009) Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Find the distance d(, ) between the points and.
More informationТЕКА Кom. Mot. i Energ. Roln. OL PAN, 2010, 10B, THEORETICAL ANALYSIS OF A CIRCLE CURRENT CONVERTER OPERATION UNDER CONTROL OF HOLLOW CYLINDERS
ТЕКА Кom. Mot. i Energ. Roln. OL PAN, 2010, 10B, 33-42 THEORETICAL ANALYSIS OF A CIRCLE CURRENT CONVERTER OPERATION UNDER CONTROL OF HOLLOW CYLINDERS Vadim Miroshnikov, Oleivi Anver Volodymyr Dal East-Ukrainian
More information2016 SEC 4 ADDITIONAL MATHEMATICS CW & HW
FEB EXAM 06 SEC 4 ADDITIONAL MATHEMATICS CW & HW Find the values of k for which the line y 6 is a tangent to the curve k 7 y. Find also the coordinates of the point at which this tangent touches the curve.
More informationPart (1) Second : Trigonometry. Tan
Part (1) Second : Trigonometry (1) Complete the following table : The angle Ratio 42 12 \ Sin 0.3214 Cas 0.5321 Tan 2.0625 (2) Complete the following : 1) 46 36 \ 24 \\ =. In degrees. 2) 44.125 = in degrees,
More informationYou can learn more about the services offered by the teaching center by visiting
MAC 232 Exam 3 Review Spring 209 This review, produced by the Broward Teaching Center, contains a collection of questions which are representative of the type you may encounter on the exam. Other resources
More information2015 Fall Startup Event Solutions
1. Evaluate: 829 7 The standard division algorithm gives 1000 + 100 + 70 + 7 = 1177. 2. What is the remainder when 86 is divided by 9? Again, the standard algorithm gives 20 + 1 = 21 with a remainder of
More informationMATH 2433 Homework 1
MATH 433 Homework 1 1. The sequence (a i ) is defined recursively by a 1 = 4 a i+1 = 3a i find a closed formula for a i in terms of i.. In class we showed that the Fibonacci sequence (a i ) defined by
More informationMAT1035 Analytic Geometry
MAT1035 Analytic Geometry Lecture Notes R.A. Sabri Kaan Gürbüzer Dokuz Eylül University 2016 2 Contents 1 Review of Trigonometry 5 2 Polar Coordinates 7 3 Vectors in R n 9 3.1 Located Vectors..............................................
More informationThe Princeton Review AP Calculus BC Practice Test 2
0 The Princeton Review AP Calculus BC Practice Test CALCULUS BC SECTION I, Part A Time 55 Minutes Number of questions 8 A CALCULATOR MAY NOT BE USED ON THIS PART OF THE EXAMINATION Directions: Solve each
More informationMATHEMATICS OF DOPING PROFILES. C(x,t) t. = D 2 C(x,t) x 2. 4Dt dx '
EE43 MATHEMATICS OF DOPING PROFILES N. Cheung The diffusion equation with constant D : has the general solution: C(x,t) = C(x,t) = D 2 C(x,t) 4πDt F(x ' ) e -(x-x' ) 2 4Dt dx ' - where F(x') is the C(x,t)
More information13 Path Planning Cubic Path P 2 P 1. θ 2
13 Path Planning Path planning includes three tasks: 1 Defining a geometric curve for the end-effector between two points. 2 Defining a rotational motion between two orientations. 3 Defining a time function
More information0, otherwise. Find each of the following limits, or explain that the limit does not exist.
Midterm Solutions 1, y x 4 1. Let f(x, y) = 1, y 0 0, otherwise. Find each of the following limits, or explain that the limit does not exist. (a) (b) (c) lim f(x, y) (x,y) (0,1) lim f(x, y) (x,y) (2,3)
More informationMu Alpha Theta National Convention 2013
Practice Round Alpha School Bowl P1. What is the common difference of the arithmetic sequence 10, 23,? P2. Find the sum of the digits of the base ten representation of 2 15. P3. Find the smaller value
More informationSenior Secondary Australian Curriculum
Senior Secondary Australian Curriculum Specialist Mathematics Glossary Unit 1 Combinatorics Arranging n objects in an ordered list The number of ways to arrange n different objects in an ordered list is
More informationAnnouncements. Topics: Homework: - sections 4.5 and * Read these sections and study solved examples in your textbook!
Announcements Topics: - sections 4.5 and 5.1-5.5 * Read these sections and study solved examples in your textbook! Homework: - review lecture notes thoroughly - work on practice problems from the textbook
More informationSec 4 Maths. SET A PAPER 2 Question
S4 Maths Set A Paper Question Sec 4 Maths Exam papers with worked solutions SET A PAPER Question Compiled by THE MATHS CAFE 1 P a g e Answer all the questions S4 Maths Set A Paper Question Write in dark
More informationb) The trend is for the average slope at x = 1 to decrease. The slope at x = 1 is 1.
Chapters 1 to 8 Course Review Chapters 1 to 8 Course Review Question 1 Page 509 a) i) ii) [2(16) 12 + 4][2 3+ 4] 4 1 [2(2.25) 4.5+ 4][2 3+ 4] 1.51 = 21 3 = 7 = 1 0.5 = 2 [2(1.21) 3.3+ 4][2 3+ 4] iii) =
More information(D) (A) Q.3 To which of the following circles, the line y x + 3 = 0 is normal at the point ? 2 (A) 2
CIRCLE [STRAIGHT OBJECTIVE TYPE] Q. The line x y + = 0 is tangent to the circle at the point (, 5) and the centre of the circles lies on x y = 4. The radius of the circle is (A) 3 5 (B) 5 3 (C) 5 (D) 5
More informationAPPM 1360 Final Exam Spring 2016
APPM 36 Final Eam Spring 6. 8 points) State whether each of the following quantities converge or diverge. Eplain your reasoning. a) The sequence a, a, a 3,... where a n ln8n) lnn + ) n!) b) ln d c) arctan
More informationLesson 5 Practice Problems
Name: Date: Lesson 5 Section 5.1: Linear Functions vs. Exponential Functions 1. Complete the table below. Function Linear or Exponential? Linear: Increasing or Decreasing? Exponential: Growth or Decay?
More information2008 Euclid Contest. Solutions. Canadian Mathematics Competition. Tuesday, April 15, c 2008 Centre for Education in Mathematics and Computing
Canadian Mathematics Competition An activity of the Centre for Education in Mathematics and Computing, University of Waterloo, Waterloo, Ontario 008 Euclid Contest Tuesday, April 5, 008 Solutions c 008
More informationMathematics, Algebra, and Geometry
Mathematics, Algebra, and Geometry by Satya http://www.thesatya.com/ Contents 1 Algebra 1 1.1 Logarithms............................................ 1. Complex numbers........................................
More informationJEE MAIN 2013 Mathematics
JEE MAIN 01 Mathematics 1. The circle passing through (1, ) and touching the axis of x at (, 0) also passes through the point (1) (, 5) () (5, ) () (, 5) (4) ( 5, ) The equation of the circle due to point
More informationSolutions to Homework 11
Solutions to Homework 11 Read the statement of Proposition 5.4 of Chapter 3, Section 5. Write a summary of the proof. Comment on the following details: Does the proof work if g is piecewise C 1? Or did
More informationChapter 3a Topics in differentiation. Problems in differentiation. Problems in differentiation. LC Abueg: mathematical economics
Chapter 3a Topics in differentiation Lectures in Mathematical Economics L Cagandahan Abueg De La Salle University School of Economics Problems in differentiation Problems in differentiation Problem 1.
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES We have: Seen how to interpret derivatives as slopes and rates of change Seen how to estimate derivatives of functions given by tables of values Learned how
More informationx n+1 = ( x n + ) converges, then it converges to α. [2]
1 A Level - Mathematics P 3 ITERATION ( With references and answers) [ Numerical Solution of Equation] Q1. The equation x 3 - x 2 6 = 0 has one real root, denoted by α. i) Find by calculation the pair
More informationUnit 8. ANALYTIC GEOMETRY.
Unit 8. ANALYTIC GEOMETRY. 1. VECTORS IN THE PLANE A vector is a line segment running from point A (tail) to point B (head). 1.1 DIRECTION OF A VECTOR The direction of a vector is the direction of the
More information6x 2 8x + 5 ) = 12x 8. f (x) ) = d (12x 8) = 12
AMS/ECON 11A Class Notes 11/6/17 UCSC *) Higher order derivatives Example. If f = x 3 x + 5x + 1, then f = 6x 8x + 5 Observation: f is also a differentiable function... d f ) = d 6x 8x + 5 ) = 1x 8 dx
More information2. Maxwell s Equations
UDC 519.63 Parallel Second Order Finite Volume Scheme for Maxwell s Equations with Discontinuous Dielectric Permittivity on Structured Meshes T. Z. Ismagilov Department of Information Technology Novosibirsk
More informationPOINT. Preface. The concept of Point is very important for the study of coordinate
POINT Preface The concept of Point is ver important for the stud of coordinate geometr. This chapter deals with various forms of representing a Point and several associated properties. The concept of coordinates
More informationHodograph transformations and generating of solutions for nonlinear differential equations
W.I. Fushchych Scientific Works 23 Vol. 5 5. Hodograph transformations and generating of solutions for nonlinear differential equations W.I. FUSHCHYCH V.A. TYCHYNIN Перетворення годографа однiєї скалярної
More informationLecture Wise Questions from 23 to 45 By Virtualians.pk. Q105. What is the impact of double integration in finding out the area and volume of Regions?
Lecture Wise Questions from 23 to 45 By Virtualians.pk Q105. What is the impact of double integration in finding out the area and volume of Regions? Ans: It has very important contribution in finding the
More informationSlide 1. Slide 2. Slide 3 Remark is a new function derived from called derivative. 2.2 The derivative as a Function
Slide 1 2.2 The derivative as a Function Slide 2 Recall: The derivative of a function number : at a fixed Definition (Derivative of ) For any number, the derivative of is Slide 3 Remark is a new function
More informationMO-ARML Practice March 17, 2007 Relay Questions
O-ARL Practice arch 17, 007 Relay Questions If a question depends on NYWR, the answer should be expressed in terms of the variable representing NYWR 1. Compute the number of distinct 11 letter palindromic
More informationNational Quali cations
H 2018 X747/76/11 National Quali cations Mathematics Paper 1 (Non-Calculator) THURSDAY, 3 MAY 9:00 AM 10:10 AM Total marks 60 Attempt ALL questions. You may NOT use a calculator. Full credit will be given
More informationQUESTION BANK ON. CONIC SECTION (Parabola, Ellipse & Hyperbola)
QUESTION BANK ON CONIC SECTION (Parabola, Ellipse & Hyperbola) Question bank on Parabola, Ellipse & Hyperbola Select the correct alternative : (Only one is correct) Q. Two mutually perpendicular tangents
More informationCLASS X FORMULAE MATHS
Real numbers: Euclid s division lemma Given positive integers a and b, there exist whole numbers q and r satisfying a = bq + r, 0 r < b. Euclid s division algorithm: This is based on Euclid s division
More informationMATHEMATICAL METHODS UNIT 1 Chapter 1 Reviewing Linear Equations Chapter 2 Coordinate geometry & linear relations
REVIEWING LINEAR EQUATIONS E da = q ε ( B da = 0 E ds = dφ. B ds = μ ( i + μ ( ε ( dφ 3 MATHEMATICAL METHODS UNIT 1 Chapter 1 Reviewing Linear Equations Chapter 2 Coordinate geometry & linear relations
More information9-12 Mathematics Vertical Alignment ( )
Algebra I Algebra II Geometry Pre- Calculus U1: translate between words and algebra -add and subtract real numbers -multiply and divide real numbers -evaluate containing exponents -evaluate containing
More informationParametric Equations and Polar Coordinates
Parametric Equations and Polar Coordinates Parametrizations of Plane Curves In previous chapters, we have studied curves as the graphs of functions or equations involving the two variables x and y. Another
More informationSection 3.7. Rolle s Theorem and the Mean Value Theorem
Section.7 Rolle s Theorem and the Mean Value Theorem The two theorems which are at the heart of this section draw connections between the instantaneous rate of change and the average rate of change of
More informationIntegration. Tuesday, December 3, 13
4 Integration 4.3 Riemann Sums and Definite Integrals Objectives n Understand the definition of a Riemann sum. n Evaluate a definite integral using properties of definite integrals. 3 Riemann Sums 4 Riemann
More informationMATH 255 Applied Honors Calculus III Winter Midterm 1 Review Solutions
MATH 55 Applied Honors Calculus III Winter 11 Midterm 1 Review Solutions 11.1: #19 Particle starts at point ( 1,, traces out a semicircle in the counterclockwise direction, ending at the point (1,. 11.1:
More informationNumerical-analytical algorithm for constructing the envelope of the projectile trajectories in midair
Numerical-analytical algorithm for constructing the envelope of the projectile trajectories in midair Peter S Chudinov Perm State Agricultural Academy, Perm, 61499, 5-th October st., 1, Russia E-mail:
More informationUnit 10 Parametric and Polar Equations - Classwork
Unit 10 Parametric and Polar Equations - Classwork Until now, we have been representing graphs by single equations involving variables x and y. We will now study problems with which 3 variables are used
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 informationChapter 1: Introduction
Chapter 1: Introduction Definition: A differential equation is an equation involving the derivative of a function. If the function depends on a single variable, then only ordinary derivatives appear and
More informationPractice Assessment Task SET 3
PRACTICE ASSESSMENT TASK 3 655 Practice Assessment Task SET 3 Solve m - 5m + 6 $ 0 0 Find the locus of point P that moves so that it is equidistant from the points A^-3, h and B ^57, h 3 Write x = 4t,
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 information1. The positive zero of y = x 2 + 2x 3/5 is, to the nearest tenth, equal to
SAT II - Math Level Test #0 Solution SAT II - Math Level Test No. 1. The positive zero of y = x + x 3/5 is, to the nearest tenth, equal to (A) 0.8 (B) 0.7 + 1.1i (C) 0.7 (D) 0.3 (E). 3 b b 4ac Using Quadratic
More informationAP Calculus (BC) Chapter 9 Test No Calculator Section Name: Date: Period:
WORKSHEET: Series, Taylor Series AP Calculus (BC) Chapter 9 Test No Calculator Section Name: Date: Period: 1 Part I. Multiple-Choice Questions (5 points each; please circle the correct answer.) 1. The
More information( ) ( ) SECTION 1.1, Page ( x 3) 5 = 4( x 5) = 7. x = = = x x+ 0.12(4000 x) = 432
CHAPTER Functions and Graphs SECTION., Page. x + x + x x x. x + x x x x x. ( x ) ( x ) x 6 x x x x x + x x 7. x + x + x + 6 8 x 8 6 x x. x x 6 x 6 x x x 8 x x 8 + x..x +..6.x. x 6 ( n + ) ( n ) n + n.
More informationKepler s laws. M. Plakhotnyk. March 2, 2017
Kepler s laws arxiv:1702.06537v2 [math.ho] 1 Mar 2017 M. Plakhotnyk March 2, 2017 Kepler s laws with introduction to differential calculus (book for scholars, who are interested in physics and mathematics)
More informationMAC Calculus II Spring Homework #6 Some Solutions.
MAC 2312-15931-Calculus II Spring 23 Homework #6 Some Solutions. 1. Find the centroid of the region bounded by the curves y = 2x 2 and y = 1 2x 2. Solution. It is obvious, by inspection, that the centroid
More informationPre-Calculus Team Questions Florida Regional Competition March C = ( )
Florida Regional Competition March 08. Given: sin ( ) sin α = for 0 < α
More informationHigher Mathematics Course Notes
Higher Mathematics Course Notes Equation of a Line (i) Collinearity: (ii) Gradient: If points are collinear then they lie on the same straight line. i.e. to show that A, B and C are collinear, show that
More informationPossible C2 questions from past papers P1 P3
Possible C2 questions from past papers P1 P3 Source of the original question is given in brackets, e.g. [P1 January 2001 Question 1]; a question which has been edited is indicated with an asterisk, e.g.
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 informationDEPARTMENT OF MATHEMATICS
DEPARTMENT OF MATHEMATICS AS level Mathematics Core mathematics 1 C1 2015-2016 Name: Page C1 workbook contents Indices and Surds Simultaneous equations Quadratics Inequalities Graphs Arithmetic series
More informationMATHEMATICS 2017 HSC Course Assessment Task 4 (Trial Examination) Thursday, 3 August 2017
MATHEMATICS 2017 HSC Course Assessment Task 4 (Trial Examination) Thursday, 3 August 2017 General instructions Working time 3 hours. (plus 5 minutes reading time) Write using blue or black pen. Where diagrams
More informationExercise: concepts from chapter 3
Reading:, Ch 3 1) The natural representation of a curve, c = c(s), satisfies the condition dc/ds = 1, where s is the natural parameter for the curve. a) Describe in words and a sketch what this condition
More informationChapter 8: Taylor s theorem and L Hospital s rule
Chapter 8: Taylor s theorem and L Hospital s rule Theorem: [Inverse Mapping Theorem] Suppose that a < b and f : [a, b] R. Given that f (x) > 0 for all x (a, b) then f 1 is differentiable on (f(a), f(b))
More informationQUESTION BANK ON STRAIGHT LINE AND CIRCLE
QUESTION BANK ON STRAIGHT LINE AND CIRCLE Select the correct alternative : (Only one is correct) Q. If the lines x + y + = 0 ; 4x + y + 4 = 0 and x + αy + β = 0, where α + β =, are concurrent then α =,
More informationCreated by T. Madas. Candidates may use any calculator allowed by the regulations of this examination.
IYGB GCE Mathematics SYN Advanced Level Snoptic Paper C Difficult Rating: 3.895 Time: 3 hours Candidates ma use an calculator allowed b the regulations of this eamination. Information for Candidates This
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 informationComplex Integration Line Integral in the Complex Plane CHAPTER 14
HAPTER 14 omplex Integration hapter 13 laid the groundwork for the study of complex analysis, covered complex numbers in the complex plane, limits, and differentiation, and introduced the most important
More informationMath-Net.Ru All Russian mathematical portal
Math-Net.Ru All Russian mathematical portal U. V. Linnik, On the representation of large numbers as sums of seven cubes, Rec. Math. [Mat. Sbornik] N.S., 1943, Volume 12(54), Number 2, 218 224 Use of the
More informationAS PURE MATHS REVISION NOTES
AS PURE MATHS REVISION NOTES 1 SURDS A root such as 3 that cannot be written exactly as a fraction is IRRATIONAL An expression that involves irrational roots is in SURD FORM e.g. 2 3 3 + 2 and 3-2 are
More informationStorm Open Library 3.0
S 50% off! 3 O L Storm Open Library 3.0 Amor Sans, Amor Serif, Andulka, Baskerville, John Sans, Metron, Ozdoby,, Regent, Sebastian, Serapion, Splendid Quartett, Vida & Walbaum. d 50% f summer j sale n
More informationIYGB. Special Paper U. Time: 3 hours 30 minutes. Created by T. Madas. Created by T. Madas
IYGB Special Paper U Time: 3 hours 30 minutes Candidates may NOT use any calculator Information for Candidates This practice paper follows the Advanced Level Mathematics Core Syllabus Booklets of Mathematical
More information2001 Solutions Euclid Contest (Grade 12)
Canadian Mathematics Competition An activity of The Centre for Education in Mathematics and Computing, University of Waterloo, Waterloo, Ontario 001 s Euclid Contest (Grade 1) for The CENTRE for EDUCATION
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 Partial Differentiation Page y f x
ENGI 344 4 Partial Differentiation Page 4-0 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can
More informationLearning Objectives These show clearly the purpose and extent of coverage for each topic.
Preface This book is prepared for students embarking on the study of Additional Mathematics. Topical Approach Examinable topics for Upper Secondary Mathematics are discussed in detail so students can focus
More informationKinematics
Classical Mechanics Kinematics Velocities Idea 1 Choose the appropriate frame of reference. Useful frames: - Bodies are at rest - Projections of velocities vanish - Symmetric motion Change back to the
More informationSection 3.5: Implicit Differentiation
Section 3.5: Implicit Differentiation In the previous sections, we considered the problem of finding the slopes of the tangent line to a given function y = f(x). The idea of a tangent line however is not
More informationCOORDINATE GEOMETRY BASIC CONCEPTS AND FORMULAE. To find the length of a line segment joining two points A(x 1, y 1 ) and B(x 2, y 2 ), use
COORDINATE GEOMETRY BASIC CONCEPTS AND FORMULAE I. Length of a Line Segment: The distance between two points A ( x1, 1 ) B ( x, ) is given b A B = ( x x1) ( 1) To find the length of a line segment joining
More informationChapter 10. Additional Topics in Trigonometry
Chapter 10 Additional Topics in Trigonometry 1 Right Triangle applications Law of Sines and Cosines Parametric equations Polar coordinates Curves in polar coordinates Summary 2 Chapter 10.1 Right Triangle
More informationUNIVERSITY OF HOUSTON HIGH SCHOOL MATHEMATICS CONTEST Spring 2018 Calculus Test
UNIVERSITY OF HOUSTON HIGH SCHOOL MATHEMATICS CONTEST Spring 2018 Calculus Test NAME: SCHOOL: 1. Let f be some function for which you know only that if 0 < x < 1, then f(x) 5 < 0.1. Which of the following
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 informationAP Calculus BC : The Fundamental Theorem of Calculus
AP Calculus BC 415 5.3: The Fundamental Theorem of Calculus Tuesday, November 5, 008 Homework Answers 6. (a) approimately 0.5 (b) approimately 1 (c) approimately 1.75 38. 4 40. 5 50. 17 Introduction In
More informationPreliminary Mathematics
NORTH SYDNEY GIRLS HIGH SCHOOL 2011 YEARLY EXAMINATION Preliminary Mathematics General Instructions Reading Time 5 minutes Working Time 2 hours Write using black or blue pen Board-approved calculators
More informationAdd Math (4047/02) Year t years $P
Add Math (4047/0) Requirement : Answer all questions Total marks : 100 Duration : hour 30 minutes 1. The price, $P, of a company share on 1 st January has been increasing each year from 1995 to 015. The
More information