Hello, I'm trying to numerically solve for the steady-state solution of the following PDE:
\frac{\partial \psi}{\partial t}=-y\frac{\partial \psi}{\partial
x}-2xy(y^2+1)\psi
It can be seen that the analytical steady state solution is:
\psi=exp(-(x^2+1)(y^2+1))
Attached is the program written in FiPy to do this. I've set the initial
function to be a Gaussian and evolved it forward in time. For the first 10
seconds of running the code, very little happens. Then the top left and
lower right portions of the grid develop non-zero values and the solution
diverges.
I'm quite sure that I'm expressing this problem in the right way (ie: using
a convection term) and that the residues and tolerances etc have been set
sufficiently small. Any idea why the solution is diverging instead of
converging to the analytic value?
Thanks for you help.
Regards,
Altan
test2.py
Description: Binary data
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
