Online Maths Assignment Help
Online Maths Assignment Help
Problems
Problem 1 : Modified inner products for column vectors
Consider the inner product (x, y) = x∗By from class (lecture 5.5
notes), where the vectors are in CN and B is an N × N Hermitian
positive-definite matrix.
(a) Show that this inner product satisfies the required properties
of inner products from class: (x, y) = (y, x), (x, x) > 0 except
for x = 0. (Linearity (x, αy+z) = α(x, y)+(x, z) is obvious from
linearity the of matrix operations; you need not show it.)
L=1
M = 100
D = diff1(M)
dx = L / (M+1)
x = dx*0.5:dx:L # sequence of x values from 0.5*dx to <= L in
steps of dx
C = ....something from c(x)...hint: use diagm...
A = -D’ * C * D / dx^2
• The rod is insulated: no heat flows out the sides or through the
ends.
Given these assumptions, you should be able to answer the
following:
(a) “Newton’s law of cooling” says that that the temperature of
an object changes at a rate (K/s) proportional to the
temperature difference with its surroundings. Derive the
equivalent here: show that our assumptions above imply that
for some constant α, for 1 < n < N.
Also give the
(slightly different) equations for n = 1
and n = N.
(b) Write your equation from the previous part in matrix form:
for some matrix A.
(c) Let T (x, t) be the temperature along the rod, and suppose
Tn(t) = T ([n − 0.5]Δx, t) (the temperature at the center of the
n-th piece). Take the limit N → ∞ (with L fixed, so that
and derive a partial differential equation
What is Aˆ? (Don’t worry about the x = 0, L ends until the next part.)
There are other ways to solve this problem of course, that are
also second-order accurate.
(b) In order to approximate Au, we did three things: compute u'
by a center-difference as in class, multiply by cm+0.5 at each
point m + 0.5, then compute the derivative by another center-
difference. The first and last steps are exactly the same center-
difference steps as in class, so they correspond as in class to
multiplying by D and −DT , respectively, where D is the (M +
1) × M matrix
(ii) We find that the dot product is ≈ 4.3 × 10−16, which is zero
up to roundoff errors (your exact value may differ, but
should be of the same order of magnitude).
Problem 3:
(a) The heat capacity equation tells us that
where dQn/dt is the rate
of change of the heat in the n-th piece. The thermal
conductivity equation tells us that dQn/dt, in turn, is equal to
the sum of the rates q at which heat flows from n + 1 and n −
1 into n:
(c) Ignoring the ends for the moment, for all the interior points
we have
which is exactly our familiar center-
difference
approximation for at the point n (x = [n − 0.5]Δx). Hence,
thus
(d) The boundary conditions are The easiest way to see this is to
observe that our heat flow q is really a first derivative, and zero
heat flow at the ends
Means zero derivatives. That is, qn+0.5 =
is really an approximate
derivative: while the flows q0.5 and
qN+0.5 to/from n = 0 and n = N + 1 is zero, and hence q0.5 =
qN+0.5 = 0 ≈
Working backwards, (setting 1 for convenience)
consider
with these boundary conditions and center-difference
approximations. We are given Tn = T([n − 0.5]Δx,t) for n =
1,...,N. First, we compute
for n = 1,...,N − 1 (−DT T using the D above). Unlike the
Dirichlet case in class, we don’t com
for 1 <n<N,
we differentiated,
In the N → ∞ limit, this gives
multiplied
where hΔy is the area of the interface between the two blocks.
Then, to convert into a rate of temperature change, we will
divide by cρhΔxΔy, where hΔxΔy is the volume of the block.
Putting this all together, we obtain:
where the thing in [··· ] is precisely the five-point stencil
approximation for ∇2 from class. Hence, we obtain
where for fun I have put the κ in the middle, which is the right
place if κ is not a constant (you were not required to do this).