Hi Daniel,
I'm getting a better grip on what's happening, and it's embarrassing I didn't
notice that before. Ex. 5.1 was solved via a "standard" FV technique, whereas
Exponential and Power Law are enhancements of such "standard" technique.
I apologize for wasting your time with this. I will be more careful next time I
post to the list.
I used both Exponential...() and PowerLaw...() methods, the former faring
slightly better than the latter in the last (the 5th) cell. Upwind...() did not
produce good results, but I'm pretty sure I'm making a mistake somewhere.
Fausto
>________________________________
> From: Daniel Wheeler <[email protected]>
>To: Multiple recipients of list <[email protected]>
>Sent: Monday, March 26, 2012 1:12:23 PM
>Subject: Re: Example 5.1 revisited
>
>
>On Mon, Mar 26, 2012 at 10:42 AM, Fausto Arinos de A. Barbuto
><[email protected]> wrote:
>
>
>>
>>Hi Daniel,
>>
>>
>>More and more I like this neat piece of software called FiPy but less and
>>less I understand
>>it. :-) I moved on to Example 5.1-ii (where the Péclet is > 2 and therefore
>>numerical
>>instabilities should pop up) but FiPy simply gave rat's ass to this all and
>>came up with
>>the most precise, almost-analytical solution I have ever seen for
>>that example.
>
>With the power law convection scheme, it should be stable (and first order
>accurate) through the whole range of Peclet numbers. Which convection scheme
>are you using?
>
>I apologize for flooding the list with such simple questions and examples, but
>before the
>>more complex CFD configurations come the simple ones. Soon I will be
>>tackling transient
>>and 2D problems, once again following V&M's tracks.
>
>It's cool. All feedback is welcome.
>
>I didn't know that there were a wiki where I might publish my findings.
>That's (good) news
>>to me.
>
>I just noticed that you can't create a wiki or blog post at the moment. I'll
>try and fix that.
>
>--
>Daniel Wheeler
>
>_______________________________________________
>fipy mailing list
>[email protected]
>http://www.ctcms.nist.gov/fipy
> [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
>
>
>
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
