Re: Problem in solving Poisson equation with 1D Cylindrical mesh

2010-05-20 Thread Daniel Wheeler
the same result... Tomorrow I'll try a completely clean install of python, fipy and all packages on another pc Thanks Date: Tue, 18 May 2010 18:46:05 -0400 From: daniel.wheel...@gmail.com To: fipy@nist.gov Subject: Re: Problem in solving Poisson equation with 1D Cylindrical mesh

RE: Problem in solving Poisson equation with 1D Cylindrical mesh

2010-05-19 Thread Eduard Manley
@nist.gov Subject: Re: Problem in solving Poisson equation with 1D Cylindrical mesh Make that branches/version-2_1 rather than trunk/. On Tue, May 18, 2010 at 5:28 PM, Daniel Wheeler daniel.wheel...@gmail.com wrote: It might be easier to do a fresh checkout of trunk rather than copying

Re: Problem in solving Poisson equation with 1D Cylindrical mesh

2010-05-19 Thread Daniel Wheeler
...@gmail.com To: fipy@nist.gov Subject: Re: Problem in solving Poisson equation with 1D Cylindrical mesh Make that branches/version-2_1 rather than trunk/. On Tue, May 18, 2010 at 5:28 PM, Daniel Wheeler daniel.wheel...@gmail.com wrote: It might be easier to do a fresh checkout of trunk rather

Re: Problem in solving Poisson equation with 1D Cylindrical mesh

2010-05-18 Thread Daniel Wheeler
@nist.gov Subject: Re: Problem in solving Poisson equation with 1D Cylindrical mesh An older version of fipy (tr...@3480, for example) gives, In [1]: from fipy import * m = CylindricalGrid1D(dx=(1.,2.)) /users/wd15/Documents/python/fipy/tr...@3480/fipy/viewers/gnuplotViewer/gnuplot1DViewer.py:76

Re: Problem in solving Poisson equation with 1D Cylindrical mesh

2010-05-17 Thread Daniel Wheeler
is wrong. I am doing some error? Thanks Date: Thu, 13 May 2010 16:51:39 -0400 From: daniel.wheel...@gmail.com To: fipy@nist.gov Subject: Re: Problem in solving Poisson equation with 1D Cylindrical mesh With any luck this has now been dealt with. Latest versions of branches/version-2_1

RE: Problem in solving Poisson equation with 1D Cylindrical mesh

2010-05-11 Thread Eduard Manley
Date: Tue, 11 May 2010 10:45:24 -0400 From: daniel.wheel...@gmail.com To: fipy@nist.gov Subject: Re: Problem in solving Poisson equation with 1D Cylindrical mesh Can you put together the simplest script possible that demonstrates the problem and I'll try and debug it? Thanks

Re: Problem in solving Poisson equation with 1D Cylindrical mesh

2010-05-06 Thread Jonathan Guyer
On May 5, 2010, at 9:53 PM, Eduard Manley wrote: I'm using a logarithmic discretization (first dr= 5e-04 and next dr increasing as 1.05)(with internal radius= 0.00125, external radius=0.03) and, for some unknown reason this create problems and wrong result with cylindrical 1D mesh. I

Re: Problem in solving Poisson equation with 1D Cylindrical mesh

2010-05-06 Thread Daniel Wheeler
Edward, This may be to do with having a very small volume (or area or line or point) for the inner most element of the domain. It should be the same whether one is using a 1D or 2D mesh. Since you are getting differences in the 1D and 2D case, it should be relatively easy to debug and figure out

RE: Problem in solving Poisson equation with 1D Cylindrical mesh

2010-05-06 Thread Eduard Manley
Subject: Re: Problem in solving Poisson equation with 1D Cylindrical mesh Edward, This may be to do with having a very small volume (or area or line or point) for the inner most element of the domain. It should be the same whether one is using a 1D or 2D mesh. Since you are getting differences

Problem in solving Poisson equation with 1D Cylindrical mesh

2010-05-05 Thread Eduard Manley
Hi, I am trying to solve an eq of this type: A(d phi/d t) = div (D grad phi) + q in cylindrical coordinates. phi is the temperature. (A= 350, D= 05, q a source(joule heating (which is as an hyperbola)) (BCs: Left - fixedflux =0, Right - FixedValue = 20.) With 2D mesh

Problem in solving Poisson equation with 1D Cylindrical mesh

2010-05-05 Thread Eduard Manley
Problem partially solved: I'm using a logarithmic discretization (first dr= 5e-04 and next dr increasing as 1.05)(with internal radius= 0.00125, external radius=0.03) and, for some unknown reason this create problems and wrong result with cylindrical 1D mesh. I tried using a uniform