Hi Dave, Adding the tag to the SOLVERINFO works (the assertion does not fail). If I try to achieve the same through adding it to the GLOBALSYSSOLNINFO, which I use, it does not help. I managed to create a quad mesh in Gmsh and use it with the coupled solver, and periodic BC. The difference was I did not use Transfinite keyword, so the mesh was slightly unstructured. What I noticed for this case, is that application of forcing through BodyForce had no effect on the flow. (This works like a charm with the VelocityCorrection). The flow was slowing down, and behaved as if no pressure gradient was present. Is it not, that since pressure is implicitly accounted for in the coupled solver, the forcing is just not applied? Not sure if this is a limitation of the coupled solver (not mentioned in the manual), or a bug. Hope this is helpful. Regards, Stan Gepner On 09.06.2015 04:49, David Moxey wrote:

Hi Stan,

Apologies for the delayed response. Sounds like you could be having some trouble with the multi-level static condensation when combined with the periodic BCs.

You could try, in your SOLVERINFO tag:

<I PROPERTY="GlobalSysSoln" VALUE="DirectStaticCond" />

and see if that works out? If that works I can try and look into why the multi-level is failing.

Cheers,

Dave

On 4 Jun 2015, at 18:45, Stanisław Gepner <sgepner@meil.pw.edu.pl> wrote:

...

Hi,

I think I found a problem when importing Gmsh meshes into Nektar. This happens if I try to use the coupled solver and the periodic BC. Works with VelocityCorrection, or if BC are standard D / N.

What I do: 1. Using the geometry in channel.geo create a mesh in Gmsh. This is a rectangular channel. I use Transfinite and Recombine to get a structured quad mesh. 2. I use MeshConvert with paralign to generate the channel.xml 3. If run with the attached cond2D.xml Nektar fails in AssemblyMap.cpp line 284, with "These values should never be the same"

The problem is persistent with other recombined/transfinite meshes.

If any of you is aware of how I could try to avoid the problem, I would be obliged.

Regards, Stan Gepner

<channel.msh><channel.xml><cond2D.xml><channel.geo>_______________________________________________ Nektar-users mailing list Nektar-users@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/nektar-users

-- David Moxey (Research Associate) d.moxey@imperial.ac.uk | www.imperial.ac.uk/people/d.moxey

Room 363, Department of Aeronautics, Imperial College London, London, SW7 2AZ, UK.