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 217 records · Page 12

How High is that Dune? A Comparison of Methods Used to Constrain the Morphometry of Aeolian Bedforms on Mars

Methods traditionally used to estimate the relative height of surface features on Mars include: photoclinometry, shadow length and stereography. The MOLA data set enables a more accurate assessment of the surface topography of Mars. However, many small-scale aeolian bedforms remain below the sample resolution of the MOLA data set. In response to this a number of research teams have adopted and refined existing methods and applied them to high resolution (2-6 m/pixel) narrow angle MOC satellite images. Collectively, the methods provide data on a range of morphometric parameters (many not previously available for dunes on Mars). These include dune height, width, length, surface area, volume, longitudinal and cross profiles). This data will facilitate a more accurate analysis of aeolian bedforms on Mars. In this paper we undertake a comparative analysis of methods used to determine the height of aeolian dunes and ripples.

Bourke, M.↗

Understanding the hydrological response of a headwater-dominated catchment by analysis of distributed surface–subsurface interactions

Abstract We computationally explore the relationship between surface–subsurface exchange and hydrological response in a headwater-dominated high elevation, mountainous catchment in East River Watershed, Colorado, USA. In order to isolate the effect of surface–subsurface exchange on the hydrological response, we compare three model variations that differ only in soil permeability. Traditional methods of hydrograph analysis that have been developed for headwater catchments may fail to properly characterize catchments, where catchment response is tightly coupled to headwater inflow. Analyzing the spatially distributed hydrological response of such catchments gives additional information on the catchment functioning. Thus, we compute hydrographs, hydrological indices, and spatio-temporal distributions of hydrological variables. The indices and distributions are then linked to the hydrograph at the outlet of the catchment. Our results show that changes in the surface–subsurface exchange fluxes trigger different flow regimes, connectivity dynamics, and runoff generation mechanisms inside the catchment, and hence, affect the distributed hydrological response. Further, changes in surface–subsurface exchange rates lead to a nonlinear change in the degree of connectivity—quantified through the number of disconnected clusters of ponding water—in the catchment. Although the runoff formation in the catchment changes significantly, these changes do not significantly alter the aggregated streamflow hydrograph. This hints at a crucial gap in our ability to infer catchment function from aggregated signatures. We show that while these changes in distributed hydrological response may not always be observable through aggregated hydrological signatures, they can be quantified through the use of indices of connectivity.

54 ENVIRONMENTAL SCIENCES↗

Transfer function models for using empirical and physics-based simulation signal response data

In many situations, real or induced flaws such as tight cracks with known morphology cannot be manufactured in part geometry specimens or in real parts. Typically, surface fatigue cracks are manufactured in simple geometry specimens such as flat plates, dog-bone shaped flat or cylindrical specimens. If a nondestructive evaluation (NDE)technique is required to provide a reliably detectable flaw size, denoted asa90/95, for induced flaws in a part, then a direct method for qualifying the NDE procedure is to use appropriate induced flaw specimens and perform NDE procedure demonstration on the specimens. Probability of detection (POD)analysis of the empirical data may provide estimation of a90/95. This approach is described as direct POD demonstration testing, which may follow guidelines of MIL-HDBK-1823. This paper considers a case, where embedded tight cracklike induced flaws are to be detected reliably using a signal response based NDE procedure. Here, it is assumed that it is not practical to make surface or embedded induced flaw specimens in part geometry or configuration. Therefore, a direct POD demonstration testing cannot be undertaken. It is also assumed that simulation of signal response is possible for both surface and embedded induced flaws in part geometry specimens using a physics-based model. The proposed approach for NDE procedure qualification uses artificial flaws in simple geometry and part geometry specimens, and induced flaws in the same type of simple geometry specimens. Signal response data is taken on all sets of artificial and induced flaws in simple geometry and part geometry specimens. Moreover, simulated signal response data is generated for surface and embedded flaws. Thus, a case of five signal response versus flaw size datasets is considered. Three of the datasets are empirical and two datasets are physics model-based simulation datasets. A method of devising and using transfer function calculation dataset blocks to estimate the either the reliably detectable flaw size or the demonstration flaw size is provided.

