Implicit Methods:
The Implicit Difference Equations:
By approximating uxx and ut at tj+1 rather than tj as well as using a backwards difference for ut, the equation ut= cuxx is approximated by:
Note that every the terms have index j + 1 except one and isolating this term leads to:
ui,j = −rui−1,j+1 + (1 + 2r)ui,j+1− rui+1,j+1 for 1 ≤ i ≤ m − 1,
Where r = ck/h2 as before.
At this time we have ujgiven in terms of uj+1. This would seem similar to a problem until you consider that the relationship is linear. Utilizing matrix notation we have:
uj= Buj+1− rbj+1,
Where bj+1 represents the boundary situation. Therefore to find uj+1 we need only solve the linear system:
Buj+1 = uj+ rbj+1,
Where uj and bj+1 are given and
Using this scheme is called the implicit method since uj+1 is defined implicitly.
Since we are solving the most significant quantity is the maximum absolute eigenvalue of B−1 which is 1 divided by the smallest ew of B. The maximum absolute ew’s of B−1 as a function of r for different size matrices. Notice that this complete maximum is always less than 1. Therefore errors are always diminished over time and thus this method is always stable. For the same reason it is as well always as accurate as the individual steps.
Both this implicit method as well as the explicit method in the previous lecture make O(h2) error in approximating uxxand O(k) error in approximating ut therefore they have total error O(h2 + k). Therefore although the stability condition permits the implicit method to use arbitrarily large k to maintain accuracy we still need k ∼h2.
Maximum complete eigenvalue EW as a function of r for the matrix B−1 from the implicit method for the heat equation calculated for matrices B of sizes m = 2 . . . 10. When EW <1 the method is stable that is, it is always stable.
Crank-Nicholson Method:
Now that we encompass two different methods for solving parabolic equation it is natural to ask ‘can we get better by taking an average of the two methods?’ The respond is yes.
We execute a weighted average of the two methods by considering an average of the approximations of uxx at j and j + 1. This direct to the equations:
The implicit process contained in these equations is called the Crank-Nicholson method. Gathering terms acquiesce the equations:
−rλui−1,j+1 + (1 + 2rλ)ui,j+1 − rλui+1,j+1 = r(1 − λ)ui−1,j+ (1 − 2r(1 − λ))ui,j+ r(1 − λ)ui+1,j.
In matrix notation this is:
Bλuj+1 = Aλuj + rbj+1,
Where
and
In this equation uj and bj+1 are known, Aujcan be compute directly as well as then the equation is solved for uj+1.
If we choose λ = 1/2 then we are in consequence doing a central difference for ut which has error O(k2). Our total error is then O(h2 + k2). With a bit of work we are able to show that the method is always stable and so we can use k ∼h without a problem.
To get best possible accuracy with a weighted average it is always necessary to use the right weights. For the Crank-Nicholson process with a given r we need to choose:
λ =(r − 1/6)/ 2r
This choice will create the method have truncation error of order O(h4 + k2) which is really excellent considering that the implicit and explicit methods each have truncation errors of order O(h2 + k). Surprisingly we are able to do even better if we also requireand consequently:
With these options the method has truncation error of order O(h6) which is absolutely amazing. To appreciate the proposition supposes that we need to solve a problem with 4 significant digits. If we utilize the explicit or implicit method alone then we will need h2≈ k ≈ 10−4. If L = 1 and T ≈ 1 then we require m ≈ 100 and n ≈ 10000. Therefore we would have a total of 1000000 grid points approximately all on the interior. This is plenty.
Next presume we solve the same problem using the optimal Crank-Nicholson method. We would need h6≈ 10−4 which would need us to take m ≈ 4.64 therefore we would take m = 5 and have h = 1/5. For k we require k = (√5/10)h2/c. If c = 1 this gives k = √5/250 ≈ 0.0089442 so we would need n ≈ 112 to get T ≈ 1. This offers us a total of 560 interior grid points or a factor of 1785 fewer than the explicit or else implicit method alone.
Latest technology based Matlab Programming Online Tutoring Assistance
Tutors, at the www.tutorsglobe.com, take pledge to provide full satisfaction and assurance in Matlab Programming help via online tutoring. Students are getting 100% satisfaction by online tutors across the globe. Here you can get homework help for Matlab Programming, project ideas and tutorials. We provide email based Matlab Programming help. You can join us to ask queries 24x7 with live, experienced and qualified online tutors specialized in Matlab Programming. Through Online Tutoring, you would be able to complete your homework or assignments at your home. Tutors at the TutorsGlobe are committed to provide the best quality online tutoring assistance for Matlab Programming Homework help and assignment help services. They use their experience, as they have solved thousands of the Matlab Programming assignments, which may help you to solve your complex issues of Matlab Programming. TutorsGlobe assure for the best quality compliance to your homework. Compromise with quality is not in our dictionary. If we feel that we are not able to provide the homework help as per the deadline or given instruction by the student, we refund the money of the student without any delay.
Refraction of Light tutorial all along with the key concepts of Refraction of Light through Rectangular Glass, Snell's law, Refraction of Light through Prism, Laws of Refraction, Real and Apparent Depths, total Internal Reflection
Theory and lecture notes of Chi-square test for independence all along with the key concepts of chi-square test for independence, homework help, assignment help. Tutorsglobe offers homework help, assignment help and tutor’s assistance on Chi-square test for independence.
www.tutorsglobe.com offers irreversible addition reactions homework help, irreversible addition reactions assignment help, online tutoring assistance, organic chemistry solutions by online qualified tutor's help.
theory and lecture notes of ttl nand gate all along with the key concepts of ttl nand gate circuit structure, logical operation and circuit analysis. tutorsglobe offers homework help, assignment help and tutor’s assistance on ttl nand gate.
Has your search for the most affordable & trustworthy Strategic Management Essay Writing Assignment Help exhausted you? Hire us for A++.
tutorsglobe.com disadvantages of capitalist economy assignment help-homework help by online capitalist economy tutors
Theory and lecture notes of ultinomial, Poisson and HyperGeometric all along with the key concepts of Multinomial Probabilities, Poisson Probabilities and Hypergeometric Probabilities. Tutorsglobe offers homework help, assignment help and tutor’s assistance on Multinomial, Poisson and HyperGeometric.
tutorsglobe.com externalities-open access and public goods assignment help-homework help by online intermediate microeconomics tutors
tutorsglobe.com evolution of parasitic association assignment help-homework help by online helminthology tutors
Manufacture of Chlor-alkali tutorial all along with the key concepts of Chlor-Alkali Processes, chlor-alkali electrolytic processes, mercury cell method, diaphragm electrolytic cell method and Manufacture of sodium carbonate
tutorsglobe.com imbibition assignment help-homework help by online absorption and movement tutors
Tutorsglobe offers computer science homework help, computer science assignment help, computer science online tutoring assistance with live online qualified tutors.
We all are well known with the everyday routine; we rinse our hair and after that automatically pick up the hairdryer, flip the switch and there are the warm air to dry fast and let us get on with the day.
One dimensional progressive waves tutorial all along with the key concepts of Waves on a Stretched String, Waves in a Fluid, Sound Waves in a Gas, Waves in a Uniform Rod, Wave motion and impedance, Impedance Offered by Strings-Transverse Waves
www.tutorsglobe.com offers development of software design model homework help, assignment help, case study, writing homework help, online tutoring assistance by computer science tutors.
1948730
Questions Asked
3689
Tutors
1469219
Questions Answered
Start Excelling in your courses, Ask an Expert and get answers for your homework and assignments!!