Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “response surface”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Investigation of Surface Skidding Response of PBX 9501 and PBX 9502

Introduction: Assess the hazards of handling polymer bonder explosives(PBX) for the case of skidding with grit particles in the frictional interface. High explosive (HE) materials tested: PBX 9501 and 9502. High melting point materials at the skidding interface between HE and a substrate is known to increase the probability of ignition by creating hot spots not subject to reduced heating from the melting of HE material. Skidding occurred at peak forces of 6-26 kN normal force, and 0.6-3.4 kN of frictional force. Velocities for skidding range from 5-30 m/s, see figure 7. Sparse amounts of 300 μm silica grit is applied loosely to the HE surface before testing. Previous works have tested ignition by grit. Testing parameters include angle of impact, drop height, grit size, grit density, substrate material. This study focuses on pre ignition conditions and detailed thermal imaging of the frictional interface. Conclusions: Temperature hot spots were observed at high points in HE surface topology. These high points may be created by grit particles. Since these points impact first high frictional work is concentrated and quickly increases temperature. Small amounts of grit do not affect the overall frictional force but do concentrate the frictional work to create hot spots.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A comparison of surrogate constitutive models for viscoplastic creep simulation of HT-9 steel

Mechanistic microstructure-informed constitutive models for the mechanical response of polycrystals are a cornerstone of computational materials science. However, as these models become increasingly more complex – often involving coupled differential equations describing the effect of specific deformation modes – their associated computational costs can become prohibitive, particularly in optimization or uncertainty quantification tasks that require numerous model evaluations. To address this challenge, surrogate constitutive models that balance accuracy and computational efficiency are highly desirable. Data-driven surrogate models, that learn the constitutive relation directly from data, have emerged as a promising solution. In this work, we develop two local surrogate models for the viscoplastic response of a steel: a piecewise response surface method and a mixture of experts model. These surrogates are designed to adapt to complex material behavior, which may vary with material parameters or operating conditions. The surrogate constitutive models are applied to creep simulations of HT-9 steel, an alloy of considerable interest to the nuclear energy sector due to its high tolerance to radiation damage, using training data generated from viscoplastic self-consistent (VPSC) simulations. In conclusion, we define a set of test metrics to numerically assess the accuracy of our surrogate models for predicting viscoplastic material behavior, and show that the mixture of experts model outperforms the piecewise response surface method in terms of accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An Investigation into the Surface Skidding Response of PBX 9501 and PBX 9502

While polymer-bonded explosives, PBX, have reduced sensitivity to ignition from mechanical shock or heating compared to conventional explosives, the investigation of the mechanical ignition mechanisms for PBX remains vital to assessing the safety during machining and general handling. Under frictional heating, hot spot generation due to high melting point contaminants in the interface is a suspected source for increased probability of detonation. Here, to investigate such frictional heating, samples of PBX 9501 and PBX 9502 were impacted and skidded against a tangentially moving target and thermally imaged. Temporally resolved temperature and forces were simultaneously measured with and without 300 µm silica grit at the frictional interface. A trend of increasing temperature was observed as the speed and tangential force on the sample increased. Grit particles in the frictional interface were found to act as frictional heat concentration spots owing to the grit's protrusion from the surface and lesser susceptibility to melting compared to the PBX. The coefficient of friction between PBX 9501 and the skidding surface was observed to be largely constant for forces and speeds at the lower end of the test range. In contrast, at high tangential speeds, the coefficient was found to be significantly lower.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Statistical Design of Experiments Enables Rapid Exploration of Perfluorobutane Sulfonate Degradation