Nondestructive evaluation↗

Chemometric optimisation of pressurised liquid extraction for the determination of alliin and S‐allyl‐cysteine in giant garlic ( Allium ampeloprasum L.) by liquid chromatography tandem mass spectrometry

Abstract Introduction Giant garlic is a functional food that contains different kinds of bioactive molecules with beneficial effects on chronic noncommunicable diseases like diabetes and cardiovascular conditions. Considering biosynthesis pathways, abundance, and biological activity, alliin and S‐allyl‐cysteine were used as chemical markers of organosulphur compounds present in giant garlic. Objective To establish a chemometric optimisation of pressurised liquid extraction for the determination of alliin and S‐allyl‐cysteine in giant garlic by liquid chromatography tandem mass spectrometry (LC–MS/MS). Methodology Samples were blanched (ca. 90°C for 10 min) to inactivate alliinase and γ‐glutamyl transpeptidase enzymes and then freeze‐dried. Chemometric optimisation was performed via response surface methodology based on central composite design (CCD). Organosulphur compound yields were determined applying a validated LC–MS/MS method in multiple reaction monitoring (MRM) mode using the following transitions: for alliin m/z 178 → 74 and for S‐allyl‐cysteine m/z 162 → 41. Results According to CCD results, under constant conditions of pressure (1500 psi) and time (20 min), the optimal conditions for pressurised liquid extraction of alliin and S‐allyl‐cysteine were 70.75 and 68.97% v/v of ethanol in extraction solvent and 76.45 and 98.88°C as extraction temperature, respectively. Multiple response optimisation for the simultaneous extraction of both organosulphur compounds was established via desirability function. Under these conditions, 2.70 ± 0.27 mg g −1 dry weight (DW) of alliin and 2.79 ± 0.22 mg g −1 DW of S‐allyl‐cysteine were extracted. Conclusions These results clearly demonstrated that pressurised liquid extraction is an efficient green technique to extract bioactive organosulphur compounds from giant garlic. Extraction yields were significantly ( p < 0.05) higher than those obtained with conventional ultra‐turrax extraction.

Peterssen‐Fonseca, Darlene↗

Secular rotational motions and the mechanical structure of a dynamical viscoelastic earth

A survey is presented of analytical methods for computing the linear responses of the rotational axis of a layered viscoelastic earth to surface loading. Theoretical research in this area is first summarized, and the differences between the mechanical boundary conditions to be applied at the interface separating the upper and lower mantles for an adiabatically and chemically stratified mantle are discussed. Some examples of polar wander and secular variation of the spin rate from glacial excitation are presented for various types of chemical and viscosity stratifications. The effects of an artificial density jump at the base of the lithosphere in models are examined, and certain issues concerning the fluid tidal Love number for different types of density stratification are addressed. The meaning of effective plate thickness over geological time scales for rotational dynamics is discussed.

Yuen, D. A.↗

Design Optimization of Composite Structures under Uncertainty

Design optimization under uncertainty is computationally expensive and is also challenging in terms of alternative formulation. The work under the grant focused on developing methods for design against uncertainty that are applicable to composite structural design 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 simultaneous design of structure and inspection periods for fail-safe structures.

Haftka, Raphael T.↗

Detecting Volcanic Co2 Emissions Hidden in Tropical Volcanic Forests Using Fixed-Wing sUAS

