Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “response surface”

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 55 records · Page 3

Evaluation of Response Surface Experiment Designs for Distributed Propulsion Aircraft Aero-Propulsive Modeling

Modern distributed hybrid and electric propulsion aircraft, including vertical, short, and conventional takeoff and landing configurations, exhibit significant aero-propulsive complexity and a large number of interacting test factors. This paper presents the development and evaluation of experiment designs for aero-propulsive characterization of distributed propulsion aircraft. Five different foundational response surface designs are evaluated to inform the development of two sequential design approaches tailored to complex aircraft aerodynamic characterization experiments. The first approach, which builds on sequential face-centered central composite designs, has been used previously to develop aero-propulsive models for complex aircraft using wind tunnel testing. The second approach is a new design strategy leveraging a regular I-optimal and nested I-optimal design that was developed for this study. The two sequential design strategies are compared for experiments with a large number of test factors using pre-experiment design evaluation metrics, as well as modeling results obtained from simulated wind tunnel data for the NASA LA-8 aircraft. The design evaluation metrics show that the sequential I-optimal base design has higher statistical power, lower correlation among candidate regressors, lower prediction variance, and more precise parameter estimates. The simulated wind tunnel experiments conducted using each design reveal that the sequential I-optimal base design has better predictive capability with fewer test points. The experiment design and evaluation procedures are described in detail to inform future aerodynamic characterization experiments for complex aircraft.

design of experiments↗

Dynamic Probabilistic Risk Assessment Based Response Surface Approach for FLEX and Accident Tolerant Fuels for Medium Break LOCA Spectrum

After the Fukushima Daiichi Accident, the safety features such as accident tolerant fuel (ATF) and diverse and flexible coping strategies (FLEX) for existing nuclear fleets are being investigated by the US Department of Energy under the Light Water Reactor Sustainability Program. This research is being conducted to quantify the risk-benefit of these safety features. Dynamic probabilistic risk assessment (DPRA)-based response-surface approach has been presented to quantify the FLEX and ATF benefits by estimating the risk associated with each option. ATFs with multilayered silicon carbide (SiC), iron-chromium-aluminum, and chromium-coated zirconium cladding were considered in this study. While these ATF candidates perform better than the current zirconium cladding (Zr), they may introduce additional failure modes in some operating conditions. The fuel failure analysis modules (FAMs) were developed to investigate ATF performance. The dynamic risk assessments were performed using RAVEN, a DPRA tool, coupled with RELAP5 and FAMs. A cumulative distribution function-based index provided a mean of comparing the benefits of safety enhancements. For medium break loss of coolant accidents, FLEX operational timing window for each fuel type was estimated. Among these ATF candidates, SiC-type ATF was the most beneficial candidate for an increased safety margin than Zr-based fuel and was found to complement FLEX strategies in terms of risk and coping time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Numbers of center points appropriate to blocked response surface experiments

Tables are given for the numbers of center points to be used with blocked sequential designs of composite response surface experiments as used in empirical optimum seeking. The star point radii for exact orthogonal blocking is presented. The center point options varied from a lower limit of one to an upper limit equal to the numbers proposed by Box and Hunter for approximate rotatability and uniform variance, and exact orthogonal blocking. Some operating characteristics of the proposed options are described.

Holms, A. G.↗

Numbers of center points appropriate to blocked response surface experiments

Tables are given for the numbers of center points to be used with blocked sequential designs of composite response surface experiments as used in empirical optimum seeking. The tables also give the star point radii for exact orthogonal blocking. The center point options vary from a lower limit of one to an upper limit equal to the numbers proposed by Box and Hunter for approximate rotatability and uniform variance, and exact orthogonal blocking. Some operating characteristics of the proposed options are described.

Holms, A. G.↗

Reliability Sensitivity Analysis and Design Optimization of Composite Structures Based on Response Surface Methodology