The phase-out of long-chain PFAS has led to the proliferation of highly recalcitrant short-chain analogs, and mineralization technologies for these short-chain PFAS are needed to mitigate their deleterious effects on environmental and human health. In this study, we utilize a statistical design of experiments, specifically, response surface methodology, to rapidly evaluate the electrochemical degradation of the short-chain PFAS, perfluorobutane sulfonate (PFBS). We evaluate the impacts of the three primary electrochemical parameters (concentration of PFBS, concentration of supporting electrolyte, and applied current) over multiple orders of magnitude on the three primary reaction outcomes of electrochemical PFBS degradation (incomplete PFBS decomposition, complete PFBS mineralization as fluorine, and anodic energy consumption). Our results correspond with literature and clearly identify the well-known tradeoff between energy consumption and complete mineralization. Intriguingly, partial PFBS decomposition and energy consumption demonstrate nonlinear dependencies in the current/supporting electrolyte concentration space and the current/PFBS concentration space, respectively. These findings highlight the utility of the response surface methodology model to efficiently interrogate a large parameter space, identifying both common results and less-obvious interactions between electrochemical parameters and their influences on reaction outcomes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Probabilistic Evaluation of Geomechanical Risks in CO2 Storage: An Exploration of Caprock Integrity Metrics Using a Multilaminate Model

The probabilistic uncertainty assessment of geomechanical risk—specifically, caprock failure—attributable to CO2 injection, as presented in a simplified hypothetical geological model, was the focus of this study. Our approach amalgamates the implementation of a multilaminate model, the creation of a response surface model in conjunction with the Box–Behnken sampling design, the execution of associated numerical modeling experiments, and the utilization of Monte Carlo simulations. Probability distributions to encapsulate the inherent variability (elastic and mechanical properties of the caprock and reservoir) and uncertainty in prediction estimates (vertical displacement, total strain, and F value) were employed. Our findings reveal that the Young modulus of the caprock is a key factor controlling equivalent total strain but is insufficient as a stand-alone indicator of caprock integrity. It is confirmed that the caprock can accommodate significant deformation without failure, if it possesses a low Young’s modulus and high mechanical strength properties, such as the friction angle and uniaxial compressive strength. Similarly, vertical displacement was found to be an unreliable indicator for caprock integrity, as caprock failure can occur across a broad spectrum of vertical displacements, particularly when both the Young modulus and mechanical strength properties have wide ranges. This study introduces the F value as the most dependable indicator for caprock failure, although it is a theoretical attribute (the shortest distance between the Mohr circle and the nearest failure envelope used to measure the sensitivity to failure) and not physically measurable in the field. Deviatoric stress levels were found to vary based on stress regimes, with the maximum levels observed under extensive and compressive stress regimes. In conjunction with the use of the response surface method, this study demonstrates the efficacy of the multilaminate framework and the Mohr–Coulomb constitutive model in providing a simplified, yet effective, probabilistic model of the mechanical behavior of caprock failure, reducing mathematical and computational complexities.

Energy & Fuels↗
DOE OSTI2023

High-moisture pelleting of corn stover using pilot- and commercial-scale systems: Impact of moisture content, L / D ratio and hammer mill screen size on pellet quality and energy consumption

Here, high-moisture pelleting of corn stover was tested in pilot- and commercial-scale systems. For the pilot-scale study, a 6.35 mm screen size hammer mill ground corn stover was reconditioned to moisture between 14 and 29.6% (w.b.), and 5, 7 and 9 L/D ratio pellet die were used. At 29% corn stover moisture content and L/D ratios of 5 and 9, the bulk density and durability of the pellets were ≤450 kg/m 3 and <90% and >540 kg/m 3 and >90%, respectively. Increasing the feedstock moisture and L/D ratio increased the pelleting energy. Increasing the hammer mill screen size to 11.11 mm reduced the bulk density but not the durability. The response surface models adequately described the pelleting process (R 2 ≥ 0.88), and analysis of variance showed a strong interaction between the process variables and pellet properties. Commercial-scale testing of high-moisture pelleting using 6.35 mm ground corn stover bales in the moisture range 24–26% (w.b.) produced pellets with a density of >675 kg/m 3 and durability of >98%, whereas increasing the hammer mill screen size to 11.11 mm reduced the bulk density by about 100 kg/m 3 but not the durability. The energy consumption of the commercial-scale unit operations for the high-moisture pelleting process was in the range 98–124 kW h/ton, which is 64–72% less energy than conventional pelleting that requires energy about 350 kW h/ton for biomass drying from 30% (w.b.) to 10% (w.b.) moisture content before pelleting. Pellets produced using high-moisture pelleting met International Organization for Standardization and Pellets Fuel Institute international standards. Published 2023. This article is a U.S. Government work and is in the public domain in the USA.

