The exact region of stability for MacCormack scheme. Hoon Hong. Research Institute for Symbolic Computation. Johannes Kepler University
|
|
- Britton Reynard Walsh
- 6 years ago
- Views:
Transcription
1 The eact region of stailit for MacCormack scheme Hoon Hong Research Institute for Smolic Computation Johannes Kepler Universit A- Linz, Austria Astract Let the two dimensional scalar avection equation @u = ^a + ^ : We prove that the stailit region @ MacCormack scheme for this equation is eactl given ^a t t + ^ 1 where t ; and are the grid distances. It is interesting to note that the stailit region is identical to the one for La-Wendro scheme proved Turkel. 1 Introduction 1 a + a + = 1=8 a = + = = 1 Stailit analsis is one of the most fundamental tasks in devising reliale numerical methods for solving partial dierential equations. In this paper, we investigate the stailit of MacCormack's scheme [8] which is widel used in solving hperolic partial dierential equations, especiall in aerodnamics. In particular we consider the following scalar two dimensional = : In [1] Wendro derived the following sucient condition for the stailit region: a + a + 1=8: where a and are respectivel ^a t and ^ t, where again t ; and are the grid distances. As it was remarked [1], numerical sampling shows that it is not a necessar condition. The main result of this paper is the derivation of the following necessar and sucient condition: a + Figure 1 compares the region otained Wendro and the one otained here. It is interesting to note that the eact stailit region of MacCormack scheme derived here is identical to that for La- Wendro scheme proved Turkel [11]. In fact, one will notice in the susequent sections that the underling proof techniques are similar. This research was done in the framework of the European project ACCLAIM Figure 1: Stailit Regions The structure of the paper is as follows. In the net section, we will give a precise description of the stailit prolem and the eact solution for MacCormack scheme. In the following section, we will give a proof. Stailit of MacCormack Scheme In this section, we recall the denitions of the La-Wendro scheme and the MacCormack scheme. Then we give a precise statement of the stailit prolem for the MacCormack scheme and the solution for it. Let T and T e the translation operators respectivel de- ned T u(; ; t) = u( + ; ; t); T u(; ; t) = u(; + ; t): Let L = I + 1 a(t? T?1 ) + 1 (T? T?1 ) + 1 a (T? I + T?1 ) + 1 (T? I + T?1 ) + 1 a(t? T?1 )(T? T?1 ): 1
2 Then the La-Wendro [6] scheme is given Let u(; ; t + t ) Lu(; ; t): G ++ = I + a(t? I) + (T? I); G?? = I + a(i? T?1 ) + (I? T?1 ); G?+ = I + a(i? T?1 ) + (T? I); G +? = I + a(t? I) + (I? T?1); M 1 = 1 (I + G ++G?? ); M = 1 (I + G?+G +? ); M = M M 1 : Then the MacCormack scheme [8] is given u(; ; t + t ) Mu(; ; t): As oserved in [1], the two operators M and L are related M = L? Q where Q = 1 a(t? I + T?1 )(T? I + T?1 ): In order to otain the amplications for the La-Wendro scheme and the MacCormack scheme, one replaces T and T with e i and e i respectivel, otaining r = 1? a (1? cos )? (1? cos )? a sin sin s = a sin + sin q = a(1? cos )(1? cos ) l = r + is m = l? q where l and m are respectivel the amplications of the La- Wendro and the MacCormack scheme. Now the stailit prolem for the MacCormack scheme is to nd a necessar and sucient condition on a and such that for all real values of and the asolute value of m is not greater than 1. More formall, the prolem is to nd a quantier-free formula F (a; ) which is equivalent to the quantied formula 8(; ) R jmj 1: Let us oserve that we can remove the trigonometric functions from the aove formula carring out the usual parameterization with u 1 = sin, v 1 = cos, u = sin and v = cos, along with the condition u 1 + v 1 = u + v = 1. The resulting formula (involving onl algeraic operations) elongs to the language of the rst order theor of real closed elds. It is well known that the quantier elimination in that theor can e carried out algorithmicall [1]. Further, there have een various reakthroughs in devising more ecient algorithms [1,,, 5,, 9]. There also eists an implementation []. Thus, in principle, the stailit prolems can e solved completel automaticall a computer program. In fact, there have een interesting and important progress in this direction as reported in [7] where Liska and Steinerg succeeded in deriving the stailit conditions for several prolems using the quantier elimination program qepcad []. But, in spite of these various progresses, the complete automation for the MacCormack scheme is not et feasile due to the huge requirement of computing times. Thus, the result reported in this paper is ased on \human-rain" proof. But it must e emphasized that during conjecturing the main result, computer plaed essential role, as will e reported in some other paper. Now we state the main result of this paper, that we will prove in the net section. Theorem 1 (Stailit of MacCormack Scheme) 8(; ) R jmj 1 a + Proof In this section, we give a proof of the theorem stated aove. The proof will consists of proofs of four lemmas. The dependenc among the lemmas is as follows. The rst lemma is used for proving the second lemma. The second lemma is used for proving the third and the fourth lemma. The third lemma proves that the proposed stailit condition is necessar. The fourth lemma proves that it is also sucient. Slightl dierent variations of the rst two lemmas were proved Turkel [11]. Thus, the reader who is familiar with [11] might want to skip reading the proofs of the two lemmas. We decided to include the complete proof here ecause our lemmas are slight dierent from Turkel's in that ours do not have asolute value function in the formulation, and thus require careful adjustments in the proofs. Lemma 1 (Univariate quantier elimination) Let P () = (a + )? a? : 8 R P () a + Proof. The proof will e case analsis. First three cases take care of the \degenerate" situations. Case: a = L.H.S. 8 R? Case: =? (? 1) L.H.S. 8 R a? a 8 R a (a? 1)
3 Case: a = 1 a (a? 1) a 1 a + 1 L.H.S. 8 R + 6 +? = (1) + 1 Case: Otherwise, that is, a 6= ; a 6= 1; 6=. Note that P () = (a? a ) + : Since (a? a ) is not zero, P is of degree in. Thus we have where L.H.S C1 ^ C C1 : a? a C : 8 R [ P () = =) P () ]: We see immediatel C1 a 1 a < 1 ecause a 6= 1. Let us now eliminate the quantier from C. C 8 R [ a(a + )? a = =) P () ] 8 R [ (a + )? a = =) P () ] 8 R [ a +? a 1 = =) P () ] 8 R [ = =) P () ] P : Finall we otained a quantier-free formula. Let us simplif this formula. P a +? a?! a 1? a? a 1? a? a? a? a? a? 1? a? a? a a? 1? a? 1? a a +? 1 : In order to link the second and the third steps from the ottom, we have used the elementar fact that (u? v ) has the same sign as (u? v). Now putting together, we nall otain L.H.S. C1 ^ C a < 1 ^?1? a a < 1 ^ a +? 1 a +? 1 a +? 1 Remark: This univariate quantier elimination prolem is small enough to e carried out a computer program. In fact, we have tried the quantier elimination program qepcad []. Given the input 8 (a + )? a? ; the program, in aout 5 seconds on a DEC station, outputs the quantier-free formula: a 6 + a + a + 6? a + 1a? + a + 1: B suitale rewriting, one see that it is equivalent to a + 1: Remark: It is interesting to note that the polnomial shown aove a 6 + a + a + 6? a + 1a? + a +? 1 is the discriminant of the input polnomial P () = (a + )? a? with respect to upto some trivial factor. Thus, for this particular case, the discriminant computation was sucient. It is ecause P has onl one real root, sa. In that case, it is well known that P is the discriminant of P () upto some constant multiple. But it is not true in general when P has more than one real root. Lemma (Bivariate quantier elimination) Let Proof. P (; ) = (a + )? a? : 8(; ) R P (; ) a + L.H.S. [ 8 R P (; ) ] ^ [ 8 R8 R? fg P (; ) ] [ 8 R (a? a ) ] ^ [ 8 R8 R? fg a +? a?! ]
4 a 1 ^ 8z R (az + )? a z? a 1 ^ a + 1 ( Lemma 1) a + 1 Lemma (Necessar condition) 8(; ) R jmj 1 =) a + Proof. The proof will e case analsis. The rst three cases will take care of \degenerate" situations. Case: a = and = We see that r = 1; s = ; q =. Thus, jmj = 1. Hence we have L.H.S. 8(; ) R R.H.S Case: a = and = 1 We see that r = cos ; s = sin ; q =. Thus, jmj = 1. Hence we have L.H.S. 8(; ) R (1) 1 R.H.S Case: = and a = 1 We see that r = cos ; s = sin ; q =. Thus, jmj = 1. Hence we have L.H.S. 8(; ) R 1 1 (1) + 1 R.H.S Case: Otherwise B epanding jmj with respect to and around (; ), we otain Let jmj = 1 + 1? (a + )? a? + O( ): P (; ) = (a + )? a? : One can easil verif that P is a non-zero polnomial (i.e. does not vanish identicall). 1 Thus, the ehaviour of jmj near the origin is determined that of P. Hence we have L.H.S. =) 8(; ) R P (; ) a + 1 ( Lemma ) R.H.S 1 In fact, the rst three cases were introduced in order to lter out the situation when P is a zero polnomial. Lemma (Sucient condition) Proof. 8(; ) R jmj 1 (= a + Let us rst note jmj = jl? q j Thus we have = (r + is)? q j = jr? s? q + rsij = (r? s? q ) + r s = r + s + q? r s? r q + s q + r s = r + s + q + r s + r q + s q? r q = (r + s + q )? r q : 8(; ) R jmj 1 8(; ) R (r + s + q )? r q 1 (= 8(; ) R (r + s + q ) 1 8(; ) R r + s + q 1: (1) A digression to tell an inside stor. Aove, I took a \old" step of ignoring the term r q. It is hard to elieve that the proof will go through without the term, ut it turns out to e so. In fact, I have spent several months, in vain, searching for proofs which take account of the term ecause I was preconditioned to the idea that the term must pla an essential role. Onl after drawing lots of pictures using plotting programs and inspecting lots of computation traces using computer algera sstems, one morning, I \accidentl" hit upon the intuition that the term might not pla an important role for stailit. Then I could reach, after onl few more hours of work, the end of the proof, to m jo! The end of the digression. Since and alwas occur as arguments of trigonometric functions, we can restrict their ranges. 8(; ) R r + s + q 1 8(; ) [?; ] r + s + q 1: Here the choice of [?; ] instead of [; ] is intentional as will e seen now. Let?1 ; 1 cos ; cos 1 = sin ; = sin : cos = p 1?, cos = p 1?. (If we had chosen [; ] as the range of and instead, we could have had dicult since we could not have determined the sign of cos and cos from and.) We rewrite r; s; and q in terms of and, otaining r = 1? a?? a p 1? p 1? s = ap 1? + p 1? q = a :
5 After some straightforward calculation, we otain w 1 (; ) = r + s + q Now oserve = 1? a? +a + +8a +16a (1?? + ) +16a p 1? p 1? +16a p 1? p 1? : 8(; ) [?; ] r + s + q 1 8(; ) [?1; 1] w 1 (; ) 1: () 1?? + = 1? (1? )? (1? ) 1: Thus we have 8(; ) [?1; 1] w 1 (; ) 1 (= 8(; ) [?1; 1] w (; ) 1 where w (; ) is otained from w 1 (; ) replacing the term (1?? + ) with 1, that is, w (; ) = 1? a? +a + +a +16a p 1? p 1? +16a p 1? p 1? : Let w (; ) e otained from w (; ) dropping the term p 1? p 1?, that is, w (; ) = 1? a? Now we claim that +a + +a +16a +16a : 8(; ) [?1; 1] w (; ) 1 (= 8(; ) [?1; 1] w (; ) 1 (5) In order to prove this claim, let us assume that 8(; ) [?1; 1] w (; ) 1: Let ( ; ) [?1; 1] e aritrar ut ed. show that w ( ; ) 1: We will do so case analsis. Case: a Clearl, we have w ( ; ) w ( ; ): we need to From the hpothesis, we also have Thus, we have Case: a < Clearl, we have w ( ; ) 1: w ( ; ) 1: w ( ; ) w ( ;? ) w ( ;? ): From the hpothesis, we also have Thus, we have w ( ;? ) 1: w ( ; ) 1: This completes the proof of the claim that we can drop the radical terms. Now note that 8(; ) [?1; 1] w (; ) 1 8(; ) [?1; 1]?a? + (a + ) (= 8(; ) R? a? + (a + ) a + 1 ( Lemma ): (6) Finall, appling the transitivit of the logical implication ((=) to the implication relations in (1),, (),, (5) and (6), we conclude that 8(; ) R jmj 1 (= a + From Lemmas and, Theorem 1 immediatel follows. Concluding Remarks In this paper, we have proved that the stailit region of the MacCormack scheme for the two dimensional scalar avection equation is eactl given a + As Turkel [11] has done for the La-Wendro scheme, it seems to e possile to generalize the result of this paper to n- valued functions where a and are commuting real smmetric matrices. This will e investigated in a susequent note. Though it is not et computationall feasile to automate the whole process a computer program, some part of the proof could e alread carried out computer programs. Furthermore various conjecturing processes were greatl aided computer programs. These eperience will also e reported in another paper. Acknowledgment: I would like to thank Stanle Steinerg and Richard Liska ecause the introduced me to the stailit prolem of MacCormack scheme and also rought Wendro's work to m attention. 5
6 References [1] G. E. Collins. Quantier elimination for the elementar theor of real closed elds clindrical algeraic decomposition. In Lecture Notes In Computer Science, pages 1{18. Springer-Verlag, Berlin, Vol.. [] G. E. Collins and H. Hong. Partial clindrical algeraic decomposition for quantier elimination. Journal of Smolic Computation, 1():99{8, Septemer [] H. Hong. An improvement of the projection operator in clindrical algeraic decomposition. In International Smposium of Smolic and Algeraic Computation ISSAC-9, pages 61{6, 199. [] H. Hong. Improvements in CAD{ased Quantier Elimination. PhD thesis, The Ohio State Universit, 199. [5] H. Hong. Simple solution formula construction in clindrical algeraic decomposition ased quantier elimination. In International Conference on Smolic and Algeraic Computation ISSAC-9, pages 177{188, 199. [6] P.D. La and B. Wendro. Dierence schemes for hperolic equations with higher orer of accurac. Comm. Pure Appl. Math, XVII:81{98, 196. [7] R. Liska and S. Steinerg. Appling Quantier Elimination to Stailit Analsis of Dierence Schemes. The Computer Journal, 6(5):97{5, 199. [8] R. MacCormack. The eect of viscosit in hpervelocit impact cratering. AIAA Paper No., (69-5), [9] J. Renegar. On the computational compleit and geometr of the rst-order theor of the reals (part III). Technical Report 856, Cornell Universit, Ithaca, New York USA, August [1] A. Tarski. A Decision Method for Elementar Algera and Geometr. Univ. of California Press, Berkele, second edition, [11] E. Turkel. Smmetric hperolic dierence schemes and matri prolems. Lin. Alg. Appl., 16:19{19, [1] B. Wendro. The stailit of MacCormack's Method for the Scalar Avection Equations. App. Math. Lett., :89{91,
Generalized Least-Squares Regressions II: Theory and Classication
Recent Advances in Intelligent Control, Modelling Comutational Science Generalized Least-Squares Regressions II Theor Classication NATANIEL GREENE Deartment of Mathematics Comuter Science Kingsorough Communit
More informationFACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures AB = BA = I,
FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 7 MATRICES II Inverse of a matri Sstems of linear equations Solution of sets of linear equations elimination methods 4
More informationMAT 1275: Introduction to Mathematical Analysis. Graphs and Simplest Equations for Basic Trigonometric Functions. y=sin( x) Function
MAT 275: Introduction to Mathematical Analsis Dr. A. Rozenblum Graphs and Simplest Equations for Basic Trigonometric Functions We consider here three basic functions: sine, cosine and tangent. For them,
More informationab is shifted horizontally by h units. ab is shifted vertically by k units.
Algera II Notes Unit Eight: Eponential and Logarithmic Functions Sllaus Ojective: 8. The student will graph logarithmic and eponential functions including ase e. Eponential Function: a, 0, Graph of an
More informationGraphs and polynomials
5_6_56_MQVMM - _t Page Frida, Novemer 8, 5 :5 AM MQ Maths Methods / Final Pages / 8//5 Graphs and polnomials VCEcoverage Areas of stud Units & Functions and graphs Algera In this chapter A The inomial
More informationGraphs and polynomials
1 1A The inomial theorem 1B Polnomials 1C Division of polnomials 1D Linear graphs 1E Quadratic graphs 1F Cuic graphs 1G Quartic graphs Graphs and polnomials AreAS of STud Graphs of polnomial functions
More informationMathematics 309 Conic sections and their applicationsn. Chapter 2. Quadric figures. ai,j x i x j + b i x i + c =0. 1. Coordinate changes
Mathematics 309 Conic sections and their applicationsn Chapter 2. Quadric figures In this chapter want to outline quickl how to decide what figure associated in 2D and 3D to quadratic equations look like.
More informationTopological structures and phases. in U(1) gauge theory. Abstract. We show that topological properties of minimal Dirac sheets as well as of
BUHEP-94-35 Decemer 1994 Topological structures and phases in U(1) gauge theory Werner Kerler a, Claudio Rei and Andreas Weer a a Fachereich Physik, Universitat Marurg, D-35032 Marurg, Germany Department
More information10 Back to planar nonlinear systems
10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form
More informationRolle s Theorem. THEOREM 3 Rolle s Theorem. x x. then there is at least one number c in (a, b) at which ƒ scd = 0.
4.2 The Mean Value Theorem 255 4.2 The Mean Value Theorem f '(c) 0 f() We know that constant functions have zero derivatives, ut could there e a complicated function, with man terms, the derivatives of
More informationOn Range and Reflecting Functions About the Line y = mx
On Range and Reflecting Functions About the Line = m Scott J. Beslin Brian K. Heck Jerem J. Becnel Dept.of Mathematics and Dept. of Mathematics and Dept. of Mathematics and Computer Science Computer Science
More informationCopyrighted by Gabriel Tang B.Ed., B.Sc. Page 111.
Algera Chapter : Polnomial and Rational Functions Chapter : Polnomial and Rational Functions - Polnomial Functions and Their Graphs Polnomial Functions: - a function that consists of a polnomial epression
More informationChapter #6 EEE State Space Analysis and Controller Design
hapter #6 EEE3-83 State Space nalsis and ontroller Design State feedack Linear Quadratic Regulator (LQR) Open and closed loop estimators Reduced Order Estimators (not assessed) racking Module Leader: Dr
More informationwe make slices perpendicular to the x-axis. If the slices are thin enough, they resemble x cylinders or discs. The formula for the x
Math Learning Centre Solids of Revolution When we rotate a curve around a defined ais, the -D shape created is called a solid of revolution. In the same wa that we can find the area under a curve calculating
More information8.1 Exponents and Roots
Section 8. Eponents and Roots 75 8. Eponents and Roots Before defining the net famil of functions, the eponential functions, we will need to discuss eponent notation in detail. As we shall see, eponents
More informationINTRODUCTION TO DIFFERENTIAL EQUATIONS
INTRODUCTION TO DIFFERENTIAL EQUATIONS. Definitions and Terminolog. Initial-Value Problems.3 Differential Equations as Mathematical Models CHAPTER IN REVIEW The words differential and equations certainl
More informationCoordinate geometry. + bx + c. Vertical asymptote. Sketch graphs of hyperbolas (including asymptotic behaviour) from the general
A Sketch graphs of = a m b n c where m = or and n = or B Reciprocal graphs C Graphs of circles and ellipses D Graphs of hperbolas E Partial fractions F Sketch graphs using partial fractions Coordinate
More informationPolynomial Degree and Finite Differences
CONDENSED LESSON 7.1 Polynomial Degree and Finite Differences In this lesson, you Learn the terminology associated with polynomials Use the finite differences method to determine the degree of a polynomial
More information5Higher-degree ONLINE PAGE PROOFS. polynomials
5Higher-degree polnomials 5. Kick off with CAS 5.2 Quartic polnomials 5.3 Families of polnomials 5.4 Numerical approimations to roots of polnomial equations 5.5 Review 5. Kick off with CAS Quartic transformations
More informationSection 3.1. ; X = (0, 1]. (i) f : R R R, f (x, y) = x y
Paul J. Bruillard MATH 0.970 Problem Set 6 An Introduction to Abstract Mathematics R. Bond and W. Keane Section 3.1: 3b,c,e,i, 4bd, 6, 9, 15, 16, 18c,e, 19a, 0, 1b Section 3.: 1f,i, e, 6, 1e,f,h, 13e,
More informationLESSON #11 - FORMS OF A LINE COMMON CORE ALGEBRA II
LESSON # - FORMS OF A LINE COMMON CORE ALGEBRA II Linear functions come in a variet of forms. The two shown below have been introduced in Common Core Algebra I and Common Core Geometr. TWO COMMON FORMS
More informationChapter 4 Analytic Trigonometry
Analtic Trigonometr Chapter Analtic Trigonometr Inverse Trigonometric Functions The trigonometric functions act as an operator on the variable (angle, resulting in an output value Suppose this process
More informationUNCORRECTED SAMPLE PAGES. 3Quadratics. Chapter 3. Objectives
Chapter 3 3Quadratics Objectives To recognise and sketch the graphs of quadratic polnomials. To find the ke features of the graph of a quadratic polnomial: ais intercepts, turning point and ais of smmetr.
More informationMath 030 Review for Final Exam Revised Fall 2010 RH/ DM 1
Math 00 Review for Final Eam Revised Fall 010 RH/ DM 1 1. Solve the equations: (-1) (7) (-) (-1) () 1 1 1 1 f. 1 g. h. 1 11 i. 9. Solve the following equations for the given variable: 1 Solve for. D ab
More informationMining the Quadratic Formula: An Integrated Approach to Quadratic Function
數學教育第二十期 (6/005) Mining the Quadratic Formula: An Integrated Approach to Quadratic Function Lee Kok Ming Li Po Chun United World College of Hong Kong 1. Introduction Mathematic lessons should not e presented
More informationThe Theory of a Scroll Profile
Purdue Universit Purdue e-pus International Compressor Engineering Conference School of Mechanical Engineering 006 The Theor of a Scroll Profile Huert Bukac Little Dnamics Follow this and additional works
More informationLESSON #42 - INVERSES OF FUNCTIONS AND FUNCTION NOTATION PART 2 COMMON CORE ALGEBRA II
LESSON #4 - INVERSES OF FUNCTIONS AND FUNCTION NOTATION PART COMMON CORE ALGEBRA II You will recall from unit 1 that in order to find the inverse of a function, ou must switch and and solve for. Also,
More informationSEPARABLE EQUATIONS 2.2
46 CHAPTER FIRST-ORDER DIFFERENTIAL EQUATIONS 4. Chemical Reactions When certain kinds of chemicals are combined, the rate at which the new compound is formed is modeled b the autonomous differential equation
More informationWeek 3 September 5-7.
MA322 Weekl topics and quiz preparations Week 3 September 5-7. Topics These are alread partl covered in lectures. We collect the details for convenience.. Solutions of homogeneous equations AX =. 2. Using
More information2.2 SEPARABLE VARIABLES
44 CHAPTER FIRST-ORDER DIFFERENTIAL EQUATIONS 6 Consider the autonomous DE 6 Use our ideas from Problem 5 to find intervals on the -ais for which solution curves are concave up and intervals for which
More informationExpansion formula using properties of dot product (analogous to FOIL in algebra): u v 2 u v u v u u 2u v v v u 2 2u v v 2
Least squares: Mathematical theory Below we provide the "vector space" formulation, and solution, of the least squares prolem. While not strictly necessary until we ring in the machinery of matrix algera,
More informationES.1803 Topic 16 Notes Jeremy Orloff
ES803 Topic 6 Notes Jerem Orloff 6 Eigenalues, diagonalization, decoupling This note coers topics that will take us seeral classes to get through We will look almost eclusiel at 2 2 matrices These hae
More informationHigher. Polynomials and Quadratics. Polynomials and Quadratics 1
Higher Mathematics Polnomials and Quadratics Contents Polnomials and Quadratics 1 1 Quadratics EF 1 The Discriminant EF Completing the Square EF Sketching Paraolas EF 7 5 Determining the Equation of a
More informationTENSOR TRANSFORMATION OF STRESSES
GG303 Lecture 18 9/4/01 1 TENSOR TRANSFORMATION OF STRESSES Transformation of stresses between planes of arbitrar orientation In the 2-D eample of lecture 16, the normal and shear stresses (tractions)
More informationTHE KINETICS OF PHASE TRANSITIONS: I. THE ROLE OF HETEROGENEITY. F. Paladi* and V. Gamurari
THE KINETICS OF PHASE TRANSITIONS: I. THE ROLE OF HETEROGENEITY F. Paladi* and V. Gamurari Department of Theoretical Phsics Moldova State Universit 60 A. Mateevici str. MD-009 Chisinau Repulic of Moldova
More informationCubic and quartic functions
3 Cubic and quartic functions 3A Epanding 3B Long division of polnomials 3C Polnomial values 3D The remainder and factor theorems 3E Factorising polnomials 3F Sum and difference of two cubes 3G Solving
More informationLESSON #12 - FORMS OF A LINE COMMON CORE ALGEBRA II
LESSON # - FORMS OF A LINE COMMON CORE ALGEBRA II Linear functions come in a variet of forms. The two shown below have been introduced in Common Core Algebra I and Common Core Geometr. TWO COMMON FORMS
More informationThe Maze Generation Problem is NP-complete
The Mae Generation Problem is NP-complete Mario Alviano Department of Mathematics, Universit of Calabria, 87030 Rende (CS), Ital alviano@mat.unical.it Abstract. The Mae Generation problem has been presented
More informationOrdinary Differential Equations of First Order
CHAPTER 1 Ordinar Differential Equations of First Order 1.1 INTRODUCTION Differential equations pla an indispensable role in science technolog because man phsical laws relations can be described mathematicall
More informationON THE INTERPRETATION OF THE LAGRANGE MULTIPLIERS IN THE CONSTRAINT FORMULATION OF CONTACT PROBLEMS; OR WHY ARE SOME MULTIPLIERS ALWAYS ZERO?
Proceedings of the ASME 214 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference IDETC/CIE 214 August 17-2, 214, Buffalo, New York, USA DETC214-3479
More informationTrigonometric Functions
TrigonometricReview.nb Trigonometric Functions The trigonometric (or trig) functions are ver important in our stud of calculus because the are periodic (meaning these functions repeat their values in a
More information5. Zeros. We deduce that the graph crosses the x-axis at the points x = 0, 1, 2 and 4, and nowhere else. And that s exactly what we see in the graph.
. Zeros Eample 1. At the right we have drawn the graph of the polnomial = ( 1) ( 2) ( 4). Argue that the form of the algebraic formula allows ou to see right awa where the graph is above the -ais, where
More information16.5. Maclaurin and Taylor Series. Introduction. Prerequisites. Learning Outcomes
Maclaurin and Talor Series 6.5 Introduction In this Section we eamine how functions ma be epressed in terms of power series. This is an etremel useful wa of epressing a function since (as we shall see)
More information0.24 adults 2. (c) Prove that, regardless of the possible values of and, the covariance between X and Y is equal to zero. Show all work.
1 A socioeconomic stud analzes two discrete random variables in a certain population of households = number of adult residents and = number of child residents It is found that their joint probabilit mass
More informationMathematics. Mathematics 2. hsn.uk.net. Higher HSN22000
Higher Mathematics UNIT Mathematics HSN000 This document was produced speciall for the HSN.uk.net website, and we require that an copies or derivative works attribute the work to Higher Still Notes. For
More information67. (a) Use a computer algebra system to find the partial fraction CAS. 68. (a) Find the partial fraction decomposition of the function CAS
SECTION 7.5 STRATEGY FOR INTEGRATION 483 6. 2 sin 2 2 cos CAS 67. (a) Use a computer algebra sstem to find the partial fraction decomposition of the function 62 63 Find the area of the region under the
More information5.4 dividing POlynOmIAlS
SECTION 5.4 dividing PolNomiAls 3 9 3 learning ObjeCTIveS In this section, ou will: Use long division to divide polnomials. Use snthetic division to divide polnomials. 5.4 dividing POlnOmIAlS Figure 1
More informationragsdale (zdr82) HW7 ditmire (58335) 1 The magnetic force is
ragsdale (zdr8) HW7 ditmire (585) This print-out should have 8 questions. Multiple-choice questions ma continue on the net column or page find all choices efore answering. 00 0.0 points A wire carring
More informationMathematics. Mathematics 2. hsn.uk.net. Higher HSN22000
hsn.uk.net Higher Mathematics UNIT Mathematics HSN000 This document was produced speciall for the HSN.uk.net website, and we require that an copies or derivative works attribute the work to Higher Still
More information#A50 INTEGERS 14 (2014) ON RATS SEQUENCES IN GENERAL BASES
#A50 INTEGERS 14 (014) ON RATS SEQUENCES IN GENERAL BASES Johann Thiel Dept. of Mathematics, New York City College of Technology, Brooklyn, New York jthiel@citytech.cuny.edu Received: 6/11/13, Revised:
More informationIdentifying second degree equations
Chapter 7 Identifing second degree equations 71 The eigenvalue method In this section we appl eigenvalue methods to determine the geometrical nature of the second degree equation a 2 + 2h + b 2 + 2g +
More informationLab 5 Forces Part 1. Physics 211 Lab. You will be using Newton s 2 nd Law to help you examine the nature of these forces.
b Lab 5 Forces Part 1 Phsics 211 Lab Introduction This is the first week of a two part lab that deals with forces and related concepts. A force is a push or a pull on an object that can be caused b a variet
More informationLab 5 Forces Part 1. Physics 225 Lab. You will be using Newton s 2 nd Law to help you examine the nature of these forces.
b Lab 5 orces Part 1 Introduction his is the first week of a two part lab that deals with forces and related concepts. A force is a push or a pull on an object that can be caused b a variet of reasons.
More informationGet Solution of These Packages & Learn by Video Tutorials on Matrices
FEE Download Stud Package from website: wwwtekoclassescom & wwwmathsbsuhagcom Get Solution of These Packages & Learn b Video Tutorials on wwwmathsbsuhagcom Matrices An rectangular arrangement of numbers
More informationPhysics Gravitational force. 2. Strong or color force. 3. Electroweak force
Phsics 360 Notes on Griffths - pluses and minuses No tetbook is perfect, and Griffithsisnoeception. Themajorplusisthat it is prett readable. For minuses, see below. Much of what G sas about the del operator
More informationPure Further Mathematics 2. Revision Notes
Pure Further Mathematics Revision Notes October 016 FP OCT 016 SDB Further Pure 1 Inequalities... 3 Algebraic solutions... 3 Graphical solutions... 4 Series Method of Differences... 5 3 Comple Numbers...
More informationPrecalculus Prerequisite Packet Paint Branch High School Math Department. Concepts To Be Assessed on the Precalculus Course Pre-assessment.
Updated /01 The problems in this packet are designed to help ou review topics from previous math courses that are important to our success in Precalculus. It is important that ou take time during summer
More informationf 0 ab a b: base f
Precalculus Notes: Unit Eponential and Logarithmic Functions Sllaus Ojective: 9. The student will sketch the graph of a eponential, logistic, or logarithmic function. 9. The student will evaluate eponential
More informationGeneralized Least-Squares Regressions III: Further Theory and Classication
Generalized Least-Squares Regressions III: Further Theor Classication NATANIEL GREENE Department of Mathematics Computer Science Kingsborough Communit College, CUNY Oriental Boulevard, Brookln, NY UNITED
More informationLecture 5. Equations of Lines and Planes. Dan Nichols MATH 233, Spring 2018 University of Massachusetts.
Lecture 5 Equations of Lines and Planes Dan Nichols nichols@math.umass.edu MATH 233, Spring 2018 Universit of Massachusetts Februar 6, 2018 (2) Upcoming midterm eam First midterm: Wednesda Feb. 21, 7:00-9:00
More informationIn this chapter a student has to learn the Concept of adjoint of a matrix. Inverse of a matrix. Rank of a matrix and methods finding these.
MATRICES UNIT STRUCTURE.0 Objectives. Introduction. Definitions. Illustrative eamples.4 Rank of matri.5 Canonical form or Normal form.6 Normal form PAQ.7 Let Us Sum Up.8 Unit End Eercise.0 OBJECTIVES In
More information39. (a) Use trigonometric substitution to verify that. 40. The parabola y 2x divides the disk into two
35. Prove the formula A r for the area of a sector of a circle with radius r and central angle. [Hint: Assume 0 and place the center of the circle at the origin so it has the equation. Then is the sum
More information10.2 The Unit Circle: Cosine and Sine
0. The Unit Circle: Cosine and Sine 77 0. The Unit Circle: Cosine and Sine In Section 0.., we introduced circular motion and derived a formula which describes the linear velocit of an object moving on
More information11.4 Polar Coordinates
11. Polar Coordinates 917 11. Polar Coordinates In Section 1.1, we introduced the Cartesian coordinates of a point in the plane as a means of assigning ordered pairs of numbers to points in the plane.
More informationRising HONORS Algebra 2 TRIG student Summer Packet for 2016 (school year )
Rising HONORS Algebra TRIG student Summer Packet for 016 (school ear 016-17) Welcome to Algebra TRIG! To be successful in Algebra Trig, ou must be proficient at solving and simplifing each tpe of problem
More informationApplications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
More informationAnalysis of Effects of Rebounds and Aerodynamics for Trajectory of Table Tennis Ball
SICE Annual Conference 21 August 18-21, 21, The Grand Hotel, Taipei, Taiwan Analsis of Effects of Reounds and Aerodnamics for Trajector of Tale Tennis Junko Nonomura Mechanical Science and Engineering,
More informationDifferentiation and applications
FS O PA G E PR O U N C O R R EC TE D Differentiation and applications. Kick off with CAS. Limits, continuit and differentiabilit. Derivatives of power functions.4 C oordinate geometr applications of differentiation.5
More informationParticular Solutions
Particular Solutions Our eamples so far in this section have involved some constant of integration, K. We now move on to see particular solutions, where we know some boundar conditions and we substitute
More information1.3 LIMITS AT INFINITY; END BEHAVIOR OF A FUNCTION
. Limits at Infinit; End Behavior of a Function 89. LIMITS AT INFINITY; END BEHAVIOR OF A FUNCTION Up to now we have been concerned with its that describe the behavior of a function f) as approaches some
More information8 Differential Calculus 1 Introduction
8 Differential Calculus Introduction The ideas that are the basis for calculus have been with us for a ver long time. Between 5 BC and 5 BC, Greek mathematicians were working on problems that would find
More informationToda s Theorem: PH P #P
CS254: Computational Compleit Theor Prof. Luca Trevisan Final Project Ananth Raghunathan 1 Introduction Toda s Theorem: PH P #P The class NP captures the difficult of finding certificates. However, in
More information8. BOOLEAN ALGEBRAS x x
8. BOOLEAN ALGEBRAS 8.1. Definition of a Boolean Algebra There are man sstems of interest to computing scientists that have a common underling structure. It makes sense to describe such a mathematical
More information10.5 Graphs of the Trigonometric Functions
790 Foundations of Trigonometr 0.5 Graphs of the Trigonometric Functions In this section, we return to our discussion of the circular (trigonometric functions as functions of real numbers and pick up where
More informationA NEW CHAOTIC SYSTEM AND BEYOND: THE GENERALIZED LORENZ-LIKE SYSTEM
Tutorials and Reviews International Journal of Bifurcation and Chaos, Vol. 14, No. 4) 17 17 c World Scientific Pulishing Compan A NEW CHAOTIC SYSTEM AND BEYOND: THE GENERALIZED LORENZ-LIKE SYSTEM JINHU
More informationMath 123 Summary of Important Algebra & Trigonometry Concepts Chapter 1 & Appendix D, Stewart, Calculus Early Transcendentals
Math Summar of Important Algebra & Trigonometr Concepts Chapter & Appendi D, Stewart, Calculus Earl Transcendentals Function a rule that assigns to each element in a set D eactl one element, called f (
More information1 HOMOGENEOUS TRANSFORMATIONS
HOMOGENEOUS TRANSFORMATIONS Purpose: The purpose of this chapter is to introduce ou to the Homogeneous Transformation. This simple 4 4 transformation is used in the geometr engines of CAD sstems and in
More informationDiagnostic Tests Study Guide
California State Universit, Sacramento Department of Mathematics and Statistics Diagnostic Tests Stud Guide Descriptions Stud Guides Sample Tests & Answers Table of Contents: Introduction Elementar Algebra
More informationa [A] +Algebra 2/Trig Final Exam Review Fall Semester x [E] None of these [C] 512 [A] [B] 1) Simplify: [D] x z [E] None of these 2) Simplify: [A]
) Simplif: z z z 6 6 z 6 z 6 ) Simplif: 9 9 0 ) Simplif: a a a 0 a a ) Simplif: 0 0 ) Simplif: 9 9 6) Evaluate: / 6 6 6 ) Rationalize: ) Rationalize: 6 6 0 6 9) Which of the following are polnomials? None
More informationMath Notes on sections 7.8,9.1, and 9.3. Derivation of a solution in the repeated roots case: 3 4 A = 1 1. x =e t : + e t w 2.
Math 7 Notes on sections 7.8,9., and 9.3. Derivation of a solution in the repeated roots case We consider the eample = A where 3 4 A = The onl eigenvalue is = ; and there is onl one linearl independent
More information* τσ σκ. Supporting Text. A. Stability Analysis of System 2
Supporting Tet A. Stabilit Analsis of Sstem In this Appendi, we stud the stabilit of the equilibria of sstem. If we redefine the sstem as, T when -, then there are at most three equilibria: E,, E κ -,,
More informationQUADRATIC FUNCTION REVIEW
Name: Date: QUADRATIC FUNCTION REVIEW Linear and eponential functions are used throughout mathematics and science due to their simplicit and applicabilit. Quadratic functions comprise another ver important
More informationLESSON #28 - POWER FUNCTIONS COMMON CORE ALGEBRA II
1 LESSON #8 - POWER FUNCTIONS COMMON CORE ALGEBRA II Before we start to analze polnomials of degree higher than two (quadratics), we first will look at ver simple functions known as power functions. The
More information1.6 ELECTRONIC STRUCTURE OF THE HYDROGEN ATOM
1.6 ELECTRONIC STRUCTURE OF THE HYDROGEN ATOM 23 How does this wave-particle dualit require us to alter our thinking about the electron? In our everda lives, we re accustomed to a deterministic world.
More informationPolynomial and Rational Functions
Polnomial and Rational Functions 5 Figure 1 35-mm film, once the standard for capturing photographic images, has been made largel obsolete b digital photograph. (credit film : modification of work b Horia
More informationNew approach to study the van der Pol equation for large damping
Electronic Journal of Qualitative Theor of Differential Equations 2018, No. 8, 1 10; https://doi.org/10.1422/ejqtde.2018.1.8 www.math.u-szeged.hu/ejqtde/ New approach to stud the van der Pol equation for
More information6 = 1 2. The right endpoints of the subintervals are then 2 5, 3, 7 2, 4, 2 9, 5, while the left endpoints are 2, 5 2, 3, 7 2, 4, 9 2.
5 THE ITEGRAL 5. Approimating and Computing Area Preliminar Questions. What are the right and left endpoints if [, 5] is divided into si subintervals? If the interval [, 5] is divided into si subintervals,
More informationManeuvering Target Tracking Method based on Unknown but Bounded Uncertainties
Preprints of the 8th IFAC World Congress Milano (Ital) ust 8 - Septemer, Maneuvering arget racing Method ased on Unnown ut Bounded Uncertainties Hodjat Rahmati*, Hamid Khaloozadeh**, Moosa Aati*** * Facult
More informationOn Information and Sufficiency
On Information and Sufficienc Huaiu hu SFI WORKING PAPER: 997-02-04 SFI Working Papers contain accounts of scientific work of the author(s) and do not necessaril represent the views of the Santa Fe Institute.
More informationExact Free Vibration of Webs Moving Axially at High Speed
Eact Free Viration of Wes Moving Aially at High Speed S. HATAMI *, M. AZHARI, MM. SAADATPOUR, P. MEMARZADEH *Department of Engineering, Yasouj University, Yasouj Department of Civil Engineering, Isfahan
More informationMATRIX TRANSFORMATIONS
CHAPTER 5. MATRIX TRANSFORMATIONS INSTITIÚID TEICNEOLAÍOCHTA CHEATHARLACH INSTITUTE OF TECHNOLOGY CARLOW MATRIX TRANSFORMATIONS Matri Transformations Definition Let A and B be sets. A function f : A B
More informationFinding Complex Solutions of Quadratic Equations
y - y - - - x x Locker LESSON.3 Finding Complex Solutions of Quadratic Equations Texas Math Standards The student is expected to: A..F Solve quadratic and square root equations. Mathematical Processes
More informationThe Force Table Introduction: Theory:
1 The Force Table Introduction: "The Force Table" is a simple tool for demonstrating Newton s First Law and the vector nature of forces. This tool is based on the principle of equilibrium. An object is
More informationFraction-free Row Reduction of Matrices of Skew Polynomials
Fraction-free Row Reduction of Matrices of Skew Polynomials Bernhard Beckermann Laboratoire d Analyse Numérique et d Optimisation Université des Sciences et Technologies de Lille France bbecker@ano.univ-lille1.fr
More informationName Please print your name as it appears on the class roster.
Berkele Cit College Practice Problems Math 1 Precalculus - Final Eam Preparation Name Please print our name as it appears on the class roster. SHORT ANSWER. Write the word or phrase that best completes
More informationAn analytic proof of the theorems of Pappus and Desargues
Note di Matematica 22, n. 1, 2003, 99 106. An analtic proof of the theorems of Pappus and Desargues Erwin Kleinfeld and Tuong Ton-That Department of Mathematics, The Universit of Iowa, Iowa Cit, IA 52242,
More informationNon-regular languages
Non-regular languages (Pumping Lemma) Fall 2017 Costas Busch - RPI 1 Non-regular languages { a n n : n 0} { vv R : v{ a, }*} Regular languages a* * c a c( a )* etc... Fall 2017 Costas Busch - RPI 2 How
More informationComputer Graphics: 2D Transformations. Course Website:
Computer Graphics: D Transformations Course Website: http://www.comp.dit.ie/bmacnamee 5 Contents Wh transformations Transformations Translation Scaling Rotation Homogeneous coordinates Matri multiplications
More informationStrong echo eect and nonlinear transient growth in. shear ows U.S.A.
Strong echo eect and nonlinear transient growth in shear ows J. Vanneste 1,P. J. Morrison and T. Warn 3 1 Dept. of Phsics, Universit oftoronto, Toronto, Canada. Dept. of Phsics and Institute for Fusion
More informationParametric Equations for Circles and Ellipses
Lesson 5-8 Parametric Equations for Circles and Ellipses BIG IDEA Parametric equations use separate functions to defi ne coordinates and and to produce graphs Vocabular parameter parametric equations equation
More information