Hello Vladimir and Robert,
Vladimir: Thanks. In the problem I'm considering, the NBC is of the form \nabla V\cdot n = V^N (i.e. without the M). This complicates things a little -- see attached pdf.
Robert: Yes, that sounds like a good name. I've attached a pdf describing a demo problem which -- I think -- requires a term of this type. If I've got my calculations correct, you'll be able to test the new term against the analytical solution.
Thanks once again,
Ben
On Friday, March 29, 2013 3:16:15 PM UTC, Robert Cimrman wrote:Hello Ben,
On 03/29/2013 03:35 PM, Ben Derrett wrote:
> Hello,
>
> In the weak formulation of my problem (a 2D diffusion problem), I have the
> term
>
> \int_{\Gamma^N} p\nabla V\cdot M n
>
> where \Gamma^N is region where Neumann BC is applied, p is a test function,
> V is the unknown function, M is a 2x2 matrix and n is the normal vector.
> Has an SfePy term been written corresponding to this?
I do not think so. However it should not be difficult to add it - could you
provide an example file using a term like that, so that I could test the
possible implementation. It should run without the non-existing term. Also,
what name would you propose? dw_surface_s_grad_s_n??
> Alternatively, is there a SfePy term corresponding to the integral
>
> \int_{\Gamma^N} p\nabla V\cdot q,
>
> where q is a vector? I saw the term dw_surface_lcouple, but it doesn't seem
> to support the weak evaluation mode.
As above, but I assume that having the above term would suffice, right?
r.