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 541 records · Page 30

Aerodynamic configuration design using response surface methodology analysis

An investigation has been conducted to determine a set of optimal design parameters for a single-stage-to-orbit reentry vehicle. Several configuration geometry parameters which had a large impact on the entry vehicle flying characteristics were selected as design variables: the fuselage fineness ratio, the nose to body length ratio, the nose camber value, the wing planform area scale factor, and the wing location. The optimal geometry parameter values were chosen using a response surface methodology (RSM) technique which allowed for a minimum dry weight configuration design that met a set of aerodynamic performance constraints on the landing speed, and on the subsonic, supersonic, and hypersonic trim and stability levels. The RSM technique utilized, specifically the central composite design method, is presented, along with the general vehicle conceptual design process. Results are presented for an optimized configuration along with several design trade cases.

Engelund, Walter C.↗

Electromagnetic Sounding of the Moon Using Apollo 16 and Lunokhod 2 Surface Magnetometer Observations (Preliminary Results)

A new technique of deep electromagnetic sounding of the Moon using simultaneous magnetic field measurements at two lunar surface sites is described. The method, used with the assumption that deep electrical conductivity is a function only of lunar radius, has the advantage of allowing calculation of the external driving field from two surface site measurements only, and therefore does not require data from a lunar orbiting satellite. A transient response calculation is presented for the example of a magnetic field discontinuity of February 13, 1973, measured simultaneously by Apollo 16 and Lunokhod 2 surface magnetometers.

Vanyan, L. L.↗

Modulating the Pseudoelastic Response of Nitinol Using Ion Implantation

This work explores whether ion beam modification can be used to modulate the austenite to martensite phase transformation in Nickel-Titanium (NiTi), thereby achieving novel or localized transformation properties in near-surface regions. We report this could provide alternatives to laser shot peening or other surface treatment methods and possibly expand applications in biomedical, aerospace, and other fields. Irradiation induces defects and internal stress that can serve as nucleation and/or pinning sites for the phase transformation. Thus, it can augment more convention- al approaches, including alloying, severe mechanical work, grain size reduction, and precipitation of coherent precipitates. A range of outcomes is possible in principle, including a shift of the critical stress or temperature for onset of the transformation, linearization, reduction of hysteresis, stabilization, and extent of transformation strain.

36 MATERIALS SCIENCE↗

Morphology, microstructure, and doping behaviour: A comparison between different deposition methods for poly-Si/SiO x passivating contacts

In this work, we study how crystallographic structures, optoelectronic properties, and nanoscale surface morphologies of ex situ phosphorus-doped polycrystalline silicon (poly-Si)/SiO x passivating contacts, formed by different deposition methods (sputtering, plasma-enhanced chemical vapour deposition [PECVD], and low-pressure chemical vapour deposition [LPCVD]), are investigated and compared. Across all these deposition technologies, we noted the same trend: higher diffusion temperatures yield films that are more crystalline but that have rougher surface morphologies due to bigger surface crystal grains. Also, the recrystallization process of the as-deposited Si films starts from the SiO x interface, rather than from the film surface and bulk. However, there are some distinct differences among these technologies. First, the LPCVD method yields the lowest deposition rate, roughest surfaces, and smallest degree of crystallinity on finished poly-Si films. In contrast, the PECVD method has the highest deposition rate and smoothest surfaces for both as-deposited Si and annealed poly-Si films. Second, as-deposited sputtered and PECVD Si films contain only an amorphous phase, whereas as-deposited LPCVD films already has some crystalline phase. Third, the LPCVD phosphorus in-diffusion into the substrate depends strongly on the initial film thickness, whereas for the other two methods, it is weakly dependent on thickness. Finally, the passivation quality of every poly-Si film type has different responses to the film thickness and diffusion temperature, suggesting that the ex situ doping optimization should be performed independently.

14 SOLAR ENERGY↗

Roll Damping Derivatives from Generalized Lifting-Surface Theory and Wind Tunnel Forced-Oscillation Tests

Improving aerodynamic models for adverse loss-of-control conditions in flight is an area being researched under the NASA Aviation Safety Program. Aerodynamic models appropriate for loss of control conditions require a more general mathematical representation to predict nonlinear unsteady behaviors. As more general aerodynamic models are studied that include nonlinear higher order effects, the possibility of measurements that confound aerodynamic and structural responses are probable. In this study an initial step is taken to look at including structural flexibility in analysis of rigid-body forced-oscillation testing that accounts for dynamic rig, sting and balance flexibility. Because of the significant testing required and associated costs in a general study, it makes sense to capitalize on low cost analytical methods where possible, especially where structural flexibility can be accounted for by a low cost method. This paper provides an initial look at using linear lifting surface theory applied to rigid-body aircraft roll forced-oscillation tests.

