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]> > 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 > >
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
