597 research outputs found
Review of Summation-by-parts schemes for initial-boundary-value problems
High-order finite difference methods are efficient, easy to program, scales
well in multiple dimensions and can be modified locally for various reasons
(such as shock treatment for example). The main drawback have been the
complicated and sometimes even mysterious stability treatment at boundaries and
interfaces required for a stable scheme. The research on summation-by-parts
operators and weak boundary conditions during the last 20 years have removed
this drawback and now reached a mature state. It is now possible to construct
stable and high order accurate multi-block finite difference schemes in a
systematic building-block-like manner. In this paper we will review this
development, point out the main contributions and speculate about the next
lines of research in this area
Spectral/hp element methods: recent developments, applications, and perspectives
The spectral/hp element method combines the geometric flexibility of the
classical h-type finite element technique with the desirable numerical
properties of spectral methods, employing high-degree piecewise polynomial
basis functions on coarse finite element-type meshes. The spatial approximation
is based upon orthogonal polynomials, such as Legendre or Chebychev
polynomials, modified to accommodate C0-continuous expansions. Computationally
and theoretically, by increasing the polynomial order p, high-precision
solutions and fast convergence can be obtained and, in particular, under
certain regularity assumptions an exponential reduction in approximation error
between numerical and exact solutions can be achieved. This method has now been
applied in many simulation studies of both fundamental and practical
engineering flows. This paper briefly describes the formulation of the
spectral/hp element method and provides an overview of its application to
computational fluid dynamics. In particular, it focuses on the use the
spectral/hp element method in transitional flows and ocean engineering.
Finally, some of the major challenges to be overcome in order to use the
spectral/hp element method in more complex science and engineering applications
are discussed
An improved local radial basis function method for solving small-strain elasto-plasticity
Strong-form meshless methods received much attention in recent years and are
being extensively researched and applied to a wide range of problems in science
and engineering. However, the solution of elasto-plastic problems has proven to
be elusive because of often non-smooth constitutive relations between stress
and strain. The novelty in tackling them is the introduction of virtual finite
difference stencils to formulate a hybrid radial basis function generated
finite difference (RBF-FD) method, which is used to solve smallstrain von Mises
elasto-plasticity for the first time by this original approach. The paper
further contrasts the new method to two alternative legacy RBF-FD approaches,
which fail when applied to this class of problems. The three approaches differ
in the discretization of the divergence operator found in the balance equation
that acts on the non-smooth stress field. Additionally, an innovative
stabilization technique is employed to stabilize boundary conditions and is
shown to be essential for any of the approaches to converge successfully.
Approaches are assessed on elastic and elasto-plastic benchmarks where
admissible ranges of newly introduced free parameters are studied regarding
stability, accuracy, and convergence rate
Differential-Algebraic Equations and Beyond: From Smooth to Nonsmooth Constrained Dynamical Systems
The present article presents a summarizing view at differential-algebraic
equations (DAEs) and analyzes how new application fields and corresponding
mathematical models lead to innovations both in theory and in numerical
analysis for this problem class. Recent numerical methods for nonsmooth
dynamical systems subject to unilateral contact and friction illustrate the
topicality of this development.Comment: Preprint of Book Chapte
The role of numerical boundary procedures in the stability of perfectly matched layers
In this paper we address the temporal energy growth associated with numerical
approximations of the perfectly matched layer (PML) for Maxwell's equations in
first order form. In the literature, several studies have shown that a
numerical method which is stable in the absence of the PML can become unstable
when the PML is introduced. We demonstrate in this paper that this instability
can be directly related to numerical treatment of boundary conditions in the
PML. First, at the continuous level, we establish the stability of the constant
coefficient initial boundary value problem for the PML. To enable the
construction of stable numerical boundary procedures, we derive energy
estimates for the variable coefficient PML. Second, we develop a high order
accurate and stable numerical approximation for the PML using
summation--by--parts finite difference operators to approximate spatial
derivatives and weak enforcement of boundary conditions using penalties. By
constructing analogous discrete energy estimates we show discrete stability and
convergence of the numerical method. Numerical experiments verify the
theoretical result
- …