Difference between revisions of "Domain decomposition"

From Serpent Wiki
Jump to: navigation, search
Line 6: Line 6:
 
# geometry-based division into sectors and a central zone
 
# geometry-based division into sectors and a central zone
  
The method is illustrated by [[Domain decomposition example 1|full-core example]], where the depletion zone division and domain decomposition is defined by input lines:
+
The method is illustrated by a [[Domain decomposition example 1|full-core example]], where the depletion zone division and domain decomposition is defined by input lines:
 
+
The following full-core example illustrates the divisions.
+
  
 
  <nowiki>
 
  <nowiki>
Line 16: Line 14:
  
 
div fuel sep 1
 
div fuel sep 1
 
% --- Plot:
 
 
plot 3 1000 1000 0 -187.6 187.6 -187.6 187.6
 
  
 
% -- Domain decomposition (1 = simple, 2 = sector, 3 = sector + center)
 
% -- Domain decomposition (1 = simple, 2 = sector, 3 = sector + center)
Line 26: Line 20:
  
 
% ------------------------------------------------------------------------------</nowiki>
 
% ------------------------------------------------------------------------------</nowiki>
 +
 +
When the calculation is run with 5 MPI tasks, it produces the following output:
 +
 +
<nowiki>
 +
sss2 -mpi 5 core
 +
 +
(...)
 +
 +
ividing materials into depletion zones...
 +
 +
Material fuel:
 +
 +
- 63624 cells
 +
- 63624 depletion zones in total
 +
 +
Sorting lists (this may take a while)...
 +
OK.
 +
 +
Counting cells...
 +
OK.
 +
 +
Decomposing 63624 divided materials into 5 domains:
 +
 +
Domain 1: 12744 materials (20.0%)
 +
Domain 2: 12720 materials (20.0%)
 +
Domain 3: 12720 materials (20.0%)
 +
Domain 4: 12720 materials (20.0%)
 +
Domain 5: 12720 materials (20.0%)
 +
 +
(...)</nowiki>
 +
 +
and a geometry plot that shows the domains in different colors:
 +
 +
[[File:Dd_mode3.png|500px]]
 +
 +
Changing the mode to 2 (sector-wise division) produces:
 +
 +
[[File:Dd_mode2.png|500px]]
 +
 +
and the simple index-based division:
 +
 +
[[File:Dd_mode1.png|500px]]

Revision as of 11:56, 16 May 2019

Domain decomposition is used for saving computer memory by spatially dividing large burnup calculation problems between multiple calculation nodes. Each node handles particle transport only within its own domain, and particles crossing the domain boundaries are sent to another node to complete the history. The implementation in Serpent relies on a collision-based approach, i.e. the transfer is invoked when the particle undergoes a collision in a material which belongs to another domain. The decomposition affects only burnable materials divided into depletion zones using the div sep option. The most typical domain decomposition case is a full-core burnup calculation, in which the fuel is decomposed into separate domains to save memory. Interaction data for the other materials is shared by all parallel tasks.

The current implementation in Serpent 2 is relatively simple, and mainly designed for reactor cores. The number of domains is determined by the number of MPI tasks, and the division using the set dd input option. There are three division types:

  1. simple, based on depletion zone index
  2. geometry-based division into sectors
  3. geometry-based division into sectors and a central zone

The method is illustrated by a full-core example, where the depletion zone division and domain decomposition is defined by input lines:

% --- Domain decomposition -----------------------------------------------------

% --- Divide fuel into depletion zones:

div fuel sep 1

% -- Domain decomposition (1 = simple, 2 = sector, 3 = sector + center)

set dd 3

% ------------------------------------------------------------------------------

When the calculation is run with 5 MPI tasks, it produces the following output:

sss2 -mpi 5 core

(...)

ividing materials into depletion zones...

Material fuel:

 - 63624 cells
 - 63624 depletion zones in total

Sorting lists (this may take a while)...
OK.

Counting cells...
OK.

Decomposing 63624 divided materials into 5 domains:

Domain 1: 12744 materials (20.0%)
Domain 2: 12720 materials (20.0%)
Domain 3: 12720 materials (20.0%)
Domain 4: 12720 materials (20.0%)
Domain 5: 12720 materials (20.0%)

(...)

and a geometry plot that shows the domains in different colors:

Dd mode3.png

Changing the mode to 2 (sector-wise division) produces:

Dd mode2.png

and the simple index-based division:

Dd mode1.png