Unformatted text preview: Lecture Notes for EGN 5422 Engineering Applications of Partial Differential Equations Summer and Spring Semesters Arthur David Snider, Professor Emeritus [email protected] The current syllabus and the takehome exams are available from the course Blackboard site. CONTENTS Generic Syllabus (Download a current version from the Blackboard site) Accurate Calculations Some Solutions Some Suggested Assignments Readings to Go with the Lectures Template for Solutions to Laplace's Equation in a Square Isotherms D'Alembert's Solution Answers for Section 5.2 Some Fourier demo sites Some Chapter 7 and Chapter 8 Solutions via USFKAD Guide to the lecture movies Lecture Notes Generic Syllabus: EGN 5422 or EEL 6935 ENGINEERING APPLICATIONS OF PARTIAL DIFFERENTIAL EQUATIONS The 5000 and 6000level sections are identical in content, but the 6000level section is graded at a higher standard. Contact the Electrical Engineering Department if you wish to enroll in a 6000level section. Instructor: Prof. A. D. (Dave) Snider, ENB 246A, 8139744785, FAX: 813/9745250 [email protected] Office hours (see Blackboard) Text: Partial Differential Equations: Sources and Solutions, A. D. Snider, Dover Publications, Mineola NY, 2006, ISBN 0486 453405, required. NOTE: Computer access (PC) is also required. Catalog Description: Power series solutions for ordinary differential equations, SturmLiouville theory, special functions. Separation of variables for partial differential equations. Green's functions. Use of USFKAD software. Course Prerequisites: MAP 2302 Differential Equations (solution methods for first order nonlinear and higher order linear equations) The lectures can be viewed (by anyone at any time) at http://netcast.usf.edu/browse.php?page=Classes/engineering/snider/PartialDiffEqs However, there is an error in the listings: there is NO movie "Lecture 48_3D Examples". Simply skip it. Welcome to my classes, USF students! You're taking a "webonly" class from me. That means you won't see me live; you'll watch lectures I taped in the past, on the course materials. I will be available in my office for consultation a few hours per week, and we can Skype. I'll answer your emails promptly, unless they are addressed in the syllabus or at the course web site in Blackboard, https://my.usf.edu. This site contains class announcements, documents, old and current exams, assignments, etc. Think about this: you're taking a course from a professor whom you'll hardly ever see. So you don't have much chance for an "A" unless you (1) buy the textbook and (2) follow the web procedures precisely. Be sure to read this longwinded syllabus entirely. Then print it out and stick it in the back of your textbook for future reference. You may have misgivings about nonlive lectures. Let me relate a few observations that I have made over the past years with this process. First of all, this procedure has proved quite successful to my surprise. I was leery about it at first, but I soon found out that the average student performance on tests was slightly better than it had been before, when I was lecturing live. I think the reasons for this are as follows: 1. You will never need to miss a lecture due to illness, conflicting appointments, being out of town, or simply being tired. The lectures are at the web site all semester long, 24/7. 2. You can rewatch any lecture, or part of a lecture, as many times as you need. 3. If you need to take a break (no student has ever fallen asleep during my lectures, of course!!!), you can stop the playback and resume when you return. The lectures are numbered in order (except for lecture 51, Numerics, which can be watched anytime after you have mastered the software and before you attempt the takehome). They were taped during different semesters, so you will have to figure out the right pace to get you through them during the time allotted in your current semester. If you can't do this, please drop my class and retake MTH 101, paying particular attention to the lessons on ratio and proportion. Any administrative instructions that I give in the lectures are probably out of date. Consult your email and Blackboard weekly for updates from me on assignments, test dates, and procedures in general. Now let me take one paragraph to try to talk you out of taking this course. Partial differential equations is a very difficult subject. In my book and lectures I have tried to make it as simple as possible; but still EGN 5422 is the hardest course I teach. Only about 5060% of the students who take it make an A; you need to be very good at mathematics to achieve an A. Don't take this course if you can't afford a B on your transcript; if you need to know the subject, consider auditing the course for noncredit. Be sure to read the document on accurate calculations. Remember these tips in all your classes. Instructions for using MATLAB to perform the integrations for the takehome assignment are covered in this document. Lecture notes replicating everything that appears on the blackboard during the lectures can be purchased from > ProCopy > 5219 E. Fowler Ave. > Tampa, Fl 33617 > [email protected]<mailto:[email protected]> > (813) 9885900 They can also be downloaded from Blackboard. (Ignore notes in the .xbk format.) If you download, print out the lecture notes before you watch the lecture, and have them in front of you.
The software USFKAD can be downloaded from Blackboard. LaTeX software is required to run USFKAD. LaTeX (the MiKTeX package) is available at www.miktex.org . Download LaTeX during the first week and read the notes at Blackboard for installing it. Experience has shown that if you wait until midsemester to start working with LaTeX, you will sabotage any chance of success in this course. Test your installation by trying to open and read the article CompPhys4.tex . Requirements and Assessment: 1. Each student must email Prof. Snider with the following data: Last name: ______ First name: _______ Class: EGN 5422, by (see Blackboard). I will not acknowledge these emails individually. In one week you will receive an acknowledgement, by email, from me (Dr. Snider) that you are in his class email address list; if you do not receive this acknowledgement, email me again until I acknowledge receipt. Thereafter each student is liable for all email notices concerning the class from Prof. Snider. Students who wish to use different email boxes should email this data from each box. Do not use one email box to request mail to a different box. 2. Each student must sign a copy of the final page of this syllabus as indicated below and submit it to Dr. Snider by (see Blackboard). You are not officially enrolled in the class until you have turned in a signed syllabus. Postalmail a hard copy to Dr. A D Snider, Dept. of Electrical Engineering, University of South Florida, 4202 East Fowler Avenue ENB 118, Tampa FL 33620; or put a copy in my EE Department mailbox. Email is not acceptable. I will not acknowledge receipt of these syllabi individually, but they can be regenerated later in the semester if a problem arises. 3. Certain homework problems will be recommended to the students, but not graded. You should regard the old tests as a prime source of homework problems; work them during the semester as the particular topic is covered in the lectures. Once you have mastered the USFKAD software, you can use it to confirm your answers. 4. A midterm examination (time and place to be arranged, based on section 5.2 only) and a final will be given. Note that answers to all the problems in Section 5.2 are given in Blackboard's "Course Documents" A time and place on the Tampa campus will be arranged for these tests, but they can be taken at remote sites and more convenient times (within limits) if a proctor agreement is worked out; I'll notify you. Additionally a takehome test, based lecture #51 "Numerics," will be assigned. The takehome is heavily computational and will be computer graded with no partial credit, but you will be allowed four attempts (with, however, different numerical parameters each time). You will submit your answers by email to the TA, who will shortly respond with your score, a tabulation of your incorrect answers, and the correct answers for the parameters you used. Your first attempt is assessed at (only) 5%, the second at 55%, the third at 20%, and the fourth at 20%. Each attempt has a deadline. If you miss a deadline, your subsequent submission will also count as the missed submission. If you stop submitting at any time, your last submission will count for the subsequent ones. The takehome tests will be available at Blackboard, and they contain more detailed instructions. Time permitting, you will be asked to present an oral defense of your takehome and exam solutions. Your final grade will be a weighted average of the midterm (20%), the takehome (40%), and the final (40%). I recommend that you take a timed midterm and a timed final from Blackboard for practice. These tests are openbook, closed notes; you are not permitted to bring old tests to the exams. (I'm not trying to be tough; experience has shown that old tests are counterproductive rather than helpful.) 5. An "incomplete" grade will be awarded if either the email, syllabus signoff, midterm, takehome, or final are not submitted. Incompletes can be made up following USF policies; you are on your own to figure out these policies. However you need to know that I am retired, and only work as an adjunct faculty, so I cannot guarantee my availability after the semester is over. If you need a flexible, cooperative professor instead of a cranky old curmudgeon, you are advised to take another class. Please mail a copy of this page to Dr. Snider. Academic Dishonesty  It is not acceptable to copy, plagiarize or otherwise make use of the work of others in completing homework, project, exam or other course assignments. The minimum penalty for doing so is an automatic zero on the assignment and an "F" in the course. I have read the syllabus for EGN 5422, Spring 2011, and agree to abide by its schedule and terms. Print name: Sign name: Accurate calculations and evaluation of integrals using MATLAB. Tips on calculating: Some tests in this course require a lot of calculation. Here are some pointers that you may not have picked up yet: 1. If you require a certain number of significant digits in a final answer, you must maintain more digits in the intermediate calculations. For instance consider the addition 1/3 + 1/3 = 2/3. The threesignificant digit expression for the answer, 2/3, is .667 . But if you first round the addends to three significant digits you get .333 + .333 = .666, which is not correct to 3 digits. A good rule of thumb is to retain at least two more digits than required, in all intermediate calculations. 2. Always introduce symbols (letters) when you evaluate a formula, especially on a calculator. For instance suppose you wanted to add 436 and 578 on a calculator. The worst way to do it is to enter 436, press enter, press +, enter 578, press enter, and press =. Because: suppose you get an answer you know is wrong, like 23. You don't know whether you entered the 436 wrong, or the 578, or the + sign. You have to reenter all the data again. The smart way is to let A equal 436, let B equal 578, and call for A+B. If the answer is absurd, you can recall A, B, and the formula; and you can correct only the one that's wrong. This is particularly significant when you're dealing with highdigit numbers, and complicated formulas that may require parentheses. Remember these tips in all your classes. 3. To add an infinite series. Unless you know a lot about the series, there is no sure way to tell when you have achieved a specified accuracy. Here's what conservative engineers do in practice: Add up a selected number of terms 5, 10, whatever. Then double the number of terms and add again. Compare the second total to the first; if they agree to within the prescribed accuracy, you're probably OK. (Naturally, take the second total as your answer.) If the required accuracy is expressed as a percentage, the difference between the two totals must be less than the specified percentage of the final total. If they don't agree sufficiently well, double the number of terms and add again, and so on until you get agreement in two consecutive sums, to within the specified accuracy. If you can't get agreement with your resources, you're going to need a professional mathematician. 4. Note: if you don't know MATLAB, you still should be able to use it to evaluate integrals; just open it (click on the MATLAB icon) and notice how functions are entered into the "@" statements below multiplication, sines, exponentials, division, CORRECT NUMBER OF PARENTHESES. 5.678 I. Onedimensional integrals. To integrate 1.234 [e 2 x cos(4 x 7) 2 2x 1 ] dx : ( x 3) ( x 1)( x 2)
2 in MATLAB type format long (Enter) a=2 (Enter) b=4 (Enter) c=7 (Enter) d=2 (Enter) e=1 (Enter) f=3 (Enter) g=1 (Enter) h=2 (Enter) lower=1.234 (Enter) upper=5.678 (Enter) accuracy = 1e8 (Enter) (Of course you may want more or less accuracy.) Now enter the formula for the integrand, REGARDING THE VARIABLE OF INTEGRATION (x) AS A MATRIX. That means inserting the period mark (.) in front of all multiplication signs *, division signs /, and exponentiation carats ^ involving x. [email protected](x)(exp(a*x.^2).*cos(b.*x+c) (d.*x+e)./((x+f).^2.*(x+g).*(x+h))) (Enter) (You'll probably get some typo error message here; retype carefully, counting parentheses and inserting dots. An error message like "??? Error using ==> mtimes Inner matrix dimensions must agree." means you left out some dots. An error message like " Error: Expression or statement is incorrectpossibly unbalanced (, {, or [." usually means your parentheses are wrong. If you change the values of any of the parameters, reenter the [email protected] statement.) Now type Q=quadl(integrand, lower, upper, accuracy) (Enter) I got an answer of Q = 0.037041299610442 . II. Twodimensional integrals. To integrate (u
4 6 5 7 2 uv) dvdu : in MATLAB type format long (Enter) lowerv=6 (Enter) upperv=7 (Enter) loweru=4 (Enter) upperu=5 (Enter) accuracy = 1e8 (Enter) [email protected](u,v)(u.^2+u.*v) (Enter) Q= dblquad(integrand, loweru,upperu,lowerv,upperv, accuracy, @quadl) (Enter) I got Q = 49.583333333333329 III. Threedimensional integrals. To integrate uvw dw dv du : 2 4 6 3 5 7 in MATLAB type format long (Enter) loww=6 (Enter) highw=7 (Enter) lowv=4 (Enter) highv=5 (Enter) lowu=2 (Enter) highu=3 (Enter) accuracy = 1e8 (Enter) [email protected](u,v,w)(u.*v.*w) (Enter) Q= triplequad(integrand, lowu,highu,lowv,highv,loww,highw,accuracy, @quadl) (Enter) I got Q = 73.124999999999986 . To be conservative, try out your code with integrands that you can do by hand, to check. Some Solutions 1 x1 + O(x0) Assignments. Don't hand in. 1. One solution of the equation y''  4y = sin x is 0.2 sin x . Use shifted cosh and sinh to find a solution meeting the initial conditions y(3) = 2, y'(3) = 1 . Hand in your answer on one sheet of paper. Answer y(x) = [2 + .2 sin 3] cosh 2(x3) + 0.5 [1 + .2 cos 3] sinh 2(x3)  .2 sin x or 2.0282 cosh 2(x3) + .401 sinh 2(x3)  .2 sin x If you got 2.0105 cosh 2(x3) + .599 sinh 2(x3)  .2 sin x , you used degrees, not radians. 2. Page 15, # 1[a,b,c,e,f], #2 Page 30 #10[a,b,c] Page 82 #2 3. Use the software to solve Example 1, p. 311 Example 1, p. 311, with the Neumann conditions changed to Dirichlet conditions Laplace's equation in a square, homogen. Dirichlet on the left and bottom, homog. Neumann on the right, nonhomog. Dirichlet on the top. Example 1, p. 333 4. #1 on page 303. (This equation is not included in USFKAD.) For those of you who are curious about the nonexistence of a steady state solution to the Neumann problem, it is discussed at length in #10, page 304. It can be really nasty in engineering situations! 5. Page 315: 1 and 3
Page 303: 4. Readings for EGN 5422, from PDEs: Sources and Solutions by A D Snider Readings to Go with the Lectures Jan 6 (Week 1A): 1.1 Jan 8 (Week 1B): 1.1 Jan 13 (Week 2A): 1.2 Jan 15 (Week 2B): 2.2, 2.3; also Nagel, Saff, and Snider, Fundamentals of DEs, "Qualitative Considerations for VariableCoefficient and Nonlinear Equations" Jan 20 (Week 3A): 2.4, 2.5 Jan 22 (Week 3B): 2.5, 3.13.4 Jan 27: 3.5, 4.3 Jan 29: 4.3 Feb 3: 4.34.5 Feb 5: 4.6 Feb 10: 5.1 Feb 12: 5.2 Feb 24: 5.2, 5.3 Feb 26: 5.3 Mar 2: 6.1 Mar 4: 6.2 Mar 16: 6.3 Mar 18: 6.2, 6.3 Mar 23: 6.3 Mar 25: 6.3, 6.4 Mar 30: 6.4 Apr 1: 7.1 Apr 6: 7.2 Apr 8: 7.3 Apr 13: 1.3, 1.4, 8.1, 8.2 Apr 15: 8.2 Template for solutions to Laplace's equation in a square, in 2 dimensions. [a1 cosh kx + a2 sinh kx] [d1 cos ky + d2 sin ky] + [b1 + b2x] [e1 + e2y] + [c1 cos kx + c2 sin kx] [f1 cosh ky + f2 sinh ky] Isotherms 50 50 25 Dave's Solution 0 10 75 100 100 125 140 150 50
60 70 Philip's Solution
80 90 25 50 0 60 75
70 80 90 100 100 125 150 D'Alembert's Solution Answers for Section 5.2 sinh 3 Some Fourier Demo Sites http://www.math.ethz.ch/~lanford/MMP/Fourier/demo3.html http://www.falstad.com/fourier/ http://www.jhu.edu/~signals/fourier2/ http://www.math.ubc.ca/~feldman/demos/demo3.html http://homepages.gac.edu/~huber/fourier/ Some Solutions from Chapters 7 and 8, via USFKAD page 439 ܵ ݷ ݵ ʼ ݾ ݵ ݵ ݵ ܷ վ ܷ վ ݷ ݷ ʼ ݵ ݵ ܷ ݵ page 444 Ƚ ʼ ֵ Ƚ ֵ ֵ ֵ Ƚ ʼ ֵ page 447 ʼ ֵ ֵ ֵ Ҵ ֵ ޢ page 449 Ƚ Ѵ page 453 ܵ ݵ ܵ ŵ ۵ ʼ ܵ ܵ ʼ ܵ ʼ page 456 ܵ page 457 Ƚ Ѵ Ƚ ֵ ֵ Ѵ ֵ page 459 Fig 6 ܵ ܵ page 459 Fig 7 ܵ ܵ page 464 ܵ ܵ page 467 ʼ ܷ ʼ page 470 ʼ ܷ ݷ ʼ page 471 ʽ ʼ ʽ ʽ ܷ ݷ ʼ ʽ page 474 page 476 Ƚ ʽ ֵ Ƚ ʽ ַ ֵ page 477 Ƚ ʽ ֵ ֵ ¾ ¾ Ƚ ʽ ֵ ֵ ַ ֽ ֵ ¾ ¾ page 478 Ƚ ֵ ֵ ޣ ֵ Ƚ ֵ ֵ ޣ ʼ ֵ page 514 ݴ ݵ ܵ ܴ ݵ ܵ ܵ ʼ ܵ ܵ ʼ ܵ ʼ page 523 Ƚ Ƚ ֵ ޣ ֵ ֵ page 527 page 529 page 530 ״ page 531 Ƚ
ֵ Խ Ƚ page 534 ʼ page 554 ܷ ʼ ݵ ʼ page 557 ʼ ݵ page 562 Ƚ ֵ Ѵ ѵ ѵ ...
View
Full Document
 Fall '11
 DavidSnider
 Partial differential equation, Arthur David Snider, ÓØ ÖÛ, ©ØÖ Ò×

Click to edit the document details