Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear dynamical systems”

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.

At least 1,315 records · Page 73

Numerical method for solution of systems of non-stationary spatially one-dimensional nonlinear differential equations

A computational scheme and a standard program is proposed for solving systems of nonstationary spatially one-dimensional nonlinear differential equations using Newton's method. The proposed scheme is universal in its applicability and its reduces to a minimum the work of programming. The program is written in the FORTRAN language and can be used without change on electronic computers of type YeS and BESM-6. The standard program described permits the identification of nonstationary (or stationary) solutions to systems of spatially one-dimensional nonlinear (or linear) partial differential equations. The proposed method may be used to solve a series of geophysical problems which take chemical reactions, diffusion, and heat conductivity into account, to evaluate nonstationary thermal fields in two-dimensional structures when in one of the geometrical directions it can take a small number of discrete levels, and to solve problems in nonstationary gas dynamics.

Morozov, S. K.↗

Chemoviscosity modeling for thermosetting resins, 2

A new analytical model for simulating chemoviscosity of thermosetting resin was formulated. The model is developed by modifying the Williams-Landel-Ferry (WLF) theory in polymer rheology for thermoplastic materials. By assuming a linear relationship between the glass transition temperature and the degree of cure of the resin system under cure, the WLF theory can be modified to account for the factor of reaction time. Temperature dependent functions of the modified WLF theory constants were determined from the isothermal cure data of Lee, Loos, and Springer for the Hercules 3501-6 resin system. Theoretical predictions of the model for the resin under dynamic heating cure cycles were shown to compare favorably with the experimental data reported by Carpenter. A chemoviscosity model which is capable of not only describing viscosity profiles accurately under various cure cycles, but also correlating viscosity data to the changes of physical properties associated with the structural transformations of the thermosetting resin systems during cure was established.

Hou, T. H.↗

Stability of a rigid rotor supported on oil-film journal bearings under dynamic load

Most published work relating to dynamically loaded journal bearings are directed to determining the minimum film thickness from the predicted journal trajectories. These do not give any information about the subsynchronous whirl stability of journal bearing systems since they do not consider the equations of motion. It is, however, necessary to know whether the bearing system operation is stable or not under such an operating condition. The stability characteristics of the system are analyzed. A linearized perturbation theory about the equilibrium point can predict the threshold of stability; however it does not indicate postwhirl orbit detail. The linearized method may indicate that a bearing is unstable for a given operating condition whereas the nonlinear analysis may indicate that it forms a stable limit cycle. For this reason, a nonlinear transient analysis of a rigid rotor supported on oil journal bearings under: (1) a unidirectional constant load, (2) a unidirectional periodic load, and (3) variable rotating load are performed. The hydrodynamic forces are calculated after solving the time-dependent Reynolds equation by a finite difference method with a successive overrelaxation scheme. Using these forces, equations of motion are solved by the fourth-order Runge-Kutta method to predict the transient behavior of the rotor. With the aid of a high-speed digital computer and graphics, the journal trajectories are obtained for several different operating conditions.

Majumdar, B. C.↗

Optimization of the structural and control system for LSS with reduced-order model

The objective is the simultaneous design of the structural and control system for space structures. The minimum weight of the structure is the objective function, and the constraints are placed on the closed loop distribution of the frequencies and the damping parameters. The controls approach used is linear quadratic regulator with constant feedback. A reduced-order control system is used. The effect of uncontrolled modes is taken into consideration by the model error sensitivity suppression (MESS) technique which modified the weighting parameters for the control forces. For illustration, an ACOSS-FOUR structure is designed for a different number of controlled modes with specified values for the closed loop damping parameters and frequencies. The dynamic response of the optimum designs for an initial disturbance is compared.

Khot, N. S.↗

Numerical modeling of flows in simulated brush seal configurations

The two-dimensional flow field in a linear brush seal configuration is modeled and the flow field determined numerically using an FEM model consisting of a matrix of 6 x 6 uniformly spaced circular bristles with a uniform flow field upstream of the simulated brush element. Results from the numerical model are compared to earlier work for a circular seal of similar configuration. Temporal disturbance of the upstream flow field is studied to elucidate system dynamics; however, only the local velocity field is investigated. The effect of tangential flow at the upstream boundary is also studied. Tangential flow produces changes in the pressure applied to the bristle but has little effect on the axial flow rates.

