Hi Ian,
One way to deal with the sharp discontinuity could be to have very
fine grid spacing near the discontinuity. The following code solves a
diffusion problem with a step in the diffusion coefficient. It only
uses 4 cells in 1D. 2 of the cells are very large and 2 are very
small. As the small
rent Diffusion Coefficients at Single FV Face
Yes!
> On Sep 16, 2016, at 2:18 PM, Raymond Smith wrote:
>
> A side note, here it may actually be more convenient to think about the D in
> terms of the volumes, so you could also define it as a cell variable with
> values specified as D1 wh
Yes!
> On Sep 16, 2016, at 2:18 PM, Raymond Smith wrote:
>
> A side note, here it may actually be more convenient to think about the D in
> terms of the volumes, so you could also define it as a cell variable with
> values specified as D1 when x<=1 and D2 when x>1, then use the
> harmonicFace
monicfacevalue ?
>
>
>
> Best Regards
>
>
>
> Krishna & Ian.
>
>
>
>
>
>
>
> *From:* fipy-boun...@nist.gov [mailto:fipy-boun...@nist.gov] *On Behalf
> Of *Raymond Smith
> *Sent:* Friday, September 16, 2016 7:18 PM
> *To:* fipy@nist.gov
> *Sub
.gov
Subject: Re: Applying Two Different Diffusion Coefficients at Single FV Face
A side note, here it may actually be more convenient to think about the D in
terms of the volumes, so you could also define it as a cell variable with
values specified as D1 when x<=1 and D2 when x>1, th
A side note, here it may actually be more convenient to think about the D
in terms of the volumes, so you could also define it as a cell variable
with values specified as D1 when x<=1 and D2 when x>1, then use the
harmonicFaceValue attribute when you put it into the governing equation to
have FiPy
Hi, Ian.
I don't think there is such a thing as having two different flux
coefficients at the same face for the same governing PDE. The flux through
a given face is calculated by the coefficient at that face times some
approximation of the gradient in a field variable at that face, like
D * grad(c