Solving modified heat equation

  • #1
psie
261
32
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
  • #2
Think about the long-term steady state. Find a function [itex]u_\infty(x)[/itex] which satisfies [tex]u_\infty'' - hu_\infty = 0[/tex] with [itex]u_\infty(0)= 0[/itex] and [itex]u_\infty(\pi) = 1[/itex]. Then [itex]f(x,t) = u(x,t) - u_\infty(x)[/itex] must satisfy [itex]f_t = f_{xx} - hf[/itex] subject to the self-adjoint boundary condition [itex]f(0,t) = f(\pi,t) = 0[/itex] and the initial condition [itex]f(x,0) = -u_\infty(x)[/itex].
 
  • Like
Likes psie
  • #3
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}.$$
 
  • #4
Do you not know how to determine the coefficients [itex]b_n[/itex]?
 
  • #5
pasmith said:
Do you not know how to determine the coefficients [itex]b_n[/itex]?
Right, I don't. I don't see any connection between ##\sinh## and ##\sin## that's useful here.
 
  • #6
If [tex] \sum_{n=1}^\infty b_n \sin nx = f(x), \quad x \in [0, \pi][/tex] then [tex]
b_n = \frac{2}{\pi} \int_0^\pi f(x) \sin nx \,dx.[/tex] This should be derived in any decent textbook on fourier series.

Integrals of the form [itex]\int \sin ax \sinh bx\,dx[/itex] can be done by integratin by parts twice or by expressing everything in terms of exponentials.
 
  • Like
Likes hutchphd, Orodruin and psie

FAQ: Solving modified heat equation

What is a modified heat equation?

A modified heat equation is a variation of the classical heat equation, which may include additional terms or modifications to account for various physical phenomena such as non-linear diffusion, variable conductivity, or external sources and sinks. These modifications can make the equation more complex and reflective of real-world scenarios.

How do you solve a modified heat equation analytically?

Solving a modified heat equation analytically typically involves methods such as separation of variables, integral transforms (like Fourier or Laplace transforms), or similarity solutions. However, due to the complexity introduced by the modifications, finding an exact analytical solution may not always be possible, and approximations or numerical methods might be required.

What numerical methods are commonly used to solve modified heat equations?

Numerical methods commonly used to solve modified heat equations include finite difference methods (FDM), finite element methods (FEM), and finite volume methods (FVM). These methods discretize the equation and solve the resulting system of algebraic equations using computational algorithms, allowing for the approximation of the solution over a defined grid or mesh.

How do boundary and initial conditions affect the solution of a modified heat equation?

Boundary and initial conditions are crucial in determining the solution of a modified heat equation. They provide the necessary constraints to ensure a unique solution. Different types of boundary conditions, such as Dirichlet, Neumann, or Robin conditions, can significantly influence the behavior of the solution. The initial condition defines the state of the system at the beginning of the observation period and sets the stage for how the solution evolves over time.

Can software tools be used to solve modified heat equations, and which ones are recommended?

Yes, software tools can be used to solve modified heat equations. Some of the recommended tools include MATLAB, COMSOL Multiphysics, ANSYS, and Mathematica. These tools offer built-in functions and solvers specifically designed for partial differential equations (PDEs), including heat equations, and provide robust environments for both analytical and numerical solutions.

Back
Top