Re: cylindrical coordinates

2018-10-10 Thread Martinus WERTS
ike to go 3D, and the symmetry of our system would allow to >> use a 2D cylindrical grid (r,z) - with zero flux at z=0 and z=L (and >> r=0), and either Dirichlet/zero flux at r=R. Looking at the mailing list >> archive and GitHub, it appears that cylindrical coordinates are at the

Re: cylindrical coordinates

2018-10-10 Thread Daniel Wheeler
hive and GitHub, it appears that cylindrical coordinates are at the > moment not working properly (missing factor) Yes, that's correct. Sorry about that. > Is this still the case? > > If so, I will start with a simple 3D Cartesian mesh, and then perhaps > move to more adapted meshes

AW: Gradient in cylindrical coordinates

2018-01-10 Thread Munoz Leyton
recipients of list <fipy@nist.gov> Betreff: Re: Gradient in cylindrical coordinates Also, note that the "leastSquaresGrad" doesn't have this issue: https://gist.github.com/wd15/fc34ccb2e57602fc6f9bea96d8160f4a#file-untitled-ipynb See, https://www.ctcms.nist.gov/fipy/fipy/generated/fip

Re: Gradient in cylindrical coordinates

2018-01-09 Thread Guyer, Jonathan E. Dr. (Fed)
Leyton - A compressible flow example contribution would be most welcome! Please submit a [pull request](https://github.com/usnistgov/fipy/pulls) and we'll work with you to get it integrated and released. [don't hesitate to ask if you need help making the pull request] - Jon > On Jan 5, 2018,

Re: Gradient in cylindrical coordinates

2018-01-08 Thread Daniel Wheeler
Fipy Team, >> >> >> >> First of all I would like to thank you for your amazing work! I love working >> with fipy. >> >> >> >> I have a problem when I calculate the gradient of my variable in cylindrical >> coordinates. You can s

Re: Gradient in cylindrical coordinates

2018-01-08 Thread Daniel Wheeler
at 7:23 PM, Munoz Leyton <leyton.mu...@hirtenberger.com> wrote: > Dear Fipy Team, > > > > First of all I would like to thank you for your amazing work! I love working > with fipy. > > > > I have a problem when I calculate the gradient of my variable in cylindrical &g

Re: 2D Cylindrical Coordinates - Pipe

2012-06-21 Thread Kendall Boniface
The Volume for a 2D mesh is the area of the cell. The Volume for a 1D mesh isthe length of the cell. The CylindricalGrid2D represents a wedge that subtends 1 rad. Ah, I see now. Thank you for the clarification! For what I want to do, I have actually found that the following few lines

Re: 2D Cylindrical Coordinates - Pipe

2012-06-20 Thread Kendall Boniface
On Tue, Jun 19, 2012 at 1:04 PM, Jonathan Guyer gu...@nist.gov wrote: On Jun 19, 2012, at 12:20 PM, Kendall Boniface wrote: I am having a bit of trouble manipulating a 2D cylindrical mesh and was wondering if anyone has any helpful advice? mesh = CylindricalGrid2D(dx=dx, dy=dy, nx=nx,

Re: 2D Cylindrical Coordinates - Pipe

2012-06-20 Thread Jonathan Guyer
On Jun 20, 2012, at 11:20 AM, Kendall Boniface wrote: I'm a bit confused about the getCellVolumes() function. Since I am using the 2D cylindrical grid, I assume it isn't actually giving me volumes so to speak. I've manually played around with some of my numbers and can't seem to figure

2D Cylindrical Coordinates - Pipe

2012-06-19 Thread Kendall Boniface
Hello again, I am having a bit of trouble manipulating a 2D cylindrical mesh and was wondering if anyone has any helpful advice? I want a mesh that is 1.5 meters in the z direction and has an inner radius of 0.0046 m and an outer radius of 0.00635 m (to model the wall of a pipe). This is what I

Re: 2D Cylindrical Coordinates - Pipe

2012-06-19 Thread Jonathan Guyer
On Jun 19, 2012, at 12:20 PM, Kendall Boniface wrote: I am having a bit of trouble manipulating a 2D cylindrical mesh and was wondering if anyone has any helpful advice? mesh = CylindricalGrid2D(dx=dx, dy=dy, nx=nx, ny=ny) + ((0.0046,),) print mesh.getCellCenters() I want the z axis

Re: 3D Transient Conduction - Cylindrical Coordinates

2012-06-05 Thread Kendall Boniface
I pasted the gmsh commands into a file and ran gmsh on it. It threw a weird error ** On entry to DGESVD parameter number 6 had an illegal value What error did you get? Did you adapt the gmsh commands from somewhere? Maybe we could work from that and check that those commands work

Re: ... cylindrical coordinates ...

2007-01-22 Thread Daniel Wheeler
On Jan 19, 2007, at 10:51 PM, Matt Koch wrote: Rejected message: sent to fipy@nist.gov by [EMAIL PROTECTED] follows. Reason for rejection: sender not subscribed. --- Hi Jonathan, I have to deal with both a

Re: ... cylindrical coordinates ...

2007-01-18 Thread Daniel Wheeler
to do it in the case of cylindrical coordinates it to avoid making changes to the discretization of the equations. I believe (needs to be confirmed) that changing the cell volumes and face areas does in fact allow you to use cartesian discretization with cylindrical coordinates. As far as my

Re: ... cylindrical coordinates ...

2007-01-18 Thread Matt Koch
are talking about a 3D model then. No need to talk about cylindrical coordinates then, simply model an arbitrary shape (such as a cylinder) in rectangular 3D space? The question then is, how well will a mesh approximate curved surfaces? Plus, no matter how thin the slice in 3D, a computational

Re: ... cylindrical coordinates ...

2007-01-18 Thread Jonathan Guyer
theta) of the slice? Exactly. We really are talking about a 3D model then. No need to talk about cylindrical coordinates then, simply model an arbitrary shape (such as a cylinder) in rectangular 3D space? Not exactly. The slice of pie will only be one cell deep. It's only a sliver of pie

Re: ... cylindrical coordinates ...

2007-01-17 Thread Daniel Wheeler
, but it is a little too advanced for me. This exchange discusses wedge shaped meshes that are a requirement for faking cylindrical coordinates using Cartesian discretization. If you want to solve a 2D cylindrical problem, you need to create a 3D wedge shaped mesh. This may be possible in gmsh

Re: ... cylindrical coordinates ...

2007-01-17 Thread Matt Koch
. This exchange discusses wedge shaped meshes that are a requirement for faking cylindrical coordinates using Cartesian discretization. If you want to solve a 2D cylindrical problem, you need to create a 3D wedge shaped mesh. This may be possible in gmsh and then you can use the Gmsh importer to read

... cylindrical coordinates ...

2007-01-16 Thread Matt Koch
in order to switch on cylindrical coordinates? From the above entry, I am almost guessing that one would have to implement their own divergences and such in cylindrical coordinates by adding the radius in proper places of the divergences and such in rectangular coordinates? That can't be right