Mullen, R. L.↗

Theoretical Study of Decomposition Pathways for HArF and HKrF

To provide theoretical insights into the stability and dynamics of the new rare gas compounds HArF and HKrF, reaction paths for decomposition processes HRgF to Rg + HF and HRgF to H + Rg + F (Rg = Ar, Kr) are calculated using ab initio electronic structure methods. The bending channels, HRgF to Rg + HF, are described by single-configurational MP2 and CCSD(T) electronic structure methods, while the linear decomposition paths, HRgF to H + Rg + F, require the use of multi-configurational wave functions that include dynamic correlation and are size extensive. HArF and HKrF molecules are found to be energetically stable with respect to atomic dissociation products (H + Rg + F) and separated by substantial energy barriers from Rg + HF products, which ensure their kinetic stability. The results are compatible with experimental data on these systems.

Chaban, Galina M.↗

Template Matching Approach to Signal Prediction

A new approach to signal prediction and prognostic assessment of spacecraft health resolves an inherent difficulty in fusing sensor data with simulated data. This technique builds upon previous work that demonstrated the importance of physics-based transient models to accurate prediction of signal dynamics and system performance. While models can greatly improve predictive accuracy, they are difficult to apply in general because of variations in model type, accuracy, or intended purpose. However, virtually any flight project will have at least some modeling capability at its disposal, whether a full-blown simulation, partial physics models, dynamic look-up tables, a brassboard analogue system, or simple hand-driven calculation by a team of experts. Many models can be used to develop a predict, or an estimate of the next day s or next cycle s behavior, which is typically used for planning purposes. The fidelity of a predict varies from one project to another, depending on the complexity of the simulation (i.e. linearized or full differential equations) and the level of detail in anticipated system operation, but typically any predict cannot be adapted to changing conditions or adjusted spacecraft command execution. Applying a predict blindly, without adapting the predict to current conditions, produces mixed results at best, primarily due to mismatches between assumed execution of spacecraft activities and actual times of execution. This results in the predict becoming useless during periods of complicated behavior, exactly when the predict would be most valuable. Each spacecraft operation tends to show up as a transient in the data, and if the transients are misaligned, using the predict can actually harm forecasting performance. To address this problem, the approach here expresses the predict in terms of a baseline function superposed with one or more transient functions. These transients serve as signal templates, which can be relocated in time and space against the signal background. One then has the ability to reconstruct a signal regardless of the precise timing of the transients. During operation, one applies the actual start times of spacecraft activities as they occur, and produces a reconstructed, accurate predict in real-time. This general approach is valid under two important conditions. First, the transients themselves must be time-invariant. Second, the transients must be reasonably consistent with respect to different operating points. Both of these assumptions are generally valid, but in the case of a complicated system with numerous types of overlapping transients, this approach may not be effective. Fortunately, there tend to be few transients in spacecraft telemetry of sensor quantities or low-level health and status information because these signals rarely reflect multiple different types of operation. Furthermore, if the predict is at least reasonably close to actual operation, the shift either in time, or in the operating point when the transient occurs - is likely to be small. The proposed approach considers three ways to recognize a transient. The first and most reliable is to use a different signal that identifies operating mode - often the transient will be correlated to a change in operating mode, and this change can usually be detected positively from discrete signals in spacecraft telemetry. If there is no useful mode signal, the second option is to identify the transient template by hand, and detect the onset of the transient using a curve-fitting approach. Finally, if one elects not to choose by hand, there is an option for automatic selection. This technique has been applied to sensor data from several JPL missions and industrial applications, demonstrating an improvement in accurate prediction of future behavior and early detection of system problems. This technique is applicable to practically any time-varying, quantitative sensor measurement.

Mackey, Ryan↗

Comparison of Aircraft Models and Integration Schemes for Interval Management in the TRACON

