Day 05

Tuesday, Feb 14, 2017

Quiz

Equilibrium temperature profile in a pipe

thick-wall-seamless-pipes-250x250.jpg

Laplacian in cylindrical coordinates is

2u = (1)/(r)()/(r)(r(u)/(r)) + (1)/(r2)(2u)/(θ2) + (2u)/(z2)

2017_02_14/20170214_131837.jpg

Actual plots of log(r ⁄ b) ⁄ log(a ⁄ b):

thick_pipe_plots.png

Ch 2 Method of separation of variables, cont'd

2.3 1D heat equation with prescribed boundary values, cont'd

From last time:

2017_02_09/20170209_135456.jpg 2017_02_14/20170214_131901.jpg 2017_02_14/20170214_131907.jpg 2017_02_14/20170214_131921.jpg 2017_02_14/20170214_135432.jpg 2017_02_14/20170214_135443.jpg 2017_02_14/20170214_135509.jpg

Example IBVP

Coefficients we computed:

Bn_for_first_heat_solution.png

Plotting solution with Python

(From 0 to 60mph in 20 minutes!)

numpy arrays, broadcasting, outer, linspace, arange, sum, base-0 counting, [1:]

functions (def), for loops

pylab, plot, xlabel, ylabel

Illustration of solution: Jupyter Notebook, HTML

Approximation by lowest-λ solution