This report discusses the development and application of two alternative strategies in the form of global and sequential local response surface (RS) techniques for the solution of reliability-based optimization (RBO) problems. The problem of a thin-walled composite circular cylinder under axial buckling instability is used as a demonstrative example. In this case, the global technique uses a single second-order RS model to estimate the axial buckling load over the entire feasible design space (FDS) whereas the local technique uses multiple first-order RS models with each applied to a small subregion of FDS. Alternative methods for the calculation of unknown coefficients in each RS model are explored prior to the solution of the optimization problem. The example RBO problem is formulated as a function of 23 uncorrelated random variables that include material properties, thickness and orientation angle of each ply, cylinder diameter and length, as well as the applied load. The mean values of the 8 ply thicknesses are treated as independent design variables. While the coefficients of variation of all random variables are held fixed, the standard deviations of ply thicknesses can vary during the optimization process as a result of changes in the design variables. The structural reliability analysis is based on the first-order reliability method with reliability index treated as the design constraint. In addition to the probabilistic sensitivity analysis of reliability index, the results of the RBO problem are presented for different combinations of cylinder length and diameter and laminate ply patterns. The two strategies are found to produce similar results in terms of accuracy with the sequential local RS technique having a considerably better computational efficiency.

Rais-Rohani, Masoud↗

Response Surface Modeling Tolerance and Inference Error Risk Specifications: Proposed Industry Standards

This paper reviews the derivation of an equation for scaling response surface modeling experiments. The equation represents the smallest number of data points required to fit a linear regression polynomial so as to achieve certain specified model adequacy criteria. Specific criteria are proposed which simplify an otherwise rather complex equation, generating a practical rule of thumb for the minimum volume of data required to adequately fit a polynomial with a specified number of terms in the model. This equation and the simplified rule of thumb it produces can be applied to minimize the cost of wind tunnel testing.

DeLoach, Richard↗

Computer program for post-flight evaluation of the control surface response for an attitude controlled missile

A FORTRAN IV coded computer program is presented for post-flight analysis of a missile's control surface response. It includes preprocessing of digitized telemetry data for time lags, biases, non-linear calibration changes and filtering. Measurements include autopilot attitude rate and displacement gyro output and four control surface deflections. Simple first order lags are assumed for the pitch, yaw and roll axes of control. Each actuator is also assumed to be represented by a first order lag. Mixing of pitch, yaw and roll commands to four control surfaces is assumed. A pseudo-inverse technique is used to obtain the pitch, yaw and roll components from the four measured deflections. This program has been used for over 10 years on the NASA/SCOUT launch vehicle for post-flight analysis and was helpful in detecting incipient actuator stall due to excessive hinge moments. The program is currently set up for a CDC CYBER 175 computer system. It requires 34K words of memory and contains 675 cards. A sample problem presented herein including the optional plotting requires eleven (11) seconds of central processor time.

Knauber, R. N.↗

Classifying Multi-Model Wheat Yield Impact Response Surfaces Showing Sensitivity to Temperature and Precipitation Change

Crop growth simulation models can differ greatly in their treatment of key processes and hence in their response to environmental conditions. Here, we used an ensemble of 26 process-based wheat models applied at sites across a European transect to compare their sensitivity to changes in temperature (minus 2 to plus 9 degrees Centigrade) and precipitation (minus 50 to plus 50 percent). Model results were analysed by plotting them as impact response surfaces (IRSs), classifying the IRS patterns of individual model simulations, describing these classes and analysing factors that may explain the major differences in model responses. The model ensemble was used to simulate yields of winter and spring wheat at four sites in Finland, Germany and Spain. Results were plotted as IRSs that show changes in yields relative to the baseline with respect to temperature and precipitation. IRSs of 30-year means and selected extreme years were classified using two approaches describing their pattern. The expert diagnostic approach (EDA) combines two aspects of IRS patterns: location of the maximum yield (nine classes) and strength of the yield response with respect to climate (four classes), resulting in a total of 36 combined classes defined using criteria pre-specified by experts. The statistical diagnostic approach (SDA) groups IRSs by comparing their pattern and magnitude, without attempting to interpret these features. It applies a hierarchical clustering method, grouping response patterns using a distance metric that combines the spatial correlation and Euclidian distance between IRS pairs. The two approaches were used to investigate whether different patterns of yield response could be related to different properties of the crop models, specifically their genealogy, calibration and process description. Although no single model property across a large model ensemble was found to explain the integrated yield response to temperature and precipitation perturbations, the application of the EDA and SDA approaches revealed their capability to distinguish: (i) stronger yield responses to precipitation for winter wheat than spring wheat; (ii) differing strengths of response to climate changes for years with anomalous weather conditions compared to period-average conditions; (iii) the influence of site conditions on yield patterns; (iv) similarities in IRS patterns among models with related genealogy; (v) similarities in IRS patterns for models with simpler process descriptions of root growth and water uptake compared to those with more complex descriptions; and (vi) a closer correspondence of IRS patterns in models using partitioning schemes to represent yield formation than in those using a harvest index. Such results can inform future crop modelling studies that seek to exploit the diversity of multi-model ensembles, by distinguishing ensemble members that span a wide range of responses as well as those that display implausible behaviour or strong mutual similarities.

