Difference between revisions of "Hyperviscosity"
From Medusa: Coordinate Free Mehless Method implementation
Zigavaupotic (talk | contribs) (Created page with "One of the most commonly mentioned drawbacks of the RBF-FD method is its stability. Furthermore, the RBF-FD method, due to scattered nodes, sometimes produces system (differen...") |
Zigavaupotic (talk | contribs) |
||
Line 1: | Line 1: | ||
− | One of the most | + | One of the most common drawbacks of the RBF-FD method is its stability. Furthermore, the RBF-FD method, due to scattered nodes, sometimes produces system (differential) matrices with spurious eigenvalues. Those are far more common than with for example FVM or FEM method. |
− | + | For such unstable systems, one usually considers numerical diffusion in the form of hyperviscosity. The hyperviscosity scheme introduces higher-order Laplacian operator (i.e. hyperviscosity operator) to the right hand side of the PDE | |
$$ | $$ | ||
− | +\gamma \Delta^\alpha u | + | +\gamma \Delta^\alpha u, |
$$ | $$ | ||
− | + | where $\gamma$ is a hyperviscosity constant and $\alpha$ is the order of hyperviscosity. |
Latest revision as of 15:24, 23 March 2024
One of the most common drawbacks of the RBF-FD method is its stability. Furthermore, the RBF-FD method, due to scattered nodes, sometimes produces system (differential) matrices with spurious eigenvalues. Those are far more common than with for example FVM or FEM method.
For such unstable systems, one usually considers numerical diffusion in the form of hyperviscosity. The hyperviscosity scheme introduces higher-order Laplacian operator (i.e. hyperviscosity operator) to the right hand side of the PDE
$$ +\gamma \Delta^\alpha u, $$
where $\gamma$ is a hyperviscosity constant and $\alpha$ is the order of hyperviscosity.