Thanks a lot for the implementation of anisotropic viscosities!
I haven't had the time to test it for now, but I'll do it as soon as possible. I will let you know the tests that I have done and how they went.
Thanks again,
Kévin
Search found 15 matches
- 13 Sep 2021, 09:23
- Forum: Software development
- Topic: GOLF rheology in ParStokes Solver?
- Replies: 8
- Views: 5344
- 02 Sep 2021, 11:16
- Forum: Software development
- Topic: GOLF rheology in ParStokes Solver?
- Replies: 8
- Views: 5344
Re: GOLF rheology in ParStokes Solver?
I have attached a modified version of the ParStokes Solver that computes an anisotropic viscosity under a (6x6) matrix form. I tried to follow the original structure of the ParStokes solver. Notably, the computation of the effective viscosity matrix is done in a separate module file, similar to orig...
- 31 Aug 2021, 10:33
- Forum: Software development
- Topic: GOLF rheology in ParStokes Solver?
- Replies: 8
- Views: 5344
Re: GOLF rheology in ParStokes Solver?
I have found one paper that uses the isotropic part of the tensorial viscosity to compute the Schur complement: https://www.sciencedirect.com/science/article/pii/S0031920114001514#b0085 . In the article they say that they use "the lumped mass matrix weighted by the inverse of the isotropic part...
- 27 Aug 2021, 15:51
- Forum: Software development
- Topic: GOLF rheology in ParStokes Solver?
- Replies: 8
- Views: 5344
Re: GOLF rheology in ParStokes Solver?
Hi Mika, Thanks a lot for your answer. I have done a bit a bibliography, and indeed there appears to be no clear and easy way to define a good approximation of the Schur complement with anisotropic viscosities. From what I understand (and this is not my original field of competence so I might be wro...
- 23 Aug 2021, 09:42
- Forum: Software development
- Topic: GOLF rheology in ParStokes Solver?
- Replies: 8
- Views: 5344
GOLF rheology in ParStokes Solver?
Hello Everyone, In Brief: I would like to use the anisotropic GOLF material law with the ParStokes solver for large simulations. It requires some modification of the ParStokes.F90 file, in order to use non scalar viscosities, but I do not know how to concretly do it. My Problem: I have tried to perf...
- 30 Jul 2021, 10:45
- Forum: Bug reports
- Topic: Linear Solver - Different behaviours on different OS
- Replies: 5
- Views: 3828
Re: Linear Solver - Different behaviours on different OS
Hi everyone, Sorry for the long response, I did not find the time until recently to test the issue a bit more. The strain-rate solver tells me that I indeed have a large variability in terms of local strain rates in my structure (basically going down to the Critical Shear Rate). This suggests that w...
- 08 Jun 2021, 16:59
- Forum: Bug reports
- Topic: Linear Solver - Different behaviours on different OS
- Replies: 5
- Views: 3828
Re: Linear Solver - Different behaviours on different OS
Thanks for the response (and sorry for the delay on my side, I wanted to do a few more tests before posting). The two Elmer installation on Ubuntu and Centos were manually compiled (no packaged binaries involved). I tried installing several version of Elmer on the Ubuntu machine (sources some from E...
- 03 Jun 2021, 12:00
- Forum: Bug reports
- Topic: Linear Solver - Different behaviours on different OS
- Replies: 5
- Views: 3828
Linear Solver - Different behaviours on different OS
Hello everyone, I am currently experiencing a strange behaviour of Elmer, and I really do not understand its origin. In a nutshell: though using the same Elmer version (9.0), same sif file and same mesh, I am getting very different results on different computers. My workflow is simply to solve the s...
- 22 Mar 2021, 19:35
- Forum: General
- Topic: Cannot read shared objects when using a User Function on HPC
- Replies: 4
- Views: 1775
Re: Cannot read shared objects when using a User Function on HPC
I recompile the USF on the HPC every time I launch a new batch, just before launching the Elmer simulation. My batch script looks like that: # Loading elmer and required libraries in the user environment module load netlib-scalapack-2.1.0-gcc-9.1.0-4dzvobu module load openmpi-4.0.5-gcc-9.1.0-erqp5sa...
- 22 Mar 2021, 18:57
- Forum: General
- Topic: Cannot read shared objects when using a User Function on HPC
- Replies: 4
- Views: 1775
Re: Cannot read shared objects when using a User Function on HPC
Hi Peter, Thanks for the very fast answer. Do you compile the .so on the HPC server? I guess then you must have access to the libraries. I did not personally compiled the .so file on the HPC, it was taken care of by the IT team. I checked the .so files on the HPC, and they have the permission rwxr-x...