CO2 emissions are among the earliest indicators of subsurface state change and reactivation of volcanic systems and these deep signals reflect in faint diffuse emissions variations on volcanic flanks. Monitoring vast areas of rainforests covering tropical volcanoes is challenging, and space-borne sensors like OCO-2 are not sensitive enough. Fortunately, trees exposed to mild enhancements of CO2 may experience fertilization and build excess CO2 into excess biomass (Cawse-Nicholson, et al. 2018, Biogeosciences). They also process excess CO2 through photosynthesis, measurable by fluorescence (Bogue et al., 2019, Biogeosciences). Furthermore, spaceborne remote sensing instruments have demonstrated increasing NDVI months to years before eruptions, unexplained by other factors (e.g., Houlie et al. 2006, EPSL; Seiler et al. 2017, PLoS One). By providing a means to map large areas of above-canopy CO2 variations over tropical rainforest on the flanks of actively degassing volcanoes, we enable to interpret the signals plants provide in response to excess CO2 with airborne and spaceborne hyperspectral remote sensing methods, including a possible future Surface Biology and Geology (SBG) NASA satellite mission.As a critical step to mature and test this concept, we deployed small unmanned aerial systems (sUAS) using a fixed-wing sUAS designed for autonomous operations and long endurance in extreme environments, without causing significant disturbance of the sampled air. We integrated an in-situ sensor to record the CO2 enhancement field above the forest canopy. Test flight results on the flanks of Turrialba volcano in Costa Rica above forest canopies covering known moderate gas seeps demonstrate the sUAS-borne detection and mapping capabilities of above-canopy elevated CO2 gradients. The strong detection capabilities and high detector signal stability resulted from key system design elements including RF shielding, mechanical stabilization, and calibration procedures. This highly robust system is readily applied for diffuse volcanic CO2 emission studies on active volcanoes covered by dense vegetation.

Volcanic↗

Dynamic Modeling using Output-Error Parameter Estimation based on Frequency Responses Estimated with Multisine Inputs

A method is developed for estimating model parameters, such as nondimensional stability and control derivatives, by fitting transfer function or state-space models to empirical frequency response data using the output-error approach. The frequency response data were computed using Fourier transforms of measured input and output data. The control surfaces were excited with periodic multisine inputs which facilitated time-ecient estimation of multiple-input multiple-output frequency responses. The method was applied to lateral data from a nonlinear flight dynamics simulation of the F-16 aircraft, and to longitudinal data from multiple repeated flight test maneuvers of the NASA T-2 subscale aircraft. Results using simulation data showed the frequency response method compared well to other standard methods for parameter estimation. In addition to including all the available inputs, outputs, and harmonic frequencies in the estimation, relatively small subsets of the measured data could also be used to focus on identifying specific parts of the model. Results from flight test data showed that parameter estimates and uncertainties determined from repeated maneuvers were accurate and in statistical agreement with each other.

Grauer, Jared A.↗

Tropospheric Response to Estimated Spectrally Discriminated Solar Forcing Over the Past 500 Years

The GISS Global Climate Middle Atmosphere Model (GCMAM) is used to investigate the effect of estimated solar irradiance changes on climate for the past 500 years. This model is employed to allow the impact of UV variations on the stratosphere to affect the troposphere via wave-mean flow interactions. Multiple experiments are done with only a total solar irradiance change (peaking at 0.2 percent from the Maunder Minimum to today); with estimated spectrally-varying irradiance changes (i.e., peak changes of 0.7 percent in the UV, 0.2 percent in the visible and near IR; and 0.07 percent in the IR greater than 1 micron); and the spectrally-varying changes in conjunction with model calculated ozone responses in the stratosphere. Results of the varying temperature patterns and radiation response will be discussed. Of interest is whether the different methods of forcing the solar-induced climate change produce different spatial surface temperature signatures, particularly ones that can be differentiated from greenhouse gas warming. In preliminary tests, spectrally-varying solar forcing with induced ozone changes for solar maximum minus solar minimum conditions results in a temperature signal that is primarily at high latitudes.The high latitude response arises due to solar/ozone-induced alterations in the stratospheric wind field that affect planetary wave propagation from the troposphere, and alter tropospheric advection patterns. In contrast, forcing by total solar irradiance changes produces significant response at low and subtropical latitudes as well, driven by water vapor and cloud feedbacks to the radiative perturbation.

Rind, David↗

Viscosity Measurement using Drop Coalescence in Microgravity

