Re: [deal.II] using VectorTools::interpolate_boundary_values :Simple verification of Stokes system with zero velocity (using step-22)

2017-02-02 Thread Jane Lee
Thanks Wolfgang, very silly of me. Have sorted the other problem too which didn't have anything to do with the way i was interpolating boundaries. On Wednesday, February 1, 2017 at 6:30:00 PM UTC, Wolfgang Bangerth wrote: > > On 02/01/2017 11:24 AM, Jane Lee wrote: > > > > where the boundary v

Re: [deal.II] using VectorTools::interpolate_boundary_values :Simple verification of Stokes system with zero velocity (using step-22)

2017-02-01 Thread Wolfgang Bangerth
On 02/01/2017 11:24 AM, Jane Lee wrote: where the boundary values are just 1.0 so that u is a constant vector field. this still gives me the bizarre 22 factor in my pressures solution (which is just linear in z - the resulting pressure solution is atleast linear), but the pressure solution aside

Re: [deal.II] using VectorTools::interpolate_boundary_values :Simple verification of Stokes system with zero velocity (using step-22)

2017-02-01 Thread Jane Lee
Hi Wolfgang, thanks for a swift reply. and apologies, I wasn't very clear. when I'm imposing dirichlet conditions everywhere, i don't have a condition for pressure (those were two separate test cases - sorry i was unclear). and i'm using: VectorTools::interpolate_boundary_values (

Re: [deal.II] using VectorTools::interpolate_boundary_values :Simple verification of Stokes system with zero velocity (using step-22)

2017-02-01 Thread Wolfgang Bangerth
I'm trying to solve a stokes system with some funny boundary conditions so I've stripped the code to make sure I was getting the correct answers to a very simple system as follows: governing equations are: div grad u + grad p = f div u = 0 i'm setting u=0 on the boundaries, giving u=0 (or close

[deal.II] using VectorTools::interpolate_boundary_values :Simple verification of Stokes system with zero velocity (using step-22)

2017-02-01 Thread Jane Lee
Hi all, This may seem like a very obvious question/mistake on my part so please do excuse me if I've overseen something very simple. I'm trying to solve a stokes system with some funny boundary conditions so I've stripped the code to make sure I was getting the correct answers to a very simpl