09 BIOMASS FUELS↗
DOE OSTI2022

Throughput Optimization of Molybdenum Carbide Nanoparticle Catalysts in a Continuous Flow Reactor Using Design of Experiments

Transition metal carbides (TMCs) have attracted significant attention because of their applications toward a wide range of catalytic transformations. However, the practicality of their synthesis is still limited because of the harsh conditions in which most TMCs are prepared. Recently, a solution-phase synthesis of phase-pure a-MoC1-x nanoparticles was presented. While this synthetic route yielded nanoparticles with exceptional catalytic performance, the reaction parameter space was not explored, and catalyst throughput was not optimized for scale-up. Continuous flow platforms coupled with statistical design of experiments (DoE) can provide a powerful method for understanding the reaction parameter space for optimizations. Here, we demonstrate the use of statistical DoE in tandem with response surface methodology for a parametric screening analysis to optimize the throughput of a MoC1-x nanoparticle synthesis utilizing a millifluidic flow reactor. A full factorial design was implemented to evaluate four input variables (reaction temperature, flow rate, solvent fraction of oleylamine, and precursor concentration) that carry statistically significant effects on three responses (throughput, residence time, and isolated yield). A Doehlert matrix was implemented to investigate each significant variable at a higher number of levels to optimize throughput. Our results give a nonintuitive set of experimental conditions that resulted in an optimized throughput of 2.2 g h-1. This translates to a 50-fold increase in throughput compared to the previously reported batch method. The catalytic performance of the MoC1-x nanoparticles produced under optimized throughput was demonstrated in the CO2 hydrogenation reaction. This DoE screening analysis and throughput optimization of MoC1-x synthesis open the door to an increased feasibility for scale-up.

design of experiments↗
DOE OSTI2026

Disentangling High Harmonic Generation from Surface and Bulk States of a Topological Insulator

The discovery of topological phases has introduced a new dimension to materials science. Three-dimensional (3D) topological insulators (TIs) are a remarkable class of matter that is insulating in the bulk while hosting conductive topological surface states (TSSs) with unique charge and spin properties. High-order harmonic generation (HHG) has emerged as a powerful tool to probe condensed matter systems by providing insights into their electronic structure and dynamic behavior. Here, we investigate HHG in the prototype 3D-TI Bi$_2$Se$_3$. We demonstrate that the contributions of bulk and surface states to the harmonic emission can be controlled by tuning the thickness of thin film samples. An ultrathin (6 nm) film substantially enhances HHG from the surface states, while the bulk states dominate HHG in a thicker (50 nm) film. By applying a quasi-static terahertz perturbing field, we disentangle the bulk and surface responses and reveal the significant impact of the surface states' shift vector and Berry curvature on HHG. Our study provides effective methods for isolating the optical responses of TSSs from those of the bulk, which opens the door to resolving an ongoing debate regarding whether it is possible to reliably extract topological signatures in HHG.

Atomic Physics (physics.atom-ph)↗
DOE OSTI2025

Surface Energy Balance Responses to Radiative Forcing in the Central Arctic From MOSAiC and Models

