Navier-Stokes solver convergence issues

Extension of Elmer in computational glaciology
Post Reply
LizR
Posts: 6
Joined: 14 Dec 2015, 08:39
Antispam: Yes
Location: ACE CRC Tasmania

Navier-Stokes solver convergence issues

Post by LizR »

Hello,

I am having some problems with the Navier-Stokes solver. I would like to perform a simple diagnostic test to apply pressure to the top of a cube of ice and get the stresses as an output. I have modified the ForceToStress test case to create a simple example. However, if I replace the ForceToStress solver with the Navier-Stokes solver it won't run. I keep getting this error - ERROR:: IterSolve: Failed convergence tolerances.

I have tried increasing the convergence tolerance but if I make it too large the answers won't mean anything and I don't think tight tolerances are the cause of the problem.

I have been trying to do essentially the same thing in a few different ways with different .sif files and different meshes but I keep hitting this problem. What am I missing?
Attachments
simple.sif
(3.12 KiB) Downloaded 561 times
cube.grd
(743 Bytes) Downloaded 559 times
raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: Navier-Stokes solver convergence issues

Post by raback »

Hi

Your problem has a large null space. Assume a solution and add constant velocity to that and it will also be a solution. Hence you cannot have a solution for a linear system. You should add some BCs which determine the velocity field uniquely.

Apart from that I don't really understand the setting. You seem to solve the N-S in a 2D plate where also the external pressure is given. It would seem more natural to give the BC for a lower dimensional entity.

-Peter
LizR
Posts: 6
Joined: 14 Dec 2015, 08:39
Antispam: Yes
Location: ACE CRC Tasmania

Re: Navier-Stokes solver convergence issues

Post by LizR »

Since the last post I have made the boundary conditions and the mesh geometry more realistic and I understand a lot more about what the boundary conditions do now.
However, I still can't seem to get the Navier-Stokes solver to do a diagnostic run. Using the ElasticSolver and equivalent BCs works just fine.
I don't understand what the difference is, and why the Navier-Stokes solution will not converge.
Any help would be much appreciated.
-Liz
Attachments
Tests.zip
Contains elastic and NS solver versions of the same problem
(422.1 KiB) Downloaded 516 times
raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: Navier-Stokes solver convergence issues

Post by raback »

Hi Liz

Elasticity solver has a Lagrangian formulation and the displacements are directly linked to the shape of the domain. One solver is therefore sufficient even with a moving surface.

The N-S solver is naturally written in Eulerien frame of reference. Velocity is not directly associated with the shape of the domain. Instead you need a free surface solver, a mesh adaptation solver, and an ALE formulation in the N-S solver. There are many examples among tests, and the Elmer/Ice community can probably provide you with real cases.

-Peter
LizR
Posts: 6
Joined: 14 Dec 2015, 08:39
Antispam: Yes
Location: ACE CRC Tasmania

Re: Navier-Stokes solver convergence issues

Post by LizR »

Hi Peter
Thank you very much for telling me this. It is good to know that the issue was things that were missing, not things that were incorrect per se. I will continue working on it.
-Liz
Post Reply