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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

CONTROLLED END GAS AUTO IGNITION WITH EXHAUST GAS RECIRCULATION ON A STOICHIOMETRIC, SPARK IGNITED, NATURAL GAS ENGINE

The goal of this study is to address fundamental limitations to achieving diesel-like efficiencies in heavy duty on-highway natural gas (NG) engines. Engine knock and misfire are barriers to pathways leading to higher efficiency engines. This study explores enabling technologies for development of high efficiency stoichiometric, spark ignited, natural gas engines. These include design strategies for fast and stable combustion and higher dilution tolerance. Additionally, advanced control methodologies are implemented to maintain stable operation between knock and misfire limits. To implement controlled end- gas autoignition (C-EGAI) strategies a Combustion Intensity Metric (CIM) is used for ignition control with the use of a Woodward large engine control module (LECM). Tests were conducted using a single cylinder, variable compression ratio, cooperative fuel research (CFR) engine with baseline conditions of 900 RPM, engine load of 800 kPa indicated mean effective pressure (IMEP), and stoichiometric air/fuel ratio. Exhaust gas recirculation (EGR) tests were performed using a custom EGR system that simulates a high pressure EGR loop and can provide a range of EGR rates from 0 to 40%. The experimental measurements included the variance of EGR rate, compression ratio, engine speed, IMEP, and CIM. These five variables were optimized through a Modified Box- Benken design Surface Response Method (RSM), with brake efficiency as the merit function. A positive linear correlation between CIM and f-EGAI was identified. Consequently, CIM was used as the feedback control parameter for C-EGAI. As such, implementation of C-EGAI effectively allowed for the utilization of high EGR rates and CRs, controlling combustion between a narrower gap between knock and lean limits. The change from fixed to parametric ignition timing with CIM targeted select values of f-EGAI with an average coefficient of variance (COV) of peak pressure of 5.4. The RSM efficiency optimization concluded with operational conditions of 1080 RPM, 1150 kPa IMEP, 10.55:1 compression ratio, and 17.8% EGR rate with a brake efficiency of 21.3%. At this optimized point of peak performance, a f-EGAI for C-EGAI was observed at 34.1% heat release due to auto ignition, a knock onset crank angle value of 10.3° aTDC and ignition timing of - 24.7° aTDC. This work has demonstrated that combustion at a fixed f-EGAI can be maintained through advanced ignition control of CIM without experiencing heavy knocking events.

Cooperative Fuel Research Engine, Auto Ignition, E↗
NASA NTRS1999

Optimization of a GO2/GH2 Swirl Coaxial Injector Element

An injector optimization methodology, method i, is used to investigate optimal design points for a gaseous oxygen/gaseous hydrogen (GO2/GH2) swirl coaxial injector element. The element is optimized in terms of design variables such as fuel pressure drop, DELTA P(sub f), oxidizer pressure drop, DELTA P(sub 0) combustor length, L(sub comb), and full cone swirl angle, theta, for a given mixture ratio and chamber pressure. Dependent variables such as energy release efficiency, ERE, wall heat flux, Q(sub w) injector heat flux, Q(sub inj), relative combustor weight, W(sub rel), and relative injector cost, C(sub rel), are calculated and then correlated with the design variables. An empirical design methodology is used to generate these responses for 180 combinations of input variables. Method i is then used to generate response surfaces for each dependent variable. Desirability functions based on dependent variable constraints are created and used to facilitate development of composite response surfaces representing some, or all, of the five dependent variables in terms of the input variables. Two examples illustrating the utility and flexibility of method i are discussed in detail. First, joint response surfaces are constructed by sequentially adding dependent variables. Optimum designs are identified after addition of each variable and the effect each variable has on the design is shown. This stepwise demonstration also highlights the importance of including variables such as weight and cost early in the design process. Secondly, using the composite response surface that includes all five dependent variables, unequal weights are assigned to emphasize certain variables relative to others. Here, method i is used to enable objective trade studies on design issues such as component life and thrust to weight ratio.

Tucker, P. Kevin↗
NASA NTRS2001

Optimization of a GO2/GH2 Impinging Injector Element