The Arctic surface energy budget (SEB) couples the atmosphere with the sea ice, making it useful for both studying surface processes as well as evaluating models. Improved understanding of atmosphere-ice interactions is required to improve models, requiring year-round observations to address seasonally dependent biases. This work uses novel observations from the MOSAiC expedition to quantify the responses of surface fluxes to radiative forcing over sea ice throughout a complete annual cycle. We identify two primary regimes of flux response: an ice growth regime in winter and an ice melt regime in summer. In the growth regime, changes in radiative forcing impact upwelling longwave, sensible heat, and subsurface heat fluxes, whereas in the melt regime changes in radiative forcing primarily alter the amount of melt and subsurface transmission because the surface temperature is fixed. These observed responses of surface fluxes to radiative forcing are used to evaluate seven weather forecast models during the ice growth regime. In most models, the responses of surface fluxes to radiative forcing do not match observations. Many models also have biased downwelling longwave. One model (the Coupled Arctic Forecast System; CAFS) adequately captures both the mean radiative forcing and the flux responses in winter. CAFS is further evaluated against observations spanning the full MOSAiC year, demonstrating sufficient agreement to provide a more generalized understanding of these SEB process relationships across the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Balancing sugar recovery and inhibitor generation during energycane processing: Coupling cryogenic grinding with hydrothermal pretreatment at low temperatures

Pretreatment of lignocellulosic biomass at high temperatures or with oxidizing chemicals generate various inhibitors that restrict the efficient bioconversion of sugars in subsequent steps. The present study systematically investigates individual and combinatorial effects of pretreatment parameters on the generation of inhibitors. Herein, a plot between pretreatment temperature and inhibitor revealed optimum pretreatment temperature for energycane bagasse i.e., 170 °C beyond which total inhibitor production increased exponentially. No inhibitor production was observed on mechanical processing i.e., disk milling/cryogenic grinding of biomass. Evaluation of response surface regression exhibited that biomass solids loading has a significant effect on inhibitor generation at higher temperatures. The concentrations of certain inhibitors such as acetic acid, furfurals, and HMF increased more than 3-folds on doubling the solids loading. Furthermore, a novel low-severity approach of low-temperature hydrothermal pretreatment coupled with cryogenic grinding for lignocellulosic biomasses has been introduced which improved sugar yields while maintaining a low inhibitor concentration.

09 BIOMASS FUELS↗
DOE OSTI2022

Algorithm 1028: VTMOP: Solver for Blackbox Multiobjective Optimization Problems

VTMOP is a Fortran 2008 software package containing two Fortran modules for solving computationally expensive bound-constrained blackbox multiobjective optimization problems. VTMOP implements the algorithm of [32], which handles two or more objectives, does not require any derivatives, and produces well-distributed points over the Pareto front. The first module contains a general framework for solving multiobjective optimization problems by combining response surface methodology, trust region methodology, and an adaptive weighting scheme. The second module features a driver subroutine that implements this framework when the objective functions can be wrapped as a Fortran subroutine. Lastly, support is provided for both serial and parallel execution paradigms, and VTMOP is demonstrated on several test problems as well as one real-world problem in the area of particle accelerator optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

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↗
DOE OSTI2020

Chemical Impacts of Potential CO2 and Brine Leakage on Groundwater Quality with Quantitative Risk Assessment: A Case Study of the Farnsworth Unit

Potential leakage of reservoir fluids is considered a key risk factor for geologic CO2 sequestration (GCS), with concerns of their chemical impacts on the quality of overlying underground sources of drinking water (USDWs). Effective risk assessment provides useful information to guide GCS activities for protecting USDWs. In this study, we present a quantified risk assessment case study of an active commercial-scale CO2-enhanced oil recovery (CO2-EOR) and sequestration field, the Farnsworth Unit (FWU). Specific objectives of this study include: (1) to quantify potential risks of CO2 and brine leakage to the overlying USDW quality with response surface methodology (RSM); and (2) to identify water chemistry indicators for early detection criteria. Results suggest that trace metals (e.g., arsenic and selenium) are less likely to become a risk due to their adsorption onto clay minerals; no-impact thresholds based on site monitoring data could be a preferable reference for early groundwater quality evaluation; and pH is suggested as an indicator for early detection of a leakage. This study may provide quantitative insight for monitoring strategies on GCS sites to enhance the safety of long-term CO2 sequestration.

