Re: Solving a convection diffusion reaction equation

2012-06-07 Thread Matej Svejda
/j7YuD This is again not the solution of \varphi for a given time but the convolution of \varphi with a weighting distribution over time. On Wed, Jun 6, 2012 at 6:01 PM, Daniel Wheeler wrote: > On Tue, Jun 5, 2012 at 12:59 PM, Matej Svejda wrote: >>> OK. My comments still stand. The d

Re: Solving a convection diffusion reaction equation

2012-06-05 Thread Matej Svejda
> OK. My comments still stand. The discrepancies do not have the appearance of > a "bad" solution, but rather that they are not solutions to the same equation > and/or boundary conditions. > > Are you solving the same equation and boundary conditions? How do you know? So the equation I am tryin

Re: Solving a convection diffusion reaction equation

2012-05-31 Thread Matej Svejda
> As I asked on May 21, 2012, at 11:08 AM: > >> It is difficult to be sure, but the green curve (Mathematica solution?) >> seems to have a small slope at alpha=1, where as the blue curve (FiPy >> solution?) looks like it has zero slope. Could there be a difference in the >> boundary conditions?

Re: Solving a convection diffusion reaction equation

2012-05-31 Thread Matej Svejda
uracy? On Wed, May 23, 2012 at 5:56 PM, Jonathan Guyer wrote: > > On May 23, 2012, at 9:28 AM, Matej Svejda wrote: > >> I want for the H(\alpha, t) function to be automatically updated as >> the Variable t changes. But I have a problem: The value of H is >> defined piec

Re: Solving a convection diffusion reaction equation

2012-05-23 Thread Matej Svejda
Jonathan, Thank you for your help. I have a few more questions: 1) I want for the H(\alpha, t) function to be automatically updated as the Variable t changes. But I have a problem: The value of H is defined piecewise, meaning that I have a treshold \epsilon. Depending on whether \alpha <= \epsilo

Re: Solving a convection diffusion reaction equation

2012-05-19 Thread Matej Svejda
First of all: Thank you so much for your help! > I would use .getFaceGrad() here, not .getFaceGradAverage(). I didn't realize > we even had such a function (for five years!), and contrary to its > documentation, I think it has to be 1st order accurate, not 2nd. For what > it's worth, it doesn't

Solving a convection diffusion reaction equation

2012-05-16 Thread Matej Svejda
Hi, I'm having some problems with solving a PDE using FiPy. The equation is of the form: C * \frac{\partial \varphi(\alpha, t)}{\partial t} = \frac{\partial^2 \varphi(\alpha, t)}{\partial \alpha^2} + \frac{\partial}{\partial \alpha} (\frac{\partial H(\alpha, t)}{\partial \alpha} * \varphi(\alpha

Re: Second derivative of function in Source

2012-05-16 Thread Matej Svejda
Thank you very much for your help :) On Mon, May 14, 2012 at 3:47 PM, Jonathan Guyer wrote: > > On May 12, 2012, at 9:44 AM, Matej Svejda wrote: > >> I'm using FiPy to solve a Convection Diffusion Reaction Equation (one >> spacial dimension, plus time). I want my so

Second derivative of function in Source

2012-05-12 Thread Matej Svejda
Hello, I'm using FiPy to solve a Convection Diffusion Reaction Equation (one spacial dimension, plus time). I want my source to: da^2(H(a, t)) * P(a, t) Where da^2(H(a,t)) is the second special derivative of a function H (a known function and NOT the function I'm solving for) and P the function