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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Evaluation of the Surface Reflectance Long-Term Data Record from AVHRR over Multiple Land Surface Types

In this work, we evaluate the performance of the AVHRR surface reflectance LTDR V4 using Landsat-5 Thematic Mapper (TM5) Collection-1 surface reflectance data over 440 globally distributed sites, which give a representative set of land surface types and climate conditions. Surface reflectance anisotropy effects were normalized using the VJB method, and spectral response differences were accounted for with spectral band adjustment factors (SBAFs) computed as a function of the Normalized Difference Vegetation Index (NDVI) from a set of over 100,000 hyperspectral spectra from Hyperion. The performance of the AVHRR record is reported in terms of the accuracy, precision, and uncertainty, as compared to the specifications of the reference product. Results show that the AVHRR record performance is close to the 0.071ρ+0.0071 specification defined in the original TM5 product evaluation.

LTDR↗
DOE OSTI2022

Estimating epistemic uncertainty in soil parameters for nonlinear site response analyses: Introducing the Latin Hypercube Sampling technique

This study quantifies the effects of epistemic uncertainty in soil parameters on nonlinear (NL) site response analysis (SRA) results, validated against the data recorded at a well-instrumented geotechnical downhole array located in Japan. To this end, a one-dimensional soil column model of the Service Hall Array (SHA) near the Kashiwazaki-Kariwa Nuclear Power Plant (KKNPP) is developed using the finite element (FE) program LS-DYNA. The dynamic stress–strain relationship is characterized by a modified two-stage hyperbolic (MTH) NL backbone curve formulation capable of capturing soil behavior at both small- and large-shear strains. The model is then validated against the ground motion recordings to capture the model bias. The uncertainties associated with the shear-wave velocity profile (a small-strain soil property) and soil shear strength (a large-strain soil property) are incorporated in NL SRA to quantify their separate and joint randomization effects on the results. This study proposes using the Latin Hypercube Sampling (LHS) method as an efficient alternative to commonly used methods, such as Standard Monte Carlo (SMC), to account for uncertainty propagation in such reliability analysis. Both low-intensity and design-level records from the recordings at the SHA are applied to study the contribution of the small- and large-strain NL dynamic soil properties. Results from 46,200 NL FE analyses (23,100 per input motion) are presented. Measured and predicted site response, using recorded ground motions at this downhole array, is compared to assess the significance of soil parameter uncertainty on the observed ground motion dispersions. It is demonstrated that increasing the number of soft realizations and implementing higher level earthquake intensity lead to higher ground motion dispersion. Unlike past studies in randomization of Vs profiles with the SMC method, the LHS method is shown to have no significant effect on the predicted median surface response spectra and amplification factors (AFs) for this case study.

Engineering↗
NASA NTRSDate not supplied

Quantifying Changes in the Arctic Shortwave Cloud Radiative Effects

The shortwave cloud radiative effect (SWCRE) is important on the Arctic surface radiation budget and the major source of inter-model spread in predictions of Arctic climate. To better understand the individual contributions of various radiative processes to changes in SWCRE, the paper presents the use of the extended APRP (Atmospheric Radiative Perturbation Potential) method. This involves adding the absorptivity for the upward beam and considering differences in reflectivity between upward and downward beams, as well as analyzing the cloud masking effect resulting from changes in surface albedo in more detail. Using data from the CMIP5 and CMIP6 climate models, the study decomposes the SWCRE over the Arctic surface and analyzes inter-model differences in quadrupled CO2 simulations. The study takes into account the fact that the response of SWCRE to Arctic warming is influenced by changes in surface albedo, cloud amount, and cloud microphysics. Results show that in the sunlight season, the reduction in surface albedo associated with sea ice loss is directly linked to strong negative SWCRE, which explains the considerable model discrepancy. Arctic clouds can hinder the positive surface albedo feedback by changing the albedo in two ways: (1) decreasing incoming shortwave radiation due to cloud reflection and (2) by decreasing the shortwave reaching the surface after being reflected by clouds. In addition, increased (decreased) cloud amount and cloud liquid water are shown to be less (more) incoming shortwave fluxes at the surface, but not dominating factors to the Arctic surface radiation budget and its inter-model variation. Overall, the extended APRP method offers a useful tool for analyzing the complex interactions between clouds and radiative process, reasonably decomposes the individual SWCRE responses at the Arctic surface, and emphasizes that considering not only the cloud amount or its properties, but also surface albedo change is critical for the prediction of SWCRE on the Arctic surface.