Pototzky, Anthony S↗

Alternative Method of On-Orbit Response-Versus-Scan-Angle Characterization for MODIS Reflective Solar Bands

The moderate resolution imaging spectroradiometer (MODIS) has 20 reflective solar bands (RSB), covering a spectral range from 0.41 to 2.2 microns, which are calibrated on-orbit using its onboard calibrators, which include a solar diffuser, a solar diffuser stability monitor, and a spectroradiometric calibration assembly. A space view (SV) port is used to provide a background reference and also facilitates near-monthly lunar observations through a spacecraft roll. In every scan, the Earth's surface, SV, and onboard calibrators are viewed via a two-sided scan mirror, the reflectance of which depends on the angle of incidence (AOI) as well as the wavelength of the incident light. Response-versus-scan-angle (RVS) is defined as a dependence function of the scan mirror's reflectance over AOI. An initial RVS for each RSB was measured prelaunch for both Terra and Aqua MODIS. Algorithms have been developed to track the on-orbit RVS variation using the measurements from the onboard calibrators, supplemented with the earth view (EV) trends from pseudoinvariant desert targets obtained at different AOI. Since the mission beginning, the MODIS characterization support team (MCST) has dedicated efforts in evaluating approaches of characterizing the on-orbit RVS. A majority of the approaches focused on fitting the data at each AOI over time and then deriving the relative change at different AOI. The current version of the on-orbit RVS algorithm, as implemented in the collection 6 (C6) level-1B (L1B), is also based on the above rationale. It utilizes the EV response trends from the pseudoinvariant Libyan desert targets to supplement the gain derived from the onboard calibrators. The primary limitation of this approach is the assumption of the temporal stability of these desert sites. Consequently, MCST developed an approach that derives the on-orbit RVS change using measurements from a single desert site, combined with the on-orbit lunar measurements. In addition, the EV and onboard responses are fit first as a function of AOI before fitting temporally in order to eliminate the dependence on the stability of the desert site. Comprehensive comparisons are performed with current C6 RVS results for both Terra and Aqua MODIS. Results demonstrate that this alternative method provides a supplemental means to monitor the on-orbit RVS for MODIS RSB.

solar bands↗

Process design tools and techno-economic analysis for capacitive deionization

Capacitive deionization (CDI) devices use cyclical electrosorption on porous electrode surfaces to achieve water desalination. Process modeling and design of CDI systems requires accurate treatment of the coupling among input electrical forcing, input flow rates, and system responses including salt removal dynamics, water recovery, energy storage, and dissipation. Techno-economic analyses of CDI further require a method to calculate and compare between a produced commodity (e.g. desalted water) versus capital and operational costs of the system. In this work, we demonstrate a new modeling and analysis tool for CDI developed as an installable Matlab program that allows direct numerical simulation of CDI dynamics and calculation of key performance and cost parameters. The program is provided for free and is used to run open-source Simulink models. The Simulink environment sends information to the program and allows for a drag and drop design space where users can connect CDI cells to relevant periphery blocks such as grid energy, battery, solar panel, waste disposal, and maintenance/labor cost streams. The program allows for simulation of arbitrary current forcing and arbitrary flow rate forcing of one or more CDI cells. We employ validated well-mixed reactor formulations together with a non-linear circuit model formulation that can accommodate a variety of electric double layer sub-models (e.g. for charge efficiency). The program includes a graphical user interface (GUI) to specify CDI plant parameters, specify operating conditions, run individual tests or parameter batch-mode simulations, and plot relevant results. The techno-economic models convert among dimensional streams of species (e.g. feed, desalted water, and brine), energy, and cost and enable a variety of economic estimates including levelized water costs.

42 ENGINEERING↗

Statistically elucidated responses from low-signal contrast mechanisms in ultrafast electron microscopy

