So, there are two boundary integrals in the formulation and I am trying to implement them for the case of fixed displacement boundary condition. The first boundary integral falls to zero for such a BC, as no displacement variation is possible in this case. So, I am implementing the second integral from the second equation. The integral can be represented as in the BC.PNG file. There are those factors in brackets, assigned displacement values and variation of the pressure multiplied together:
*local_rhs*(i) += -porosity*(1.+Q[k]/R[k]) *d_boundary_values[q] *fe_face_values[pressure].value(i,q) *fe_face_values.JxW(q); But the way it is written now it cannot work, as the *d_boundary_values* term is a vector of vectors (as we have three displacement components). So, the question is, whether there is a way to treat the displacement components seperately in this *d_boundary_values* term? Or maybe even a smarter way to make it work? Please, let me know. Thank you a lot! четверг, 30 сентября 2021 г. в 14:28:26 UTC+2, Wolfgang Bangerth: > On 9/30/21 2:48 AM, Мария Бронзова wrote: > > > > Thank you a lot for the reply! I am trying to implement the Biot > equations for > > the propagation of acoustic waves in porous media. The equations are > attached > > in the "Biot .PNG" file. There are two boundary integrals there, > containing > > information on work, done by external forces on the fluid and skeleton > part of > > the media. I want to implement fixed displacement boundary conditions > > ("BC.PNG") on a part of the geometry with ID = 2, so I need to set the > > displacement components to zero. I was wondering, whether there is a way > to do > > it similarly to the way it"s done in step-20 example for the pressure > boundary > > condition (so, in the local_rhs vector). > > Maria: so which of the terms in that formula are you interested in > implementing? > > Best > W. > > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/7f753205-6fe6-4da2-93bb-340c87144688n%40googlegroups.com.