In this article,discrete variants of several results from vector calculus are studied for clas-sical finite difference summation by parts operators in two and three space dimensions.It is shown that existence theorems...In this article,discrete variants of several results from vector calculus are studied for clas-sical finite difference summation by parts operators in two and three space dimensions.It is shown that existence theorems for scalar/vector potentials of irrotational/solenoidal vector fields cannot hold discretely because of grid oscillations,which are characterised explicitly.This results in a non-vanishing remainder associated with grid oscillations in the discrete Helmholtz Hodge decomposition.Nevertheless,iterative numerical methods based on an interpretation of the Helmholtz Hodge decomposition via orthogonal projections are pro-posed and applied successfully.In numerical experiments,the discrete remainder vanishes and the potentials converge with the same order of accuracy as usual in other first-order partial differential equations.Motivated by the successful application of the Helmholtz Hodge decomposition in theoretical plasma physics,applications to the discrete analysis of magnetohydrodynamic(MHD) wave modes are presented and discussed.展开更多
We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusi...We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusing on discontinuous spectral element semidis-cretizations,we design new controllers for existing methods and for some new embedded Runge-Kutta pairs.We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice.We compare a wide range of error-control-based methods,along with the common approach in which step size con-trol is based on the Courant-Friedrichs-Lewy(CFL)number.The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances,while additionally providing control of the temporal error at tighter tolerances.The numerical examples include challenging industrial CFD applications.展开更多
文摘In this article,discrete variants of several results from vector calculus are studied for clas-sical finite difference summation by parts operators in two and three space dimensions.It is shown that existence theorems for scalar/vector potentials of irrotational/solenoidal vector fields cannot hold discretely because of grid oscillations,which are characterised explicitly.This results in a non-vanishing remainder associated with grid oscillations in the discrete Helmholtz Hodge decomposition.Nevertheless,iterative numerical methods based on an interpretation of the Helmholtz Hodge decomposition via orthogonal projections are pro-posed and applied successfully.In numerical experiments,the discrete remainder vanishes and the potentials converge with the same order of accuracy as usual in other first-order partial differential equations.Motivated by the successful application of the Helmholtz Hodge decomposition in theoretical plasma physics,applications to the discrete analysis of magnetohydrodynamic(MHD) wave modes are presented and discussed.
基金Open Access funding enabled and organized by Projekt DEAL.
文摘We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusing on discontinuous spectral element semidis-cretizations,we design new controllers for existing methods and for some new embedded Runge-Kutta pairs.We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice.We compare a wide range of error-control-based methods,along with the common approach in which step size con-trol is based on the Courant-Friedrichs-Lewy(CFL)number.The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances,while additionally providing control of the temporal error at tighter tolerances.The numerical examples include challenging industrial CFD applications.
