Dear all

I have included below an example code which converges with

commit ed1bdc6fc1931a664bc30920f66f054907a49cce

but fails to converge with

commit d1e3ffcaf6733e7d8509d93785b1f9e2b16a977e

In this code, I solve the Winslow equations using continuous finite
elements. I am using curved boundaries and I could reproduce the failure
only with a complicated manifold.

The manifold I use is a NACA0012 airfoil for which I have implemented a
*project_to_manifold* function, see *winslow/parallel/naca.h*. I project
the candidate point onto the true boundary in the normal direction to the
boundary edge.

The code is here

https://bitbucket.org/cpraveen/deal_ii/src/master/winslow/parallel/?at=master

To get the code
git clone https://cprav...@bitbucket.org/cpraveen/deal_ii.git

The code is in

deal_ii/winslow/parallel

Steps to run the code

cd deal_ii/winslow/parallel
cmake .
make
cd ../run
gmsh -2 naca_struct.geo
../parallel/main 3

I realise that this is a bit complicated example but I dont know what to
remove to simplify it. Any tips will be useful.

This problem could be related to what Luca wrote here

https://groups.google.com/forum/#!searchin/dealii/SphericalManifold$20and$
20accuracy$20problem%7Csort:relevance/dealii/o33ApWoDnPQ/4MGY6SjMBQAJ

Best
praveen

On Tue, Aug 23, 2016 at 10:36 PM, Praveen C <cprav...@gmail.com> wrote:

> So it looks like the problem is with
>
> commit d1e3ffcaf6733e7d8509d93785b1f9e2b16a977e
> Author: Luca Heltai <luca.hel...@sissa.it>
> Date:   Tue Jul 12 12:39:30 2016 +0200
>
>     Added get_new_point with two points and a weight.
>
> I will try to send my code which is affected by the "bug", but it may not
> be very simple.
>
> Best
> praveen
>

-- 
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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply via email to