Search found 15 matches

by kfourteau
13 Sep 2021, 09:23
Forum: Software development
Topic: GOLF rheology in ParStokes Solver?
Replies: 8
Views: 5251

Re: GOLF rheology in ParStokes Solver?

Thanks a lot for the implementation of anisotropic viscosities! :D

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
by kfourteau
02 Sep 2021, 11:16
Forum: Software development
Topic: GOLF rheology in ParStokes Solver?
Replies: 8
Views: 5251

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...
by kfourteau
31 Aug 2021, 10:33
Forum: Software development
Topic: GOLF rheology in ParStokes Solver?
Replies: 8
Views: 5251

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...
by kfourteau
27 Aug 2021, 15:51
Forum: Software development
Topic: GOLF rheology in ParStokes Solver?
Replies: 8
Views: 5251

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...
by kfourteau
23 Aug 2021, 09:42
Forum: Software development
Topic: GOLF rheology in ParStokes Solver?
Replies: 8
Views: 5251

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...
by kfourteau
30 Jul 2021, 10:45
Forum: Bug reports
Topic: Linear Solver - Different behaviours on different OS
Replies: 5
Views: 3763

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...
by kfourteau
08 Jun 2021, 16:59
Forum: Bug reports
Topic: Linear Solver - Different behaviours on different OS
Replies: 5
Views: 3763

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...
by kfourteau
03 Jun 2021, 12:00
Forum: Bug reports
Topic: Linear Solver - Different behaviours on different OS
Replies: 5
Views: 3763

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...
by kfourteau
22 Mar 2021, 19:35
Forum: General
Topic: Cannot read shared objects when using a User Function on HPC
Replies: 4
Views: 1769

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...
by kfourteau
22 Mar 2021, 18:57
Forum: General
Topic: Cannot read shared objects when using a User Function on HPC
Replies: 4
Views: 1769

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...