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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Probing topological phase transitions using high-harmonic generation

We report the prediction and realization of topological insulators have sparked great interest in experimental approaches to the classification of materials. The phase transition between non-trivial and trivial topological states is important, not only for basic materials science but also for next-generation technology, such as dissipation-free electronics. It is therefore crucial to develop advanced probes that are suitable for a wide range of samples and environments. Here we demonstrate that circularly polarized laser-field-driven high-harmonic generation is distinctly sensitive to the non-trivial and trivial topological phases in the prototypical three-dimensional topological insulator bismuth selenide. The phase transition is chemically initiated by reducing the spin–orbit interaction strength through the substitution of bismuth with indium atoms. We find strikingly different high-harmonic responses of trivial and non-trivial topological surface states that manifest themselves as a conversion efficiency and elliptical dichroism that depend both on the driving laser ellipticity and the crystal orientation. The origins of the anomalous high-harmonic response are corroborated by calculations using the semiconductor optical Bloch equations with pairs of surface and bulk bands. As a purely optical approach, this method offers sensitivity to the electronic structure of the material, including its nonlinear response, and is compatible with a wide range of samples and sample environments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Transient Lattice Response upon Photoexcitation in CuInSe 2 Nanocrystals with Organic or Inorganic Surface Passivation

CuInSe 2 nanocrystals offer promise for optoelectronics including thin-film photovoltaics and printed electronics. Additive manufacturing methods such as photonic curing controllably sinter particles into quasi-continuous films and offer improved device performance. To gain understanding of nanocrystal response under such processing conditions, we investigate impacts of photoexcitation on colloidal nanocrystal lattices via time-resolved X-ray diffraction. We probe three sizes of particles and two capping ligands (oleylamine and inorganic S 2- ) to evaluate resultant crystal lattice temperature, phase stability, and thermal dissipation. Elevated fluences produce heating and loss of crystallinity, the onset of which exhibits particle size dependence. We find size-dependent recrystallization and cooling lifetimes ranging from 90 to 200 ps with additional slower cooling on the nanosecond time scale. Sulfide-capped nanocrystals show faster recrystallization and cooling compared to oleylamine-capped nanocrystals. Using these lifetimes, we find interfacial thermal conductivities from 3 to 28 MW/(m 2 K), demonstrating that ligand identity strongly influences thermal dissipation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

An Iterative Response-Surface-Based Approach for Chance-Constrained AC Optimal Power Flow Considering Dependent Uncertainty

A modern power system is characterized by a stochastic variation of the loads and an increasing penetration of renewable energy generation, which results in large uncertainties in its states. These uncertainties bring formidable challenges to the power system planning and operation process. To address these challenges, we propose a cost-effective, iterative response-surface-based approach for the chance-constrained AC optimal power-flow problem that aims to ensure the secure operation of the power systems considering dependent uncertainties. Starting from a stochastic-sampling-based framework, we first utilize the copula theory to simulate the dependence among multivariate uncertain inputs. Then, to reduce the prohibitive computational time required in the traditional Monte-Carlo method, we propose, instead of using the original complicated power-system model, to rely on a polynomial-chaos-based response surface. This response surface allows us to efficiently evaluate the time-consuming power-system model at arbitrary distributed sampled values with a negligible computational cost. This further enables us to efficiently conduct an online stochastic testing for the system states that not only screens out the statistical active constraints, but also assists in a better design of the tightened bounds without using any Gaussian or symmetric assumption. Finally, an iterative procedure is executed to fine-tune the optimal solution that better satisfies a predefined probability. The simulations conducted in multiple test systems demonstrate the excellent performance of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

The role of surfaces in flexoelectricity

We develop the relationship between the strain derivative of the mean-inner potential and surface contributions to flexoelectricity, identifying the true surface-specific component of the flexoelectric response of finite samples. Density functional theory calculations on a range of experimentally observed, low energy SrTiO 3 , MgO, and Si surfaces demonstrate that the mean-inner potential and its contributions to flexoelectricity are sensitive to small differences in surface structure, chemistry, and adsorbates. We also introduce a method to estimate mean-inner potential contributions to flexoelectricity using electron scattering factors and use this approximation to predict total flexoelectric responses for a variety of insulators. Strategies to experimentally disentangle bulk and surface flexoelectric terms are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Micropolar Elastoplasticity Using a Fast Fourier Transform‐Based Solver

