Undergrad Problem getting the coefficients of a non-homogeneous PDE using the Fourier method

Click For Summary
SUMMARY

The discussion focuses on solving a non-homogeneous partial differential equation (PDE) using the Fourier method. The equation under consideration is \(u_t = ku_{xx} + h\) with specified initial and boundary conditions. The solution is split into two parts: a steady-state solution \(U(x)\) and a transient solution \(v(x,t)\). The coefficients \(a_n\) for the Fourier series expansion are derived but do not match the expected solution, particularly for \(n=1\), due to the orthogonality of sine and cosine functions.

PREREQUISITES
  • Understanding of non-homogeneous PDEs
  • Familiarity with Fourier series and orthogonality principles
  • Knowledge of boundary value problems
  • Proficiency in separation of variables technique
NEXT STEPS
  • Study the derivation of Fourier series coefficients in non-homogeneous PDEs
  • Learn about the orthogonality of sine and cosine functions in Fourier analysis
  • Explore the method of separation of variables for solving PDEs
  • Investigate the implications of boundary conditions on the solution of PDEs
USEFUL FOR

Mathematicians, physicists, and engineers dealing with heat transfer, wave equations, or any field requiring the solution of non-homogeneous PDEs using Fourier methods.

Phys pilot
Messages
28
Reaction score
0
Hello, I have to solve this problem:

$$u_t=ku_{xx}+h \; \;\; \; \; 0<x<1 \; \; \,\; t>0$$
$$u(x,0)=u_0(1-\cos{\pi x}) \; \;\; \; \; 0\leq x \leq 1$$
$$u(0,t)=0 \; \;\; \; \; u(1,t)=2u_0 \; \;\; \; \; t\geq0$$

So I know that I can split the solution in two (I don't know the reason. I would appreciate a short explanation because we introduce a new function depending just on x and because the sum of the both solutions are also a solution):
$$u(x,t)=v(x,t)+U(x)$$

This give us two different problems.
$$kU''(x)+h=0$$
and
$$v_t=kv_{xx}$$
$$\text{I.C}\begin{cases} v(x,0)=u_0(1-\cos{\pi x})-U(x)\end{cases}$$
$$\text{B.C}\begin{cases} v(0,t)=0 \\ v(1,t)=0 \end{cases}$$

First I can sol the U(x) which is giving me:
$$kU''(x)+h=0$$
$$U(x)=-\frac{hx^2}{2k}+\left(2u_0+\frac{h}{2k} \right)x$$

Now I solve the PDE with the new boundary conditions using the separation of variables and I get:
$$X_n(x)=B_n\sin{n\pi x}$$
$$T_n(t)=C_ne^{-kn^2\pi^2 t}$$
$$v(x,t)=\sum_{n=1}^\infty{=a_ne^{-kn^2\pi^2 t}\sin{n\pi x}}$$

If I apply the initial condition:

$$v(x,0)=\sum_{n=1}^\infty{=a_n\sin{n\pi x}}=u_0(1-cos(\pi x))-U(x)$$

So finally the coefficients should be given by the expression:

$$a_n=2\int_0^1 (u_0(1-cos(\pi x))-U(x))\sin(n \pi x)dx$$

But I'm not gettin the supposed solution which is:
$$a_n=\frac{2u_0}{n\pi}[1+(-1)^n]+\frac{2u_0n}{(n^2-1)\pi}[1-(-1)^n]+\frac{2h}{k\pi^3 n^3}[(-1)^n-1]$$

Actually the supposed solution to the problem is:

$$u(x,t)=-\frac{hx^2}{2k}+\left(2u_0+\frac{h}{2k}\right)x-\frac{4h}{k\pi}e^{-k\pi^2 t}\sin{\pi x}+ \sum_{n=2}^\infty{a_n e^{-kn^2\pi^2 t}\sin{n\pi x}} $$
As you can see the expansion series starts at 2 instead of one so how am I supposed to get the term for n=1 because you can not get it from a_n
 
Last edited:
Physics news on Phys.org
Recall that \cos \pi x and \sin \pi x are orthogonal, so there is no n=1 term in the expansion of \cos \pi x, although there may be from u_0 - U(x).
 

Similar threads

  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 36 ·
2
Replies
36
Views
5K
  • · Replies 4 ·
Replies
4
Views
3K
Replies
6
Views
2K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 11 ·
Replies
11
Views
4K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K