The emergence of ultrafast electron microscopy (UEM) has enabled the discovery of strongly correlated dynamic mechanisms, including electron–phonon coupling, structural phase transitions, thermal transport, and electromagnetic deflection. Most UEM systems operate stroboscopically, meaning that the technique is susceptible to artifacts, mistakes, and misinterpretation of the data due to extensive experimental effort. In contrast to the ultrafast designation, data acquisition is extraordinarily slow because the electron beam has significantly reduced signal compared to traditional transmission electron microscopy due to pulsing the electron beam. Consequently, the sample may drift, tilt, or undergo irreversible structural changes that are independent of the time-resolved dynamics throughout the experimental time frame. Furthermore, these datasets require significant user interpretation that can be problematic when proper controls are not implemented thoroughly. Here, we demonstrate a new algorithm designed to separate ultrafast structural dynamics from long-term artifacts using a LiNbO 3 sample experiencing electrically driven surface acoustic wave propagation. Additionally, we provide examples of the impact of user bias when analyzing the data and provide a methodology, which enables the extraction of time-resolved responses when the image signal is extraordinarily low. Overall, the goal of this publication is to provide methods that validate the experimental results and reduce researcher biases during UEM data interpretation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Separation of lymphocytes by electrophoresis under terrestrial conditions and at zero gravity

Electrophoretic mobility (EPM) of human peripheral lymphocytes were examined with the following objectives: To determine differences in EPM of lymphocytes under immuno-stimulated and immuno-suppressed states. To define the conditions necessary for the separation of lymphocyte sub-populations in normal and pathological conditions; To investigate immunological active, charged chemical groups on lymphocyte surfaces; and to investigate pathophysiological mechanisms of immune responsiveness, as reflected by alterations in EPM. To evaluate the potential of lymphocyte electrophoresis as: (1) a means of monitoring the immune status of kidney transplant recipients, (2) in predicting the outcome of kidney transplants, and (3) as a method for separation of lymphocyte sub-populations, the EPM was studied for unfractionated human peripheral lymphocytes and of populations enriched with T and "B" cells from normal adults, hemodialysis patients and kidney transplant recipients.

Rubin, A. L.↗

Shape Measurement of Large Aerospace Structures Using Sensitivity Electrical TDR Distributed Strain Sensor

Electrical time domain reflectometry (ETDR) sensing technique can be best described as "closed-loop radar," where the information is derived from the reflections of a voltage pulse sent through a transmission medium. The ETDR sensing technique is a well-developed method and has been widely used to locate and evaluate discontinuities in long coaxial power transmission cables. The ETDR technique provides a true distributed sensing capability which can not only sense the distributed loading condition of the structure but also can pin-point the location of disturbance, such as the locations of stress concentration and structural damages. Proof-of-concept experiments have been conducted using photoelastic specimens with embedded commercial coaxial cables, i.e., RG85/U and RG174, to demonstrate the stress/strain sensing capability of ETDR sensors for structural health monitoring application. Although the test results showed that the ETDR sensor signals capture specimen deformation pattern both in bending and tension and indicate the location and type of crack damages of the photoelastic specimen; yet, the low signal-to-noise ratio of the sensor signal smears the details of the strain measurement that the ETDR signals can convey. A high-sensitivity ETDR coaxial strain sensor prototype newly developed at Clark Atlanta University will be presented. The construction of the prototype sensing cable as well as its electrical properties relevant to distributed strain sensing application will be shown in details. Test results of the sensitivity and tension responses of the ETDR signal of the prototype sensor will be presented and compared with those of commercial coaxial cables. Promising potentials of the ETDR distributed strain sensing method for shape measurement application of large aerospace structures will also be demonstrated using long slender beam with surface-bonded ETDR distributed strain sensor.

Lin, Mark W.↗

Boundary layer receptivity and control