ABSTRACT This work presents a micromechanical spectral formulation for obtaining the full‐field and homogenized response of elastoplastic micropolar composites. A closed‐form radial‐return mapping is derived from thermodynamics‐based micropolar elastoplastic constitutive equations to determine the increment of plastic strain necessary to return the generalized stress state to the yield surface, and the algorithm implementation is verified using the method of numerically manufactured solutions. Then, size‐dependent material response and micro‐plasticity are shown as features that may be efficiently simulated in this micropolar elastoplastic framework. The computational efficiency of the formulation enables the generation of large datasets in reasonable computing times.

42 ENGINEERING↗

DOE OSTI2025

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↗

DOE OSTI2020

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↗

DOE OSTI2025

Surface Trap Dynamics and Gamma‐Ray Detection Enhancement in Passivated FAPbBr 3 Perovskite Crystals

Perovskite‐based direct radiation detectors offer a compelling platform for room‐temperature gamma spectroscopy due to their high sensitivity, tunable optoelectronic properties, and compatibility with solution processing. While prior studies on formamidinium lead tribromide (FAPbBr 3 ) have shown that surface passivation mitigates deep trap state impact on the charge collection efficiency (CCE), the influence of surface‐based shallow traps on time‐resolved detector response remains largely unexplored. In this article, a digital signal processing (DSP) method is applied to analyze waveform rising edge dynamics in FAPbBr 3 single‐crystal detectors, identifying prompt and delayed charge collection components. Using the second derivative of the collected charge signal, the onset of shallow trap reemission is isolated, occurring 1.1 µs after carrier generation. Surface passivation suppressed deep trap states, shifting carrier occupancy toward shallower traps and extending the temporal window of delayed charge collection. This delayed contribution correlates with measurable changes in the temporal and amplitude distributions of detector events. By leveraging time‐domain filtering informed by these dynamics, a selective enhancement in signal quality and gamma peak‐to‐Compton ratio (PCR) is achieved, while preserving over 90% of photoelectric events. An improvement in the PCR to 1.80, from the original 1.09, is enabled by focusing on high collection efficiency events corresponding to reemitted carriers in the 1.1 to 3 µs range after the generation event. These results demonstrate a pathway for leveraging shallow trap dynamics to enhance spectral fidelity in time‐resolved perovskite detectors.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Understanding the hydrological response of a headwater-dominated catchment by analysis of distributed surface–subsurface interactions

Abstract We computationally explore the relationship between surface–subsurface exchange and hydrological response in a headwater-dominated high elevation, mountainous catchment in East River Watershed, Colorado, USA. In order to isolate the effect of surface–subsurface exchange on the hydrological response, we compare three model variations that differ only in soil permeability. Traditional methods of hydrograph analysis that have been developed for headwater catchments may fail to properly characterize catchments, where catchment response is tightly coupled to headwater inflow. Analyzing the spatially distributed hydrological response of such catchments gives additional information on the catchment functioning. Thus, we compute hydrographs, hydrological indices, and spatio-temporal distributions of hydrological variables. The indices and distributions are then linked to the hydrograph at the outlet of the catchment. Our results show that changes in the surface–subsurface exchange fluxes trigger different flow regimes, connectivity dynamics, and runoff generation mechanisms inside the catchment, and hence, affect the distributed hydrological response. Further, changes in surface–subsurface exchange rates lead to a nonlinear change in the degree of connectivity—quantified through the number of disconnected clusters of ponding water—in the catchment. Although the runoff formation in the catchment changes significantly, these changes do not significantly alter the aggregated streamflow hydrograph. This hints at a crucial gap in our ability to infer catchment function from aggregated signatures. We show that while these changes in distributed hydrological response may not always be observable through aggregated hydrological signatures, they can be quantified through the use of indices of connectivity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Chemometric optimisation of pressurised liquid extraction for the determination of alliin and S‐allyl‐cysteine in giant garlic ( Allium ampeloprasum L.) by liquid chromatography tandem mass spectrometry