We present in here details of a new method, using drop coalescence, for application in microgravity environment for determining the viscosity of highly viscous undercooled liquids. The method has the advantage of eliminating heterogeneous nucleation at container walls caused by crystallization of undercooled liquids during processing. Also, due to the rapidity of the measurement, homogeneous nucleation would be avoided. The technique relies on both a highly accurate solution to the Navier-Stokes equations as well as on data gathered from experiments conducted in near zero gravity environment. The liquid viscosity is determined by allowing the computed free surface shape relaxation time to be adjusted in response to the measured free surface velocity of two coalescing drops. Results are presented from two validation experiments of the method which were conducted recently on board the NASA KC-135 aircraft. In these tests the viscosity of a highly viscous liquid, such as glycerine at different temperatures, was determined to reasonable accuracy using the liquid coalescence method. The experiments measured the free surface velocity of two glycerine drops coalescing under the action of surface tension alone in low gravity environment using high speed photography. The free surface velocity was then compared with the computed values obtained from different viscosity values. The results of these experiments were found to agree reasonably well with the calculated values.

Antar, Basil N.↗

Fatigue Reliability of Gas Turbine Engine Structures

The results of an investigation are described for fatigue reliability in engine structures. The description consists of two parts. Part 1 is for method development. Part 2 is a specific case study. In Part 1, the essential concepts and practical approaches to damage tolerance design in the gas turbine industry are summarized. These have evolved over the years in response to flight safety certification requirements. The effect of Non-Destructive Evaluation (NDE) methods on these methods is also reviewed. Assessment methods based on probabilistic fracture mechanics, with regard to both crack initiation and crack growth, are outlined. Limit state modeling techniques from structural reliability theory are shown to be appropriate for application to this problem, for both individual failure mode and system-level assessment. In Part 2, the results of a case study for the high pressure turbine of a turboprop engine are described. The response surface approach is used to construct a fatigue performance function. This performance function is used with the First Order Reliability Method (FORM) to determine the probability of failure and the sensitivity of the fatigue life to the engine parameters for the first stage disk rim of the two stage turbine. A hybrid combination of regression and Monte Carlo simulation is to use incorporate time dependent random variables. System reliability is used to determine the system probability of failure, and the sensitivity of the system fatigue life to the engine parameters of the high pressure turbine. 'ne variation in the primary hot gas and secondary cooling air, the uncertainty of the complex mission loading, and the scatter in the material data are considered.

Cruse, Thomas A.↗

Laser Doppler vibrometry for piezoelectric coefficient ($d_{33}$) measurements in irradiated aluminum nitride

Sensors used for experiments in advanced reactors must survive in harsh environments. Few material systems can be used to construct sensors viable for extreme conditions. Aluminum nitride (AlN) is one such material because it has high thermal stability and radiation resistance and sustains good piezoelectric and dielectric properties at high temperatures. In this work, the piezoelectric coefficient $d_{33}$ of the AlN single-crystal with thermal and irradiation damage was investigated using an indirect method with a laser Doppler vibrometer (LDV). Surface electrodes were deposited on the AlN samples, and the vibration response of the samples to an applied voltage was monitored using the LDV as a function of the excitation frequency. The $d_{33}$ estimation was based on the excitation voltage and the thickness-mode displacement extracted from the LDV measurements. Further, six AlN substrates were irradiated with 8 MeV Al 2+ at three fluences (10 15 , 10 16 , and 10 17 ions/cm 2 ) and two temperatures (300 °C and 500 °C). The $d_{33}$ for the six irradiated samples and one pristine sample were measured, and the measurement uncertainty was estimated based on five repeated tests. All samples were also measured by a commercial piezometer for comparison. The experimental results demonstrate that the piezoelectric coefficients obtained by LDV were about 0.8–1.16 pm/V lower than those obtained by the piezometer. With the compensation of the clamping effect, the corrected LDV values are similar to the piezometer results. Both show similar trends in all samples, which validates the feasibility of the proposed method for $d_{33}$ measurement. Based on the LDV results, the irradiated samples show a 12%–22% decrease in $d_{33}$ compared with the pristine samples. The samples irradiated under the same fluence at a higher temperature (500 °C) demonstrated a lower $d_{33}$ than those at 300 °C. The effect of the retro-reflective tape, sample temperature, and sample size on the $d_{33}$ measurement were also studied.