Receptivity processes initiate natural instabilities in a boundary layer. The instabilities grow and eventually break down to turbulence. Consequently, receptivity questions are a critical element of the analysis of the transition process. Success in modeling the physics of receptivity processes thus has a direct bearing on technological issues of drag reduction. The means by which transitional flows can be controlled is also a major concern: questions of control are tied inevitably to those of receptivity. Adjoint systems provide a highly effective mathematical method for approaching many of the questions associated with both receptivity and control. The long term objective is to develop adjoint methods to handle increasingly complex receptivity questions, and to find systematic procedures for deducing effective control strategies. The most elementary receptivity problem is that in which a parallel boundary layer is forced by time-harmonic sources of various types. The characteristics of the response to such forcing form the building blocks for more complex receptivity mechanisms. The first objective of this year's research effort was to investigate how a parallel Blasius boundary layer responds to general direct forcing. Acoustic disturbances in the freestream can be scattered by flow non-uniformities to produce Tollmien-Schlichting waves. For example, scattering by surface roughness is known to provide an efficient receptivity path. The present effort is directed towards finding a solution by a simple adjoint analysis, because adjoint methods can be extended to more complex problems. In practice, flows are non-parallel and often three-dimensional. Compressibility may also be significant in some cases. Recent developments in the use of Parabolized Stability Equations (PSE) offer a promising possibility. By formulating and solving a set of adjoint parabolized equations, a method for mapping the efficiency with which external forcing excites the three-dimensional motions of a non-parallel boundary layer was developed. The method makes use of the same computationally efficient formulation that makes the PSE currently so appealing. In the area of flow control, adjoint systems offer a powerful insight into the effect of control forces. One of the simplest control strategies for boundary layers involves the application of localized mean wall suction.

Hill, D. C.↗

Comparison of theory with the experimental characterization of the spatial frequency response of interferometers using a binary pseudo-random array sample

Experimental evaluations of the surface height response of an interference microscope using a binary pseudo-random array test sample are compared with a theory based on a Fourier optics model. Measurements of key instrument characteristics, including the illumination, imaging, and obscuring apertures of three different Mirau objectives, support the theoretical calculations. Agreement between experimental and theoretical modeling confirms the predictability of the spatial frequency response for the purpose of specification and optimization of instrument configuration for specific metrology tasks. The results also provide confidence in methods of compensating for the decrease in instrument response with spatial frequency.

Calibration↗

Post-Buckling Analysis of Curved Honeycomb Sandwich Panels Containing Interfacial Disbonds

A numerical study on the effect of facesheet-core disbonds on the post-buckling response of curved honeycomb sandwich panels is presented herein. This work was conducted as part of the development of a damage tolerance plan for the next-generation Space Launch System heavy lift launch vehicle payload fairing. As such, the study utilized full-scale fairing barrel segments as the structure of interest. The panels were composed of carbon fiber reinforced polymer facesheets and aluminum honeycomb core. The panels were analyzed numerically using the finite element method incorporating geometric nonlinearity. In a predetermined circular region, facesheet and core nodes were detached to simulate a disbond, between the outer mold line facesheet and honeycomb core, induced via low-speed impact. Surface-to-surface contact in the disbonded region was invoked to prevent interpenetration of the facesheet and core elements and obtain realistic stresses in the core. The diameter of this disbonded region was varied and the effect of the size of the disbond on the post-buckling response was observed. Significant changes in the slope of the edge load-deflection response were used to determine the onset of global buckling and corresponding buckling load. Finally, several studies were conducted to determine the sensitivity of the numerical predictions to refinement in the finite element mesh.

elastic buckling↗

Shadowgraph Study of Gradient Driven Fluctuations

