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 181 records · Page 10

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

Use of Subsonic Kernel Function in an Influence-Coefficient Method of Aeroelastic Analysis and some Comparisons with Experiment

This paper illustrates the development and application of an influence-coefficient method of analysis for calculating the response of a flexible wing in an airstream to an oscillating disturbing force and for treating such aeroelastic instabilities as flutter and divergence. Aerodynamic coefficients are derived on the basis of lifting - surface theory for subsonic compressible flow by use of the method presented in NASA Technical Report R-48. Application of the analysis is made to a uniform cantilever wing- tip tank configuration for which responses to a sinusoidal disturbing force and flutter speeds were measured over a range of subsonic Mach numbers and densities. Calculated responses and flutter speeds based on flexibility influence coefficients measured at nine stations are in good agreement with experiment, provided the aerodynamic load is distributed over the wing so that local centers of pressure very nearly coincide with these nine influence stations. The use of experimental values of bending and torsional structural damping coefficients in the analysis generally improved the agreement between calculated and experimental responses. Some calculations were made to study the effects on density on responses near the flutter conditions, and linear response trends were obtained over a wide range of densities.

Sewall, John L.↗

Effect of Non-Equilibrium Surface Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using non-equilibrium gas/surface interaction model provides time-accurate solutions for the 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 mass 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 the 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 nonequilibrium gas/surface interaction chemistry between air and the 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 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-Kanq↗

Micropolar Elastoplasticity Using a Fast Fourier Transform‐Based Solver

ABSTRACT This work presents a micromechanical spectral formulation for obtaining the full‐field and homogenized response of elastoplastic micropolar composites. A closed‐form radial‐return mapping is derived from thermodynamics‐based micropolar elastoplastic constitutive equations to determine the increment of plastic strain necessary to return the generalized stress state to the yield surface, and the algorithm implementation is verified using the method of numerically manufactured solutions. Then, size‐dependent material response and micro‐plasticity are shown as features that may be efficiently simulated in this micropolar elastoplastic framework. The computational efficiency of the formulation enables the generation of large datasets in reasonable computing times.

42 ENGINEERING↗

Non-Contact Thermal Properties Measurement with Low-Power Laser and IR Camera System

As shown by the Phoenix Mars Lander's Thermal and Electrical Conductivity Probe (TECP), contact measurements of thermal conductivity and diffusivity (using a modified flux-plate or line-source heat-pulse method) are constrained by a number of factors. Robotic resources must be used to place the probe, making them unavailable for other operations for the duration of the measurement. The range of placement is also limited by mobility, particularly in the case of a lander. Placement is also subject to irregularities in contact quality, resulting in non-repeatable heat transfer to the material under test. Most important from a scientific perspective, the varieties of materials which can be measured are limited to unconsolidated or weakly-cohesive regolith materials, rocks, and ices being too hard for nominal insertion strengths. Accurately measuring thermal properties in the laboratory requires significant experimental finesse, involving sample preparation, controlled and repeatable procedures, and, practically, instrumentation much more voluminous than the sample being tested (heater plates, insulation, temperature sensors). Remote measurements (infrared images from orbiting spacecraft) can reveal composite properties like thermal inertia, but suffer both from a large footprint (low spatial resolution) and convolution of the thermal properties of a potentially layered medium. In situ measurement techniques (the Phoenix TECP is the only robotic measurement of thermal properties to date) suffer from problems of placement range, placement quality, occupation of robotic resources, and the ability to only measure materials of low mechanical strength. A spacecraft needs the ability to perform a non-contact thermal properties measurement in situ. Essential components include low power consumption, leveraging of existing or highly-developed flight technologies, and mechanical simplicity. This new in situ method, by virtue of its being non-contact, bypasses all of these problems. The use of photons to both excite and measure the thermal response of any surface material to a high resolution (estimated footprint = 10 square centimeters) is a generational leap in physical properties measurements. The proposed method consists of spot-heating the surface of a material with a low (less than 1 W) power laser. This produces a moderate (5-10 K) temperature increase in the material.

Hudson, Troy L.↗

A method of computing the transient temperature of thick walls from arbitrary variation of adiabatic-wall temperature and heat-transfer coefficient

A method of calculating the temperature of thick walls has been developed in which the time series and the response to a unit triangle variation of surface temperature concepts are used, together with essentially standard formulas for transient temperature and heat flow into thick walls. The method can be used without knowledge of the mathematical tools of its development. The method is particularly suitable for determining the wall temperature in one-dimensional thermal problems in aeronautics where there is a continuous variation of the heat-transfer coefficient and adiabatic-wall temperature. The method also offers a convenient means for solving the inverse problem of determining the heat-flow history when temperature history is known.

