WOLFRAM NOTEBOOK

WOLFRAM|DEMONSTRATIONS PROJECT

One-Dimensional Heat Conduction with Temperature-Dependent Conductivity

trial solution order
2
3
4
exact and trial solutions
trial solution:1.32624x-0.326238
2
x
Consider the following nonlinear boundary value problem
d
dz
k(T)
dT
dz
=0
, with
T(z=0)=
T
0
,
T(z=d)=
T
1
, and
0zd
. This equation with the boundary conditions (BCs) describes the steady-state behavior of the temperature of a slab with a temperature-dependent heat conductivity given by
k(T)=
k
0
+α(T-
T
0
)
.
Introducing nondimensional temperature and position using
θ=
T-
T
0
T
1
-
T
0
and
x=
z
d
, the governing equation and BCs become
d
dx
(1+aθ)
dθ
dx
=0
, with
a=α/
k
0
,
θ(x=0)=0
, and
θ(x=1)=1
. Without loss of generality, set
a=1
. Then this problem admits an analytical solution given by
θ(x)=-1+
1+3x
.
We use Galerkin's method to find an approximate solution in the form
θ
N
(x)=x+
N
i=1
A
i
(
i+1
x
-x)
. The unknown coefficients
A
i
of the trial solution are determined using the residual
R(
θ
N
,x)=(1+a
θ
N
(x))
2
d
θ
N
d
2
x
+a
2
d
θ
N
dx
and setting the integrals
1
0
R(
T
N
,x)(
j
x
-x)dx=0
for
j=2,,N+1
. You can vary the order of the trial solution,
N+1
. The Demonstration plots the analytical solution (in gray) as well as the approximate solution (in dashed cyan). The polynomial expression of the trial solution is displayed on the same plot. In addition, the difference between the exact and approximate solutions,
e(x)
, is plotted. This error is seen to decrease as the order of the trial solution,
N+1
, is increased.
Wolfram Cloud

You are using a browser not supported by the Wolfram Cloud

Supported browsers include recent versions of Chrome, Edge, Firefox and Safari.


I understand and wish to continue anyway »

You are using a browser not supported by the Wolfram Cloud. Supported browsers include recent versions of Chrome, Edge, Firefox and Safari.