On the analysis of linear and nonlinear dynamical systems from transient-response data
Explore the source record for details and available documents.
Shinbrot, Marvin↗
SEARCH · Engineering Papers
Results for “Linear dynamics”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Is Large-Scale Low-Frequency Variability Satisfying Linear Barotropic Vorticity Dynamics?
What is the dynamics governing the sea level variabilities at spatial scales larger than the mesoscale and periods longer than a few weeks? This question is investigated using the TOPEX/POSEIDON altimetry data. For periods shorter than a decade, it is believed that the large-scale variabilities in the open ocean should be barotropic and can be described by the barotropic vorticity equation. To avoid the error-sensitive Laplacian operator, the above equation was integrated over a large area. Preliminary analysis was performed in the northeastern Pacific, where the eddy energy is relatively low and the bathymetry is relatively smooth. Good correlation between the two sides of the equation was obtained at periods longer than 60 days. The largest error in the data is suspected to be the ocean tides. Empirical correction for the ocean tides will be performed for further analysis. Preliminary results of a global calculation will be presented.
Fu, Lee-Lueng↗
Controllability of higher order linear systems.
Dynamical system with time evolution described by second order linear differential equation in complex Banach space, discussing controllability
Fattorini, H. O.↗
Fixed order dynamic compensation for multivariable linear systems
This paper considers the design of fixed order dynamic compensators for multivariable time invariant linear systems, minimizing a linear quadratic performance cost functional. Attention is given to robustness issues in terms of multivariable frequency domain specifications. An output feedback formulation is adopted by suitably augmenting the system description to include the compensator states. Either a controller or observer canonical form is imposed on the compensator description to reduce the number of free parameters to its minimal number. The internal structure of the compensator is prespecified by assigning a set of ascending feedback invariant indices, thus forming a Brunovsky structure for the nominal compensator.
Kramer, F. S.↗
Classical linear-control analysis applied to business-cycle dynamics and stability
Linear control analysis is applied as an aid in understanding the fluctuations of business cycles in the past, and to examine monetary policies that might improve stabilization. The analysis shows how different policies change the frequency and damping of the economic system dynamics, and how they modify the amplitude of the fluctuations that are caused by random disturbances. Examples are used to show how policy feedbacks and policy lags can be incorporated, and how different monetary strategies for stabilization can be analytically compared. Representative numerical results are used to illustrate the main points.
Wingrove, R. C.↗
Determination of dynamic processes in piecewise linear automatic control systems <opredeleniye dinamicheskikm protsessov v kusochno-lineynykh sistemakh avtomaticheskogo regulirovaniya<
Dynamic processes in piecewise linear automatic control systems
AUTOMATIC CONTROL↗
A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems
This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.
97 MATHEMATICS AND COMPUTING↗
Streaming Instability with Multiple Dust Species – II. Turbulence and Dust–Gas Dynamics at Non-Linear Saturation
The streaming instability is a fundamental process that can drive dust–gas dynamics and ultimately planetesimal formation inprotoplanetary discs. As a linear instability, it has been shown that its growth with a distribution of dust sizes can be classifiedinto two distinct regimes, fast- and slow-growth, depending on the dust-size distribution and the total dust-to-gas density ratio. Using numerical simulations of an unstratified disc, we bring three cases in different regimes into non-linear saturation. Wefind that the saturation states of the two fast-growth cases are similar to its single-species counterparts. The one with maximumdimensionless stopping timeτs,max=0.1 and=2 drives turbulent vertical dust–gas vortices, while the other withτs,max=2and=0.2 leads to radial traffic jams and filamentary structures of dust particles. The dust density distribution for the former isflat in low densities, while the one for the latter has a low-end cut-off. By contrast, the one slow-growth case results in a virtuallyquiescent state. Moreover, we find that in the fast-growth regime, significant dust segregation by size occurs, with large particlesmoving towards dense regions while small particles remain in the diffuse regions, and the mean radial drift of each dust speciesis appreciably altered from the (initial) drag-force equilibrium. The former effect may skew the spectral index derived frommultiwavelength observations and change the initial size distribution of a pebble cloud for planetesimal formation. The latteralong with turbulent diffusion may influence the radial transport and mixing of solid materials in young protoplanetary discs.
Chao-chin Yang↗
Quantum Algorithm for Linear Non-unitary Dynamics with Near-Optimal Dependence on All Parameters
We introduce a family of identities that express general linear non-unitary evolution operators as a linear combination of unitary evolution operators, each solving a Hamiltonian simulation problem. This formulation can exponentially enhance the accuracy of the recently introduced linear combination of Hamiltonian simulation (LCHS) method [An, Liu, and Lin, Physical Review Letters, 2023]. For the first time, this approach enables quantum algorithms to solve linear differential equations with both optimal state preparation cost and near-optimal scaling in matrix queries on all parameters.
Applied Dynamical Systems↗
Non-linear controls influence functions in an aircraft dynamics simulator
In the development and testing of novel structural and controls concepts, such as morphing aircraft wings, appropriate models are needed for proper system characterization. In most instances, available system models do not provide the required additional degrees of freedom for morphing structures but may be modified to some extent to achieve a compatible system. The objective of this study is to apply wind tunnel data collected for an Unmanned Air Vehicle (UAV), that implements trailing edge morphing, to create a non-linear dynamics simulator, using well defined rigid body equations of motion, where the aircraft stability derivatives change with control deflection. An analysis of this wind tunnel data, using data extraction algorithms, was performed to determine the reference aerodynamic force and moment coefficients for the aircraft. Further, non-linear influence functions were obtained for each of the aircraft s control surfaces, including the sixteen trailing edge flap segments. These non-linear controls influence functions are applied to the aircraft dynamics to produce deflection-dependent aircraft stability derivatives in a non-linear dynamics simulator. Time domain analysis of the aircraft motion, trajectory, and state histories can be performed using these nonlinear dynamics and may be visualized using a 3-dimensional aircraft model. Linear system models can be extracted to facilitate frequency domain analysis of the system and for control law development. The results of this study are useful in similar projects where trailing edge morphing is employed and will be instrumental in the University of Maryland s continuing study of active wing load control.
Guerreiro, Nelson M.↗
Human response in closed loop control of dynamic systems
Linear transfer function for describing human response to aircraft control
Adams, J. J.↗
Essential uncontrollability of discrete linear, time-invariant, dynamical systems
The concept of a 'best approximating m-dimensional subspace' for a given set of vectors in n-dimensional whole space is introduced. Such a subspace is easily described in terms of the eigenvectors of an associated Gram matrix. This technique is used to approximate an achievable set for a discrete linear time-invariant dynamical system. This approximation characterizes the part of the state space that may be reached using modest levels of control. If the achievable set can be closely approximated by a proper subspace of the whole space then the system is 'essentially uncontrollable'. The notion finds application in studies of failure-tolerant systems, and in decoupling.
Cliff, E. M.↗
Identification of a linear model of rotor-fuselage dynamics from nonlinear simulation data
Linear regression techniques are used to obtain 9- and 12-degree-of-freedom linear rotorcraft models from the input-output data generated by a nonlinear, blade-element rotorcraft simulation in hover. The resulting models are used to evaluate the coupling of the fuselage modes with the rotor flapping and lead-lag modes at various frequencies. New techniques are proposed and evaluated to improve the identification process, including a method of verifying the assumed model structure by using data sets generated at different input frequencies.
Duval, R. W.↗
Data-Driven Distributed Algorithms for Estimating Eigenvalues and Eigenvectors of Interconnected Dynamical Systems
Here, the paper presents data-driven algorithms to estimate in a distributed manner the eigenvalues, right and left eigenvectors of an unknown linear (or linearized) interconnected dynamic system. In particular, the proposed algorithms do not require the identification of the system model in advance before performing the estimation. As a first step, we consider interconnected dynamical system with distinct eigenvalues. The proposed strategy first estimates the eigenvalues using the well-known Prony method. The right and left eigenvectors are then estimated by solving distributively a set of linear equations. One important feature of the proposed algorithms is that the topology of communication network used to perform the distributed estimation can be chosen arbitrarily, given that it is connected, and is also independent of the structure or sparsity of the system (state) matrix. The proposed distributed algorithms are demonstrated via a numerical example.
97 MATHEMATICS AND COMPUTING↗
Impact of Cross-Axis Structural Dynamics on Validation of Linear Models for Space Launch System
A feasibility study was performed to examine the advisability of incorporating a set of Programmed Test Inputs (PTIs) during the Space Launch System (SLS) vehicle flight. The intent of these inputs is to provide validation to the preflight models for control system stability margins, aerodynamics, and structural dynamics. During October 2009, Ares I-X program was successful in carrying out a series of PTI maneuvers which provided a significant amount of valuable data for post-flight analysis. The resulting data comparisons showed excellent agreement with the preflight linear models across the frequency spectrum of interest. However unlike Ares I-X, the structural dynamics associated with the SLS boost phase configuration are far more complex and highly coupled in all three axes. This presents a challenge when implementing this similar system identification technique to SLS. Preliminary simulation results show noticeable mismatches between PTI validation and analytical linear models in the frequency range of the structural dynamics. An alternate approach was examined which demonstrates the potential for better overall characterization of the system frequency response as well as robustness of the control design.
Pei, Jing↗
Polariton spectra under the collective coupling regime. I. Efficient simulation of linear spectra and quantum dynamics
We outline two general theoretical techniques to simulate polariton quantum dynamics and optical spectra under the collective coupling regimes described by a Holstein–Tavis–Cummings (HTC) model Hamiltonian. The first one takes advantage of sparsity of the HTC Hamiltonian, which allows one to reduce the cost of acting polariton Hamiltonian onto a state vector to the linear order of the number of states, instead of the quadratic order. The second one is applying the well-known Chebyshev series expansion approach for quantum dynamics propagation and to simulate the polariton dynamics in the HTC system; this approach allows us to use a much larger time step for propagation and only requires a few recursive operations of the polariton Hamiltonian acting on state vectors. These two theoretical approaches are general and can be applied to any trajectory-based non-adiabatic quantum dynamics methods. We apply these two techniques with our previously developed Lindblad-partially linearized density matrix approach to simulate the linear absorption spectra of the HTC model system, with both inhomogeneous site energy disorders and dipolar orientational disorders. Our numerical results agree well with the previous analytic and numerical work.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A NASTRAN DMAP alter for linear buckling analysis under dynamic loading
A unique modification to the NASTRAN solution sequence for transient analysis with direct time integration (COSMIC NASTRAN rigid format 9) was developed and incorporated into a DMAP alter. This DMAP alter calculates the buckling stability of a dynamically loaded structure, and is used to predict the onset of structural buckling under stress wave loading conditions. The modified solution sequence incorporates the linear buckling analysis capability (rigid format 5) of NASTRAN into the existing Transient solution rigid format in such a way as to provide a time dependent eigensolution which is used to assess the buckling stability of the structure as it responds to the impulsive load. As a demonstration of the validity of this modified solution procedure, the dynamic buckling of a prismatic bar subjected to an impulsive longitudinal compression is analyzed and compared to the known theoretical solution. In addition, a dynamic buckling analysis is performed for the analytically less tractable problem of the localized dynamic buckling of an initially flawed composite laminate under transverse impact loading. The addition of this DMAP alter to the transient solution sequence in NASTRAN facilitates the prediction of both time and mode of buckling.
Aiello, Robert A.↗
A NASTRAN DMAP alter for linear buckling analysis under dynamic loading
A modification to the NASTRAN solution sequence for transient analysis with direct time integration (COSMIC NASTRAN rigid format 9) was developed and incorporated into a DMAP alter. This DMAP alter calculates the buckling stability of a dynamically loaded structure, and is used to predict the onset of structural buckling under stress-wave loading conditions. The modified solution sequence incorporates the linear buckling analysis capability (rigid format 5) of NASTRAN into the existing Transient solution rigid format in such a way as to provide a time dependent eigensolution which is used to assess the buckling stability of the structure as it responds to the impulsive load. As a demonstration of the validity of this modified solution procedure, the dynamic buckling of a prismatic bar subjected to an impulsive longitudinal compression is analyzed and compared to the known theoretical solution. In addition, a dynamic buckling analysis is performed for the analytically less tractable problem of the localized dynamic buckling of an initially flawed composite laminate under transverse impact loading. The addition of this DMAP alter to the transient solution sequence in NASTRAN facilitates the computational prediction of both the time at which the onset of dynamic buckling occurs in an impulsively loaded structure, and the dynamic buckling mode shapes of that structure.
Aiello, Robert A.↗