On Mar 28, 2016, at 3:04 PM, Djoeke Schoonenberg
wrote:
> Hmm, when I do as you suggest:
>
> xface = mesh.faceCenters[0]
>
> Eq = TransientTerm(var=S_vap) == \
> ConvectionTerm(coeff=vgas(xface)+p*Dif/xface, var=S_vap)\
> + ImplicitDiffusionTerm(coeff=Dif, var=S_vap)\
>
> I
On 28 Mar 2016, at 19:19, Guyer, Jonathan E. Dr. (Fed)
wrote:
>
> On Mar 28, 2016, at 12:20 PM, Djoeke Schoonenberg
> wrote:
>
>> So it seems that I’m implementing the variable speed in a wrong way. What I
>> did was ConvectionTerm(FaceVariable(mesh=mesh, value=vgas(x)), var=S_vap),
>> wh
Thanks, If you think DFT in FiPy is strange, you should see how I make
isosurfaces in POV-Ray. :)
That said, If anyone is tinkering with my example and has questions, I
am glad to answer them!
Cheers,
-Mike Waters
On 03/28/2016 01:46 PM, Guyer, Jonathan E. Dr. (Fed) wrote:
> Mike, thanks for
Mike, thanks for the example, and for the rather perverse application of FiPy!
I'll fiddle with this and see what I get.
- Jon
On Mar 25, 2016, at 7:16 PM, Michael Waters wrote:
> Hello,
>
> I still have a large memory leak when using Trilinos. I am not sure where to
> start looking so I mad
On Mar 28, 2016, at 12:20 PM, Djoeke Schoonenberg
wrote:
> So it seems that I’m implementing the variable speed in a wrong way. What I
> did was ConvectionTerm(FaceVariable(mesh=mesh, value=vgas(x)), var=S_vap),
> where vgas(x) is the speed as function of position. How should this be done
>
On Thu, Mar 24, 2016 at 6:36 PM, Djoeke Schoonenberg
wrote:
>
> I have a couple of questions/issues:
>
> 1) The analytical solution for S_vap (for a constant mass flux) is a 1/r
> profile —> its gradient is equal to - S_vap / r. However, setting this as a
> boundary condition at r1 (which updates
