Re: [deal.II] Re: Replacement of FE_Q with FE_Bernstein fails

2019-01-16 Thread Wolfgang Bangerth
On 1/17/19 12:43 AM, 'Maxi Miller' via deal.II User Group wrote: > For convergence tests I wanted to replace the commonly used elements FE_Q > with > FE_Bernstein. Unfortunately apparently it is not possible to simply replace > all occurences. When running my program, I get the following error:

Re: [deal.II] How to implement GridGenerator::channel_with_cylinder?

2019-01-16 Thread Jean-Paul Pelteret
Dear Nick, Great, thanks for being willing to add this feature! To get you started, here are some guidelines as to how to use git to create a feature branch and make a pull request to the deal.II github repository. https://github.com/dealii/dealii/wiki/Contributing

[deal.II] Re: Replacement of FE_Q with FE_Bernstein fails

2019-01-16 Thread 'Maxi Miller' via deal.II User Group
Update: It looks as if this issue: https://github.com/dealii/dealii/issues/2237 might be related to the currently occuring problem? Am Donnerstag, 17. Januar 2019 08:42:01 UTC+1 schrieb Maxi Miller: > > For convergence tests I wanted to replace the commonly used elements FE_Q > with

[deal.II] Replacement of FE_Q with FE_Bernstein fails

2019-01-16 Thread 'Maxi Miller' via deal.II User Group
For convergence tests I wanted to replace the commonly used elements FE_Q with FE_Bernstein. Unfortunately apparently it is not possible to simply replace all occurences. When running my program, I get the following error: An error

Re: [deal.II] How to implement GridGenerator::channel_with_cylinder?

2019-01-16 Thread nboychenko1
Thank you David, Thank you Jean-Paul for the hints. Fro the contribution of such functionality, I will be glad to contribute and will really appreciate your help. Nick. On Wednesday, January 16, 2019 at 11:52:09 AM UTC-5, Jean-Paul Pelteret wrote: > > Dear Nick, > > Indeed, it would be quite

Re: [deal.II] Re: Problem with Block Preconditioner and Fully Periodic Boundary Conditions

2019-01-16 Thread Wolfgang Bangerth
On 1/16/19 10:19 AM, ntets...@berkeley.edu wrote: > > thanks for the reply. I'm working on Immersed Boundary Method and flagella > swimming, and, more specifically, I'm trying to model a unit square domain, > which is infinitely replicated in two dimensions, governed by the > incompressible

Re: [deal.II] apply Pouseuille BC at multiple inlets with different velocity magnitudes

2019-01-16 Thread Lailai Zhu
Thanks very much, Jean-Paul, during the time i was waiting for the answer, i kind of figured out a way similar to the first one, now since you have also proposed a more elegant way, i will look into it also. Thanks, best, lailai On Wednesday, January 16, 2019 at 12:08:18 PM UTC-5, Jean-Paul

Re: [deal.II] Re: Problem with Block Preconditioner and Fully Periodic Boundary Conditions

2019-01-16 Thread ntetsika
Hi Wolfgang, thanks for the reply. I'm working on Immersed Boundary Method and flagella swimming, and, more specifically, I'm trying to model a unit square domain, which is infinitely replicated in two dimensions, governed by the incompressible time-dependent stokes equations for low Reynolds

Re: [deal.II] apply Pouseuille BC at multiple inlets with different velocity magnitudes

2019-01-16 Thread Jean-Paul Pelteret
Dear Laila, I have two suggestions that may help you here: 1/ You could add the boundary ID as an argument to the PoiseuilleBoundaryValues class: template class PoiseuilleBoundaryValues : public Function { public: PoiseuilleBoundaryValues (const typename types::boundary_id b_id)

Re: [deal.II] How to implement GridGenerator::channel_with_cylinder?

2019-01-16 Thread Jean-Paul Pelteret
Dear Nick, Indeed, it would be quite straightforward to extend the GridGenerator::plate_with_a_hole() function to one that has a filled inclusion using the steps that David has listed.

[deal.II] Re: How to implement GridGenerator::channel_with_cylinder?

2019-01-16 Thread David Wells
Hi Nick, I think the following would work: 1. Create a triangulation via GridGenerator::plate_with_a_hole. 2. Create a second Triangulation via GridGenerator::hyper_ball. If you are working in 3D then you will have to extrude this. Use GridTools::transform to recenter this Triangulation in the

[deal.II] Re: Implement of adaptive mesh refinement cannot work in time-dependent problem.

2019-01-16 Thread Daniel Arndt
Chucui, Thanks for your quick answer! In my deal.ii (version 8.5.1), the code also > fails in the 9th time step. Shall I need to reinstall my deal.ii into the > newest version? > Normally, we only fix bugs in the library in the developer version which is also very stable unless you are

[deal.II] Re: Implement of adaptive mesh refinement cannot work in time-dependent problem.

2019-01-16 Thread chucui1016
Dear Prof. Arndt, Thanks for your quick answer! In my deal.ii (version 8.5.1), the code also fails in the 9th time step. Shall I need to reinstall my deal.ii into the newest version? For the shortest test code, please give me some time. Because I have deleted about 2000 lines. If I delete

[deal.II] Re: Implement of adaptive mesh refinement cannot work in time-dependent problem.

2019-01-16 Thread Daniel Arndt
Chucui, with a recent developer version your code fails for me in the 9th time step when the solver fails to converge within 1000 iterations. Hence, I can't reproduce the problem with a newer deal.II version. Please provide a minimal example by removing all functionality that is not necessary