Arctic cloud↗
NASA NTRS1992

A boundary integral method for an inverse problem in thermal imaging

An inverse problem in thermal imaging involving the recovery of a void in a material from its surface temperature response to external heating is examined. Uniqueness and continuous dependence results for the inverse problem are demonstrated, and a numerical method for its solution is developed. This method is based on an optimization approach, coupled with a boundary integral equation formulation of the forward heat conduction problem. Some convergence results for the method are proved, and several examples are presented using computationally generated data.

Bryan, Kurt↗
DOE OSTI2020

UQ Toolkit v 2.0

The Uncertainty Quantification (UQ) Toolkit is a software library for the characterizaton and propagation of uncertainties in computational models. For the characterization of uncertainties, Bayesian inference tools are provided to infer uncertain model parameters, as well as Bayesian compressive sensing methods for discovering sparse representations of high-dimensional input-output response surfaces, and also Karhunen-Loève expansions for representing stochastic processes. Uncertain parameters are treated as random variables and represented with Polynomial Chaos expansions (PCEs). The library implements several spectral basis function types (e.g. Hermite basis functions in terms of Gaussian random variables or Legendre basis functions in terms of uniform random variables) that can be used to represent random variables with PCEs. For propagation of uncertainty, tools are provided to propagate PCEs that describe the input uncertainty through the computational model using either intrusive methods (Galerkin projection of equations onto basis functions) or non-intrusive methods (perform deterministic operation at sampled values of the random values and project the obtained results onto basis functions).

Safta, Cosmin↗
NASA NTRS2014

Determination of Flaw Size from Thermographic Data

Conventional methods for reducing the pulsed thermographic responses of delaminations tend to overestimate the size of the flaw. Since the heat diffuses in the plane parallel to the surface, the resulting temperature profile over the flaw is larger than the flaw. A variational method is presented for reducing the thermographic data to produce an estimated size for the flaw that is much closer to the true size of the flaw. The size is determined from the spatial thermal response of the exterior surface above the flaw and a constraint on the length of the contour surrounding the flaw. The technique is applied to experimental data acquired on a flat bottom hole composite specimen.

Winfree, William P.↗
NASA NTRS1972

Excitation, response, and fatigue life estimation methods for the structural design of externally blown flaps

Means for predicting the fluctuating pressures acting on externally blown flap surfaces are developed on the basis of generalizations derived from non-dimensionalized empirical data. Approaches for estimation of the fatigue lives of skin-stringer and honeycomb-core sandwich flap structures are derived from vibration response analyses and panel fatigue data. Approximate expressions for fluctuating pressures, structural response, and fatigue life are combined to reveal the important parametric dependences. The two-dimensional equations of motion of multi-element flap systems are derived in general form, so that they can be specialized readily for any particular system. An introduction is presented of an approach to characterizing the excitation pressures and structural responses which makes use of space-time spectral concepts and promises to provide useful insights, as well as experimental and analytical savings.

Ungar, E. E.↗
NASA NTRS1984

Turbine Engine Hot Section Technology, 1984

Presentations were made concerning the hot section environment and behavior of combustion liners, turbine blades, and waves. The presentations were divided into six sessions: instrumentation, combustion, turbine heat transfer, structural analysis, fatigue and fracture, and surface properties. The principal objective of each session was to disseminate research results to date, along with future plans. Topics discussed included modeling of thermal and fluid flow phenomena, structural analysis, fatigue and fracture, surface protective coatings, constitutive behavior, stress-strain response, and life prediction methods.

Source record↗
NASA NTRS2018

Ascent Aerodynamic Force and Moment Database Development for the Space Launch System

The Space Launch System Aerodynamics Task Team is responsible for delivering aerodynamic force and moment databases from liftoff through ascent until the rocket leaves the Earth’s atmosphere. The process for developing the ascent portion of this database is described in the current paper. The data used to develop the database were generated using a combination of wind tunnel testing and CFD simulations. The details of the wind tunnel testing performed at the NASA Ames Unitary Plan Wind Tunnel and CFD simulations performed using FUN3D at wind tunnel and flight conditions are discussed, and comparisons of these data sets are provided. The methods used for converting the source data into the final database response surfaces with corresponding uncertainty are also detailed.

