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 91 records · Page 5

Wireless Measurement of Contact and Motion Between Contact Surfaces

This method uses a magnetic-field- response contact sensor that is designed to identify surface contact and motion between contact locations. The sensor has three components: (1) a capacitor-inductor circuit with two sets of electrical contact pads, (2) a capacitor with a set of electrical contact pads, and (3) an inductor with a set of electrical contact pads. A unique feature of this sensor is that it is inherently multifunctional. Information can be derived from analyzing such sensor response attributes as amplitude, frequency, and bandwidth. A change in one attribute can be due to a change in a physical property of a system. A change in another attribute can be due to another physical property, which has no relationship to the first one.

Woodard, Stanley E.↗

Improved Sizing of Impact Damage in Composites Based on Thermographic Response

Impact damage in thin carbon fiber reinforced polymer composites often results in a relatively small region of damage at the front surface, with increasing damage near the back surface. Conventional methods for reducing the pulsed thermographic responses of the composite tend to underestimate the size of the back surface damage, since the smaller near surface damage gives the largest thermographic indication. A method is presented for reducing the thermographic data to produce an estimated size for the impact damage that is much closer to the size of the damage estimated from other NDE techniques such as microfocus x-ray computed tomography and pulse echo ultrasonics. Examples of the application of the technique to experimental data acquired on specimens with impact damage are presented. The method is also applied to the results of thermographic simulations to investigate the limitations of the technique.

Winfree, William P.↗

Optimization of Car Body under Constraints of Noise, Vibration, and Harshness (NVH), and Crash

To be competitive on the today's market, cars have to be as light as possible while meeting the Noise, Vibration, and Harshness (NVH) requirements and conforming to Government-man dated crash survival regulations. The latter are difficult to meet because they involve very compute-intensive, nonlinear analysis, e.g., the code RADIOSS capable of simulation of the dynamics, and the geometrical and material nonlinearities of a thin-walled car structure in crash, would require over 12 days of elapsed time for a single design of a 390K elastic degrees of freedom model, if executed on a single processor of the state-of-the-art SGI Origin2000 computer. Of course, in optimization that crash analysis would have to be invoked many times. Needless to say, that has rendered such optimization intractable until now. The car finite element model is shown. The advent of computers that comprise large numbers of concurrently operating processors has created a new environment wherein the above optimization, and other engineering problems heretofore regarded as intractable may be solved. The procedure, shown, is a piecewise approximation based method and involves using a sensitivity based Taylor series approximation model for NVH and a polynomial response surface model for Crash. In that method the NVH constraints are evaluated using a finite element code (MSC/NASTRAN) that yields the constraint values and their derivatives with respect to design variables. The crash constraints are evaluated using the explicit code RADIOSS on the Origin 2000 operating on 256 processors simultaneously to generate data for a polynomial response surface in the design variable domain. The NVH constraints and their derivatives combined with the response surface for the crash constraints form an approximation to the system analysis (surrogate analysis) that enables a cycle of multidisciplinary optimization within move limits. In the inner loop, the NVH sensitivities are recomputed to update the NVH approximation model while keeping the Crash response surface constant. In every outer loop, the Crash response surface approximation is updated, including a gradual increase in the order of the response surface and the response surface extension in the direction of the search. In this optimization task, the NVH discipline has 30 design variables while the crash discipline has 20 design variables. A subset of these design variables (10) are common to both the NVH and crash disciplines. In order to construct a linear response surface for the Crash discipline constraints, a minimum of 21 design points would have to be analyzed using the RADIOSS code. On a single processor in Origin 2000 that amount of computing would require over 9 months! In this work, these runs were carried out concurrently on the Origin 2000 using multiple processors, ranging from 8 to 16, for each crash (RADIOSS) analysis. Another figure shows the wall time required for a single RADIOSS analysis using varying number of processors, as well as provides a comparison of 2 different common data placement procedures within the allotted memories for each analysis. The initial design is an infeasible design with NVH discipline Static Torsion constraint violations of over 10%. The final optimized design is a feasible design with a weight reduction of 15 kg compared to the initial design. This work demonstrates how advanced methodology for optimization combined with the technology of concurrent processing enables applications that until now were out of reach because of very long time-to-solution.

Kodiyalam, Srinivas↗

Parallel Study of HEND, RAD, and DAN Instrument Response to Martian Radiation and Surface Conditions

