Consider the Initial Value Problem du/dt = 50/u-50u, 0...
Fantastic news! We've Found the answer you've been seeking!
Question:
Transcribed Image Text:
Consider the Initial Value Problem du/dt = 50/u-50u, 0<t<1, u(0) = √2. We will use Implicit Euler and Implicit Trapezoid to solve this nonlin- ear problem. (a) First, write out each of the schemes for this problem (the equation relating ui+1 to uż). Formulate each as an equation of the form g(i+1, Ui) = 0, a root finding problem. (b) In this particular case, one could solve the root finding problems analytically. However, we will not do so. Instead, we will use Newton's method. As always, this is done by making an initial ,(0) guess for what ₁+1 might be, call this u, then generating a (k) sequence of successive guesses ut based on the standard New- ton's method formula. For each of the two schemes, write out ,,(k+1) the equation for generating the next guess ut from the prior (k) guess u via Newton's method. (c) Now, implement the above in code, and solve the system with each method using N = 10 equal width subdivisions of the inter- val 0 < t < 1. For your root finding, use a stopping tolerance of € = 10-6. For your initial Newton's method guess u at each ,(0) timestep, use the solution you found for the prior timestep, ₁. Plot the two solutions. (d) The exact solution to this problem is u(t)=√1+e-100€. Given this and your numerical solutions above, plot the error E₁ |u(ti) - ui as a function of t for each method, and comment on what you find. = Consider the Initial Value Problem du/dt = 50/u-50u, 0<t<1, u(0) = √2. We will use Implicit Euler and Implicit Trapezoid to solve this nonlin- ear problem. (a) First, write out each of the schemes for this problem (the equation relating ui+1 to uż). Formulate each as an equation of the form g(i+1, Ui) = 0, a root finding problem. (b) In this particular case, one could solve the root finding problems analytically. However, we will not do so. Instead, we will use Newton's method. As always, this is done by making an initial ,(0) guess for what ₁+1 might be, call this u, then generating a (k) sequence of successive guesses ut based on the standard New- ton's method formula. For each of the two schemes, write out ,,(k+1) the equation for generating the next guess ut from the prior (k) guess u via Newton's method. (c) Now, implement the above in code, and solve the system with each method using N = 10 equal width subdivisions of the inter- val 0 < t < 1. For your root finding, use a stopping tolerance of € = 10-6. For your initial Newton's method guess u at each ,(0) timestep, use the solution you found for the prior timestep, ₁. Plot the two solutions. (d) The exact solution to this problem is u(t)=√1+e-100€. Given this and your numerical solutions above, plot the error E₁ |u(ti) - ui as a function of t for each method, and comment on what you find. =
Expert Answer:
Related Book For
Calculus For Scientists And Engineers Early Transcendentals
ISBN: 9780321849212
1st Edition
Authors: William L Briggs, Bernard Gillett, Bill L Briggs, Lyle Cochran
Posted Date:
Students also viewed these accounting questions
-
The statement of cash flows for Picture Perfect Photography reported the following for the year ended December 31, 2014: Cash flows from financing activities: Dividends (declared and) paid $(8,000)...
-
Mr. Johnson is worried about his wife because she had a heart attack and has been in the hospital for over a month now. Over stimulation of the CRH --> ACTH--> -> Cortisol pathway by intense chronic...
-
A cardioid can be described by the polar cquation r = 2+2 cos 0 Determine the gradient of this curve when 0 = T/2 making clear your workings. 5
-
Express the following sum in closed form. 2 k 2 (4 + 3 - 4 ) = n k 1
-
A river flowing steadily at a rate of 240 m3/s is considered for hydroelectric power generation. It is determined that a dam can be built to collect water and release it from an elevation difference...
-
Is Bi-Economy entitled to consequential damages for the destruction of its business?
-
Which corporate characteristic is a disadvantage? a. Double taxation b. Mutual agency c. Limited liability d. None of the above
-
A major focus of the trial in this case was BDO Seidmans consideration of, and response to, the red flags apparent during the 1995 HMI audit. Define or describe the phrase red flags. Explain the...
-
Find the derivative of f(x)=6e5+e-x f'(x) = 9.
-
Hastings Company bases its variable overhead performance report on the actual direct labour- hours of the period. Data concerning the most recent year, which ended on December 31, are as follows:...
-
The chief cashier of a bank must manage the cash holding so as to gain interest by investing excess cash where possible, but at the same time keep sufficient cash on hand to meet the bank's needs. A...
-
Describe the four types of audit opinions.
-
Why would studying and evaluating each of the client's controls not necessarily be the most efficient way of determining the scope of related substantive tests?
-
Describe the narrowly focused matters that the auditor must consider at the account balance level when deciding on the scope of the audit work necessary to provide reasonable assurance of detecting...
-
A survey of 32 medical doctors shows the number of days that they worked in a recent year. Construct a frequency distribution for the data using six classes. Summarize the results. 240 232 280 261...
-
What four steps should the auditor take if he or she determines that a difference between amounts on the client's accounting records and audited amounts is fraud that could be material to the...
-
Veto Company bought equipment on January 1 for $99,000. The expected life is 10 years, and the residual value is $11,000. Based on three acceptable depreciation methods, the company computes annual...
-
A container holds 2.0 mol of gas. The total average kinetic energy of the gas molecules in the container is equal to the kinetic energy of an 8.0 10-3-kg bullet with a speed of 770 m/s. What is the...
-
Given the following vector fields and oriented curves C, evaluate c F T ds. F = (x, y) on the curve r(t) = (1, 312), for 1 t2 (x + y)/2
-
Give the antiderivatives of x p . For what values of p does your answer apply?
-
Determine whether the following vector fields are conservative on the specified region. If so, determine a potential function. Let R* and D* be open regions of R 2 and R 3 , respectively, that do not...
-
Steam flows steadily and isentropically through a nozzle. At an upstream section where the speed is negligible, the temperature and pressure are \(450^{\circ} \mathrm{C}\) and \(6 \mathrm{MPa}\)...
-
Space debris impact is a real concern for spacecraft. If a piece of space debris were to create a hole of \(0.001 \mathrm{in} .^{2}\) area in the hull of the International Space Station (ISS), at...
-
Oxygen discharges from a tank through a convergent nozzle. The temperature and velocity in the jet are \(-20^{\circ} \mathrm{C}\) and \(270 \mathrm{~m} / \mathrm{s}\), respectively. What is the...
Study smarter with the SolutionInn App