Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Response surface model”

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 217 records · Page 12

Effects of surface anomalies on certain model generated meteorological histories

The Mintz-Arakawa 2-level general circulation model has been used in a series of experiments to compute the response of the model atmosphere to: (1) a positive sea surface temperature anomaly in the North Pacific Ocean in summer and in winter, (2) an identical anomaly in the South Pacific Ocean in the Southern Hemisphere winter, and (3) anomalous northward and southward displacements of the Northern Hemisphere snow line over the continents. In each case computations were carried out for 90 forecast days. Results are shown in terms of the differences between anomaly and control histories. Time series of certain regional response indices, including area-average sea level pressure and 600 mb circulation indices, as well as 30-day mean sea level pressure maps are used in the analysis. Of particular interest is the evidence of significant interhemisphere influence.

Spar, J.↗

Transfer function models for using empirical and physics-based simulation signal response data

In many situations, real or induced flaws such as tight cracks with known morphology cannot be manufactured in part geometry specimens or in real parts. Typically, surface fatigue cracks are manufactured in simple geometry specimens such as flat plates, dog-bone shaped flat or cylindrical specimens. If a nondestructive evaluation (NDE)technique is required to provide a reliably detectable flaw size, denoted asa90/95, for induced flaws in a part, then a direct method for qualifying the NDE procedure is to use appropriate induced flaw specimens and perform NDE procedure demonstration on the specimens. Probability of detection (POD)analysis of the empirical data may provide estimation of a90/95. This approach is described as direct POD demonstration testing, which may follow guidelines of MIL-HDBK-1823. This paper considers a case, where embedded tight cracklike induced flaws are to be detected reliably using a signal response based NDE procedure. Here, it is assumed that it is not practical to make surface or embedded induced flaw specimens in part geometry or configuration. Therefore, a direct POD demonstration testing cannot be undertaken. It is also assumed that simulation of signal response is possible for both surface and embedded induced flaws in part geometry specimens using a physics-based model. The proposed approach for NDE procedure qualification uses artificial flaws in simple geometry and part geometry specimens, and induced flaws in the same type of simple geometry specimens. Signal response data is taken on all sets of artificial and induced flaws in simple geometry and part geometry specimens. Moreover, simulated signal response data is generated for surface and embedded flaws. Thus, a case of five signal response versus flaw size datasets is considered. Three of the datasets are empirical and two datasets are physics model-based simulation datasets. A method of devising and using transfer function calculation dataset blocks to estimate the either the reliably detectable flaw size or the demonstration flaw size is provided.

Nondestructive evaluation↗

Radiant Heating Tests on Meteoritic Material

A new effort geared toward modeling the physics of meteor entry and break-up is underway at NASA Ames Research Center. This is part of a broader interdisciplinary effort on providing physics-based risk assessment models for potentially hazardous objects. As part of the entry modeling task we are seeking to improve our understanding of, among other things, the ablation of meteoric material during high speed entry into earths atmosphere. Meteoroid entry differs greatly in some key respects from spacecraft entry modeling. First, the aerothermal environment at these high velocities (18 km/s) is dominated by radiation. Second, meteoroids less than, say, 50m in size, will likely lose a significant portion of their mass during the high-speed atmospheric entry process due to vaporization as well as melting and spallation. The mass of the object, in turn, directly affects the amount of energy deposited in the atmosphere, and therefore the amount of damage done. Thus it is important for us to understand and be able to model this process in greater detail in order to assess the hazard posed by these objects.In this presentation, we first give an overview of the simple ablation models that are typically used in meteor entry calculations. Additionally, we will present a new model which utilizes a similar approach to what is typically done for spacecraft TPS response modeling. This uses an equilibrium assumption near the surface to compute the ablation rate, as is done in heritage material response codes. Next we describe the radiant heating experiment which uses a high-powered laser to emulate the radiation dominated heating environment experienced by meteoroids during atmospheric entry. The facility the Laser Hardened Material Evaluation Laboratory (LHMEL) permits us to expose samples of meteoritic material to heating rates in excess of 100kW/sq.cm. An overview of the experimental set-up and test plan for the initial exploratory campaign at this facility will be given. Then we present both qualitative and quantitative results from this initial test series. Comparisons between predicted and measured ablation rates suggest that there may be significant blockage of the incident beam by the ablation plume. Furthermore, melt is shown to be a significant ablation mechanics, even at high heating rates (16 kW/sq.cm). Finally, comparisons between the phenomenology of the ablation of terrestrial rocks namely, basalt -- to that of meteorites show very different behavior. This is shown to likely be due in part to the effect of composition on the melt viscosity.