Classification; Climate change; Crop model; Ensemb↗

Development of a Response Surface Thermal Model for Orion Mated to the International Space Station

A study was performed to determine if a Design of Experiments (DOE)/Response Surface Methodology could be applied to on-orbit thermal analysis and produce a set of Response Surface Equations (RSE) that accurately predict vehicle temperatures. The study used an integrated thermal model of the International Space Station and the Orion Outer mold line model. Five separate factors were identified for study: yaw, pitch, roll, beta angle, and the environmental parameters. Twenty external Orion temperatures were selected as the responses. A DOE case matrix of 110 runs was developed. The data from these cases were analyzed to produce an RSE for each of the temperature responses. The initial agreement between the engineering data and the RSE predictions was encouraging, although many RSEs had large uncertainties on their predictions. Fourteen verification cases were developed to test the predictive powers of the RSEs. The verification showed mixed results with some RSE predicting temperatures matching the engineering data within the uncertainty bands, while others had very large errors. While this study to not irrefutably prove that the DOE/RSM approach can be applied to on-orbit thermal analysis, it does demonstrate that technique has the potential to predict temperatures. Additional work is needed to better identify the cases needed to produce the RSEs

Miller, Stephen W.↗

Improvements to a Response Surface Thermal Model for Orion

A study was performed to determine if a Design of Experiments (DOE)/Response Surface Methodology could be applied to on-orbit thermal analysis and produce a set of Response Surface Equations (RSE) that predict Orion vehicle temperatures within 10 F. The study used the Orion Outer Mold Line model. Five separate factors were identified for study: yaw, pitch, roll, beta angle, and the environmental parameters. Twenty-three external Orion components were selected and their minimum and maximum temperatures captured over a period of two orbits. Thus, there are 46 responses. A DOE case matrix of 145 runs was developed. The data from these cases were analyzed to produce a fifth order RSE for each of the temperature responses. For the 145 cases in the DOE matrix, the agreement between the engineering data and the RSE predictions was encouraging with 40 of the 46 RSEs predicting temperatures within the goal band. However, the verification cases showed most responses did not meet the 10 F goal. After reframing the focus of the study to better align the RSE development with the purposes of the model, a set of RSEs for both the minimum and maximum radiator temperatures was produced which predicted the engineering model output within +/-4 F. Therefore, with the correct application of the DOE/RSE methodology, RSEs can be developed that provide analysts a fast and easy way to screen large numbers of environments and assess proposed changes to the RSE factors.

Miller, Stephen W.↗

Improvements to a Response Surface Thermal Model for Orion Mated to the International Space Station

