Benchmarks made in the PLEIADES platform

Fuel performance codes based on the Cast3M finite element solver

Numerous performances assessments were made within the PLEIADES platform : replacing fortran implementations by their MFront counterparts led to significant improvements, from \(30\%\) to \(50\%\) of the total computational times of some fuel performance codes developed in the platform.

This improvements were mainly due to fact that the behaviour integration schemes changed from explicit Runge-Kutta schemes to implicit ones. The main benefit of MFront was to grant users an easier access to the implicit schemes.

Cyrano3 fuel performance code

Relying on the specific modelling hypotheses supported by this code, namely axisymmetrical generalised plane stress and axisymmetrical generalised plane strain, highly specialised and efficient implementations of mechanical behaviours were developed in Cyrano3 fuel performance code [see 1] for both isotropic and orthotropic materials : numerical integration boils down to solving a scalar non linear equation in both cases and provides the consistent tangent operator.

The figure below compares the total computational times of a native implementation of a cladding behaviour to its equivalent MFront implementation: the last one appears to be competitive with the native implementation (the average computational time using the MFront implementation is sightly lower than the average computational time using the native implementation).

Comparing the total computational times of a native implementation of a cladding behaviour available in Cyrano3 to its equivalent MFront implementation

Benchmarks based on the Code-Aster finite element solver

Some benchmarks comparing the implementation generated through MFront to the native implementation provided by the Code-Aster finite element solver. Graphical illustrations shows that the results obtained with both implementations are indistinguishable.
Behaviour and test description Algorithm Total computational times (Code-Aster vs MFront) Graphical illustration
Visco-plastic and damaging for steel [see 2, 3] - 3D Notched specimen implying large deformation Implicit \(17mn 43s\) vs \(7mn 58s\)
Damaging for concrete [see 4, 5], 3D beam bending Default \(45mn\) vs \(63mn\)
Generic Single crystal viscoplasticity [see 7, 8], 3D aggregate, 300 grains Implicit \(28mn\) vs \(24mn\)
FCC single crystal viscoplasticity [8] , 2D specimen with displacement boundary conditions from EBSD experiment Ìmplicit \(33m54s\) vs \(29m30s\)
FCC homogeneized polycrystals 30 grains [8, see 10], unit testing Runge-Kutta 4/5 \(9s67\) vs \(8s22\)
Anisotropic creep with phase transformation, 3D pipe [see 11] Implicit \(180s\) vs \(171s\)

Developers of the Code-Aster general purpose finite element solver, made independent extensive tests, comparing their own native implementations to the ones generated with MFront, generally using an implicit scheme in both cases. Without discussing the very details of each test performed, several general conclusions can be drawn:

For a given behaviour, the development time was found significantly lower with MFront.

References

1.
Thouvenin, Gilles, Baron, Daniel, Largenton, Nathalie, Largenton, Rodrigue and Thevenin, Philippe. EDF CYRANO3 code, recent innovations. In : LWR Fuel Performance Meeting/TopFuel/WRFPM. Orlando, Florida, USA, September 2010.
2.
Mustata, R. and Hayhurst, D. R. Creep constitutive equations for a 0.5Cr 0.5 Mo 0.25V ferritic steel in the temperature range 565 °C - 675 °C. International Journal of Pressure Vessels and Piping. May 2005. Vol. 82, no. 5, p. 363–372. DOI 10.1016/j.ijpvp.2004.11.002. Available from: http://www.sciencedirect.com/science/article/pii/S0308016105000037
3.
EDF. R5.03.13 révision : 8886: Comportement viscoplastique avec endommagement de Hayhurst. Référence du Code Aster. EDF-R&D/AMA, 2012. Available from: http://www.code-aster.org
4.
Mazars, Jacky, Hamon, François and Grange, Stéphane. A new 3D damage model for concrete under monotonic, cyclic and dynamic loadings. Materials and Structures. October 2014. P. 1–15. DOI 10.1617/s11527-014-0439-8. Available from: http://link.springer.com/article/10.1617/s11527-014-0439-8
5.
EDF. R7.01.08 révision : 10461: Modèle d’endommagement de Mazars. Référence du Code Aster. EDF-R&D/AMA, 2013. Available from: http://www.code-aster.org
7.
Méric, L. and Cailletaud, Georges. Single crystal modelling for structural calculations. Journal of Engineering Material and Technology. January 1991. Vol. 113, p. 171–182.
7.
Méric, L. and Cailletaud, Georges. Single crystal modelling for structural calculations. Journal of Engineering Material and Technology. January 1991. Vol. 113, p. 171–182.
8.
EDF. R5.03.11 révision : 10623: Comportements élastoviscoplastiques mono et polycristallins. Référence du Code Aster. EDF-R&D/AMA, 2013. Available from: http://www.code-aster.org
9.
Monnet, G., Naamane, S. and Devincre, B. Orowan strengthening at low temperatures in bcc materials studied by dislocation dynamics simulations. Acta Materialia. January 2011. Vol. 59, no. 2, p. 451–461. DOI 10.1016/j.actamat.2010.09.039. Available from: http://www.sciencedirect.com/science/article/pii/S1359645410006166
10.
Berveiller, M. and Zaoui, A. An extension of the self-consistent scheme to plastically-flowing polycrystals. Journal of the Mechanics and Physics of Solids. October 1978. Vol. 26, no. 5–6, p. 325–344. DOI 10.1016/0022-5096(78)90003-0. Available from: http://www.sciencedirect.com/science/article/pii/0022509678900030
11.
EDF. R4.04.05: Modèle de comportement élasto-visqueux META_LEMA_ANI avec prise en compte de la métallurgie pour les tubes de gaine du crayon combustible. Documentation du Code-Aster. EDF-R&D/AMA, 2013. Available from: http://www.code-aster.org
12.
Chaboche, J. L. and Cailletaud, G. Integration methods for complex plastic constitutive equations. Computer method in applied mechanics and engineering. 1996. Vol. 133, p. 125–155.
13.
Brent, Richard P. Algorithms for Minimization Without Derivatives. Dover Publications, 1973. ISBN 9780486419985.