Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “response surface methods”

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 109 records · Page 6

Lapse-Rate Feeback: The Key to Reconciling TOA and Surface Attributions of Surface Warming

The impact of climate feedbacks on surface warming is conventionally evaluated using a decomposition of the top-of-atmosphere (TOA) energy budget. Alternatively, the climate feedback analysis can also be carried out using the surface energy budget. However, the two perspectives do not provide the same interpretation of process contributions to surface warming, particularly when executing a spatial analysis. The TOA energy budget is equal to the sum of the surface and atmospheric energy budgets. Using the CMIP5 RCP 8.5 model projections, we show that the major discrepancies between the surface and TOA climate feedback attributions of surface warming are due to non-negligible changes in the atmospheric energy budget that differ from their counterparts at the surface. Individual radiative and non-radiative processes cause vertically non-uniform energy flux perturbations, in response to an external forcing, that naturally lead to a vertically non-uniform temperature change. The TOA lapse-rate feedback is the manifestation of these multiple processes that produce a vertically non-uniform warming response such that it accounts for the asymmetry between the changes in the atmospheric and surface energy budgets. Using the climate feedback-response analysis method, we can decompose the lapse-rate feedback into contributions by individual processes. The negative lapse-rate feedback in the tropics is due to greater moist convection and condensational heating along with ocean heat storage. On the other hand, the positive lapse-rate feedback in polar regions is primarily due to the surface albedo and water vapor feedbacks. Combining the process contributions that are hidden within the lapse-rate feedback with their respective direct impacts on the TOA energy budget allows for a very consistent picture of process contributions to surface warming and its inter-model spread as that given by the surface energy budget approach. With the modified approach, both perspectives indicate water vapor feedback is the largest contributor in the tropics, while surface albedo feedback is the greatest contributor in polar regions. Dynamics and ocean heat storage are the main suppressors of surface warming, except over Antarctica. Both perspectives show there is large inter-model uncertainty in the contributions of water vapor, clouds, albedo, and dynamics plus ocean heat storage.

Sergio A Sejas↗

Toward Understanding Wake Vortices and Atmospheric Turbulence Interactions Using Large-Eddy Simulation

The vortices produced by an aircraft in flight are a complex phenomena created from a 'sheet of vorticity' leaving the trailing edge of the aircraft surfaces. This sheet tends to roll-up into two counter-rotating vortices. After a few spans downstream of the aircraft, the roll-up process is complete and the vortex pair may be characterized in a simple manner for modeling purposes. Our research will focus on what happens to these post roll-up vortices in the vicinity of an airport terminal. As the aircraft wake vortices descend, they are transported by the air mass which they are embedded and are decayed by both internal and external processes. In the vicinity of the airport, these external influences are usually due to planetary boundary layer (PBL) turbulence. Using large-eddy simulation (LES), one may simulate a variety of PBL conditions. In the LES method, turbulence is generated in the PBL as a response to surface heat flux, horizontal pressure gradient, wind shear, and/or stratification, and may produce convective or unstably stratified, neutral, or stably stratified PBL's. Each of these PBL types can occur during a typical diurnal cycle of the PBL. Thus it is important to be able to characterize these conditions with the LES method. Once this turbulent environment has been generated, a vortex pair will be introduced and the interactions are observed. The objective is to be able to quantify the PBL turbulence vortex interaction and be able to draw some conclusions of vortex behavior from the various scale interactions. This research is ongoing, and we will focus on what has been accomplished to date and the future direction of this research. We will discuss the model being used, show results that validate its use in the PBL, and present a nested-grid method proposed to analyze the entire PBL and vortex pair simultaneously.

DeCroix, David↗

A Survey of Patch Methods

Patch methods are someshow a response to the fact that surface geometry is local, that is, only small parts of a surface are created at a time. The two categories of patches, transfinite patches and finite dimensional patches are examined and a discussion of trivariate patches is presented.

Barnhill, R. E.↗

Effects of Microstructural Variability on Thermo-Mechanical Properties of a Woven Ceramic Matrix Composite

The objectives of this paper include identifying important architectural parameters that describe the SiC/SiC five-harness satin weave composite and characterizing the statistical distributions and correlations of those parameters from photomicrographs of various cross sections. In addition, realistic artificial cross sections of a 2D representative volume element (RVE) are generated reflecting the variability found in the photomicrographs, which are used to determine the effects of architectural variability on the thermo-mechanical properties. Lastly, preliminary information is obtained on the sensitivity of thermo-mechanical properties to architectural variations. Finite element analysis is used in combination with a response surface and it is shown that the present method is effective in determining the effects of architectural variability on thermo-mechanical properties.