This study is an extension of previous work to evaluate the applicability of Design of Experiments (DOE)/Response Surface Methodology to on-orbit thermal analysis. The goal was to determine if the methodology could produce a Response Surface Equation (RSE) that predicted the thermal model temperature results within +/-10 F. An RSE is a polynomial expression that can then be used to predict temperatures for a defined range of factor combinations. Based on suggestions received from the previous work, this study used a model with simpler geometry, considered polynomials up to fifth order, and evaluated orbital temperature variations to establish a minimum and maximum temperature for each component. A simplified Outer Mold Line (OML) thermal model of the Orion spacecraft was used in this study. The factors chosen were the vehicle's Yaw, Pitch, and Roll (defining the on-orbit attitude), the Beta angle (restricted to positive beta angles from 0 to 75), and the environmental constants (varying from cold to hot). All factors were normalized from their native ranges to a non-dimensional range from -1.0 to 1.0. Twenty-three components from the OML were chosen and the minimum and maximum orbital temperatures were calculated for each to produce forty-six responses for the DOE model. A customized DOE case matrix of 145 analysis cases was developed which used analysis points at the factor corners, mid-points, and center. From this data set, RSE s were developed which consisted of cubic, quartic, and fifth order polynomials. The results presented are for the fifth order RSE. The RSE results were then evaluated for agreement with the analytical model predictions to produce a +/-3(sigma) error band. Forty of the 46 responses had a +/-3(sigma) value of 10 F or less. Encouraged by this initial success, two additional sets of verification cases were selected. One contained 20 cases, the other 50 cases. These cases were evaluated both with the fifth order RSE and with the analytical model. For the maximum temperature predictions, 12 of the 23 components had all predictions within +/-10 F and 17 were within +/-20 F. For the minimum temperature predictions, only 4 of the 23 components (the four radiator temperatures), were within the 10 F goal. The maximum temperature RSEs were then run through 59,049 screening cases. The RSE predictions were then filtered to find 55 cases that produced the hottest temperatures. These 55 cases were then analyzed using the thermal model and the results compared against the RSE predictions. As noted earlier, 12 of the 23 responses were within +/-10 F at 17 within +/-20 F. These results demonstrate that if properly formulated, an RSE can provide a reliable, fast temperature prediction. Despite this progress, additional work is needed to determine why the minimum temperatures responses and 6 of the hot temperature responses did not produce reliable RSEs. Recommend focus areas are the model itself (arithmetic vs. diffusion nodes) and seeking consultations with statistical application experts.

Miller, StephenW.↗

Assessment of Response Surface Models using Independent Confirmation Point Analysis

This paper highlights various advantages that confirmation-point residuals have over conventional model design-point residuals in assessing the adequacy of a response surface model fitted by regression techniques to a sample of experimental data. Particular advantages are highlighted for the case of design matrices that may be ill-conditioned for a given sample of data. The impact of both aleatory and epistemological uncertainty in response model adequacy assessments is considered.

DeLoach, Richard↗

Application of Neural Networks to Wind tunnel Data Response Surface Methods

The integration of nonlinear neural network methods with conventional linear regression techniques is demonstrated for representative wind tunnel force balance data modeling. This work was motivated by a desire to formulate precision intervals for response surfaces produced by neural networks. Applications are demonstrated for representative wind tunnel data acquired at NASA Langley Research Center and the Arnold Engineering Development Center in Tullahoma, TN.

Lo, Ching F.↗

Shape Optimization of Supersonic Turbines Using Response Surface and Neural Network Methods

Turbine performance directly affects engine specific impulse, thrust-to-weight ratio, and cost in a rocket propulsion system. A global optimization framework combining the radial basis neural network (RBNN) and the polynomial-based response surface method (RSM) is constructed for shape optimization of a supersonic turbine. Based on the optimized preliminary design, shape optimization is performed for the first vane and blade of a 2-stage supersonic turbine, involving O(10) design variables. The design of experiment approach is adopted to reduce the data size needed by the optimization task. It is demonstrated that a major merit of the global optimization approach is that it enables one to adaptively revise the design space to perform multiple optimization cycles. This benefit is realized when an optimal design approaches the boundary of a pre-defined design space. Furthermore, by inspecting the influence of each design variable, one can also gain insight into the existence of multiple design choices and select the optimum design based on other factors such as stress and materials considerations.

Papila, Nilay↗

Distinct Surface Response to Black Carbon Aerosols

