Re: cylindrical coordinates

2018-10-10 Thread Martinus WERTS
o 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 m

cylindrical coordinates

2018-09-20 Thread Martinus WERTS
conditions). Now I would like 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

AW: Gradient in cylindrical coordinates

2018-01-10 Thread Munoz Leyton
recipients of list 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/fipy

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
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 see the code I used below this line. >> >> &g

Re: Gradient in cylindrical coordinates

2018-01-08 Thread Daniel Wheeler
On Fri, Jan 5, 2018 at 7:23 PM, Munoz Leyton 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 > coordinates.

Gradient in cylindrical coordinates

2018-01-05 Thread Munoz Leyton
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 coordinates. You can see the code I used below this line. from fipy import * from fipy import CellVariable

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

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 > figu

Re: 2D Cylindrical Coordinates - Pipe

2012-06-20 Thread Kendall Boniface
On Tue, Jun 19, 2012 at 1:04 PM, Jonathan Guyer 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, ny=ny)

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 a

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 h

Re: 3D Transient Conduction - Cylindrical Coordinates

2012-06-06 Thread Daniel Wheeler
On Tue, Jun 5, 2012 at 1:08 PM, Kendall Boniface wrote: > > Hi Daniel, > I have set the edge lengths to 1 and placed the gmsh commands in a .geo file > as you suggested. I've attached the .geo file to this message. I also > removed the 4 lines and the last line of the gmsh commands as I think they

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: 3D Transient Conduction - Cylindrical Coordinates

2012-06-05 Thread Daniel Wheeler
On Mon, Jun 4, 2012 at 2:01 PM, Kendall Boniface wrote: > I have tried to set up a 3D mesh using Gmsh, but have not been very > successful. Could anyone offer some insight as to how I can set this up > successfully? I pasted the gmsh commands into a file and ran gmsh on it. It threw a weird erro

3D Transient Conduction - Cylindrical Coordinates

2012-06-04 Thread Kendall Boniface
conductivity of that material changes as a function of the temperature. My previous post was just looking at a simple 1D problem (just straight conduction through a slab of the material). Now I would like to look at it in 3D cylindrical coordinates. We have an insulated 1.5m long pipe with an inner and

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 cylinder

Re: ... cylindrical coordinates ...

2007-01-19 Thread Matt Koch
, if I am not too mistaken, I might as well view this like modelling a slice of pie in cartesian (rectangular) 3D space, and applying zero flux, etc., boundary conditions to the two flanks (planes of constant theta) of the slice? Exactly. We really are talking about a 3D model then. No need

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 o

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
; in this case, I'm going to attempt to clarify (but will probably fail miserably). Let's take our cylindrical coordinates to be r: radial distance from axis of cylinder theta: azimuthal angle around axis of cylinder z: distance along length of cylinder It is

Re: ... cylindrical coordinates ...

2007-01-18 Thread Daniel Wheeler
ell? No. 2D mesh objects are really 2D. The reason you would have 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 c

Re: ... cylindrical coordinates ...

2007-01-17 Thread Matt Koch
eems to speak to the subject, 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. Th

Re: ... cylindrical coordinates ...

2007-01-17 Thread Daniel Wheeler
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 and

... cylindrical coordinates ...

2007-01-16 Thread Matt Koch
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