An injector optimization methodology, method i, is used to investigate optimal design points for a gaseous oxygen/gaseous hydrogen (GO2/GH2) impinging injector element. The unlike impinging element, a fuel-oxidizer- fuel (F-O-F) triplet, is optimized in terms of design variables such as fuel pressure drop, (Delta)P(sub f), oxidizer pressure drop, (Delta)P(sub o), combustor length, L(sub comb), and impingement half-angle, alpha, for a given mixture ratio and chamber pressure. Dependent variables such as energy release efficiency, ERE, wall heat flux, Q(sub w), injector heat flux, Q(sub inj), relative combustor weight, W(sub rel), and relative injector cost, C(sub rel), are calculated and then correlated with the design variables. An empirical design methodology is used to generate these responses for 163 combinations of input variables. Method i is then used to generate response surfaces for each dependent variable. Desirability functions based on dependent variable constraints are created and used to facilitate development of composite response surfaces representing some, or all, of the five dependent variables in terms of the input variables. Three examples illustrating the utility and flexibility of method i are discussed in detail. First, joint response surfaces are constructed by sequentially adding dependent variables. Optimum designs are identified after addition of each variable and the effect each variable has on the design is shown. This stepwise demonstration also highlights the importance of including variables such as weight and cost early in the design process. Secondly, using the composite response surface which includes all five dependent variables, unequal weights are assigned to emphasize certain variables relative to others. Here, method i is used to enable objective trade studies on design issues such as component life and thrust to weight ratio. Finally, specific variable weights are further increased to illustrate the high marginal cost of realizing the last increment of injector performance and thruster weight.

Tucker, P. Kevin↗
NASA NTRS2000

An Optimization-Based Approach to Injector Element Design

An injector optimization methodology, method i, is used to investigate optimal design points for gaseous oxygen/gaseous hydrogen (GO2/GH2) injector elements. A swirl coaxial element and an unlike impinging element (a fuel-oxidizer-fuel triplet) are used to facilitate the study. The elements are optimized in terms of design variables such as fuel pressure drop, APf, oxidizer pressure drop, deltaP(sub f), combustor length, L(sub comb), and full cone swirl angle, theta, (for the swirl element) or impingement half-angle, alpha, (for the impinging element) at a given mixture ratio and chamber pressure. Dependent variables such as energy release efficiency, ERE, wall heat flux, Q(sub w), injector heat flux, Q(sub inj), relative combustor weight, W(sub rel), and relative injector cost, C(sub rel), are calculated and then correlated with the design variables. An empirical design methodology is used to generate these responses for both element types. Method i is then used to generate response surfaces for each dependent variable for both types of elements. Desirability functions based on dependent variable constraints are created and used to facilitate development of composite response surfaces representing the five dependent variables in terms of the input variables. Three examples illustrating the utility and flexibility of method i are discussed in detail for each element type. First, joint response surfaces are constructed by sequentially adding dependent variables. Optimum designs are identified after addition of each variable and the effect each variable has on the element design is illustrated. This stepwise demonstration also highlights the importance of including variables such as weight and cost early in the design process. Secondly, using the composite response surface that includes all five dependent variables, unequal weights are assigned to emphasize certain variables relative to others. Here, method i is used to enable objective trade studies on design issues such as component life and thrust to weight ratio. Finally, combining results from both elements to simulate a trade study, thrust-to-weight trends are illustrated and examined in detail.

Tucker, P. Kevin↗
NASA NTRS2013

Bayesian Revision of Residual Detection Power

This paper addresses some issues with quality assessment and quality assurance in response surface modeling experiments executed in wind tunnels. The role of data volume on quality assurance for response surface models is reviewed. Specific wind tunnel response surface modeling experiments are considered for which apparent discrepancies exist between fit quality expectations based on implemented quality assurance tactics, and the actual fit quality achieved in those experiments. These discrepancies are resolved by using Bayesian inference to account for certain imperfections in the assessment methodology. Estimates of the fraction of out-of-tolerance model predictions based on traditional frequentist methods are revised to account for uncertainty in the residual assessment process. The number of sites in the design space for which residuals are out of tolerance is seen to exceed the number of sites where the model actually fails to fit the data. A method is presented to estimate how much of the design space in inadequately modeled by low-order polynomial approximations to the true but unknown underlying response function.

