# Discussion Forum

## computing accurate fluxes with periodic boundary conditions?

RSS feed   |   Turn on email notifications   |   2 Replies   Last post: January 26, 2013 12:42am UTC

Sabine

January 24, 2013 6:42pm UTC

computing accurate fluxes with periodic boundary conditions?

Hi everybody,

does anyone know if I can compute the fluxes accurately, as described in the user guide on page 308 ff. and file 0973 from the knowdledge database etc, when using periodic boundary conditions on the respective boundaries? In the descriptions, it is explained only for Dirichlet boundary conditions :( And if you think it should work, should I then use the weak constraint option? (and, how, if you have time to explain?) I model in the PDE mode of version 4.1. with matlab.

All best,
Sabine

Mattia Montanari

January 25, 2013 3:13pm UTC in response to Sabine

Re: computing accurate fluxes with periodic boundary conditions?

Hi Sabine,

well.. if you know (at least you know where to look for) how to deal with Diriclet BC's, then you have just 2 options left:

- Neumann BC's: in this case what you look for is exactly what you prescribed (trivial)

- Robin BC's: write as follow: ( a*x + b*d(x) ) = f(x)
this is nothing but a combination of the two previous cases. of course you know somehow a, b, f(x). On the boundary you know x from the solution vector, so immediately comes d(x) which is the flux that you want.

since you're working with Matlab this shouldn't require a big effort. I didn't read thourgh the manual, but this is something that you can write yourself directly in Matlab!

Hope it helps!
Mattia

Ivar Kjelberg

January 26, 2013 12:42am UTC in response to Sabine

Re: computing accurate fluxes with periodic boundary conditions?

Hi

another point I have noticed, is that to get good flux values along boundaries, you should ensure your mesh is fine while getting close to this boundary, by using edge pre-seeds, or even a boundary mesh (sometimes the latter requires further tweaking)

--
Good luck
Ivar

Rules and guidelines