MATLAB Answers

Green's Function Solution in Matlab

113 views (last 30 days)
yusuf
yusuf on 4 Feb 2014
Answered: Ted Shultz on 16 Oct 2019
I wrote a code for stationary temperature distribution which is;
f = 0.09;
b = 0.0044;
q = 3.73e-9;
L = 1;
Tw = 250;
Tam = 27;
syms c x g
T = 2*c*cosh(x*((f-b*g)/q)^0.5)+g/(f-b*q);
c = solve(subs(T,'x',L/2)==0,c);
z = simplify(int(subs(T,'c',c),x,-L/2,L/2));
g = solve(z==L*(Tw-Tam),g)
T
Now I try to use perturbation method and find T1 temperature distribution. My equation is;
q*T1''(x)-T1(x)*(f-b*g+i*w*p)=T*(f1-b*g1)-g1
And my prof. say that can be solved by Green's function G(x,y), where the G(x,y) is soluton of;
q*G''(x,y)-G(x,y)*(f-b*g+i*w*p)=Diract(x-y)
And
G(L/2)=0 , G(-L/2)=0
if I can find the G(x,y), I will get the solution of T1.
Can anyone help me please?
Thank you for everything
Yusuf

Answers (1)

Ted Shultz
Ted Shultz on 16 Oct 2019
The person who asked this question has some (slow) sample code that will solve Green's function.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by