Re: [deal.II] Petrov-Galerkin stabilization approach to solve Navier-Stoke equations

2017-10-20 Thread Wolfgang Bangerth
On 10/20/2017 10:36 AM, Jie Cheng wrote: Jie -- I think this would be a very interesting program for others as well. Would you be interested in contributing it to the code gallery? I'd love to contribute after I work it out! Jie, that would be fantastic -- please do let us know

Re: [deal.II] Petrov-Galerkin stabilization approach to solve Navier-Stoke equations

2017-10-20 Thread Jie Cheng
Hi Martin and Wolfgang Thank you very much for the helpful comments and references. I'll start to read the works you mentioned. > Jie -- I think this would be a very interesting program for others as > well. > Would you be interested in contributing it to the code gallery? > I'd love

Re: [deal.II] Petrov-Galerkin stabilization approach to solve Navier-Stoke equations

2017-10-20 Thread Wolfgang Bangerth
On 10/19/2017 10:46 PM, Jie Cheng wrote: but if you treat the nonlinearity implicitly, and time discretized model will remain nonlinear and needs to be solved with a Newton scheme or similar. The appropriate solver for that system then remains FGMRES. Thank you for your

Re: [deal.II] Petrov-Galerkin stabilization approach to solve Navier-Stoke equations

2017-10-20 Thread Wolfgang Bangerth
On 10/20/2017 03:42 AM, Martin Kronbichler wrote: Then, there are also pressure reaction-convection-diffusion operators that also take convection into account. I have not seen them a lot for the full block system in the time-dependent case, probably because the range where they are really better

Re: [deal.II] Re: Calculating lift and drag

2017-10-20 Thread Martin Kronbichler
Hi Jie, One needs to move some of the vertices around. I have some code that does generate this geometry in one of my projects: https://github.com/kronbichler/adaflo/blob/master/tests/flow_past_cylinder.cc (check the create_triangulation function) The mesh that comes out from this code is

Re: [deal.II] Re: Calculating lift and drag

2017-10-20 Thread Martin Kronbichler
Hi Jie, One needs to move some of the vertices around. I have some code that does generate this geometry in one of my projects: https://github.com/kronbichler/adaflo/blob/master/tests/flow_past_cylinder.cc (check the create_triangulation function) The mesh that comes out from this code is

Re: [deal.II] Petrov-Galerkin stabilization approach to solve Navier-Stoke equations

2017-10-20 Thread Martin Kronbichler
Hi Jie, > Thank you for your hints on my previous question. Now I have > extended step-57 to be time-dependent. I added a u,t term which is > approximated using backward Euler, and all the remaining terms are > treated implicitly. Newton iteration is still used at every time step. > The system

[deal.II] Re: Signal: Segmentation fault in parallel::distributed::Triangulation

2017-10-20 Thread Daniel Arndt
Yiayng, Am Freitag, 20. Oktober 2017 01:42:09 UTC+2 schrieb Yiyang Zhang: > > Hello Daniel, > > I am using the following code on my macbook. > > mpirun -n 1 xterm -e lldb ./test : -n 1 ./test > So you are running into the error with only one MPI process as well? Then, you should be able to omit