Steady state performance of WhitneyAVHarmonicSolver

General discussion about Elmer

Steady state performance of WhitneyAVHarmonicSolver

Postby vencels » 16 Aug 2017, 16:03


I am doing coupled Elmer and OpenFOAM simulation using EOF-Library.

Both meshes have ~1M elements/cells. Simulation is running on 8+8 MPI processes. For every Elmer steady state computation (it takes 170sec) there are ~8 OpenFOAM timesteps (all combined ~30sec). Basically, Elmer is taking about 85% of computational time.

I am looking for a way to improve Elmer's steady state performance, something like disabling unnecessary tasks. Considering that the problem is changing very little between two consecutive computations, I would assume that Elmer takes result from previous steady state, then does a couple more iterations until converges up to tolerance, calculates fields (current, B, J. etc..).

Any ideas how to boost performance?

Case and log file attached.
Case file
(3.42 KiB) Downloaded 14 times
Log file
(643.87 KiB) Downloaded 9 times
Posts: 22
Joined: 20 Sep 2016, 17:05
Location: Latvia

Re: Steady state performance of WhitneyAVHarmonicSolver

Postby raback » 16 Aug 2017, 23:49

Hi vencels

The assembly could be improved but most of that requires altering the code. Your material parameter is not constant so you really need to recompute the matrix every iteration. I would write F90 routines for the current densities. MATC really takes quite a lot of time. This should take easily some ~20% of the assembly time away.

Now for the linear system there is not that much that can be done. I trust you've tested different strategies. The matrix is really quite a beast and has a lot of connections. It is huge compared to the sparse matrices coming from similar mesh from the splitted schemes in FVM method.

What you could try is to have sloppier criteria at the start, say
Code: Select all
Linear System Convergence Tolerance = Variable "coupled iter"
    1.0 1.0e-3
    2.0 1.0e-5
    3.0 1.0e-7
    4.0 1.0e-7

Or you could say "Linear System Residual Mode = Logical True" and set the tolerance to ~1.0e-3 which is then an incremental increase in accuracy. If the latter works I would prefer that.

Also between OpenFOAM and Elmer you might require some relaxation. Also that could depend on iteration as above. It could be 1.0 the 1st iteration, and ~0.7 after that. Loosely coupled schemes often overshoot. This is of course very problem dependent...

Site Admin
Posts: 3188
Joined: 22 Aug 2009, 11:57
Location: Espoo, Finland

Return to General

Who is online

Users browsing this forum: No registered users and 3 guests