1
$\begingroup$

I am using the FFT method to solve the two dimensional PDE on a $2\pi\times 2\pi$ square. I have boundary conditions that $u(x,0,t)=0$ and $u(x,2\pi,t)=0$. How do I go about setting these boundaries? I can't just expand in different basis functions because of the way my fft function works.

Thanks in advance.

$\endgroup$
1
  • $\begingroup$ Sounds like a job for the discrete sine transform. Any reference on FFT should include something on how to get a canned FFT routine to do DST. Alternatively, you might want to look at FFTW. $\endgroup$ Oct 29, 2010 at 15:04

1 Answer 1

1
$\begingroup$

If you don't mind wasting some resources, extend the $y$ coordinate to the interval $[0,4\pi]$. Then only look at solutions that are odd (anti-symmetric) under the flip $y\mapsto 4\pi-y$. In particular, if your PDE has source terms, they have to be anti-symmetrically mirrored from the $[0,2\pi]$ to the $[2\pi,4\pi]$ half of the extended interval. When restricted to $y\in[0,2\pi]$ a periodic solution on the extended interval will satisfy the desired boundary conditions.

$\endgroup$
1
  • $\begingroup$ This is essentially what a discrete sine transform does, only with less waste of space due to the exploitation of the antisymmetry. $\endgroup$ Oct 31, 2010 at 10:32

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.