Hello,
To solve a steady-state diffusion-convection problem such as the one presented
in
pg. 83 of FiPy Manual Release 2.1.3 (examples.convection.exponential1D.mesh1D)
a tridiagonal and a column matrix must be put together. I was wondering if there
isn't a way to display these matrices using some FiPy built-in function. Is
that
possible? I did a quick (and unsuccessful) search in the manual but didn't find
anything specifically related to that.
I ask that because I've recently started to learn FiPy (and I'm very well
impressed
so far) and decided to solve the worked-out examples in Versteeg &
Malalasekera's
great book, "An Introduction to Computational Fluid Dynamics". I know, the
examples are very simple but, as I said, I've just started with FiPy and the
idea is
to move on to more complex situations later. I successfully solved the first
two pure
diffusion examples (4.1 and 4.2) and decided to tackle a diffusion-convection
one,
namely Example 5.1 in page 137. The results I found with FiPy are very close
to the
figures presented in the book, but aren't quite the same. That is the reason
why I
want to inspect the matrices -- I want to compare them with the ones presented
in
the book. I'm sure this is not a FiPy problem, it's my code. Hopefully the
matrices
will tell me where I'm goofing.
Thanks,
Fausto
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
