3,879 research outputs found
An Improved High Order Finite Difference Method for Non-conforming Grid Interfaces for the Wave Equation
This paper presents an extension of a recently developed high order finite difference method for the wave equation on a grid with non-conforming interfaces. The stability proof of the existing methods relies on the interpolation operators being norm-contracting, which is satisfied by the second and fourth order operators, but not by the sixth order operator. We construct new penalty terms to impose interface conditions such that the stability proof does not require the norm-contracting condition. As a consequence, the sixth order accurate scheme is also provably stable. Numerical experiments demonstrate the improved stability and accuracy property
An Improved High Order Finite Difference Method for Non-conforming Grid Interfaces for the Wave Equation
This paper presents an extension of a recently developed high order finite difference method for the wave equation on a grid with non-conforming interfaces. The stability proof of the existing methods relies on the interpolation operators being norm-contracting, which is satisfied by the second and fourth order operators, but not by the sixth order operator. We construct new penalty terms to impose interface conditions such that the stability proof does not require the norm-contracting condition. As a consequence, the sixth order accurate scheme is also provably stable. Numerical experiments demonstrate the improved stability and accuracy property
High-order numerical methods for 2D parabolic problems in single and composite domains
In this work, we discuss and compare three methods for the numerical
approximation of constant- and variable-coefficient diffusion equations in both
single and composite domains with possible discontinuity in the solution/flux
at interfaces, considering (i) the Cut Finite Element Method; (ii) the
Difference Potentials Method; and (iii) the summation-by-parts Finite
Difference Method. First we give a brief introduction for each of the three
methods. Next, we propose benchmark problems, and consider numerical tests-with
respect to accuracy and convergence-for linear parabolic problems on a single
domain, and continue with similar tests for linear parabolic problems on a
composite domain (with the interface defined either explicitly or implicitly).
Lastly, a comparative discussion of the methods and numerical results will be
given.Comment: 45 pages, 12 figures, in revision for Journal of Scientific Computin
A Comparison of Two Shallow Water Models with Non-Conforming Adaptive Grids: classical tests
In an effort to study the applicability of adaptive mesh refinement (AMR)
techniques to atmospheric models an interpolation-based spectral element
shallow water model on a cubed-sphere grid is compared to a block-structured
finite volume method in latitude-longitude geometry. Both models utilize a
non-conforming adaptation approach which doubles the resolution at fine-coarse
mesh interfaces. The underlying AMR libraries are quad-tree based and ensure
that neighboring regions can only differ by one refinement level.
The models are compared via selected test cases from a standard test suite
for the shallow water equations. They include the advection of a cosine bell, a
steady-state geostrophic flow, a flow over an idealized mountain and a
Rossby-Haurwitz wave. Both static and dynamics adaptations are evaluated which
reveal the strengths and weaknesses of the AMR techniques. Overall, the AMR
simulations show that both models successfully place static and dynamic
adaptations in local regions without requiring a fine grid in the global
domain. The adaptive grids reliably track features of interests without visible
distortions or noise at mesh interfaces. Simple threshold adaptation criteria
for the geopotential height and the relative vorticity are assessed.Comment: 25 pages, 11 figures, preprin
- …