Nonlinearly-Stable High-Order Methods on Simplices with Improved Efficiency - presented by David Zingg

Nonlinearly-Stable High-Order Methods on Simplices with Improved Efficiency

David Zingg

DZ
Slide at 22:00
Tensor-Product Split-Simplex Operators:
Error convergence comparisons
SBP-E
SBP-E p=1 0p=2-0p=3*p=4-0p=5
TPSS p=1-p=2-p=3-p=1- -p=5
TPSS
p=1 p=2 p=3 p=4 p=5
10-2
10-1
10-3
10-2
10-3
10-5
10-4
10-5
10-5
10-4
10-8
10-7
10-6
10-5
10-4
1/dof
1/dof
Two dimensions
Three dimensions
Compared to SBP-E operator for isentropic vortex problem
1D operators based on Legendre-Gauss-Lobatto quadrature rules
Experiments on highly distorted meshes show comparable efficiency gains
Share slide
Summary (AI generated)

Here are various nodal distributions that illustrate the underlying tensor product structures across three quadrants. To implement this effectively, your one-dimensional operator must be one order higher than the desired final multidimensional operator, specifically D minus one orders higher.

In the results presented, we compare 2D representations on the left with 3D representations on the right. The notation used includes TPSS, which stands for Tensor Product Split Simplex, and SVP, which refers to the Split Boundary Operator.

The open symbols represent the multidimensional operator, while the solid symbols denote the new operators. By matching colors, you can observe a significant advantage in freedom and solution error.

This example focuses on an isotropic vortex problem. Notably, we applied extensive warping and skewing to the meshes. Interestingly, despite these distortions, we are not encountering the expected issues typically associated with highly distorted meshes. Additionally, we present normalized comparisons of residual computation times.