RE: casting implicit Boundary Conditions in FiPy

Hi Raymond, Sorry, it was a typo. Yes, It is indeed d (phi)/dx, the spatial derivative BC. I shall try setting phi.faceGrad.constrain([k*phi], mesh.facesRight), and see if it will work. Thanks for pointing this out. Krishna From: fipy-boun...@nist.gov [mailto:fipy-boun...@nist.gov] On

Re: how to create large grid with holes and odd geometry

Thank you. Because of the regular nature of the original data, it is easy to make it all triangles. Cheers, Jamie On Wed, Jun 8, 2016 at 2:10 PM, Guyer, Jonathan E. Dr. (Fed) < jonathan.gu...@nist.gov> wrote: > Meshes with holes are not a problem for FiPy. Daniel will be happy to help > you

Re: how to create large grid with holes and odd geometry

Meshes with holes are not a problem for FiPy. Daniel will be happy to help you create a Mesh2D from the output of the triangle package. Basically, you need a list of vertex coordinates, a list of vertex IDs that make up faces, and a list of faces that make up cells. Having all triangles should

Re: how to create large grid with holes and odd geometry

Hi Jamie, The simplest way might be just to take a basic grid and just "switch off" the simulation outside of the black zone. This can be achieved by including a transient term in the equation and setting the coefficient to be large in the white zone and zero or something small in the black zone.

Re: how to create large grid with holes and odd geometry

If the domain were not so large and so sparse, I'd be inclined to create a simple, rectilinear Grid2D of the full extent and then use known coefficients to mask out (set B to zero?) the solution where you don't know/care. Assuming the axes are labeled in grid spacings (?), then your mesh would

Re: Matrix coefficient for FiPy Transient Term

There's no need to create, extract, and reshape the result. Just use the mesh's y coordinates: TransientTerm(coeff=mesh.y**2) == DiffusionTerm(coeff=D) should work. > On Jun 8, 2016, at 9:01 AM, Gopalakrishnan, Krishnakumar > wrote: > > src="

Matrix coefficient for FiPy Transient Term

https://cdn.rawgit.com/google/code-prettify/master/loader/run_prettify.js";> Hello, The following PDE is defined in a non-square 2D (x-y) cartesian grid (20 x 10) domain $$ \frac{\partial}{\partial t} (y^2 \phi(x,y,t)) = D \nabla^2 \phi(x,y,t) $$ wherein $y^2$ is the 'y' co-ordinate of