A fluid or fluid mixture, subjected to a vertical temperature and/or concentration gradient in a gravitational field, exhibits greatly enhanced light scattering at small angles. This effect is caused by coupling between the vertical velocity fluctuations due to thermal energy and the vertically varying refractive index. Physically, small upward or downward moving regions will be displaced into fluid having a refractive index different from that of the moving region, thus giving rise to the enhanced scattering. The scattered intensity is predicted to vary with scattering wave vector q, as q(sup -4), for sufficiently large q, but the divergence is quenched by gravity at small q. In the absence of gravity, the long wavelength fluctuations responsible for the enhanced scattering are predicted to grow until limited by the sample dimensions. It is thus of interest to measure the mean-squared amplitude of such fluctuations in the microgravity environment for comparison with existing theory and ground based measurements. The relevant wave vectors are extremely small, making traditional low-angle light scattering difficult or impossible because of stray elastically scattered light generated by optical surfaces. An alternative technique is offered by the shadowgraph method, which is normally used to visualize fluid flows, but which can also serve as a quantitative tool to measure fluctuations. A somewhat novel shadowgraph apparatus and the necessary data analysis methods will be described. The apparatus uses a spatially coherent, but temporally incoherent, light source consisting of a super-luminescent diode coupled to a single-mode optical fiber in order to achieve extremely high spatial resolution, while avoiding effects caused by interference of light reflected from the various optical surfaces that are present when using laser sources. Results obtained for a critical mixture of aniline and cyclohexane subjected to a vertical temperature gradient will be presented. The sample was confined between two horizontal parallel sapphire plates with a vertical spacing of 1 mm. The temperatures of the sapphire plates were controlled by independent circulating water loops that used Peltier devices to add or remove heat from the room air as required. For a mixture with a temperature gradient, two effects are involved in generating the vertical refractive index gradient, namely thermal expansion and the Soret effect, which generates a concentration gradient in response to the applied temperature gradient. For the aniline/cyclohexane system, the denser component (aniline) migrates toward the colder surface. Consequently, when heating from above, both effects result in the sample density decreasing with altitude and are stabilizing in the sense that no convective motion occurs regardless of the magnitude of the applied temperature gradient. The Soret effect is strong near a binary liquid critical point, and thus the dominant effect is due to the induced concentration gradient. The results clearly show the divergence at low q and the predicted gravitational quenching. Results obtained for different applied temperature gradients at varying temperature differences from the critical temperature, clearly demonstrate the predicted divergence of the thermal diffusion ratio. Thus, the more closely the critical point is approached, the smaller becomes the temperature gradient required to generate the same signal. Two different methods have been used to generate pure concentration gradients. In the first, a sample cell was filled with a single fluid, ethylene glycol, and a denser miscible fluid, water, was added from below thus establishing a sharp interface to begin the experiment. As time went on the two fluids diffused into each other, and large amplitude fluctuations were clearly observed at low q. The effects of gravitational quenching were also evident. In the second method, the aniline/cyclohexane sample was used, and after applying a vertical temperature gradient for several hours, the top and bottom temperatures were set equal and the thermal gradient died on a time scale of seconds, leaving the Soret induced concentration gradient in place. Again, large-scale fluctuations were observed and died away slowly in amplitude as diffusion destroyed the initial concentration gradient.

Cannell, David↗

Lunar soil properties and soil mechanics. Flow in porous media under rarefied gas conditions. Research phase: Fluid conductivity of lunar surface materials

Theoretical and experimental research on fluid conductivity of lunar surface materials is summarized. Theoretical methods were developed for the analysis of transitional and free-molecular flows, and for analysis of lunar permeability probe data in general. Experimental studies of rarefied flows under conditions of a large pressure gradient show flows in the continuum regime to be responsible for the largest portion of the pressure drop between source and sink for one dimensional flow, provided the entrance Knudsen number is sufficiently small. The concept of local similarity leading to a universal nondimensional function of Knudsen number was shown to have approximate validity; flows in all regimes may be described in terms of an area fraction and a single length parameter. Synthetic porous media prepared from glass beads exhibited flow behavior similar in many regards to that of a natural sandstone; studies using artificial stones with known pore configurations may lead to new insight concerning the structure of natural materials. The experimental method involving the use of segmented specimens of large permeability is shown to be fruitful.

Hurlbut, F. C.↗

Electromagnetic sounding of the moon using Apollo 16 and Lunokhod 2 surface magnetometer observations /preliminary results/

A technique of deep electromagnetic sounding of the moon using simultaneous magnetic-field measurements at two lunar surface sites is described. The method, used with the assumption that deep electrical conductivity is a function only of lunar radius, has the advantage of allowing calculation of the external driving field from two surface-site measurements only and therefore does not require data from a lunar orbiting satellite. A transient-response calculation is presented for the example of a magnetic-field discontinuity, measured simultaneously by Apollo 16 and Lunokhod 2 surface magnetometers.

Vanian, L. L.↗

Investigations of plasma response associated with resonant magnetic perturbation fields using perturbation method in KSTAR H-mode plasmas

Abstract The plasma response associated with the resonant magnetic perturbation (RMP) field was investigated using the small edge perturbations induced by a modulated supersonic molecular beam injection (SMBI) in KSTAR. The modulated SMBI provides a time-varying perturbation of the plasma density source in the region just inside the last closed flux surface and a modulated flow damping rate. Radial propagation of the toroidal rotation perturbation induced by SMBI from the q = 3 surface to the q = 2 surface was observed. Theoretical analysis using the general perturbed equilibrium code of the RMP intensity profiles of the RMP field is consistent with the phase profile of the toroidal rotation perturbation.

Physics↗

Tangle-Free Mesh Motion for Ablation Simulations

Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and correctness. The focus of the presentation will be on the algorithm; specifics on how the techniques may be used in spacecraft design will be not discussed.

Droba, Justin↗