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 127 records · Page 7

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.↗

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↗

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↗

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↗

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.↗

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.↗

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↗

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.↗

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↗

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↗

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.↗

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↗

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↗

Mesh Deformation Boundary Conditions for Three-Dimensional Ablation Solvers

A method for determining boundary node displacements for three-dimensional ablation solvers is developed. The approach is applicable to a wide range of mesh motion algorithms used to determine internal node displacements. The method leverages radial basis functions to analytically define the geometry to which the boundary nodes should be constrained. Constrained optimization, elasticity, and spring analogy solvers are utilized to redistribute the nodes on the evolved geometry surface while maintaining mesh quality. The method is implemented in the CHarring Ablator Response (CHAR) code, and a realistic demonstration case is presented on the Boundary Layer Transition (BOLT) flight experiment configuration. The method is shown to be robust for highly-curved geometries with large domain deformations.

Adam J Amar↗

On Calculation Methods and Results for Straight Cylindrical Roller Bearing Deflection, Stiffness, and Stress

The purpose of this study was to assess some calculation methods for quantifying the relationships of bearing geometry, material properties, load, deflection, stiffness, and stress. The scope of the work was limited to two-dimensional modeling of straight cylindrical roller bearings. Preparations for studies of dynamic response of bearings with damaged surfaces motivated this work. Studies were selected to exercise and build confidence in the numerical tools. Three calculation methods were used in this work. Two of the methods were numerical solutions of the Hertz contact approach. The third method used was a combined finite element surface integral method. Example calculations were done for a single roller loaded between an inner and outer raceway for code verification. Next, a bearing with 13 rollers and all-steel construction was used as an example to do additional code verification, including an assessment of the leading order of accuracy of the finite element and surface integral method. Results from that study show that the method is at least first-order accurate. Those results also show that the contact grid refinement has a more significant influence on precision as compared to the finite element grid refinement. To explore the influence of material properties, the 13-roller bearing was modeled as made from Nitinol 60, a material with very different properties from steel and showing some potential for bearing applications. The codes were exercised to compare contact areas and stress levels for steel and Nitinol 60 bearings operating at equivalent power density. As a step toward modeling the dynamic response of bearings having surface damage, static analyses were completed to simulate a bearing with a spall or similar damage.

Krantz, Timothy L.↗

An Approximate Ablative Thermal Protection System Sizing Tool for Entry System Design

A computer tool to perform entry vehicle ablative thermal protection systems sizing has been developed. Two options for calculating the thermal response are incorporated into the tool. One, an industry-standard, high-fidelity ablation and thermal response program was integrated into the tool, making use of simulated trajectory data to calculate its boundary conditions at the ablating surface. Second, an approximate method that uses heat of ablation data to estimate heat shield recession during entry has been coupled to a one-dimensional finite-difference calculation that calculates the in-depth thermal response. The in-depth solution accounts for material decomposition, but does not account for pyrolysis gas energy absorption through the material. Engineering correlations are used to estimate stagnation point convective and radiative heating as a function of time. The sizing tool calculates recovery enthalpy, wall enthalpy, surface pressure, and heat transfer coefficient. Verification of this tool is performed by comparison to past thermal protection system sizings for the Mars Pathfinder and Stardust entry systems and calculations are performed for an Apollo capsule entering the atmosphere at lunar and Mars return speeds.

Dec, John A.↗

An Approximate Ablative Thermal Protection System Sizing Tool for Entry System Design

A computer tool to perform entry vehicle ablative thermal protection systems sizing has been developed. Two options for calculating the thermal response are incorporated into the tool. One, an industry-standard, high-fidelity ablation and thermal response program was integrated into the tool, making use of simulated trajectory data to calculate its boundary conditions at the ablating surface. Second, an approximate method that uses heat of ablation data to estimate heat shield recession during entry has been coupled to a one-dimensional finite-difference calculation that calculates the in-depth thermal response. The in-depth solution accounts for material decomposition, but does not account for pyrolysis gas energy absorption through the material. Engineering correlations are used to estimate stagnation point convective and radiative heating as a function of time. The sizing tool calculates recovery enthalpy, wall enthalpy, surface pressure, and heat transfer coefficient. Verification of this tool is performed by comparison to past thermal protection system sizings for the Mars Pathfinder and Stardust entry systems and calculations are performed for an Apollo capsule entering the atmosphere at lunar and Mars return speeds.

Dec, John A.↗