Re: [deal.II] Issue with convergence of iterative linear solver for system matrix in modified step-57 with no-normal flux constraints

2018-02-18 Thread Timo Heister
Are you setting the conditions correctly for the zero_constraints and the nonzero_constraints and are you sure you are not applying other boundary conditions on that wall (for example with interpolate_boundary_values())? On Sat, Feb 17, 2018 at 5:46 PM, Bruno Blais wrote: > Hello everyone, > > I

Re: [deal.II] Issue with convergence of iterative linear solver for system matrix in modified step-57 with no-normal flux constraints

2018-02-18 Thread Bruno Blais
Hello, My gmsh mesh produces 3 physical groups which are independent. Say my simulation is the flow around a cylinder. The cylinder is a physical group (0), the inlet (say left wall) is a class inheriting from the Function class which sets the X component to 1 and the Y component to 0 and uses t

Re: [deal.II] Issue with convergence of iterative linear solver for system matrix in modified step-57 with no-normal flux constraints

2018-02-19 Thread Timo Heister
> Does the order in which I apply the nonzero and zero constraints matter? These are two independent objects, so no. > Currently I apply the inlet and then the no-slip in the nonzero_constraints, > thus the bottom and top wall appear after the inlet. Afterward the cylinder > is put in the zero co

Re: [deal.II] Issue with convergence of iterative linear solver for system matrix in modified step-57 with no-normal flux constraints

2018-02-19 Thread Bruno Blais
Hello, Sorry I feel I have not explained myself correctly. Here is a drawing of the case: With Ux may be a profile or a constant. Initially I had set-up that cas

Re: [deal.II] Issue with convergence of iterative linear solver for system matrix in modified step-57 with no-normal flux constraints

2018-02-22 Thread Bruno Blais
I have found my problem. A boundary condition in my mesh was ill-defined and this lead to my error. Sorry for the confusion. Thanks On Monday, 19 February 2018 13:41:35 UTC-5, Bruno Blais wrote: > > Hello, > Sorry I feel I have not explained myself correctly. Here is a drawing of > the case: >