Shea, Patrick R.↗
NASA NTRS2014

A Modular Approach to Model Oscillating Control Surfaces Using Navier Stokes Equations

The use of active controls for rotorcraft is becoming more important for modern aerospace configurations. Efforts to reduce the vibrations of helicopter blades with use of active-controls are in progress. Modeling oscillating control surfaces using the linear aerodynamics theory is well established. However, higher-fidelity methods are needed to account for nonlinear effects, such as those that occur in transonic flow. The aeroelastic responses of a wing with an oscillating control surface, computed using the transonic small perturbation (TSP) theory, have been shown to cause important transonic flow effects such as a reversal of control surface effectiveness that occurs as the shock wave crosses the hinge line. In order to account for flow complexities such as blade-vortex interactions of rotor blades higher-fidelity methods based on the Navier-Stokes equations are used. Reference 6 presents a procedure that uses the Navier-Stokes equations with moving-sheared grids and demonstrates up to 8 degrees of control-surface amplitude, using a single grid. Later, this procedure was extended to accommodate larger amplitudes, based on sliding grid zones. The sheared grid method implemented in EulerlNavier-Stokes-based aeroelastic code ENS AERO was successfully applied to active control design by industry. Recently there are several papers that present results for oscillating control surface using Reynolds Averaged Navier-Stokes (RANS) equations. References 9 and 10 report 2-D cases by filling gaps with overset grids. Reference 9 compares integrated forces with the experiment at low oscillating frequencies whereas Ref. 10 reports parametric studies but with no validation. Reference II reports results for a 3D case by modeling the gap region with a deformed grid and compares force results with the experiment only at the mid-span of flap. In Ref. II grid is deformed to match the control surface deflections at the section where the measurements are made. However, there is no indication in Ref. II that the gaps are explicitly modeled as in Ref. 6. Computations using overset grids are reported in Ref. 12 for a case by adding moving control surface to an existing blade but with no validation either with an experiment or another computation.

Navier↗
DOE OSTI2022

Probing topological phase transitions using high-harmonic generation

We report the prediction and realization of topological insulators have sparked great interest in experimental approaches to the classification of materials. The phase transition between non-trivial and trivial topological states is important, not only for basic materials science but also for next-generation technology, such as dissipation-free electronics. It is therefore crucial to develop advanced probes that are suitable for a wide range of samples and environments. Here we demonstrate that circularly polarized laser-field-driven high-harmonic generation is distinctly sensitive to the non-trivial and trivial topological phases in the prototypical three-dimensional topological insulator bismuth selenide. The phase transition is chemically initiated by reducing the spin–orbit interaction strength through the substitution of bismuth with indium atoms. We find strikingly different high-harmonic responses of trivial and non-trivial topological surface states that manifest themselves as a conversion efficiency and elliptical dichroism that depend both on the driving laser ellipticity and the crystal orientation. The origins of the anomalous high-harmonic response are corroborated by calculations using the semiconductor optical Bloch equations with pairs of surface and bulk bands. As a purely optical approach, this method offers sensitivity to the electronic structure of the material, including its nonlinear response, and is compatible with a wide range of samples and sample environments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Transient Lattice Response upon Photoexcitation in CuInSe 2 Nanocrystals with Organic or Inorganic Surface Passivation

CuInSe 2 nanocrystals offer promise for optoelectronics including thin-film photovoltaics and printed electronics. Additive manufacturing methods such as photonic curing controllably sinter particles into quasi-continuous films and offer improved device performance. To gain understanding of nanocrystal response under such processing conditions, we investigate impacts of photoexcitation on colloidal nanocrystal lattices via time-resolved X-ray diffraction. We probe three sizes of particles and two capping ligands (oleylamine and inorganic S 2- ) to evaluate resultant crystal lattice temperature, phase stability, and thermal dissipation. Elevated fluences produce heating and loss of crystallinity, the onset of which exhibits particle size dependence. We find size-dependent recrystallization and cooling lifetimes ranging from 90 to 200 ps with additional slower cooling on the nanosecond time scale. Sulfide-capped nanocrystals show faster recrystallization and cooling compared to oleylamine-capped nanocrystals. Using these lifetimes, we find interfacial thermal conductivities from 3 to 28 MW/(m 2 K), demonstrating that ligand identity strongly influences thermal dissipation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