Reusable models of common elements for communication, computation, decision and control in air traffic management are necessary in order to enable simulation, analysis and assurance of emergent properties, such as safety and stability, for a given operational concept. Uncertainties due to faults, such as dropped messages, along with non-linearities and sensor noise are an integral part of these models, and impact emergent system behavior. Flight control algorithms designed using a linearized version of the flight mechanics will exhibit error due to model uncertainty, and may not be stable outside a neighborhood of the given point of linearization. Moreover, the communication mechanism by which the sensed state of an aircraft is fed back to a flight control system (such as an ADS-B message) impacts the overall system behavior; both due to sensor noise as well as dropped messages (vacant samples). Additionally simulation of the flight controller system can exhibit further numerical instability, due to selection of the integration scheme and approximations made in the flight dynamics. We examine the theoretical and numerical stability of a speed controller under the Euler and Runge-Kutta schemes of integration, for the Maintain phase for a Mid-Term (2035-2045) Interval Management (IM) Operational Concept for descent and landing operations. We model uncertainties in communication due to missed ADS-B messages by vacant samples in the integration schemes, and compare the emergent behavior of the system, in terms of stability, via the boundedness of the final system state. Any bound on the errors incurred by these uncertainties will play an essential part in a composable assurance argument required for real-time, flight-deck guidance and control systems,. Thus, we believe that the creation of reusable models, which possess property guarantees, such as safety and stability, is an innovative and essential requirement to assessing the emergent properties of novel airspace concepts of operation.

Neogi, Natasha↗

Thermal Fault Tolerance Analysis of Carbon Fiber Rope Barrier Systems for Use in the Reusable Solid Rocket Motor ( RSRM) Nozzle Joints

Carbon Fiber Rope (CFR) thermal barrier systems are being considered for use in several RSRM (Reusable Solid Rocket Motor) nozzle joints as a replacement for the current assembly gap close-out process/design. This study provides for development and test verification of analysis methods used for flow-thermal modeling of a CFR thermal barrier subject to fault conditions such as rope combustion gas blow-by and CFR splice failure. Global model development is based on a 1-D (one dimensional) transient volume filling approach where the flow conditions are calculated as a function of internal 'pipe' and porous media 'Darcy' flow correlations. Combustion gas flow rates are calculated for the CFR on a per-linear inch basis and solved simultaneously with a detailed thermal-gas dynamic model of a local region of gas blow by (or splice fault). Effects of gas compressibility, friction and heat transfer are accounted for the model. Computational Fluid Dynamic (CFD) solutions of the fault regions are used to characterize the local flow field, quantify the amount of free jet spreading and assist in the determination of impingement film coefficients on the nozzle housings. Gas to wall heat transfer is simulated by a large thermal finite element grid of the local structure. The employed numerical technique loosely couples the FE (Finite Element) solution with the gas dynamics solution of the faulted region. All free constants that appear in the governing equations are calibrated by hot fire sub-scale test. The calibrated model is used to make flight predictions using motor aft end environments and timelines. Model results indicate that CFR barrier systems provide a near 'vented joint' style of pressurization. Hypothetical fault conditions considered in this study (blow by, splice defect) are relatively benign in terms of overall heating to nozzle metal housing structures.

Clayton, J. Louie↗

Efficient Computational Model of Hysteresis

A recently developed mathematical model of the output (displacement) versus the input (applied voltage) of a piezoelectric transducer accounts for hysteresis. For the sake of computational speed, the model is kept simple by neglecting the dynamic behavior of the transducer. Hence, the model applies to static and quasistatic displacements only. A piezoelectric transducer of the type to which the model applies is used as an actuator in a computer-based control system to effect fine position adjustments. Because the response time of the rest of such a system is usually much greater than that of a piezoelectric transducer, the model remains an acceptably close approximation for the purpose of control computations, even though the dynamics are neglected. The model (see Figure 1) represents an electrically parallel, mechanically series combination of backlash elements, each having a unique deadband width and output gain. The zeroth element in the parallel combination has zero deadband width and, hence, represents a linear component of the input/output relationship. The other elements, which have nonzero deadband widths, are used to model the nonlinear components of the hysteresis loop. The deadband widths and output gains of the elements are computed from experimental displacement-versus-voltage data. The hysteresis curve calculated by use of this model is piecewise linear beyond deadband limits.

Shields, Joel↗

Systems engineering studies of on-orbit assembly operation