Abstract Introduction Giant garlic is a functional food that contains different kinds of bioactive molecules with beneficial effects on chronic noncommunicable diseases like diabetes and cardiovascular conditions. Considering biosynthesis pathways, abundance, and biological activity, alliin and S‐allyl‐cysteine were used as chemical markers of organosulphur compounds present in giant garlic. Objective To establish a chemometric optimisation of pressurised liquid extraction for the determination of alliin and S‐allyl‐cysteine in giant garlic by liquid chromatography tandem mass spectrometry (LC–MS/MS). Methodology Samples were blanched (ca. 90°C for 10 min) to inactivate alliinase and γ‐glutamyl transpeptidase enzymes and then freeze‐dried. Chemometric optimisation was performed via response surface methodology based on central composite design (CCD). Organosulphur compound yields were determined applying a validated LC–MS/MS method in multiple reaction monitoring (MRM) mode using the following transitions: for alliin m/z 178 → 74 and for S‐allyl‐cysteine m/z 162 → 41. Results According to CCD results, under constant conditions of pressure (1500 psi) and time (20 min), the optimal conditions for pressurised liquid extraction of alliin and S‐allyl‐cysteine were 70.75 and 68.97% v/v of ethanol in extraction solvent and 76.45 and 98.88°C as extraction temperature, respectively. Multiple response optimisation for the simultaneous extraction of both organosulphur compounds was established via desirability function. Under these conditions, 2.70 ± 0.27 mg g −1 dry weight (DW) of alliin and 2.79 ± 0.22 mg g −1 DW of S‐allyl‐cysteine were extracted. Conclusions These results clearly demonstrated that pressurised liquid extraction is an efficient green technique to extract bioactive organosulphur compounds from giant garlic. Extraction yields were significantly ( p < 0.05) higher than those obtained with conventional ultra‐turrax extraction.

Peterssen‐Fonseca, Darlene↗

DOE OSTI2022

Laser Doppler vibrometry for piezoelectric coefficient ($d_{33}$) measurements in irradiated aluminum nitride

Sensors used for experiments in advanced reactors must survive in harsh environments. Few material systems can be used to construct sensors viable for extreme conditions. Aluminum nitride (AlN) is one such material because it has high thermal stability and radiation resistance and sustains good piezoelectric and dielectric properties at high temperatures. In this work, the piezoelectric coefficient $d_{33}$ of the AlN single-crystal with thermal and irradiation damage was investigated using an indirect method with a laser Doppler vibrometer (LDV). Surface electrodes were deposited on the AlN samples, and the vibration response of the samples to an applied voltage was monitored using the LDV as a function of the excitation frequency. The $d_{33}$ estimation was based on the excitation voltage and the thickness-mode displacement extracted from the LDV measurements. Further, six AlN substrates were irradiated with 8 MeV Al 2+ at three fluences (10 15 , 10 16 , and 10 17 ions/cm 2 ) and two temperatures (300 °C and 500 °C). The $d_{33}$ for the six irradiated samples and one pristine sample were measured, and the measurement uncertainty was estimated based on five repeated tests. All samples were also measured by a commercial piezometer for comparison. The experimental results demonstrate that the piezoelectric coefficients obtained by LDV were about 0.8–1.16 pm/V lower than those obtained by the piezometer. With the compensation of the clamping effect, the corrected LDV values are similar to the piezometer results. Both show similar trends in all samples, which validates the feasibility of the proposed method for $d_{33}$ measurement. Based on the LDV results, the irradiated samples show a 12%–22% decrease in $d_{33}$ compared with the pristine samples. The samples irradiated under the same fluence at a higher temperature (500 °C) demonstrated a lower $d_{33}$ than those at 300 °C. The effect of the retro-reflective tape, sample temperature, and sample size on the $d_{33}$ measurement were also studied.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Functional stimuli-responsive polymers on micro- and nano-patterned interfaces

Micro- and nano-patterned surfaces offer precise control over morphology and chemical composition, enhancing the stability, durability, and functionality of coating materials. When combined with stimuli-responsive polymers, these surfaces gain dynamic adaptability, enabling reversible binding, reusable sensing, and selective molecular capture. Furthermore, while recent review articles have explored various aspects of stimuli-responsive materials, from hydrogel patterns for bioanalytical applications to shape-morphing hydrogels for soft robotics and sensors, a comprehensive review focused on the integration of smart polymers with micro- or nano-patterned interfaces remains absent. This review addresses key surface patterning techniques, including soft lithography, colloidal lithography, and polymer brush photolithography, as well as advances in surface-initiated polymerization methods, such as surface-initiated controlled radical polymerization (SI-CRP). In addition, we discuss recent progress in integrating stimuli-responsive polymers with patterned surfaces to create advanced, functional materials.

Colloidal lithography↗

DOE OSTI2024

Efficient learning of accurate surrogates for simulations of complex systems

Machine learning methods are increasingly deployed to construct surrogate models for complex physical systems at a reduced computational cost. However, the predictive capability of these surrogates degrades in the presence of noisy, sparse or dynamic data. Here, we introduce an online learning method empowered by optimizer-driven sampling that has two advantages over current approaches: it ensures that all local extrema (including endpoints) of the model response surface are included in the training data, and it employs a continuous validation and update process in which surrogates undergo retraining when their performance falls below a validity threshold. We find, using benchmark functions, that optimizer-directed sampling generally outperforms traditional sampling methods in terms of accuracy around local extrema even when the scoring metric is biased towards assessing overall accuracy. Finally, the application to dense nuclear matter demonstrates that highly accurate surrogates for a nuclear equation-of-state model can be reliably autogenerated from expensive calculations using few model evaluations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Efficient Learning of Accurate Surrogates for Simulations of Complex Systems

