Parametrically-defined Approximate Control Variates
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Multi Model Monte Carlo with Python (\mxmc {}) is a software package developed as a general capability for computing the statistics of outputs from an expensive, high-fidelity model by leveraging faster, low-fidelity models for speedup. Motivated by uncertainty propagation problems where classical Monte Carlo (MC) simulation is computationally intractable, various multi-model MC approaches have recently emerged that yield unbiased estimators with significantly reduced variance relative to MC for the same cost. These existing methods include multi-level Monte Carlo (MLMC), multi-fidelity Monte Carlo (MFMC), and approximate control variates (ACV). Given a fixed computational budget and a collection of models with varying cost/accuracy, each method seeks a sample allocation strategy across the models that results in an estimator with optimal variance reduction. \mxmc {} is a versatile tool that enables convenient access to many existing multi-model MC approaches within one modular and extensible package. With \mxmc {}, users can easily compare existing methods to determine the best choice for their particular problem, while developers have a basis for implementing and sharing new variance reduction approaches. This report introduces the \mxmc {} software, providing a summary of the problem-solving workflow for users as well as a brief overview of the code layout for developers.
Multi-model Monte Carlo methods have been demonstrated to be an efficient and accurate alternative to standard Monte Carlo (MC) in the model-based propagation of uncertainty in entry, descent, and landing (EDL) applications. These multi-model MC methods fuse predictions from low-fidelity models with the high-fidelity EDL model of interest to produce unbiased statistics with a fraction of the computational cost. The accuracy and efficiency of the multi-model MC methods are dependent upon the magnitude of correlations of the low-fidelity models with the high-fidelity model, but also upon the correlation among the low-fidelity models, and their relative computational cost. Because of this layer of complexity, the question of how to optimally select the set of low-fidelity models has remained open. In this work, methods for optimal model construction and tuning are investigated as a means to increase the speed and precision of trajectory simulation for EDL. Specifically, the focus is on the inclusion of low-fidelity model tuning within the sample allocation optimization that accompanies multi-model MC methods. Preliminary results indicate that low-fidelity model tuning can significantly improve efficiency and precision of trajectory simulations and provide an increased edge to multi-model MC methods when compared to standard MC. The challenges and potential benefits to exploring a fully iterative and comprehensive optimization strategy in future work are highlighted.
Optimal trajectory and control approximate solution by asymptotic expansion, considering Mayer variational problem solution
The results of flight tests to determine flying qualities of a Chance Vought F4U-4 airplane are presented and discussed herein. In addition to comprehensive measurements at low altitude (about 8000 ft), tests of limited scope were made at high altitude (about 25,000 ft). The more important characteristics, based on a comparison of the test results and opinions of the pilots with the Navy requirements, can be summarized as follows: 1. The short-period control-free oscillations of the elevator angle and the normal acceleration were satisfactorily damped. 2. The most rearward center-of-gravity locations for satisfactory static longitudinal stability with power on, as determined by the control-force variations, were approximately 30 and 27 percent M.A.C. with flaps and gear up and down, respectively. 3. In maneuvering flight the conditions for which control-force gradients of satisfactory magnitude were obtained were seriously limited by sizable changes in the gradient with center-of-gravity location, airspeed, altitude, acceleration factor, and direction of turn. 4. The elevator and rudder controls were satisfactory for landings and take-offs. 5. The trim tabs were sufficiently effective for all controls. 6. The directional and lateral dynamic stability was positive, but the rudder oscillation did not damp within one cycle. The airplane oscillation damped sufficiently at low altitude but not at high altitude. 7. Both rudder-fixed and rudder-free static directional stability were positive over a sideslip range of +/-15 deg. However, the rudder force tended to reverse at high angles of right sideslip with flaps and gear up, power on, at low speeds. 8. The stick-fixed static lateral stability (dihedral effect) was positive in all conditions, but the stick-free dihedral effect was neutral at low speeds with flap and gear down, power on. 9. The yaw due to abrupt full aileron deflection at low speed was mot excessive, and the rudder control was adequate to hold trim sideslip. 10. In abrupt rudder-fixed aileron rolls in the clean configuration the maximum pb/2V for full aileron deflection at low and normal speeds was only 0.064. 11. The stalling characteristics were considered unsatisfactory in all configurations in both straight and turning flight due to inadequate stall warning. The motions in the stalls were not unduly severe, and recovery could be effected promptly by normal use of the controls.
The effects of control mode and rate on the fracture toughness of ceramics were measured by using chevron-notched flexure specimen in accordance with ASTM C1421. The use of stroke control gave consistent results with about 2% variation in measured fracture toughness for a very wide range of rates (0.005 to 0.5 mm/min). Use of strain or CMOD control gave an approximately 5% variation over a very wide range of rates, with the measurements being a function of rate. However, the effect was eliminated by use of dry nitrogen, implying a stress corrosion effect rather than a stability effect. With the use of nitrogen for strain control, fracture toughness values were within about 1% over a wide range of rates (1 to 80 micro-epsilon/sec.). CMOD or strain control did allow stable crack extension well past maximum load, and thus is preferred for energy calculations. The effort is being used to confirm recommendations in ASTM International standard C1421 on fracture toughness measurement.
Optimum design of dominant type adaptive control systems with large parameter variations, using fourth order approximation
This paper presents a neural network (NN) approximator-based online reinforcement learning (ORL) controller design for Mars Ascent Vehicle (MAV) under parametric variation and significant external disturbances. The ORL controller, which does not require any offline training, involves two NNs where an action NN produces optimal short-term control performance while a critic NN evaluates the performance of the action NN using an approximated cost function. The simulation example with comparisons against baseline Proportional-Integral-Derivative (PID) and gain scheduled pole-placement PID (GS-PP-PID) controllers show the proposed controller’s effectiveness and robustness under parametric variation and high external disturbances.
The paper describes the development of a robust control methodology that can tolerate large uncertainties in plant dynamics in terms of proportional gain variations as well as frequencies and damping. For this to be accomplished, the methodology starts by designing a feedback controller based on a loop shaping technique. The design is then improved by employing a root locus control design approach. Finally, the response of the feedback system is further improved by augmenting the control design using model reference control. The design approach also takes into consideration the speed of the process actuator for a fast response and for high disturbance rejection. Results from two control design examples covered in the paper show that this control approach can tolerate plant proportional gain uncertainty of approximately ± two orders of magnitude with relatively large variations in plant frequencies and damping.
Candidate autopilot control laws that control the engine failure transient sink rates by demonstrating the engineering application of modern state variable control theory were defined. The results of approximate modal analysis were compared to those derived from full state analyses provided from computer design solutions. The aircraft was described, and a state variable model of its longitudinal dynamic motion due to engine and control variations was defined. The classical fast and slow modes were assumed to be sufficiently different to define reduced order approximations of the aircraft motion amendable to hand analysis control definition methods. The original state equations of motion were also applied to a large scale state variable control design program, in particular OPTSYS. The resulting control laws were compared with respect to their relative responses, ease of application, and meeting the desired performance objectives.
An organic Rankine cycle engine is used as part of a solar thermal power conversion assembly (PCA). The PCA, including a direct-heated cavity receiver and a shaft-mounted alternator, is mounted at the focal point of a parabolic dish concentrator. The engine controls are required to maintain approximately constant values of turbine inlet temperature and shaft speed, despite variation in the concentrated solar power input to the receiver. The controls design approach, system models, and initial stability and performance analysis results are presented herein.
A multi-objective optimal control technique is modified to accommodate changing cost function weights and is used to control a flexible wing aircraft model. Variation of the weights is used to adjust the relative importance of each objective according to either a prescribed function of time or of the state. Several techniques for obtaining a practical approximation to the optimal control solution are presented, and stability of a specific weight structure with the optimal controller is demonstrated. Functionality of the multi-objective control design with weight variation is demonstrated in simulation of a flexible wing transport aircraft and is shown to improve performance over the fixed weight version both at a constant flight condition and across changing flight conditions.
Lower stratospheric air trajectories entering the region over Alaska at the approximately 125 mb level during late May, 1975 indicate a substantial shift in the geographical source regions for the air masses present during that time. This shift coincides with an approximately 25% decrease in the observed halocarbon mixing ratios at the 125 mb level as determined from a daily sequence of halocarbon profiles. Since the halocarbon species measured are essentially chemically inactive at this level, the observed variation is linked to the changing meteorological pattern.
Decoupling theory and robust servomechanism theory are applied to the design of linear multivariable systems with large parameter variations. In addition to being approximately decoupled in the transient period, the over-all system achieves tracking and disturbance rejection robustly in the steady state. An example in flight control system is given.
Regionally confined interannual perturbations in the tropical radiation budget associated with east-west circulations are examined along with the forcing mechanisms of these perturbations in relation to surface and cloudiness variability. A zonal harmonic analysis of emitted longwave radiation emphasizes that these variations are largely controlled at the planetary wave scale. The overall effect leads to an approximately 50 deg/yr propagation phase speed that is considerably slower than the oceanic Kelvin wave capable of driving east-west longwave (LW) anomalies through SST feedback. Examination of the time-dependent radiative energetics over the tropics reveals that the aforementioned anomaly LW propagation is mainly due to cloud forcing associated with east-west circulation changes. Diabatic heating associated with coupled ocean-atmosphere feedback appears to be larger responsible for the LW anomaly propagation.
An electron emission control circuit of the switching regulator type operating at 100 kHz has been developed which maintains a constant emission current within 0.1% for a cathode power demand variation of approximately 100%. The power output stage has an efficiency of 67%, and the overall efficiency is 45% when driving a thoria-coated iridium cathode having a nominal resistance at operating temperature of 2.5 ohms. Under optimum conditions, the bus power demand is 1.75 W. The circuit is useful in controlling the electron emission current of ion sources in applications which involve a substantial variation of the cathode work function, such as oxygen partial pressure measurements over a large dynamic range.
Chemical data for the parental magmas of the nine known VLT Array I pyroclastic glasses show statistically significant trends due to olivine (approximately Fo83) control during their formation by partial melting. The compositional scatter is largely due to compositional variations in the source regions on the + or - 1 percent level. This compositional scatter is small when one considers that the scale of the source region is up to 1000 km, but is sufficient to make positive identification of the residual phase(s) in the source regions difficult. Nevertheless, when the effects of the scatter are properly modeled, it is relatively clear that olivine is the residual phase in the source region. Hence these data and additional constraints indicate that the source regions are at shallow depths in the moon.