Solving modified heat equation

Click For Summary
The discussion focuses on solving a modified heat equation with non-constant boundary conditions. The transformation using u(x,t) = v(x,t)e^{-ht} simplifies the problem to the standard heat equation, but the boundary conditions for v(x,t) complicate the solution. A steady-state function u_\infty(x) is identified, satisfying the required conditions, and the transformation f(x,t) = e^{-ht}v(x,t) is suggested to revert to the heat equation. The challenge lies in determining the coefficients b_n for the Fourier series representation of the initial condition, with a method involving integration by parts or exponential expressions proposed. The conversation emphasizes the need for a solid understanding of Fourier series to derive the coefficients effectively.
psie
Messages
315
Reaction score
40
Homework Statement
Find a solution of the following problem \begin{align} u_t&= u_{xx} - hu,\qquad &0<x<\pi, \ t>0; \\ u(0,t)&=0,u(\pi,t)=1,\qquad &t>0; \\ u(x,0)&=0,\qquad &0<x<\pi.\end{align} Here ##h>0## is a constant.
Relevant Equations
The heat equation in the form ##u_t= u_{xx}## and its solution ##u(x,t)=\sum_{n=1}^\infty b_n e^{-n^2t}\sin nx##.
In my Fourier analysis book, the author introduces some basic PDE problems and how one can solve these using Fourier series. I know how to solve basic heat equation problems, but the above one is different from the previous problems I've worked in terms of the boundary conditions. Using ##u(x,t)=v(x,t)e^{-ht}## I can transform the equation into the heat equation, i.e. ##v_t= v_{xx}## , however, the boundary conditions become $$v(0,t)=0,\quad v(\pi,t)=e^{ht}.$$ I don't know how to deal with non-constant boundary conditions...any ideas on how to proceed?
 
Physics news on Phys.org
Think about the long-term steady state. Find a function u_\infty(x) which satisfies u_\infty&#039;&#039; - hu_\infty = 0 with u_\infty(0)= 0 and u_\infty(\pi) = 1. Then f(x,t) = u(x,t) - u_\infty(x) must satisfy f_t = f_{xx} - hf subject to the self-adjoint boundary condition f(0,t) = f(\pi,t) = 0 and the initial condition f(x,0) = -u_\infty(x).
 
Ok. I found $$u_{\infty}(x)=\frac{\sinh \sqrt{h}x}{\sinh \sqrt{h}\pi}.$$ I still think I have to use the trick ##f(x,t)=e^{-ht}v(x,t)## to turn ##f_t = f_{xx} - hf## into ##v_t = v_{xx}##. The solution to the latter will be $$v(x,t)=\sum_{n=1}^\infty b_n e^{-n^2t}\sin nx.$$ But I'm stuck at how to solve for ##b_n## in $$f(x,0)=v(x,0)=\sum_{n=1}^\infty b_n \sin nx=-\frac{\sinh \sqrt{h}x}{\sinh \sqrt{h}\pi}.$$
 
Do you not know how to determine the coefficients b_n?
 
pasmith said:
Do you not know how to determine the coefficients b_n?
Right, I don't. I don't see any connection between ##\sinh## and ##\sin## that's useful here.
 
If \sum_{n=1}^\infty b_n \sin nx = f(x), \quad x \in [0, \pi] then <br /> b_n = \frac{2}{\pi} \int_0^\pi f(x) \sin nx \,dx. This should be derived in any decent textbook on fourier series.

Integrals of the form \int \sin ax \sinh bx\,dx can be done by integratin by parts twice or by expressing everything in terms of exponentials.
 
  • Like
Likes hutchphd, Orodruin and psie
First, I tried to show that ##f_n## converges uniformly on ##[0,2\pi]##, which is true since ##f_n \rightarrow 0## for ##n \rightarrow \infty## and ##\sigma_n=\mathrm{sup}\left| \frac{\sin\left(\frac{n^2}{n+\frac 15}x\right)}{n^{x^2-3x+3}} \right| \leq \frac{1}{|n^{x^2-3x+3}|} \leq \frac{1}{n^{\frac 34}}\rightarrow 0##. I can't use neither Leibnitz's test nor Abel's test. For Dirichlet's test I would need to show, that ##\sin\left(\frac{n^2}{n+\frac 15}x \right)## has partialy bounded sums...