skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: An algorithm for computing robust forward invariant sets of two dimensional nonlinear systems
Summary Robustness of nonlinear systems can be analyzed by computing robust forward invariant sets (RFISs). Knowledge of the smallest RFIS of a system, can help analyze system performance under perturbations. A novel algorithm is developed to compute an approximation of the smallest RFIS for two‐dimensional nonlinear systems subjected to a bounded additive disturbance. The problem of computing an RFIS is formulated as a path planning problem, and the algorithm developed plans a path which iteratively converges to the boundary of an RFIS. Rigorous mathematical analysis shows that the proposed algorithm terminates in a finite number of iterations, and that the output of the proposed algorithm is an RFIS. Simulations are presented to illustrate the proposed algorithm, and to support the mathematical results. This work may aid future development, for use with higher dimensional systems.  more » « less
Award ID(s):
1849228 1828678
PAR ID:
10450067
Author(s) / Creator(s):
 ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Asian Journal of Control
Volume:
23
Issue:
5
ISSN:
1561-8625
Page Range / eLocation ID:
p. 2403-2419
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract A technique based on the Wiener path integral (WPI) is developed for determining the stochastic response of diverse nonlinear systems with fractional derivative elements. Specifically, a reduced-order WPI formulation is proposed, which can be construed as an approximation-free dimension reduction approach that renders the associated computational cost independent of the total number of stochastic dimensions of the problem. In fact, the herein developed technique can determine, directly, any lower-dimensional joint response probability density function corresponding to a subset only of the response vector components. This is done by utilizing an appropriate combination of fixed and free boundary conditions in the related variational, functional minimization, problem. Notably, the reduced-order WPI formulation is particularly advantageous for problems where the interest lies in few only specific degrees-of-freedom whose stochastic response is critical for the design and optimization of the overall system. An indicative numerical example is considered pertaining to a stochastically excited tuned mass-damper-inerter nonlinear system with a fractional derivative element. Comparisons with relevant Monte Carlo simulation data demonstrate the accuracy and computational efficiency of the technique. 
    more » « less
  2. Abstract We propose easy‐to‐implement heuristics for time‐constrained applications of a problem referred to in the literature as the facility location problem with immobile servers, stochastic demand, and congestion, the service system design problem, or the immobile server problem (ISP). The problem is typically posed as one of allocating capacity to a set of M/M/1 queues to which customers with stochastic demand are assigned with the objective of minimizing a cost function composed of a fixed capacity‐acquisition cost, a variable customer‐assignment cost, and an expected‐waiting‐time cost. The expected‐waiting‐time cost results in a nonlinear term in the objective function of the standard binary programming formulation of the problem. Thus, the solution approaches proposed in the literature are either sophisticated linearization or relaxation schemes, or metaheuristics. In this study, we demonstrate that an ensemble of straightforward, greedy heuristics can rapidly find high‐quality solutions. In addition to filling a gap in the literature on ISP heuristics, new stopping criteria for an existing cutting plane algorithm are proposed and tested, and a new mixed‐integer linear model requiring no iterating algorithm is developed. In many cases, our heuristic approach finds solutions of the same or better quality than those found by exact methods implemented with expensive, state‐of‐the‐art mathematical programming software, in particular a commercial nonlinear mixed‐integer linear programming solver, given a five‐minute time limit. 
    more » « less
  3. Abstract We study a second order ensemble method for fast computation of an ensemble of magnetohydrodynamics flows at small magnetic Reynolds number. Computing an ensemble of flow equations with different input parameters is a common procedure for uncertainty quantification in many engineering applications, for which the computational cost can be prohibitively expensive for nonlinear complex systems. We propose an ensemble algorithm that requires only solving one linear system with multiple right‐hands instead of solving multiple different linear systems, which significantly reduces the computational cost and simulation time. Comprehensive stability and error analyses are presented proving conditional stability and second order in time convergent. Numerical tests are provided to illustrate theoretical results and demonstrate the efficiency of the proposed algorithm. 
    more » « less
  4. Abstract This article introduces a computational design framework for obtaining three‐dimensional (3D) periodic elastoplastic architected materials with enhanced performance, subject to uniaxial or shear strain. A nonlinear finite element model accounting for plastic deformation is developed, where a Lagrange multiplier approach is utilized to impose periodicity constraints. The analysis assumes that the material obeys a von Mises plasticity model with linear isotropic hardening. The finite element model is combined with a corresponding path‐dependent adjoint sensitivity formulation, which is derived analytically. The optimization problem is parametrized using the solid isotropic material penalization method. Designs are optimized for either end compliance or toughness for a given prescribed displacement. Such a framework results in producing materials with enhanced performance through much better utilization of an elastoplastic material. Several 3D examples are used to demonstrate the effectiveness of the mathematical framework. 
    more » « less
  5. This paper considers the single source shortest path (SSSP) problem, which is the key for many applications such as navigation, mapping, routing, and social networking. Existing SSSP algorithms are designed mostly for shared-memory systems. Nevertheless, with the prevalence of diverse smart devices like drones, there is a growing interest in deploying SSSP algorithms over distributed computing systems so that they can run efficiently onboard of smart devices via Mobile Ad Hoc Computing or at the network edges via Mobile Edge Computing. In this paper, we introduce a communication-efficient ∆-stepping algorithm for distributed computing systems. The proposed algorithm is featured by 1) a message coordination architecture for reducing message exchanges between workers, 2) a pruning technique for reducing redundant computations, and 3) an aggregation technique for further reducing message exchanges when communication delay is significant. Theoretical analyses and experimental studies on real-world graph datasets demonstrate the promising performance of proposed algorithm. 
    more » « less