Goldsmith, Marlana B.↗

Dynamic Response of Control Servo System Installed in NAES-Equipped SB2C-5 Airplane (BuAer No. 83135)

Dynamic--response measurements for various conditions of displacement and rate signal input, sensitivity setting, and simulated hinge moment were made of the three control-surface servo systems of an NAES-equipped remote-controlled airplane while on the ground. The basic components of the servo systems are those of the General Electric Company type G-1 autopilot using electrical signal. sources, solenoid-operated valves, and hydraulic pistons. The test procedures and difficulties are discussed, Both frequency and transient-response data, are presented and comparisons are made. The constants describing the servo system, the undamped natural frequency, and the damping ratio, are determined by several methods. The response of the system with the addition of airframe rate signal is calculated. The transfer function of the elevator surface, linkage, and cable system is obtained. The agreement between various methods of measurement and calculation is considered very good. The data are complete enough and in such form that they may be used directly with the frequency-response data of an airplane to predict the stability of the autopilot-airplane combination.

Smaus, Louis H.↗

Hybrid NN/SVM Computational System for Optimizing Designs

A computational method and system based on a hybrid of an artificial neural network (NN) and a support vector machine (SVM) (see figure) has been conceived as a means of maximizing or minimizing an objective function, optionally subject to one or more constraints. Such maximization or minimization could be performed, for example, to optimize solve a data-regression or data-classification problem or to optimize a design associated with a response function. A response function can be considered as a subset of a response surface, which is a surface in a vector space of design and performance parameters. A typical example of a design problem that the method and system can be used to solve is that of an airfoil, for which a response function could be the spatial distribution of pressure over the airfoil. In this example, the response surface would describe the pressure distribution as a function of the operating conditions and the geometric parameters of the airfoil. The use of NNs to analyze physical objects in order to optimize their responses under specified physical conditions is well known. NN analysis is suitable for multidimensional interpolation of data that lack structure and enables the representation and optimization of a succession of numerical solutions of increasing complexity or increasing fidelity to the real world. NN analysis is especially useful in helping to satisfy multiple design objectives. Feedforward NNs can be used to make estimates based on nonlinear mathematical models. One difficulty associated with use of a feedforward NN arises from the need for nonlinear optimization to determine connection weights among input, intermediate, and output variables. It can be very expensive to train an NN in cases in which it is necessary to model large amounts of information. Less widely known (in comparison with NNs) are support vector machines (SVMs), which were originally applied in statistical learning theory. In terms that are necessarily oversimplified to fit the scope of this article, an SVM can be characterized as an algorithm that (1) effects a nonlinear mapping of input vectors into a higher-dimensional feature space and (2) involves a dual formulation of governing equations and constraints. One advantageous feature of the SVM approach is that an objective function (which one seeks to minimize to obtain coefficients that define an SVM mathematical model) is convex, so that unlike in the cases of many NN models, any local minimum of an SVM model is also a global minimum.

Rai, Man Mohan↗

Effects of Spaceflight Relevant Carbon Dioxide Levels on Pathogenesis Related Microbial Characteristics

BACKGROUND Multiple stressors in the spaceflight environment have the potential to alter microbial pathogenesis and virulence characteristics, including microgravity, radiation, and hostile/confined conditions. This combination of factors may lead to a stacking of risks, potentially creating unexpected increase in the risk of infectious disease. One understudied variable is the higher CO2 atmospheric concentration. The current Spaceflight Maximum Allowance Concentration (SMAC) for 24-hour average CO2 is 0.4% (3 mm Hg). These elevated ambient CO2 levels aboard the ISS could potentially influence the diversity and phenotypic responses of the resident microbial communities from both the spacecraft environment (air, surface, water) and crewmembers (gut, nasal, skin microbiomes). METHODS For these studies, we will evaluate the response of Staphylococcus aureus, Streptococcus pneumoniae and Enterobacter aerogenes, which are opportunistic pathogens that are of medical significance and have been or are likely to be found aboard spacecraft. These microorganisms will be cultured in both spaceflight analog and control conditions at CO2levels representing terrestrial atmospheric concentration of 0.04% (0.3 mmHg), ISS elevated concentration of 0.4% (3 mm Hg), and an unexpected elevated concentration of 1 % (7.6 mmHg) to represent a worst case scenario. DISCUSSION Overwhelming evidence from studies performed in food microbiology, marine biology and terrestrial environmental biology suggest altered gene expression, selective bacterial inhibition, increased growth and diversity, and increased antibiotic resistance of bacterial communities individually and in biofilm formation when exposed to increased levels of CO2. While many of these findings are based on levels of CO2 higher than would be found during spaceflight exploration missions, the potential for subtle changes in CO2 levels to exacerbate infectious disease risks warrants an evaluation of pathogen responses to spaceflight CO2 conditions. Our study will provide answers to spaceflight relevant atmospheric conditions, as described in the Risk of Adverse Health Effects Due to Host-Microorganism Interaction (GapMicro 103).

