RE: Applying Bounding Limits to PDE Solution Variables

2016-08-02 Thread Campbell, Ian
Hi Daniel, Thank you very much for your prompt reply. I'm manually trying to implement the bound-checking feature immediately following a sweep, as you've suggested. As an initial step towards this, I need to manually compute my own residual and compare to that computed via '.sweep' in FiPy.

Re: Applying Bounding Limits to PDE Solution Variables

2016-08-02 Thread Daniel Wheeler
Hi Ian, Sorry, there isn't any way to do this as an implicit command in FiPy, whether as a simple reset during sweeping or as a fully implicit constraint in the equations. Of course, you can always try these techniques explicitly in the Python code or in the design of the equations. I agree that I

Applying Bounding Limits to PDE Solution Variables

2016-08-02 Thread Campbell, Ian
Hi All, I would like to prevent the CellVariables I am solving for from ever reaching physically-meaningless and incorrect values, such as concentrations becoming negative. Is there a way to do this? The most natural place I could think of it being is an argument in CellVariable declaration,

Re: reaction-advection-diffusion; problems making the solution converge

2016-08-02 Thread Daniel Wheeler
Hi Nils, Could you possible post your code so I can play around with it? I can't think offhand what could be the issue or exactly what the error is that you're calculating, it seems like the normalized residual or something. This is a linear equation so I'm not entirely sure why the non-linear res

Re: error in gradient even with orthogonal grid

2016-08-02 Thread James Pringle
Thank you -- since my problem solves for a stream function, this was a very helpful tip, and has made my solutions appear much better. Jamie On Mon, Aug 1, 2016 at 4:49 PM, Daniel Wheeler wrote: > Hi Jamie, > > Sorry for the delayed reply. > > On Thu, Jul 21, 2016 at 11:31 AM, James Pringle wr