SEARCH · Search NASA
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.
System identification techniques are used to represent the dynamic behavior of the Space Shuttle Main Engine. The transfer function matrices of the linearized models of both the closed loop and the open loop system are obtained by using the recursive maximum likelihood method.
System identification techniques are used to represent the dynamic behavior of the Space Shuttle Main Engine. The transfer function matrices of the linearized models of both the closed loop and the open loop system are obtained by using the recursive maximum likelihood method.
A covariance analysis of the performance and sensitivity of the attitude determination Extended Kalman Filter (EKF) used by the On Board Computer (OBC) of the Extreme Ultra Violet Explorer (EUVE) spacecraft is presented. The linearized dynamics and measurement equations of the error states are derived which constitute the truth model describing the real behavior of the systems involved. The design model used by the OBC EKF is then obtained by reducing the order of the truth model. The covariance matrix of the EKF which uses the reduced order model is not the correct covariance of the EKF estimation error. A true covariance analysis has to be carried out in order to evaluate the correct accuracy of the OBC generated estimates. The results of such analysis are presented which indicate both the performance and the sensitivity of the OBC EKF.
This paper presents a covariance analysis of the performance and sensitivity of the attitude determination Extended Kalman Filter (EKF) used by the On Board Computer (OBC) of the Extreme Ultra Violet Explorer (EUVE) spacecraft. The linearized dynamics and measurement equations of the error states are used in formulating the 'truth model' describing the real behavior of the systems involved. The 'design model' used by the OBC EKF is then obtained by reducing the order of the truth model. The covariance matrix of the EKF which uses the reduced order model is not the correct covariance of the EKF estimation error. A 'true covariance analysis' has to be carried out in order to evaluate the correct accuracy of the OBC generated estimates. The results of such analysis are presented which indicate both the performance and the sensitivity of the OBC EKF.
This paper presents a covariance analysis of the performance and sensitivity of the attitude determination Extended Kalman Filter (EKF) used by the onboard computer (OBC) of the Extreme Ultra Violet Explorer (EUVE) spacecraft. The linearized dynamics and measurement equations of the error states are used in formulating the 'truth model' describing the order of the systems involved. The 'design model' used by the OBC EKF is then obtained by reducing the order of the truth model. The covariance matrix of the EKF which uses the reduced order model is not the correct covariance of the EKF estimation error. A 'true covariance analysis' has to be carried out in order to evaluate the correct accuracy of the OBC generated estimates. The results of such analysis are presented which indicate both the performance and the sensitivity of the OBC EKF.
During this period of performance (1 Jan. 1992 - 31 Dec. 1992), the investigation of the following tasks was accomplished: (1) application of the Non-linear Force-free Field (NLFF) Model to study the active region evolution; (2) coronal dynamical responses due to emerging flux including the transition region; (3) loss of MHD equilibrium due to footpoint motions; and (4) two-dimensional MHD global coronal model: steady-state streamers. The detailed description of these studies are included.
A linear, time-dependent, primitive-equation, 3D numerical model is used to investigate the dynamical coupling between the troposphere and stratosphere through the vertical propagation of planetry waves. Attention is given to the roles of the tropopause in affecting the propagation of planetary waves and to the extent to which the stratosphere can determine its own wave driving, i.e., the sensitivity of the generation of the wave activity in the troposphere to the change of zonal wind in the stratosphere. Two maxima in Eliassen-Palm (EP) flux convergence are found for wave 1 in high latitudes, one in the upper stratosphere and the other in the upper troposphere. For wave 2 there is a third maximum north of the zero wind line in the tropical stratosphere. The formation of the maximum in EP-flux convergence in the upper troposphere is related to the dissipation of non-WKBJ waves near a local maximum of the refractive index. The existence of the local maximum in refractive index at the tropopause is a result of the great gradient of the buoyancy frequency as well as the vertical shear of the zonal flow there.
A damped linear oscillator (LRC circuit) predicts a significant amount of the AL index variations when it is driven by the IMF's B(South) component. The model parameters are determined from fitting the AL time series and show a small, apparently nonsystematic variation with the activity level. Following an isolated external disturbance the LRC model's response peaks at 15 min - 0.5 h and decays after 1-2 more hours. The time scales and the correlation between observed AL and the circuit output (60-80 percent) are in good agreement with earlier linear prediction filter results. The results support the modeling of the global magnetospheric behavior with dynamical systems of a few degrees of freedom.
The objective of system identification is to determine reliable dynamical models of a structure by systematically using its measured excitation and response. It brings together in an integrated fashion, experimental, analytical, and computational techniques in structural dynamics. Areas of application for system identification include the following: (1) Model Evaluation--assessing assumptions (linearity and equivalent viscous damping) and techniques (finite-element modeling) used to construct theoretical models of a structure; (2) Model Improvement--updating of a theoretical model to enable more accurate response predictions for possible future loads on the structure, or for control of the structure; (3) Empirical Modelling--developing empirical relationships (nonlinear models) or empirical parameter values (modal damping) because the present state of the art does not provide theoretical results; and (4) Damage Detection and Assessment--continual or episodic updating of a structural model through vibration monitoring to detect and locate any structural damage. It can be argued that since the construction or modification of models using test data is subject to inherent uncertainties, the above problems should be properly treated within a Bayesian probabilistic framework. Such a methodology is presented which allows the precision of the estimates of the model parameters to be computed. It also leads to a guiding principle in applications. Namely, when selecting a single model from a given class of models, one should take the most probable model in the class based on the experimental data. Practical applications of this principle are given which are based on the utilization of measured seismic motions in large civil structures. Examples include the application of a computer program MODE-ID to identify modal properties directly from seismic excitation and response time histories from a nine-story steel-frame building at JPL and from a freeway overpass bridge.
The Annular Suspension and Pointing System (ASPS) developed by the Flight System division of Sperry Corporation is a six-degree of freedom payload pointing system designed for use with the space shuttle. This magnetic suspension and pointing system provides precise controlled pointing in six-degrees of freedom, isolation of payload-carrier disturbances, and end mount controlled pointing. Those are great advantages over the traditional mechanical joints for space applications. In this design, we first analyzed the assumed model of the single degree ASPS bearing actuator and obtained the plant dynamics equations. By linearizing the plant dynamics equations, we designed the cascade and feedback compensators such that a stable and satisfied result was obtained. The specified feedback compensator was computer simulated with the nonlinearized plant dynamics equations. The results indicated that an unstable output occurred. In other words, the designed feedback compensator failed. The failure of the design is due to the Taylor's series expansion not converging.
The dynamics of baroclinic wave growth in a saturated environment is examined using linear and nonlinear models employing a parameterization of latent heat release that assumes all rising air is saturated, and saturation equivalent potential temperature is conserved on ascent. Piecewise potential vorticity (PV) diagnostics are used to interpret the results. When the stability to vertical displacements in saturated air is allowed to increase with height, as it must in an atmosphere with a constant, positive lapse rate of potential temperature, the growth rates of the most unstable modes of the Eady problem grow only marginally faster than the modes of the dry problem. The vertical variation of moist static stability produces a gradient of moist potential vorticity in the rising air, eliminating the short-wave cutoff present in the dry Eady problem. The destabilization of the short waves is shown to be associated with the interaction between surface potential temperature anomalies and diabatically generated lower-tropospheric potential vorticity anomalies. Nonlinear primitive equation simulations, starting from normal-mode initial conditions, show that while the dry wave grows at nearly the linear growth rate until maximum amplitude is reached, the moist wave grows significantly faster than the linear growth rate at finite amplitude. This enhanced growth is associated with the rapid amplification of a mesoscale PV anomaly generated by latent heat release at the warm front. The rapid amplification of the surface cyclone results from the superposition of the circulation associated with this misoscale PV anomaly upon the circulation associated with the surface and upper boundary potential temperature anomalies. Additional integrations with finite-amplitude initial conditions more typical of atmospheric conditions exhibit similar behavior. It is suggested that many of the rapid cyclogenesis events that occur as upper-troposheric PV anomalies cross the east coasts of continents may arise from the rapid generation of PV anomalies by condensational heating in the moist maritime lower troposphere.
This study examines the uncertainty in forecasts of the January-February-March (JFM) mean extratropical circulation, and how that uncertainty is modulated by the El Nino/Southern Oscillation (ENSO). The analysis is based on ensembles of hindcasts made with an Atmospheric General Circulation Model (AGCM) forced with sea surface temperatures observed during; the 1983 El Nino and 1989 La Nina events. The AGCM produces pronounced interannual differences in the magnitude of the extratropical seasonal mean noise (intra-ensemble variability). The North Pacific, in particular, shows extensive regions where the 1989 seasonal mean noise kinetic energy (SKE), which is dominated by a "PNA-like" spatial structure, is more than twice that of the 1983 forecasts. The larger SKE in 1989 is associated with a larger than normal barotropic conversion of kinetic energy from the mean Pacific jet to the seasonal mean noise. The generation of SKE due to sub-monthly transients also shows substantial interannual differences, though these are much smaller than the differences in the mean flow conversions. An analysis of the Generation of monthly mean noise kinetic energy (NIKE) and its variability suggests that the seasonal mean noise is predominantly a statistical residue of variability resulting from dynamical processes operating on monthly and shorter times scales. A stochastically-forced barotropic model (linearized about the AGCM's 1983 and 1989 base states) is used to further assess the role of the basic state, submonthly transients, and tropical forcing, in modulating the uncertainties in the seasonal AGCM forecasts. When forced globally with spatially-white noise, the linear model generates much larger variance for the 1989 base state, consistent with the AGCM results. The extratropical variability for the 1989 base state is dominanted by a single eigenmode, and is strongly coupled with forcing over tropical western Pacific and the Indian Ocean, again consistent with the AGCM results. Linear calculations that include forcing from the AGCM variance of the tropical forcing and submonthly transients show a small impact on the variability over the Pacific/North American region compared with that of the base state differences.
A series of Reynolds-averaged Navier-Stokes calculations were employed to study the performance of rocket-based combined-cycle systems operating in an all-rocket mode. This parametric series of calculations were executed within a statistical framework, commonly known as design of experiments. The parametric design space included four geometric and two flowfield variables set at three levels each, for a total of 729 possible combinations. A D-optimal design strategy was selected. It required that only 36 separate computational fluid dynamics (CFD) solutions be performed to develop a full response surface model, which quantified the linear, bilinear, and curvilinear effects of the six experimental variables. The axisymmetric, Reynolds-averaged Navier-Stokes simulations were executed with the NPARC v3.0 code. The response used in the statistical analysis was created from Isp efficiency data integrated from the 36 CFD simulations. The influence of turbulence modeling was analyzed by using both one- and two-equation models. Careful attention was also given to quantify the influence of mesh dependence, iterative convergence, and artificial viscosity upon the resulting statistical model. Thirteen statistically significant effects were observed to have an influence on rocket-based combined-cycle nozzle performance. It was apparent that the free-expansion process, directly downstream of the rocket nozzle, can influence the Isp efficiency. Numerical schlieren images and particle traces have been used to further understand the physical phenomena behind several of the statistically significant results.
A comparison was made between two approaches to predict acoustic emission waveforms in thin plates. A normal mode solution method for Mindlin plate theory was used to predict the response of the flexural plate mode to a point source, step-function load, applied on the plate surface. The second approach used a dynamic finite element method to model the problem using equations of motion based on exact linear elasticity. Calculations were made using properties for both isotropic (aluminum) and anisotropic (unidirectional graphite/epoxy composite) materials. For simulations of anisotropic plates, propagation along multiple directions was evaluated. In general, agreement between the two theoretical approaches was good. Discrepancies in the waveforms at longer times were caused by differences in reflections from the lateral plate boundaries. These differences resulted from the fact that the two methods used different boundary conditions. At shorter times in the signals, before reflections, the slight discrepancies in the waveforms were attributed to limitations of Mindlin plate theory, which is an approximate plate theory. The advantages of the finite element method are that it used the exact linear elasticity solutions, and that it can be used to model real source conditions and complicated, finite specimen geometries as well as thick plates. These advantages come at a cost of increased computational difficulty, requiring lengthy calculations on workstations or supercomputers. The Mindlin plate theory solutions, meanwhile, can be quickly generated on personal computers. Specimens with finite geometry can also be modeled. However, only limited simple geometries such as circular or rectangular plates can easily be accommodated with the normal mode solution technique. Likewise, very limited source configurations can be modeled and plate theory is applicable only to thin plates.
The Prelaunch loads analysis traditionally captures the vehicle and Mobile Launcher (ML) loads associated with stacking in the Vehicle Assembly Building (VAB), rollout to the pad, and high-speed winds while out on the pad prior to liftoff. The following paper details the methodology used in this analysis for loads for the Artemis Flight Readiness assessment. This paper focuses on the stacking and assembly loads, the loads induced by fueling, and the load seen during rollout to the pad. Methods for integrating multiple structural models in dynamic simulations using interface forces, constraint modes, and residual vectors is discussed. Both linear and nonlinear responses are addressed. Simulation results are compared to measured data from the Dynamic Rollout Test (DRT) and Wet Dress Rehearsal (WDR).
We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.
A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5\% to 5\%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.
A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5% to 5%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.