Nuclear detection methods are being used to understand the radiation environment at Mars. JPL (Jet Propulsion Laboratory) assets on Mars include: Orbiter -2001 Mars Odyssey [High Energy Neutron Detector (HEND)]; Mars Science Laboratory Rover -Curiosity [(Radiation Assessment Detector (RAD); Dynamic Albedo Neutron (DAN))]. Spacecraft have instruments able to detect ionizing and non-ionizing radiation. Instrument response on orbit and on the surface of Mars to space weather and local conditions [is discussed] - Data available at NASA-PDS (Planetary Data System).

Mars↗

Progressive Response Surfaces

Response surface functions are often used as simple and inexpensive replacements for computationally expensive computer models that simulate the behavior of a complex system over some parameter space. Progressive response surfaces are ones that are built up progressively as global information is added from new sample points in the parameter space. As the response surfaces are globally upgraded based on new information, heuristic indications of the convergence of the response surface approximation to the exact (fitted) function can be inferred. Sampling points can be incrementally added in a structured fashion, or in an unstructured fashion. Whatever the approach, at least in early stages of sampling it is usually desirable to sample the entire parameter space uniformly. At later stages of sampling, depending on the nature of the quantity being resolved, it may be desirable to continue sampling uniformly over the entire parameter space (Progressive response surfaces), or to switch to a focusing/economizing strategy of preferentially sampling certain regions of the parameter space based on information gained in early stages of sampling (Adaptive response surfaces). Here we consider Progressive response surfaces where a balanced indication of global response over the parameter space is desired.We use a variant of Moving Least Squares to fit and interpolate structured and unstructured point sets over the parameter space. On a 2-D test problem we compare response surface accuracy for three incremental sampling methods: Progressive Lattice Sampling; Simple-Random Monte Carlo; and Halton Quasi-Monte-Carlo sequences. We are ultimately after a system for constructing efficiently upgradable response surface approximations with reliable error estimates.

Romero, V. J.↗

The Photogrammetry Cube

We can determine distances between objects and points of interest in 3-D space to a useful degree of accuracy from a set of camera images by using multiple camera views and reference targets in the camera s field of view (FOV). The core of the software processing is based on the previously developed foreign-object debris vision trajectory software (see KSC Research and Technology 2004 Annual Report, pp. 2 5). The current version of this photogrammetry software includes the ability to calculate distances between any specified point pairs, the ability to process any number of reference targets and any number of camera images, user-friendly editing features, including zoom in/out, translate, and load/unload, routines to help mark reference points with a Find function, while comparing them with the reference point database file, and a comprehensive output report in HTML format. In this system, scene reference targets are replaced by a photogrammetry cube whose exterior surface contains multiple predetermined precision 2-D targets. Precise measurement of the cube s 2-D targets during the fabrication phase eliminates the need for measuring 3-D coordinates of reference target positions in the camera's FOV, using for example a survey theodolite or a Faroarm. Placing the 2-D targets on the cube s surface required the development of precise machining methods. In response, 2-D targets were embedded into the surface of the cube and then painted black for high contrast. A 12-inch collapsible cube was developed for room-size scenes. A 3-inch, solid, stainless-steel photogrammetry cube was also fabricated for photogrammetry analysis of small objects.

Source record↗

Advances in Aero-Propulsive Modeling for Fixed-Wing and eVTOL Aircraft Using Experimental Data

Small unmanned aircraft and electric vertical takeoff and landing (eVTOL) aircraft have recently emerged as vehicles able to perform new missions and stimulate future air transportation methods. This dissertation presents several system identification research advancements for these modern aircraft configurations enabling accurate mathematical model development for flight dynamics simulations based on wind-tunnel and flight-test data. The first part of the dissertation focuses on advances in flight-test system identification methods using small, fixed-wing, remotely-piloted, electric, propeller-driven aircraft. A generalized approach for flight dynamics model development for small fixed-wing aircraft from flight data is described and is followed by presentation of novel flight-test system identification applications, including: aero-propulsive model development for propeller aircraft and nonlinear dynamic model identification without mass properties. The second part of the dissertation builds on established fixed-wing and rotary-wing aircraft system identification methods to develop modeling strategies for transitioning, distributed propulsion, eVTOL aircraft. Novel wind-tunnel experiment designs and aero-propulsive modeling approaches are developed using a subscale, tandem tilt-wing, eVTOL aircraft, leveraging design of experiments and response surface methodology techniques. Additionally, a method applying orthogonal phase-optimized multisine input excitations to aircraft control effectors in wind-tunnel testing is developed to improve test efficiency and identified model utility. Finally, the culmination of this dissertation is synthesis of the techniques described throughout the document to form a flight-test system identification approach for eVTOL aircraft that is demonstrated using a high-fidelity flight dynamics simulation. The research findings highlighted throughout the dissertation constitute substantial progress in efficient empirical aircraft modeling strategies that are applicable to many current and future aeronautical vehicles enabling accurate flight simulation development, which can subsequently be used to foster advancement in many other pertinent technology areas.