DeLoach, Richard↗
DOE OSTI2020

Dynamic PRA-Based Estimation of PWR Coping Time Using a Surrogate Model for Accident Tolerant Fuel

In this study, we propose an interpolation-based response surface surrogate methodology to manage a large number of scenarios in dynamic probabilistic risk assessment. It adopts the shape Dynamic Time Warping algorithm to cluster the interpolation neighborhood from time series sample data. The interpolation method was adapted from Taylor Kriging to allow a reduced-order model of the Taylor series. In order to demonstrate its applicability to complex issues in risk assessment for nuclear engineering, an example risk response surface to estimate emergency core cooling system (ECCS) criteria for triplex silicon carbide (SiC) accident-tolerant fuel was constructed. The response surface was exploited to estimate the cumulative failure probability of the fuel cladding structure due to the uncertainties in operator actions and safety systems. The functional failures were assessed based on a combination of individual layer failures computed by coupling Risk Analysis Virtual Environment software with a pressurized water reactor 1000-MW(electric) RELAP5 model and the in-house fuel performance assessment module. Results showed that SiC cladding failure probability spiked less than 1 min after a large-break loss-of- coolant accident whenever the current ECCS criteria for Zircaloy-4 (Zr-4) cladding was used. However, it still provides an increased safety margin of three orders of magnitude compared to Zr-4. This positive margin could be utilized to relax active ECCS requirements by allowing deviations of up to 450 s in its actuation time. The proposed surrogate methodology generated a response surface of SiC cladding failure probability reasonably well, with a significant savings of computation time. This methodology is expected to be useful in the analysis of system response with complex uncertainty sources.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2002

Variable Complexity Optimization of Composite Structures

The use of several levels of modeling in design has been dubbed variable complexity modeling. The work under the grant focused on developing variable complexity modeling strategies with emphasis on response surface techniques. Applications included design of stiffened composite plates for improved damage tolerance, the use of response surfaces for fitting weights obtained by structural optimization, and design against uncertainty using response surface techniques.

Haftka, Raphael T.↗
NASA NTRS2002

Probabilistic Design of a Plate-Like Wing to Meet Flutter and Strength Requirements

An approach is presented for carrying out reliability-based design of a metallic, plate-like wing to meet strength and flutter requirements that are given in terms of risk/reliability. The design problem is to determine the thickness distribution such that wing weight is a minimum and the probability of failure is less than a specified value. Failure is assumed to occur if either the flutter speed is less than a specified allowable or the stress caused by a pressure loading is greater than a specified allowable. Four uncertain quantities are considered: wing thickness, calculated flutter speed, allowable stress, and magnitude of a uniform pressure load. The reliability-based design optimization approach described herein starts with a design obtained using conventional deterministic design optimization with margins on the allowables. Reliability is calculated using Monte Carlo simulation with response surfaces that provide values of stresses and flutter speed. During the reliability-based design optimization, the response surfaces and move limits are coordinated to ensure accuracy of the response surfaces. Studies carried out in the paper show the relationship between reliability and weight and indicate that, for the design problem considered, increases in reliability can be obtained with modest increases in weight.

Stroud, W. Jefferson↗
NASA NTRS2008

Reducing Surface Clutter in Cloud Profiling Radar Data