36 MATERIALS SCIENCE↗

Functional stimuli-responsive polymers on micro- and nano-patterned interfaces

Micro- and nano-patterned surfaces offer precise control over morphology and chemical composition, enhancing the stability, durability, and functionality of coating materials. When combined with stimuli-responsive polymers, these surfaces gain dynamic adaptability, enabling reversible binding, reusable sensing, and selective molecular capture. Furthermore, while recent review articles have explored various aspects of stimuli-responsive materials, from hydrogel patterns for bioanalytical applications to shape-morphing hydrogels for soft robotics and sensors, a comprehensive review focused on the integration of smart polymers with micro- or nano-patterned interfaces remains absent. This review addresses key surface patterning techniques, including soft lithography, colloidal lithography, and polymer brush photolithography, as well as advances in surface-initiated polymerization methods, such as surface-initiated controlled radical polymerization (SI-CRP). In addition, we discuss recent progress in integrating stimuli-responsive polymers with patterned surfaces to create advanced, functional materials.

Colloidal lithography↗

initMIP-Antarctica: an Ice Sheet Model Initialization Experiment of ISMIP6

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

ice sheet model↗

Efficient learning of accurate surrogates for simulations of complex systems

Machine learning methods are increasingly deployed to construct surrogate models for complex physical systems at a reduced computational cost. However, the predictive capability of these surrogates degrades in the presence of noisy, sparse or dynamic data. Here, we introduce an online learning method empowered by optimizer-driven sampling that has two advantages over current approaches: it ensures that all local extrema (including endpoints) of the model response surface are included in the training data, and it employs a continuous validation and update process in which surrogates undergo retraining when their performance falls below a validity threshold. We find, using benchmark functions, that optimizer-directed sampling generally outperforms traditional sampling methods in terms of accuracy around local extrema even when the scoring metric is biased towards assessing overall accuracy. Finally, the application to dense nuclear matter demonstrates that highly accurate surrogates for a nuclear equation-of-state model can be reliably autogenerated from expensive calculations using few model evaluations.

79 ASTRONOMY AND ASTROPHYSICS↗

Efficient Learning of Accurate Surrogates for Simulations of Complex Systems

Machine learning methods are increasingly deployed to construct surrogate models for complex physical systems at a reduced computational cost. However, the predictive capability of these surrogates degrades in the presence of noisy, sparse, or dynamic data. We introduce an online learning method empowered by optimizer-driven sampling that has two advantages over current approaches: it ensures that all local extrema (including endpoints) of the model response surface are included in the training data and it employs a continuous validation and update process in which surrogates undergo retraining when their performance falls below a validity threshold. We find, using benchmark functions, that optimizer-directed sampling generally outperforms traditional sampling methods in terms of accuracy around local extrema even when the scoring metric is biased towards assessing overall accuracy. Finally, the application to dense nuclear matter demonstrates that highly accurate surrogates for a nuclear equation of state model can be reliably auto-generated from expensive calculations using few model evaluations.

97 MATHEMATICS AND COMPUTING↗

A new response surface approach for structural reliability analysis

This paper describes a new approach for computing structural reliability by post-processing previously computed probabilistic results for stress and strength. The objective is to provide an accurate method whereby independent probabilistic analyses for stress and strength functions can be performed independently and combined at a later time to compute probability of failure. The method provides a capability for testing different strength measures without the need for re-computing the probabilistic stress response. The proposed approach takes full account of the basic random variables effecting both stress and strength, and the failure region in the variable space identified during separate stress/strength probabilistic analyses. A simple closed-form example and a more complex analysis of a turbine blade subject to creep rupture is used to illustrate the method.

Thacker, B. H.↗

Calibrating hypersonic turbulence flow models with the HIFiRE-1 experiment using data-driven machine-learned models.

In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model’s input parameters. Here we develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model’s response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.

42 ENGINEERING↗