Lemmy: Bestiverse
  • Communities
  • Create Post
  • Create Community
  • heart
    Support Lemmy
  • search
    Search
  • Login
  • Sign Up
RSS BotMB to Hacker NewsEnglish · 13 days ago

Non-uniform finite-element meshes defined by ray dynamics for Helmholtz problems

arxiv.org

external-link
message-square
0
fedilink
4
external-link

Non-uniform finite-element meshes defined by ray dynamics for Helmholtz problems

arxiv.org

RSS BotMB to Hacker NewsEnglish · 13 days ago
message-square
0
fedilink
The $h$-version of the finite-element method ($h$-FEM) applied to the high-frequency Helmholtz equation has been a classic topic in numerical analysis since the 1990s. It is now rigorously understood that (using piecewise polynomials of degree $p$ on a mesh of a maximal width $h$) the conditions "$(hk)^p ρ$ sufficiently small" and "$(hk)^{2p} ρ$ sufficiently small" guarantee, respectively, $k$-uniform quasioptimality (QO) and bounded relative error (BRE), where $ρ$ is the norm of the solution operator with $ρ\sim k$ for non-trapping problems. Empirically, these conditions are observed to be optimal in the context of $h$-FEM with a uniform mesh. This paper demonstrates that QO and BRE can be achieved using certain non-uniform meshes that violate the conditions above on $h$ and involve coarser meshes away from trapping and in the perfectly matched layer (PML). The main theorem details how varying the meshwidth in one region affects errors both in that region and elsewhere. One notable consequence is that, for any scattering problem (trapping or nontrapping), in the PML one only needs $hk$ to be sufficiently small; i.e. there is no pollution in the PML. The motivating idea for the analysis is that the Helmholtz data-to-solution map behaves differently depending on the locations of both the measurement and data, in particular, on the properties of billiards trajectories (i.e. rays) through these sets. Because of this, it is natural that the approximation requirements for finite-element spaces in a subset should depend on the properties of billiard rays through that set. Inserting this behaviour into the latest duality arguments for the FEM applied to the high-frequency Helmholtz equation allows us to retain detailed information about the influence of $\textit{both}$ the mesh structure $\textit{and}$ the behaviour of the true solution on local errors in FEM.

Comments

alert-triangle
You must log in or register to comment.

Hacker News

hackernews

Subscribe from Remote Instance

You are not logged in. However you can subscribe from another Fediverse account, for example Lemmy or Mastodon. To do this, paste the following into the search field of your instance: !hackernews@lemmy.bestiver.se
lock
Community locked: only moderators can create posts. You can still comment on posts.

Posts from the RSS Feed of HackerNews.

The feed sometimes contains ads and posts that have been removed by the mod team at HN.

Visibility: Public
globe

This community can be federated to other instances and be posted/commented in by their users.

  • 381 users / day
  • 1.64K users / week
  • 3.45K users / month
  • 9K users / 6 months
  • 2 local subscribers
  • 1.94K subscribers
  • 25.3K Posts
  • 9.17K Comments
  • Modlog
  • mods:
  • patrick
  • RSS Bot
  • BE: 0.19.5
  • Modlog
  • Instances
  • Docs
  • Code
  • join-lemmy.org