A A Medina-Colorado↗

Graphite/copper alloy interfacial energies determined using the sessile drop method

Graphite surfaces are not wet by pure copper. This lack of wetting is responsible for a debonding phenomenon that was found in continuous graphite fiber/copper matrix composites materials subjected to elevated temperatures. By suitably alloying copper, its capability to wet graphite surfaces can be enhanced. In situ measurements of graphite/copper alloy wetting angles were made using the sessile drop method. Interfacial energy values were calculated based upon these measurements.

Devincent, Sandra M.↗

Analysis of Mars 2020 Entry with Coupled Material Response and CFD

When computing the response of an ablating thermal protection system during entry, the aerothermal environment and material response are traditionally computed independently, with the introduction of a blowing correction term in the material response model to account for the outgassing of char and pyrolysis gases. This study presents an approach, where the pyrolysis blowing gases, calculated within the PATO material response code, are integrated into the DPLR hypersonic CFD code via a blowing boundary condition. This methodology employs an iterative process, whereby the blown pyrolysis gas products from PATO are incorporated into DPLR, refining surface heating predictions. The NEQAIR code is used to compute radiative heating. The method is applied in computing 3D material response of the Mars 2020 entry. When compared with uncoupled material response, the coupled simulations show a lower surface heat flux initially and a higher heat flux at peak heating as shown in Fig. 1. This is in agreement with previous work using a sphere case with environments from MSL.

Heat Transfer↗

Uncertainty Analysis of Decomposing Polyurethane Foam

Sensitivity/uncertainty analyses are necessary to determine where to allocate resources for improved predictions in support of our nation's nuclear safety mission. Yet, sensitivity/uncertainty analyses are not commonly performed on complex combustion models because the calculations are time consuming, CPU intensive, nontrivial exercises that can lead to deceptive results. To illustrate these ideas, a variety of sensitivity/uncertainty analyses were used to determine the uncertainty associated with thermal decomposition of polyurethane foam exposed to high radiative flux boundary conditions. The polyurethane used in this study is a rigid closed-cell foam used as an encapsulant. Related polyurethane binders such as Estane are used in many energetic materials of interest to the JANNAF community. The complex, finite element foam decomposition model used in this study has 25 input parameters that include chemistry, polymer structure, and thermophysical properties. The response variable was selected as the steady-state decomposition front velocity calculated as the derivative of the decomposition front location versus time. An analytical mean value sensitivity/uncertainty (MV) analysis was used to determine the standard deviation by taking numerical derivatives of the response variable with respect to each of the 25 input parameters. Since the response variable is also a derivative, the standard deviation was essentially determined from a second derivative that was extremely sensitive to numerical noise. To minimize the numerical noise, 50-micrometer element dimensions and approximately 1-msec time steps were required to obtain stable uncertainty results. As an alternative method to determine the uncertainty and sensitivity in the decomposition front velocity, surrogate response surfaces were generated for use with a constrained Latin Hypercube Sampling (LHS) technique. Two surrogate response surfaces were investigated: 1) a linear surrogate response surface (LIN) and 2) a quadratic response surface (QUAD). The LHS techniques do not require derivatives of the response variable and are subsequently relatively insensitive to numerical noise. To compare the LIN and QUAD methods to the MV method, a direct LHS analysis (DLHS) was performed using the full grid and timestep resolved finite element model. The surrogate response models (LIN and QUAD) are shown to give acceptable values of the mean and standard deviation when compared to the fully converged DLHS model.

Hobbs, Michael L.↗

Analysis of Mars 2020 Entry with Coupled Material Response and CFD