Stern, Eric↗

The average impulse response of a rough surface and its applications

This paper is concerned with the theoretical model for short pulse scattering from a statistically random planar surface with particular application to current state of the art radar altimetry. A short review of the assumptions inherent in the convolutional model is presented. Simplified expressions are obtained for both the impulse response and the average backscattered power for near normal incidence under the assumptions common to satellite radar altimetry systems. In particular, it is shown that the conventional two-dimensional surface integration can be reduced to a closed form solution. Two applications of these results are presented relative to radar altimetry, namely, radar antenna pointing angle determination and altitude bias correction for pointing angle and surface roughness effects. It is also shown that these results have direct application to the analysis of the two frequency system proposed by Weissman, and a possible combined long pulse altimeter and two frequency system is suggested.

Brown, G. S.↗

Numerical simulation of the atmospheric response to equatorial Pacific sea surface temperature anomalies

There exists observational and modeling evidence to the effect that sea surface temperature (SST) anomalies in the equatorial Pacific have a significant effect upon wintertime climate in the Pacific/North American sector of the hemisphere. In connection with the emergence of a more definitive observational basis for designing general circulation model (GCM) experiments to investigate the atmospheric response to equatorial SST anomalies, it was recommended to conduct a series of GCM experiments. The present study involves a series of extended (1-2 month) integrations.Attention is given to a review of previous GCM experiments, the plan for the present study, the general circulation model and its simulation of the winter and summer circulation described by Shukla et al. (1981), tropical precipitation, upper level flow, low-level circulation and temperature, and the tropical tropospheric temperature.

Shukla, J.↗

Comparison of Transmission Line Methods for Surface Acoustic Wave Modeling

Surface Acoustic Wave (SAW) technology is low cost, rugged, lightweight, extremely low power and can be used to develop passive wireless sensors. For these reasons, NASA is investigating the use of SAW technology for Integrated Vehicle Health Monitoring (IVHM) of aerospace structures. To facilitate rapid prototyping of passive SAW sensors for aerospace applications, SAW models have been developed. This paper reports on the comparison of three methods of modeling SAWs. The three models are the Impulse Response Method (a first order model), and two second order matrix methods; the conventional matrix approach, and a modified matrix approach that is extended to include internal finger reflections. The second order models are based upon matrices that were originally developed for analyzing microwave circuits using transmission line theory. Results from the models are presented with measured data from devices. Keywords: Surface Acoustic Wave, SAW, transmission line models, Impulse Response Method.

Wilson, William↗

Dynamically morphing microchannels in liquid crystal elastomer coatings containing disclinations

Liquid crystal elastomers (LCEs) hold a major promise as a versatile material platform for smart soft coatings, since their orientational order can be predesigned to program a desired dynamic profile. In this work, we introduce temperature-responsive dynamic coatings based on LCEs with arrays of singular defects-disclinations that run parallel to the surface. The disclinations form in response to antagonistic patterns of the molecular orientation at the top and bottom surfaces, imposed by the plasmonic mask photoalignment. Upon heating, an initially flat LCE coating develops linear microchannels located above each disclination. The stimulus that causes a non-flat profile of LCE coatings upon heating is the activation force induced by the gradients of molecular orientation around disclinations. To describe the formation of microchannels and their thermal response, we adopt a Frank-Oseen model of disclinations in patterned director field and propose a linear elasticity theory to connect the complex spatially-varying molecular orientation to the displacements of the LCE. The thermo-responsive surface profiles predicted by the theory and by the finite element modeling are in good agreement with the experimental data; in particular, higher gradients of molecular orientation produce a stronger modulation of the coating profile. The elastic theory and the finite element simulations allow us to estimate the material parameter that characterizes the elastomer coating’s response to the thermal activation. The disclinations containing LCEs show potential as soft dynamic coatings with a predesigned responsive surface profile.