An algorithm has been devised to reduce ground clutter in the data products of the CloudSat Cloud Profiling Radar (CPR), which is a nadir-looking radar instrument, in orbit around the Earth, that measures power backscattered by clouds as a function of distance from the instrument. Ground clutter contaminates the CPR data in the lowest 1 km of the atmospheric profile, heretofore making it impossible to use CPR data to satisfy the scientific interest in studying clouds and light rainfall at low altitude. The algorithm is based partly on the fact that the CloudSat orbit is such that the geodetic altitude of the CPR varies continuously over a range of approximately 25 km. As the geodetic altitude changes, the radar timing parameters are changed at intervals defined by flight software in order to keep the troposphere inside a data-collection time window. However, within each interval, the surface of the Earth continuously "scans through" (that is, it moves across) a few range bins of the data time window. For each radar profile, only few samples [one for every range-bin increment ((Delta)r = 240 m)] of the surface-clutter signature are available around the range bin in which the peak of surface return is observed, but samples in consecutive radar profiles are offset slightly (by amounts much less than (Delta)r) with respect to each other according to the relative change in geodetic altitude. As a consequence, in a case in which the surface area under examination is homogenous (e.g., an ocean surface), a sequence of consecutive radar profiles of the surface in that area contains samples of the surface response with range resolution (Delta)p much finer than the range-bin increment ((Delta)p << r). Once the high-resolution surface response has thus become available, the profile of surface clutter can be accurately estimated by use of a conventional maximum-correlation scheme: A translated and scaled version of the high-resolution surface response is fitted to the observed low-resolution profile. The translation and scaling factors that optimize the fit in a maximum-correlation sense represent (1) the true position of the surface relative to the sampled surface peak and (2) the magnitude of the surface backscatter. The performance of this algorithm has been tested on CloudSat data acquired over an ocean surface. A preliminary analysis of the test data showed a surface-clutter-rejection ratio over flat surfaces of >10 dB and a reduction of the contaminated altitude over ocean from about 1 km to about 0.5 km (over the ocean). The algorithm has been embedded in CloudSat L1B processing as of Release 04 (July 2007), and the estimated flat surface clutter is removed in L2B-GEOPROF product from the observed profile of reflectivity (see CloudSat product documentation for details and performance at http://www.cloudsat.cira.colostate.edu/ dataSpecs.php?prodid=1).

Tanelli, Simone↗
DOE OSTI2022

Neutral Mode Dominates the Forced Global and Regional Surface Temperature Response in the Past and Future

Abstract Using a large suite of Green's function perturbation experiments, we construct the emergent dynamical operator for surface temperature and extract the most excitable mode. The leading mode turns out to be the most dominant mode excited by the climate forcing of doubling CO 2 , alone capturing 56% of the total spatial variances of the surface temperature response. The pattern of the leading mode is partly organized by the atmospheric annular modes in both hemispheres, and further modulated by the feedbacks from ocean circulations. Though derived from a single model, the leading mode can capture the spatio‐temporal variabilities of the global mean surface temperature over the past century in both the CMIP6 model simulations and observations. Moreover, the leading mode is most efficiently excited by radiative forcings from the midlatitude bands centered around 45°N and 25°S, where most human activities take place, suggesting the potency of human influence on the climate change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2018

Use of the Ames Check Standard Model for the Validation of Wall Interference Corrections

The new check standard model of the NASA Ames 11-ft Transonic Wind Tunnel was chosen for a future validation of the facility's wall interference correction system. The chosen validation approach takes advantage of the fact that test conditions experienced by a large model in the slotted part of the tunnel's test section will change significantly if a subset of the slots is temporarily sealed. Therefore, the model's aerodynamic coefficients have to be recorded, corrected, and compared for two different test section configurations in order to perform the validation. Test section configurations with highly accurate Mach number and dynamic pressure calibrations were selected for the validation. First, the model is tested with all test section slots in open configuration while keeping the model's center of rotation on the tunnel centerline. In the next step, slots on the test section floor are sealed and the model is moved to a new center of rotation that is 33 inches below the tunnel centerline. Then, the original angle of attack sweeps are repeated. Afterwards, wall interference corrections are applied to both test data sets and response surface models of the resulting aerodynamic coefficients in interference-free flow are generated. Finally, the response surface models are used to predict the aerodynamic coefficients for a family of angles of attack while keeping dynamic pressure, Mach number, and Reynolds number constant. The validation is considered successful if the corrected aerodynamic coefficients obtained from the related response surface model pair show good agreement. Residual differences between the corrected coefficient sets will be analyzed as well because they are an indicator of the overall accuracy of the facility's wall interference correction process.

wind tunnel wall interferenc↗
NASA NTRS1975

Early detection of disease program: Evaluation of the cellular immune response

Surfaces of normal, cultured, and mitogen-stimulated mouse lymphoid cells were examined by scanning electron microscopy (SEM). Lymphocytes with smooth, highly villous and intermediate surfaces were observed in cell suspensions from both spleens and thymuses of normal mice and from spleens of congenitally athymic (nude) mice. Several strain-specific surface features were noted, including the spine-like appearance of microvilli on C57B1/6 lymphocytes. Although thymus cell suspensions contained somewhat more smooth cells than did spleen cell preparations, lymphocyte derivation could not be inferred from SEM examination. Studies of cells stimulated with mitogenic agents for thymus-derived lymphocytes (concanavalin A) or for bone marrow-derived lymphocytes (lipopolysaccharide) suggested that, in the mouse, development of a complex villous surface is a general concomitant of lymphocyte activation and transformation.

Criswell, B. S.↗
NASA NTRS1997

Timescales of Land Surface Evapotranspiration Response

Soil and vegetation exert strong control over the evapotranspiration rate, which couples the land surface water and energy balances. A method is presented to quantify the timescale of this surface control using daily general circulation model (GCM) simulation values of evapotranspiration and precipitation. By equating the time history of evaporation efficiency (ratio of actual to potential evapotranspiration) to the convolution of precipitation and a unit kernel (temporal weighting function), response functions are generated that can be used to characterize the timescales of evapotranspiration response for the land surface model (LSM) component of GCMS. The technique is applied to the output of two multiyear simulations of a GCM, one using a Surface-Vegetation-Atmosphere-Transfer (SVAT) scheme and the other a Bucket LSM. The derived response functions show that the Bucket LSM's response is significantly slower than that of the SVAT across the globe. The analysis also shows how the timescales of interception reservoir evaporation, bare soil evaporation, and vegetation transpiration differ within the SVAT LSM.

Scott, Russell↗
DOE OSTI2022

Aggregation in bottom-up vulnerability assessments and equity implications: The case of Jordanian households’ water supply

Bottom-up methods for water resources modeling rely on acceptability thresholds to find, through a response surface, which deeply uncertain futures lead to system failure. They commonly treat water users as aggregate actors, which may preclude analysis of the equity impacts of interventions. This paper explores how aggregation choices for large groups of water users lead to different policy recommendations in response surface assessments. Herein, two aggregation methods with varying parameters are considered: percentile satisfaction targets and generalized mean. A 2-dimensional stress-test assessment across groundwater availability and population is applied to household water supply in Jordan. The study compares six different policies covering supply enhancement and rebalancing, using a country-wide multi-agent model that characterizes households across socioeconomic strata. For different aggregation levels, policies are ordered by their associated robustness index. Results show that aggregation choices may modify response surfaces as much as policy changes and strongly determine policy preference. Modifying allocation rules can substantially reduce the disparity in household vulnerability. Preferences defined by aggregation intervals provide a finer understanding of trade-offs among water users and may improve deliberation over equity under deep uncertainty.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Ocean response to surface heat anomalies

An ocean general circulation model (OGCM) is used to study the response of ocean heat and mass transport to positive and negative heat flux anomalies at the ocean surface. As expected, tropical and low-latitude mixed layers respond rapidly (e-folding time about 50-70 years) to external forcing, while the response of the high-latitude mixed layer, especially the Southern Ocean and northern North Atlantic, is very slow (e-folding time greater than 300 yr). The overall response is faster for negative than positive heat flux anomaly at the surface. The meridional heat transport changes by 15% in the first 50 yr in the southern high latitudes. Surprisingly, for the next 400-500 yr the change is very small. The analysis shows that the meridional mass transport intensifies in response to a negative surface heat flux anomaly but weakens in response to a positive heat flux anomaly. For example, at model year 100 the North Atlantic Deep Water (NADW) is reduced from about 18 Sv to about 10 Sv for the positive heat flux experiment but increased to about 26 Sv for the negative heat flux experiment.

Jiang, Xingjian↗
DOE OSTI2021

Phenotypic plasticity of fungal traits in response to moisture and temperature

Phenotypic plasticity of traits is commonly measured in plants to improve understanding of organismal and ecosystem responses to climate change but is far less studied for microbes. Specifically, decomposer fungi are thought to display high levels of phenotypic plasticity and their functions have important implications for ecosystem dynamics. Assessing the phenotypic plasticity of fungal traits may therefore be important for predicting fungal community response to climate change. Here, we assess the phenotypic plasticity of 15 fungal isolates (12 species) from a Southern California grassland. Fungi were incubated on litter at five moisture levels (ranging from 4–50% water holding capacity) and at five temperatures (ranging from 4–36 °C). After incubation, fungal biomass and activities of four extracellular enzymes (cellobiohydrolase (CBH), β -glucosidase (BG), β -xylosidase (BX), and N-acetyl- β - D -glucosaminidase (NAG)) were measured. We used response surface methodology to determine how fungal phenotypic plasticity differs across the moisture-temperature gradient. We hypothesized that fungal biomass and extracellular enzyme activities would vary with moisture and temperature and that the shape of the response surface would vary between fungal isolates. We further hypothesized that more closely related fungi would show more similar response surfaces across the moisture-temperature gradient. In support of our hypotheses, we found that plasticity differed between fungi along the temperature gradient for fungal biomass and for all the extracellular enzyme activities. Plasticity also differed between fungi along the moisture gradient for BG activity. These differences appear to be caused by variation mainly at the moisture and temperature extremes. We also found that more closely related fungi had more similar extracellular enzymes activities at the highest temperature. Altogether, this evidence suggests that with global warming, fungal biodiversity may become increasingly important as functional traits tend to diverge along phylogenetic lines at higher temperatures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Revisiting the equatorial Pacific sea surface temperature response to global warming

The relative roles of the oceanic and atmospheric processes in the pattern formation of the equatorial Pacific sea surface temperature (SST) response to global warming is investigated using a set of climate model experiments embedded with a novel partial coupling technique. Here, the modeling results show that the SST response experiences a transition from a La Niña-like warming pattern at the initial stage to an El Niño-like warming pattern at the quasi-equilibrium stage. By decomposing anomalous equatorial Pacific SST into atmosphere-forced passive component and ocean dynamics-induced active component, it is found that the SST warming pattern at both stages is entirely induced by its active component. Specifically, the meridional and vertical ocean circulation changes play a dominant role in forming the La Niña-like SST warming pattern at the initial stage, and the zonal and meridional ocean circulation changes are responsible for the formation of the El Niño-like SST warming pattern at the quasi-equilibrium stage. In contrast, the passive SST at both stages is characterized by a zonally uniform warming along the equator, which can be explained by a balance between the cooling effect associated with mean upwelling and the warming effect due to surface passive heat flux change. In addition, this study finds that it is the slowdown of the Pacific subtropical cells during the transition period that controls the evolution of the equatorial SST warming pattern by changing the meridional and vertical ocean heat transports.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Dual-Responsive Macromolecular Surfaces with Binary Patterns

The development of nanopatterned interfaces incorporating multiple stimuli-responsive polymers has great potential for advancing smart sensors and molecular capture devices, yet it remains challenging due to limitations in patterning techniques and “grafting from” methods. Here, we present a binary-patterned surface featuring honeycomb-shaped photothermal polypyrrole (PPy) and thermally responsive poly(ethylene glycol methyl ether acrylate-co-poly(ethylene glycol) methyl ether acrylate) [poly(EGMEA-co-PEGMEA)] brushes within the cavities, fabricated using colloidal lithography and surface-initiated photoinduced electron transfer-reversible addition–fragmentation chain transfer (SI-PET-RAFT) polymerization. The morphology of the surface can be precisely tuned by adjusting the height of each domain. This design integrates unique topographical features with the stimuli-responsiveness of the polymer, enabling the collapse of the copolymer brushes through direct heating or photothermal conversion of PPy under near-IR light. In conclusion, this innovative system offers potential capture-release functionality, providing versatility for diverse applications and enhancing adaptability across different functional scenarios.

Copolymers↗