Radial Basis Function¶

This interpolation filter interpolates field data from a source mesh to a target mesh based on radial basis functions (RBFs), where the local Wendland kernel together with a modified Shepard's method was chosen. Thereby, node-based data can be interpolated either to the cell centroids or the nodes of the source mesh, whereas element-based data can currently only be interpolated to the cell-centroids. The theoretical background of the implemented interpolation scheme is published in ¹. When using the conservative filters please provide a citation in you publication:

Schoder, Stefan, et al. "Aeroacoustic source term computation based on radial basis functions." International Journal for Numerical Methods in Engineering (2020).

Like for all interpolation filters, the input and output quantitiy, the source and target regions, and the target mesh have to be defined within the xml-scheme. By using the optional tag noSlipWall, the output quantity of the nodes on the specified wall region (boundary surface) are set to zero. This tag is intended to consider no-slip-walls for the interpolation of the flow velocity but can theoretically also be used for other quantities.

The field data is interpolated to the nodes of the target mesh by default. By setting useElemAsTarget true, the field data is interpolated to the cell centroids of the target mesh. The tag interpolationExponent handles the locality of the approximation on the target mesh and is considered within the weight function (for details see ¹). The larger it is, the more local the approach making it less accurate but capable of resolving stronger gradients. Furthermore, globalFactor allows a scaling of the output by the defined factor.

The number of neighbour points $N_q$ (numNeighbours) and influence points $N_w$ (numNeighbours_weight) are by default set to $N_q=18$ and $N_w=13$ , which is a compromise of accuracy and numerical efficiency for 2D and 3D interpolations. The interopolation result can be enhanced by tuning these parameters specifically for the involved source and target mesh. By setting the tag useCGAL4RBF to true, the nearest neighbor search will be carried out be by CGAL, which requires the CGAL 3.3 build (cmake flag must be set in the cmake build properties).

<interpolation type="FieldInterpolation_RBF" id="interpolationRBCF" inputFilterIds="input">
    <IntSchemeRBF numNeighbours="..." numNeighbours_weight="..." globalFactor="..." useCGAL4RBF="false"/>
        <useElemAsTarget>false</useElemAsTarget>
    <noSlipWall name="noSlipWall"/>
    <targetMesh>
        <hdf5 fileName="pathToMeshFile/meshFile.cfs"/>
    </targetMesh>
    <singleResult>
                <inputQuantity resultName="inputQuantity"/>
                <outputQuantity resultName="outputQuantity"/>
    </singleResult>
    <regions>
        <sourceRegions>
            <region name="sourceRegionR" />
        </sourceRegions>
        <targetRegions>
                    <region name="targetRegion"/>
                </targetRegions>
    </regions>
</interpolation>

Acknowledgement¶

Please provide an acknowledgement at the end of your publication using this software part for simulations

The computational results presented have been achieved [in part] using the software openCFS [Radial Basis Function].

References¶

Stefan Schoder, Klaus Roppert, Michael Weitz, Clemens Junger, and Manfred Kaltenbacher. Aeroacoustic source term computation based on radial basis functions. International Journal for Numerical Methods in Engineering, 121(9):2051–2067, 2020. ↩↩