private function spongeLayer_radial_polyn5_for_treeids(me, treeIds, tree, nDim, n) result(res)
This function calculates the sigma for the radial spongelayer
polynomial order 5 from treeids.
Sponge profile:
where, \sigma_A - sponge strength, L - thickness, x0 - start of sponge.
Profile is taken from:
Xu, Hui; Sagaut, Pierre (2013): Analysis of the absorbing layers for the
weakly-compressible lattice Boltzmann methods. In Journal of Computational
Physics 245, pp. 14-42. DOI: 10.1016/j.jcp.2013.02.051.
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed
arrows point from an interface to procedures which implement that interface.
This could include the module procedures in a generic interface or the
implementation in a submodule of an interface in a parent module.
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed
arrows point from an interface to procedures which implement that interface.
This could include the module procedures in a generic interface or the
implementation in a submodule of an interface in a parent module.