An official website of the United States government
Abstract This work focuses on the development of a novel, strongly-coupled, second-order partitioned method for fluid–poroelastic structure interaction. The flow is assumed to be viscous and incompressible, and the poroelastic material is described using the Biot model. To solve this problem, a numerical method is proposed, based on Robin interface conditions combined with the refactorization of the Cauchy’s one-legged ‘ϑ-like’ method. This approach allows the use of the mixed formulation for the Biot model. The proposed algorithm consists of solving a sequence of Backward Euler–Forward Euler steps. In the Backward Euler step, the fluid and poroelastic structure problems are solved iteratively until convergence. Then, the Forward Euler problems are solved using equivalent linear extrapolations. We prove that the iterative procedure in the Backward Euler step is convergent, and that the converged method is stable whenϑ∈ [1/2, 1]. Numerical examples are used to explore convergence rates with varying parameters used in our scheme, and to compare our method to a monolithic method based on Nitsche’s coupling approach.
more »
« less
Simultaneous mode, input and state estimation for switched linear stochastic systems
Summary In this paper, we propose a filtering algorithm for simultaneously estimating the mode, input and state of hidden mode switched linear stochastic systems with unknown inputs. Using a multiple‐model approach with a bank of linear input and state filters for each mode, our algorithm relies on the ability to find the most probable model as a mode estimate, which we show is possible with input and state filters by identifying a key property, that a particular residual signal we callgeneralized innovationis a Gaussian white noise. We also provide an asymptotic analysis for the proposed algorithm and provide sufficient conditions forasymptoticallyachieving convergence to the true model (consistency), or to the “closest” model according to an information‐theoretic measure (convergence). A simulation example of intention‐aware vehicles at an intersection is given to demonstrate the effectiveness of our approach.
more »
« less
- PAR ID:
- 10455032
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal of Robust and Nonlinear Control
- Volume:
- 31
- Issue:
- 2
- ISSN:
- 1049-8923
- Page Range / eLocation ID:
- p. 640-661
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Finite-time Sample Complexity Analysis of Least Square Identifying Stochastic Switched Linear SystemIn this paper, we examine the high-probability finite-time theoretical guarantees of the least squares method for system identification of switched linear systems with process noise and without control input. We consider two scenarios: one in which the switching is i.i.d., and the other in which the switching is according to a Markov process. We provide concentration inequalities using a martingale-type argument to bound the identification error at each mode, and we use concentration lemmas for the switching signal. Our bound is in terms of state dimension, trajectory length, finite-time gramian, and properties of the switching signal distribution. We then provide simulations to demonstrate the accuracy of the identification. Additionally, we show that the empirical convergence rate is consistent with our theoretical bound.more » « less
-
Abstract We study the problem of fitting a piecewise affine (PWA) function to input–output data. Our algorithm divides the input domain into finitely many regions whose shapes are specified by a user-provided template and such that the input–output data in each region are fit by an affine function within a user-provided error tolerance. We first prove that this problem is NP-hard. Then, we present a top-down algorithmic approach for solving the problem. The algorithm considers subsets of the data points in a systematic manner, trying to fit an affine function for each subset using linear regression. If regression fails on a subset, the algorithm extracts a minimal set of points from the subset (an unsatisfiable core) that is responsible for the failure. The identified core is then used to split the current subset into smaller ones. By combining this top-down scheme with a set-covering algorithm, we derive an overall approach that provides optimal PWA models for a given error tolerance, where optimality refers to minimizing the number of pieces of the PWA model. We demonstrate our approach on three numerical examples that include PWA approximations of a widely used nonlinear insulin–glucose regulation model and a double inverted pendulum with soft contacts.more » « less
-
In this paper, we develop a distributed consensus algorithm for agents whose states evolve on a manifold. This algorithm is complementary to traditional consensus, predominantly developed for systems with dynamics on vector spaces. We provide theoretical convergence guarantees for the proposed manifold consensus provided that agents are initialized within a geodesically convex (g-convex) set. This required condition on initialization is not restrictive as g-convex sets may be comparatively “large” for relevant Riemannian manifolds. Our approach to manifold consensus builds upon the notion of Riemannian Center of Mass (RCM) and the intrinsic structure of the manifold to avoid projections in the ambient space. We first show that on a g-convex ball, all states coincide if and only if each agent’s state is the RCM of its neighbors’ states. This observation facilitates our convergence guarantee to the consensus submanifold. Finally, we provide simulation results that exemplify the linear convergence rate of the proposed algorithm and illustrates its statistical properties over randomly generated problem instances.more » « less
-
Abstract This paper studies the linear convergence of the subspace constrained mean shift (SCMS) algorithm, a well-known algorithm for identifying a density ridge defined by a kernel density estimator. By arguing that the SCMS algorithm is a special variant of a subspace constrained gradient ascent (SCGA) algorithm with an adaptive step size, we derive the linear convergence of such SCGA algorithm. While the existing research focuses mainly on density ridges in the Euclidean space, we generalize density ridges and the SCMS algorithm to directional data. In particular, we establish the stability theorem of density ridges with directional data and prove the linear convergence of our proposed directional SCMS algorithm.more » « less