system identification↗

Electromagnetic Response of Superconducting RF Cavities

Recently reported anomalies in frequency shift of order δƒ ~ 0.1–10 kHz for Niobium SRF cavities in a narrow temperature region near T c ≃ 9 K are sensitive to the surface preparation and Nitrogen doping. We developed methods for calculating the surface current response of Nb SRF cavities, as well as the resonant frequency shift and quality factor, as functions of temperature, frequency and disorder based on the Keldysh formulation of the theory for nonequilibrium superconductivity coupled with Maxwell’s theory and boundary conditions for the response functions and electromagnetic field. We show that the anomaly is sensitive to impurity disorder. Our results for the anomaly in the frequency shift are in good agreement with experimental results for Nb with an anisotropic gap on the Fermi surface and inhomogeneous non-magnetic disorder. We also show that the quality factor as a function of the impurity scattering rate is maximum for in SRF cavities with intermediate disorder, ℏ/τΔ~O(1) with the maximum $Q$ decreasing with increasing frequency.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optimization of a Cyclone Using Multiphase Flow Computational Fluid Dynamics

The U.S. Department of Energy National Energy Technology Laboratory's (NETL) 50 kW th chemical looping reactor (CLR) has an underperforming cyclone, which was designed using empirical correlations. To improve the performance of this cyclone using computational fluid dynamics (CFD)-based modeling simulations, four critical design parameters including the vortex tube radius and length, barrel radius, and the inlet width and height were optimized. Here, NETL's open source multiphase flow with interphase exchange (MFiX) CFD code has been used to model a series of cyclones by systematically varying the geometric design parameters. To perform the optimization process, the surrogate modeling and sensitivity analysis followed by the optimization capability in nodeworks was used. The basic methodology for the process is to employ a statistical design of experiments (DOE) method to generate sampling simulations that fill the design space. Corresponding CFD models are then created, executed, and postprocessed. A response surface is created to characterize the relationship between input parameters and the quantities of interest (QoI). Finally, the CFD-surrogate is used by an optimization method to find the optimal design condition based on the objective and constraints prescribed. The resulting optimal cyclone has a larger diameter and longer vortex tube, a larger diameter barrel, and a taller and narrower solids inlet. The improved design has a predicted pressure drop 11 times lower than the original design while reducing the mass loss by a factor of 2.3.

42 ENGINEERING↗

Rotational symmetry relation for efficient response function generation in the coarse mesh transport method COMET

The coarse mesh transport code COMET is a continuous energy hybrid stochastic-deterministic neutronics solver with high fidelity and formidable computation speed in solving reactor core problems. Its method is based on the incident flux expansion theory. In this work, we take advantage of the local geometric symmetry in many reactor cores lattices (e.g., fuel lattices and reflector blocks) to develop relations among the flux response expansion coefficients for symmetric surfaces to further improve the computational efficiency of the COMET response function generation tool (method). This is done by a rigorous derivation of the transformation matrices for the angular and spatial expansion moments resulting from a rotation of a coarse mesh by an arbitrary angle. The relations for the response coefficients for the symmetric surfaces can be then written as the Kronecker product of those transformation matrices. The method is implemented into COMET and tested on two advanced high temperature reactor (AHTR) full-length single assembly benchmark problems. The COMET results using the response function library based on the symmetry relations were compared to those using the library directly generated by continuous energy Monte Carlo for all surfaces. It was found that the eigenvalues and stripe-wise fission densities using the two libraries are in statistical agreement as expected. This indicates that the new method maintains the high fidelity of the original COMET method while improving the computational efficiency in the response function generation by 270% to 400%, depending on the local geometric symmetry. This method also reduces the size of the response function library by the same magnitude (270% to 400%). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

First Principles Nonadiabatic Excited-State Molecular Dynamics in NWChem

