A modified Legendre-Gauss-Radau collocation method is developed for solving optimal control problems whose solutions contain a nonsmooth optimal control. The method includes an additional variable that defines the location of nonsmoothness. In addition, collocation constraints are added at the end of a mesh interval that defines the location of nonsmoothness in the solution on each differential equation that is a function of control along with a control constraint at the endpoint of this same mesh interval. The transformed adjoint system for the modified Legendre-Gauss-Radau collocation method along with a relationship between the Lagrange multipliers of the nonlinear programming problem and a discrete approximation of the costate of the optimal control problem is then derived. Finally, it is shown via example that the new method provides an accurate approximation of the costate.
more »
« less
This content will become publicly available on January 15, 2026
A new nonsmooth optimal control framework for wind turbine power systems
Optimal control theory extending from the calculus of variations has not been used to study the wind turbine power system (WTPS) control problem, which aims at achieving two targets: (i) maximizing power generation in lower wind speed conditions; and (ii) maintaining the output power at the rated level in high wind speed conditions. A lack of an optimal control framework for the WTPS (i.e., no access to actual optimal control trajectories) reduces optimal control design potential and prevents competing control methods of WTPSs to have a reference control solution for comparison. In fact, the WTPS control literature often relies on reduced and linearized models of WTPSs, and avoids the nonsmoothness present in the system during transitions between different conditions of operation. In this paper, we introduce a novel optimal control framework for the WTPS control problem. We use in our formulation a recent accurate, nonlinear differential–algebraic equation (DAE) model of WTPSs, which we then generalize over all wind speed ranges using nonsmooth functions. We also use developments in nonsmooth optimal control theory to take into account nonsmoothness present in the system. We implement this new WTPS optimal control approach to solve the problem numerically, including (i) different wind speed profiles for testing the system response; (ii) real-world wind data; and (iii) a comparison with smoothing and naive approaches. Results show the effectiveness of the proposed approach.
more »
« less
- PAR ID:
- 10582560
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- Journal of the Franklin Institute
- Volume:
- 362
- Issue:
- 3
- ISSN:
- 0016-0032
- Page Range / eLocation ID:
- 107498
- Subject(s) / Keyword(s):
- Wind turbine Optimal control Nonsmooth systems Power systems Generalized derivatives
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)A new method is developed for solving optimal control problems whose solutions are nonsmooth. The method developed in this paper employs a modified form of the Legendre–Gauss–Radau orthogonal direct collocation method. This modified Legendre–Gauss–Radau method adds two variables and two constraints at the end of a mesh interval when compared with a previously developed standard Legendre– Gauss–Radau collocation method. The two additional variables are the time at the interface between two mesh intervals and the control at the end of each mesh inter- val. The two additional constraints are a collocation condition for those differential equations that depend upon the control and an inequality constraint on the control at the endpoint of each mesh interval. The additional constraints modify the search space of the nonlinear programming problem such that an accurate approximation to the location of the nonsmoothness is obtained. The transformed adjoint system of the modified Legendre–Gauss–Radau method is then developed. Using this transformed adjoint system, a method is developed to transform the Lagrange multipliers of the nonlinear programming problem to the costate of the optimal control problem. Fur- thermore, it is shown that the costate estimate satisfies one of the Weierstrass–Erdmann optimality conditions. Finally, the method developed in this paper is demonstrated on an example whose solution is nonsmooth.more » « less
-
A new method is developed for solving optimal control problems whose solutions contain a nonsmooth optimal control. The method developed in this paper employs a modified form of the Legendre-Gauss-Radau (LGR) orthogonal direct collocation method in which an additional variable and two additional constraints are included at the end of a mesh interval. The additional variable is the switch time where a discontinuity occurs. The two additional constraints are a collocation condition on each differential equation that is a function of control along with a control constraint at the endpoint of the mesh interval that defines the location of the nonsmoothness. These additional constraints modify the search space of the NLP in a manner such that an accurate approximation to the location of the nonsmoothness is obtained. An example with a nonsmooth solution is used throughout the paper to illustrate the improvement of the method over the standard Legendre-Gauss-Radau collocation method.more » « less
-
The data-driven newsvendor problem with features has recently emerged as a significant area of research, driven by the proliferation of data across various sectors such as retail, supply chains, e-commerce, and healthcare. Given the sensitive nature of customer or organizational data often used in feature-based analysis, it is crucial to ensure individual privacy to uphold trust and confidence. Despite its importance, privacy preservation in the context of inventory planning remains unexplored. A key challenge is the nonsmoothness of the newsvendor loss function, which sets it apart from existing work on privacy-preserving algorithms in other settings. This paper introduces a novel approach to estimating a privacy-preserving optimal inventory policy within the f-differential privacy framework, an extension of the classical [Formula: see text]-differential privacy with several appealing properties. We develop a clipped noisy gradient descent algorithm based on convolution smoothing for optimal inventory estimation to simultaneously address three main challenges: (i) unknown demand distribution and nonsmooth loss function, (ii) provable privacy guarantees for individual-level data, and (iii) desirable statistical precision. We derive finite-sample high-probability bounds for optimal policy parameter estimation and regret analysis. By leveraging the structure of the newsvendor problem, we attain a faster excess population risk bound compared with that obtained from an indiscriminate application of existing results for general nonsmooth convex loss. Our bound aligns with that for strongly convex and smooth loss function. Our numerical experiments demonstrate that the proposed new method can achieve desirable privacy protection with a marginal increase in cost. This paper was accepted by J. George Shanthikumar, data science. Funding: This work was supported by the National Science Foundation [Grants DMS-2113409 and DMS 2401268 to W.-X. Zhou, and FRGMS-1952373 to L. Wang]. Supplemental Material: The online appendix and data files are available at https://doi.org/10.1287/mnsc.2023.01268 .more » « less
-
In this article, we present a constraint-driven optimal control framework that achieves emergent cluster flocking within a constrained 2D environment. We formulate a decentralized optimal control problem that includes safety, flocking, and predator avoidance constraints. We explicitly derive conditions for constraint compatibility and propose an event-driven constraint relaxation scheme. We map this to an equivalent switching system that intuitively describes the behavior of each agent in the system. Instead of minimizing control effort, as it is common in the ecologically-inspired robotics literature, in our approach, we minimize each agent’s deviation from their most efficient locomotion speed. Finally, we demonstrate our approach in simulation both with and without the presence of a predator.more » « less