36 MATERIALS SCIENCE↗

Critique of Macro Flow/Damage Surface Representations for Metal Matrix Composites Using Micromechanics

Guidance for the formulation of robust, multiaxial, constitutive models for advanced materials is provided by addressing theoretical and experimental issues using micromechanics. The multiaxial response of metal matrix composites, depicted in terms of macro flow/damage surfaces, is predicted at room and elevated temperatures using an analytical micromechanical model that includes viscoplastic matrix response as well as fiber-matrix debonding. Macro flow/damage surfaces (i.e., debonding envelopes, matrix threshold surfaces, macro 'yield' surfaces, surfaces of constant inelastic strain rate, and surfaces of constant dissipation rate) are determined for silicon carbide/titanium in three stress spaces. Residual stresses are shown to offset the centers of the flow/damage surfaces from the origin and their shape is significantly altered by debonding. The results indicate which type of flow/damage surfaces should be characterized and what loadings applied to provide the most meaningful experimental data for guiding theoretical model development and verification.

Lissenden, Cliff J.↗

A physics-driven Σ-Y atomization model for heavy-duty engine simulations

The atomization of a liquid jet is a multi-scale and multi-physics problem of interest for many engineering applications. Particularly, it drives the fuel-oxidizer mixing and dictates the efficiency of combustion engines. High-fidelity multi-phase simulations remain challenging due to the excessive computational cost required to capture all the atomization scales. Thus, atomization models are necessary to represent sub-grid liquid structures. Here, in this work, a modified Σ-Y model in the context of Eulerian-Lagrangian Spray Atomization (ELSA) is used to transport the surface area density of the spray. The model's predictive performance is assessed under various operating conditions relevant to heavy-duty engines using the Engine Combustion Network (ECN) Spray C and Spray D research-grade injectors. The classic droplet collision formulation of the Σ-Y model alone does not replicate the response of the measured spray surface area to changes in injector, ambient pressure, and injection pressure, requiring individual tuning of the model's parameters. Instead, a transition between dense and dilute spray breakup mechanisms is proposed in terms of the average droplet spacing. The collision breakup mechanism represents the dilute spray, whereas the droplet size in the dense spray is driven by a competition between the integral scale of turbulence and the balance between the turbulent kinetic energy and the surface energy of the droplets. Such an approach minimizes the requirements for model tuning. Moreover, the role of the constants of a compressible standard k-ε RANS model is assessed in the injector's internal flow and external spray simulation framework, and an updated set is proposed. The results are validated against x-ray radiography and Ultra-Small Angle X-ray Scattering (USAXS) data and highlight the predictive capabilities of the proposed physics-driven Σ-Y model, which is compatible with engine simulation turnaround times.

42 ENGINEERING↗

A Monte Carlo/response surface strategy for sensitivity analysis: application to a dynamic model of vegetative plant growth

We describe the application of a strategy for conducting a sensitivity analysis for a complex dynamic model. The procedure involves preliminary screening of parameter sensitivities by numerical estimation of linear sensitivity coefficients, followed by generation of a response surface based on Monte Carlo simulation. Application is to a physiological model of the vegetative growth of soybean plants. The analysis provides insights as to the relative importance of certain physiological processes in controlling plant growth. Advantages and disadvantages of the strategy are discussed.

Non-NASA Center↗

Sensitivity Analysis of Particle-In-Cell Modeling Parameters in Settling Bed, Bubbling Fluidized Bed and Circulating Fluidized Bed