When computing the response of an ablating thermal protection system during entry, the aerothermal environment and material response are traditionally computed independently, with the introduction of a blowing correction term in the material response model to account for the outgassing of char and pyrolysis gases. This study presents an approach, where the pyrolysis blowing gases, calculated within the PATO material response code [1], are integrated into the DPLR hypersonic CFD code [2] via a blowing boundary condition. This methodology employs an iterative process, whereby the blown pyrolysis gas products from PATO are incorporated into DPLR, refining surface heating predictions. The NEQAIR code is used to compute radiative heating [3]. The method is applied in computing 3D material response of the Mars 2020 entry. When compared with uncoupled material response, the coupled simulations show a lower surface heat flux initially and a higher heat flux at peak heating as shown in Fig. 1. This is in agreement with previous work using a sphere case with environments from MSL [4].

Heat Transfer↗

Acoustic fatigue - A Monte Carlo approach

The present paper reviews the basic concepts of nonlinear structural response to surface flow and acoustic pressure inputs. A time domain solution and the Monte Carlo method are utilized to estimate the required statistics of the nonlinear response. Simulation procedures of multi-dimensional and multi-variate random processes are discussed. Application of these procedures to acoustic fatigue of flight structures are considered. Numerical examples are included to illustrate the applicability of this method for the solution of nonlinear problems.

Vaicaitis, R.↗

Control Surface Interaction Effects of the Active Aeroelastic Wing Wind Tunnel Model

This paper presents results from testing the Active Aeroelastic Wing wind tunnel model in NASA Langley s Transonic Dynamics Tunnel. The wind tunnel test provided an opportunity to study aeroelastic system behavior under combined control surface deflections, testing for control surface interaction effects. Control surface interactions were observed in both static control surface actuation testing and dynamic control surface oscillation testing. The primary method of evaluating interactions was examination of the goodness of the linear superposition assumptions. Responses produced by independently actuating single control surfaces were combined and compared with those produced by simultaneously actuating and oscillating multiple control surfaces. Adjustments to the data were required to isolate the control surface influences. Using dynamic data, the task increases, as both the amplitude and phase have to be considered in the data corrections. The goodness of static linear superposition was examined and analysis of variance was used to evaluate significant factors influencing that goodness. The dynamic data showed interaction effects in both the aerodynamic measurements and the structural measurements.

Heeg, Jennifer↗

Comparison of Response Surface and Kriging Models in the Multidisciplinary Design of an Aerospike Nozzle

The use of response surface models and kriging models are compared for approximating non-random, deterministic computer analyses. After discussing the traditional response surface approach for constructing polynomial models for approximation, kriging is presented as an alternative statistical-based approximation method for the design and analysis of computer experiments. Both approximation methods are applied to the multidisciplinary design and analysis of an aerospike nozzle which consists of a computational fluid dynamics model and a finite element analysis model. Error analysis of the response surface and kriging models is performed along with a graphical comparison of the approximations. Four optimization problems are formulated and solved using both approximation models. While neither approximation technique consistently outperforms the other in this example, the kriging models using only a constant for the underlying global model and a Gaussian correlation function perform as well as the second order polynomial response surface models.

Simpson, Timothy W.↗

Accuracy Assessment of Response Surface Approximations for Supersonic Turbine Design

There is a growing trend to employ CFD tools to supply the necessary information for design optimization of fluid dynamics components/systems. Such results are prone to uncertainties due to reasons including discretization. errors, incomplete convergence of computational procedures, and errors associated with physical models such as turbulence closures. Based on this type of information, gradient-based optimization algorithms often suffer from the noisy calculations, which can seriously compromise the outcome. Similar problems arise from the experimental measurements. Global optimization techniques, such as those based on the response surface (RS) concept are becoming popular in part because they can overcome some of these barriers. However, there are also fundamental issues related to such global optimization technique such as RS. For example, in high dimensional design spaces, typically only a small number of function evaluations are available due to computational and experimental costs. On the other hand, complex features of the design variables do not allow one to model the global characteristics of the design space with simple quadratic polynomials. Consequently a main challenge is to reduce the size of the region where we fit the RS, or make it more accurate in the regions where the optimum is likely to reside. Response Surface techniques using either polynomials or and Neural Network (NN) methods offer designers alternatives to conduct design optimization. The RS technique employs statistical and numerical techniques to establish the relationship between design variables and objective/constraint functions, typically using polynomials. In this study, we aim at addressing issues related to the following questions: (1) How to identify outliers associated with a given RS representation and improve the RS model via appropriate treatments? (2) How to focus on selected design data so that RS can give better performance in regions critical to design optimization? (3) How to combine NN and polynomial techniques for improving the accuracy of the RS model?

Papila, Nilay↗

Carbon-Temperature-Water Change Analysis for Peanut Production Under Climate Change: A Prototype for the AgMIP Coordinated Climate-Crop Modeling Project (C3MP)

