39,524 research outputs found
Positive trigonometric polynomials for strong stability of difference equations
We follow a polynomial approach to analyse strong stability of linear
difference equations with rationally independent delays. Upon application of
the Hermite stability criterion on the discrete-time homogeneous characteristic
polynomial, assessing strong stability amounts to deciding positive
definiteness of a multivariate trigonometric polynomial matrix. This latter
problem is addressed with a converging hierarchy of linear matrix inequalities
(LMIs). Numerical experiments indicate that certificates of strong stability
can be obtained at a reasonable computational cost for state dimension and
number of delays not exceeding 4 or 5
Robust Stability Analysis of Nonlinear Hybrid Systems
We present a methodology for robust stability analysis of nonlinear hybrid systems, through the algorithmic construction of polynomial and piecewise polynomial Lyapunov-like functions using convex optimization and in particular the sum of squares decomposition of multivariate polynomials. Several improvements compared to previous approaches are discussed, such as treating in a unified way polynomial switching surfaces and robust stability analysis for nonlinear hybrid systems
Three-dimensional central-moments-based lattice Boltzmann method with external forcing: A consistent, concise and universal formulation
The cascaded or central-moments-based lattice Boltzmann method (CM-LBM) is a
robust alternative to the more conventional BGK-LBM for the simulation of
high-Reynolds number flows. Unfortunately, its original formulation makes its
extension to a broader range of physics quite difficult. To tackle this issue,
a recent work [A. De Rosis, Phys. Rev. E 95, 013310 (2017)] proposed a more
generic way to derive concise and efficient three-dimensional CM-LBMs. Knowing
the original model also relies on central moments that are derived in an adhoc
manner, i.e., by mimicking those of the Maxwell-Boltzmann distribution to
ensure their Galilean invariance a posteriori, a very recent effort [A. De
Rosis and K. H. Luo, Phys. Rev. E 99, 013301 (2019)] was proposed to further
generalize their derivation. The latter has shown that one could derive
Galilean invariant CMs in a systematic and a priori manner by taking into
account high-order Hermite polynomials in the derivation of the discrete
equilibrium state. Combining these two approaches, a compact and mathematically
sound formulation of the CM-LBM with external forcing is proposed. More
specifically, the proposed formalism fully takes advantage of the D3Q27
discretization by relying on the corresponding set of 27 Hermite polynomials
(up to the sixth order) for the derivation of both the discrete equilibrium
state and the forcing term. The present methodology is more consistent than
previous approaches, as it properly explains how to derive Galilean invariant
CMs of the forcing term in an a priori manner. Furthermore, while keeping the
numerical properties of the original CM-LBM, the present work leads to a
compact and simple algorithm, representing a universal methodology based on CMs
and external forcing within the lattice Boltzmann framework.Comment: Published in Phys. Fluids as Editor's Pic
Stability of quaternionic linear systems
The main goal of this paper is to characterize stability and bounded-input-bounded-output (BIBO)-stability of quaternionic dynamical systems. After defining the quaternion skew-field, algebraic properties of quaternionic polynomials such as divisibility and coprimeness are investigated. Having established these results, the Smith and the Smith-McMillan forms of quaternionic matrices are introduced and studied. Finally, all the tools that were developed are used to analyze stability of quaternionic linear systems in a behavioral framework
Multi-patch discontinuous Galerkin isogeometric analysis for wave propagation: explicit time-stepping and efficient mass matrix inversion
We present a class of spline finite element methods for time-domain wave
propagation which are particularly amenable to explicit time-stepping. The
proposed methods utilize a discontinuous Galerkin discretization to enforce
continuity of the solution field across geometric patches in a multi-patch
setting, which yields a mass matrix with convenient block diagonal structure.
Over each patch, we show how to accurately and efficiently invert mass matrices
in the presence of curved geometries by using a weight-adjusted approximation
of the mass matrix inverse. This approximation restores a tensor product
structure while retaining provable high order accuracy and semi-discrete energy
stability. We also estimate the maximum stable timestep for spline-based finite
elements and show that the use of spline spaces result in less stringent CFL
restrictions than equivalent piecewise continuous or discontinuous finite
element spaces. Finally, we explore the use of optimal knot vectors based on L2
n-widths. We show how the use of optimal knot vectors can improve both
approximation properties and the maximum stable timestep, and present a simple
heuristic method for approximating optimal knot positions. Numerical
experiments confirm the accuracy and stability of the proposed methods
- âŠ