The objective of the work presented is to perform a preliminary sensitivity analysis of particle-in-cell (PIC) model parameters when applied to settling bed, bubbling fluidized bed, and circulating fluidized bed simulations. These examples correspond to widely different flow conditions commonly seen in chemical engineering applications. Simulations were performed using the PIC method in the open-source software Multiphase Flow with Interphase eXchanges (MFiX) developed by the National Energy Technology Laboratory (NETL). As part of the non-intrusive uncertainty quantification (UQ) analysis, simulation campaigns were generated using Nodeworks. Sampling locations or settings for PIC model parameters were determined using the Latin Hypercube method. Response surfaces were created using radial basis functions (RBF), and Sobol’ indices were estimated to quantify the influence of model parameters on the quantities of interest (QoI). This study marks a first step towards systematically determining optimal ranges for model parameters used in MFiX-PIC. Based on limited experience, it is expected that these values would depend strongly on flow conditions. Given the complexity of the multiphase flow systems under analysis, a non-intrusive UQ based approach is used to identify the most influential parameters in each case. This prior knowledge will help in proposing an effective design of experiments (DoE) and determine optimal parameters through techniques such as deterministic or Bayesian calibration, which will be pursued in the future.

42 ENGINEERING↗

QuadTune version 1: a regional tuner for global atmospheric models

When a new, better-formulated physical parameterization is introduced into a global atmospheric model, aspects of the global model solutions are sometimes degraded. Then, in order to use the new global model to address science questions, there is an incentive to restore its accuracy. Oftentimes this restoration is achieved by tuning of model parameter values. Unfortunately, the retuning process is expensive because characterizing the parameter dependence requires numerous time-consuming global simulations. To reduce the cost of tuning, this manuscript introduces a “poor man's” model tuner, “QuadTune”. QuadTune carves the globe into regions and approximates the model parameter dependence through the use of an uncorrelated quadratic emulator (i.e., response surface). The simplicity of the emulator reduces the required number of global model simulations and aids explainability of tuner behavior. Tuning removes parametric error but leaves behind model structural error. Structural error manifests itself as regional residual biases, such as stubborn biases and tuning trade-offs. To visualize these residual biases, QuadTune's software includes a set of diagnostic plots. This paper illustrates the use of the plots for characterizing residual biases with an example tuning problem.

Environmental sciences↗

The sensitivity of a one-dimensional thermodynamic sea ice model to changes in cloudiness

A thermodynamic sea ice-lead model is used to assess the importance of cloud cover changes to modeled ice thickness. For regions of either permanent multiyear ice or seasonal sea ice, the cloud amount variations have relatively little impact. However, for regions where the presence of summer ice is variable from year to year, the predicted ice thickness is strongly dependent on cloud cover. In general, with a snow covered surface, decreased cloud leads to surface cooling while increased cloud gives rise to a surface warming. For a melting bare ice surface, the reverse occurs. The ice model response time is too long for interannual variations in cloud amount to explain interannual variations in ice thickness and extent. Nevertheless, the implication of the results is that numerical modeling of sea ice distribution requires accurate cloud data or cloud prediction and that trends in cloud cover may lead to significant perturbations in sea ice extent and thickness.

Shine, K. P.↗

Conjugate (solid/fluid) computational fluid dynamics analysis of the space shuttle solid rocket motor nozzle/case and case field joints

Three-dimensional, conjugate (solid/fluid) heat transfer analyses of new designs of the Solid Rocket Motor (SRM) nozzle/case and case field joints are described. The main focus was to predict the consequences of multiple rips (or debonds) in the ambient cure adhesive packed between the nozzle/case joint surfaces and the bond line between the mating field joint surfaces. The models calculate the transient temperature responses of the various materials neighboring postulated flow/leakpaths into, past, and out from the nozzle/case primary O-ring cavity and case field capture O-ring cavity. These results were used to assess if the design was failsafe (i.e., no potential O-ring erosion) and reusable (i.e., no excessive steel temperatures). The models are adaptions and extensions of the general purpose PHOENICS fluid dynamics code. A non-orthogonal coordinate system was employed and 11,592 control cells for the nozzle/case and 20,088 for the case field joints are used with non-uniform distribution. Physical properties of both fluid and solids are temperature dependent. A number of parametric studies were run for both joints with results showing temperature limits for reuse for the steel case on the nozzle joint being exceeded while the steel case temperatures for the field joint were not. O-ring temperatures for the nozzle joint predicted erosion while for the field joint they did not.

