> I am coding for my graduation design, and hoping to achieve a coupling > geomechanics and fluid model. > I use FE_RT<dim> for velocities and FE_Q<dim> for displacement of the rock. > > In order to valid my model and programming, I use 1 D Terzaghi. > But my boundary conditions appears have some problems.
I assume you then expect that the vertical velocity should be a 1d function independent of x and only dependent on y? The picture you attach clearly does not satisfy this. I don't know why this isn't so, but here are a few things you can check: * The vertical velocity you show is zero at the left and right sides of the domain. You say that you are applying "no flux" at boundaries 1, 2, 3 (which sides of the rectangle are these?) but that does not seem to be the case. I think that the function you are using, 'VectorTools::project_boundary_values_div_conforming', should do this correctly, but you may want to look at a small testcase with just the 1x10 mesh and output the constraints you get. There is a function in DoFTools that also lets you output the physical location of every degree of freedom so that you can correlate what you get from the constraints object to where these DoFs are located. Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
