salomemed
6.5.0

For fields with polynomial representation on each cell, the components of the discretized field on the source side can be expressed as linear combinations of the components of the discretized field on the target side, in terms of a matrixvector product:
is called the interpolation matrix. The objectives of interpolators is to compute the matrix W depending on their physical properties (Overview: intensive and extensive field) and their mesh discretisation (P0, P1,...).
At the basis of many CFD numerical schemes is the fact that physical quantities such as density, momentum per unit volume or energy per unit volume obey some balance laws that should be preserved at the discrete level on every cell.
It is therefore often desired that the process interpolation preserve the integral of on any domain. At the discrete level, for any target cell , the following general interpolation formula has to be satisfied :
This equation is used to compute , based on the fields representation ( P0, P1, P1d etc..) and the geometry of source and target mesh cells.
Another important property of the interpolation process is the maximum principle: the field values resulting from the interpolation should remain between the upper and lower bounds of the original field. When interpolation is performed between a source mesh S and a target mesh T the aspect of overlapping is important. In fact if any cell of of S is fully overlapped by cells of T and inversely any cell of T is fully overlapped by cells of S that is
then the meshes S and T are said to be overlapping and all the algorithms will return the same results.
The ideal interpolation algorithm should be conservative and respect the maximum principle. However such an algorithm can be impossible to design if the two meshes do not overlap. When the meshes do not overlap, using either or in the formula one obtains an algorithm that respects either conservativity either the maximum principle (see the nature of field summary table).
We assume that the field is represented by a vector with a discrete value on each cell. This value can represent either
For intensive fields such as mass density or power density, the left hand side in the general interpolation equation becomes :
Here Vol represents the volume when the mesh dimension is equal to 3, the area when mesh dimension is equal to 2, and length when mesh dimension is equal to 1.
In the general interpolation equation the right hand side becomes :
As the field values are constant on each cell, the coefficients of the linear remapping matrix are given by the formula :
In code coupling from neutronics to hydraulics, extensive field of power is exchanged and the total power should remain the same. The discrete values of the field represent the total power contained in the cell. Hence in the general interpolation equation the left hand side becomes :
while the right hand side is now :
The coefficients of the linear remapping matrix are then given by the formula :
In the case of fields with P0 representation, if the meshes do not overlap the scheme is either conservative or maximum preserving (not both) and depending on the prioritised property and the Nature of a field the interpolation coefficients take the following value
Intensive  extensive  
Conservation 
RevIntegral 
IntegralGlobConstraint 
Maximum principle 
ConservativeVolumic 
Integral 
Let's consider the following case with a source mesh containing two cells and a target mesh containing one cell. Let's consider a field FS on cells on the source mesh that we want to interpolate on the target mesh.
The value of FS on the cell#0 is 4 and the value on the cell#1 is 100.
The aim here is to compute the interpolated field FT on the target mesh of field FS depending on the nature of the field.
The first step of the interpolation leads to the following M1 matrix :
If we apply the formula above it leads to the following matrix :
As we can see here the maximum principle is respected.This nature of field is particulary recommended to interpolate an intensive field such as temperature or pression.
If we apply the formula above it leads to the following matrix :
This type of interpolation is typically recommended for the interpolation of power (NOT power density !) for a user who wants to conserve the quantity only on the intersecting part of the source mesh (the green part on the example)
This type of interpolation is equivalent to the computation of followed by a multiplication by where is given by :
In the particular case treated here, it means that only a power of 25.055 W is intercepted by the target cell !
So from the 104 W of the source field , only 25.055 W are transmited in the target field using this nature of field. In order to treat differently a power field, another policy, integral global constraint nature is available.
If we apply the formula above it leads to the following matrix :
This type of interpolation is typically recommended for the interpolation of power (NOT power density !) for a user who wants to conserve all the power in its source field. Here we have 104 W in source field, we have 104 W too, in the output target interpolated field.
BUT, As we can see here, the maximum principle is not respected here, because the target cell #0 has a value higher than the two intercepted source cells.
If we apply the formula above it leads to the following matrix :
This type of nature is particulary recommended to interpolate an intensive density field (moderator density, power density). The difference with conservative volumic seen above is that here the target field is homogeneized to the whole target cell. It explains why this nature of field does not follow the maximum principle.
To illustrate the case, let's consider that is a power density field in . With this nature of field the target cell #0 cumulates 0.125*4=0.5 W of power from the source cell #0 and 0.75*100=75 W of power from the source cell #1. It leads to 75.5 W of power on the whole target cell #0. So, the final power density is equal to 75.5/1.5=50.333 W/m^2.