For the radiative impact of individual climate forcings, most previous studies focused on the global mean values at the top of the atmosphere (TOA), and less attention has been paid to surface processes, especially for black carbon (BC) aerosols. In this study, the surface radiative responses to five different forcing agents were analyzed by using idealized model simulations. Our analyses reveal that for greenhouse gases, solar irradiance, and scattering aerosols, the surface temperature changes are mainly dictated by the changes of surface radiative heating, but for BC, surface energy redistribution between different components plays a more crucial role. Globally, when a unit BC forcing is imposed at TOA, the net shortwave radiation at the surface decreases by −5.87±0.67 W m −2 (W m −2 ) −1 (averaged over global land without Antarctica), which is partially offset by increased downward longwave radiation (2.32±0.38 W m −2 (W m −2 ) −1 from the warmer atmosphere, causing a net decrease in the incoming downward surface radiation of −3.56±0.60 W m −2 (W m −2 ) −1 . Despite a reduction in the downward radiation energy, the surface air temperature still increases by 0.25±0.08 K because of less efficient energy dissipation, manifested by reduced surface sensible (−2.88±0.43 W m −2 (W m −2 ) −1 ) and latent heat flux (−1.54±0.27 W m −2 (W m −2 ) −1 ), as well as a decrease in Bowen ratio (−0.20±0.07 (W m −2 ) −1 ). Such reductions of turbulent fluxes can be largely explained by enhanced air stability (0.07±0.02 K (W m −2 ) −1 ), measured as the difference of the potential temperature between 925 hPa and surface, and reduced surface wind speed (−0.05±0.01 m s −1 (W m −2 ) −1 ). The enhanced stability is due to the faster atmospheric warming relative to the surface, whereas the reduced wind speed can be partially explained by enhanced stability and reduced Equator-to-pole atmospheric temperature gradient. These rapid adjustments under BC forcing occur in the lower atmosphere and propagate downward to influence the surface energy redistribution and thus surface temperature response, which is not observed under greenhouse gases or scattering aerosols. Our study provides new insights into the impact of absorbing aerosols on surface energy balance and surface temperature response.

radiative impact↗

Aerodynamic configuration design using response surface methodology analysis

An investigation has been conducted to determine a set of optimal design parameters for a single-stage-to-orbit reentry vehicle. Several configuration geometry parameters which had a large impact on the entry vehicle flying characteristics were selected as design variables: the fuselage fineness ratio, the nose to body length ratio, the nose camber value, the wing planform area scale factor, and the wing location. The optimal geometry parameter values were chosen using a response surface methodology (RSM) technique which allowed for a minimum dry weight configuration design that met a set of aerodynamic performance constraints on the landing speed, and on the subsonic, supersonic, and hypersonic trim and stability levels. The RSM technique utilized, specifically the central composite design method, is presented, along with the general vehicle conceptual design process. Results are presented for an optimized configuration along with several design trade cases.

Engelund, Walter C.↗

Multiobjective Aerodynamic Shape Optimization Using Pareto Differential Evolution and Generalized Response Surface Metamodels

Differential Evolution (DE) is a simple, fast, and robust evolutionary algorithm that has proven effective in determining the global optimum for several difficult single-objective optimization problems. The DE algorithm has been recently extended to multiobjective optimization problem by using a Pareto-based approach. In this paper, a Pareto DE algorithm is applied to multiobjective aerodynamic shape optimization problems that are characterized by computationally expensive objective function evaluations. To improve computational expensive the algorithm is coupled with generalized response surface meta-models based on artificial neural networks. Results are presented for some test optimization problems from the literature to demonstrate the capabilities of the method.

Madavan, Nateri K.↗

A Response Surface Methodology for Bi-Level Integrated System Synthesis (BLISS)

The report describes a new method for optimization of engineering systems such as aerospace vehicles whose design must harmonize a number of subsystems and various physical phenomena, each represented by a separate computer code, e.g., aerodynamics, structures, propulsion, performance, etc. To represent the system internal couplings, the codes receive output from other codes as part of their inputs. The system analysis and optimization task is decomposed into subtasks that can be executed concurrently, each subtask conducted using local state and design variables and holding constant a set of the system-level design variables. The subtasks results are stored in form of the Response Surfaces (RS) fitted in the space of the system-level variables to be used as the subtask surrogates in a system-level optimization whose purpose is to optimize the system objective(s) and to reconcile the system internal couplings. By virtue of decomposition and execution concurrency, the method enables a broad workfront in organization of an engineering project involving a number of specialty groups that might be geographically dispersed, and it exploits the contemporary computing technology of massively concurrent and distributed processing. The report includes a demonstration test case of supersonic business jet design.

Altus, Troy David↗