Effect of Surface Nonequilibrium Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using finite-rate gas/surface interaction model provides time-accurate solutions for multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal Response and Ablation Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas momentum conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of gas/surface interaction chemistry between air and carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was a Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kang↗
NASA NTRS1988

A study of rain effects on radar scattering from water waves

Results are presented from a laboratory investigation of microwave power return due to rain-generated short waves on a wind wave surface. The wind wave tank, sensor, and data processing methods used in the study are described. The study focuses on the response of a 36-GHz radar system, orientated 30 deg from nadir and pointing upwind, to surface waves generated by various combinations of rain and wind. The results show stronger radar signal levels due to short surface waves generated by rain impacting the wind wave surface, supporting the results of Moore et al. (1979) for a 14-GHz radar.

Bliven, Larry F.↗
DOE OSTI2021

An Iterative Response-Surface-Based Approach for Chance-Constrained AC Optimal Power Flow Considering Dependent Uncertainty

A modern power system is characterized by a stochastic variation of the loads and an increasing penetration of renewable energy generation, which results in large uncertainties in its states. These uncertainties bring formidable challenges to the power system planning and operation process. To address these challenges, we propose a cost-effective, iterative response-surface-based approach for the chance-constrained AC optimal power-flow problem that aims to ensure the secure operation of the power systems considering dependent uncertainties. Starting from a stochastic-sampling-based framework, we first utilize the copula theory to simulate the dependence among multivariate uncertain inputs. Then, to reduce the prohibitive computational time required in the traditional Monte-Carlo method, we propose, instead of using the original complicated power-system model, to rely on a polynomial-chaos-based response surface. This response surface allows us to efficiently evaluate the time-consuming power-system model at arbitrary distributed sampled values with a negligible computational cost. This further enables us to efficiently conduct an online stochastic testing for the system states that not only screens out the statistical active constraints, but also assists in a better design of the tightened bounds without using any Gaussian or symmetric assumption. Finally, an iterative procedure is executed to fine-tune the optimal solution that better satisfies a predefined probability. The simulations conducted in multiple test systems demonstrate the excellent performance of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

The role of surfaces in flexoelectricity

We develop the relationship between the strain derivative of the mean-inner potential and surface contributions to flexoelectricity, identifying the true surface-specific component of the flexoelectric response of finite samples. Density functional theory calculations on a range of experimentally observed, low energy SrTiO 3 , MgO, and Si surfaces demonstrate that the mean-inner potential and its contributions to flexoelectricity are sensitive to small differences in surface structure, chemistry, and adsorbates. We also introduce a method to estimate mean-inner potential contributions to flexoelectricity using electron scattering factors and use this approximation to predict total flexoelectric responses for a variety of insulators. Strategies to experimentally disentangle bulk and surface flexoelectric terms are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2004

High-Order Shock-Capturing Methods for Modeling Dynamics of the Solar Atmosphere

We use one-dimensional high-order central shock capturing numerical methods to study the response of various model solar atmospheres to forcing at the solar surface. The dynamics of the atmosphere is modeled with the Euler equations in a variable-sized flux tube in the presence of gravity. We study dynamics of the atmosphere suggestive of spicule formation and coronal oscillations. These studies are performed on observationally-derived model atmospheres above the quiet sun and above sunspots. To perform these simulations, we provide a new extension of existing second- and third- order shock-capturing methods to irregular grids. We also solve the problem of numerically maintaining initial hydrostatic balance via the introduction of new variables in the model equations and a careful initialization mechanism. We find several striking results: all model atmospheres respond to a single impulsive perturbation with several strong shock waves consistent with the rebound-shock model. These shock waves lift material and the transition region well into the initial corona, and the sensitivity of this lift to the initial impulse depends non-linearly on the details of the atmosphere model. We also reproduce an observed 3-minute coronal oscillation above sunspots compared to 5-minute oscillations above the quiet sun.

Bryson, Steve↗
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↗