Hello Tariq,
I am doing similar things for error estimation: I have to compute the
difference of the jumps across a face. Therefore I just created a second
FEFaceValues object and initialize it to the neighboring cell's face.
Would this also work for you?
Best Regards,
Markus
Am 11.01.11 14:50, schrieb [email protected]:
Dear all,
I am trying to solve a problem from structural mechanics involving two
(finite) elastic domains with different elastic constants.
In order for this problem to be well posed I have to specify internal
boundary conditions, in my case a vanishing "jump" of a function of
the solution at the boundary.
As a consequence the variational formulation of my problem contains
terms that involve objects on different sides of the internal boundary:
more precisely boundary integrals of functions assembled in a cell on
one side multiplied by test functions on the other side.
The best way to solve this problem I can think of is to figure out the
test/shape function in the same cell that corresponds to the one I
really want - which should be possible since
both cells are images of the same unit cell .
Has anybody experiences with something like this or maybe an idea that
solves the problem more elegantly?
Thanks in advance,
tariq
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