58 GEOSCIENCES↗
DOE OSTI2020

Modeling, analysis, and code/data validation of DIII-D tokamak divertor experiments on ELM and non-ELM plasma tungsten sputtering erosion

Abstract We analyzed recent DIII-D tokamak tungsten divertor probe experiments using advanced, coupled, sputter erosion/redeposition, plasma, and surface response code packages. Modeling is done for ELMing H-mode, and L-mode plasmas, impinging on various size tungsten deposits on Divertor Material Evaluation System (DiMES) carbon probes. The simulations compute 3D, full kinetic, sub-gyromotion, impurity sputtering and transport, including changes in tungsten surface composition and response due to mixed deuterium and carbon ions irradiation. Per our analysis, ELM (edge localized mode) plasma sputtering in DIII-D mostly involves free-streaming high energy (∼500–1000 eV) D + and C +6 ions, with high near-surface plasma density. L-Mode sputtering is due to impurity sputtering (C, W) only, with lower density. All cases show complete redeposition of tungsten on the divertor, with significant redeposition on the tungsten spots themselves, and low self-sputtering. Comparison of ELM plasma gross tungsten erosion simulation results with in-situ spectroscopic data is good, as are code/data comparisons of net erosion using post-exposure Rutherford backscattering (RBS) data for the L-mode probes. The analysis, extrapolated to a full tungsten divertor, implies low net erosion and negligible plasma contamination from sputtering. These results support the use of high-Z plasma facing surfaces in ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Catalytic Activity and Water Stability of the MgO(111) Surface for 2-Pentanone Condensation

Nanomaterials derived from earth-abundant metal oxides have gained tremendous interest as catalysts; although, water stability remains a challenge. This study examines MgO(111) surfaces for 2-pentanone condensation and their evolution during D2O hydration. Catalyst screening confirmed the high activity of fresh MgO(111) for 2-pentanone condensation relative to conventionally prepared MgO(100). Computational modeling suggests that the (111) surface is readily hydroxylated, and that surface hydroxyls help stabilize the surface and reduce the barrier for 2-pentanone condensation. Vapor-phase D2O hydration after 3 min increased MgO(111) hydroxyls and retained surface area and activity; however, after 1 h, deuteroxide formation reduced the surface area and activity by >30 %. After 24 h, deuteroxide growth slowed down, and surface area and activity remained stable. This suggests MgO(111)-derived hydroxide may be the dominant surface responsible for 2-pentanone condensation following water exposure. Thermal regeneration of the 24-h sample restored 86 % of the surface area and 94 % of the activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Risk Assessment and Management Workflow—An Example of the Southwest Regional Partnership

This paper summarizes the risk assessment and management workflow developed and applied to the Southwest Regional Partnership on Carbon Sequestration (SWP) Phase III Demonstration Project. The risk assessment and management workflow consists of six primary tasks, including management planning, identification, qualitative analysis, quantitative analysis, response planning, and monitoring. Within the workflow, the SWP assembled and iteratively updated a risk registry that identifies risks for all major activities of the project. Risk elements were ranked with respect to the potential impact to the project and the likelihood of occurrence. Both qualitative and quantitative risk analyses were performed. To graphically depict the interactions among risk elements and help building risk scenarios, process influence diagrams were used to represent the interactions. The SWP employed quantitative methods of risk analysis including Response Surface Method (RSM), Polynomial Chaos Expansion (PCE), and the National Risk Assessment Partnership (NRAP) toolset. The SWP also developed risk response planning and performed risk control and monitoring to prevent the risks from affecting the project and ensure the effectiveness of risk management. As part of risk control and monitoring, existing and new risks have been tracked and the response plan was subsequently evaluated. Findings and lessons learned from the SWP’s risk assessment and management efforts will provide valuable information for other commercial geological CO2 storage projects.

Lee, Si-Yong↗