While the practice of construction has a long history, the underlying theory of construction is relatively young. Very little has been documented as to techniques of logistic support, construction planning, construction scheduling, construction testing, and inspection. The lack of 'systems approaches' to construction processes is certainly one of the most serious roadblocks to the construction of space structures. System engineering research efforts at CSC are aimed at developing concepts and tools which contribute to a systems theory of space construction. The research is also aimed at providing means for trade-offs of design parameters for other research areas in CSC. Systems engineering activity at CSC has divided space construction into the areas of orbital assembly, lunar base construction, interplanetary transport vehicle construction, and Mars base construction. A brief summary of recent results is given. Several models for 'launch-on-time' were developed. Launch-on-time is a critical concept to the assembly of such Earth-orbiting structures as the Space Station Freedom, and to planetary orbiters such as the Mars transfer vehicle. CSC has developed a launch vehicle selection model which uses linear programming to find optimal combinations of launch vehicles of various sizes (Atlas, Titan, Shuttles, HLLV's) to support SEI missions. Recently, the Center developed a cost trade-off model for studying on orbit assembly logistics. With this model it was determined that the most effective size of the HLLV would be in the range of 120 to 200 metric tons to LEO, which is consistent with the choices of General Stafford's Synthesis Group Report. A second-generation Dynamic Construction Activities Model ('DYCAM') process model has been under development, based on our past results in interruptability and our initial DYCAM model. This second-generation model is built on the paradigm of knowledge-based expert systems. It is aimed at providing answers to two questions: (1) what are some necessary or sufficient conditions for judging conceptual designs of spacecraft?, and (2) can a methodology be formulated such that these conditions may be used to provide computer-aided tools for evaluating conceptual designs and planning for space assembly sequences? Early simulation results indicate that the DYCAM model has a clear ability to emulate and simulate human orbital construction processes.

Morgenthaler, George W.↗

A simple Fourier analysis technique for measuring the dynamic response of manual control systems.

From a practical standpoint most dynamic response analysis techniques are complicated, expensive, and time consuming to mechanize and operate. Adaptive parameter tracking techniques partially overcome these faults, but they have prescribed forms which is an undesirable restriction for basic research. To overcome these drawbacks, an on-line Fourier analysis technique has been developed which deals only with the input and error signals in a manual feedback control system. Several sine waves are used for the input, and the on-line data measurement includes the error variance and a simple means for obtaining the sine and cosine transforms of the error signal at each of the input frequencies. From this relatively small amount of data, we show simple ways to compute the open- and closed-loop dynamic response at input frequencies and the relative amount of linearly correlated power in the error signal.

Allen, R. W.↗

Shape determination and control for large space structures

An integral operator approach is used to derive solutions to static shape determination and control problems associated with large space structures. Problem assumptions include a linear self-adjoint system model, observations and control forces at discrete points, and performance criteria for the comparison of estimates or control forms. Results are illustrated by simulations in the one dimensional case with a flexible beam model, and in the multidimensional case with a finite model of a large space antenna. Modal expansions for terms in the solution algorithms are presented, using modes from the static or associated dynamic mode. These expansions provide approximated solutions in the event that a used form analytical solution to the system boundary value problem is not available.

Weeks, C. J.↗

Feedback control laws for highly maneuverable aircraft

The results of a study of the application of H infinity and mu synthesis techniques to the design of feedback control laws for the longitudinal dynamics of the High Angle of Attack Research Vehicle (HARV) are presented. The objective of this study is to develop methods for the design of feedback control laws which cause the closed loop longitudinal dynamics of the HARV to meet handling quality specifications over the entire flight envelope. Control law designs are based on models of the HARV linearized at various flight conditions. The control laws are evaluated by both linear and nonlinear simulations of typical maneuvers. The fixed gain control laws resulting from both the H infinity and mu synthesis techniques result in excellent performance even when the aircraft performs maneuvers in which the system states vary significantly from their equilibrium design values. Both the H infinity and mu synthesis control laws result in performance which compares favorably with an existing baseline longitudinal control law.

Garrard, William L.↗

Linear and nonlinear trending and prediction for AVHRR time series data

The variability of AVHRR calibration coefficient in time was analyzed using algorithms of linear and non-linear time series analysis. Specifically we have used the spline trend modeling, autoregressive process analysis, incremental neural network learning algorithm and redundancy functional testing. The analysis performed on available AVHRR data sets revealed that (1) the calibration data have nonlinear dependencies, (2) the calibration data depend strongly on the target temperature, (3) both calibration coefficients and the temperature time series can be modeled, in the first approximation, as autonomous dynamical systems, (4) the high frequency residuals of the analyzed data sets can be best modeled as an autoregressive process of the 10th degree. We have dealt with a nonlinear identification problem and the problem of noise filtering (data smoothing). The system identification and filtering are significant problems for AVHRR data sets. The algorithms outlined in this study can be used for the future EOS missions. Prediction and smoothing algorithms for time series of calibration data provide a functional characterization of the data. Those algorithms can be particularly useful when calibration data are incomplete or sparse.

Smid, J.↗

Newton-Euler Dynamic Equations of Motion for a Multi-body Spacecraft

The Magnetospheric MultiScale (MMS) mission employs a formation of spinning spacecraft with several flexible appendages and thruster-based control. To understand the complex dynamic interaction of thruster actuation, appendage motion, and spin dynamics, each spacecraft is modeled as a tree of rigid bodies connected by spherical or gimballed joints. The method presented facilitates assembling by inspection the exact, nonlinear dynamic equations of motion for a multibody spacecraft suitable for solution by numerical integration. The building block equations are derived by applying Newton's and Euler's equations of motion to an "element" consisting of two bodies and one joint (spherical and gimballed joints are considered separately). Patterns in the "mass" and L'force" matrices guide assembly by inspection of a general N-body tree-topology system. Straightforward linear algebra operations are employed to eliminate extraneous constraint equations, resulting in a minimum-dimension system of equations to solve. This method thus combines a straightforward, easily-extendable, easily-mechanized formulation with an efficient computer implementation.

Stoneking, Eric↗

Simulation of Dynamics of a Flexible Miniature Airplane

A short report discusses selected aspects of the development of the University of Florida micro-aerial vehicle (UFMAV) basically, a miniature airplane that has a flexible wing and is representative of a new class of airplanes that would operate autonomously or under remote control and be used for surveillance and/or scientific observation. The flexibility of the wing is to be optimized such that passive deformation of the wing in the presence of aerodynamic disturbances would reduce the overall response of the airplane to disturbances, thereby rendering the airplane more stable as an observation platform. The aspect of the development emphasized in the report is that of computational simulation of dynamics of the UFMAV in flight, for the purpose of generating mathematical models for use in designing control systems for the airplane. The simulations are performed by use of data from a wind-tunnel test of the airplane in combination with commercial software, in which are codified a standard set of equations of motion of an airplane, and a set of mathematical routines to compute trim conditions and extract linear state space models.

Waszak, Martin R.↗

Signal-Induced Noise Effects in a Photon Counting System for Stratospheric Ozone Measurement

Signal-induced noise (SIN) is a common effect resulting when a photomultiplier tube (PMT) is saturated, for a brief moment, with a high intensity light pulse. After the laser pulse is sent into the atmosphere a very large light return, from either the near-field or a cloud, causes the PMT to momentarily saturate. The PMT is gated off at this time so no signal is seen at the anode. When the PMT gate is turned on, the far-field light return from the atmosphere is observed. This signal is distorted, however because of the addition of SIN to the received light signal causing a slower than expected decay of the atmospheric signal return. We have characterized SIN responses to varying parameters of the incident light on the PMT. These varied parameters included incident wavelength, PMT voltage, incident intensity, and tube type. We found that only the amplitude of the SIN was effected by varying PMT voltages and light intensities. The amplitude increased linearly as input light intensity increased. Different incident wavelengths at the same intensity did not effect the amplitude or the temporal behavior of the SIN response. Finally, different PMT tubes with similar physical structures exhibited similar SIN responses although with different amplitudes. The different amplitudes can be attributed to the different gains and operating voltages of each tube. These results suggest that SIN is caused by photocathode electron dynamics such as charge accumulation on internal PMT surfaces. These surfaces then emit the electrons slowly resulting in a long decay noise signal. With the SIN responses characterized we can now try to develop a method to reduce or eliminate SIN in DIAL systems.

Harper, David B.↗