Machine learning methods are increasingly deployed to construct surrogate models for complex physical systems at a reduced computational cost. However, the predictive capability of these surrogates degrades in the presence of noisy, sparse, or dynamic data. We introduce an online learning method empowered by optimizer-driven sampling that has two advantages over current approaches: it ensures that all local extrema (including endpoints) of the model response surface are included in the training data and it employs a continuous validation and update process in which surrogates undergo retraining when their performance falls below a validity threshold. We find, using benchmark functions, that optimizer-directed sampling generally outperforms traditional sampling methods in terms of accuracy around local extrema even when the scoring metric is biased towards assessing overall accuracy. Finally, the application to dense nuclear matter demonstrates that highly accurate surrogates for a nuclear equation of state model can be reliably auto-generated from expensive calculations using few model evaluations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Calibrating hypersonic turbulence flow models with the HIFiRE-1 experiment using data-driven machine-learned models.

In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model’s input parameters. Here we develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model’s response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.

42 ENGINEERING↗

DOE OSTI2021

Design and Optimization of a Gas-Cooled, Airfoil Fin Microchannel Heat Exchanger

High-performance microchannel heat exchangers are needed to supply heat for power conversion for nuclear microreactors. An airfoil fin microchannel design, constructed of Alloy 617 with helium as the working fluid, was analyzed and optimized using a design of experiments with artificial intelligence and machine learning techniques. The use of airfoil fins offers the potential to reduce pressure drop across the heat exchanger, as compared to other types of channel configurations. A framework for topology optimization of airfoil fin PCHEs has been developed that can be readily extended to different fin sizes and shapes, as well as different inlet and operating conditions, materials of construction, and working fluids. An optimization procedure was developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. STAR-CCM+ was used to analyze a simplified two-channel configuration where five parameters were varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods were compared for generating surrogate models – a 5D polynomial and a regression neural network. A response surface approximation was created from the surrogate models and input to a genetic algorithm. The genetic algorithm identified a set of optimal points on the Pareto front. The optimal geometry was found across six channel Reynolds numbers ranging from 1000 to 5000 to analyze how varying inlet conditions affects the optimal design. A set of optimal designs that maximizes heat transfer and minimizes pressure drop was identified, and a thermal stress analysis was performed on the optimal design. This work has developed a digital framework for the expedient topology design and evaluation of PCHE designs for gas-cooled microreactor applications. Correlations for the Nusselt number and Darcy friction factor were developed that can be useful for thermal hydraulic analyses using system codes. A thermal stress analysis was conducted and a brief discussion of the status of code cases of PCHEs for nuclear applications is given. Testing and thermomechanical modeling is needed to facilitate future code compliance of PCHEs for high pressure and high temperature applications.

42 ENGINEERING↗

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 OSTI2020

Multi-objective optimization of peel and shear strengths in ultrasonic metal welding using machine learning-based response surface methodology

Ultrasonic metal welding (UMW) is a solid-state joining technique with varied industrial applications. Despite of its numerous advantages, UMW has a relative narrow operating window and is sensitive to variations in process conditions. As such, it is imperative to quantitatively characterize the influence of welding parameters on the resulting joint quality. The quantification model can be subsequently used to optimize the parameters. Conventional response surface methodology (RSM) usually employs linear or polynomial models, which may not be able to capture the intricate, nonlinear input-output relationships in UMW. Furthermore, some UMW applications call for simultaneous optimization of multiple quality indices such as peel strength, shear strength, electrical conductivity, and thermal conductivity. To address these challenges, this paper develops a machine learning (ML)-based RSM to model the input-output relationships in UMW and jointly optimize two quality indices, namely, peel and shear strengths. The performance of various ML methods including spline regression, Gaussian process regression (GPR), support vector regression (SVR), and conventional polynomial regression models with different orders is compared. A case study using experimental data shows that GPR with radial basis function (RBF) kernel and SVR with RBF kernel achieve the best prediction accuracy. The obtained response surface models are then used to optimize a compound joint strength indicator that is defined as the average of normalized shear and peel strengths. In addition, the case study reveals different patterns in the response surfaces of shear and peel strengths, which has not been systematically studied in the literature. While developed for the UMW application, the method can be extended to other manufacturing processes.

42 ENGINEERING↗