Doran, D.↗

A Computational Approach for Probabilistic Analysis of Water Impact Simulations

NASA's development of new concepts for the Crew Exploration Vehicle Orion presents many similar challenges to those worked in the sixties during the Apollo program. However, with improved modeling capabilities, new challenges arise. For example, the use of the commercial code LS-DYNA, although widely used and accepted in the technical community, often involves high-dimensional, time consuming, and computationally intensive simulations. The challenge is to capture what is learned from a limited number of LS-DYNA simulations to develop models that allow users to conduct interpolation of solutions at a fraction of the computational time. This paper presents a description of the LS-DYNA model, a brief summary of the response surface techniques, the analysis of variance approach used in the sensitivity studies, equations used to estimate impact parameters, results showing conditions that might cause injuries, and concluding remarks.

Horta, Lucas G.↗

Damage Propagation Modeling for Aircraft Engine Prognostics

This paper describes how damage propagation can be modeled within the modules of aircraft gas turbine engines. To that end, response surfaces of all sensors are generated via a thermo-dynamical simulation model for the engine as a function of variations of flow and efficiency of the modules of interest. An exponential rate of change for flow and efficiency loss was imposed for each data set, starting at a randomly chosen initial deterioration set point. The rate of change of the flow and efficiency denotes an otherwise unspecified fault with increasingly worsening effect. The rates of change of the faults were constrained to an upper threshold but were otherwise chosen randomly. Damage propagation was allowed to continue until a failure criterion was reached. A health index was defined as the minimum of several superimposed operational margins at any given time instant and the failure criterion is reached when health index reaches zero. Output of the model was the time series (cycles) of sensed measurements typically available from aircraft gas turbine engines. The data generated were used as challenge data for the Prognostics and Health Management (PHM) data competition at PHM 08.

Saxena, Abhinav↗

Adaptive Aft Signature Shaping of a Low-Boom Supersonic Aircraft Using Off-Body Pressures

The design and optimization of a low-boom supersonic aircraft using the state-of-the- art o -body aerodynamics and sonic boom analysis has long been a challenging problem. The focus of this paper is to demonstrate an e ective geometry parameterization scheme and a numerical optimization approach for the aft shaping of a low-boom supersonic aircraft using o -body pressure calculations. A gradient-based numerical optimization algorithm that models the objective and constraints as response surface equations is used to drive the aft ground signature toward a ramp shape. The design objective is the minimization of the variation between the ground signature and the target signature subject to several geometric and signature constraints. The target signature is computed by using a least-squares regression of the aft portion of the ground signature. The parameterization and the deformation of the geometry is performed with a NASA in- house shaping tool. The optimization algorithm uses the shaping tool to drive the geometric deformation of a horizontal tail with a parameterization scheme that consists of seven camber design variables and an additional design variable that describes the spanwise location of the midspan section. The demonstration cases show that numerical optimization using the state-of-the-art o -body aerodynamic calculations is not only feasible and repeatable but also allows the exploration of complex design spaces for which a knowledge-based design method becomes less effective.

Ordaz, Irian↗

Effects of a transient sea surface temperature anomaly on the energetics of the Mintz-Arakawa model atmosphere

The possible response of the atmosphere, as simulated by the two level Mintz-Arakawa global general circulation model, to a transient North Pacific sea surface temperature anomaly is investigated in terms of the energetics both in the spatial and wave number domains. Results indicate that the transient SST variations of reasonable magnitude in the North Pacific Ocean can induce a disturbing effect on the global energetics both in the spatial and wave number domains. The ability of the two level Mintz-Arakawa model to simulate the atmospheric energetics is also examined. Except in the tropics, the model exhibits a reasonable and realistic energy budget.

Chow, S. H.↗