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 199 records · Page 11

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↗

Implicit Coupling Approach for Simulation of Charring Carbon Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver with nonequilibrium gas/surface interaction for simulation of charring carbon ablators can be performed using an implicit approach. The material thermal response code used in this study is the three-dimensional version of Fully Implicit Ablation and Thermal response program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation method. Coupling between the material response and flow codes is performed by solving the surface mass balance in flow solver and the surface energy balance in material response code. Thus, the material surface recession is predicted in flow code, and the surface temperature and pyrolysis gas injection rate are computed in material response code. It is demonstrated that the time-lagged explicit approach is sufficient for simulations at low surface heating conditions, in which the surface ablation rate is not a strong function of the surface temperature. At elevated surface heating conditions, the implicit approach has to be taken, because the carbon ablation rate becomes a stiff function of the surface temperature, and thus the explicit approach appears to be inappropriate resulting in severe numerical oscillations of predicted surface temperature. Implicit coupling for simulation of arc-jet models is performed, and the predictions are compared with measured data. Implicit coupling for trajectory based simulation of Stardust fore-body heat shield is also conducted. The predicted stagnation point total recession is compared with that predicted using the chemical equilibrium surface assumption

Ablation↗

The GLAS physical inversion method for analysis of TIROS N data

The Goddard Laboratory for Atmospheric Sciences approach to temperature sounding from satellite observations is fundamentally different from that of NESS in that heavy reliance is placed on the ability to model accurately the instrumental response to atmospheric and surface conditions, while no use is made of statistical relationships between satellite observations and atmospheric temperature profiles. The method involves starting with a guess set of atmospheric and surface conditions from which expected brightness temperatures for the satellite observations are computed. Then, iterative relaxation of atmospheric and surface conditions is performed according to the difference between observed and computed brightness temperatures until sufficient agreement is reached. In the absence of sufficient agreement, no retrieval is produced for that location.

Susskind, J.↗

Climate-focused Life Cycle Assessments of Biochar Production by an ARTi Pyrolysis Reactor and an Air Burners CharBoss® Air Curtain Incinerator