Hill, P R↗

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine a reduced-order model (ROM) of a computational fluid dynamics (CFD) solver in support of linear aeroservoelastic model development and feedback control design. The approach was applied to the FUN3D code for the half-span wind tunnel test article used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. In a transonic flow condition, multiple inputs (11 structural mode displacements and 3 control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines while multiple outputs (the corresponding 14 generalized aerodynamic forces) were recorded. From these recorded times series, the matrix of frequency responses was computed and subsequently fit using rational function approximations (RFAs). It was found that the entire (14 x 14) matrix of frequency responses could be determined from a single CFD run and that results generally followed trends predicted using other methods. Differences were attributed to the modeling fidelity and nonlinearities from structural mode and control surface interactions at higher reduced frequencies. More accurate fits of the RFAs to the frequency response data were obtained by making two CFD runs, one with only structural mode excitations and one with only control surface excitations, which reduced the degree of nonlinearity in the modeling data.

Aeroservoelasticity↗

MEDLI2: MISP Inferred Aerothermal Environment and Flow Transition Assessment

The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite on the Mars2020 mission contained multiple sensors on the aeroshell to measure the aerothermal environment during entry into the Martian atmosphere. These sensors performed superbly and successfully returned forebody and aftbody heating measurements. Analysis of MEDLI2 data indicated flow transitioning from a laminar to turbulent state on the heatshield. No evidence of flow transition was observed on the backshell. Two methods were used to estimate flow transition times on the heatshield: (1) temperature gradient of near-surface thermocouple data and (2) heat flux gradient from an inverse reconstruction approach using thermocouple data and material response modeling. Both methods produced similar transition times with an estimated accuracy of ±1 s. To assess various transition criteria, transition parameters were evaluated at each sensor location using flow field solutions from computational fluid dynamics (CFD) simulations. The idea was to use conservative values inferred from MEDLI2 data as transition criteria for other Mars missions. To test this hypothesis, MEDLI data from the Mars Science Laboratory (MSL) mission was used to compare predicted vs. actual flow transition times. The comparisons suggest smooth wall transition criteria are not well-suited in modeling the rapid progression of a turbulent transition front. Transition criteria containing a roughness element parameter agreed better with the flight data. In summary, critical transition values derived from MEDLI2 data may be used as a starting point in constructing a flow transition model for future Mars missions.

Chun Y Tang↗

Interpretation of scanning electron microscope measurements of minority carrier diffusion lengths in semiconductors

In scanning electron microscope (SEM) injection measurements of minority carrier diffusion lengths some uncertainties of interpretation exist when the response current is nonlinear with distance. This is significant in epitaxial layers where the layer thickness is not large in relation to the diffusion length, and where there are large surface recombination velocities on the incident and contact surfaces. An image method of analysis is presented for such specimens. A method of using the results to correct the observed response in a simple convenient way is presented. The technique is illustrated with reference to measurements in epitaxial layers of GaAs. Average beam penetration depth may also be estimated from the curve shape.

Flat, A.↗

A Reduced-frequency Approach for Calculating Dynamic Derivatives

Computational Fluid Dynamics (CFD) is increasingly being used to both augment and create an aerodynamic performance database for aircraft configurations. This aerodynamic database contains the response of the aircraft to varying flight conditions and control surface deflections. The current work presents a novel method for calculating dynamic stability derivatives which reduces the computational cost over traditional unsteady CFD approaches by an order of magnitude, while still being applicable to arbitrarily complex geometries over a wide range of flow regimes. The primary thesis of this work is that the response to a forced motion can often be represented with a small, predictable number of frequency components without loss of accuracy. By resolving only those frequencies of interest, the computational effort is significantly reduced so that the routine calculation of dynamic derivatives becomes practical. The current implementation uses this same non-linear, frequency-domain approach and extends the application to the 3-D Euler equations. The current work uses a Cartesian, embedded-boundary method to automate the generation of dynamic stability derivatives.

Murman, Scott M.↗

Stereology techniques in radiation biology

Clinicians involved in conventional radiation therapy are very concerned about the dose-response relationships of normal tissues. Before proceeding to new clinical protocols, radiation biologists involved with conformal proton therapy believe it is necessary to quantify the dose response and tolerance of the organs and tissues that will be irradiated. An important focus is on the vasculature. This presentation reviews the methodology and format of using confocal microscopy and stereological methods to quantify tissue parameters, cell number, tissue volume and surface area, and vessel length using the microvasculature as a model tissue. Stereological methods and their concepts are illustrated using an ongoing study of the dose response of the microvessels in proton-irradiated hemibrain. Methods for estimating the volume of the brain and the brain cortex, the total number of endothelial cells in cortical microvessels, the length of cortical microvessels, and the total surface area of cortical microvessel walls are presented step by step in a way understandable for readers with little mathematical background. It is shown that stereological techniques, based on a sound theoretical basis, are powerful and reliable and have been used successfully.

Non-NASA Center↗

Closed Loop Active Flow Separation Detection and Control in a Multistage Compressor

Active closed loop flow control was successfully demonstrated on a full annulus of stator vanes in a low speed axial compressor. Two independent methods of detecting separated flow conditions on the vane suction surface were developed. The first technique detects changes in static pressure along the vane suction surface, while the second method monitors variation in the potential field of the downstream rotor. Both methods may feasibly be used in future engines employing embedded flow control technology. In response to the detection of separated conditions, injection along the suction surface of each vane was used. Injected mass flow on the suction surface of stator vanes is known to reduce separation and the resulting limitation on static pressure rise due to lowered diffusion in the vane passage. A control algorithm was developed which provided a proportional response of the injected mass flow to the degree of separation, thereby minimizing the performance penalty on the compressor system.

Bright, Michelle M.↗

OSEE responses on D6AC steel due to sample preparation

A cleaning procedure and a contamination method developed at NASA Langley to contaminate sample surfaces with a certain thickness of HD-2 grease are demonstrated. Varying levels of contamination are obtained. Optically stimulated electron emission (OSEE) measurements were performed using the parallel electric field cell and the Pat scan device for clean and contaminated samples with a known amount of common contaminant. The OSEE response is found to be more sensitive for the clean sample than for the contaminated one. It is concluded that the OSEE measurements obtained using these techniques are sensitive to surface conditions and the contaminant thickness variations. It is also shown that the wavelength of illuminating light can have a profound effect on OSEE current levels.

Abedin, M. N.↗

Turbine Engine Hot Section Technology, 1987

Presentations were made concerning the development of design analysis tools for combustor liners, turbine vanes, and turbine blades. Presentations were divided into six sections: instrumentation, combustion, turbine heat transfer, structural analysis, fatigue and fracture, surface protective coatings, constitutive behavior of materials, stress-strain response and life prediction methods.

Source record↗