AMS 212B Perturbation Methods Lecture 05 Copyright by Hongyun Wang, UCSC
|
|
- Nelson Allen
- 5 years ago
- Views:
Transcription
1 AMS B Pertrbation Methods Lectre 5 Copright b Hongn Wang, UCSC Recap: we discssed bondar laer of ODE Oter epansion Inner epansion Matching: ) Prandtl s matching ) Matching b an intermediate variable (Skip it) Consider the range = η ϵ α where η > and / < α < are fied as +. Epand ( ot )( ) =η α Epand ( inn ) =η α Match (ot) and (inn) Note: This matching method does not prodce a composite epansion. We appl matching b an intermediate variable to or eample + + = =, ( ) = ( ot) inn = e + ( ), + ( ) = c ( e )+ c ( e ) c ( + e ) Epand ( ot )( ) and ( inn) =η α in and keep terms p to O(). =η α ot ( ) = =η α e e ηα + η α +! + = e η α +! - -
2 AMS B Pertrbation Methods innt = e( η α + )+! = c ( e ηα )+ c ( e ηα ) c ( + e ηα )η α +! =η α = c + c c η α +! = c c η α + c +! Matching condition: ( ot ) = ( inn ) =η α =η α ==> e( η α + )+!= c c η α + c +! ==> c = e, c = e Again, this matching method does not prodce a composite epansion. 3) Van Dke s matching We etend Prandtl s matching b considering more sophisticated limits, respectivel for small and for large. Inner limit of the oter epansion (ot)(inn) (): To calclate (ot)(inn) (), we do a) epand ( ot )( ) in the Talor series arond =, and then b) sbstitte = to write it in terms of. Oter limit of the inner epansion (inn)(ot) (): To calclate (inn)(ot) (), we treat e - as zero and keep other terms intact. The accrac of this epression improves eponentiall as increases. Note: Both (ot)(inn) () and (inn)(ot) () are fnctions, not scalars. Van Dke s matching condition: (ot)(inn) () = (inn)(ot) () Composite epansion: To constrct the composite epansion, we define - -
3 AMS B Pertrbation Methods ( ) ( inn )( ot) m = ot ( inn) Note: (m) () is a fnction, not a scalar. Using (m) () we write ot the composite epansion: ( ) = ( ot )( )+ ( inn )( ) ( m ) c ( ) = Remark: ) In the composite epansion, it is more convenient if we first calclate ( inn )( ) ( m )( ) and then evalate it at =. ) The composite epansion is one epression valid for the whole interval [, ]. We appl Van Dke s matching to or eample. + + = =, ( ) = ( ot) inn = e + ( ), + ( ) = c ( e )+ c ( e ) c ( + e ) The oter limit of inner epansion: ( inn )( ot) ( ) = c + c c The inner limit of oter epansion: We do Talor epansion arond =, and keep terms p to O(). Note that = ϵ = O()). ( ot )( ) = e e + e ( ) = e + +O ==> ( ot )( inn) Matching condition: ( inn )( ot) ( ) = ( ot )( inn) ==> c + c c ( ) = e + = e
4 AMS B Pertrbation Methods ==> c = e, c = e ==> ( m ) = e + e e Using c = e and c = e, we write ot (inn) () ( inn )( ) = e( e )+ e( e ) e( + e ) We calclate (inn) () - (m) () in preparation for writing ot the composite epansion ( inn) ( ) ( m) The composite epansion is c = e e + e = e + + ( ) = ( ) = ( ot )( )+ ( inn )( ) ( m ) = e + e + + The figre below compares the eact soltion and composite epansions for = =. Eact soltion asmptotic with term asmptotic with terms Qestion: Where is the bondar laer? We first consider the simple case of constant coefficients with a. We hope to find some gidelines from this simple case
5 AMS B Pertrbation Methods +a + b = = α, ( ) = β Characteristic eqation: λ +aλ + b =, a, At =, it becomes a linear eqation, aλ + b =. ==> λ ( ) = b a For, the characteristic eqation has two roots: λ λ ( ) b a ( )λ = b ==> λ ( ) a A general soltion of the ODE is = c e λ + c e λ (the reglar root) (the singlar root) The bondar laer is cased b the term e λ. To determine the location of bondar laer, we std Behavior of e λ when λ is large: For λ, inepect the behavior of mae λ eλ = e λ. << ==> For λ, the bondar laer is at =. For λ +, inepect the behavior of Conclsion:. mae λ eλ = e λ << ==> For λ +, the bondar laer is at =. When λ ( ) a, the bondar laer is at =
6 AMS B Pertrbation Methods When λ ( ) a +, the bondar laer is at =. Eample: + + = As +, λ ( ) a. ==> the bondar laer is at =. As -, λ ( ) a +, ==> the bondar laer is at =. Remark: The sign of a is determined b and coefficient a. It is not affected b fnction term (b ) in ODE ϵ + a + b =. We can tr to appl the reslt to an ODE with non-linear fnction term. Eample (nonlinear ODE with bondar laer at = ) + + = =, ( ) =, As -, a = + Highlights of the soltion: ==> the bondar laer is at = ) The ODE is non-linear. There is no eact soltion! ) Since the bondar laer is at =, the inner variable is selected as = 3) After scaling, ODE ( )+ ( )+ f ( ) = becomes - 6 -
7 ( ) ( )+ f AMS B Pertrbation Methods ( ) = This confirms that the fnction term ƒ(()) in ODE does not affect the leading term behavior of bondar laer (inclding the location and width). 4) For (, ), as ϵ - we have (instead of + ). In Van Dke matching, when calclating ( inn )( ot) ( ), we treat e as zero. 5) The composite epansion is ( ) = c ( +) e + e + log ( +)+ log ( + ) (See Appendi A for the derivation) e 8 e e ( ) Qestion: What happens if a =? Answer: Two possibilities: *) bondar laers; *) no soltion. Eample (Two bondar laers) = sin ( ) =, ( ) = Claim: It has two bondar laers., + We can arrive at this conclsion in two was. View #: (a simple view) At =, the eqation becomes = sin ==> it does not need or accommodate an bondar condition. Ths, we need two bondar laers to accommodate the two bondar conditions. View #: (a more detailed view) - 7 -
8 AMS B Pertrbation Methods Characteristic eqation of the ODE is λ = ==> λ =, λ = Recall the behavior of e λ. λ = + ==> a bondar laer is at =. λ = ==> a bondar laer is at =. Therefore, we have two bondar laers. The epression of λ and λ also sggests that both bondar laers have the width ϵ. Now we derive the epansions from the differential eqation assming that both bondar laers have the width ϵ. Later on, we will discss how to find the width of a bondar laer. Oter epansion: We seek an epansion of the form ( ot )( ) = a ( )+O( ) Note: no bondar condition is imposed on (ot). Sbstitting into eqation ields ( ot) = sin Inner epansion at = : Let = be the inner variable. We have d d = d d, d d = d d =, sin sin - 8 -
9 AMS B Pertrbation Methods In terms of inner variable, the eqation becomes ( ) = sin = sin cos +O( ) We seek an epansion of the form ( innl )( ) = a ( )+ a ( )+O Bondar condition at = : () = =, = ϵ. ==> a ( )+ a ( )+!= ==> a() =, a() = Note: Onl the bondar condition at = is imposed on the inner epansion at =. Sbstitting into eqation ields ( a ( )+ a ( ) ) a + a ( ) = sin cos +O( ) ==> a ( ) a ( ) sin + Both coefficients mst be zero. a a ( )+cos +O = : a a = sin a ( ) = ==> a = sin ( e )+ c e e (See Appendi B for the derivation) / : a a = cos a ( ) = ==> a = cos + c e e (See Appendi C for the derivation) The inner epansion at = is - 9 -
10 AMS B Pertrbation Methods ( ) = sin ( e )+ c ( e e )+ cos + c ( e e ) innl Matching at = : Recall the oter epansion: ( ) = sin ot The inner limit of oter epansion: We do Talor epansion arond =, and keep terms p to O( ϵ).. Note that = = O ( ot )( innl) ( ) = sin = sin +cos +! = sin + cos +! The oter limit of inner epansion: As +, = + and we treat e- as zero. ( innl )( ot) ( ) = sin + cos c e c e +! Van Dke matching condition: ( innl )( ot) ( ) = ( ot )( innl) ( ) ==> c =, c = ==> ( ml ) ( innl) = sin = sin ==> ( innl) ( ) ( ml) + cos ( e )+ cos = sin e - -
11 AMS B Pertrbation Methods Inner epansion at = : Note: We select the inner variable sch that ) v = at = and ) v + as moves ot of the bondar laer and into the oter epansion region. Let v = ( v) v be the inner variable. The ODE becomes = sin v + The inner epansion at = is innr = sin +cos v +O( ) ( v) = sin ( e v )+ d ( e v e v )+ cos v + d e v e v Matching at = : Recall the oter epansion: ( ) = sin ot The inner limit of oter epansion: We do Talor epansion arond =. ( ot )( innr) ( v) = sin = sin ( ) = sin cos ( )+! = sin cos v +! ( innr) ( v) ( mr) v = sin e v The composite epansion (inclding two bondar laers) is ( ) = ( ot )( )+ ( innl )( ) ( ml ) c ( ) = ( ) v= + ( innr )( v) ( mr ) v - -
12 AMS B Pertrbation Methods ( ) = sin +sin e sin e c Remark: In the composite epansion above, the ϵ term disappears. This is attribted to two factors: The roots of the characteristic eqation have the form of λ = c + c +! with c =. There is no ϵ term in the oter epansion If we change the problem to + = =, ( ) =, + there will be ϵ term in the composite epansion. The figre below compares the eact soltion and the composite epansion for = = Eact soltion asmptotic with term Eample (Asmptotic epansion does not eist) - -
13 AMS B Pertrbation Methods = sin ( ) =, ( ) =, We epand both () and sin in Forier sine series. (The se of sine series is dictated b the bondar conditions () = () =.) sin = k= a sin kπ k = α k sin( kπ) k= Sbstitting into the ODE ields + α k kπ sin kπ = a k sin kπ k= a ==> α k = k kπ + k= When =, we have ϵ(k π) + = and the BVP has no soltion. ( kπ) As ϵ -, it will enconter -/(k π) infinitel man times. Therefore, this problem does not have an asmptotic epansion. Appendi A Derivation of a two-term epansion for the bondar vale problem + + = =, ( ) =, As -, we have a = +. Oter epansion: ==> The bondar laer is at =
14 AMS B Pertrbation Methods We seek an epansion of the form ( ot )( ) = a ( )+ a ( )+! Bondar condition at = : () = = ==> a ( )+ a ( )+!= ==> a() =, a() = Note: When the bondar laer is at =, onl the bondar condition at = is imposed on the oter epansion. Sbstitting into eqation ( a +! )+( a + a +! )+ a + a +! = ==> a + a + a +a a + a +!= : a + a = a = ==> a a = ==> = a ==> = + c ==> a a ( ) = + c Using the bondar condition a() =, we obtain a ( ) = + : a +a a = a = a ( ) = ( +) 3 ==> a + + a = + 3 ==> + a +( +)a = ( +) - 4 -
15 AMS B Pertrbation Methods a ==> + ==> a = ( +) = + log + + ==> ( +) a = log( +)+ c c ( +) Using the bondar condition a() =, we obtain a ( ) = log + ( +) The oter epansion is ( ) = ot ( +) log + + Inner epansion: Let = be the inner variable. We have d d = d d, d d = d d Sbstitting into eqation, we have d d d d + = ==> ( ) ( )+ ( ) = This is the ODE in terms of the inner variable after scaling. We seek an epansion of the form ( inn )( ) = a ( )+ a ( )+! Bondar condition at = : () = =, ( ) = ϵ. ==> a ( )+ a ( )+!= ==> a() =, a() = Note: Onl the bondar condition at = is imposed on the inner epansion at = Sbstitting into eqation - 5 -
16 AMS B Pertrbation Methods ( a + a +! ) ( a + a +! )+ a +! ==> a a + a a + a +!= = : : a a = a = ==> a ( ) = c ( e ) a a = a = c ( e ) = c e + e a ( ) = Let A s = L a ( ) and c = a (). Taking Laplace transform of both sides, we have s A( s) sa ( ) a ( ) sa( s) a ==> ( s s)a s ==> A s = c = c c = c s s + s s s + s s c s s s s s + s = c 5 c s + c s s + s c s s Renaming c = c 5 c and taking inverse transform, we get a ( ) = L A( s) = c ( e )+ c ( +e ) c e e In the above, we have sed L e a = s + a and Ths, the inner epansion is L e a = s + a - 6 -
17 AMS B Pertrbation Methods ( ) = c ( e )+ c ( e )+ c ( +e ) c ( e e ) inn where coefficients c and c are to be determined in matching. Matching: Recall the oter epansion: ( ) = ot The oter limit of inner epansion: ( +) log + + For (, ), as ϵ - we have = (-)/ϵ. Ths, we treat e as zero. ( inn )( ot) ( ) = c + c + c The inner limit of oter epansion: We do Talor epansion arond = and keep terms p to O(). Note that (-) = ϵ = O()). ( ot )( inn) ( ) = ( ) = + 4 log ( ) Van Dke s Matching condition: ( inn )( ot) ( ) = ( ot )( inn) ( ) ==> c + c + c = + 4 log ==> c =, c = log ==> ( m ) inn = + 4 log ( ) ( ) = ( e )+ log ( ) = ( + ) log ( e )+ ( 4 +e ) ( 8 e e ) 4 log ( ) - 7 -
18 AMS B Pertrbation Methods ( ) ( m ) inn The composite epansion is c = e + log ( ) ( ) = ( ot )( )+ ( inn )( ) ( m ) = ( +) e ( ) = + e + log ( +)+ log ( + ) e + e 8 e e e 8 e e ( ) The figre below compares the eact soltion and composite epansions for =...8 = Eact soltion asmptotic with term asmptotic with terms Appendi B Soltion of Let A s a a = sin a ( ) = = L a ( ) and c = a (). Taking Laplace transform of both sides, we have - 8 -
19 AMS B Pertrbation Methods s A( s) sa ( ) a ( ) A s ==> ( s )A s ==> A( s) = c ==> A( s) = c ==> A s = c = c +sin = sin s s +sin s s s s + s + sin s + + s s + sin s + s +sin s + s Renaming c = c + sin and taking inverse transform, we obtain a ( ) = L A( s) = c ( e e )+sin e In the above, we have sed L e a = s + a s Appendi C Soltion of Let A s a a = cos a ( ) = = L a ( ) and c = a (). Taking Laplace transform of both sides, we have s A( s) sa ( ) a ( ) A s ==> ( s )A s = c cos = cos s s - 9 -
20 AMS B Pertrbation Methods ==> A( s) = c ==> A s = cos s cos s s c s + s +cos s Renaming c = cos c and taking inverse transform, we arrive at a ( ) = L A( s) = c ( e e )+cos - -
Chapter 2 Difficulties associated with corners
Chapter Difficlties associated with corners This chapter is aimed at resolving the problems revealed in Chapter, which are cased b corners and/or discontinos bondar conditions. The first section introdces
More information3 2D Elastostatic Problems in Cartesian Coordinates
D lastostatic Problems in Cartesian Coordinates Two dimensional elastostatic problems are discssed in this Chapter, that is, static problems of either plane stress or plane strain. Cartesian coordinates
More informationChapter 6 Momentum Transfer in an External Laminar Boundary Layer
6. Similarit Soltions Chapter 6 Momentm Transfer in an Eternal Laminar Bondar Laer Consider a laminar incompressible bondar laer with constant properties. Assme the flow is stead and two-dimensional aligned
More informationEE2 Mathematics : Functions of Multiple Variables
EE2 Mathematics : Fnctions of Mltiple Variables http://www2.imperial.ac.k/ nsjones These notes are not identical word-for-word with m lectres which will be gien on the blackboard. Some of these notes ma
More informationComplex Variables. For ECON 397 Macroeconometrics Steve Cunningham
Comple Variables For ECON 397 Macroeconometrics Steve Cnningham Open Disks or Neighborhoods Deinition. The set o all points which satis the ineqalit
More information10.4 Solving Equations in Quadratic Form, Equations Reducible to Quadratics
. Solving Eqations in Qadratic Form, Eqations Redcible to Qadratics Now that we can solve all qadratic eqations we want to solve eqations that are not eactl qadratic bt can either be made to look qadratic
More informationSection 7.4: Integration of Rational Functions by Partial Fractions
Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to
More informationMean Value Formulae for Laplace and Heat Equation
Mean Vale Formlae for Laplace and Heat Eqation Abhinav Parihar December 7, 03 Abstract Here I discss a method to constrct the mean vale theorem for the heat eqation. To constrct sch a formla ab initio,
More informationcalled the potential flow, and function φ is called the velocity potential.
J. Szantr Lectre No. 3 Potential flows 1 If the flid flow is irrotational, i.e. everwhere or almost everwhere in the field of flow there is rot 0 it means that there eists a scalar fnction ϕ,, z), sch
More informationSecond-Order Wave Equation
Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order
More informationPartial Differential Equations with Applications
Universit of Leeds MATH 33 Partial Differential Eqations with Applications Eamples to spplement Chapter on First Order PDEs Eample (Simple linear eqation, k + = 0, (, 0) = ϕ(), k a constant.) The characteristic
More information1 Differential Equations for Solid Mechanics
1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to
More informationPhysicsAndMathsTutor.com
C Integration - By sbstittion PhysicsAndMathsTtor.com. Using the sbstittion cos +, or otherwise, show that e cos + sin d e(e ) (Total marks). (a) Using the sbstittion cos, or otherwise, find the eact vale
More informationBLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students
BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they
More informationIncompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System
International Jornal of Compter Applications (97 8887) Volme 79 No., October Incompressible Viscoelastic Flow of a Generalised Oldroed-B Flid throgh Poros Medim between Two Infinite Parallel Plates in
More information4 Exact laminar boundary layer solutions
4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved
More informationPulses on a Struck String
8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a
More informationEXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L.
.4 WAVE EQUATION 445 EXERCISES.3 In Problems and solve the heat eqation () sbject to the given conditions. Assme a rod of length.. (, t), (, t) (, ),, > >. (, t), (, t) (, ) ( ) 3. Find the temperatre
More information1. State-Space Linear Systems 2. Block Diagrams 3. Exercises
LECTURE 1 State-Space Linear Sstems This lectre introdces state-space linear sstems, which are the main focs of this book. Contents 1. State-Space Linear Sstems 2. Block Diagrams 3. Exercises 1.1 State-Space
More informationChapter 1: Differential Form of Basic Equations
MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)
More informationLinear System Theory (Fall 2011): Homework 1. Solutions
Linear System Theory (Fall 20): Homework Soltions De Sep. 29, 20 Exercise (C.T. Chen: Ex.3-8). Consider a linear system with inpt and otpt y. Three experiments are performed on this system sing the inpts
More informationA Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane
Adv. Theor. Appl. Mech., Vol. 1, 8, no. 1, 9 A Decomposition Method for Volme Flx and Average Velocit of Thin Film Flow of a Third Grade Flid Down an Inclined Plane A. Sadighi, D.D. Ganji,. Sabzehmeidani
More informationLecture 17 Errors in Matlab s Turbulence PSD and Shaping Filter Expressions
Lectre 7 Errors in Matlab s Trblence PSD and Shaping Filter Expressions b Peter J Sherman /7/7 [prepared for AERE 355 class] In this brief note we will show that the trblence power spectral densities (psds)
More informationMECHANICS OF SOLIDS COMPRESSION MEMBERS TUTORIAL 2 INTERMEDIATE AND SHORT COMPRESSION MEMBERS
MECHANICS O SOIDS COMPRESSION MEMBERS TUTORIA INTERMEDIATE AND SHORT COMPRESSION MEMBERS Yo shold jdge yor progress by completing the self assessment exercises. On completion of this ttorial yo shold be
More informationHomotopy Perturbation Method for Solving Linear Boundary Value Problems
International Jornal of Crrent Engineering and Technolog E-ISSN 2277 4106, P-ISSN 2347 5161 2016 INPRESSCO, All Rights Reserved Available at http://inpressco.com/categor/ijcet Research Article Homotop
More informationLinear Strain Triangle and other types of 2D elements. By S. Ziaei Rad
Linear Strain Triangle and other tpes o D elements B S. Ziaei Rad Linear Strain Triangle (LST or T6 This element is also called qadratic trianglar element. Qadratic Trianglar Element Linear Strain Triangle
More informationSetting The K Value And Polarization Mode Of The Delta Undulator
LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions
More informationGeometry of Span (continued) The Plane Spanned by u and v
Geometric Description of Span Geometr of Span (contined) 2 Geometr of Span (contined) 2 Span {} Span {, } 2 Span {} 2 Geometr of Span (contined) 2 b + 2 The Plane Spanned b and If a plane is spanned b
More informationMomentum Equation. Necessary because body is not made up of a fixed assembly of particles Its volume is the same however Imaginary
Momentm Eqation Interest in the momentm eqation: Qantification of proplsion rates esign strctres for power generation esign of pipeline systems to withstand forces at bends and other places where the flow
More informationBertrand s Theorem. October 8, µr 2 + V (r) 0 = dv eff dr. 3 + dv. f (r 0 )
Bertrand s Theorem October 8, Circlar orbits The eective potential, V e = has a minimm or maximm at r i and only i so we mst have = dv e L µr + V r = L µ 3 + dv = L µ 3 r r = L µ 3 At this radis, there
More informationChapter 3. Preferences and Utility
Chapter 3 Preferences and Utilit Microeconomics stdies how individals make choices; different individals make different choices n important factor in making choices is individal s tastes or preferences
More informationdenote the space of measurable functions v such that is a Hilbert space with the scalar products
ω Chebyshev Spectral Methods 34 CHEBYSHEV POLYOMIALS REVIEW (I) General properties of ORTHOGOAL POLYOMIALS Sppose I a is a given interval. Let ω : I fnction which is positive and continos on I Let L ω
More informationDifferentiation of Exponential Functions
Differentiation of Eponential Fnctions The net derivative rles that o will learn involve eponential fnctions. An eponential fnction is a fnction in the form of a constant raised to a variable power. The
More informationFormal Methods for Deriving Element Equations
Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact
More informationSUBJECT:ENGINEERING MATHEMATICS-I SUBJECT CODE :SMT1101 UNIT III FUNCTIONS OF SEVERAL VARIABLES. Jacobians
SUBJECT:ENGINEERING MATHEMATICS-I SUBJECT CODE :SMT0 UNIT III FUNCTIONS OF SEVERAL VARIABLES Jacobians Changing ariable is something e come across er oten in Integration There are man reasons or changing
More informationThe New (2+1)-Dimensional Integrable Coupling of the KdV Equation: Auto-Bäcklund Transformation and Non-Travelling Wave Profiles
MM Research Preprints, 36 3 MMRC, AMSS, Academia Sinica No. 3, December The New (+)-Dimensional Integrable Copling of the KdV Eqation: Ato-Bäcklnd Transformation and Non-Traelling Wae Profiles Zhena Yan
More informationIT is well known that searching for travelling wave solutions
IAENG International Jornal of Applied Mathematics 44: IJAM_44 2 A Modification of Fan Sb-Eqation Method for Nonlinear Partial Differential Eqations Sheng Zhang Ao-Xe Peng Abstract In this paper a modification
More informationWall treatment in Large Eddy Simulation
Wall treatment in arge Edd Simlation David Monfort Sofiane Benhamadoche (ED R&D) Pierre Sagat (Université Pierre et Marie Crie) 9 novembre 007 Code_Satrne User Meeting Wall treatment in arge Edd Simlation
More informationConcept of Stress at a Point
Washkeic College of Engineering Section : STRONG FORMULATION Concept of Stress at a Point Consider a point ithin an arbitraril loaded deformable bod Define Normal Stress Shear Stress lim A Fn A lim A FS
More informationA Single Species in One Spatial Dimension
Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,
More informationThis Topic follows on from Calculus Topics C1 - C3 to give further rules and applications of differentiation.
CALCULUS C Topic Overview C FURTHER DIFFERENTIATION This Topic follows on from Calcls Topics C - C to give frther rles applications of differentiation. Yo shold be familiar with Logarithms (Algebra Topic
More informationLecture 11: Fourier Cosine Series
Introductory lecture notes on Partial Differential Equations - c Anthony Peirce Not to be copied, used, or revised without eplicit written permission from the copyright owner ecture : Fourier Cosine Series
More informationMath 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation:
Math 263 Assignment #3 Soltions 1. A fnction z f(x, ) is called harmonic if it satisfies Laplace s eqation: 2 + 2 z 2 0 Determine whether or not the following are harmonic. (a) z x 2 + 2. We se the one-variable
More informationMAT389 Fall 2016, Problem Set 6
MAT389 Fall 016, Problem Set 6 Trigonometric and hperbolic fnctions 6.1 Show that e iz = cos z + i sin z for eer comple nmber z. Hint: start from the right-hand side and work or wa towards the left-hand
More informationThe Dual of the Maximum Likelihood Method
Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All
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 informationPrimary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector
Chapter 4: Flids Kinematics 4. Velocit and Description Methods Primar dependent ariable is flid elocit ector V V ( r ); where r is the position ector If V is known then pressre and forces can be determined
More informationA NEW EVALUATION PLATFORM FOR NAVIGATION SYSTEMS
A NEW EVALUATION PLATFORM FOR NAVIGATION SYSTEMS Thomas Hanefeld Sejerøe Niels Kjølstad polsen Ole Ravn Ørsted DTU, Atomation, The Technical Universit of Denmar, Bilding 6, DK-8 Kgs. Lngb, Denmar E-mail:9977@stdent.dt.d,
More informationMotion in an Undulator
WIR SCHAFFEN WISSEN HEUTE FÜR MORGEN Sven Reiche :: SwissFEL Beam Dnamics Grop :: Pal Scherrer Institte Motion in an Undlator CERN Accelerator School FELs and ERLs On-ais Field of Planar Undlator For planar
More informationOptimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications
Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications Navin Khaneja lectre notes taken by Christiane Koch Jne 24, 29 1 Variation yields a classical Hamiltonian system Sppose that
More informationNonlocal Symmetries and Interaction Solutions for Potential Kadomtsev Petviashvili Equation
Commn. Theor. Phs. 65 (16) 31 36 Vol. 65, No. 3, March 1, 16 Nonlocal Smmetries and Interaction Soltions for Potential Kadomtsev Petviashvili Eqation Bo Ren ( ), Jn Y ( ), and Xi-Zhong Li ( ) Institte
More information4753 Mark Scheme June 2017 Question Answer Marks Guidance. x x x A1 correct expression, allow. isw
475 Mark Scheme Jne 7 Qestion Answer Marks Gidance (5 ) (5 ) d Chain rle on (5 ) (5 ) ( 6 )( )(5 ) crect epression, allow ( 6 )( ) o.e. (5 ) isw cao isw (5 ) inverted v shape throgh (, ), (, ) and (, )
More information10.2 Solving Quadratic Equations by Completing the Square
. Solving Qadratic Eqations b Completing the Sqare Consider the eqation ( ) We can see clearl that the soltions are However, What if the eqation was given to s in standard form, that is 6 How wold we go
More informationHigher Maths A1.3 Recurrence Relations - Revision
Higher Maths A Recrrence Relations - Revision This revision pack covers the skills at Unit Assessment exam level or Recrrence Relations so yo can evalate yor learning o this otcome It is important that
More informationApproximate Solution of Convection- Diffusion Equation by the Homotopy Perturbation Method
Gen. Math. Notes, Vol. 1, No., December 1, pp. 18-114 ISSN 19-7184; Copyright ICSRS Pblication, 1 www.i-csrs.org Available free online at http://www.geman.in Approximate Soltion of Convection- Diffsion
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationUncertainties of measurement
Uncertainties of measrement Laboratory tas A temperatre sensor is connected as a voltage divider according to the schematic diagram on Fig.. The temperatre sensor is a thermistor type B5764K [] with nominal
More informationSTABILIZATIO ON OF LONGITUDINAL AIRCRAFT MOTION USING MODEL PREDICTIVE CONTROL AND EXACT LINEARIZATION
8 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES STABILIZATIO ON OF LONGITUDINAL AIRCRAFT MOTION USING MODEL PREDICTIVE CONTROL AND EXACT LINEARIZATION Čeliovsý S.*, Hospodář P.** *CTU Prage, Faclty
More informationShock wave structure for Generalized Burnett Equations
Shock wave strctre for Generalized Brnett Eqations A.V. Bobylev, M. Bisi, M.P. Cassinari, G. Spiga Dept. of Mathematics, Karlstad University, SE-65 88 Karlstad, Sweden, aleander.bobylev@ka.se Dip. di Matematica,
More informationConditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane
Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching
More information4.4 Moment of a Force About a Line
4.4 Moment of a orce bot a Line 4.4 Moment of a orce bot a Line Eample 1, page 1 of 3 1. orce is applied to the end of gearshift lever DE. Determine the moment of abot shaft. State which wa the lever will
More informationGeometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018
Geometric Image Maniplation Lectre 4 Wednesda, Janar 4, 08 Programming Assignment Image Maniplation: Contet To start with the obvios, an image is a D arra of piels Piel locations represent points on the
More informationIntegration of Basic Functions. Session 7 : 9/23 1
Integration o Basic Fnctions Session 7 : 9/3 Antiderivation Integration Deinition: Taking the antiderivative, or integral, o some nction F(), reslts in the nction () i ()F() Pt simply: i yo take the integral
More informationA Computational Study with Finite Element Method and Finite Difference Method for 2D Elliptic Partial Differential Equations
Applied Mathematics, 05, 6, 04-4 Pblished Online November 05 in SciRes. http://www.scirp.org/jornal/am http://d.doi.org/0.46/am.05.685 A Comptational Stdy with Finite Element Method and Finite Difference
More informationLab Manual for Engrd 202, Virtual Torsion Experiment. Aluminum module
Lab Manal for Engrd 202, Virtal Torsion Experiment Alminm modle Introdction In this modle, o will perform data redction and analsis for circlar cross section alminm samples. B plotting the torqe vs. twist
More informationAn alternative approach to evaluate the average Nusselt number for mixed boundary layer conditions in parallel flow over an isothermal flat plate
An alternative approach to evalate the average Nsselt nber for ied bondary layer conditions in parallel flo over an isotheral flat plate Viacheslav Stetsyk, Krzysztof J. Kbiak, ande i and John C Chai Abstract
More informationρ u = u. (1) w z will become certain time, and at a certain point in space, the value of
THE CONDITIONS NECESSARY FOR DISCONTINUOUS MOTION IN GASES G I Taylor Proceedings of the Royal Society A vol LXXXIV (90) pp 37-377 The possibility of the propagation of a srface of discontinity in a gas
More informationRange of validity and intermittent dynamics of the phase of oscillators with nonlinear self-excitation
Range of validit and intermittent dnamics of the phase of oscillators with nonlinear self-ecitation D.V. Strnin 1 and M.G. Mohammed Comptational Engineering and Science Research Centre, Faclt of Health,
More informationSTEP Support Programme. STEP III Hyperbolic Functions: Solutions
STEP Spport Programme STEP III Hyperbolic Fnctions: Soltions Start by sing the sbstittion t cosh x. This gives: sinh x cosh a cosh x cosh a sinh x t sinh x dt t dt t + ln t ln t + ln cosh a ln ln cosh
More informationThe Scalar Conservation Law
The Scalar Conservation Law t + f() = 0 = conserved qantity, f() =fl d dt Z b a (t, ) d = Z b a t (t, ) d = Z b a f (t, ) d = f (t, a) f (t, b) = [inflow at a] [otflow at b] f((a)) f((b)) a b Alberto Bressan
More informationMEG 741 Energy and Variational Methods in Mechanics I
MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)
More informationOn the scaling of entropy viscosity in high order methods
On te scaling of entropy viscosity in ig order metods Adeline Kornels and Daniel Appelö Abstract In tis work, we otline te entropy viscosity metod and discss ow te coice of scaling inflences te size of
More information3.4-Miscellaneous Equations
.-Miscellaneos Eqations Factoring Higher Degree Polynomials: Many higher degree polynomials can be solved by factoring. Of particlar vale is the method of factoring by groping, however all types of factoring
More informationDiscontinuous Fluctuation Distribution for Time-Dependent Problems
Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation
More informationA Macroscopic Traffic Data Assimilation Framework Based on Fourier-Galerkin Method and Minimax Estimation
A Macroscopic Traffic Data Assimilation Framework Based on Forier-Galerkin Method and Minima Estimation Tigran T. Tchrakian and Sergiy Zhk Abstract In this paper, we propose a new framework for macroscopic
More informationBoundary layer develops in the flow direction, δ = δ (x) τ
58:68 Trblent Flos Handot: Bondar Laers Differences to Trblent Channel Flo Bondar laer develops in the flo direction, not knon a priori Oter part of the flo consists of interittent trblent/non-trblent
More informationfür Mathematik in den Naturwissenschaften Leipzig
Ma-Planck-Institt für Mathematik in den Natrwissenschaften Leipzig Nmerical simlation of generalized KP type eqations with small dispersion by Christian Klein, and Christof Sparber Preprint no.: 2 26 NUMERICAL
More informationElastico-Viscous MHD Free Convective Flow Past an Inclined Permeable Plate with Dufour Effects in Presence of Chemical Reaction
International Jornal of Engineering and Technical Research (IJETR) ISSN: 2321-0869 (O) 2454-4698 (P), Volme-3, Isse-8, Agst 2015 Elastico-Viscos MHD Free Convective Flow Past an Inclined Permeable Plate
More information5.5 U-substitution. Solution. Z
CHAPTER 5. THE DEFINITE INTEGRAL 22 5.5 U-sbstittion Eample. (a) Fin the erivative of sin( 2 ). (b) Fin the anti-erivative cos( 2 ). Soltion. (a) We se the chain rle: sin(2 )=cos( 2 )( 2 ) 0 =cos( 2 )2
More informationTheorem (Change of Variables Theorem):
Avance Higher Notes (Unit ) Prereqisites: Integrating (a + b) n, sin (a + b) an cos (a + b); erivatives of tan, sec, cosec, cot, e an ln ; sm/ifference rles; areas ner an between crves. Maths Applications:
More informationUNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL
8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a
More information4.7. Newton s Method. Procedure for Newton s Method HISTORICAL BIOGRAPHY
4. Newton s Method 99 4. Newton s Method HISTORICAL BIOGRAPHY Niels Henrik Abel (18 189) One of the basic problems of mathematics is solving equations. Using the quadratic root formula, we know how to
More informationSolving the Lienard equation by differential transform method
ISSN 1 746-7233, England, U World Jornal of Modelling and Simlation Vol. 8 (2012) No. 2, pp. 142-146 Solving the Lienard eqation by differential transform method Mashallah Matinfar, Saber Rakhshan Bahar,
More informationModelling by Differential Equations from Properties of Phenomenon to its Investigation
Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University
More informationLecture Notes: Finite Element Analysis, J.E. Akin, Rice University
9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)
More information{ } is an asymptotic sequence.
AMS B Perturbation Methods Lecture 3 Copyright by Hongyun Wang, UCSC Recap Iterative method for finding asymptotic series requirement on the iteration formula to make it work Singular perturbation use
More informationAlgebraic Multigrid. Multigrid
Algebraic Mltigrid We re going to discss algebraic ltigrid bt irst begin b discssing ordinar ltigrid. Both o these deal with scale space eaining the iage at ltiple scales. This is iportant or segentation
More informationThe Oscillatory Stable Regime of Nonlinear Systems, with two time constants
6th WSES International Conference on CIRCUITS SYSTEMS ELECTRONICSCONTROL & SIGNL PROCESSING Cairo Egpt Dec 9-3 7 5 The Oscillator Stable Regime of Nonlinear Sstems with two time constants NUŢU VSILE *
More informationPolymer confined between two surfaces
Appendix 4.A 15 Polymer confined between two srfaces In this appendix we present in detail the calclations of the partition fnction of a polymer confined between srfaces with hard wall bondary conditions.
More informationAssignment Fall 2014
Assignment 5.086 Fall 04 De: Wednesday, 0 December at 5 PM. Upload yor soltion to corse website as a zip file YOURNAME_ASSIGNMENT_5 which incldes the script for each qestion as well as all Matlab fnctions
More informationFLUID MECHANICS D203 SAE SOLUTIONS TUTORIAL 2 APPLICATIONS OF BERNOULLI SELF ASSESSMENT EXERCISE 3
FLUI MEHNIS 0 SE SOLUTIONS TUTORIL PPLITIONS OF ERNOULLI SELF SSESSMENT EXERISE Take the density of water to be 997 kg/m throghot nless otherwise stated.. Ventri meter is 50 mm bore diameter at inlet and
More informationChaotic and Hyperchaotic Complex Jerk Equations
International Jornal of Modern Nonlinear Theory and Application 0 6-3 http://dxdoiorg/0436/ijmnta000 Pblished Online March 0 (http://wwwscirporg/jornal/ijmnta) Chaotic and Hyperchaotic Complex Jerk Eqations
More informationApproximate Solution for the System of Non-linear Volterra Integral Equations of the Second Kind by using Block-by-block Method
Astralian Jornal of Basic and Applied Sciences, (1): 114-14, 008 ISSN 1991-8178 Approximate Soltion for the System of Non-linear Volterra Integral Eqations of the Second Kind by sing Block-by-block Method
More informationAdvanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski
Advanced topics in Finite Element Method 3D trss strctres Jerzy Podgórski Introdction Althogh 3D trss strctres have been arond for a long time, they have been sed very rarely ntil now. They are difficlt
More informationLecture 3Section 7.3 The Logarithm Function, Part II
Lectre 3Section 7.3 The Logarithm Fnction, Part II Jiwen He Section 7.2: Highlights 2 Properties of the Log Fnction ln = t t, ln = 0, ln e =. (ln ) = > 0. ln(y) = ln + ln y, ln(/y) = ln ln y. ln ( r) =
More informationA Radial Basis Function Method for Solving PDE Constrained Optimization Problems
Report no. /6 A Radial Basis Fnction Method for Solving PDE Constrained Optimization Problems John W. Pearson In this article, we appl the theor of meshfree methods to the problem of PDE constrained optimization.
More informationMath 273b: Calculus of Variations
Math 273b: Calcls of Variations Yacob Kreh Homework #3 [1] Consier the 1D length fnctional minimization problem min F 1 1 L, or min 1 + 2, for twice ifferentiable fnctions : [, 1] R with bonary conitions,
More information1 Undiscounted Problem (Deterministic)
Lectre 9: Linear Qadratic Control Problems 1 Undisconted Problem (Deterministic) Choose ( t ) 0 to Minimize (x trx t + tq t ) t=0 sbject to x t+1 = Ax t + B t, x 0 given. x t is an n-vector state, t a
More informationHomogeneous Liner Systems with Constant Coefficients
Homogeneos Liner Systems with Constant Coefficients Jly, 06 The object of stdy in this section is where A is a d d constant matrix whose entries are real nmbers. As before, we will look to the exponential
More informationIntroduction to Differential Equations
Math0 Lecture # Introduction to Differential Equations Basic definitions Definition : (What is a DE?) A differential equation (DE) is an equation that involves some of the derivatives (or differentials)
More information