Computational simulation of non-adiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons, to name a few. Here, we report an implementation of the fewest-switches surface hopping algorithm in the NWChem computational chemistry program. Here the surface hopping method is combined with linear-response time-dependent density functional theory calculations of adiabatic excited state potential energy surfaces. To treat quantum transitions between arbitrary electronic Born{Oppenheimer states, we have implemented both numerical and analytical differentiation schemes for derivative non-adiabatic couplings. A numerical approach for the time-derivative non-adiabatic couplings together with an analytical method for calculating non-adiabatic coupling vectors is an efficient combination for surface hopping approaches. Additionally, electronic decoherence schemes and a state reassigned unavoided crossings algorithm are also implemented to improve the accuracy of the simulated dynamics and to handle trivial unavoided crossings. We apply our code to study the ultrafast decay of photoexcited benzene, including a detailed analysis of the potential energy surface, population decay time scales, and vibrational coordinates coupled to the excitation dynamics. The development presented in this work is a baseline for future implementations of more sophisticated frameworks for simulating of non-adiabatic molecular dynamics in NWChem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Processing Near-Infrared Imagery of the Orion Heatshield During EFT-1 Hypersonic Reentry

The Scientifically Calibrated In-Flight Imagery (SCIFLI) team captured high-resolution, calibrated, near-infrared imagery of the Orion capsule during atmospheric reentry of the EFT-1 mission. A US Navy NP-3D aircraft equipped with a multi-band optical sensor package, referred to as Cast Glance, acquired imagery of the Orion capsule's heatshield during a period when Orion was slowing from approximately Mach 10 to Mach 7. The line-of-sight distance ranged from approximately 65 to 40 nmi. Global surface temperatures of the capsule's thermal heatshield derived from the near-infrared intensity measurements complemented the in-depth (embedded) thermocouple measurements. Moreover, these derived surface temperatures are essential to the assessment of the thermocouples' reliance on inverse heat transfer methods and material response codes to infer the surface temperature from the in-depth measurements. The paper describes the image processing challenges associated with a manually-tracked, high-angular rate air-to-air observation. Issues included management of significant frame-to-frame motions due to both tracking jerk and jitter as well as distortions due to atmospheric effects. Corrections for changing sky backgrounds (including some cirrus clouds), atmospheric attenuation, and target orientations and ranges also had to be made. The image processing goal is to reduce the detrimental effects due to motion (both sensor and capsule), vibration (jitter), and atmospherics for image quality improvement, without compromising the quantitative integrity of the data, especially local intensity (temperature) variations. The paper will detail the approach of selecting and utilizing only the highest quality images, registering several co-temporal image frames to a single image frame to the extent frame-to-frame distortions would allow, and then co-adding the registered frames to improve image quality and reduce noise. Using preflight calibration data, the registered and averaged infrared intensity images were converted to surface temperatures on the Orion capsule's heatshield. Temperature uncertainties will be discussed relative to uncertainties of surface emissivity and atmospheric transmission loss. Comparison of limited onboard surface thermocouple data to the image derived surface temperature will be presented.

Spisz, Thomas S.↗

Wavelet analysis introduction and application to radar scattering from water waves

Double-parameter expansions, leading to an instantaneous multiscale analysis, are used to obtain time-frequency graphs from time-series data for interpretation of scatterometer response to sea surface geometry. This method is used to define, as well as possible, one-dimensional signal input (elevation and slope). However, some important fundamental parameters, such as the asymmetric development of short scales in crests and troughs, the long-wave amplitude modulation, the different zero-crossings, and so forth, are easily extracted from the wavelet images.

Chapron, B.↗

The Launch Systems Operations Cost Model

One of NASA's primary missions is to reduce the cost of access to space while simultaneously increasing safety. A key component, and one of the least understood, is the recurring operations and support cost for reusable launch systems. In order to predict these costs, NASA, under the leadership of the Independent Program Assessment Office (IPAO), has commissioned the development of a Launch Systems Operations Cost Model (LSOCM). LSOCM is a tool to predict the operations & support (O&S) cost of new and modified reusable (and partially reusable) launch systems. The requirements are to predict the non-recurring cost for the ground infrastructure and the recurring cost of maintaining that infrastructure, performing vehicle logistics, and performing the O&S actions to return the vehicle to flight. In addition, the model must estimate the time required to cycle the vehicle through all of the ground processing activities. The current version of LSOCM is an amalgamation of existing tools, leveraging our understanding of shuttle operations cost with a means of predicting how the maintenance burden will change as the vehicle becomes more aircraft like. The use of the Conceptual Operations Manpower Estimating Tool/Operations Cost Model (COMET/OCM) provides a solid point of departure based on shuttle and expendable launch vehicle (ELV) experience. The incorporation of the Reliability and Maintainability Analysis Tool (RMAT) as expressed by a set of response surface model equations gives a method for estimating how changing launch system characteristics affects cost and cycle time as compared to today's shuttle system. Plans are being made to improve the model. The development team will be spending the next few months devising a structured methodology that will enable verified and validated algorithms to give accurate cost estimates. To assist in this endeavor the LSOCM team is part of an Agency wide effort to combine resources with other cost and operations professionals to support models, databases, and operations assessments.

Prince, Frank A.↗

Bayesian Inference in the Modern Design of Experiments

This paper provides an elementary tutorial overview of Bayesian inference and its potential for application in aerospace experimentation in general and wind tunnel testing in particular. Bayes Theorem is reviewed and examples are provided to illustrate how it can be applied to objectively revise prior knowledge by incorporating insights subsequently obtained from additional observations, resulting in new (posterior) knowledge that combines information from both sources. A logical merger of Bayesian methods and certain aspects of Response Surface Modeling is explored. Specific applications to wind tunnel testing, computational code validation, and instrumentation calibration are discussed.

DeLoach, Richard↗

Effects of Hybrid Flow Control on a Normal Shock Boundary-Layer Interaction

Hybrid flow control, a combination of micro-ramps and steady micro-jets, was experimentally investigated in the 15x15 cm Supersonic Wind Tunnel at the NASA Glenn Research Center. A central composite design of experiments method, was used to develop response surfaces for boundary-layer thickness and reversed-flow thickness, with factor variables of inter-ramp spacing, ramp height and chord length, and flow injection ratio. Boundary-layer measurements and wall static pressure data were used to understand flow separation characteristics. A limited number of profiles were measured in the corners of the tunnel to aid in understanding the three-dimensional characteristics of the flowfield.

Hirt, Stefanie M.↗

Statistical Methodologies to Integrate Experimental and Computational Research

Development of advanced algorithms for simulating engine flow paths requires the integration of fundamental experiments with the validation of enhanced mathematical models. In this paper, we provide an overview of statistical methods to strategically and efficiently conduct experiments and computational model refinement. Moreover, the integration of experimental and computational research efforts is emphasized. With a statistical engineering perspective, scientific and engineering expertise is combined with statistical sciences to gain deeper insights into experimental phenomenon and code development performance; supporting the overall research objectives. The particular statistical methods discussed are design of experiments, response surface methodology, and uncertainty analysis and planning. Their application is illustrated with a coaxial free jet experiment and a turbulence model refinement investigation. Our goal is to provide an overview, focusing on concepts rather than practice, to demonstrate the benefits of using statistical methods in research and development, thereby encouraging their broader and more systematic application.

Parker, P. A.↗

Alternating Conditional Expectations: Introducing a Non‐Parametric Statistical Method to Interpret Long‐Term Greenhouse Gas Flux Measurements Over Semi‐Arid and Wetland Ecosystems

Abstract We explore the potential of using a non‐parametric statistical method called Alternating Conditional Expectations, ACE, to quantify functional relationships in biogeosciences. Here, ACE is used to quantify the non‐linear and multi‐faceted responses of greenhouse gas fluxes to a set of biophysical forcings, when the shapes of those response surfaces are unknown. We evaluated the statistical method over two contrasting ecosystems and two contrasting time steps. One case involved quantifying the biophysical controls of water vapor and carbon dioxide (CO 2 ) fluxes over a semi‐arid oak savanna using daily integrated fluxes. The other case evaluated the responses of CO 2 and methane (CH 4 ) flux measurements to a set of biophysical forcings at a restored tidal wetland using thirty‐minute averages. The statistical model, based on 4 independent variables, explained up over 90% of the variation in daily integrated flux densities of water vapor and net carbon dioxide exchange at the savanna site. This fit was defined by distinct non‐linear responses to such drivers as gross primary production, photosynthetically active radiation, air temperature, vapor pressure deficit and soil moisture. At the tidal wetland site, we evaluated net carbon dioxide and methane fluxes with short‐term measurements to capture the influence of rising and falling tides and seasonality in biological activity. The statistical model defined the shape of the forcing of fluxes due to the roles of carbon exudates, water table depth, oxygen level in the water column, temperature and vegetation status. The statistical fits of the greenhouse gas fluxes were less precise than the savanna case. The fetch varies on a run‐to‐run basis as it is comprised of a heterogeneous mosaic of open water and vegetation. Furthermore, it is difficult to monitor the environmental conditions of the archaea and bacteria in the sediments that produce methane and carbon dioxide.

Environmental Sciences & Ecology↗