Difference between revisions of "Hyperviscosity"

From Medusa: Coordinate Free Mehless Method implementation
Jump to: navigation, search
(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...")
 
 
Line 1: Line 1:
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 (differential) matrices with spurious eigenvalues. Those are far more common than with for example FVM or FEM method. The instabilities arise fastly even with linear PDEs.
+
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.  
  
With meshless methods (especially RBF-FD) one usually stabilises such schemes by adding a higher-order Laplacian (hyperviscosity) operator on the right-hand side.
+
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,
 
$$
 
$$
  
The final scheme is referred to as hyperviscosity scheme
+
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.