Tuesday, January 1, 2019

Kreiss-Oliger Spatial Filtering


Because of high frequency numerical noise, this version also includes Kriess-Oliger spatial filtering. This filter is implemented as:

Note that, since we are using  a 3-dimensional spatial grid, this filter is taken in all three spatial dimensions. That is to say, we perform the filter, for example, first on the x-dimension, then in the y-dimension, and then in the z-dimension. (Obviously, the order is not important here.)



No comments:

Post a Comment

Overview -- Numerical Relativity Using CUDA C/C++ Is Easier Than You Think!

Simulation results of a binary black hole head-on collision run on a GPU based home gaming computer Not too long ago, the only way ...