- #1

- 212

- 0

I'm using Green's Functions for heat conduction problems, and I'm trying to solve the following integral:

http://img28.imageshack.us/img28/4923/026307b169b04faa8364086.png [Broken]

Where:

http://img820.imageshack.us/img820/3742/6332938c445f4b9e8da8ba5.png [Broken]

N/A

I obliviously broke this up into two integrals which I'll call Integral #1 and Integral #2.

For Integral #1 I used some algebraic simplification and a u-substitution (u=4α(t-τ)) to get the following:

http://img600.imageshack.us/img600/2892/121092a52813459b8280e05.png [Broken]

Using Mathmatica, I got:

http://img856.imageshack.us/img856/719/92de0714ca1a4ec8abd7c29.png [Broken]

For which I get the following after putting in the limits of the integral and evaluate some limits:

http://img195.imageshack.us/img195/6397/26da0028e4824870aa77365.png [Broken]

For Integral #2 I again did some algebraic simplification and a u-substitution (u=α(t-τ)) to get:

http://img838.imageshack.us/img838/7506/fb9afa29d9e44617861b1a6.png [Broken]

Using Mathmatica again:

http://img194.imageshack.us/img194/2729/25901d5d37a644cbab647ec.png [Broken]

Where a=(x/2) and b=(h/k).

After doing some algebra cleanup, substituting back in the values of a and b, and applying the limits of integration:

http://img571.imageshack.us/img571/7754/c206229ff433444c9d7b31d.png [Broken]

After combining the solutions for Integrals #1 and #2 with the constants that were pulled at the start:

http://img191.imageshack.us/img191/9814/b5fd0b3ed9294586bb55efc.png [Broken]

Where f

Now, when I did dimensional analysis for x=t=0, I got units of temperature; however, when I tried to plot my solution in Matlab, I end up getting garbage. I've doubled checked my code in Matlab to make sure the solution is typed in correctly (and I am in the process of triple checking) and I have double checked my integration. So far, I can't find a mistake.

I'd prefer to do the integrals without Mathmatica, but I haven't been able to find an appropriate table of integrals. I THINK I can do Integral #1 entirely by hand using the definition of the error function, but that still leaves Integral #2, which I haven't been able to crack.

Any thoughts/help/advice?

## Homework Statement

http://img28.imageshack.us/img28/4923/026307b169b04faa8364086.png [Broken]

Where:

http://img820.imageshack.us/img820/3742/6332938c445f4b9e8da8ba5.png [Broken]

## Homework Equations

N/A

## The Attempt at a Solution

I obliviously broke this up into two integrals which I'll call Integral #1 and Integral #2.

For Integral #1 I used some algebraic simplification and a u-substitution (u=4α(t-τ)) to get the following:

http://img600.imageshack.us/img600/2892/121092a52813459b8280e05.png [Broken]

Using Mathmatica, I got:

http://img856.imageshack.us/img856/719/92de0714ca1a4ec8abd7c29.png [Broken]

For which I get the following after putting in the limits of the integral and evaluate some limits:

http://img195.imageshack.us/img195/6397/26da0028e4824870aa77365.png [Broken]

For Integral #2 I again did some algebraic simplification and a u-substitution (u=α(t-τ)) to get:

http://img838.imageshack.us/img838/7506/fb9afa29d9e44617861b1a6.png [Broken]

Using Mathmatica again:

http://img194.imageshack.us/img194/2729/25901d5d37a644cbab647ec.png [Broken]

Where a=(x/2) and b=(h/k).

After doing some algebra cleanup, substituting back in the values of a and b, and applying the limits of integration:

http://img571.imageshack.us/img571/7754/c206229ff433444c9d7b31d.png [Broken]

After combining the solutions for Integrals #1 and #2 with the constants that were pulled at the start:

http://img191.imageshack.us/img191/9814/b5fd0b3ed9294586bb55efc.png [Broken]

Where f

_{i}(τ)=hT_{∞}Now, when I did dimensional analysis for x=t=0, I got units of temperature; however, when I tried to plot my solution in Matlab, I end up getting garbage. I've doubled checked my code in Matlab to make sure the solution is typed in correctly (and I am in the process of triple checking) and I have double checked my integration. So far, I can't find a mistake.

I'd prefer to do the integrals without Mathmatica, but I haven't been able to find an appropriate table of integrals. I THINK I can do Integral #1 entirely by hand using the definition of the error function, but that still leaves Integral #2, which I haven't been able to crack.

Any thoughts/help/advice?

Last edited by a moderator: