Georgios, I am working on a combined stokes laplace problem (I also attach the code). > Essentially I am solving stokes in one part of the domain and laplace in > another. To this extend, I want my velocity to be both continuous and and > have continuous first derivatives everywhere in my domain (including on > between interfaces). My implementation is similar to step 46. I am > choosing my spaces like I show below (in the stokes part and in the laplace > part of my domain respectively): > > stokes_fe (FE_Q<dim>(stokes_degree+2), dim, > > FE_Q<dim>(stokes_degree+1), 1), > > laplace_fe (FE_Q<dim>(stokes_degree+2), dim, > > FE_Nothing<dim>(), 1), > > > > spaces with the intention of getting zero jump in the velocity gradient > between the interfaces, however, this does not happen. I have taken care > to choose my pressure in L^2_0 on both parts of the domain. Any suggestion > on how I can enforce continuity of velocity gradients across the interface? > The (discrete) ansatz spaces you have chosen only imply continuity and continuity of the gradients in tangential direction across faces. Hence, you propably would need to modify your ansatz spaces to get continuity in normal direction. >From a mathematical point of view, I would expect weak differentiability to be sufficient though. What is the (continuous) solution space for the continuous formulation of your problem?
In the end, the question is: Why do you need continuous gradients? Best, Daniel -- 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