Climate change is projected to push the limits of cropping systems and has the potential to disrupt the agricultural sector from local to global scales. This article introduces the Coordinated Climate-Crop Modeling Project (C3MP), an initiative of the Agricultural Model Intercomparison and Improvement Project (AgMIP) to engage a global network of crop modelers to explore the impacts of climate change via an investigation of crop responses to changes in carbon dioxide concentration ([CO2]), temperature, and water. As a demonstration of the C3MP protocols and enabled analyses, we apply the Decision Support System for Agrotechnology Transfer (DSSAT) CROPGRO-Peanut crop model for Henry County, Alabama, to evaluate responses to the range of plausible [CO2], temperature changes, and precipitation changes projected by climate models out to the end of the 21st century. These sensitivity tests are used to derive crop model emulators that estimate changes in mean yield and the coefficient of variation for seasonal yields across a broad range of climate conditions, reproducing mean yields from sensitivity test simulations with deviations of ca. 2% for rain-fed conditions. We apply these statistical emulators to investigate how peanuts respond to projections from various global climate models, time periods, and emissions scenarios, finding a robust projection of modest (<10%) median yield losses in the middle of the 21st century accelerating to more severe (>20%) losses and larger uncertainty at the end of the century under the more severe representative concentration pathway (RCP8.5). This projection is not substantially altered by the selection of the AgMERRA global gridded climate dataset rather than the local historical observations, differences between the Third and Fifth Coupled Model Intercomparison Project (CMIP3 and CMIP5), or the use of the delta method of climate impacts analysis rather than the C3MP impacts response surface and emulator approach.

climate change↗

Shear-Sensitive Liquid Crystal Coating Method: Surface-Inclination Effects on Shear Vector Measurements

The shear-sensitive liquid crystal coating (SSLCC) method is an image-based technique for both visualizing dynamic surface-flow phenomena, such as transition and separation, and for measuring the continuous shear-stress vector distribution acting on an aerodynamic surface. Under proper lighting and viewing conditions (discussed below), the coating changes color in response to an applied aerodynamic shear. This color-change response is continuous and reversible, with a response time of milliseconds, and is a function of both the shear magnitude and the shear vector orientation relative to the observer. The liquid crystal phase of matter is a weakly-ordered, viscous, non-Newtonian fluid state that exists between the nonuniform liquid phase and the ordered solid phase of certain organic compounds. Cholesteric liquid crystal compounds possess a helical molecular arrangement that selectively scatters white light, incident along the helical axis, as a three-dimensional spectrum. This property is linked to the helical pitch length, which is within the range of wavelengths in the visible spectrum. The pitch length, and hence the wavelength of the scattered light, is influenced by shear stress normal to the helical axis. This unique optical property produces a measurable color change in response to an applied shearing force. The full-surface shear stress vector measurement method, developed at NASA-Ames, is schematically illustrated. As with the visualization method, the coated test surface is illuminated from the normal direction with white light and the camera is positioned at an above-plane view angle of approximately 30 deg. Experiments have been initiated at NASA Ames to begin the process of quantifying surface-inclination (surface-curvature) effects on shear vector measurement accuracy. In preliminary experiments, surface-inclination angles theta(sub x), theta(sub y) of 0, +/-5, +/-10, and +/-15 deg were employed. In this arrangement, white-light illumination was positioned normal to the untilted test surface, and the camera above-plane view angle was set at 30 deg relative to the untilted test surface. As can be seen, vector-aligned lambda(sub d) values showed no dependence on theta(sub x) or theta(sub y) for absolute values of these tilt angles is less than or equal to 15 deg. Acquisition and analyses of full-surface color images are presently underway to definitively document the insensitivity limits of the shear vector measurement methodology to surface-slope variations.

Reda, Daniel C.↗

Measurement of Flaw Size From Thermographic Data

Simple methods for reducing the pulsed thermographic responses of delaminations tend to overestimate the size of the delamination, since the heat diffuses in the plane parallel to the surface. The result is a temperature profile over the delamination which is larger than the delamination size. A variational approach is presented for reducing the thermographic data to produce an estimated size for a flaw that is much closer to the true size of the delamination. The method is based on an estimate for the thermal response that is a convolution of a Gaussian kernel with the shape of the flaw. The size is determined from both the temporal and spatial thermal response of the exterior surface above the delamination and constraints on the length of the contour surrounding the delamination. Examples of the application of the technique to simulation and experimental data are presented to investigate the limitations of the technique.

Winfree, William P.↗