Hello Tariq,
thank you for the additional information. I am doing such things the
following way:
Loop over all cells
if (cell->material_id () == A)
Loop over all faces
if (cell->neighbor (face)->material_id () == B)
reinit fe_face_values_A
get function values for face of cell A
reinit fe_face_values_B
get function values for face of cell B
compute integral
Best Regards,
Markus
Am 12.01.11 16:22, schrieb Tariq Baig:
Thanks for your advice,
it took me a while to think of an example that explains better what my
problem is.
I have two domains named A and B separated by a boundary G, by the way
suggested by Martin and Markus I can determine if my cell lies at G
and on which side.
The functional I want to assemble contains something like :
\int_G (F_A-F_B)*ðF_A dS
where F_B is the solution on domain B , F_A the solution on domain A
and ðF_A denotes the variation of the solution on domain A.
In the simplest case ðF_A is just a shape function, however the
problem of multiplying things from different sides of the boundary
remains.
greetings tariq
@Martin I am solving the problem of an elastic inclusion in a
mooney-rivlin elastic material
2011/1/11 Markus Bürg <[email protected] <mailto:[email protected]>>
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]:
<mailto:[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 listhttp://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