This report presents a limited, dynamic, consequential life cycle assessment (LCA) to compare the climate impacts of two biochar production methods using wood as a feedstock. The two methods are a pyrolysis reactor supplied by ARTi (Des Moines, IA, https://www.arti.com/) and a T26 CharBoss® air curtain incinerator supplied by Air Burners, Inc. (Palm City, FL, https://airburners.com/). The underlying LCA methodology is described in a chapter by Singh et al. (2024) and implemented in the form of a workbook freely available as online Supplementary Material for the chapter. For the convenience of the reader, a pre-print version of the relevant portions of Singh et al. (2024) is attached as Appendix A to this report. Specific assumptions and calculations to obtain the input parameters used in this LCA for each production method are described in the Methodology section below. This implementation of the LCA considers emissions associated with biomass loading, comminution and conversion, biochar decay in soil, and the production and use of bioenergy generated during the conversion process. The LCA is “limited” in that upstream emissions associated with biomass production, harvest, transportation, and land-use change, as well as embodied emissions in equipment and facilities are not considered. Similarly, downstream emissions from biochar transport and incorporation into soil (i.e., tillage), and the impact of biochar soil amendments on soil greenhouse gas (GHG) emissions (other than CO 2 from biochar decay in soil), soil organic carbon stocks, crop response, and surface albedo are not considered. As the intent is to compare different biochar production methods in a simple unbiased manner, the primary alternative biomass pathway for the LCA is immaculate combustion, which is the hypothetical instantaneous and complete conversion of carbon in the biomass to CO 2 at time zero without generation of any other greenhouse gases or aerosols (GHGAs) or any useful bioenergy. Use of this pathway provides relative values for the production methods and, when the embodied emissions are similar and the same feedstock is used, these relative values are reasonable approximations for those attained with a full LCA.

09 BIOMASS FUELS↗

Long-time dynamic compatibility of two ethylene propylene elastomers with hydrazine

A test method is described for predicting the long-time survivability of elastomers in hydrazine under dynamic stressing conditions. The method selected is based upon the existence and approximate invariance of a 'physical property surface' relating the mechanical response of the elastomer in terms of stress, strain, time, and temperature. The property surface was generated for two selected elastomers (EPT-10 and AF-E-332). A novel carousel testing tank was designed to allow sequential testing of eight tensile specimens immersed in liquid hydrazine within a constant-temperature water bath. The test procedure and data reduction methods used to generate the property surface are described. The utility and validity of these results applied to fatigue and flexure loading to these elastomeric materials over long-time periods are discussed.

Coulbert, C. D.↗

Finite Element and Plate Theory Modeling of Acoustic Emission Waveforms

A comparison was made between two approaches to predict acoustic emission waveforms in thin plates. A normal mode solution method for Mindlin plate theory was used to predict the response of the flexural plate mode to a point source, step-function load, applied on the plate surface. The second approach used a dynamic finite element method to model the problem using equations of motion based on exact linear elasticity. Calculations were made using properties for both isotropic (aluminum) and anisotropic (unidirectional graphite/epoxy composite) materials. For simulations of anisotropic plates, propagation along multiple directions was evaluated. In general, agreement between the two theoretical approaches was good. Discrepancies in the waveforms at longer times were caused by differences in reflections from the lateral plate boundaries. These differences resulted from the fact that the two methods used different boundary conditions. At shorter times in the signals, before reflections, the slight discrepancies in the waveforms were attributed to limitations of Mindlin plate theory, which is an approximate plate theory. The advantages of the finite element method are that it used the exact linear elasticity solutions, and that it can be used to model real source conditions and complicated, finite specimen geometries as well as thick plates. These advantages come at a cost of increased computational difficulty, requiring lengthy calculations on workstations or supercomputers. The Mindlin plate theory solutions, meanwhile, can be quickly generated on personal computers. Specimens with finite geometry can also be modeled. However, only limited simple geometries such as circular or rectangular plates can easily be accommodated with the normal mode solution technique. Likewise, very limited source configurations can be modeled and plate theory is applicable only to thin plates.

Prosser, W. H.↗

Response of stiffened sandwich panels

An increasingly important part of the design of skin-stringer systems, such as those found in flight vehicles, is the ability to reduce surface vibration levels when they are subjected to the dynamic loads of jet, or propeller noise. One common method of reducing the response of these structures is by layered construction using viscoelastic materials. Thus analytical methods to describe such systems are essential.

Slazak, M.↗

Calibration to surface reflectance of terrestrial imaging spectrometry data: Comparison of methods

Many algorithms for spectral analysis of imaging spectroscopy data of the Earth's surface require that the data be calibrated to surface reflectance. Calibration requires removing instrumental response, solar irradiance, atmospheric transmittance, and atmospheric scattering from the radiance detected at the sensor. Depending on the amount of support data, this can be a formidable task. This paper examines four methods of calibration: (1) a radiative transfer model from the University of Colorado (ATREM: Gao and Goetz, 1990; Gao et al., 1992), (2) a MODTRAN-based method developed at the Jet Propulsion Lab by Green et al., (1191), (3) a ground calibration using known sites as standards, and (4) a combined approach using radiative transfer methods and ground calibration. Data from the Airborne Visual and Infra-Red Imaging Spectrometer (AVIRIS) instrument were evaluated from data sets obtained over multiple years and multiple sites.

Clark, Roger N.↗

Dynamic PRA-Based Estimation of PWR Coping Time Using a Surrogate Model for Accident Tolerant Fuel

In this study, we propose an interpolation-based response surface surrogate methodology to manage a large number of scenarios in dynamic probabilistic risk assessment. It adopts the shape Dynamic Time Warping algorithm to cluster the interpolation neighborhood from time series sample data. The interpolation method was adapted from Taylor Kriging to allow a reduced-order model of the Taylor series. In order to demonstrate its applicability to complex issues in risk assessment for nuclear engineering, an example risk response surface to estimate emergency core cooling system (ECCS) criteria for triplex silicon carbide (SiC) accident-tolerant fuel was constructed. The response surface was exploited to estimate the cumulative failure probability of the fuel cladding structure due to the uncertainties in operator actions and safety systems. The functional failures were assessed based on a combination of individual layer failures computed by coupling Risk Analysis Virtual Environment software with a pressurized water reactor 1000-MW(electric) RELAP5 model and the in-house fuel performance assessment module. Results showed that SiC cladding failure probability spiked less than 1 min after a large-break loss-of- coolant accident whenever the current ECCS criteria for Zircaloy-4 (Zr-4) cladding was used. However, it still provides an increased safety margin of three orders of magnitude compared to Zr-4. This positive margin could be utilized to relax active ECCS requirements by allowing deviations of up to 450 s in its actuation time. The proposed surrogate methodology generated a response surface of SiC cladding failure probability reasonably well, with a significant savings of computation time. This methodology is expected to be useful in the analysis of system response with complex uncertainty sources.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hierarchical Modeling and Robust Synthesis for the Preliminary Design of Large Scale Complex Systems

Large-scale complex systems are characterized by multiple interacting subsystems and the analysis of multiple disciplines. The design and development of such systems inevitably requires the resolution of multiple conflicting objectives. The size of complex systems, however, prohibits the development of comprehensive system models, and thus these systems must be partitioned into their constituent parts. Because simultaneous solution of individual subsystem models is often not manageable iteration is inevitable and often excessive. In this dissertation these issues are addressed through the development of a method for hierarchical robust preliminary design exploration to facilitate concurrent system and subsystem design exploration, for the concurrent generation of robust system and subsystem specifications for the preliminary design of multi-level, multi-objective, large-scale complex systems. This method is developed through the integration and expansion of current design techniques: Hierarchical partitioning and modeling techniques for partitioning large-scale complex systems into more tractable parts, and allowing integration of subproblems for system synthesis; Statistical experimentation and approximation techniques for increasing both the efficiency and the comprehensiveness of preliminary design exploration; and Noise modeling techniques for implementing robust preliminary design when approximate models are employed. Hierarchical partitioning and modeling techniques including intermediate responses, linking variables, and compatibility constraints are incorporated within a hierarchical compromise decision support problem formulation for synthesizing subproblem solutions for a partitioned system. Experimentation and approximation techniques are employed for concurrent investigations and modeling of partitioned subproblems. A modified composite experiment is introduced for fitting better predictive models across the ranges of the factors, and an approach for constructing partitioned response surfaces is developed to reduce the computational expense of experimentation for fitting models in a large number of factors. Noise modeling techniques are compared and recommendations are offered for the implementation of robust design when approximate models are sought. These techniques, approaches, and recommendations are incorporated within the method developed for hierarchical robust preliminary design exploration. This method as well as the associated approaches are illustrated through their application to the preliminary design of a commercial turbofan turbine propulsion system. The case study is developed in collaboration with Allison Engine Company, Rolls Royce Aerospace, and is based on the Allison AE3007 existing engine designed for midsize commercial, regional business jets. For this case study, the turbofan system-level problem is partitioned into engine cycle design and configuration design and a compressor modules integrated for more detailed subsystem-level design exploration, improving system evaluation. The fan and low pressure turbine subsystems are also modeled, but in less detail. Given the defined partitioning, these subproblems are investigated independently and concurrently, and response surface models are constructed to approximate the responses of each. These response models are then incorporated within a commercial turbofan hierarchical compromise decision support problem formulation. Five design scenarios are investigated, and robust solutions are identified. The method and solutions identified are verified by comparison with the AE3007 engine. The solutions obtained are similar to the AE3007 cycle and configuration, but are better with respect to many of the requirements.

Koch, Patrick N.↗

Remote sensing of ice phenomena from orbit by signal correlation of multiple receiver responses

The method of signal correlation of microwave responses as applied to the measurement of Earth-surface ice temperatures from orbit is explained and summarized. Ice temperatures are estimated by a correlation function that is derived from the processes of a forward stepwise correlator. Subsets of the post-detected outputs of microwave receiving channels are combined in a multivariate cross-correlation function which operates as a spatial filter and serves to improve the spatial resolution of the thermal gradients in ice structures. The correlator is designed to selectively identify the correlative components among the microwave responses and to strongly suppress or cancel the non-correlative components appearing in the post-detected outputs.

Stacey, J. M.↗

Turbine Engine Hot Section Technology 1986

On October 21 and 22, 1986, the Turbine Engine Hot Section Technology (HOST) Project of the NASA Lewis Research Center sponsored a workshop to discuss current research pertinent to turbine engine durability problems. Presentations were made concerning the hot-section environment and the behavior of combustion liners, turbine blades, and turbine vanes. The presentations were divided Into six sessions: Instrumentation, Combustion, Turbine Heat Transfer, Structural Analysis, Fatigue and Fracture, and Surface Protection. The principal objective of each session was to disseminate research results to date. Topics discussed included modeling of thermal and fluid-flow phenomena, structural analysis, fatigue and fracture, surface protective coatings, constitutive behavior of materials, stress-strain response, and life-prediction methods. Researchers from industry, academia, and government presented results of their work sponsored by the HOST Project. This publication contains extended abstracts and visual material presented during the workshop.

Turbine Engine Hot Section Technology (HOST)↗

Universal first-order reliability concept applied to semistatic structures

A reliability design concept was developed for semistatic structures which combines the prevailing deterministic method with the first-order reliability method. The proposed method surmounts deterministic deficiencies in providing uniformly reliable structures and improved safety audits. It supports risk analyses and reliability selection criterion. The method provides a reliability design factor derived from the reliability criterion which is analogous to the current safety factor for sizing structures and verifying reliability response. The universal first-order reliability method should also be applicable for air and surface vehicles semistatic structures.

Verderaime, V.↗

Terahertz NDE for Under Paint Corrosion Detection and Evaluation

Corrosion under paint is not visible until it has caused paint to blister, crack, or chip. If corrosion is allowed to continue then structural problems may develop. Identifying corrosion before it becomes visible would minimize repairs and costs and potential structural problems. Terahertz NDE imaging under paint for corrosion is being examined as a method to inspect for corrosion by examining the terahertz response to paint thickness and to surface roughness.

Anastasi, Robert F.↗

Viscosity Measurement Using Drop Coalescence in Microgravity

We present in here validation studies of a new method for application in microgravity environment which measures the viscosity of highly viscous undercooled liquids using drop coalescence. The method has the advantage of avoiding heterogeneous nucleation at container walls caused by crystallization of undercooled liquids during processing. Homogeneous nucleation can also be avoided due to the rapidity of the measurement using this method. The technique relies on measurements from experiments conducted in near zero gravity environment as well as highly accurate analytical formulation for the coalescence process. The viscosity of the liquid is determined by allowing the computed free surface shape relaxation time to be adjusted in response to the measured free surface velocity for two coalescing drops. Results are presented from two sets of validation experiments for the method which were conducted on board aircraft flying parabolic trajectories. In these tests the viscosity of a highly viscous liquid, namely glycerin, was determined at different temperatures using the drop coalescence method described in here. The experiments measured the free surface velocity of two glycerin drops coalescing under the action of surface tension alone in low gravity environment using high speed photography. The liquid viscosity was determined by adjusting the computed free surface velocity values to the measured experimental data. The results of these experiments were found to agree reasonably well with the known viscosity for the test liquid used.

Antar, Basil N.↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