Oplossingen warmtevergelijking

Moderators: dirkwb, Xilvo

Reageer
Gebruikersavatar
Berichten: 524

Oplossingen warmtevergelijking

Hallo allemaal,

De warmtevergelijking
\(u_t = u_{xx}\)
heeft voor elke
\(a \in \mathbb{R}\)
oplossingen van de vorm:
\(u(x,t) = \frac{1}{t^a}f(y)\)
, met
\(y = \frac{x}{\sqrt{t}}, t > 0\)
.

Laat zien dat
\(f(y)\)
voldoet aan:
\(f''(y) + \frac{1}{2}yf'(y) + af(y) = 0\)
.[/b]

---

Volgens mij is dit simpel aan te tonen, maar ik weet niet goed hoe. Ik dacht eraan dat
\(f(y) = u(x,t)t^a\)
, maar dat lijkt me niet de bedoeling.

Alvast bedankt!

- Fruitschaal.

Gebruikersavatar
Berichten: 10.179

Re: Oplossingen warmtevergelijking

Vul je functie u(x, t) gewoon eens in in je PDE?
Zoek je graag naar het meest interessante wetenschapsnieuws? Wij zoeken nog een vrijwilliger voor ons nieuwspostteam.

Berichten: 4.246

Re: Oplossingen warmtevergelijking

Is ' afgeleide naar de tijd? En met Driekse's hint kom je er trouwens.
Quitters never win and winners never quit.

Gebruikersavatar
Berichten: 524

Re: Oplossingen warmtevergelijking

dirkwb schreef: di 09 apr 2013, 18:56
Is ' afgeleide naar de tijd? En met Driekse's hint kom je er trouwens.
Gewoon de afgeleide van
\(f\)
, want die heeft als enige variabele
\(y\)
. Die variabele is wel gelijk aan een uitdrukking met twee variabelen, maar dat maakt hier volgens mij niet. Voor
\(u_t\)
differentieer je f(y) naar t en voor
\(u_{xx}\)
tweemaal naar x.

Ondertussen was ik al op de oplossing gekomen. Ik had wat ik dacht te snel afgeschreven, maar dat leek toch de manier te zijn.
\(u(x,t) = \frac{1}{t^a}f(\frac{x}{\sqrt{t}})\)
\(u_t = \frac{1}{t^a}\left(-\frac{x}{2t\sqrt{t}}f'\left(\frac{x}{\sqrt{t}}\right) - \frac{a}{t}f\left(\frac{x}{\sqrt{t}}\right)\right) = \frac{1}{t^{a+1}}\left(-\frac{1}{2}yf'(y) - af(y)\right)\)
\(u_x = \frac{1}{t^a\sqrt{t}}f'(y)\)
\(u_{xx} = \frac{1}{t^{a+1}}f''(y)\)
Er geldt:
\(u_t = u_{xx} \Rightarrow u_{xx} - u_t = 0\)
Dus:
\(\frac{1}{t^{a+1}}f''(y) - \frac{1}{t^{a+1}}\left(-\frac{1}{2}yf'(y) - af(y)\right) = 0\)
\(\frac{1}{t^{a+1}}\left(f''(y) + \frac{1}{2}yf'(y) + af(y)\right) = 0\)
\(\frac{1}{t^{a+1}} \neq 0\)
, dus:
\(f''(y) + \frac{1}{2}yf'(y) + af(y) = 0\)
En dit moest aangetoond worden.

Bedankt voor het bevestigen :)

-----

Bepaal de oplossing van
\(u_t = u_{xx}\)
die voldoet aan
\(u(0,t) = 0\)
voor
\(t > 0\)
en
\(u(x,0) = 1\)
voor
\(x > 0\)
. Hint: Stel
\(a = 0\)
.[/b]

Als
\(a = 0\)
, dan
\(u(x,t) = f(y)\)
. Dus dat zou betekenen dat als ik
\(f''(y) + \frac{1}{2}yf'(y) + af(y) = 0\)
kan oplossen voor
\(f(y)\)
, ik ook
\(u(x,t)\)
gevonden heb, toch? Dan de beginvoorwaarden invullen en dat moet de oplossing geven.

Als dit het geval is: hoe los ik
\(f''(y) + \frac{1}{2}yf'(y) + af(y) = 0\)
op?

Gebruikersavatar
Berichten: 10.179

Re: Oplossingen warmtevergelijking

Schrijf je vergelijking als
\(f''(y) + \frac{1}{2} (y f(y))' = 0\)
of
\(f'(y) + \frac{1}{2} (y f(y)) = cte\)


Edit: hmm, die a over het hoofd gezien... Even denken dus ;) .
Zoek je graag naar het meest interessante wetenschapsnieuws? Wij zoeken nog een vrijwilliger voor ons nieuwspostteam.

Gebruikersavatar
Berichten: 524

Re: Oplossingen warmtevergelijking

Hoe moet het dan met die a? :P

Gebruikersavatar
Berichten: 524

Re: Oplossingen warmtevergelijking

Drieske schreef: di 09 apr 2013, 20:52
Schrijf je vergelijking als
\(f''(y) + \frac{1}{2} (y f(y))' = 0\)
of
\(f'(y) + \frac{1}{2} (y f(y)) = cte\)
Edit: hmm, die a over het hoofd gezien... Even denken dus ;) .
Oh, die
\(af(y)\)
valt natuurlijk weg, want
\(a = 0\)
.

Dus dan moet ik
\(f''(y) + \frac{1}{2}yf'(y) = 0\)
oplossen.

Gebruikersavatar
Berichten: 10.179

Re: Oplossingen warmtevergelijking

Ik snap eerlijk gezegd echt niet vanwaar die hint "stel a=0" komt eerlijk gezegd. Wat is het nut daarvan? Persoonlijk vind a=1/2 veel makkelijker en logischer :P . Ben je zeker dat die hint klopt?
Zoek je graag naar het meest interessante wetenschapsnieuws? Wij zoeken nog een vrijwilliger voor ons nieuwspostteam.

Gebruikersavatar
Berichten: 524

Re: Oplossingen warmtevergelijking

Jep, die klopt.

Het is toch juist voordelig als die hele f(y) wegvalt daardoor?

Moet ik nu stellen dat g(y) = f'(y) zodat de vergelijking g'(y) = -(1/2)y*g(y) wordt?

Gebruikersavatar
Berichten: 10.179

Re: Oplossingen warmtevergelijking

Ja natuurlijk, dom van mij. Ja, dat werkt (komt btw quasi hetzelfde uit met a=1/2). Kun je die vgl oplossen?
Zoek je graag naar het meest interessante wetenschapsnieuws? Wij zoeken nog een vrijwilliger voor ons nieuwspostteam.

Gebruikersavatar
Berichten: 524

Re: Oplossingen warmtevergelijking

Ik dacht dat ik hier al geantwoord had.
\(g'(y) = -(1/2)*y*g(y)\)
\(g(y) = Ce^{-\frac{1}{4}y^2}\)
Dus
\(f'(y) = Ce^{-\frac{1}{4}y^2}\)
Dat houdt in dat:
\(u(x,t) = \int f'(y) dy = \int Ce^{-\frac{1}{4}y^2} dy = C \int e^{-\frac{x}{4t}} d\left(\frac{x}{\sqrt{t}}\right)\)
Als tip heb ik meegekregen:
\(\int^{\infty}_{-\infty} e^{-\alpha x^2} dx = \sqrt{\frac{\pi}{\alpha}}\)
,
\(\alpha > 0\)
In dit geval is het duidelijk dat
\(\alpha = \frac{1}{4}\)
, dus
\(\int_{-\infty}^{\infty} f'(y) dy = C\sqrt{4\pi}\)
Het probleem is dat ik niet weet hoe ik:
\(u(x,t) = C \int e^{-\frac{x}{4t}} d\left(\frac{x}{\sqrt{t}}\right)\)
moet oplossen.

Ik dacht wel aan iets als:
\(u(0,t) = 0 \Rightarrow C \int e^0 d(0) = 0\)
. Geen flauw idee wat hier staat.
\(u(x,0) = 1 \Rightarrow C \int e^{-\infty} d(\infty) = 1\)
. Wederom geen idee.

Reageer