kill(all)$

maxima: heat:heat conduction (1D),variable : k,version 1.0,part1(plane wall)
peter.vlasschaert@gmail.com,20/07/2021

variable thermalbconductivity =k
=========================
depends on the conducting medium

model :plane wall
=====

→ ode :

pa:'diff(k*'diff(T,x),x)=0;

pb:k=ko*(1+α*x^2);

depends(T,x);

pc:subst([pb],pa);

≥ final ode 

pc1:ev(pc,diff);

→ bc : (b)oundary (c)onditions (L =length of the wall)

≥ x=0 (begin : wall) ,T=T1

bc1:T(0) = T1;

≥ x=L (end : wall),T=T2

bc2:T(L)=T2;

→ heat flow 'through the wall'

pd:q=-k*A*diff(T,x);

pd1:subst([pb],pd);

solve : model ode plane wall .
=======================
use :pc1

assume(α>0);

m:ode2(pc1,T,x);

listofvars(m);

use :{bc1,bc2}

→ bc1

bc1a:subst([x=0,T=T1],m);

→ bc2

bc2:subst([x=L,T=T2],m);

→ find

p1:solve([bc1a,bc2],[%k1,%k2]);

fnal solution : ode+bc
================

p2:subst(flatten(p1),m);

when α=0

p2a:'limit(p2,α,0)=rhs(ratsimp(limit(p2,α,0)));

heat tranfer to wall:q

first step :differentation p2

p3:'diff(T,x)=diff(rhs(p2),x);

second step: heat flow 'through the wall'

p4:ev(pd1);

p4a:subst([p3],p4);

suppose : A=1, 'surface area plane wall'

A:1;

p4b:ev(p4a);

p4b1:facsum(p4b);