Suggested Problems Set #8 ECH 6847, Fall 2016 Completion Date: Wednesday December 14th
Suggested Problems Set #8 ECH 6847, Fall 2016 Completion Date: Wednesday December 14th
Suggested Problems Set #8 ECH 6847, Fall 2016 Completion Date: Wednesday December 14th
Take logarithms of the equation and write a model that is linear in the parameters
ln(k0 ) and E/R. Summarize the data and model with the linear algebra problem
A x = b,
in which x contains the parameters of the least squares problem,
ln(k0 )
x =
E/R
b)
c)
Euler scheme. Both methods are first order accurate, as both can be developed from firstorder Taylor series approximations. The explicit method estimates the solution for x(t + t)
using known information at the current time t. The implicit method uses the value of f at the
new point, where the value is not known, requiring the solution to a potentially complicated
algebraic equation at each step.
So why might one use the implict method instead of the simpler explicit method? The
implict Euler method is unconditionally stable, whereas the explict method is only conditionally stable. To demonstrate, compute x at t = 14.7 using t = 2.1 using both methods
and compare to the exact solution when f (x, t) = x; plot your results as a function of
time. You will see that the explicit method generates a solution that alternates about the
real solution and that the implicit method does not. The first is unstable, the second is
stable.
A criteria for stability is that |G| < 1 for these methods, where
x(t + t) = Gx(t).
Identify the range of time steps at which the explicit method is stable and unstable. Also
demonstrate that the implicit method is stable, regardless of the time step.
Problem #2: Temperature measurement by analyis of an extended surface
A long rod within an oven having an air temperature of T = 400oC is pressed firmly onto
the surface of a billet. Thermocouples imbedded in the rod at locations 25 and 100 mm from
the billet register temperatures of 325 and 375oC, respectively. What is the temperature of
the billet?
Note that the temperature in the rod is governed by the differential equation,
d2 T (z)
4h
(T (z) T ) = 0
2
dz
kD
where z is the coordinate along the rod, z = 0 at the billet, and the temperature approaches
4h
, can be determined from the provided
400o C as z goes to infinity. The coefficient, kD
information.