1.3. Spectral Solver


The present implementation is based on the works by P. Suquet [3] and R.A. Lebensohn [4]. The »basic scheme« included in DAMASK is described in [1] and [2]. As an faster and more robust variant of the originally proposed fix-point algorithm, DAMASK includes also a variant using any suitable non-linear solver provided by PETSc. As an second alternative, an version using the Polarization field is included.


Solver variants




Convergence criterion

Since the iterative »basic scheme« (using either the fix-point algorithm or a non-linear solver provided by PETSc) strictly maintains compatibility, it is only necessary to achieve (static) equilibrium. As suitable convergence criterion Moulinec and Suquet [3] suggested

\begin{alignat}{1} \varepsilon_\text{eq} \ge \frac{\operatorname{RMS}(\operatorname{Div} \textbf{P}(\textbf x)) }{|| \bar{\textbf{P}}||} \times \text{m} \label{eq:convergenceCriterion} \end{alignat}

based on the the root mean square (RMS) value of $\operatorname{Div}(\tnsr P)$, which can be conveniently calculated in Fourier space employing Parseval's theorem.

Corrections for physical size and grid resolution

The magnitude of the divergence depends on the distance between two neighboring Fourier points, i.e. for the same number of Fourier points but different physical dimension the iterative scheme will stop after different numbers of iterations. To correct for a dependence on the size and resolution of the volume element (VE), the spectral solver can be instructed to rescale all grid point distances by the same factor such that for the dimension with the medium pixel resolution a grid point distance of

  1. as before (uncorrected)
  2. 1 meter / grid point number (fixed dimension)
  3. 1 meter (fixed grid point distance, default)
  4. 1 meter / square root of grid point number (dimension and square root resolution corrected)

is obtained.

Fixed grain structure at varying VE size and resolution

As example we take a cubic grain structure (included as example in DAMASK) at either 1 or 10 m side length and discretized by either 163, 323, or 643 Fourier points (FPs).

(a) 163

(b) 323

(c) 643
Figure 1: Periodic crystal of 20 grains (included in DAMASK as example) discretized at different resolutions.

Figure 2 presents the normalized divergence (see eq. \eqref{eq:convergenceCriterion}) for all four above rescaling options.

(a) uncorrected

(b) fixed dimension

(c) fixed grid point distance

(d) dimension and square root of resolution corrected
Figure 2: Evolution of convergence criterion with iterations for the exemplary polycrystal at various resolutions (basic fix-point rev. 2232, first step of tensile test load case).

Constant Grains per Fourier Point

(a) original

(b) 33 copies

(c) 36 copies
Figure 3: Polycrystal of 20 grains at different resolution periodic copies

(b) fixed dimension

(c) fixed grid point distance

(d) dimension and square root of resolution corrected
Figure 4: Divergence over iterations for periodic copies of VE (basic fix-point rev. 2232, first step of tensile test load case).

Convergence behavior

(a) Polycrystal, colors indicate orientation

(b) Dual-phase steel, martensite in dark blue, ferrite grains indicated by other colors
Figure 5: Input Data

(a) Polycrystal

(b) Dual-phase steel
Figure 6: Evolution of divergence and plastic strain (basic PETSc rev. 2232)



P. Eisenlohr, M. Diehl, R.A. Lebensohn, F. Roters
A spectral method solution to crystal elasto-viscoplasticity at finite strains
International Journal of Plasticity 46 (2013) 37–53
Online version

M. Diehl
A spectral method using fast Fourier transform to solve elastoviscoplastic mechanical boundary value problems
Diploma Thesis, TU München (2010)
Download here

H. Moulinec, P. Suquet
A numerical method for computing the overall response of nonlinear composites with complex microstructure
Computer Methods in Applied Mechanics and Engineering 157 (1998) 69–94
Online version

R.A. Lebensohn
N-site modeling of a 3D viscoplastic polycrystal using Fast Fourier Transform
Acta Materialia 49 (2001) 2723–2737
Online version

P. Shanthraj, P. Eisenlohr, M. Diehl, F. Roters
Numerically robust spectral methods for crystal plasticity simulations of heterogeneous materials
International Journal of Plasticity 66 (2015) 31–45
Online version

Topic revision: r41 - 03 May 2018, JaeyongJung

  • News
26 Mar 2019
DREAM.3D 6.5.119
(released 2019/03/22) comes with a DAMASK export filter
25 Mar 2019
Release of version v2.0.3
21 Jan 2019
DAMASK overview paper finally published with full citation information available
01 Dec 2018
DAMASK overview paper now online
17 Sep 2018
CMCn2018 & DAMASK user meeting to be hosted at Max-Planck-Institut für Eisenforschung
22 May 2018
Release of version v2.0.2
01 Sep 2016
CMCn2016 & DAMASK user meeting to be hosted at Max-Planck-Institut für Eisenforschung
25 Jul 2016
Release of version v2.0.1
08 Mar 2016
Release of version v2.0.0
22 Feb 2016
New webserver up and running
09 Feb 2016
Migrated code repository from Subversion to GitLab
17 Dec 2014
Release of revision 3813
14 May 2014
Release of revision 3108
02 Apr 2014
Release of revision 3062
16 Oct 2013
Release of revision 2689
15 Jul 2013
Release of revision 2555
15 Feb 2013
Release of revision 2174
13 Feb 2013
Doxygen documentation
16 Dec 2012
Powered by MathJax rendering
23 Nov 2012
Release of revision 1955
15 Nov 2012
Release of revision 1924
01 Nov 2012
Updated sidebar
30 Oct 2012
Significant website updates and content extensions

This site is powered by FoswikiCopyright by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding DAMASK? Send feedback
§ Imprint § Data Protection